Tuning the Particle Sizes in Spherical Agglomeration - Crystal Growth

The drag coefficient CD can be estimated as a function of the Reynolds ... the drag force FD increases with the droplet size ddrop, whereas the interf...
0 downloads 0 Views 2MB Size
Subscriber access provided by University of South Dakota

Article

Tuning the Particle Sizes in Spherical Agglomeration Pawel M. Orlewski, Byeongho Ahn, and Marco Mazzotti Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.8b01134 • Publication Date (Web): 05 Sep 2018 Downloaded from http://pubs.acs.org on September 7, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

Tuning the Particle Sizes in Spherical Agglomeration Pawel M. Orlewski, Byeongho Ahn, Marco Mazzotti* Institute of Process Engineering, ETH Zurich, 8092 Zurich, Switzerland E-mail: marco.mazzotti@ipe.mavt.ethz.ch Phone: +41 44 632 24 56. Fax: +41 632 11 41

Abstract

This paper presents an experimental technique for tuning the size of spherical agglomerates together with a simple model that provides qualitative information on the agglomerate size. Spherical agglomerates of benzoic acid of varying sizes using toluene as the bridging liquid were obtained with a very good reproducibility. The bridging liquid was injected into the crystal suspension via a capillary to control the initial size of the bridging liquid droplets. Experimental results obtained at different operating conditions show a clear correlation between the initial droplet size and the agglomerate size. Agglomerate size noticeably decreases with reducing the capillary size. Horizontal injection of the bridging liquid produces smaller agglomerates as compared to those produced with a vertical injection, while the size of agglomerates is insensitive to the different injection positions used in this study. Agglomerate size decreases

ACS Paragon Plus Environment

1

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 34

along with the increasing stirring rate. The presented model adequately captures the behavior of the system and gives a qualitative agreement with the experimental results.

1

1. INTRODUCTION

2

Many chemical and pharmaceutical products are prepared as powders, often consisting of

3

particles made of functional materials. The morphology of those particles plays a crucial role, as

4

it can have an impact on the efficiency of downstream product handling (e.g., filtration and

5

drying) and on the biopharmaceutical properties of particles (e.g., dissolution and

6

bioavailability)1. Accordingly, a great deal of attention has been paid to modifying the crystal

7

morphology during the crystallization phase by controlling the process parameters, e.g.

8

supersaturation, cooling rate, stirring rate or type of additives. However, sometimes undesired

9

shapes or size distributions are inevitable, thus often forcing manufacturers to deal with micron-

10

sized needle-shaped or platelet-shaped particles that exhibit poor flowability. This leads to a

11

considerable decrease in the efficiency of the downstream processes such as filtration, blending

12

and compression (tableting).

13

Spherical crystallization has recently gained interest in the chemical and pharmaceutical

14

industries to address the flowability-related problems by producing crystals in a spherical form

15

with better micrometric properties such as flowability, packability, and compressibility compared

16

to the needle- or platelet-shaped crystals2. There are two types of spherical crystallization: quasi-

17

emulsion solvent diffusion (QESD) and spherical agglomeration. In the case of the former,

18

quasi-emulsion droplets of solvent in anti-solvent are formed, within which crystallization occurs

19

due to the counter diffusion of solvent out and anti-solvent into the droplets. The latter on the

20

other hand is induced by feeding a liquid binder (often called a bridging liquid), i.e. an

21

immiscible fluid that preferentially wets the crystals, into the suspension, often after

ACS Paragon Plus Environment

2

Page 3 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

1

crystallization, but also possibly during crystallization3. In such a process, the liquid binder is

2

dispersed into a suspension of primary crystals in the form of droplets, which capture primary

3

crystals by capillary action4. Due to its higher flexibility and simplicity, preferential

4

crystallisation is considered more promising than QESD3.

5

Numerous studies have examined the effect of various operating parameters on the spherical

6

agglomeration process: stirring rate5,6, stirring time4,6–8, initial size of primary crystals5,9, the

7

amount1,4,7,9 and the injection rate5,10 of the liquid binder. Among these parameters, often the

8

amount of the liquid binder, together with its injection rate, has been shown to be the most

9

critical parameter influencing the size of agglomerates6,8,11. Surprisingly though, the influence of

10

the initial binder dispersion, namely the initial size of binder droplets, on the size of spherical

11

agglomerates has not been studied so far.

