Next Article in Journal
Mathematical Tools that Connect Different Indexing Analyses
Previous Article in Journal
In Vitro Bioaccessibility of Citrus Pomace Compounds Possessing Health Promoting Properties with Potential to Reduce the Risk of Diabetes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Systematic Model-Based Steady State and Dynamic Optimization of Combined Cooling and Antisolvent Multistage Continuous Crystallization Processes †

Department of Chemical Engineering, Loughborough University, Loughborough LE11 3TU, UK
*
Author to whom correspondence should be addressed.
Presented at the 2nd International Electronic Conference on Crystals, 10–20 November 2020; Available online: https://iocc_2020.sciforum.net/.
Proceedings 2020, 62(1), 7; https://doi.org/10.3390/proceedings2020062007
Published: 31 December 2020
(This article belongs to the Proceedings of The 2nd International Online Conference on Crystals)

Abstract

:
Currently, one of the key challenges in the pharmaceutical industry is the transformation of traditional batch production methods into robust continuous processes with the intention of reducing manufacturing costs and time and improving product quality. Crystallization is by far the most important purification technology in Pharma, as more than 80% of the active pharmaceutical ingredients (API) require at least one crystallization step. A successful crystallization process requires tight control over crystal size, shape and polymorphic purity. A rigorous and systematic methodology is presented to design and optimize multistage combined cooling and antisolvent continuous (mixed-suspension, mixed-product removal- MSMPR) crystallizers. The crystallization of acetylsalicylic acid (API) in ethanol (solvent) and water (anti-solvent) is used as a case study. A predictable and validated mathematical model of the system, which consists of a one-dimensional population balance model, was used to develop several optimizations strategies. Firstly, the attainable region of the mean particle size was determined for both minimum and maximum attainable crystal size. The method helped identify the most suitable number of stages and total residence time or volume for a cascade of continuous crystallizers. This was followed by a steady state optimization which helped determine the optimal operating temperatures and antisolvent flowrates. To minimize the startup time, a series of dynamic optimization strategies were implemented, assuming starting from empty vessels. The optimal dynamic profiles of the temperature and antisolvent flow rate, at different crystallization steps, were identified using a systematic and rigorous approach allowing a reduction in the startup time by 31%.

1. Introduction

Continuous manufacturing is increasingly seen as the most flexible and viable option for the pharmaceutical industry [1]. However, despite the strides forward, many technical challenges are still to be overcome. Continuous pharmaceutical campaigns are anticipated to have short operating windows [2]. This makes the impact of startup and shut down extremely important in both cost of production and environmental impact. To address some of these key issues, the next-generation pharmaceutical plants require systematic, rigorous and robust optimal strategies for process design, operation, and control for single processes and integrated plants [2].
Crystallization is widely used in the pharmaceutical industry for its outstanding efficiency and economic performance in purifying active pharmaceutical ingredients (API). The quality of the final product is commonly determined by the crystal size distribution, shape, and purity, which impact the downstream processability and drug safety and efficacy. Commonly, the higher average crystal size and lower coefficient of variation (CV) mean higher quality in the pharmaceutical industry [3].
The fundamental driving force for crystallization is supersaturation. Typically, the supersaturation is generated by cooling, solvent evaporation, or the addition of an antisolvent. Most of the crystallization literature is focused on the design, optimization, and control of batch processes. Various model-based and model-free techniques have been extensively implemented and thoroughly discussed in the literature. However only a few decision and control variables are available to optimally design and operate a batch crystallization process. The most common approach is to use the temperature or/and antisolvent ratio to optimize the crystal size distribution, CV, and batch time [4].
Over the last decade, continuous crystallization of API has received a growing interest. To date, most of the literature has focused on three main continuous crystallizers: mixed suspension mixed product removal (MSMPR), plug flow crystallizers and continuous oscillatory baffled crystallizers. Currently, the most dominant type of crystallizer in the pharmaceutical industry is based on a stirred tank design which exploits well-established theory and know-how built over decades of batch crystallization. Consequently, many experimental and modelling efforts have been devoted to continuous MSMPR crystallizers. Several studies focused on the optimization of single, multistage MSMPR and integrated end-to-end continuous pharmaceutical plants with a series of MSMPR crystallizers [1,5,6,7]. Despite the importance of the startup and its anticipated impact on cost and environmental footprint, only a few investigations are reported in the literature, which identified mainly the optimal dynamic profile of the cooling temperature and antisolvent flowrate [8]. However, the MSMPR vessels were assumed to be prefilled with a solution with the steady state composition. This approach is unrealistic and disregards the impact of filling the vessels during the startup, which may extend the startup time even further. To date, only Benyahia and coworkers showed the potential impact of different startup scenarios, including starting from empty MSMPR vessels, on continuous pharmaceutical campaigns [2]. The latter laid the ground for the current investigation, which consists of a systematic optimization approach starting from steady state optimization, aimed at the identification of the attainable mean particle size, followed by different startup optimization strategies including startup from empty vessels.

