xuebaunion@vip.163.com

3551 Trousdale Rkwy, University Park, Los Angeles, CA

留学生论文指导和课程辅导

无忧GPA：https://www.essaygpa.com

工作时间：全年无休-早上8点到凌晨3点

微信客服：xiaoxionga100

微信客服：ITCS521

matlab代写-CGP006

时间：2021-03-14

CGP006- Advanced Process Design and Optimisation Loughborough University

1

B Benyahia

Advanced Process Design and Optimisation

To be submitted electronically via Learn by Friday 26th March 2021

before 16:00

Crystallisation is a key separation and purification technology in the pharmaceutical industry.

The driving force for crystallisation is supersaturation which is commonly generated by

cooling, addition of antisolvent or solvent evaporation. A batch cooling crystallisation process

is available to recover paracetamol from water. The initial concentration of paracetamol is

0.0256 g/g solvent and the initial temperature is 315 K. Supersaturation is generated by

cooling. Figure 1 illustrates the principle of cooling crystallisation. The system, which is initially

undersaturated at point A, is cooled down and crosses the solubility curve at point B. In fact,

nucleation does not occur at this stage and further cooling is required to reach the limit of the

metastable zone (point C), where the nucleation occurs spontaneously. The nuclei undergo

growth and as a result, supersaturation drops gradually bringing the system closer to the

solubility curve. Further cooling is then required to maintain higher supersaturation to help

the crystals grow larger while maximising solute (paracetamol) recovery. At the end, the batch

process is terminated at point D close to the solubility curve.

Figure 1. Phase diagram and principle of cooling crystallisation

The mathematical model of Paracetamol batch crystallisation is described below. The

standard method of moment is used for the population balance which captures all variations

in the population of crystals du to nucleation and growth. This population balance is critical

to predict the crystal size distribution over time.

0

= , (1)

1

= 0, (2)

CGP006- Advanced Process Design and Optimisation Loughborough University

2

B Benyahia

2

= 21, (3)

3

= 32 (4)

= − 32 (5)

Where is the ith moment and C the concentration of Paracetamol in water (g/g water).

is the shape factor ( = 0.24 ) and the density of the crystals ( = 1.296e6 g −3).

and are respectively growth and nucleation kinetics given by

= ( − ) (6)

= ( − ) (7)

Where , , , are the kinetic parameters and ( − ) is the supersaturation.

The concentration at saturation, (solubility curve) is given by the polynomial function

below:

= 1.5846e − 5 × 2 − 9.0567e − 3 × + 1.3066 (8)

where T is the process (solution) temperature.

The mathematical model above allows the prediction of the key features of the products

obtained. For example, the number based mean crystal size (d) and the coefficient of

variation (CV) can be obtained by

= 1

0

(9)

= �2012 − 1 (10)

Several experimental measurements of the mean crystal size and concertation are available

in table 1 below.

Table1. Experimental measurements of the mean crystal size and concentration

Time (min) 0 20 32 45 60 80

Mean size (m) 0 1.91E-05 2.71E-05 3.02E-05 3.15E-05 3.30E-05

Concentration

(g/g solvent) 0.0256 0.0221 0.0175 0.0149 0.0132 0.0125

CGP006- Advanced Process Design and Optimisation Loughborough University

3

B Benyahia

1. Simulate the process using random values of ( ∈ [2.018 , 5.020]) and

( ∈ [0.001 , 0.5]). Show and comment on your results.

[10 marks]

2. Estimate the kinetic parameters and using the experimental data available in

Table 1. Show that the prediction capability of the model is acceptable and use a

quantitative criterion to prove that the model fits well with the data.

We assume that the exponent parameters are given by = . and = .

[20 marks]

3. Develop a dynamic optimisation strategy to maximize the mean crystal size. Justify the

algorithm/method used and show and comment on your results.

Tips: to address the optimisation problem,

• Use at least 5 linear intervals for the dynamic cooling profile and use either the

temperature or the cooling rate as the vector of decision variables

• The final temperature of the batch crystallisation process must be 283 K

• Only cooling is allowed (no heating)

• The cooling rate bounds are: – 2 ℃/min and 0 ℃/min

[20 marks]

4. Using a similar approach as question 3, develop a dynamic optimisation strategy to

minimise the batch time.

[20 marks]

5. Develop a dynamic multiobjective optimisation (MOO) strategy by considering two

cost functions: maximising the mean crystal size and maximising yield. Formulate and

solve the MOO problem. Justify the algorithm/method used and show and comment

on the Pareto front. Suggest an operating point from the Pareto front and justify your

choice. Is it possible to consider an additional objective function? Justify your answer.

[30 marks]

Requirements: A written report with the worked solutions must be submitted electronically via

Learn by the submission deadline stated on CASPA. The report must include ORIGINAL figures

properly labelled and numbered (and cited in the text) with the mathematical formulation of all

optimisation problems. The minimum admissible font size is Times New Roman 11 or Arial 10,

single paragraph, minimum margins of 20 mm all round. There is a strict page limit of nine (8) A4

pages for the full report – any content beyond 8 pages will NOT be marked. Submit the report with

a front cover page including your student ID (which does not count to page limit). Please save your

file as a PDF file as “StudentID-CGP0062021” – this will help us saving time during marking. You

should submit all MATLAB codes as a zipped file under the same name “StudentID-CGP0062021”

学霸联盟