Search

Droplet dispersion in liquid–liquid systems is a crucial step in many unit operations throughout the chemical, food, and pharmaceutic industries, where improper operation causes billions of dollars of loss annually. A theoretical background for the description of droplet breakup has been established, but many assumptions are still unconfirmed by experimental observations. In this investigation, a von Kármán swirling flow device was used to produce homogeneous, low-intensity turbulence suitable for carrying out droplet breakage experiments using optical image analysis. Individual droplets of known, adjustable, and repeatable sizes were introduced into an isotropic turbulent flow field providing novel control over two of the most important factors impacting droplet breakage: turbulence dissipation rate and parent droplet size. Introducing droplets one at a time, large data sets were gathered using canola, safflower, and sesame oils for the droplet phase and water as the continuous phase. Automated image analysis was used to determine breakage time, breakage probability, and child droplet size distribution for various turbulence intensities. Breakage time and breakage probability were observed to increase with increasing parent droplet size, consistent with the classic and widely used Coulaloglou–Tavlarides breakage model (C–T model). The shape of the child drop size distribution function was found to depend upon the size of the parent droplet. [ABSTRACT FROM AUTHOR]

Improving the efficiency of the production of uniformly sized functional nanoparticles for pharmaceutical and agricultural applications has been a problem of great interest. The macroscale multi-inlet vortex Reactor (macro-MIVR) could potentially be used for this purpose due to its ability to achieve the rapid mixing necessary for the flash nanoprecipitation nanoparticle fabrication technique. In the presented work, the coherent structures, a key contributor to the turbulent mixing, were investigated for the turbulent swirling flow within the MIVR. The two-point spatial correlations of velocity and concentration fluctuations at various basepoints were measured from instantaneous velocity and concentration fields obtained using simultaneous stereoscopic particle image velocimetry and planar laser-induced fluorescence. The basepoint locations were chosen as the middle and at the edge of the partially mixed concentration spiral arms. The correlations were found to be elliptical in shape, inclined, and peaked at the basepoints. A region near the basepoint was positively correlated and was surrounded by negatively correlated regions. Autocorrelations of concentration were also elliptical and curved toward the center of the reactor. The linear stochastic estimation was used to interpret the coherent structure features that would result in the observed spatial correlations. The linear stochastic estimates of the velocity fields were computed directly from the cross correlations of the tangential velocity fluctuations with the concentration fluctuations. The estimated conditional velocity fields revealed obliquely oriented counter-rotating vortical structures that stir the fluid from high-concentration regions to low-concentration regions, and the orientation of these vortical structures depended on the local concentration gradient. [ABSTRACT FROM AUTHOR]

Modeling particle-laden turbulent flows at high volume fractions requires accounting for the coupling between phases. The latter is often a sensitive point, and proper closure of the exchange and production terms due to the presence of particles is not straightforward. In the present work, a Lagrangian probability-density-function model developed for homogeneous cluster-induced turbulence is extended to a channel flow. The key features are consistent two-way coupling and the decomposition of the particle velocity into spatially correlated and uncorrelated components, which is crucial for dense flows and which allows dealing with collisions from a statistical point of view. A numerical scheme for the coupled solution of the stochastic differential equations for the particles and a Reynolds-stress model for the fluid is developed. Tests with tracer particles without two-way coupling are done to assess the validity and the consistency of the numerical scheme. Finally, two sets of numerical simulations with particles with different diameters in a turbulent channel flow at a shear Reynolds of R e τ = 300 are reported. The effect of two-way coupling by varying the mass loading of the dispersed phase in the mass-loading range φ = 0–2 is analyzed, and the results are compared to previous Eulerian–Lagrangian and Eulerian–Eulerian direct-numerical simulation (DNS) studies. Mean velocities and turbulent kinetic energy show good agreement with DNS, especially regarding the trend with respect to mass loading. Consistent with prior work, increased mass loading causes a drastic reduction of turbulent kinetic energy in the range φ = 0–2. [ABSTRACT FROM AUTHOR]