2. Materials and Methods

The combined cooling and antisolvent crystallization of aspirin (acetylsalicylic acid, ASA) in ethanol (solvent) and water (antisolvent) was considered as a case study. A dynamic mathematical model of the continuous crystallization process (MSMPR) was derived in conjunction with the literature [8], based on several key assumptions summarized below:
  • All vessels are assumed to be well-mixed.
  • The total volume of the solvent and antisolvent is additive.
  • The secondary nucleation, breakage and agglomeration are negligible.
The population balance equations (PBE) of the first and i t h   stages can be written as follows:
V i n i ( L , t ) t = F e e d   n f e e d ( L , t ) F o u t , i n i ( L , t ) V i G i n i ( L , t ) L     i = 1
V i n i ( L , t ) t = F o u t , i 1 n i 1 ( L , t ) F o u t , i n i ( L , t ) V i G i n i ( L , t ) L     i = 2 , , N
where V i refers to the volume of i t h stage vessel, feed is the inlet flow rate to the first crystallizer, n f e e d ( L , t ) is the number density function of the crystals in the feed solution, F o u t , i is the outlet flow rate of i t h crystallization stage, and G i is the growth rate of API crystals.
By using the standard method of moment, the   j t h moment at the i t h stage is defined by:
μ j , i ( t ) = 0 L j n i ( L , t ) d L ,     j = 0 , 1 , ,
Assuming that the growth rate G i is size-independent and the volume V i is constant. Equations (1) and (2) yield:
V i d μ j , i d x = F o u t , i 1 μ j , i 1 F o u t , i μ j , i + V i G i j μ j 1 , i ,     j 1
d μ j , i d x = F o u t , i 1 V i μ j , i 1 F o u t , i V i μ j , i + G i j μ j 1 , i ,     j 1
In this case study, the growth rate and nucleation rate kinetics, G and B respectively, are adopted from the literature [9] with the kinetic parameters outlined in Table 1.
G = k G 1 e x p ( k G 2 R T ) ( c * ( S 1 ) ) k G 3
B = k B 1 e x p ( k B 2 R T ) e x p ( k B 3 l n 2 S )
The supersaturation (S) is defined by Equation (8), where c is the concentration of the API at a given time expressed as mass ratio (Equation (9)) and c * is the concentration at saturation calculated by a polynomial model adopted from the literature [10], and the coefficients are shown in Table 2.
The solubility curve of the API (acetylsalicylic acid) is described by a polynomial equation:
c * =   P 11 + P 21 × W + P 12 × T + P 31 × W 2 + P 22 × W T + P 13 × T 2 +   P 41 × W 3 + P 32 × W 2 T + P 23 × W T 2 + P 51 × W 4 + P 42 × W 3 T +   P 33 × W 2 T 2 + P 61 × W 5 + P 52 × W 4 T + P 43 × W 3 T 2
              S = c c *
