Conditional statistics of inert droplet effects on turbulent combustion in reacting mixing layers

Direct numerical simulation (DNS) of turbulent reacting mixing layers laden with evaporating inert droplets is used to assess the droplet effects in the context of the conditional moment closure (CMC) for multiphase turbulent combustion. The temporally developing mixing layer has an initial Reynolds number of 1000 based on the vorticity thickness with more than 16 million Lagrangian droplets traced. An equivalent mixture fraction incorporating the inert vapour mass fractions clearly demonstrates the effects of vapour dilution on the mixture. Instantaneous fields and conditional statistics, such as the singly conditioned scalar dissipation rate, the gas temperature 〈 T _g|η 〉, conditional variance of the gas temperature 〈 T _g ^”2|η 〉 and conditional covariance between the fuel mass fraction and gas temperature 〈 Y _f ^” T _g ^”|η 〉 show considerable droplet effects. Comparison between the droplet-free and droplet-laden reacting mixing layer cases suggests significant extinction in the latter case. The resulting large conditional fluctuations around the conditional means contradict the basic assumption behind the first-order singly conditioned CMC. More sophisticated CMC approaches, such as doubly conditioned or second-order CMCs are, in principle, better able to incorporate the effects of evaporating droplets, but significant modelling challenges exist. The scalar dissipation rate doubly conditioned on the mixture fraction and a normalized gas temperature 〈 χ | η, ζ 〉 exemplifies the modelling complexity in the CMC of multiphase combustion.     

Modeling of a two-regime crystallization in a multicomponent lipid system under shear flow

The kinetics of phase transitions of milk fat triacylglycerols, as model multicomponent lipid systems, were studied under shear in a Couette cell at 17 °C, 17.5 °C and 20 °C under shear rates ranging from 0 to 2880s^-1 using synchrotron X-ray diffraction. Two-dimensional diffraction patterns were captured during the crystallization process. No effect of shear on onset time for phase α from the liquid was observed. Afterwards a two-regime crystallization process was observed. During the first regime, as observed in other systems, shear reduced the onset time of the phase transition from phase α to 2880s⁻. The model previously developed for palm oil (ODE model) worked well to describe this regime, confirming the general value of the proposed ODE model. However, the ODE model did not satisfactorily describe the second regime. We found that, as the system gets closer to equilibrium, the growth regime becomes controlled by diffusion, manifested by the kinetics following a dependence. This regime was found to be consistent with a mechanism combining step growth at a kink with progressive selection of the crystallizing moieties. This mechanism is in agreement with the displacement of the diffraction peak positions, which revealed how increased shear rate promotes the crystallization of the higher melting fraction affecting the composition of the crystallites.     

微液滴动力学特性的耗散粒子动力学模拟

对传统的耗散粒子动力学方法进行了改进.改进的耗散粒子动力学方法采用了包含远程吸引力和近距排斥力的保守力势函数,从而使得用耗散粒子动力学方法模拟多相流动成为可能.应用改进的耗散粒子动力学方法,对微尺度下液滴的形成及液滴在微重力下的大幅度振荡变形进行了数值模拟.计算结果表明,改进的耗散粒子动力学(DPD)方法能够有效地描述微尺度下液滴的动力学特性,对研究复杂流体多相流动有着重要的意义. 关键词： 多相流 微液滴 耗散粒子动力学(DPD)方法 保守力势函数     

Simulating bubble dynamics in a buoyant system