In the context of modeling turbulent scalar mixing using probability density function (PDF) methods, the treatment of molecular mixing is of paramount importance. The conditional moment closure (CMC) offers a high-fidelity description for molecular mixing in nonpremixed flows. Recent work has demonstrated that first-order CMC can be implemented numerically using the moments of the conditioning variable and first-order joint moments of the scalar of interest. When solving the CMC using, for example, quadrature-based moment methods (QBMM), a functional form must be chosen for the conditional scalar dissipation rate (CSDR) of the conditioning variable. In prior work, the CSDR was chosen to produce a β-PDF for the conditioning variable (mixture fraction) at steady state. This choice has the advantage that the system of moment equations used in QBMM-CMC can be written in closed form. In this work, an alternative choice for the CSDR is investigated, namely, the amplitude mapping closure (AMC). With the AMC, the moment equations can be closed using the quadrature method of moments incorporated into a realizable ordinary differential equation solver. Results are compared with the β-CSDR closure for binary, passive scalar mixing in homogeneous single- and disperse-phase turbulent flows. It is also demonstrated that the moment formulation of CMC provides a straightforward method for modeling the effect of differential diffusion in the context of CMC. [ABSTRACT FROM AUTHOR]

The conditional velocity time averages (⟨Ui|ξ⟩) and conditional mixture fraction time averages (⟨Φ|ωi⟩) were computed based on the Eulerian approach from the experimental data measured in a macroscale multi-inlet vortex chemical reactor. The conditioning events were determined by equally sized intervals of the sample space variable for the mixture fraction (ξ) and the velocity vector (ωi). The experimental data, which consisted of instantaneous velocities and concentration fields for two Reynolds numbers (Re = 3250 and 8125), were acquired using the simultaneous stereoscopic particle image velocimetry (stereo-PIV) and planar laser induced fluorescence techniques. Two mathematical models, the linear approximation and probability density function (PDF) gradient diffusion, were validated by experimental results. The results of the velocity conditioned on the mixture fraction demonstrated that the linear model works well in a low turbulence region away from the reactor center. Near the reactor center, high velocity gradients coupled with low concentration gradients reduce the accuracy of the linear model predictions. Nevertheless, an excellent agreement was found for the conditional events within ±2Φrms (mixture fraction root mean square). Due to lower concentration gradient in the tangential direction, the linear model better predicted the tangential velocity component for all locations investigated. The PDF model with an isotropic turbulent diffusivity performed inadequately for the tangential and axial velocity components. A modified version of the PDF model that considers the three components of the turbulent diffusivity produced a better agreement with the experimental data especially in the spiral arms regions of significant concentration gradients. Furthermore, the mixture fraction conditioned on the velocity vector components showed a more linear behavior near the reactor center, where the PDF of the mixture fraction is a Gaussian distribution. As the concentration gradients became prominent away from the reactor, ⟨Φ|ωi⟩ also deviated from the linear pattern. This was especially remarkable for the mixture fraction conditioned on the tangential velocity. The overall prediction of ⟨Φ|ωi⟩ improves at higher Reynolds number as the fluid mixing is enhanced. [ABSTRACT FROM AUTHOR]

In Part I, simulations of strongly coupled fluid-particle flow in a vertical channel were performed with the purpose of understanding, in general, the fundamental physics of wall-bounded multiphase turbulence and, in particular, the roles of the spatially correlated and uncorrelated components of the particle velocity. The exact Reynoldsaveraged (RA) equations for high-mass-loading suspensions were presented, and the unclosed terms that are retained in the context of fully developed channel flow were evaluated in an Eulerian-Lagrangian (EL) framework. Here, data from the EL simulations are used to validate a multiphase Reynolds-stress model (RSM) that predicts the wall-normal distribution of the two-phase, one-point turbulence statistics up to second order. It is shown that the anisotropy of the Reynolds stresses both near the wall and far away is a crucial component for predicting the distribution of the RA particle-phase volume fraction. Moreover, the decomposition of the phase-average (PA) particle-phase fluctuating energy into the spatially correlated and uncorrelated components is necessary to account for the boundary conditions at the wall. When these factors are properly accounted for in the RSM, the agreement with the EL turbulence statistics is satisfactory at first order (e.g., PA velocities) but less so at second order (e.g., PA turbulent kinetic energy). Finally, an algebraic stress model for the PA particle-phase pressure tensor and the Reynolds stresses is derived from the RSM using the weak-equilibrium assumption. [ABSTRACT FROM AUTHOR]

