On the Comparison between Presumed and Full PDF Methods for Turbulent Precipitation
|
|
- Lily McDaniel
- 5 years ago
- Views:
Transcription
1 Chemical and Biological Engineering Publications Chemical and Biological Engineering 2001 On the Comparison between Presumed and Full PDF Methods for Turbulent Precipitation Daniele L. Marchisio Politecnico di Torino Rodney O. Fox Iowa State University, Antonello A. Barresi Politecnico di Torino See next page for additional authors Follow this and additional works at: Part of the Process Control and Systems Commons The complete bibliographic information for this item can be found at cbe_pubs/100. For information on how to cite this item, please visit howtocite.html. This Article is brought to you for free and open access by the Chemical and Biological Engineering at Iowa State University Digital Repository. It has been accepted for inclusion in Chemical and Biological Engineering Publications by an authorized administrator of Iowa State University Digital Repository. For more information, please contact
2 On the Comparison between Presumed and Full PDF Methods for Turbulent Precipitation Abstract Turbulent precipitation is an important topic in chemical reaction engineering because of its numerous industrial applications. Several modeling approaches have been used in the past, but in recent years, computational fluid dynamics (CFD) coupled with micromixing models has been successfully applied to predict the influence of mixing on the crystal size distribution (CSD). The micromixing model is generally based on the presumed probability density function (PDF) approach, such as finite-mode PDF or beta PDF, and the aim of this work is to compare presumed PDF predictions and full PDF predictions with experimental data. The experimental data were obtained from a tubular reactor in which turbulent precipitation of barium sulfate is carried out from aqueous solutions of barium chloride and sodium sulfate. The implementation of the presumed PDF model was done using FLUENT user-defined subroutines, whereas the full PDF calculations were carried out with an in-house code based on Monte Carlo methods using the flow field prediction from FLUENT. Disciplines Chemical Engineering Process Control and Systems Comments This article is from Industrial & Engineering Chemistry Research 40 (2001): , doi: / ie Posted with permission. Authors Daniele L. Marchisio, Rodney O. Fox, Antonello A. Barresi, and Giancarlo Baldi This article is available at Iowa State University Digital Repository:
3 5132 Ind. Eng. Chem. Res. 2001, 40, On the Comparison between Presumed and Full PDF Methods for Turbulent Precipitation Daniele L. Marchisio,*, Rodney O. Fox, Antonello A. Barresi, and Giancarlo Baldi Dipartimento Scienza dei Materiali ed Ingegneria Chimica, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Torino, Italy, and Department of Chemical Engineering, Iowa State University, 2114 Sweeney Hall, Ames, Iowa Turbulent precipitation is an important topic in chemical reaction engineering because of its numerous industrial applications. Several modeling approaches have been used in the past, but in recent years, computational fluid dynamics (CFD) coupled with micromixing models has been successfully applied to predict the influence of mixing on the crystal size distribution (CSD). The micromixing model is generally based on the presumed probability density function (PDF) approach, such as finite-mode PDF or beta PDF, and the aim of this work is to compare presumed PDF predictions and full PDF predictions with experimental data. The experimental data were obtained from a tubular reactor in which turbulent precipitation of barium sulfate is carried out from aqueous solutions of barium chloride and sodium sulfate. The implementation of the presumed PDF model was done using FLUENT user-defined subroutines, whereas the full PDF calculations were carried out with an in-house code based on Monte Carlo methods using the flow field prediction from FLUENT. 1. Introduction The main advantage of using probability density function (PDF) methods lies in the fact that the chemical reaction term can be treated exactly, whereas most of the other methods require a closure assumption. PDF methods can be classified as full and presumed. Several differences between the two types of models exist, in terms of the discretization of the composition space, the spatial transport, and the numerical approach adopted. For full PDF methods, for example, spatial transport is treated as a random diffusion process, whereas presumed PDF methods are based on an assumption concerning scalar conditional mean concentrations. In the case of an inert scalar with uniform mean gradients, the full PDF method would predict that the joint scalar PDF is Gaussian (as indicated by direct numerical simulation results, as seen in the Experimental Setup and Numerical Simulation section), whereas the presumed PDF method cannot predict the correct limiting PDF. Moreover, unlike presumed PDF methods, full PDF methods do not require a priori knowledge of the shape of the joint composition PDF, but they are more CPU-intensive. In recent years, the presumed PDF method coupled with computational fluid dynamics (CFD) has been used to investigate turbulent precipitation in different devices. 1,2 The investigation of turbulent precipitation requires an adequate description of the population balance for the solid phase. This can be done by using discretization-based methods 3 or the moment method. 4 Unlike the former, the moment method requires only a small number of scalars, and recent results 5 show that is also applicable to complex problems such as sizedependent growth. * Corresponding author. marchis@athena.polito.it. Fax: Politecnico di Torino. Iowa State University. In general, using CFD coupled with a presumed PDF, good agreement with experimental data is found, but the uncertainty in kinetic expressions and the poor reproducibility of the experimental data highlight the need for other validation techniques. Thus, a comparison with full PDF predictions provides a useful means of examining the approximation of choosing the shape of the PDF a priori and its influence on the final crystal size distribution (CSD). Previous results in a Taylor- Couette reactor 6 and in stirred vessels 7-9 showed that mixing has a strong influence on the development of supersaturation and thus on particle nucleation and growth. In this work, barium sulfate precipitation from aqueous solutions in a tubular reactor with coaxial feeds is investigated. 2. Full and Presumed PDF Methods The full PDF description involves a set of notional particles that obey stochastic differential equations and move in the computational domain mimicking PDF transport in physical and composition space. A detailed description of PDF formulation and solution using Monte Carlo methods can be found in ref 10, and here, we limit the discussion to the transport equation of the composition PDF f φ t + u j f φ + [S x j ψ R (ψ)f φ ] ) R - [ u x i ψ f φ ] + [ D i ψ R 2 φ R ψ f φ ] (1) R where f φ (ψ;x,t) is the joint composition PDF of all scalars, u i is the fluid velocity, φ R represents the scalars under consideration (i.e., concentration, temperature, moments of the CSD), ψ R represents the corresponding variables in phase space, D R is the molecular diffusivity of φ R, S R (ψ) is the corresponding source term, and indicates the Reynolds average /ie CCC: $ American Chemical Society Published on Web 06/19/2001
4 Ind. Eng. Chem. Res., Vol. 40, No. 23, The first term on the right-hand side of eq 1 represents the convective transport due to the scalar conditioned velocity fluctuations. It can be closed by using a gradient-diffusion model that introduces a turbulent diffusion term u i ψ f φ )-Γ t f φ x i (2) where Γ t is the turbulent diffusivity defined as Γ t ) C µ k 2 Sc t ɛ Here, C µ is a constant equal to 0.09, Sc t is the turbulent Schmidt number equal to 0.7, k is the turbulent kinetic energy, and ɛ is the turbulent energy dissipation. The second term on the right-hand side of eq 1 represents the micromixing term that can be closed by using different approaches. In the case of nonreacting scalars, it has been shown 10 that the micromixing model must fulfill some constraints concerning the shape of the PDF: the mean values must be constant, the variance decay must be exponential, and the asymptotic shape of the PDF must be Gaussian. In this work, a simple interaction-with-the-mean (IEM) model was used. This model is also known as the linear mean square estimation (LMSE) model and is based on a linear relaxation of the scalars toward their mean values 11 (3) D R 2 φ R ψ ) 1 2τ φ ( φ R - ψ R ) (4) where φ R is the local scalar mean and τ φ is the local micromixing time. The IEM approach gives poor performance in the so-called partially segregated CSTR, where the combined effects of macro-, meso-, and micromixing are described by the micromixing model alone. 12 However, when the IEM model is properly used as a model for only micromixing and macro- and mesomixing due to spatial transport are handled separately, it yields acceptable predictions, once the micromixing time is chosen correctly. Several approaches have been proposed for evaluating this characteristic time. Baldyga 13 proposed a cascade model for the variance decay prediction, whereas more recently, Fox 14 proposed a spectral relaxation model in which the micromixing time was evaluated by solving the transport equations of the scalar variance and the scalar variance dissipation rate. However, these multiscale models are based on the idea of fully developed turbulence spectra (high Reynolds numbers), and to our knowledge, no fully validated models exist for lower Reynolds numbers. As explained below, this work was carried out at moderate Reynolds numbers, relatively close to the transition region. Under these conditions, the use of a simple large-scale-motion-dominated model, 10 seems to be more reasonable 1 ) C φ ɛ τ φ k (5) where C φ is a constant that can be adjusted to take into account the above-mentioned moderate Reynolds number effects. In the presumed PDF approach, the functional form of the PDF is assumed a priori. In this work, we make use of the finite-mode-pdf (FM-PDF) model. Using this model, the PDF is represented by a finite set of delta functions 15 where p n (x,t) is the probability of mode n, φ R n (x,t) is the value of scalar R corresponding to mode n, N e is the total number of modes, and m is the total number of scalars. The total number of modes affects the model s ability to approximate the real PDF. When this number is increased, the description becomes more accurate, but the computational cost increases. The advantage of this model is its ability to accurately describe nonpremixed reacting flows with a small number of modes. Previous results 16 showed that N e ) 3 is sufficient to obtain good accuracy. A detailed description of this approach can be found in ref 2. Without loss of generality, we take the first scalar to be the mixture fraction. The evolution of the mode probabilities is defined once a transport equation for the probability vector p n (x,t) has been provided. To ensure the conservation of probability (probabilities sum to unity) and the conservation of mass, the transport equations can be written in terms of p n (x,t) and the weighted mixture fraction in mode 3 ( s ξ 3 ) ξ 3 p 3 )as and N e m f φ (ψ;x,t) p n (x,t) δ(ψ R - φ R n (x,t)) (6) n)1 R)1 p 1 t + ( u x i p 1 ) ) i x i( Γ p 1 t x i) + γ s p 3 - γp 1 (1 - p 1 ) (7) p 2 t + ( u x i p 2 ) ) i x i( Γ p 2 t x i) + γ s p 3 - γp 2 (1 - p 2 ) (8) s ξ 3 t + ( u x i s ξ 3 ) ) i x i ( Γ s ξ 3 t x i ) - γ s p 3 ( ξ 1 + ξ 2 ) + γp 1 (1 - p 1 ) ξ 1 + γp 2 (1 - p 2 ) ξ 2 (9) In these equations, p 3 ) 1 - p 1 - p 2 ; Γ t is the turbulent diffusivity; γ s and γ are the spurious dissipation rate and the micromixing rate, respectively; and ξ 1 and ξ 2 are the mixture fraction values in modes 1 and 2, respectively, which are constant and equal to their inlet values (i.e., ξ 1 ) 1 and ξ 2 ) 0). Note that the terms involving γ correspond to molecular mixing and, thus, generate the scalar dissipation rate in the transport equation for the mixture fraction variance, whereas the terms involving γ s are required to eliminate spurious scalar dissipation resulting from the finite-mode representation of turbulent diffusion. When the full PDF method is used, this spurious dissipation term is eliminated at the price of requiring a large ensemble of notional particles, as will be explained below. In the FM-PDF model, this term can be eliminated by setting γ s as follows γ s ) 2Γ t ξ 3 ξ 3 (10) 1-2 ξ 3 (1 - ξ 3 ) x i x i The term γ can be calculated by forcing the mixture
5 5134 Ind. Eng. Chem. Res., Vol. 40, No. 23, 2001 Table 1. Source Terms for Each Scalar mixture fraction (ξ) 0 jth moment of the CSD (m j) (0) j B(c A,c B) + jg(c A,c B)m j-1 reaction progress variable (Y) F3k vm 2G(c A,c B)/[ξ stc AoM] fraction variance to follow the exact transport equation. Then, according to eq 5 γ ) C φ ɛ k ξ 2 [p 1 (1 - p 1 )(1 - ξ 3 ) 2 + p 2 (1 - p 2 ) ξ 3 2 ] (11) where ξ 2 is the Reynolds-averaged mixture fraction variance and C φ was kept equal to 1. Concerning the choice of C φ, the effect of this constant on the CSD was discussed in ref 16, and this value was chosen with particular regard to the problem of dealing with nonfully developed turbulence. Moreover, it should be highlighted that, because the purpose of this work is the comparison of the two approaches, it is important to use the same value of C φ in both models. Concerning the other scalars, generally for mode n, it is possible to define a local concentration denoted by φ R n, whereas the weighted concentration s R n is defined by s R n ) φ R n p n (12) The transport equation for the weighted concentration of scalar R for mode 3 is s R 3 + ( u t x i s R 3 ) ) i x i( Γ s R 3 t x i ) - γ t p 3 ( φ R 1 + φ R 2 ) + γ s p 1 (1 - p 1 ) φ R 1 + γ s p 2 (1 - p 2 ) φ R 2 + p 3 S R ( φ 3 ) (13) where φ R 1 and φ R 2 are the local concentrations in modes 1 and 2, respectively, and S R ( φ 3 ) is the source term for scalar R. 3. Moment Method Approach As highlighted in the Introduction, the population balance must be solved in order to describe particle nucleation and growth. Using the moment method, the system can be solved in terms of the moments of the CSD. The first five moments are of particular interest, as they are related to the total particle number density (N t ) m 0 ), the total particle area (A t ) k a m 2 ), and the total solid volume (V t ) k v m 3 ) by shape factors (k a, k v ) that depend on particle morphology. Using this approach, the mean crystal size (L 43 ) and the solid concentration (c C ) can be written as L 43 ) m 4, c m C ) Fk v m 3 3 M (14) where F is the crystal density, k v is the volume shape factor, and M is the molecular weight of the crystal. Using this approach, a set of five transport equations must be solved. Depending on the approach used (i.e., full or presumed PDF), the implementation is slightly different. In Table 1, the source term that should be used for each scalar is reported. These terms must be included in eq 1 or eq 13. Note that, for the FM-PDF model, the transport equation for the mixture fraction has already been written (see eq 9). In Table 1, the nucleation rate [B(c A,c B )] and the growth rate [G(c A,c B )] appear, and they can be calculated once the reactant concentrations are known using the equations where c A c Ao ) ξ - ξ s Y, c B c Bo ) (1 - ξ) - (1 - ξ s )Y (15) ξ s ) c Bo c Ao + c Bo (16) and Y is the reaction progress variable whose source term is defined in Table 1. Note that Y ) 0 when ξ ) 0 and ξ ) 1. The nucleation rate is taken from literature 7 { c [1/(m 3 s)] for c e 10 mol/m 3 (heterogeneous) B(c A,c B ) ) c 15 [1/(m 3 s)] for c > 10 mol/m 3 (homogeneous) (17) where c ) c A c B - k s and k s is the solubility product of barium sulfate (at room temperature, k s ) mol 2 /m 6 ). For the growth rate, a two-step diffusion-adsorption model 17 was used, which gives the following expression G(c A,c B ) ) c s 2 ) k d (c A - c As ) ) k d (c B - c Bs ) (m/s) (18) where c s ) c As c Bs - k s and k d is the mass transfer coefficient. According to ref 18, this coefficient is sizedependent but remains nearly constant for particles smaller than 10 µm, notwithstanding the difference in solute. As a consequence, k d was fixed to fall in the range between 10-7 and 10-8 (m/s)(m 3 /mol). 7 More recent results 19 have shown that the limit of the Sherwood number Sh ) 2 is also valid for microparticles. Accordingly, the mass transfer coefficient can be calculated from k d ) ShDM L 43 F (19) For average particle sizes (between reactor inlet and outlet) of about 0.1 and 1.0 µm, k d is equal to 10-6 and 10-7 (m/s)(m 3 /mol), respectively. Concerning the shape factors, for the operating conditions under consideration, the crystals are flat, and the mean crystal size measured by a laser particle sizer is very close to the width of the crystal. In previous work, 20,21 shape factors using the width as the characteristic dimension were determined, and their effects on the CSD were studied. Nevertheless, the population balance works in terms of a unique characteristic dimension (i.e., L 43 ), and thus, equivalent shape factors for a system with equidimensional growth along each particle axis must be used. 22 Using this approach, it was found that k v ) 5.0 and k a ) Because the comparison with experimental data must be made in consistent units, the characteristic dimension of the population balance has to be converted to the dimension measured by the particle sizer, that is, the characteristic
6 Ind. Eng. Chem. Res., Vol. 40, No. 23, dimension used in the former definition of shape factors, by using d 43 ) φ c L 43 (20) where φ c is a function of shape factors and was found to be equal to Experimental Setup and Numerical Simulation Simulations and experiments were carried out in a single-jet tubular reactor consisting of a small tube centered on the reactor axis with a 1-mm inner diameter, a 1.5-mm outer diameter, and a m length. The tubular reactor has an inner diameter of 1 cm and a length of 2.1 m. The velocity ratio (VR) between the jet and the coaxial flow was kept equal to 1, and the macroscopic Reynolds number was near Reactants were fed at several concentrations, with the reactant concentration ratio (R )c Bo /c Ao ) varied in the range for c Ao ) mol/m 3. The index A denotes that the reactant fed in the inner tube for which ξ ) 1. Likewise, B denotes the reactant fed in the annular region for which ξ ) 0. A detailed description of the experimental setup can be found in ref 23, and here, we limit our discussion to the final results. Barium chloride and sodium sulfate were alternatively fed in the inner tube and in the annular region, keeping the concentrations constant. Hereinafter, these two configurations are denoted by AB and BA, respectively. The flow field was modeled by using the k - ɛ model 24 in FLUENT (version 5.2), and standard wall functions were used for the wall boundary layer. The grid for the flow field simulation had 131 nodes in the axial direction and 35 nodes in the radial direction. This grid included the annular region before the injection zone (x o ) 168 mm with 20 nodes) so that the velocity profile in the inlet zone could be predicted correctly. This region was not considered in the full PDF simulation because no reaction occurs. Thus, the final grid for the full PDF simulation had 111 nodes in the axial direction and 35 nodes in the radial direction. The full PDF code was initialized by reading the axial and radial coordinates used in FLUENT and the steadystate solution of the flow field, including the mean velocities (u x and u y ), k, ɛ, the mean pressure, and the fluid density. A number of particles equal to n p n cell, where n p ) 100, were positioned in the computational domain, and each particle s flow properties were fixed by linear interpolation depending on its location. Particles were initialized at the inlet and were moved according to the velocity vector and the turbulent diffusivity assigned to them. Then, particles interacted with each other as a result of micromixing (see eq 4), where the local mean values were derived by taking the ensemble average of the scalar in each cell. (For details, see ref 25.) The changes in scalar composition due to chemical reaction, nucleation, and crystal growth were calculated for each particle by direct integration using a routine based on a fourth-order Runge-Kutta method. Concerning the FM-PDF simulations, transport equations for the mode probabilities p 1 and p 2, the weighted mixture fraction s ξ 3, the reaction progress variable s Y 3, and the five moments of the CSD s mj 3 (j ) 0,..., 4) were introduced in the CFD code as user-defined scalars. Once the flow field had converged, this set of equations was solved. Appropriate under-relaxation Figure 1. Full PDF prediction of concentration of reactant A (c Ao ) mol/m 3, R)1). Figure 2. Full PDF prediction of concentration of reactant B (c Ao ) mol/m 3, R)1). factors were used depending on the stiffness of the problem (i.e., initial reactant concentration). Note that, unlike the full PDF model (which is based on Monte Carlo methods), the FM-PDF model was solved in the CFD code itself based on a finite-volume discretization. With this approach, the grid had to be fine enough to avoid numerical diffusion. To ensure that the solution was grid-independent, different grids and different discretization schemes were tested, but for these operating conditions, no substantial differences were detected. 5. Results and Discussion In Figures 1 and 2, the reactant concentrations (c A, c B ) obtained with the full PDF code are reported. As one can see, reactant A is fed in the small tube, whereas reactant B is fed in the annular region. As a result of turbulent diffusion, the reactants mix together and react. Reactant B is in great excess, and thus, its concentration remains almost constant, whereas reactant A disappears, first because of mixing and later because of chemical reaction. The plot of the nucleation rate (see Figure 3) shows a peak near the injection point, which is due to very high supersaturation in this region.
7 5136 Ind. Eng. Chem. Res., Vol. 40, No. 23, 2001 Figure 3. Full PDF prediction of nucleation rate (c Ao ) mol/m 3, R)1). Figure 5. Comparison between the two PDF approaches for the mean mixture fraction (ξ) and the intensity of segregation I s at (x - x o)/d ) 0.0. Figure 4. Full PDF prediction of mean crystal size (c Ao ) mol/m 3, R)1). The nucleation rate reaches very high values in this region and then decreases rapidly. The plot of the mean crystal size (L 43 ) is shown in Figure 4. As can be seen, the mean crystal size increases until the final value is reached. It should be highlighted that the wavelike fluctuations in the scalar fields are caused by statistical errors, which can be improved by using more notional particles per cell or by averaging over several time steps. To compare results provided by different models that are based on different methods, it is useful to compare the mixture fractions and the intensities of the segregation profiles in different sections of the reactor. The mixture fraction is a nonreacting scalar, and thus, its mean value ( ξ ) and its variance ( ξ 2 ) are affected only by convection and turbulent diffusion. The intensity of segregation I s is defined by I s ) ξ 2 ξ (1 - ξ ) (21) and is an index of the degree of mixing; in fact, it is equal to 1 when the system is perfectly segregated and 0 when the system is perfectly micromixed. In Figures 5-7, the mean mixture fraction and the intensity of segregation versus the radial direction at several axial Figure 6. Comparison between the two PDF approaches for the mean mixture fraction (ξ) and the intensity of segregation I s at (x - x o)/d ) 1.0. positions [(x - x o )/D] are reported. A comparison is made between the time-averaged full PDF model predictions and the FM-PDF model predictions, and as can be seen, the predictions are quite similar. Concerning the reacting scalars, contour plots of the mean reactant concentrations, nucleation rate, and mean crystal size look similar, but differences were detected. For example, the FM-PDF model seems to underpredict the nucleation rate (as compared to the full PDF results) and, thus, the final particle number density. The value of the total number density is extremely sensitive to numerics because of the strong nonlinearity of the nucleation rate. However, as mentioned before, by testing different grids, the solution was shown to be grid-independent. Likewise, by using different discretization schemes, the approach used was shown to be adequate to handle the stiffness of the problem. In Figure 8, the total particle number density along the reactor axis is reported for c Ao ) mol/ m 3 and R)1. The peak of m 0 after the injection zone predicted by all models is initially very high but decreases sharply. The comparison of the mean crystal sizes (see Figure 9) highlights the fact that, without any closure, this property is detectably underestimated.
8 Ind. Eng. Chem. Res., Vol. 40, No. 23, Figure 7. Comparison between the two PDF approaches for the mean mixture fraction (ξ) and the intensity of segregation I s at (x - x o)/d ) Figure 9. Comparison between model predictions for mean crystal size along the center line (c Ao ) mol/m 3, R)1). Figure 8. Comparison between model predictions for total number density along the center line (c Ao ) mol/m 3, R) 1). Notice that the full PDF predictions presented in the last two figures are time-averaged values. However, for other values of R, agreement between the FM-PDF and full PDF predictions was better. Radial-averaged results at the end of the reactor and a comparison with experimental data are presented below. Concerning the comparison with experimental data, the choice of the mass transfer coefficient (k d ) is crucial. In Figure 10, the FM-PDF predictions with different k d values are compared with experimental data for the two configurations (AB and BA). The values of k d reported in the legend can be obtained by using eq 19 and assuming average mean crystal sizes (i.e., L 43 ) of about 1, 0.1, and 0.01 µm, respectively. As one can see, the two configurations (AB and BA) result in different experimental CSDs. This has been explained by some authors in terms of activity coefficients 26 that should be included in the growth and nucleation rates or by invoking different growth mechanisms, as suggested in ref 27. It should be also highlighted that the models presented here are perfectly symmetric and, thus, are unable to recognize the difference between excess Figure 10. Comparison between FM-PDF predictions and the experimental mean crystal size (c Ao ) mol/m 3 ). barium versus excess sulfate. Although the results show good agreement, they highlight the need for using a sizedependent mass transfer coefficient. In fact, the best agreement is obtained by using a different k d value for each point depending on the corresponding mean crystal size. In Figure 11, model predictions for the three models are compared with experimental data. As can be seen, the overall agreement is quite good, but in general, for these operating conditions, the effect of micromixing is small. At higher concentrations, the effect of micromixing should be greater, but aggregation begins to play an important role, at which point all of the models definitively underpredict the mean crystal size. However, at higher concentrations, the use of a micromixing model should be even more important, as it is needed for correctly predicting the local interactions between microcrystals close to the injection zone because of local segregation. 6. Conclusion Turbulent precipitation of barium sulfate in a tubular reactor was studied. Reactants were fed in a coaxial single-jet configuration. The effect of the reactant
9 5138 Ind. Eng. Chem. Res., Vol. 40, No. 23, 2001 Figure 11. Comparison between FM-PDF predictions and the experimental mean crystal size [c Ao ) mol/m 3, k d ) 10-6 (m/s)(m 3 /mol)]. concentration ratio on the CSD was studied, and the predictions of the FM-PDF and full PDF approaches were compared with experimental data. Results show that the population balance treatment needs to be improved; in fact, the importance of using a sizedependent mass transfer coefficient was highlighted (i.e., size-dependent growth rate.) This can be handled by using the quadrature of moments method as explained in ref 5. Moreover, at higher concentrations both models fail in predicting the final mean crystal size because aggregation is neglected. However, the ability of this approach to describe turbulent precipitation appears to be extremely promising, and the use of a complete description of solid evolution (including aggregation) linked with the PDF description represents the natural development of this work. Acknowledgment This research was partially supported by an Italian National Research Project (MURST 40%sMultiphase reactors: Hydrodynamics analysis and solid-liquid analysis). Notation B ) nucleation rate, number/(m 3 s) c A ) barium chloride concentration, mol/m 3 c B ) sodium sulfate concentration, mol/m 3 c C ) solid concentration, kg/m 3 C µ ) turbulent constant C φ ) micromixing constant d 43 ) mean crystal size, m f φ (ψ;x,t) ) joint probability density function G ) crystal growth rate, m/s k ) turbulent kinetic energy, m 2 /s 2 k a ) surface shape factor k d ) mass transfer coefficient, (m/s)(m 3 /mol) k s ) barium sulfate solubility product, mol 2 /m 6 k v ) volume shape factor L 43 ) ratio of the fourth and third moments of the CSD, m m j ) jth moment of the CSD, m j-3 M ) barium sulfate molecular weight, kg/mol n cell ) number of cells in the computational domain n p ) number of notational particles per cell N e ) number of modes for the description of the FM-PDF model p n ) probability of mode n s ξ n ) weighted mixture fraction in mode n s R n ) weighted concentration of scalar R in mode n Sc t ) turbulent Schimdt number u i ) Reynolds averaged velocity in the ith direction, m/s Y ) reaction progress variable Greek Letters R)reactant concentration ratio Γ t ) turbulent diffusivity, m 2 /s γ ) micomixing rate, 1/s γ s ) spurious dissipation rate, 1/s ɛ ) turbulent dissipation rate, m 2 s 3 ɛ ξ ) scalar dissipation rate, 1/s ν ) kinematic viscosity, m 2 /s ξ n ) local mixture fraction in mode n ξ 2 ) mixture fraction variance ξ s ) stoichiometric mixture fraction F)crystal density, kg/m 3 τ φ ) micromixing time, s φ R n ) local concentration of scalar R in mode n φ c ) length shape factor Subscript o ) value in the inlet stream Operators ) expected value ) conditioned expected value Literature Cited (1) Baldyga, J.; Orciuch, W. Closure method for precipitation in inhomogeneous turbulence. Proceedings of the 14th Symposium on Industrial Crystallization, Sept 12-16, 1999; Institution of Chemical Engineers: Rugby, U.K., 1999; Paper 86, p (2) Marchisio, D. L.; Fox, R. O.; Barresi, A. A. Simulation of turbulent precipitation in a Taylor-Couette reactor using CFD. AIChE J. 2001, 47, 664. (3) Vanni, M. Approximate population balance equations for aggregation-breakage processes. J. Colloid Interface Sci. 2000, 221, 143. (4) Rivera, T.; Larson, M. A. A model for the precipitation of pentaerythritol tetranitrate (PETN). Ind. Eng. Chem. Process Des. Dev. 1978, 17, 183. (5) McGraw, R. Description of aerosol dynamics by quadrature methods of moments. Aerosol Sci. Technol. 1997, 27, 255. (6) Barresi, A. A.; Marchisio, D.; Baldi, G. On the role of microand mesomixing in a continuous Couette-type precipitator. Chem. Eng. Sci. 1999, 54, (7) Baldyga, J.; Podgorska, W.; Pohorecky, H. Mixing precipitation model with application to double feed semibatch precipitation. Chem. Eng. Sci. 1995, 50, (8) Fitchett, D. E.; Tarbell, J. M. Effect of mixing on the precipitation of barium sulphate in an MSMPR reactor. AIChE J. 1990, 36, 511. (9) Phillips, R.; Rohani, S.; Baldyga, J. Micromixing in a singlefeed semi-batch precipitation process. AIChE J. 1999, 45, 82. (10) Pope, S. Pdf methods for turbulent reacting flows. Prog. Energy Combust. Sci. 1985, 11, 119. (11) Dopazo, C. Probability density function approach for turbulent axisymmetric heated jet. Centerline evolution. Phys. Fluids 1975, 18, 397. (12) Baldyga, J.; Bourne, J. R. Comparison of the engulfment and the interaction-with-the-mean micromixing models. Chem. Eng. J. 1990, (13) Baldyga, J. Turbulent mixer model with application to homogeneous, instantaneous chemical reactions. Chem. Eng. Sci. 1989, 44, (14) Fox, R. O. The Lagrangian spectral relaxation model of the scalar dissipation in homogeneous turbulence. Phys. Fluids 1997, 9, 2364.
10 Ind. Eng. Chem. Res., Vol. 40, No. 23, (15) Fox, R. O. On the relationship between Lagrangian micromixing models and computational fluid dynamics. Chem. Eng. Proc. 1998, 37, 521. (16) Piton, D.; Fox, R. O.; Marcant, B. Simulation of fine particle formation by precipitation using computational fluid dynamics. Can. J. Chem. Eng , 983. (17) Nielsen, A. E. Electrolyte crystal growth mechanisms. J. Crystal Growth 1984, 67, 289. (18) Nagata, S. Mixing Principles and Applications; Kodansha: Tokyo, Japan, (19) Armenante, P. M.; Kirwan, D. J. Mass transfer to microparticles in agitated system. Chem. Eng. Sci. 1989, 44, (20) Pagliolico, S.; Marchisio, D.; Barresi, A. A. Influence of operating conditions on BaSO 4 crystal size and morphology in a continuous Couette-type precipitator. J. Therm. Anal. Calorim. 1999, 56, (21) Marchisio, D. L.; Barresi, A. A.; Fox, R. O.; Baldi, G. Comparison of different modelling approaches to turbulent precipitation. In Proceedings of the 10th European Conference on Mixing, Delft, The Netherlands, July 2-5, 2000; van den Akker, H. E. A., Derksen, J. J., Eds.; Elsevier Science: New York, 2000; p 77. (22) Randolph, A. D.; Larson, M. A. Theory of Particulate Processes, 2nd ed.; Academic Press: New York, (23) Marchisio, D. L.; Fox, R. O.; Barresi, A. A.; Garbero, M.; Baldi, G. On the simulation of turbulent precipitation in a tubular reactor via Computational Fluid Dynamics. Proceedings of the 4th International Symposium on Mixing in Industrial Processes (ISMIP4), CD-ROM ed., Toulouse, France, May 14-16, 2001, paper 55; submitted to Chem. Eng. Res. Dev. (24) Fluent Inc. FLUENT 5 User s Guide; Fluent Inc.: Lebanon, NH, (25) Tsai, K.; Fox, R. O. PDF simulation of turbulent seriesparallel reaction in an axisymmetric reactor. Chem. Eng. Sci. 1994, 49, (26) Wei, H.; Garside, J. Application of CFD modelling of precipitation system. Chem. Eng. Res. Dev. 1997, 75A, 219. (27) Wong, D. C. Y.; Jaworski, J.; Nienow, A. W. Effect of ion excess on particle size and morphology during barium sulphate precipitation: An experimental study. Chem. Eng. Sci. 2001, 56, 1. Received for review December 1, 2000 Revised manuscript received April 4, 2001 Accepted April 26, 2001 IE
STUDY OF AGGREGATION IN BARIUM SULPHATE PRECIPITATION
15 th International Symposium on Industrial Crystallization Sorrento, Italy, 15-18 September 2002 paper 171 STUDY OF AGGREGATION IN BARIUM SULPHATE PRECIPITATION D. L. Marchisio 1, A. A. Barresi, M. Garbero,
More informationEXPERIMENTAL INVESTIGATION AND CFD MODELING OF MICROMIXING OF A SINGLE-FEED SEMI-BATCH PRECIPITATION PROCESS IN A LIQUID-LIQUID STIRRED REACTOR
14 th European onference on Mixing Warszawa, 10-13 September 2012 EXPERIMENTAL INVESTIGATION AND FD MODELING OF MIROMIXING OF A SINGLE-FEED SEMI-BATH PREIPITATION PROESS IN A LIQUID-LIQUID STIRRED REATOR
More informationprecipitation in a Confined Impinging Jets Reactor by means of
Modeling of nanoparticles precipitation in a Confined Impinging Jets Reactor by means of Computational Fluid Dynamics E. Gavi, D.L. Marchisio, A.A. Barresi Politecnico di Torino, Department of Material
More informationConditional statistics for passive-scalar mixing in a confined rectangular turbulent jet
Chemical and Biological Engineering Publications Chemical and Biological Engineering 2007 Conditional statistics for passive-scalar mixing in a confined rectangular turbulent jet H. Feng Iowa State University
More informationModeling Industrial Crystallizers
Modeling Industrial Crystallizers Kumar Dhanasekharan, Ph.D. Fluent Inc., 10 Cavendish Court, Lebanon, NH 03766 Phone: (603) 643-2600 x323; Email: kd@fluent.com One of the main challenges in industrial
More informationSimulation of mixing and precipitation of nanoparticles for pharmaceutical applications with CFD and MD
14th European Conference on Mixing Warszawa, 10-13 September 2012 Simulation of mixing and precipitation of nanoparticles for pharmaceutical applications with CFD and MD N. Di Pasquale 1, P. Carbone 2,
More informationEffect of Feed-Stream Configuration on Gas-Phase Chlorination Reactor Performance
Chemical and Biological Engineering Publications Chemical and Biological Engineering 2003 Effect of Feed-Stream Configuration on Gas-Phase Chlorination Reactor Performance Venkatramanan Raman Iowa State
More informationS. T. Smith Iowa State University. Rodney O. Fox Iowa State University,
Chemical and Biological Engineering Publications Chemical and Biological Engineering 2007 A term-by-term direct numerical simulation validation study of the multi-environment conditional probability-density-function
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationCFD ANALYSIS OF TURBULENCE EFFECT ON REACTION IN STIRRED TANK REACTORS
CFD ANALYSIS OF TURBULENCE EFFECT ON REACTION IN STIRRED TANK REACTORS Udaya Bhaskar Reddy R*, Gopalakrishnan S, Ramasamy E Department of Chemical Engineering, Coimbatore Institute of Technology, Coimbatore-
More informationHybrid CFD-Multizonal Modelling of Polymorphs and Agglomeration Phenomena in Crystallisation Processes
Hybrid CFD-Multizonal Modelling of Polymorphs and Agglomeration Phenomena in Crystallisation Processes Sean K. Bermingham 1, Stefania Trivellato 1,2, Maddalena Vernier 1,2, Fabrizio Bezzo 2,* 1 Process
More informationCFD Simulation of Binary Fluidized Mixtures: Effects of Restitution Coefficient and Spatial Discretization Methods
Engineering Conferences International ECI Digital Archives The 14th International Conference on Fluidization From Fundamentals to Products Refereed Proceedings 2013 CFD Simulation of Binary Fluidized Mixtures:
More informationAggregation and Breakage of Nanoparticle Dispersions in Heterogeneous Turbulent Flows
ggregation and reakage of Nanoparticle Dispersions in Heterogeneous Turbulent Flows M. Soos,, D. Marchisio*, J. Sefcik, M. Morbidelli Swiss Federal Institute of Technology Zürich, Z CH *Politecnico di
More informationSoot formation in turbulent non premixed flames
Soot formation in turbulent non premixed flames A. Cuoci 1, A. Frassoldati 1, D. Patriarca 1, T. Faravelli 1, E. Ranzi 1, H. Bockhorn 2 1 Dipartimento di Chimica, Materiali e Ing. Chimica, Politecnico
More informationAlgebraic Turbulence-Chemistry Interaction Model
50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 09-12 January 2012, Nashville, Tennessee AIAA 2012-0183 50th AIAA Aerospace Sciences Meeting, January 9-12,
More informationPolitecnico di Torino. Porto Institutional Repository
Politecnico di Torino Porto Institutional Repository [Proceeding] Simulation of macromolecule self-assembly in solution: multiscale approach a Original Citation: Lavino, A; Di Pasquale, N; Carbone, P;
More informationNumerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure
21 st ICDERS July 23-27, 27 Poitiers, France Numerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure Ahmad S. El Sayed, Cécile B. Devaud Department of Mechanical
More informationOn the Comparison between Population Balance Models for CFD Simulation of Bubble Columns
Chemical and Biological Engineering Publications Chemical and Biological Engineering 2005 On the Comparison between Population Balance Models for CFD Simulation of Bubble Columns Jay Sanyal Fluent Inc.
More informationComputational Fluid Dynamics Simulation of Chemical Reactors: Application of in Situ Adaptive Tabulation to Methane Thermochlorination Chemistry
Chemical and Biological Engineering Publications Chemical and Biological Engineering 1999 Computational Fluid Dynamics Simulation of Chemical Reactors: Application of in Situ Adaptive Tabulation to Methane
More informationTransported PDF Calculations of Combustion in Compression- Ignition Engines
International Multidimensional Engine Modeling User s Group Meeting at the SAE Congress Detroit, MI 15 April 2013 Transported PDF Calculations of Combustion in Compression- Ignition Engines V. Raj Mohan
More informationTurbulent Mixing in a Tubular Reactor: Assessment of an FDF/LES Approach
Turbulent Mixing in a Tubular Reactor: Assessment of an FDF/LES Approach E. van Vliet, J. J. Derksen and H. E. A. van den Akker Kramers Laboratorium voor Fysische Technologie, Faculty of Applied Sciences,
More informationCOMPARISON OF ANALYTICAL SOLUTIONS FOR CMSMPR CRYSTALLIZER WITH QMOM POPULATION BALANCE MODELING IN FLUENT
CHINA PARTICUOLOGY Vol. 3, No. 4, 13-18, 5 COMPARISON OF ANALYTICAL SOLUTIONS FOR CMSMPR CRYSTALLIZER WITH QMOM POPULATION BALANCE MODELING IN FLUENT Bin Wan 1, Terry A. Ring 1, *, Kumar M. Dhanasekharan
More informationOverview of Reacting Flow
Overview of Reacting Flow Outline Various Applications Overview of available reacting flow models Latest additions Example Cases Summary Reacting Flows Applications in STAR-CCM+ Chemical Process Industry
More informationCFD and Kinetic Analysis of Bluff Body Stabilized Flame
CFD and Kinetic Analysis of Bluff Body Stabilized ame A. Dicorato, E. Covelli, A. Frassoldati, T. Faravelli, E. Ranzi Dipartimento di Chimica, Materiali e Ingegneria Chimica, Politecnico di Milano, ITALY
More informationLiquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses
Liquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses by W. A. Sirignano Mechanical and Aerospace Engineering University of California, Irvine Collaborators:
More informationflame synthesis of inorganic nanoparticles
flame synthesis of inorganic nanoparticles Shraddha Shekar and Jethro Akroyd Markus KRAFT Karlsruhe 12 July 2010 Part 1 Introduction Precursor (TES) Mesoporous silica nanoparticles Aim: To answer the following
More informationEVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS
EVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS A Aroussi, S Kucukgokoglan, S.J.Pickering, M.Menacer School of Mechanical, Materials, Manufacturing Engineering and
More informationAnovel type of continuous mixing system is shown in
Mixing in an Agitated Tubular Reactor J. J. Derksen* School of Engineering, University of Aberdeen, Aberdeen, UK We analyze, through numerical simulations, the single-phase liquid flow and associated passive
More informationCRYSTALLIZATION AND PRECIPITATION ENGINEERING
CRYSTALLIZATION AND PRECIPITATION ENGINEERING Alan Jones, Rudi Zauner and Stelios Rigopoulos Department of Chemical Engineering University College London, UK www.chemeng.ucl.ac.uk Acknowledgements to:
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationINVESTIGATIONS OF MASS TRANSFER AND MICROMIXING EFFECTS IN TWO-PHASE LIQUID-LIQUID SYSTEMS WITH CHEMICAL REACTION
14 th European Conference on Mixing Warszawa, 10-13 September 20 INVESTIGATIONS OF MASS TRANSFER AND MICROMIXING EFFECTS IN TWO-PHASE LIQUID-LIQUID SYSTEMS WITH CHEMICAL REACTION M. Jasińska a, J. Bałdyga
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationTHE INFLUENCE OF MIXING ON FAST PRECIPITATION PROCESSES A COUPLED CFD- PBE APPROACH USING THE DIRECT QUADRATURE METHOD OF MOMENTS (DQMOM)
Eleventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 7-9 December 2015 THE INFLUENCE OF MIXING ON FAST PRECIPITATION PROCESSES A COUPLED CFD- PBE
More informationStudy on residence time distribution of CSTR using CFD
Indian Journal of Chemical Technology Vol. 3, March 16, pp. 114-1 Study on residence time distribution of CSTR using CFD Akhilesh Khapre*, Divya Rajavathsavai & Basudeb Munshi Department of Chemical Engineering,
More informationDEVELOPMENT OF A NUMERICAL APPROACH FOR SIMULATION OF SAND BLOWING AND CORE FORMATION
TMS (The Minerals, Metals & Materials Society), DEVELOPMENT OF A NUMERICAL APPROACH FOR SIMULATION OF SAND BLOWING AND CORE FORMATION G.F. Yao, C. W. Hirt, and
More informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationLecture (9) Reactor Sizing. Figure (1). Information needed to predict what a reactor can do.
Lecture (9) Reactor Sizing 1.Introduction Chemical kinetics is the study of chemical reaction rates and reaction mechanisms. The study of chemical reaction engineering (CRE) combines the study of chemical
More informationAerosol nucleation and growth in a turbulent jet using the Stochastic Fields method
Loughborough University Institutional Repository Aerosol nucleation and growth in a turbulent jet using the Stochastic Fields method This item was submitted to Loughborough University's Institutional Repository
More informationCFD Modelling of Turbulent Mass Transfer in a Mixing Channel
CFD Modelling of Turbulent Mass Transfer in a Mixing Channel Lene K. Hjertager Osenbroch, Bjørn H. Hjertager and Tron Solberg Aalborg University Esbjerg Esbjerg, Denmark Homepage: hugin.aue.auc.dk Prepared
More informationNUMERICAL SIMULATION OF THREE DIMENSIONAL GAS-PARTICLE FLOW IN A SPIRAL CYCLONE
Applied Mathematics and Mechanics (English Edition), 2006, 27(2):247 253 c Editorial Committee of Appl. Math. Mech., ISSN 0253-4827 NUMERICAL SIMULATION OF THREE DIMENSIONAL GAS-PARTICLE FLOW IN A SPIRAL
More informationLES of the Sandia Flame D Using an FPV Combustion Model
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 82 (2015 ) 402 409 ATI 2015-70th Conference of the ATI Engineering Association LES of the Sandia Flame D Using an FPV Combustion
More informationA MODIFIED MODEL OF COMPUTATIONAL MASS TRANSFER FOR DISTILLATION COLUMN
A MODIFIED MODEL OF COMPUTATIONAL MASS TRANSFER FOR DISTILLATION COLUMN Z. M. Sun, X. G. Yuan, C. J. Liu, K. T. Yu State Key Laboratory for Chemical Engineering (Tianjin University) and School of Chemical
More informationBasic Concepts in Reactor Design
Basic Concepts in Reactor Design Lecture # 01 KBK (ChE) Ch. 8 1 / 32 Introduction Objectives Learning Objectives 1 Different types of reactors 2 Fundamental concepts used in reactor design 3 Design equations
More informationIMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE FLUCTUATIONS
Proceedings of the Combustion Institute, Volume 9, 00/pp. 1971 1978 IMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE
More informationFlow and added small-scale topologies in a turbulent premixed flame
Flow and added small-scale topologies in a turbulent premixed flame L. Cifuentes*, A. Kempf* and C. Dopazo** luis.cifuentes@uni-due.de *University of Duisburg-Essen, Chair of Fluid Dynamics, Duisburg -
More informationIntroduction. Mathematical model and experimental
COMPARISON OF THE PERFORMANCE OF A REVERSE FLOW REACTOR AND NETWORKS OF NON-STATIONARY CATALYTIC REACTORS FOR CATALYTIC COMBUSTION OF METHANE IN LEAN MIXTURES Miguel A. G. Hevia 1, Davide Fissore 2, Salvador
More informationCHAPTER 4 EXPERIMENTAL METHODS
47 CHAPTER 4 EXPERIMENTAL METHODS 4.1 INTRODUCTION The experimental procedures employed in the present work are oriented towards the evaluation of residence time distribution in a static mixer and the
More informationAPPLICATION OF THE CONDITIONAL QUADRATURE METHOD OF MOMENTS FOR THE SIMULATION OF COALESCENCE, BREAKUP AND MASS TRANSFER IN GAS-LIQUID STIRRED TANKS
14 th European Conference on Mixing Warszawa, 10-13 Septemer 2012 APPLICATION OF THE CONDITIONAL QUADRATURE METHOD OF MOMENTS FOR THE SIMULATION OF COALESCENCE, BREAKUP AND MASS TRANSFER IN GAS-LIQUID
More informationPrediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor
Journal of Engineering Science, Vol. 10, 95 105, 2014 Prediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor Vahid Akbari and Mohd.
More informationchromatographic reactors
This is an electronic version (author's version) of the paper: Fissore D., Barresi A. A., Botar-Jid C. C. (2006). NOx removal in forced unsteady-state chromatographic reactors. Chemical Engineering Science
More informationPDF modeling and simulation of premixed turbulent combustion
Monte Carlo Methods Appl. Vol. 4 No. 4 (8), pp. 343 377 DOI. / MCMA.8.6 c de Gruyter 8 PDF modeling and simulation of premixed turbulent combustion Michael Stöllinger and Stefan Heinz Abstract. The use
More informationAssessing Mixing Sensitivities for Scale-up. Matthew Jörgensen Atlanta, 10/20/2014
Assessing Mixing Sensitivities for Scale-up Matthew Jörgensen Atlanta, 10/20/2014 Nalas: Solutions that Scale! Engineering Services Broad range of services & capabilities supporting chemical process development
More informationTHE VERIFICATION OF THE TAYLOR-EXPANSION MOMENT METHOD IN SOLVING AEROSOL BREAKAGE
144 THERMAL SCIENCE, Year 1, Vol. 16, No. 5, pp. 144-148 THE VERIFICATION OF THE TAYLOR-EXPANSION MOMENT METHOD IN SOLVING AEROSOL BREAKAGE by Ming-Zhou YU a,b * and Kai ZHANG a a College of Science, China
More informationA NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER
SCIENTIFIC RESEARCH AND EDUCATION IN THE AIR FORCE-AFASES 2016 A NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER Alexandru DUMITRACHE*, Florin FRUNZULICA ** *Institute of
More informationSimulation of Turbulent Lifted Flames and their Transient Propagation
25 th ICDERS August 2-7th, 2015 Leeds, UK Simulation of Turbulent Lifted Flames and their Transient Propagation S. Ruan, Z. Chen, N. Swaminathan University of Cambridge Cambridge, UK 1 Introduction Turbulent
More informationNumerical Simulations of Hydrogen Auto-ignition in a Turbulent Co-flow of Heated Air with a Conditional Moment Closure
Numerical Simulations of Hydrogen Auto-ignition in a Turbulent Co-flow of Heated Air with a Conditional Moment Closure I. Stanković*, 1, A. Triantafyllidis, E. Mastorakos, C. Lacor 3 and B. Merci 1, 4
More informationVERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION
VERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION Cheng Z, Agranat V.M. and Tchouvelev A.V. A.V.Tchouvelev & Associates, Inc., 659 Spinnaker Circle, Mississauga, Ontario, Canada L5W R Hydrogenics
More informationPublished in: Proceedings of the Fluent Benelux User Group Meeting, 6-7 October 2005, Wavre, Belgium
The application of Flamelet Generated Manifolds in partailly-premixed flames Ramaekers, W.J.S.; Albrecht, B.A.; van Oijen, J.A.; de Goey, L.P.H.; Eggels, R.L.G.M. Published in: Proceedings of the Fluent
More informationEffect of External Recycle on the Performance in Parallel-Flow Rectangular Heat-Exchangers
Tamkang Journal of Science and Engineering, Vol. 13, No. 4, pp. 405 412 (2010) 405 Effect of External Recycle on the Performance in Parallel-Flow Rectangular Heat-Exchangers Ho-Ming Yeh Energy and Opto-Electronic
More informationOverview of Turbulent Reacting Flows
Overview of Turbulent Reacting Flows Outline Various Applications Overview of available reacting flow models LES Latest additions Example Cases Summary Reacting Flows Applications in STAR-CCM+ Ever-Expanding
More informationSIMULATION OF FLOW IN A RADIAL FLOW FIXED BED REACTOR (RFBR)
SIMULATION OF FLOW IN A RADIAL FLOW FIXED BED REACTOR (RFBR) Aqeel A. KAREERI, Habib H. ZUGHBI, *, and Habib H. AL-ALI * Ras Tanura Refinery, SAUDI ARAMCO, Saudi Arabia * Department of Chemical Engineering,
More informationAn unrestricted algorithm for accurate prediction of multiple-input multiple-output (MIMO) Wiener processes
Chemical and Biological Engineering Publications Chemical and Biological Engineering 2004 An unrestricted algorithm for accurate prediction of multiple-input multiple-output (MIMO Wiener processes Swee-Teng
More informationEvaluation of a Scalar Probability Density Function Approach for Turbulent Jet Nonpremixed Flames
THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN, KOREA, OCTOBER 7 30, 003 Evaluation of a Scalar Probability Density Function Approach for Turbulent Jet Nonpremixed Flames S. Noda, Y. Horitani and K.
More informationModelling and Simulation of a Batch Poly(Vinyl Chloride) Reactor
769 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 32, 23 Chief Editors: Sauro Pierucci, Jiří J. Klemeš Copyright 23, AIDIC Servizi S.r.l., ISBN 978-88-968-23-; ISSN 974-979 The Italian Association
More informationSTABILITY CONDITIONS AND OBSERVER DESIGN FOR A CONTINUOUS CRYSTALLIZER
STABILITY CONDITIONS AND OBSERVER DESIGN FOR A CONTINUOUS CRYSTALLIZER Juan Du and B. Erik Ydstie* Carnegie Mellon University Pittsburgh, PA 15213 Abstract The population balance equation is used to describe
More informationFlow Generated by Fractal Impeller in Stirred Tank: CFD Simulations
Flow Generated by Fractal Impeller in Stirred Tank: CFD Simulations Gunwant M. Mule and Amol A. Kulkarni* Chem. Eng. & Proc. Dev. Division, CSIR-National Chemical Laboratory, Pune 411008, INDIA *Corresponding
More informationPDF Modeling and Simulation of Premixed Turbulent Combustion
Monte Carlo Methods Appl. Vol. No. (), pp. 43 DOI 5 / MCMA.7. c de Gruyter PDF Modeling and Simulation of Premixed Turbulent Combustion Michael Stöllinger and Stefan Heinz Abstract. The use of probability
More informationFlamelet modelling of non-premixed turbulent combustion with local extinction and re-ignition
INSTITUTE OF PHYSICS PUBLISHING Combust. Theory Modelling 7 (2003) 317 332 COMBUSTION THEORY AND MODELLING PII: S1364-7830(03)38343-3 Flamelet modelling of non-premixed turbulent combustion with local
More informationTurbulent Flows. g u
.4 g u.3.2.1 t. 6 4 2 2 4 6 Figure 12.1: Effect of diffusion on PDF shape: solution to Eq. (12.29) for Dt =,.2,.2, 1. The dashed line is the Gaussian with the same mean () and variance (3) as the PDF at
More informationKinetic Modeling of Nanoprecipitation using CFD Coupled with a Population Balance
Chemical and Biological Engineering Publications Chemical and Biological Engineering 2010 Kinetic Modeling of Nanoprecipitation using CFD Coupled with a Population Balance Janine C. Cheng Iowa State University
More informationAn investigation of the performance of turbulent mixing models
Combustion and Flame 136 (24) 28 216 www.elsevier.com/locate/jnlabr/cnf An investigation of the performance of turbulent mixing models Zhuyin Ren and Stephen B. Pope Sibley School of Mechanical and Aerospace
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More informationConsistent turbulence modeling in a hybrid LES/RANS PDF method for non-premixed flames
Consistent turbulence modeling in a hybrid LES/RANS PDF method for non-premixed flames F. Ferraro, Y. Ge and M. Pfitzner Institut für Thermodynamik, Fakultät für Luft- und Raumfahrrttechnik Universität
More informationB.1 NAVIER STOKES EQUATION AND REYNOLDS NUMBER. = UL ν. Re = U ρ f L μ
APPENDIX B FLUID DYNAMICS This section is a brief introduction to fluid dynamics. Historically, a simplified concept of the boundary layer, the unstirred water layer, has been operationally used in the
More informationINTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING---
ILASS-Europe 2002 Zaragoza 9 11 September 2002 INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING--- Z. Tapia, A. Chávez e-mail: ztapia@imp.mx Instituto Mexicano del Petróleo Blvd. Adolfo Ruiz Cortines
More informationNUMERICAL PREDICTIONS OF DEPOSTION WITH A PARTICLE CLOUD TRACKING TECHNIQUE
Committed Individuals Solving Challenging Problems NUMERICAL PREDICTIONS OF DEPOSTION WITH A PARTICLE CLOUD TRACKING TECHNIQUE by James R. Valentine Reaction Engineering International Philip J. Smith Department
More informationEffect of near-wall treatments on airflow simulations
Proceedings of 29 International Conference on Computational Methods for Energy Engineering and Environment: ICCM3E. Sousse, Tunisia, 2-22 November, 29, pp. Effect of near-wall treatments on airflow simulations
More informationBASIC DESIGN EQUATIONS FOR MULTIPHASE REACTORS
BASIC DESIGN EQUATIONS FOR MULTIPHASE REACTORS Starting Reference 1. P. A. Ramachandran and R. V. Chaudhari, Three-Phase Catalytic Reactors, Gordon and Breach Publishers, New York, (1983). 2. Nigam, K.D.P.
More informationOutline. Definition and mechanism Theory of diffusion Molecular diffusion in gases Molecular diffusion in liquid Mass transfer
Diffusion 051333 Unit operation in gro-industry III Department of Biotechnology, Faculty of gro-industry Kasetsart University Lecturer: Kittipong Rattanaporn 1 Outline Definition and mechanism Theory of
More informationDetailed Chemical Kinetics in Multidimensional CFD Using Storage/Retrieval Algorithms
13 th International Multidimensional Engine Modeling User's Group Meeting, Detroit, MI (2 March 23) Detailed Chemical Kinetics in Multidimensional CFD Using Storage/Retrieval Algorithms D.C. Haworth, L.
More informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of
More informationNumerical Modeling of Anode Baking Process with COMSOL Multiphysics
Numerical Modeling of Anode Baking Process with COMSOL Multiphysics P. Nakate 1,2, W. Voogd de 3, D.J.P. Lahaye 1, C. Vuik 1, M. Talice 4 1. Department of applied mathematics, Delft University of technology,
More informationModeling of dispersed phase by Lagrangian approach in Fluent
Lappeenranta University of Technology From the SelectedWorks of Kari Myöhänen 2008 Modeling of dispersed phase by Lagrangian approach in Fluent Kari Myöhänen Available at: https://works.bepress.com/kari_myohanen/5/
More informationThe Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water Interface
Turbulence, Heat and Mass Transfer 4 K. Hanjalic, Y. Nagano and M. Tummers (Editors) 3 Begell House, Inc. The Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water
More informationA new jet-stirred reactor for chemical kinetics investigations Abbasali A. Davani 1*, Paul D. Ronney 1 1 University of Southern California
0 th U. S. National Combustion Meeting Organized by the Eastern States Section of the Combustion Institute pril 23-26, 207 College Park, Maryland new jet-stirred reactor for chemical kinetics investigations
More informationModeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet Model with Heat-Loss Effects
9 th U. S. National Combustion Meeting Organized by the Central States Section of the Combustion Institute May 17-20, 2015 Cincinnati, Ohio Modeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet
More informationConditional-moment Closure with Differential Diffusion for Soot Evolution in Fire
Chemical and Biological Engineering Conference Presentations and Proceedings Chemical and Biological Engineering 26 Conditional-moment Closure with Differential Diffusion for Soot Evolution in Fire J.
More informationHydrogen addition to the Andrussow process for HCN synthesis
Applied Catalysis A: General 201 (2000) 13 22 Hydrogen addition to the Andrussow process for HCN synthesis A.S. Bodke, D.A. Olschki, L.D. Schmidt Department of Chemical Engineering and Materials Science,
More informationWhere ω is the mass fraction of species k, Ω
THERMOCOUPLES POSITIONING FOR EARLY-WARNING DETECTION OF THERMAL RUNAWAY 1 MATTEO RIZZOTTO, 2 FEDERICO FLORIT, 3 RENATO ROTA, 4 VALENTINA BUSINI Politecnico di Milano, Dipartimento di Chimica, Materiali
More informationCOMSOL Multiphysics Simulation of 3D Single- hase Transport in a Random Packed Bed of Spheres
COMSOL Multiphysics Simulation of 3D Single- hase Transport in a Random Packed Bed of Spheres A G. Dixon *1 1 Department of Chemical Engineering, Worcester Polytechnic Institute Worcester, MA, USA *Corresponding
More informationNUMERICAL AND EXPERIMENTAL STUDIES OF LIQUID-LIQUID MIXING IN A KENICS STATIC MIXER
14 th European Conference on Mixing Warszawa, 10-13 September 2012 NUMERICAL AND EXPERIMENTAL STUDIES OF LIQUID-LIQUID MIXING IN A KENICS STATIC MIXER Zdzislaw Jaworski a, Halina Murasiewicz a, a Institute
More informationLES modeling of heat and mass transfer in turbulent recirculated flows E. Baake 1, B. Nacke 1, A. Umbrashko 2, A. Jakovics 2
MAGNETOHYDRODYNAMICS Vol. 00 (1964), No. 00, pp. 1 5 LES modeling of heat and mass transfer in turbulent recirculated flows E. Baake 1, B. Nacke 1, A. Umbrashko 2, A. Jakovics 2 1 Institute for Electrothermal
More informationErratum to: High speed mixture fraction and temperature imaging of pulsed, turbulent fuel jets auto igniting in high temperature, vitiated co flows
DOI 10.1007/s00348-015-2101-9 ERRATUM Erratum to: High speed mixture fraction and temperature imaging of pulsed, turbulent fuel jets auto igniting in high temperature, vitiated co flows Michael J. Papageorge
More informationEffect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process
, pp. 844 850 Effect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process Baokuan LI and Fumitaka TSUKIHASHI 1) Department of Thermal Engineering, The School
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationA copy can be downloaded for personal non-commercial research or study, without prior permission or charge
Oni, T. O., and Paul, M. C. (2014) Numerical simulation of turbulent heat transfer and fluid flow in different tube designs. In: World Congress on Engineering 2014, 2-4 July 2014, London, UK. Copyright
More informationPRESENTATION SLIDES: Hydrodynamic Scale-Up of Circulating Fluidized Beds
Refereed Proceedings The 12th International Conference on Fluidization - New Horizons in Fluidization Engineering Engineering Conferences International Year 2007 PRESENTATION SLIDES: Hydrodynamic Scale-Up
More informationMOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow
TRANSPORT PHENOMENA MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow Introduction to Turbulent Flow 1. Comparisons of laminar and turbulent flows 2. Time-smoothed equations of change for incompressible
More informationA First Course on Kinetics and Reaction Engineering Unit 12. Performing Kinetics Experiments
Unit 12. Performing Kinetics Experiments Overview Generating a valid rate expression for a reaction requires both a reactor and and an accurate mathematical model for that reactor. Unit 11 introduced the
More informationChemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore
Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Lecture No. #40 Problem solving: Reactor Design Friends, this is our last session
More information