c = M A P I M s o l v e n t + M a n t i s o l v e n t
The mass balance of API, solvent and antisolvent are shown in Equations (11)–(16), where ρ f e e d is the density of feed, ρ i is the density of the solution density at the i t h stage, ρ c is the density of the crystal, which is 1400 g/kg, and k v is the shape factor ( π / 6 ) .
d M A P I , i d t = F e e d   ρ f e e d ( ω s , f e e d + ω a s , f e e d ) c F e e d F o u t , 1 ρ 1 ( ω s ,   1 + ω a s , 1 ) c 1   3 ρ c k v G 1 μ 2 , 1 V 1 , i = 1  
d M A P I , i d t = F o u t , i 1 ρ i 1 ( ω s ,   i 1 + ω a s , i 1 ) c i 1 F o u t , i , ρ i ( ω s ,   i + ω a s , i ) c i   3 ρ c k v G i μ 2 , i V i ,     i = 2 , 3
d M s , i d t = F e e d ρ f e e d ω s , f e e d F o u t , i ρ i ω s ,   i ,     i = 1
d M s , i d t = F o u t , i 1 ρ i 1 ω s ,   i 1 F o u t , i ρ i ω s ,   i ,     i = 2 , 3
d M a s , i d t = F e e d ρ f e e d ω a s , f e e d F o u t , i ρ i ω a s ,   i + F a s , i ρ a s ,     i = 1
d M a s , i d t = F o u t , i 1 ρ i 1 ω a s ,   i 1 F o u t , i ρ i ω a s ,   i + F a s , i ρ a s ,     i = 2 , 3

3. Results and Discussion

3.1. Steady State Optimization

The objective of the steady state optimization is to generate the optimal temperature ( T i ), antisolvent flow rate ( F a s , i ) and residence time ( R i ) at each stage, which maximize the average crystal size ( L n ). A three-stage continuous crystallization platform is considered where the supersaturation is created by the combination of cooling and antisolvent addition (Figure 1)
The feed temperature is 40 °C and the minimum and final temperature are both set at 25 °C to generate enough supersaturation. This is captured by the constraints C1 and C7. The concentration of antisolvent is subject to an upper bound which depends on the solubility performance. For the current system, ASA in ethanol (solvent) and water (antisolvent), the solubility of the API decreases as the antisolvent (water) weight ratio increases from 25% to 70%. The maximum antisolvent flow rate in each stage is set at 18.25 g/min, as expressed by the constraint C5. Similarly, the total residence time should not exceed 30 min, while the minimum residence time in each stage is not less than 1 min (constraints C2 and C6). In addition, the constraint C4 ensures that the temperature in the ( i + 1 ) th stage is always lower or equal to the temperature of the i th stage, which prevents undersaturated conditions and helps avoid heating which may lead to the dissolving the API crystals.
The mathematical formulation of the steady state single objective optimization problem discussed above is described below:
M a x T i , R i , F a s i d = μ 1 , 3 μ 0 , 3   i = 1 , 2 , 3
s . t .   M o d e l   e q u a t i o n s   5 16
C 1 :   25 T i 40
C 2 :   1 R i 30
C 3 :   0   F a s , i 18.5
C 4 :   T i T i 1
C 5 :   i = 1 3   F a s , i = 18.25
C 6 :   i = 1 3 R i = 30
C 7 :   T 3 = 25
where T i ,   R i ,   F a s , i are the temperature, residence time and antisolvent flow rate of the i t h stage, respectively.
This optimization problem was solved by a genetic algorithm using MATLAB with a population size of 50 and maximum number of generations of 500. The maximum average crystal size obtained was 392   μ m , which corresponds to the optimal temperatures, antisolvent flowrates, and residence times outlined in Table 3. The performance of the genetic algorithm can be evaluated by using the difference between the best solution and the mean value in one generation, as shown in Figure 2a. This approach also provides a convergence (stop) criterion for the genetic algorithm.
As shown in Figure 2.b, the dynamic profiles of the mean crystal size at different crystallization stages shows that the crystal size nearly doubled at the last stage compared to the first stage, due the high supersaturation generated in stages 2 and 3. The first crystallizer has a short residence time and was dominated by nucleation, making this stage ideal for seed generation, despite some growth exhibited during this stage. Although the dynamic mathematical model was used for the sake of steady state optimization, the results also show the intrinsic startup dynamics during which the crystals produced are off-specification and will be considered as waste. This wasted material represents a significant loss, particularly when considering a short-operation steady state window as anticipated in the pharmaceutical industry [6,11].