Simulations of strongly coupled (i.e., high-mass-loading) fluid-particle flows in vertical channels are performed with the purpose of understanding the fundamental physics of wall-bounded multiphase turbulence. The exact Reynolds-averaged (RA) equations for high-mass-loading suspensions are presented, and the unclosed terms that are retained in the context of fully developed channel flow are evaluated in an Eulerian-Lagrangian (EL) framework for the first time. A key distinction between the RA formulation presented in the current work and previous derivations of multiphase turbulence models is the partitioning of the particle velocity fluctuations into spatially correlated and uncorrelated components, used to define the components of the particle-phase turbulent kinetic energy (TKE) and granular temperature, respectively. The adaptive spatial filtering technique developed in our previous work for homogeneous flows [J. Capecelatro, O. Desjardins, and R. O. Fox, "Numerical study of collisional particle dynamics in cluster-induced turbulence," J. Fluid Mech. 747, R2 (2014)] is shown to accurately partition the particle velocity fluctuations at all distances from the wall. Strong segregation in the components of granular energy is observed, with the largest values of particle-phase TKE associated with clusters falling near the channel wall, while maximum granular temperature is observed at the center of the channel. The anisotropy of the Reynolds stresses both near the wall and far away is found to be a crucial component for understanding the distribution of the particle-phase volume fraction. In Part II of this paper, results from the EL simulations are used to validate a multiphase Reynolds-stress turbulence model that correctly predicts the wall-normal distribution of the two-phase turbulence statistics. [ABSTRACT FROM AUTHOR]

There was also an error in Table V. Table V appears correct below: TABLE V. First-order temporal discretization (Euler scheme). Publisher's Note: "A Lagrangian probability-density-function model for turbulent particle-laden channel flow in the dense regime" [Phys. [Extracted from the article]

A technique for visualizing and quantifying reactive mixing for laminar and turbulent flow in a microscale chemical reactor using confocal-based microscopic laser induced fluorescence (confocal μ-LIF) was demonstrated in a microscale multi-inlet vortex nanoprecipitation reactor. Unlike passive scalar μ-LIF, the reactive μ-LIF technique is able to visualize and quantify micromixing effects. The confocal imaging results indicated that the flow in the reactor was laminar and steady for inlet Reynolds numbers of 10, 53, and 93. Mixing and reaction were incomplete at each of these Reynolds numbers. The results also suggested that although mixing by diffusion was enhanced near the midplane of the reactor at Rej = 53 and 93 due to very thin bands of acidic and basic fluid forming as the fluid spiraled towards the center of the reactor, near the top, and bottom walls of the reactor, the lower velocities due to fluid friction with the walls hindered the formation of these thin bands, and, thus, resulted in large regions of unmixed and unreacted fluid. At Rej = 240, the flow was turbulent and unsteady. The mixing and reaction processes were still found to be incomplete even at this highest Reynolds number. At the reactor midplane, the flow images at Rej = 240 showed unmixed base fluid near the center of the reactor, suggesting that just as in the Rej = 53 and 93 cases, lower velocities near the top and bottom walls of the reactor hinder the mixing and rection of the acidic and basic streams. Ensemble averages of line-scan profiles for the Rej = 240 were then calculated to provide statistical quantification of the microscale mixing in the reactor. These results further demonstrate that even at this highest Reynolds number investigated, mixing and reaction are incomplete. Visualization and quantification of micromixing using this reactive μ-LIF technique can prove useful in the validation of computational fluid dynamics models of micromixing within microscale chemical reactors. [ABSTRACT FROM AUTHOR]