12

The goal of this work is to understand the effect of the initial droplet size on the size of

13

spherical agglomerates by exploiting the advantages offered by a microfluidic device, which

14

enables the generation of a large quantity of binder droplets, whose size can be carefully tuned. It

15

also allows decoupling the features of the initial binder dispersion from other critical operating

16

parameters, thus enabling a better understanding of the role of the liquid binder in the spherical

17

agglomeration process.

18 19

2. TUNING THE SIZE OF DROPLETS: THEORY

20

Advances in microfluidics have facilitated droplet generation in a controlled manner, as

21

summarized in a recent review paper12. Among the numerous techniques for droplet generation,

22

we have employed the simplest method, i.e. injecting the dispersed phase into the continuous

23

(immiscible) liquid phase via a cylindrical microcapillary. In this method, under certain flow

ACS Paragon Plus Environment

3

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 34

1

conditions, the droplet forms at the capillary tip and further grows until it reaches a critical size,

2

after which the droplet starts to detach thus eventually leaving the capillary (see a high-speed

3

camera video of the whole process in the supplementary material available on-line). The size of

4

the detached droplet can vary depending on the capillary size, the flow field around the capillary

5

tip, and the properties of the fluids involved and of the capillary13. Since the droplet size is at the

6

heart of this work, the following section presents a simple model that has been used to estimate

7

such droplet size.

8

2.1 Droplet Size Estimation

9

There are four main forces acting on a droplet forming at the tip of a capillary, namely the

10

buoyancy force  , the kinetic force  , the drag force  , and the interfacial tension force  .

11

The aforementioned forces can be classified into two groups: lifting forces, that act to separate

12

the droplet from the tip, and restraining forces, that act to keep it attached to the tip14. If the flow

13

directions of the continuous and dispersed phase are assumed to be parallel (Figure 1), then the

14

only force keeping the droplet attached to the nozzle is the interfacial tension force, whereas the

15

remaining three forces act as lifting forces. As long as the restraining forces are larger than the

16

lifting forces, the droplet will stay attached to the tip and continue growing. As soon as the lifting

17

forces exceed the restraining forces, the droplet will start clearing the tip of the capillary14. The

18

size of droplets at the breaking point corresponds to the size of a droplet, for which the

19

restraining forces are balanced by the lifting forces. Below, the definitions of the four forces

20

acting on the droplet are provided.

21



The interfacial tension force  , that keeps the droplet at the capillary tip, is given by  =

 +   # 1 2

ACS Paragon Plus Environment

4

Page 5 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

1

where is the interfacial tension between the two fluids, and  and  is the inner and

2

outer diameter of the capillary respectively. Here, the perimeter of the triple interface has

3

been defined as  +  /2 instead of the usual  , because the actual perimeter of

4

the triple-phase contact line among the two fluids and the capillary tip is somewhere in

5

between  and  , and well approximated by their arithmetic average15.

6 7



8

The inertial force  , that is generated due to the momentum of the dispersed phase inside the capillary and that pushes the droplet in the streamwise direction, is defined as  =    # 2

where  is the density of the dispersed phase,  is its average velocity inside the

9

capillary, and  is its volumetric flowrate.

10 11



The drag force, which also acts to detach the droplet from the capillary, is given by ! "#     = # 3 2

12

where  is the drag coefficient (dimensionless),  is the density of the continuous phase,

13

 is the relative velocity of the continuous phase with respect to the droplet, and "# is

14

the reference area. The drag coefficient  can be estimated as a function of the Reynolds

15