Multiphase flows are critical components of many physical systems; however, numerical models of multiphase flows with large parameter gradients can be challenging. Here, two different numerical methods, volume of fluid (VOF) and smoothed particle hydrodynamics (SPH), are used to model the buoyant rise of isolated gas bubbles through quiescent fluids for a range of Bond and Reynolds numbers. The VOF is an Eulerian grid–based method, whereas the SPH is Lagrangian and mesh free. Each method has unique strengths and weaknesses, and a comparison of the two approaches as applied to multiphase phenomena has not previously been performed. The VOF and SPH simulations are compared, verified, and validated. Results using two-dimensional VOF and SPH simulations are similar to each other and are able to reproduce numerical benchmarks and experimental results for sufficiently large Morton and Reynolds numbers. It is also shown that at low Reynolds numbers, the two methods, SPH and VOF, diverge in the transient regime of the bubble rise. Regimes that require simulations capable of representing three-dimensional drag are identified as well as regimes in which results from VOF and SPH diverge.     

Effects of magnetic fluids on crystallization characterizations in a multi-component and multiphase system

In this study, experiments are carried out on the effects of magnetic fluids on the crystallization char- acterizations in a multi-component and multiphase system, which contains the liquid and the vapor of HCFC141b, water, water vapor, and gas hydrates. The mass transfer phenomena between the phase interfaces of water-HCFC141b and water-vapor are also researched. The experimental results show that in the presence of a rotary magnetic field, magnetic fluids can remarkably enhance the heat and mass transfer between phase interfaces and, therefore, improve the performance of crystallization, especially in improving the formation temperature and velocity.     

螺旋轴流式多相泵外特性实验研究

本文在建立多相混输泵实验装置的基础上,以自行设计的螺旋轴流式多相泵原理机为实验研究对象,分别以纯水、空气—水为实验介质,系统地研究了多相泵的泵轴转速、进口压力、混合物特性等对多相输送性能的影响,为多相泵的进一步开发研制及完善设计提供了试验资料。     

Computational fluid dynamics simulations of aerosol behavior in a high-speed (Heubach) rotating drum dustiness tester

Potential exposure from hazardous dust may be assessed by evaluating the dustiness of the powders being handled. Dustiness is the tendency of a powder to aerosolize with a given input of energy. Previously we used computational fluid dynamics (CFD) to numerically investigate the flow inside the European Standard (EN15051) rotating drum dustiness tester during its operation. The present work extends those CFD studies to the widely used Heubach rotating drum. Air flow characteristics are investigated within the Abe-Kondoh-Nagano k-epsilon turbulence model; the aerosol is incorporated via a Euler-Lagrangian multiphase approach. The air flow inside these drums consists of a well-defined axial jet penetrating relatively quiescent air. The spreading of the Heubach jet results in a fraction of the jet recirculating as back-flow along the drum walls; at high rotation rates, the axial jet becomes unstable. This flow behavior qualitatively differs from the stable EN15051 flow pattern. The aerodynamic instability promotes efficient mixing within the Heubach drum, resulting in higher particle capture efficiencies for particle sizes d < 80 μm.     

Application of the full approximation storage method to the numerical simulation of two-dimensional steady incompressible viscous multiphase flows

In recent years multigrid algorithms have been applied to increasingly difficult systems of partial differential equations and major improvements in both speed of convergence and robustness have been achieved. Problems involving several interacting fluids are of great interest in many industrial applications, especially in the process and petro-chemical sectors. However, the multifluid version of the Navier–Stokes equations is extremely complex and represents a challenge to advanced numerical algorithms. In this paper, we describe an extension of the full approximation storage (FAS) multigrid algorithm to the multifluid equations. A number of special issues had to be addressed. The first was the development of a customised, non-linear, coupled relaxation scheme for the smoothing step. Automatic differentiation was used to facilitate the coding of a robust, globally convergent quasi-Newton method. It was also necessary to use special inter-grid transfer operators to maintain the realisability of the solution. Algorithmic details are given and solutions for a series of test problems are compared with those from a widely validated, commercial code. The new approach has proved to be robust; it achieves convergence without resorting to specialised initialisation methods. Moreover, even though the rate of convergence is complex, the method has achieved very good reduction factors: typically five orders of magnitude in 50 cycles. © 1998 John Wiley & Sons, Ltd.