Experimental results for conditional statistics in a confined turbulent wake were compared to the linear model and the transported probability density function (PDF) model for predicting the velocity conditioned on scalar mixture fraction and to a linear model for predicting the scalar mixture fraction conditioned on velocity. For velocity conditioned on scalar mixture fraction, the linear model was found to accurately predict the conditional transverse velocity, but for the conditional streamwise velocity, agreement between the experiments and the linear model was poor. The transverse conditional velocity was accurately predicted by the transported PDF model, but for the streamwise conditional velocity, a modified transported PDF model that incorporated two components of the scalar diffusivity tensor had to be used. Finally, it was found that a linear model for predicting the fluctuating scalar conditioned on velocity performed poorly due to the strongly non-Gaussian scalar PDF at the downstream locations investigated. [ABSTRACT FROM AUTHOR]

Attempts to combine the generalized IEM (GIEM) model with the mapping closure to model the molecular mixing terms in the probability distribution function balance equation. Assumption of a linear relationship between the rates of change of the reactive scalars and an inert scalar during the mixing step.

Mass production of functional nanoparticles may be realized through flash nanoprecipiation in microscale reactors such as the multi-inlet vortex reactor (MIVR). A comprehensive understanding of mixing in the MIVR is required for process control and reactor design. Mixing in the MIVR is studied using a technique coupling laser induced fluorescence with confocal laser scanning microscopy. It is shown to provide meaningful qualitative and statistical data of the scalar field for analysis and comparison with numerical simulations. Data were collected for four flow rates, showing that mixing is incomplete even at the highest flow rate. [ABSTRACT FROM AUTHOR]

Microscale reactors capable of generating turbulent flow are used in Flash NanoPrecipitation, an approach to produce functional nanoparticles with unique optical, mechanical and chemical properties. Microreactor design and optimization could be greatly enhanced by developing reliable computational models of the nanoprecipitation process. A microscale multi-inlet vortex nanoprecipitation reactor was investigated using microscopic particle image velocimetry and computational fluid dynamics. Velocity data such as the mean velocity and turbulent kinetic energy displayed good agreement between experiment and simulation over flow conditions ranging from fully laminar to turbulent, demonstrating the accuracy of the simulation model over the entire turbulent transition range. [ABSTRACT FROM AUTHOR]

A technique for visualizing and quantifying reactive mixing for laminar and turbulent flow in a microscale chemical reactor using confocal-based microscopic laser induced fluorescence (confocal μ-LIF) was demonstrated in a microscale multi-inlet vortex nanoprecipitation reactor. Unlike passive scalar μ-LIF, the reactive μ-LIF technique is able to visualize and quantify micromixing effects. The confocal imaging results indicated that the flow in the reactor was laminar and steady for inlet Reynolds numbers of 10, 53, and 93. Mixing and reaction were incomplete at each of these Reynolds numbers. The results also suggested that although mixing by diffusion was enhanced near the midplane of the reactor at Rej = 53 and 93 due to very thin bands of acidic and basic fluid forming as the fluid spiraled towards the center of the reactor, near the top, and bottom walls of the reactor, the lower velocities due to fluid friction with the walls hindered the formation of these thin bands, and, thus, resulted in large regions of unmixed and unreacted fluid. At Rej = 240, the flow was turbulent and unsteady. The mixing and reaction processes were still found to be incomplete even at this highest Reynolds number. At the reactor midplane, the flow images at Rej = 240 showed unmixed base fluid near the center of the reactor, suggesting that just as in the Rej = 53 and 93 cases, lower velocities near the top and bottom walls of the reactor hinder the mixing and rection of the acidic and basic streams. Ensemble averages of line-scan profiles for the Rej = 240 were then calculated to provide statistical quantification of the microscale mixing in the reactor. These results further demonstrate that even at this highest Reynolds number investigated, mixing and reaction are incomplete. Visualization and quantification of micromixing using this reactive μ-LIF technique can prove useful in the validation of computational fluid dynamics models of micromixing within microscale chemical reactors.