A collision algorithm for anisotropic disperse flows based on ellipsoidal parcel representations

In this study, a new approach to the prediction of collision probabilities is proposed, which is specifically derived for Lagrangian Monte Carlo simulations of highly anisotropic disperse flows. The collision probability formulation eliminates the influence of the collision mesh by introduction of a parcel diameter. Two non-parametric parcel diameter estimators are given, an isotropic parcel diameter estimator for homogeneous droplet clouds, and an anisotropic parcel diameter estimator for inhomogeneous droplet clouds of reduced dimensionality, which obtains an ellipsoidal parcel representation from a weighted principal component analysis of the parcel population. The new formulation effectively reduces the parameters for spatial resolution from two (mesh size, number of parcels) to one (number of parcels only). For validation, simple synthetic test cases are run that allow a comparison of the collision probability formulation to purely stochastic and purely deterministic schemes. Dependence on numerical parameters is tested exemplary with a highly anisotropic hollow-cone spray for gasoline direct injection.     

Effects of Spray and Turbulence Modelling on the Mixing and Combustion Characteristics of an n-heptane Spray Flame Simulated with Dynamic Adaptive Chemistry

Accurate modelling of spray combustion process is essential for efficiency improvement and emissions reduction in practical combustion engines. In this work, both unsteady Reynolds-averaged Navier-Stokes (URANS) simulations and large eddy simulations (LES) are performed to investigate the effects of spray and turbulence modelling on the mixing and combustion characteristics of an n-heptane spray flame in a constant volume chamber at realistic conditions. The non-reacting spray process is first simulated with URANS to investigate the effects of entrainment gas-jet model on the penetration characteristics and fuel vapor distributions. It is found that the droplet motion near the nozzle has significant influence on the fuel vapor distribution, while the liquid penetration length is controlled by the evaporation process and insensitive to gas-jet model. For the case considered, both URANS with the gas-jet model and large eddy simulations can properly predict the vapor penetration. For the combustion characteristics, it is found that LES yields better predictions in the global combustion characteristics. The URANS with gas jet model yields a comparable flame length and lift-off-length (LOL) to LES, but results in a larger ignition delay time compared to the experimental data. Another focus of this work is to qualify the convergence characteristics of the dynamic adaptive chemistry (DAC) method in these transient combustion simulations, where DAC is applied to reduce the mechanism locally and on-the-fly to accelerate chemistry calculations. The instantaneous flame structures and global combustion characteristics such as ignition delay time, flame lift-off length and emissions are compared between simulations with and without DAC. For URANS, good agreements are observed both on instantaneous flame structures and global characteristics. For LES, it is shown that the errors incurred by DAC are small for scatter distributions in composition space and global combustion characteristics, while they may significantly affect instantaneous flame structures in physical space. The study reveals that for DAC application in transient simulations, global or statistic information should be used to assess the accuracy, such as manifolds in composition space, conditional quantities and global combustion characteristics. For the cases investigated, a speed-up factor of more than two is achieved by DAC with a 92-species skeletal mechanism with less than 0.2 % and 3.0 % discrepancy in ignition delay and LOL, respectively.     

A Comprehensive Study of Effects of Mixing and Chemical Kinetic Models on Predictions of n-heptane Jet Ignitions with the PDF Method

The composition Probability Density Function (PDF) model is coupled with a Reynolds-averaged k???ε turbulence model and three computationally efficient, yet widely used chemical mechanisms to simulate transient n-heptane spray injection and ignition in a high temperature and high density ambient fluid. Molecular diffusion is modelled by three mixing models, namely the interaction by exchange with the mean (IEM), modified Curl (MC) and Euclidean minimum spanning trees (EMST) models. The liquid phase is modelled by a discrete phase model (DPM). This represents among the first applications of the PDF method in practical diesel engine conditions. A non-reacting case is first considered, with the focus on the ability of the model to capture the spray structure, e.g., vapour penetration and liquid length, fuel mixture fraction and its variance. Reacting cases are then investigated to compare and evaluate the three different chemical mechanisms and the three mixing models. It is concluded that the EMST mixing model in conjunction with a reduced chemical kinetic model (Lu et al., Combust Flame 156(8):1542–1551, 2009) performs the best among the options considered. The sensitivity of the results to the choice of the mixing constant is also studied to understand its effect on the flame ignition and stabilisation. Finally, the PDF model is compared to a well-mixed model that assumes turbulent fluctuations are negligible, which has been widely used in the diesel spray combustion community. Significant structural differences in the modelled flame are revealed comparing the PDF method with the well-mixed model. Quantitatively, the PDF model exhibits excellent agreement with the measurements and shows much better results than the well-mixed model in all ambient O₂ and temperature conditions.     