3.2. Startup Optimization

As discussed earlier, operating at the steady state optimal conditions (i.e., optimal temperatures and antisolvent flow rates) will ensure the attainment of the maximum mean crystal size at steady state. However, this approach provides no means of control over the intrinsic dynamics exhibited during the startup period during which large amounts of off-specification (waste) material is produced. To minimize the startup time and maximize productivity and reduce wastes, a dynamic optimization approach was developed as described in this section. To allow the dynamic optimization of a process, it is common to optimize piecewise constant and piecewise continuous profiles of the control variables. The approach has been widely used in polymer engineering and related research and industrial sectors [12]. However, only a few startup optimization approaches have been applied to the continuous crystallization processes, as reported in the literature [8]. The reported work focused mainly on one or two crystallization stages and considered the crystallizers as being entirely prefilled. Here, a more systematic and rigorous approach was developed to optimize the startup time of a three-stage MSMPR using piecewise constant profiles of antisolvent flowrate and piecewise continuous profiles of temperatures. The approach also addresses the startup optimization with initially empty vessels. For the sake of brevity, only the optimization results associated with the antisolvent flowrates at different stages are presented here.
The dynamic optimization aims at minimizing the t 95 , which represents the time to reach the steady state mean crystal size (obtained from steady state optimization) within a relative error band of ± 5 % . The mathematical formulation of the optimization problem is described below.
M i n   F a s , i , k t 95 ,   i = 1 ,   2 ,   3 ,   k = 1 , 2 , 5  
s . t .   C 1 :   0   F a s , 1 , k 4 ,   k = 1 ,   2 , 5   C 2 :   0   F a s , 2 , k 15.2 ,   k = 1 , 2 , 5   C 3 :   0   F a s , 3 , k 20 ,   k = 1 , 2 , 5
M o d e l l   e q u a t i o n s   5 16
The piecewise constant (discrete) profile of the antisolvent flowrate consists of five equal time intervals of 10 min. Only the results associated with the antisolvent flow rate, used as decision variables for the startup optimization, are presented here. The three vessels were all kept empty at the beginning. For each stage, the flow rate of antisolvent depends on the feed and the outlet flow rates.
In the optimization problem above,   F a s , i , k is the antisolvent flow rate of i t h stage and k t h time intervals. The antisolvent flowrate upper and lower bounds were set in C1–3. The optimization problem was solved using a genetic algorithm in MATLAB. The population and other settings are set the same as for the steady state optimization. The results are summarized in Table 4 below.
The optimal dynamic profiles of the average crystal size at different stages are shown in Figure 3. With the optimized antisolvent profile, the startup time is shortened to 63.43 min. Compared to the unoptimized startup (startup time of 86.62 min), the optimization saved 23.2 min (26.77% performance improvement) based on a startup from empty vessels and using equally spaced piece-wise antisolvent flowrate profiles. When the time interval number is increased to 10, the performance becomes even better. The startup time is shortened to 61.09 min (31.3% performance improvement). t95 was used as an example of the objective function. However, more tight quality requirements can be achieved by using t99. This is particularly useful when the purity is considered as one the critical quality attributes alongside mean crystal size.

4. Conclusions

A systematic and rigorous model-based optimization approach was implemented to the combined cooling and antisolvent crystallization of ASA in continuous MSMPR cascade crystallizers. The steady state optimization approach, based on the attainable regions, made it possible to design a cascade of three MSMPR crystallizers. In addition, it was possible to identify the optimal temperature, antisolvent flowrate and residence time for each MSMPR.
The minimization of startup is a major challenge in continuous pharmaceutical manufacturing due to the anticipated short operating windows. In order to minimize the startup time of the cascade of MSMPR crystallizers, which in turn maximizes the efficiency by reducing waste and increasing on-specification products, several dynamic optimization strategies were implemented. More realistic scenarios were considered by starting from empty MSMPR vessels. To address the dynamic optimization problem, several decision/control variables and discretization rules were adopted to enhance the available degrees of freedom and achieve a more effective optimization. The optimal dynamic profile of cooling temperature, antisolvent flowrate, and seed addition were identified and compared. Only the optimal profile of antisolvent flowrates is presented here for the sake of brevity. It was shown that the startup optimization improves the startup time by 26.77% and 31.3%, which represents a significant gain in terms of production and waste minimization.

