Regular Article
Assessment of LESSTRIP approach for modeling of droplet dispersion in diesellike sprays
^{1}
Laboratoire de Mécanique des Fluides et d’Acoustique, École Centrale de Lyon, 36 Avenue Guy de Collongue, 69134 Écully Cedex, France
^{2}
Volvo Group Trucks Technology, 99 Route de Lyon, 69806 Saint Priest Cedex, France
^{*} Corresponding author: surya.oruganti@eclyon.fr
Received:
16
December
2018
Accepted:
11
April
2019
In this paper, the stochastic equations of droplet motion in turbulent flow, proposed recently by Gorokhovski and Zamansky (2018, Phys. Rev. Fluids 3, 3, 034602), are assessed for turbulent spray dispersion in diesel like conditions along with Large Eddy Simulation (LES) for the gaseous flow. For droplets above the Kolmogorov length scale, this model introduces the concept of the stochastic drag, independently of laminar viscosity. For droplets below the Kolmogorov length scale, the model equation does depend on the laminar viscosity through the Stokes drag but the particle motion is stochastically forced. Both the stochastic drag and the stochastic forcing of the Stokes drag equation are based on the simple lognormal stochastic process for the viscous dissipation (ϵ) “seen” along the droplet trajectory. In this paper, this model is applied in the framework of twoway coupling, wherein the turbulence generated by the spray inturn controls the spray dispersion. The criterion for the choice of one of the approaches, i.e., the stochastic drag or the stochastic forcing, follows the classical condition for drag coefficient based on the droplet Reynolds number (Re _{p}). The nonvaporizing spray experiments from Engine Combustion Network (ECN) are used as test cases. In addition to the comparison of the spray penetration length, spreading angle and spray structure with the experimental data, a qualitative analysis of the statistics of the droplet acceleration and gas phase velocities is presented. It was shown that the new approach is much more effective in modeling the spray dynamics on relatively coarser mesh. Consequently, the new approach in the framework of twoway coupling may predict the preferential concentration effects better, which is important for spray combustion.
© S. Kaundinya Oruganti et al., published by IFP Energies nouvelles, 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
With increasing restrictions on emissions, there has been a continuous evolution towards high pressure fuel injection systems, resulting in highly turbulent incylinder flows.
When the spray combustion takes place in highly turbulent conditions and when the chemical reactions are fast, it was shown in [1, 2] that the dispersion of vaporizing droplets is a strong factor controlling the combustion rate. When modeling this dispersion, the effect of turbulent structures of very different length scales on the droplet dynamics have to be considered. Direct Numerical Simulation (DNS), where the entire spectrum of turbulent scales are resolved, have been used to study the sprayturbulence interactions in simple flow configurations. Elghobashi et al. [3–5] have studied the mechanisms of twoway coupling in particle laden isotropic turbulence. They have characterized the modification of the turbulence by inertial particles depending on the magnitude of particle response time (τ _{p}) in comparison to Kolmogorov time scale (τ _{η}). Squires and Eaton [6] and Rouson and Eaton [7] have studied the preferential concentration of different particles i.e., accumulation of dense particles because of their inertia in regions of low vorticity and high strain. Toschi et al. [8] and Cencini et al. [9] have characterized the effects of intermittency at small scales on the acceleration statistics for inertial particles. They showed the presence of highly nonGaussian acceleration Probability Density Functions (PDF) with broad tails indicating high acceleration events. But DNS can be computationally very expensive and not feasible for real engine applications. Over the last decade, there has been a significant increase in application of Large Eddy Simulations (LES) in Internal Combustion Engines (ICE) because of its ability to directly resolve the large flow structures using spatially filtered equations. Usually, the small unresolved scales, often referred to as SubGrid Scales (SGS) are modeled using SGS turbulence models [10–13]. As stated by Rutland [14], it is a common practice in the context of Lagrangian modeling in LES to either use laminar correlations or extend submodels from RANS for modeling the interaction of SGS with the droplets. These approaches usually require very high grid resolutions. Hence Rutland [14], has emphasized on the need to develop SGS models specific to the context of LES modeling of liquid sprays.
To this end, Bharadwaj et al. [15, 16], have developed a model to account for spray induced turbulence, by reconstructing the SGS gas velocity at droplet positions using by deconvolution of the resolved velocity field. And the effect of unresolved scales on droplet dispersion is modeled decomposing the instantaneous gas velocity into the resolved velocity and the velocity at the SGS. The later is calculated from the subgrid kinetic energy (k _{sgs}) and assumed to be piecewise continuous in time [17]. Tsang et al. [18] further extended the model by decomposing the SGS velocity into deterministic and stochastic components, where the former is evaluated from deconvolution [15] and the later from the dispersion model [17]. The performance of the SGS model under diesel spray conditions has been evaluated. It was concluded that the dispersion only influences the spatial distribution of liquid droplets and does not have any effect on statistics of the carrier gas phase. Pozorski and Apte [19] have modeled the SGS velocity seen by the droplets using a Langevin stochastic equation to study preferential concentration in particle laden flows. This model has been successful in modeling the dynamics of only largeinertia particles but did not work well with small inertial particles which are concentrated in the small scales usually unresolved by LES. Bini and Jones [20, 21] have modeled the particle velocity increments using the Langevin equation with the variable intensity of the Wiener process. The model has been shown to capture the broad tails in the velocity increment distribution.
All the above proposed SGS models are based on reconstructing the fluctuating flow velocity along the particle trajectories from the subgrid kinetic energy (k _{sgs}). On the other hand, it was shown in the book of Kuznetzov and Sabel’nikov [22] that the statistics of the droplet velocity relative to the flow velocity is also controlled significantly by the mean viscous dissipation. The latter characterizes the smallscale dynamics in turbulence. This was also confirmed in DNS by Zamansky et al. [23, 24]. Consequently, in the droplet acceleration models, the instantaneous dissipation rate seen by the droplet along its trajectory is to be considered as an important parameter. However in LES this parameter is usually underresolved, thereby the droplet acceleration is underresolved as well. This motivated Gorokhovski and Zamansky [25, 26] to propose a new approach, in which the inertial particle acceleration on SGS is simulated by a stochastic process for viscous dissipation rate seen by the droplet along its trajectory. The parameters of those processes are defined from the locally resolved velocity field. Two SGS models are proposed, one for a particle above the Kolmogorov scale (referred to as “stochastic drag” model), and another one for the particle below the Kolmogorov scale as described in [25, 26]. This approach is assessed by DNS [27] for heavy particles immersed in the homogeneously sheared turbulence, and is referred to as LESSTRIP (Stochastic Response of Inertial Particle). Another parameter of this model is the stochastic direction of the droplet acceleration. The model for orientation vector of droplet acceleration is based on random walk over a unit sphere. While in [28], the stochastic equation for the orientation vector on the unit diffusion sphere is based on OrnsteinUhlenbeck (OU) process, this approach is retained in our current study.
In the previous works of Gorokhovski et al. [25–28] particle laden flows were studied with oneway coupling i.e., dynamics of particles do not influence the turbulence. But in direct injection engines, the momentum transfer from spray to the ambient gas flow contributes significantly to turbulence production. This motivated us to assess the LESSTRIP approach in diesellike conditions by applying twoway coupling between the spray and the surrounding gas flow. This paper is organized as follows: In Section 2, the computational framework i.e., the governing equations of both EulerianLagrangian phases and the standard dispersion model are described. In Section 3, the LESSTRIP model is recalled in the modified formulation. In Section 4, the experimental test cases and related numerical modeling issues are listed. In Section 5, validation of the dispersion model using the temporal evolution of spray structure is presented. Additionally, an analysis of the statistics of droplet and gas flow fields is provided to illustrate the capabilities of the new model to account for the twoway coupling of sprayturbulence interaction.
2 Governing equations
In this study the twophase problem is addressed by computing the carrier gas phase using LES and solving the individual droplet dynamics using the Lagrangian approach. The interaction between the two phases due to mass and momentum transfer are modeled by the source terms in the governing equations of the gas phase.
2.1 Gas phase equations – LES
In LES, the instantaneous velocity (u) of ambient gas phase is decomposed into filtered () and SGS components (u′):
(1)where, is solved by the continuity and momentum equation obtained by filtering the Navier Stokes equations of the instantaneous velocity field:
In equation (3), τ _{ij} is the viscous stress tensor, S _{liq} is the source term accounting for the interaction between the gas and liquid phases. On the other hand Γ_{ij}, is the subgrid shear stress tensor which is modeled using subgrid scale eddy viscosity (ν _{sgs}):
(4)where S _{ij} is the mean strain rate. In this study, the oneequation eddy viscosity model [9] is used, where SGS viscosity (ν _{sgs}) is calculated from the subgrid kinetic energy (k _{sgs}):
(5)where Δ is the filter width and k _{sgs} is the subgrid kinetic energy. The k _{sgs} is solved using a transport equation accounting for its production, dissipation and convection:
(6)where is the rate of dissipation at the resolved scales.
2.2 Liquid phase equations – Lagrangian approach
Usually the liquid phase is modeled as discrete parcels, where each parcel represents a ensemble of droplets with same properties. The dynamics of these parcels are obtained by individually solving the momentum equation of each parcel:
(9)where a _{p} is the droplet acceleration, u _{p} is the droplet velocity, τ _{p} is the droplet response time to fluid solicitations and C _{d} is the drag coefficient. The drag coefficient C _{d} is a function of droplet Reynolds number (). For (Re _{p} < 1000), it is assumed that the boundary layer around the droplet experiences a transition from laminar to turbulent and the drag coefficient is strongly dependent on Re:
For much higher velocities (Re _{p} > 1000), the boundary layer is highly chaotic with vortices of many different scales being shed from the droplet. In this regime the drag coefficient is roughly constant i.e., C _{d} = 0.42.
2.3 Source terms – twoway coupling
In the case of high pressure diesel sprays, the relative velocities of the two phases is significantly high (usually in the order of 100–400 m/s). So the impact of the momentum transfer from the liquid droplets on the dynamics of the carrier gas phase has to be considered. To this end, a source term S _{liq} is added to the momentum equation of the gas. In case of nonvaporizing sprays, it accounts only for the contribution of the drag force (F _{p} = m _{p}a _{p}) between the two phases. For a given computational cell, the net momentum transfer is the volume average of the contributions of all the droplets in the cell:
(11)where n _{dm} is number of droplets in a given parcel m, V _{cell} is volume of the computational cell and F _{dm} is the drag force acting on a single droplet in a given parcel m.
2.4 Standard SGS model for dispersion
Most commonly used SGS dispersion model for diesel spray studies [15, 16] is extended from RANS [17], where each component of SGS gas velocity is randomly sampled from a Gaussian distribution:
(12)where , k _{sgs} is interpreted as the subgrid kinetic energy. And the fluctuating SGS velocity u′_{p} is sampled once over a turbulent time scale, t _{turb}. The sampling time t _{turb} usually represents the time taken by the droplet to traverse through an eddy.
3 LESSTRIP approach for droplet acceleration
The main proposal in LESSTRIP approach [25, 26] is as follows. The acceleration of a finite size particle with size larger than the Kolmogorov length scale but smaller than the filter width is given by:
(14)where τ _{p,t}i.e., the droplet response time is a fluctuating variable due to the turbulence “seen” at the length scale of its diameter.
Using the Kolmogorov scaling for turbulent viscosity, , the droplet acceleration in equation (14) is characterized by a “fluctuating drag”. The statistics of the latter is conditioned by random viscous dissipation (ϵ) seen by a droplet along its trajectory:
Introducing the unit vector of the droplet acceleration direction,
(17)and from the Kolmogorov scaling, , the droplet acceleration in equation (16) can be expressed by the product of two independent random values i.e., the dissipation rate (ϵ) and the acceleration direction vector (). Each of these random variables is represented by a stochastic process in the LESSTRIP approach:
It is worthwhile to note that the experimental studies on the droplet acceleration statistics [29, 30] have shown that the variance of acceleration scales with , as it is the case in equation (18).
For small droplets with size below the Kolmogorov length scale, the LESSTRIP approach [25, 26] is based on the decomposition of the instantaneous acceleration of droplet into two parts i.e., response to resolved flow and the contribution on SGS:
(19)where the residual part is given by:
Equations (18)–(20) require stochastic equations for unit vector of the droplet acceleration direction () and for powers of viscous dissipation i.e., and . The OrnsteinUhlenbeck process for and the method of its integration is taken from [28]. While the stochastic equations for and are derived by Ito transformation of stochastic lognormal process [31] as described in [32].
(21)where dW(t) is the increment of standard Brownian process i.e., 〈dW(t)〉 = 0 and 〈dW(t)^{2}〉 = dt. In equation (21), ϵ _{sgs} represents the locally resolved dissipation rate, and the parameters and T _{χ} are chosen as functions of the mean eddyviscosity (ν _{t}), the Kolmogorov length scale η and the spatial filter width (Δ):
3.1 LESSTRIP model for diesel sprays
Because of the high injection velocities in diesel sprays, the flow in the nearnozzle region is strongly sheared. This results in very high levels of viscous dissipation along with very small magnitudes of local Kolmogorov’s dissipative scales. At the same time, the liquid is not yet well atomized, with liquid elements typically larger than the local Kolmogorov length scale. Equation (18) assumes the drag to be independent of the laminar viscosity. It is coherent with the classical observation for the sphere drag if Re > 1000. Therefore the condition for equation (18) is Re > 1000. Finally, the reformulated LESSTRIP model for diesel spray applications is given by:
For Re _{p} > 1000 and Δ > d _{p} > η
Otherwise for Re _{p} < 1000
4 Experimental and numerical setup
4.1 Experimental data
The nonvaporizing spray data from Engine Combustion Network, ECN [33, 34] is used . Three different test cases with dodecane (C_{12}H_{26}) as surrogate of the diesel fuel are used in this study. The parameters for different test cases are listed in Table 1. Apart from the comparison of the spray penetration length and spreading angle, the schlieren images of the spray evolution for test case3, as shown in Figure 1 are used to validate the LESSTRIP model.
Fig. 1 Schlieren images of spray structure at different times for test case3 (approx. 0.5 ms and 1.0 ms). 
Nonvaporizing spray test conditions.
4.2 Numerical models
The simulations are performed in OpenFOAM, using a twoway coupled EulerianLagrangian approach with exchange of mass and momentum between the phases. In Lagrangian spray simulations, the dynamics of the droplets representing the liquid fuel are tracked in a fixed Eulerian framework. A constant volume, cylindrical vessel with a radius of 25 mm and length of 100 mm is used in the simulations as shown in Figure 2. The vessel has a rectangular cross section (region 1) of 10 mm close to the injector to discretize the mesh with a fine grid to better resolve the spray dynamics. Initially a uniform grid size of 0.25 mm (referred to as Grid A) is used in both the regions to discretize the geometry. An additional fine mesh with grid size of 0.125 mm in region1 and a nonuniform grid with cell size geometrically expanding from 0.125 mm to 0.25 mm in region2 (referred to as Grid B) is simulated for studying the mesh sensitivity. The parcels are injected using a solidcone injection method, where the parcels are sampled randomly from the userdefined initial size distribution. In this study RosinRammler distribution model is used to represent the initial size distribution of injected parcels. The parcel injection velocity is calculated from the mass flow rate profile and its initial direction is randomly sampled to lie within a userspecified cone angle. Because of the high injection velocities, the secondary breakup occurs close to the injector and thereby the choice of an appropriate breakup model significantly effects the predictions of the spray evolution. In several studies [35–37], KHRT breakup model [38] has indicated its superior performance in predicting the spray structure under high pressure diesel conditions. The KHRT model includes two modes of breakup: KH breakup, accounts for the growth of unstable waves on the liquid jet due to differences in velocity between the ambient gas and the liquid; and RT breakup, accounts for growth of waves on the droplet’s surface due to aerodynamic drag at the dropletgas interface. An Euler scheme for temporal discretion and a second order scheme for spatial discretization are used in the simulation. The injection and breakup model constants have been selected based on [39, 40], where influence of these parameters on the spray results have been reported in detail.
Fig. 2 Computational geometry representing a constant volume spray chamber. 
5 Results
5.1 Global spray characteristics
In order to evaluate the dispersion models, a quantitative comparison of the liquid penetration lengths, spreading angles for the different experimental conditions is presented in this section. In this study, the liquid penetration length is defined as the distance where the accumulated liquid droplet mass reaches 95% of the total liquid mass injected at any given instance of time. A comparison of liquid penetration length predicted for different test cases using Grid A is shown in Figures 3–5. The experimental penetration lengths presented in this study have an error of about 1–2%.
Fig. 3 Temporal evolution of liquid penetration length – test case1 with Grid A. 
Fig. 4 Temporal evolution of liquid penetration length – test case2 with Grid A. 
Fig. 5 Temporal evolution of liquid penetration length – test case3 with Grid A. 
Initial penetration is strongly controlled by high momentum and secondary breakup of liquid droplets. So the initial underprediction of the penetration lengths by both models in Figures 3–5 are overlooked in this study. But further downstream (i.e., after an axial distance of 20 mm away from injector), the LESSTRIP model gives more accurate prediction of the evolution of penetration length. Also it is clearly evident from Figures 3–5, that as the injection pressure increases the deviation in the penetration length predicted by the standard SGS model from the experiment increases significantly. Figures 6 and 7 present the grid density analysis of the SGS models for test case3 using two different filter sizes i.e., Grid A and B.
Fig. 6 Grid sensitivity of the standard SGS model for test case3. 
Fig. 7 Grid sensitivity of the LESSTRIP model for test case3. 
Figure 6 shows that the standard SGS model requires much finer grid resolution to capture the evolution of penetration length. As grid resolution is increased, more and more scales are directly resolved by LES through . Therefore, the effect of SGS velocity component ( is negligible, which clearly implies that the effect of unresolved scales is not well accounted by the model. On the other hand, from Figure 7 it can be seen that the STRIP model captures the evolution of penetration length even on coarse grid size. This is because, the formulation of STRIP model is based on modeling the droplet acceleration in terms of the ϵ which is a characteristic of the inertial scales. The model parameters of the stochastic equation for ϵ in equation (22) account for the effects of grid size on the evolution of ϵ, thereby making the dispersion model less sensitive to the grid size. The variation in the penetration length for different grid sizes as shown in Figure 7, is attributed primarily to the grid sensitivity of the KHRT breakup model. In addition to the penetration lengths, the spreading angle is a most commonly used parameter to quantify spray dispersion. The spreading angle is calculated from the isosceles triangle built considering 75% of the spray penetration, similar to [41]. For a much simplified analysis, an average spreading angle is calculated by averaging only the angles obtained when the penetration was greater than 35 mm. For all the simulations, the droplet parcels are injected within a cone angle of 20. From Figure 8 it can be seen that the LESSTRIP model consistently predicts much larger spreading angle (~4–5) than the standard SGS model for all the experimental test conditions.
Fig. 8 Comparison of the time averaged spray spreading angle for different test conditions using Grid A. 
Finally a direct comparison of the spray structure with the experimental schlieren images of test case3 is presented in Figures 9–11.
Fig. 9 Experimental schlieren images of spray for test case3 at three different time instants. (a) t = 0.5 ms, (b) t = 1.0 ms, (c) t = 1.5 ms. 
Fig. 10 Instantaneous spray structure from LESSTRIP using Grid A at three different time instants. (a) t = 0.5 ms, (b) t = 1.0 ms, (c) t = 1.5 ms. 
Fig. 11 Instantaneous spray structure from standard SGS using Grid A at three different time instants (a) t = 0.5 ms, (b) t = 1.0 ms, (c) t = 1.5 ms. 
Figure 9 shows the schlieren images of spray structure for test case3 at different time instances. Figure 10 shows the spray structure for the corresponding time instances obtained using LESSTRIP model on Grid A. Even though dispersion is not evident close to the injector (upto 20 mm), overall spray structure is in good agreement with the experiment. On the other hand, Figure 11 represents the spray structure predicted by standard SGS model on Grid A. It can be noticed that the standard SGS model predicts a highly penetrative spray core surrounded by large number of parcels dispersed very close to the injector.
5.2 Droplet acceleration statistics
In order to explain the underlying physical differences in the results predicted by the two approaches, a more detailed analysis of the droplet acceleration statistics are presented in this section. All the results presented in this section are for the test case3 using Grid A. In recent experiments in homogeneous statistically stationary turbulence [42–47], it was shown that the norm of acceleration is correlated on large times comparable with the integral time, while the direction is correlated on short times of order of the Kolmogorov time (τ _{η}). In this way the intermittency is manifested i.e., the vortical filaments of small scales are as much energetic as the largescale turbulent structures. So a comparison of the autocorrelation function ( of the droplet acceleration (a _{p}) for the two dispersion models is presented in Figures 12–13. The autocorrelation function of the droplet acceleration is calculated using equation (25):
(25)where k = 1, 2 represents the axial and the radial components of acceleration respectively. The correlation time (τ) is normalized by the the Kolmogorov time scale (τ _{η}).
Fig. 12 Comparison of autocorrelation of the droplet acceleration norm, axial and radial components for standard SGS model. 
Fig. 13 Comparison of autocorrelation of the droplet acceleration norm, axial and radial components for LESSTRIP model. 
In the case of standard SGS model (Fig. 12), the axial component of droplet acceleration vector has much longer correlation (~8τ _{η}) than its radial component and norm. This shows that the droplets retain their axial acceleration component for a longer duration implying less radial dispersion. This is the consequence of the model employed in equation (12), where the norm of the acceleration obtained from sampling of the turbulent velocity (u′) is noncorrelated in time. At the same time the direction of the droplet acceleration is retained on times defined by large turbulent structures. Therefore, a large number of droplets in the spray core only penetrate axially without any radial dispersion. Moreover, only the secondary droplets from the KHRT breakup which acquire an additional radial velocity component are dispersed radially as shown in Figure 11. On the other hand, in the case of LESSTRIP model (Fig. 13) both the axial and radial components of droplet acceleration correlate with the Kolmogorov time scales (τ _{η}), whereas the norm has much longer time correlation (5τ _{η}). The shorter time correlations of the components is because of the stochastic model, where the frequency of the fluctuations in orientation vector of droplet acceleration are scaled in the order of Kolmogorov time scales. Therefore, the LESSTRIP model is more effective in modeling the effects of intermittency at small scales on droplet dispersion.
A comparison of the PDF of the droplet acceleration at two different time instances i.e., t = 1.0 and 1.5 ms is shown in Figures 14 and 15. With increasing injection pressure, the probability of the high acceleration events increases. But as seen from Figures 14 and 15, the standard SGS model usually underresolves the droplet acceleration (a _{p}). On the other hand, the acceleration PDF of LESSTRIP model has much broader tails indicating better resolution of the large acceleration events resulting from interaction of droplet with inertial scales.
Fig. 14 Comparison of droplet acceleration PDF at t = 1.0 ms. The acceleration a _{p} is normalized by the root mean square value a _{rms}. 
Fig. 15 Comparison of droplet acceleration PDF at t = 1.5 ms. The acceleration a _{p} is normalized by the root mean square value a _{rms.} 
5.3 Impact on the gas phase dynamics
Finally, the impact of the SGS dispersion models on resolved gas phase dynamics is investigated. A qualitative comparison is made by plotting isosurfaces of the Qcriterion to visualize the vortex structures. A positive value of Q representing the local rotational motion of the fluid was chosen for plotting the isosurfaces to visualize coherent vortices. Figure 16 shows the isosurfaces of the Q criterion predicted by the different SGS models. Similar to the notable difference of the liquid spray structures as shown in Figures 7 and 8 the structure of the gas jets predicted by the SGS models are clearly different.
Fig. 16 Comparison of isosurface of Qcriterion Q = 1.0e + 08 s^{−2} colored by the gas velocities at time t = 1.5 ms. (a) Standard model, (b) LESSTRIP model. 
Next a comparison of radial profiles of the gas phase velocity at two different crosssections i.e., 20 mm and 30 mm downstream are presented in Figures 17 and 18. The results represent an ensemble average of five different realizations of each model. The average velocity profiles from simulations without any dispersion model and standard SGS model are of the same order of magnitude. On the other hand the average velocity profiles from the LESSTRIP simulations differ by an order of magnitude. This qualitatively implies that in the case of LESSTRIP approach the spray dynamics and gasflow field are strongly coupled. Since, LESSTRIP model directly models the droplet drag force, the effect of spray dispersion on the ambient gas dynamics is implicitly modeled through the source term (S _{liq}) represented by equation (11).
Fig. 17 Ensemble average gas velocity profile at axial distance of 20 mm from the injector. 
Fig. 18 Ensemble average gas velocity profile at axial distance of 30 mm from the injector. 
6 Conclusion
In this paper, the recently proposed stochastic equations (i.e., LESSTRIP) for the droplet motion in turbulent flows are used to model spray dispersion in diesel like conditions. The model uses the concept of stochastic drag for droplets larger than Kolmogorov scale and large droplet Reynolds numbers. On the other hand, for droplets smaller than Kolmogorov length scales the concept of stochastic forcing on droplet motion is employed. A detailed analysis of the performance of LESSTRIP in comparison with the most commonly used SGS dispersion model for three different nonvaporizing spray conditions is performed. In comparison with the standard dispersion model, LESSTRIP model gave good liquid spray penetration predictions for varying parameters such as injection pressure, and nozzle size. The gridsensitivity analysis has shown that the LESSTRIP model is less sensitive to the grid size and can resolve the spray dynamics even on relatively coarse grid. In order to quantify spray dispersion, a comparison of the time averaged spreading angles for different test cases is presented. In comparison with the standard dispersion model, LESSTRIP has consistently predicted larger spreading angles. The higher dispersion in case of LESSTRIP is explained using the temporal correlation of the droplet acceleration. The shorter correlation of droplet acceleration components represent large fluctuations in orientation of droplet acceleration due to intermittency effects of small scales. Accounting for the intermittency effects by LESSTRIP results in larger spray dispersion. On the other hand, the higher spray penetration with standard dispersion model is attributed to longer correlation of the axial component of the droplet acceleration. Finally, to study the impact of the dispersion on twoway coupling of sprayturbulence interaction, a comparison of the gas flow velocities is presented. It was found that, while the results predicted without any dispersion model and with standard model are of same order of magnitude, the results of LESSTRIP differ by an order of magnitude. This qualitatively shows that the LESSTRIP model is more effective in modeling the twoway coupling of sprayturbulence interaction.
Acknowledgments
The authors would like to thank the Association Nationale de la Recherche et de la Technologie (ANRT) and the Volvo group for sponsoring the project. All the simulations in this study are performed on internal clusters of Volvo group and PMCS2I cluster at Ecole Centrale de Lyon. The authors would like to express their gratitude to Engine combustion network and Professor Raul Payri from CMT – Motores Térmicos, Universitat Politècnica de València for allowing access to some of their experimental data base.
References
 Durand P., Gorokhovski M., Borghi R. (1999) An application of the probability density function model to diesel engine combustion, Combust. Sci. Technol. 144, 1–6, 47–78. [Google Scholar]
 Sabel'nikov V., Gorokhovski M., Baricault N. (2006) The extended iem mixing model in the framework of the composition pdf approach: Applications to diesel spray combustion, Combust. Theory Modell. 10, 1, 155–169. [CrossRef] [Google Scholar]
 Elghobashi S., Truesdell G.C. (1993) On the twoway interaction between homogeneous turbulence and dispersed solid particles. I: Turbulence modification, Phys. Fluids A: Fluid Dyn. 5, 7, 1790–1801. [CrossRef] [Google Scholar]
 Ferrante A., Elghobashi S. (2003) On the physical mechanisms of twoway coupling in particleladen isotropic turbulence, Phys. Fluids 15, 2, 315–329. [CrossRef] [Google Scholar]
 Truesdell G.C., Elghobashi S. (1994) On the twoway interaction between homogeneous turbulence and dispersed solid particles. II. Particle dispersion, Phys. Fluids 6, 3, 1405–1407. [CrossRef] [Google Scholar]
 Squires K.D., Eaton J.K. (1991) Preferential concentration of particles by turbulence, Phys. Fluids A: Fluid Dyn. 3, 5, 1169–1178. [NASA ADS] [CrossRef] [Google Scholar]
 Rouson D.W., Eaton J.K. (2001) On the preferential concentration of solid particles in turbulent channel flow, J. Fluid Mech. 428, 149–169. [Google Scholar]
 Toschi F., Bec J., Biferale L., Boffetta G., Celani A., Cencini M., Lanotte A.S., Musacchio S. (2008) Acceleration statistics of inertial particles from high resolution DNS turbulence, in: IUTAM Symposium on Computational Physics and New Perspectives in Turbulence, Springer, Dordrecht, pp. 73–78. [Google Scholar]
 Cencini M., Bec J., Biferale L., Boffetta G., Celani A., Lanotte A.S., Toschi F. (2006) Dynamics and statistics of heavy particles in turbulent flows, J. Turbul. 7, N36. [CrossRef] [Google Scholar]
 Smagorinsky J. (1963) General circulation experiments with the primitive equations: I. The basic experiment, Mon. Weather Rev. 91, 3, 99–164. [Google Scholar]
 Yoshizawa A., Horiuti K. (1985) A statisticallyderived subgridscale kinetic energy model for the largeeddy simulation of turbulent flows, J. Phys. Soc. Jpn. 54, 8, 2834–2839. [CrossRef] [Google Scholar]
 Germano M., Piomelli U., Moin P., Cabot W.H. (1991) A dynamic subgridscale eddy viscosity model, Phys. Fluids A: Fluid Dyn. 3, 7, 1760–1765. [NASA ADS] [CrossRef] [Google Scholar]
 Chumakov S.G., Rutland C.J. (2005) Dynamic structure subgridscale models for large eddy simulation, Int. J. Num. Methods Fluids 47, 8–9, 911–923. [CrossRef] [Google Scholar]
 Rutland C.J. (2011) Largeeddy simulations for internal combustion engines – A review, Int. J. Engine Res. 12, 5, 421–451. [CrossRef] [EDP Sciences] [Google Scholar]
 Bharadwaj N., Rutland C.J., Chang S.M. (2009) Large eddy simulation modelling of sprayinduced turbulence effects, Int. J. Engine Res. 10, 2, 97–119. [CrossRef] [Google Scholar]
 Bharadwaj N., Rutland C.J. (2010) A largeeddy simulation study of subgrid twophase interaction in particleladen flows and diesel engine sprays, Atomization Sprays 20, 8. [CrossRef] [Google Scholar]
 Amsden A.A., Butler T.D., O’Rourke P.J. (1987) The KIVAII computer program for transient multidimensional chemically reactive flows with sprays, SAE Trans. 373–383. [Google Scholar]
 Tsang C.W., Kuo C.W., Trujillo M., Rutland C. (2018) Evaluation and validation of largeeddy simulation subgrid spray dispersion models using highfidelity volumeoffluid simulation data and engine combustion network experimental data, Int. J. Engine Res., 1468087418772219. [Google Scholar]
 Pozorski J., Apte S.V. (2009) Filtered particle tracking in isotropic turbulence and stochastic modeling of subgridscale dispersion, Int. J. Multiphase Flow 35, 2, 118–128. [CrossRef] [Google Scholar]
 Bini M., Jones W.P. (2007) Particle acceleration in turbulent flows: A class of nonlinear stochastic models for intermittency and heavy tailed pdfs, Phys. Fluids 19, 3, 035104. [CrossRef] [Google Scholar]
 Bini M., Jones W.P. (2008) Large eddy simulation of particle laden turbulent flows, J. Fluid Mech. 614, 207–252. [Google Scholar]
 Kuznetsov V.R., Sabel’nikov V.A. (1990) Turbulence and combustion, Hemispher Publishing Coorporation, Washington, DC. [Google Scholar]
 Zamansky R., Vinkovic I., Gorokhovski M. (2011) Acceleration statistics of solid particles in turbulent channel flow, Phys. Fluids 23, 11, 113304. [CrossRef] [Google Scholar]
 Zamansky R., Vinkovic I., Gorokhovski M. (2013) Acceleration in turbulent channel flow: Universalities in statistics, subgrid stochastic models and an application, J. Fluid Mech. 721, 627–668. [Google Scholar]
 Gorokhovski M., Zamansky R. (2018) Modeling the effects of small turbulent scales on the drag force for particles below and above the Kolmogorov scale, Phys. Rev. Fluids 3, 3, 034602. [Google Scholar]
 Gorokhovski M., Zamansky R. (2014) Lagrangian simulation of large and small inertial particles in a high Reynolds number flow: Stochastic simulation of subgrid turbulence/particle interactions, in: Center for Turbulence Research, Proceedings of the Summer Program, pp. 37–46. [Google Scholar]
 Barge A., Gorokhovski M. (2019) Accelerations in stationary turbulence under homogeneous shear: DNS, subsequent subgrid stochastic models and application for inertial particle dynamics, J. Fluid Mech. (in preparation). [Google Scholar]
 Sabelnikov V., Barge A., Gorokhovski M. (2019) Stochastic modeling of fluid acceleration on residual scales and dynamics of suspended inertial particles in turbulence, Phys. Rev. Fluids 4, 4, 044301. [Google Scholar]
 Qureshi N.M., Arrieta U., Baudet C., Cartellier A., Gagne Y., Bourgoin M. (2008) Acceleration statistics of inertial particles in turbulent flow, Eur. Phys. J. B 66, 4, 531–536. [EDP Sciences] [Google Scholar]
 Qureshi N.M., Bourgoin M., Baudet C., Cartellier A., Gagne Y. (2007) Turbulent transport of material particles: An experimental study of finite size effects, Phys. Rev. Lett. 99, 18, 184502. [CrossRef] [PubMed] [Google Scholar]
 Pope S.B., Chen Y.L. (1990) The velocitydissipation probability density function model for turbulent flows, Phys. Fluids A: Fluid Dyn. 2, 8, 1437–1449. [CrossRef] [Google Scholar]
 Sabelnikov V., ChtabDesportes A., Gorokhovski M. (2011) New subgrid stochastic acceleration model in LES of highReynoldsnumber flows, Eur. Phys. J. B 80, 177. [EDP Sciences] [Google Scholar]
 Pickett L.M., Genzale C.L., Bruneaux G., Malbec L.M., Hermant L., Christiansen C., Schramm J. (2010) Comparison of diesel spray combustion in different hightemperature, highpressure facilities, SAE Int. J. Engines 3, 2, 156–181. [CrossRef] [Google Scholar]
 Payri R., Viera J.P., Gopalakrishnan V., Szymkowicz P.G. (2016) The effect of nozzle geometry over internal flow and spray formation for three different fuels, Fuel 183, 20–33. [CrossRef] [Google Scholar]
 Kitaguchi K., Fujii T., Hatori S., Hori T., Senda J. (2014) Effect of breakup model on largeeddy simulation of diesel spray evolution under high back pressures, Int. J. Engine Res. 15, 5, 522–538. [CrossRef] [Google Scholar]
 Fujimoto H., Tsukasa H.O.R.I., Senda J. (2009) Effect of breakup model on diesel spray structure simulated by large eddy simulation (No. 2009240024). SAE Technical Paper. [Google Scholar]
 Wehrfritz A., Vuorinen V., Kaario O., Larmi M. (2013) Large eddy simulation of highvelocity fuel sprays: Studying mesh resolution and breakup model effects for spray A, Atomization Sprays 23, 5, 419–442. [CrossRef] [Google Scholar]
 Patterson M.A., Reitz R.D. (1998) Modeling the effects of fuel spray characteristics on diesel engine combustion and emission (No. 980131). SAE Technical Paper. [Google Scholar]
 Hori T., Kuge T., Senda J., Fujimoto H. (2008) Effect of convective schemes on LES of fuel spray by use of KIVALES (No. 2008010930). SAE Technical Paper. [Google Scholar]
 Vogiatzaki K., Crua C., Morgan R., Heikal M. (2017) A study of the controlling parameters of fuel air mixture formation for ECN Spray A. [Google Scholar]
 Naber J.D., Siebers D. (1996) Effects of gas density and vaporization on penetration and dispersion of Diesel sprays, SAE Technical Paper 960034 105, 412, 82–111. [Google Scholar]
 Mordant N., Delour J., Léveque E., Arnéodo A., Pinton J.F. (2002) Long time correlations in Lagrangian dynamics: A key to intermittency in turbulence, Phys. Rev. Lett. 89, 25, 254502. [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Mordant N., Crawford A.M., Bodenschatz E. (2004) Threedimensional structure of the Lagrangian acceleration in turbulent flows, Phys. Rev. Lett. 93, 21, 214501. [CrossRef] [PubMed] [Google Scholar]
 Voth G.A., La Porta A., Crawford A.M., Alexander J., Bodenschatz E. (2002) Measurement of particle accelerations in fully developed turbulence, J. Fluid Mech. 469, 121–160. [Google Scholar]
 Yeung P.K., Pope S.B. (1989) Lagrangian statistics from direct numerical simulations of isotropic turbulence, J. Fluid Mech. 207, 531–586. [Google Scholar]
 Yeung P.K., Pope S.B., Lamorgese A.G., Donzis D.A. (2006) Acceleration and dissipation statistics of numerically simulated isotropic turbulence, Phys. Fluids 18, 065103. [CrossRef] [Google Scholar]
 Toschi F., Bodenschatz E. (2009) Lagrangian properties of particles in turbulence, Ann. Rev. Fluid Mech. 41, 1, 375–404. [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1 Schlieren images of spray structure at different times for test case3 (approx. 0.5 ms and 1.0 ms). 

In the text 
Fig. 2 Computational geometry representing a constant volume spray chamber. 

In the text 
Fig. 3 Temporal evolution of liquid penetration length – test case1 with Grid A. 

In the text 
Fig. 4 Temporal evolution of liquid penetration length – test case2 with Grid A. 

In the text 
Fig. 5 Temporal evolution of liquid penetration length – test case3 with Grid A. 

In the text 
Fig. 6 Grid sensitivity of the standard SGS model for test case3. 

In the text 
Fig. 7 Grid sensitivity of the LESSTRIP model for test case3. 

In the text 
Fig. 8 Comparison of the time averaged spray spreading angle for different test conditions using Grid A. 

In the text 
Fig. 9 Experimental schlieren images of spray for test case3 at three different time instants. (a) t = 0.5 ms, (b) t = 1.0 ms, (c) t = 1.5 ms. 

In the text 
Fig. 10 Instantaneous spray structure from LESSTRIP using Grid A at three different time instants. (a) t = 0.5 ms, (b) t = 1.0 ms, (c) t = 1.5 ms. 

In the text 
Fig. 11 Instantaneous spray structure from standard SGS using Grid A at three different time instants (a) t = 0.5 ms, (b) t = 1.0 ms, (c) t = 1.5 ms. 

In the text 
Fig. 12 Comparison of autocorrelation of the droplet acceleration norm, axial and radial components for standard SGS model. 

In the text 
Fig. 13 Comparison of autocorrelation of the droplet acceleration norm, axial and radial components for LESSTRIP model. 

In the text 
Fig. 14 Comparison of droplet acceleration PDF at t = 1.0 ms. The acceleration a _{p} is normalized by the root mean square value a _{rms}. 

In the text 
Fig. 15 Comparison of droplet acceleration PDF at t = 1.5 ms. The acceleration a _{p} is normalized by the root mean square value a _{rms.} 

In the text 
Fig. 16 Comparison of isosurface of Qcriterion Q = 1.0e + 08 s^{−2} colored by the gas velocities at time t = 1.5 ms. (a) Standard model, (b) LESSTRIP model. 

In the text 
Fig. 17 Ensemble average gas velocity profile at axial distance of 20 mm from the injector. 

In the text 
Fig. 18 Ensemble average gas velocity profile at axial distance of 30 mm from the injector. 

In the text 