Numerical modeling and experimental measurements of a high speed solid-cone water spray for use in fire suppression applications

Experimental measurements and numerical simulations of a high-speed water spray are presented. The numerical model is based on a stochastic separated flow technique that includes submodels for droplet dynamics, heat and mass transfer, and droplet–droplet collisions. Because the spray characteristics near the nozzle are difficult to ascertain, a new method for initialization of particle diameter size is developed that assumes a Rosin–Rammler distribution for droplet size, which correctly reproduces experimentally measured Sauter and arithmetic mean diameters. By relating the particle initialization to lower moments of the droplet statistics, it is possible to take advantage of measurements without substantial penalties associated with the greater experimental uncertainty of individual droplet measurements. Overall, very good agreement is observed in the comparisons of experimental measurements to computational predictions for the streamwise development of mean drop size and velocity. In addition, the importance of modeling droplet–droplet collisions is highlighted with comparison of selected droplet–droplet collision models.     

Effects of turbulence enhancement on combustion process using a double injection strategy in direct-injection spark-ignition (DISI) gasoline engines

Direct-injection spark-ignition (DISI) gasoline engines have been spotlighted due to their high thermal efficiency. Increase in the compression ratio that result from the heat absorption effect of fuel vaporization induces higher thermal efficiency than found in port fuel injection (PFI) engines. Since fuel is injected at the cylinder directly, various fuel injection strategies can be used. In this study, turbulent intensity was improved by a double injection strategy while maintaining mixture homogeneity. To analyze the turbulence enhancement effects using the double injection strategy, a side fuel injected, homogeneous-charge-type DISI gasoline engine with a multi-hole-type injector was utilized. The spray model was evaluated using experimental data for various injection pressures and the combustion model was evaluated for varied ignition timing. First and second injection timing was swept by 20 degree interval. The turbulent kinetic energy and mixture inhomogeneity index were mapped. First injection at the middle of the intake stroke and second injection early in the compression stroke showed improved turbulent characteristics that did not significantly decrease with mixture homogeneity. A double injection case that showed improved turbulent intensity while maintaining an adequate level of mixture homogeneity and another double injection case that showed significantly improved turbulent intensity with a remarkable decrease in mixture homogeneity were considered for combustion simulation. We found that the improved turbulent intensity increased the flame propagation speed. Also, the mixture homogeneity affected the pressure rise rate.     

Adaptive collision meshing and satellite droplet formation in spray simulations

Improved numerical methods and physical models have been applied to droplet collision modeling. Numerically, an adaptive collision mesh method is developed to calculate collision rate. This method produces a collision mesh that is independent of the gas phase mesh and adaptively refined according to local parcel number density. An existing model describing the satellite droplet formation during the collision process is improved to reflect the experimental findings that the satellite droplets are much smaller than the parent droplets. The adaptive collision mesh and the improved satellite model have been used to simulate three impinging spray experiments. The model was able to qualitatively predict the occurrence of small satellite drops and bi-modal post-collision drop size distributions. The effect of the collision mesh and the satellite droplet model on a high-speed non-evaporating diesel spray is also assessed.     

Modeling droplet dispersion and interphase turbulent kinetic energy transfer using a new dual-timescale Langevin model

Dispersion of spray droplets and the modulation of turbulence in the ambient gas by the dispersing droplets are two coupled phenomena that are closely linked to the evolution of global spray characteristics, such as the spreading rate of the spray and the spray cone angle. Direct numerical simulations (DNS) of turbulent gas flows laden with sub-Kolmogorov size particles, in the absence of gravity, report that dispersion statistics and turbulent kinetic energy (TKE) evolve on different timescales. Furthermore, each timescale behaves differently with Stokes number, a non-dimensional flow parameter (defined in this context as the ratio of the particle response time to the Kolmogorov timescale of turbulence) that characterizes how quickly a particle responds to turbulent fluctuations in the carrier or gas phase. A new dual-timescale Langevin model (DLM) composed of two coupled Langevin equations for the fluctuating velocities, one for each phase, is proposed. This model possesses a unique feature that the implied TKE and velocity autocorrelation in each phase evolve on different timescales. Consequently, this model has the capability of simultaneously predicting the disparate Stokes number trends in the evolution of dispersion statistics, such as velocity autocorrelations, and TKE in each phase. Predictions of dispersion statistics and TKE from the new model show good agreement with published DNS of non-evaporating and evaporating droplet-laden turbulent flow.     