References

  1. Mascia, S.; Heider, P.L.; Zhang, H.; Lakerveld, R.; Benyahia, B.; Barton, P.I.; Braatz, R.D.; Cooney, C.L.; Evans, J.M.B.; Jamison, T.F.; et al. End-to-End Continuous Manufacturing of Pharmaceuticals: Integrated Synthesis, Purification, and Final Dosage Formation. Angew. Chem. Int. Ed. 2013, 52, 12359–12363. [Google Scholar] [CrossRef] [PubMed]
  2. Benyahia, B.; Lakerveld, R.; Barton, P.I. A Plant-Wide Dynamic Model of a Continuous Pharmaceutical Process. Ind. Eng. Chem. Res. 2012, 51, 15393–15412. [Google Scholar] [CrossRef]
  3. Fysikopoulos, D.; Benyahia, B.; Borsos, A.; Nagy, Z.K.; Rielly, C.D. A Framework for Model Reliability and Estimability Analysis of Crystallization Processes with Multi-Impurity Multi-Dimensional Population Balance Models. Comput. Chem. Eng. 2019, 122, 275–292. [Google Scholar] [CrossRef]
  4. Nagy, Z.K.; Fujiwara, M.; Braatz, R.D. Modelling and Control of Combined Cooling and Antisolvent Crystallization Processes. J. Process Control 2008, 18, 856–864. [Google Scholar] [CrossRef]
  5. Lakerveld, R.; Benyahia, B.; Heider, P.L.; Zhang, H.; Wolfe, A.; Testa, C.J.; Ogden, S.; Hersey, D.R.; Mascia, S.; Evans, J.M.B.; et al. The Application of an Automated Control Strategy for an Integrated Continuous Pharmaceutical Pilot Plant. Org. Process Res. Dev. 2015, 19, 1088–1100. [Google Scholar] [CrossRef]
  6. Benyahia, B. Applications of a Plant-Wide Dynamic Model of an Integrated Continuous Pharmaceutical Plant: Design of the Recycle in the Case of Multiple Impurities, 1st ed.; Elsevier: Amsterdam, The Netherlands, 2018; Volume 41. [Google Scholar]
  7. Su, Q.; Benyahia, B.; Nagy, Z.K.; Rielly, C.D. Mathematical Modeling, Design, and Optimization of a Multisegment Multiaddition Plug-Flow Crystallizer for Antisolvent Crystallizations. Org. Process Res. Dev. 2015, 19, 1859–1870. [Google Scholar] [CrossRef]
  8. Yang, Y.; Nagy, Z.K. Combined Cooling and Antisolvent Crystallization in Continuous Mixed Suspension, Mixed Product Removal Cascade Crystallizers: Steady-State and Startup Optimization. Ind. Eng. Chem. Res. 2015, 54, 5673–5682. [Google Scholar] [CrossRef]
  9. Lindenberg, C.; Krättli, M.; Cornel, J.; Mazzotti, M.; Brozio, J. Design and Optimization of a Combined Cooling/Antisolvent Crystallization Process. Cryst. Growth Des. 2008, 9, 1124–1136. [Google Scholar] [CrossRef]
  10. Barik, K.; Prusti, P.; Mohapatra, S.S. Single- and Multi-Objective Optimisation for a Combined Cooling and Antisolvent Semi-Batch Crystallisation Process with an ACADO Toolkit. Indian Chem. Eng. 2020, 62, 287–300. [Google Scholar] [CrossRef]
  11. Parekh, R.; Benyahia, B.; Rielly, C.D. A Global State Feedback Linearization and Decoupling MPC of a MIMO Continuous MSMPR Cooling Crystallization Process; Elsevier: Amsterdam, The Netherlands, 2018; Volume 43. [Google Scholar]
  12. Benyahia, B.; Latifi, M.A.; Fonteix, C.; Pla, F. Multicriteria dynamic optimization of an emulsion copolymerization reactor. Comput. Chem. Eng. 2012, 35, 2886–2895. [Google Scholar] [CrossRef]