number Re as follows16: 9.06 !  = 0.292 (1 + - # 4 √+,

16

Here, the Reynolds number Re is defined as17 +, =

  / 0 −  2 # 5 3

17

where 0 is the diameter of the droplet and 3 is the dynamic viscosity of the

18

continuous phase. The term 0 −  reflects the fact that the capillary shields the

ACS Paragon Plus Environment

5

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 34

1

inner side of the droplet from the continuous phase, as shown in Figure 1. The relative

2

velocity of the continuous phase with respect to the droplet  is  − 0 , where  is

3

the average velocity of the continuous phase, 0 is the velocity of the droplet in the

4

! streamwise direction and it is given by  / 0 . The reference area "# is typically

5

defined as the area of an object perpendicular to the flow direction. In this work, it is

6

defined as17 "# =

! ! −  2 / 0

4

# 6

7

which accounts for the fact that the droplet is partially shielded from the continuous phase

8

flow by the capillary itself, as illustrated in Figure 1.

9 10

The buoyancy force  , that is generated due to the density difference between the two phases, is defined as  = 5 gΔ# 7

11

where 5 is the volume of the droplet, g is the acceleration of gravity, and Δ is the

12

difference between the densities of the two phases. Under the investigated conditions, the

13

buoyancy force is negligible in comparison with the other three forces, and therefore has

14

been neglected; this is often the case in liquid-liquid drop formation13,17–19.

15 16

By balancing the three horizontal forces, the size of droplets at the breaking point can be calculated:  =  +  0 ,  # 8

17

where only the drag force  increases with the droplet size 0 , whereas the interfacial

18

tension force  and the inertial force  are independent of 0 . Here, the average velocity of

19

the continuous phase  used in the drop size estimation is obtained from computational fluid

20

dynamics (CFD) simulations, as explained in Section 2.2.

ACS Paragon Plus Environment

6

Page 7 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

1

There are two main limitations to this model. First, the flow directions of the continuous and

2

dispersed phase are assumed to be parallel (referred in this work as the co-flow configuration). In

3

the case of cross-flow, the droplet breaks off from the external surface of the capillary rather than

4

from the tip, thus providing a larger contact area between the droplet and the capillary, and a

5

larger shielding effect from the continuous phase as compared to the co-flow case (Figure 3).

6

This increases the effect of the surface tension force and decreases the drag force, thus resulting

7

in a larger droplet size at break-off. Unfortunately, the actual break-off point is not known and

8

therefore the droplet size cannot be precisely estimated. Since developing a new model for

9

droplet formation lays beyond the scope of this paper, we have decided to use the simple model.

10

Second, the droplets form in the “dripping regime” where the dispersed phase drips from the

11

capillary tip in the form of droplets, thus being not applicable to the “jetting regime” drop

12

formation where the dispersed phase forms a liquid jet and brakes up into droplets due to

13

Rayleigh instabilities20.

14

To check the validity of the second assumption, a theoretical approach to characterize the flow

15

regime of droplet formation has been employed by comparing the kinetic energy per unit length

16

with the interfacial tension energy per unit length17. When the kinetic energy term is smaller than

17

the interfacial energy term, droplets form under the dripping regime. The kinetic energy of pipe

18

! flow per unit length is  !  /8 , where  is the density of the dispersed phase,  is the

19

average velocity of the dispersed phase inside the capillary, and  is the inner diameter of the

20

capillary. The interfacial tension energy per unit length is  +  /2. In this work, the

21

interfacial tension energy term is always larger than the kinetic energy term, thus implying that

22

the dripping regime governs the drop formation. It is, however, important to note that this

ACS Paragon Plus Environment

7

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 34

1

approach does not consider the energy transferred from the co-flowing continuous phase to the

2

droplet, which can influence the drop formation mechanism13.

3 4

5 6

Figure 1. Schematic drawing of the droplet formation at the capillary tip in the dripping regime.

7

The dispersed phase is injected through a capillary into a co-flowing immiscible outer fluid.

8

Three major forces acting on the droplets are displayed with solid arrows: interfacial tension

9

force  , inertial force  , and drag force  .

10

2.2 CFD Model

11

In order to get a better understanding of the flow conditions in the reactor and to obtain the local

12

fluid velocities needed for the estimation of the size of the generated droplets, a CFD model of

13

the crystallizer has been developed using the commercial CFD software Fluent 17.2. As a

14

significant vortex formation due to the intensity of mixing has been observed during the

15

experiments, a two-phase model was used to capture this effect. To solve the Navier-Stokes

16

equations, a standard Reynolds stress model (RSM) with linear pressure-strain has been used

17

together with the volume of fluid (VOF) multiphase model where the liquid phase has been

18

modeled as pure water and the gas phase as air. The simulations were run at steady-state using

19

the in-built multi reference frame (MRF) model for the description of the motion of the stirrer.

20

To reduce the computational burden, periodic boundary conditions have been used with only a

ACS Paragon Plus Environment

8

Page 9 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

1

quarter of the reactor (having a volume of 500 mL, see section 3.2) being sufficient to describe

2

the whole system. The grid has been refined until a mesh-independent solution was obtained

3

with the final grid consisting of 237,000 cells.

4 5

3. EXPERIMENTAL

6

Spherical agglomeration of benzoic acid crystals using toluene as a bridging liquid has been

7

investigated in this work, whose main goal is to study the effect of the initial droplet distribution

8

of the bridging liquid on the final particle size distribution (PSD) of the spherical agglomerates.

9

Therefore, the experiments have been designed (1) to study the effect of the initial droplet size

10

distribution on the final size distribution of the spherical agglomerates (series ;< ), (2) to establish

11

the effect of feeding position on the final size of agglomerates (series ;! ), and (3) to determine

12

which effect dominates in determining the final size of agglomerates – either the agglomeration

13

kinetics or the initial dispersion of the bridging liquid (series ;= ). Each experiment consisted of

14

two stages: first the unseeded generation of primary crystals (stage I), carried out under the same

15

conditions for all the experiments, followed by the agglomeration of the primary crystals into

16

spheres (stage II), performed using different operating conditions to explore the effects described

17

above.

18

3.1 Materials

19

In this work, benzoic acid (99.5% purity, Sigma-Aldrich, Switzerland) has been used as solute,

20

ethanol (absolute, Fischer Chemicals, Switzerland) as solvent, ultrapure deionized water (18.2

21

MΩ cm, Millipore, Switzerland) as antisolvent, and toluene (HPLC grade, Fischer Chemicals,

22

Switzerland) as bridging liquid. The physical properties of the liquids at room temperature are

23

summarised in Table 1.

ACS Paragon Plus Environment

9

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 34

1

Table 1. Physical properties of liquids used in this work at room temperature. Mixture

2

proportions are given in w/w Phase

Composition

 [kg/m= ]

3 [mPa s]

Continuous

Water/Ethanol (90/10 w/w)

963a

1.86a

Toluene

865

0.56

Dispersed 3

[mN/m] 20b

a21, b22

4

3.2 Experimental Setup

5

All experiments have been carried out as a semi-batch process using a 500 mL curved-bottom

6

cylindrical vessel with a diameter F=10 cm (GlasKeller, Germany) (see Figure 2). The

7

temperature inside the reactor was controlled using a Huber Ministat CC3 thermostat (Huber,

8

Germany). The solution was agitated using an overhead, 4-blade glass impeller with a diameter

9

of 75 mm and 40° inclined blades. The impeller was powered by a Eurostar Digital stirrer (IKA,

10

Germany) and positioned 25 mm above the reactor bottom with the shaft turning counter-

11

clockwise, thus resulting, in this case, in an upwards conveying flow direction. The lid of the

12

reactor was made of glass and had ports for a thermocouple, a liquid inlet used for antisolvent

13

crystallization, and a liquid inlet used for spherical agglomeration.

14

ACS Paragon Plus Environment

10

Page 11 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

1 2

Figure 2. Overview of the experimental setup

3

Stage I of the process – unseeded antisolvent crystallization – was performed by feeding an

4

undersaturated solution of benzoic acid in the solvent (ethanol) to the reactor filled with the

5

antisolvent (water). The stock solution was placed on a balance (Mettler Toledo, Switzerland)

6

and fed to the reactor using a LabAlliance Series II HPLC pump (Scientific Systems, USA), thus

7

enabling a precise dosage of the solution through the process. The feeding tube was located 1 cm

8

above the stirrer and 2 cm away from the shaft.

9

Stage II of the process – spherical agglomeration – was carried out by feeding the bridging

10

liquid (toluene) to the reactor containing crystals suspended in a saturated solution. The bridging

11

liquid was introduced using a PHD4400 syringe pump (Harvard Apparatus, USA). Due to the

12

high pressure drop in the system, a 20 mL stainless steel syringe (Harvard Apparatus, USA) has

13

been used to ensure precise dosage for all the process conditions investigated in this work. As

14

feeding tubes, fused silica capillaries (BGB Analytik AG, Switzerland) with four different sizes

15

have been selected with different inner and outer diameters,  and  , respectively (see Table

16

2). To protect them from mechanical damage and to allow for a precise positioning of the inlet

ACS Paragon Plus Environment

11

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 34

1

tip in the reactor, the feeding tubes have been guided within the reactor through a metal housing

2

(see Figure 3).

3

Table 2. Inner and outer diameter of capillaries used for feeding bridging liquid Capillary

H