Numerical and experimental investigation on droplet dynamics and dispersion of a jet engine injector

In the case of turbine combustors operating with liquid fuel the combustion process is governed by the liquid fuel atomization and its dispersion in the combustion chamber. By highly unsteady flow field conditions the transient interaction between the liquid and the gaseous phase is of interest, because it results in a temporal variation of air–fuel ratio which leads to a fluctuating temperature distribution. The objective of this research was the investigation of transient flow field phenomena (e.g. large coherent structures) on droplet dynamics and dispersion of an isothermal flow (of inert water droplets) as a necessary first step towards a full analysis of spray combustion in real-life devices. The advanced injector system for lean jet engine combustors PERM (Partial Evaporated Rapid Mixing) was applied, generating a dilute polydispersed spray in a swirled flow field. Experiments were performed using Phase Doppler Anemometry (PDA) and a patternator to determine the droplet polydispersity, concentration maps, and velocity profiles in the flow. An important finding is the effect of large-scale coherent structures due mainly to the precessing of the vortex core (PVC) of the swirling air jet on the particle dispersion patterns. The experimental results then serve as reference data to assess the accuracy of the Eulerian–Lagrangian computations using a Large Eddy Simulation (LES), a Unsteady Reynolds-Average Navier–Stokes Simulation (URANS) and two simplified (steady-state) simulations. There, a simplified droplet injection model was used and the required boundary conditions of injected droplet sizes were obtained from measurements. Important transient effects of deterministic droplet separation observed during experiments, could be perfectly replicated with this injection model. It is convincingly shown, through extensive computations, that the resolution of instantaneous vortical structures is indeed crucial; hence the LES, or a reasonably-well resolved URANS are preferred over the steady-state solutions with additional, stochastic-type, turbulent dispersion models.     

Nonlinear elasticity modeling of biogels

A second-order nonlinear elastic model is developed for biogels, whereby the elastic constants are estimated from the free energy densities for polymer networks and polymer-water mixtures. For typical values of the constants, important physical insights are revealed: a biogel stiffens under tension, contracts longitudinally under torsion (inverse Poynting effect), and experiences a negative normal stress under simple shear. This approach of extracting nonlinear elastic constants works for any constitutive law, and has the advantage of exploring physical phenomena without recourse to computationally intensive simulations.     

Experimental investigation on splashing and nonlinear fingerlike instability of large water drops

The fluid physics of the splashing and spreading of a large-scale water drop is experimentally observed and investigated. New phenomena of drop impact that differ from the conventional Rayleigh–Taylor instability theory are reported. Our experimental data shows good agreement with previous work at low Weber number but the number of fingers or instabilities begins to deviate from the R–T equation of Allen at high Weber numbers. Also observed were multiple waves (or rings) on the spreading liquid surface induced from pressure bouncing (or pulsation) within the impacting liquid. The first ring is transformed into a radially ejecting spray whose initial speed is accelerated to a velocity of 4–5 times that of the impacting drop. This first ring is said to be “splashing,” and its structure is somewhat chaotic and turbulent, similar to a columnar liquid jet surrounded by neighboring gas jets at relatively high impact speed. At lower impact speeds, splashing occurs as a crown-shaped cylindrical sheet. A second spreading ring is observed that transforms into fingers in the circumferential direction during spreading. At higher Weber number, the spreading of a third ring follows that of the second. This third ring, induced by the pressure pulsation, overruns and has fewer fingers than the second, which is still in a transitional spreading stage. Several important relationships between the drop impact speed, the spray ejection speed of the first ring, and the number of fingers of the second and third rings are presented, based on data acquired during a set of drop impact experiments. Issues related to the traditional use of the R–T instability are also addressed.     

An algebraic closure for the DNS of fiber-induced turbulent drag reduction in a channel flow