Figure 1. Cascade of three MSMPR crystallizers.
Figure 1. Cascade of three MSMPR crystallizers.
Proceedings 62 00007 g001
Figure 2. (a) Performance of the genetic algorithm showing the mean and best results of the objective function at different generations. (b) Optimal mean crystal size at different MSMPR stages.
Figure 2. (a) Performance of the genetic algorithm showing the mean and best results of the objective function at different generations. (b) Optimal mean crystal size at different MSMPR stages.
Proceedings 62 00007 g002
Figure 3. Dynamic profiles of the mean crystal size at different crystallization stages during startup.
Figure 3. Dynamic profiles of the mean crystal size at different crystallization stages during startup.
Proceedings 62 00007 g003
Table 1. kGi and kBi are empirical parameters.
Table 1. kGi and kBi are empirical parameters.
kGkB
1 ( 3.21 ± 0.18 ) × 10 4   ms 1 ( 1.15 ± 0.51 ) × 10 21   m 3 s 1
2 ( 2.58 ± 0.14 ) × 10 4   Jmol 1 ( 7.67 ± 0.11 ) × 10 4   Jmol 1
3 1.00 ± 0.01 0.16 ± 0.01
Table 2. Coefficients of the solubility polynomial equation
Table 2. Coefficients of the solubility polynomial equation
Pij123456
139.94.851−0.17040.007273−0.0001213 5.949 × 10 7
23.7140.2022−0.01658 7.666 × 10 5 1.14 × 10 6
30.1693−0.0036340.0003186 3.706 × 10 6
Table 3. Optimal temperatures, antisolvent flowrates, and residence times obtained in the case of the maximization of the average crystal size.
Table 3. Optimal temperatures, antisolvent flowrates, and residence times obtained in the case of the maximization of the average crystal size.
Stage 1Stage 2Stage 3
Temperature (°C)25.1125.0925.00
Antisolvent flow rate (g/min)2.057.598.61
Residence time (min)3.7710.2016.03
Table 4. Optimal antisolvent flow rate profiles obtained by the dynamic optimization approach.
Table 4. Optimal antisolvent flow rate profiles obtained by the dynamic optimization approach.
  F a s , 1 ( g / min )   F a s , 2 ( g / min )   F a s , 3 ( g / min )
0–10 min0.031.5216.47
10–20 min0.209.4813.06
20–30 min0.395.6219.80
30–40 min0.149.8519.88
40–50 min0.259.6219.32
50–end (steady state setting)2.057.598.61
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, J.; Benyahia, B. Systematic Model-Based Steady State and Dynamic Optimization of Combined Cooling and Antisolvent Multistage Continuous Crystallization Processes. Proceedings 2020, 62, 7. https://doi.org/10.3390/proceedings2020062007

AMA Style

Liu J, Benyahia B. Systematic Model-Based Steady State and Dynamic Optimization of Combined Cooling and Antisolvent Multistage Continuous Crystallization Processes. Proceedings. 2020; 62(1):7. https://doi.org/10.3390/proceedings2020062007

Chicago/Turabian Style

Liu, Jiaxu, and Brahim Benyahia. 2020. "Systematic Model-Based Steady State and Dynamic Optimization of Combined Cooling and Antisolvent Multistage Continuous Crystallization Processes" Proceedings 62, no. 1: 7. https://doi.org/10.3390/proceedings2020062007

APA Style

Liu, J., & Benyahia, B. (2020). Systematic Model-Based Steady State and Dynamic Optimization of Combined Cooling and Antisolvent Multistage Continuous Crystallization Processes. Proceedings, 62(1), 7. https://doi.org/10.3390/proceedings2020062007

Article Metrics

Back to TopTop