An algebraic closure for the non-Newtonian Navier–Stokes equations is presented which accounts for the effect of a dilute fiber suspension. The model is intended to be used in simulations of turbulent drag reduction by fiber additives, and can be considered as a computationally efficient alternative to the existing rheological models for fiber suspensions in turbulent wall-bounded flows. It is based on the assumption that the suspended elongated particles are aligned with the local velocity fluctuation vector. The model is proved to be Galilean invariant. One-way coupled simulations and comparison with a direct solution of the underlying Fokker–Planck equation show a considerable improvement over an existing and comparable model. Finally, two-way coupled simulations demonstrate that the model predicts flow statistics that are in very good agreement with those obtained by the moment approximation approach. Interestingly, the model is realistic in terms of the polymer concentration. Using the proposed model, the cost of simulating a drag-reduced flow in terms of CPU-time is slightly more than that of a Newtonian flow.     

Turbulent collision of inertial particles: Point-particle based,hybrid simulations and beyond

Point-particle based direct numerical simulation (PPDNS) has been a productive research tool for studying both single-particle and particle-pair statistics of inertial particles suspended in a turbulent carrier flow. Here we focus on its use in addressing particle-pair statistics relevant to the quantification of turbulent collision rate of inertial particles. PPDNS is particularly useful as the interaction of particles with small-scale (dissipative) turbulent motion of the carrier flow is mostly relevant. Furthermore, since the particle size may be much smaller than the Kolmogorov length of the background fluid turbulence, a large number of particles are needed to accumulate meaningful pair statistics. Starting from the relative simple Lagrangian tracking of so-called ghost particles, PPDNS has significantly advanced our theoretical understanding of the kinematic formulation of the turbulent geometric collision kernel by providing essential data on dynamic collision kernel, radial relative velocity, and radial distribution function. A recent extension of PPDNS is a hybrid direct numerical simulation (HDNS) approach in which the effect of local hydrodynamic interactions of particles is considered, allowing quantitative assessment of the enhancement of collision efficiency by fluid turbulence. Limitations and open issues in PPDNS and HDNS are discussed. Finally, on-going studies of turbulent collision of inertial particles using large-eddy simulations and particle-resolved simulations are briefly discussed.     

Modeling of binary droplet collisions for application to inter-impingement sprays

In this paper, we focused on modeling the collision phenomenon between two liquid droplets for application in spray simulations. It has been known that the existing O’Rourke collision model widely used in CFD codes is inaccurate in determining collision outcomes and droplet behavior. In addition, since the collision probability of the model follows a statistical approach involving computational cell geometry, the prediction results should be strongly dependent on the cell size. As a result, to more accurately calculate droplet collisions, the technique for predicting the droplet velocity and its direction after collision must be extended for use in spray modeling. Further, it is also necessary to consider all the possible collision outcomes, such as bouncing, stretching separation, reflexive separation and coalescence. Therefore, this paper describes the appropriateness of a composite concept for modeling collision outcomes and the implementation of deterministic collision algorithms into a multidimensional CFD code for the calculation of post-collisional droplet movements. Furthermore, the existing model does not consider the formation of satellite droplets. For this reason, our present modeling concept includes a fragmenting droplet collision model. Using the present model, we have validated the collision interactions between liquid droplets under high Weber number conditions by comparing our calculations with experimental results from a binary droplet collision. This paper also deals with the application of the model to inter-impingement sprays by analyzing the atomization characteristics, such as mean droplet size and velocity, spray tip penetrations and spray-shapes of the impinging spray using the suggested collision algorithms and then comparing the results with available experimental data.     

Counter-current three-phase fluidization in a turbulent contact absorber: A CFD simulation

A computational fluid dynamics study of three-phase counter-current fluidization occurring in a turbulent contact absorber was performed. A two-dimensional, transient Eulerian multi-fluid model was used, in which the dispersed solid phase was modeled employing a kinetic theory of granular flow. The grid independence of the model, the effect of wall boundary conditions, the choice of granular temperature model, the effects of order of discretization scheme and drag models were studied for a base case setting. The results of simulations were validated against experimental results obtained from the literature. Once the model settings were finalized, simulations were performed for different gas and liquid velocities to predict the hydrodynamics of the absorber. Computed bed expansion and pressure drop were compared with experimental data. Good agreement between the two was found for low velocities of gas and liquid.     

海上溢油迁移转化的双层数学模型

海上溢油的双层数学模型将泄漏在海域中的溢油考虑为表层油膜和分布在整个水体中的悬浮 油滴两层组成. 采用Lagrangian追踪法模拟油膜的输移扩散. 在此方法中，油膜可用大量的小油滴表示，对每个油滴都规定随时间而变化的坐标系，油滴 的运动受到风、潮流和周围油的浓度影响. 油的迁移过程包括对流、扩展、湍动扩 散、附着在岸边以及沉降到海底等过程. 转化过程包括挥发、溶解、乳化等. 此外，光化学反应、生物降解能够改变油的特征以及减小油的污染. 该模型可以用于瞬时和连续溢油的情况，不仅可以用于模拟油，也可以模拟其它与油的密度 相近的有害物质的泄漏. 水动力学模型采用美国普林斯顿海洋模式(POM). 该溢油模型已应用于渤海海峡突发性溢油事故的模拟与预报.     

Spray analysis of a gasoline direct injector by means of two-phase PIV

The hollow-cone spray of a high-pressure swirl injector for a direct-injection spark-ignition (DISI) engine was investigated inside a pressure vessel by means of particle image velocimetry (PIV). As the interaction between the spray droplets and the ambient air is of particular interest for the mixture preparation process, two-phase PIV techniques were applied. To allow phase discrimination, fluorescent seeding particles were used to trace the gas phase. Because of the periodicity of piston engine injection, a statistical evaluation of ensemble-averaged fields to reduce cycle-to-cycle variations and to provide more general information about the two-phase flow was performed. Besides the general spray/air interaction process the investigation of the spray collapse at elevated ambient pressures was the main focus of the study. Future investigations of transient interaction processes require simultaneous techniques in combination with a high-speed camera to resolve the transient interaction phenomena. Therefore, optical filters that attenuate Mie-scattered light and transmit fluorescent light were used to collect both phases on the same image. Consequently, phase separation techniques were employed for data analysis. A masking and a peak separation technique are described and a comparison between the results of an instantaneous two-phase flow field in the spray cone of a DISI injector is presented in the paper.     

Hybrid turbulence simulation of spray impingement cooling: The effect of vortex motion on turbulent heat flux

Spray impingement has been a major interest of researchers in the areas of spray cooling, internal combustion, fire suppression and spray cooling, etc. for a long time. Numerous studies have been done in the area of spray cooling. Spray cooling with phase change takes advantage of relatively large amounts of latent heat and is capable of removing high heat fluxes from the surface, which has generated the interest of many researchers. In this paper, the turbulent characteristics of vapor formed during the spray impingement are studied. Water and gasoline are used in the numerical analysis of the two‐phase spray impingement on a heated wall. Hybrid turbulence modeling was used for the analysis where the subgrid scale model was employed away from the wall and k–ε model was used near the wall. Gasoline, at 298 K, was sprayed on the heated wall, kept constant at 650 K. The surrounding temperature was maintained at 400 K at the start of the simulation. In case of water and gasoline at Reynolds number 2750, the heated wall was kept constant at 400 K and the surrounding temperature was maintained at 298 K at the start of the simulations. The nozzle diameter of 100µm was used for this study, with the nozzle plate spacing ratio at 10. The spray was impinged on the flat plate at angles of 0, 15, and 30^°. Root mean‐squared velocities and turbulent heat flux were plotted in the water spray impingement for the different angles of impingement. The effect of turbulence on the heat transfer was observed. The effect of vortex motion on the turbulent heat flux values was analyzed using different Reynolds numbers of impingement and at different angles in case of gasoline. The turbulent heat flux attained the maximum values with high vortex formation. Upwash of fluid transported heat away from the wall, producing higher heat flux values in the region. Copyright © 2008 John Wiley & Sons, Ltd.     

Analysis of Mixing Models for use in Simulations of Turbulent Spray Combustion

The present study concerns the investigation of different mixing models for use in the transported probability density function (PDF) modeling of turbulent (reacting) spray flows. The modeling of the turbulent mixing and other characteristic scalar variables such as gas enthalpy using transported (joint) PDFs has become an important method to describe turbulent (reacting) spray flows since the evaporation process causes the PDF of the mixture fraction to deviate from the widely used β function, which is typically used in models for turbulent gas flows. In the PDF transport equation, the molecular mixing does not appear in closed form so that modeling strategies are required. For gas combustion, the interaction-by-exchange-with-the-mean (IEM) model, the modified Curl (MC) model, and the Euclidean minimum spanning tree (EMST) models are used. More recently, a new mixing model, the PSP model, which is based on parameterized scalar profiles has been developed. The present study focuses on the use and analysis of the IEM, MC and PSP models for turbulent spray flames. For this purpose, the models are reconsidered with respect to the evaporation process that must be included and evaluated if spray combustion is considered. For model evaluation, turbulent ethanol/air spray flames are simulated, and the results are compared to experimental data by A. Masri, University of Sydney, Australia.