Large-Eddy Simulation of Atomizing Spray With Stochastic Modeling of Secondary Breakup

Author(s):  
Sourabh V. Apte ◽  
Mikhael Gorokhovski ◽  
Parviz Moin

Large-eddy simulation (LES) of reacting multi-phase flows in practical combustor geometries is essential to accurately predict complex physical phenomena of turbulent mixing and combustion dynamics. This necessitates use of Lagrangian particle-tracking methodology for liquid phase in order to correctly capture the droplet evaporation rates in the sparse-liquid regime away from the fuel injector. Our goal in the present work is to develop a spray-atomization methodology which can be used in conjuction with the standard particle-tracking schemes and predict the droplet-size distribution accurately. The intricate process of primary atomization and lack of detailed experimental observations close to the injector requires us to model its global effects and focus on secondary breakup to capture the evolution of droplet sizes. Accordingly, a stochastic model for LES of atomizing spray is developed. Following Kolmogorov’s idea of viewing solid particle-breakup as a discrete random process, atomization of liquid blobs at high relative liquid-to-gas velocity is considered in the framework of uncorrelated breakup events, independent of the initial droplet size. Kolmogorov’s discrete model of breakup is represented by Fokker-Planck equation for the temporal and spatial evolution of droplet radius distribution. The parameters of the model are obtained dynamically by relating them to the local Weber number. A novel hybrid-approach involving tracking of individual droplets and a group of like-droplets known as parcels is developed to reduce the computational cost and maintain the essential features and dynamics of spray evolution. The present approach is shown to capture the complex fragmentary process of liquid atomization in idealized and realistic Diesel and gas-turbine combustors.

Author(s):  
S. Gallot-Lavallée ◽  
W. P. Jones ◽  
A. J. Marquis

AbstractA computational investigation of three configurations of the Delft Spray in Hot-diluted Co-flow (DSHC) is presented. The selected burner comprises a hollow cone pressure swirl atomiser, injecting an ethanol spray, located in the centre of a hot co-flow generator, with the conditions studied corresponding to Moderate or Intense Low-oxygen Dilution (MILD) combustion. The simulations are performed in the context of Large Eddy Simulation (LES) in combination with a transport equation for the joint probability density function (pdf) of the scalars, solved using the Eulerian stochastic field method. The liquid phase is simulated by the use of a Lagrangian point particle approach, where the sub-grid-scale interactions are modelled with a stochastic approach. Droplet breakup is represented by a simple primary breakup model in combination with a stochastic secondary breakup formulation. The approach requires only a minimal knowledge of the fuel injector and avoids the need to specify droplet size and velocity distributions at the injection point. The method produces satisfactory agreement with the experimental data and the velocity fields of the gas and liquid phase both averaged and ‘size-class by size-class’ are well depicted. Two widely accepted evaporation models, utilising a phase equilibrium assumption, are used to investigate the influence of evaporation on the evolution of the liquid phase and the effects on the flame. An analysis on the dynamics of stabilisation sheds light on the importance of droplet size in the three spray flames; different size droplets play different roles in the stabilisation of the flames.


2017 ◽  
Author(s):  
Mikko Auvinen ◽  
Leena Järvi ◽  
Antti Hellsten ◽  
Üllar Rannik ◽  
Timo Vesala

Abstract. Conventional footprint models cannot account for the heterogeneity of the urban landscape imposing a pronounced uncertainty on the spatial interpretation of eddy-covariance (EC) flux measurements in urban studies. This work introduces a computational methodology that enables the generation of detailed footprints in arbitrarily complex urban flux measurements sites. The methodology is based on conducting high-resolution large-eddy simulation (LES) and Lagrangian stochastic (LS) particle analysis on a model that features a detailed topographic description of a real urban environment. The approach utilizes an arbitrarily sized target volume set around the sensor in the LES domain, to collect a dataset of LS particles which are seeded from the potential source-area of the measurement and captured at the sensor site. The urban footprint is generated from this dataset through a piecewise post-processing procedure, which divides the footprint evaluation into multiple independent processes that each yield an intermediate result that are ultimately selectively combined to produce the final footprint. The strategy reduces the computational cost of the LES-LS simulation and incorporates techniques to account for the complications that arise when the EC sensor is mounted on a building instead of a conventional flux tower. The presented computational framework also introduces a result assessment strategy which utilizes the obtained urban footprint together with a detailed land cover type dataset to estimate the potential error that may arise if analytically derived footprint models were employed instead. The methodology is demonstrated with a case study that concentrates on generating the footprint for a building-mounted EC measurement station in downtown Helsinki, Finland, under neutrally stratified atmospheric boundary layer.


Author(s):  
Y. See ◽  
M. Wang ◽  
J. Bohbot ◽  
O. Colin

Abstract The Species-Based Extended Coherent Flamelet Model (SB-ECFM) was developed and previously validated for 3D Reynolds-Averaged Navier-Stokes (RANS) modeling of a spark-ignited gasoline direct injection engine. In this work, we seek to extend the SB-ECFM model to the large eddy simulation (LES) framework and validate the model in a homogeneous charge spark-ignited engine. In the SB-ECFM, which is a recently developed improvement of the ECFM, the progress variable is defined as a function of real species instead of tracer species. This adjustment alleviates discrepancies that may arise when the numerical treatment of real species is different than that of the tracer species. Furthermore, the species-based formulation also allows for the use of second-order numeric, which can be necessary in LES cases. The transparent combustion chamber (TCC) engine is the configuration used here for validating the SB-ECFM. It has been extensively characterized with detailed experimental measurements and the data are widely available for model benchmarking. Moreover, several of the boundary conditions leading to the engine are also measured experimentally. These measurements are used in the corresponding computational setup of LES calculations with SB-ECFM. Since the engine is spark ignited, the Imposed Stretch Spark Ignition Model (ISSIM) is utilized to model this physical process. The mesh for the current study is based on a configuration that has been validated in a previous LES study of the corresponding motored setup of the TCC engine. However, this mesh was constructed without considering the additional cells needed to sufficiently resolve the flame for the fired case. Thus, it is enhanced with value-based Adaptive Mesh Refinement (AMR) on the progress variable to better capture the flame front in the fired case. As one facet of model validation, the ensemble average of the measured cylinder pressure is compared against the LES/SB-ECFM prediction. Secondly, the predicted cycle-to-cycle variation by LES is compared with the variation measured in the experimental setup. To this end, the LES computation is required to span a sufficient number of engine cycles to provide statistical convergence to evaluate the coefficient of variation (COV) in peak cylinder pressure. Due to the higher computational cost of LES, the runtime required to compute a sufficient number of engine cycles sequentially can be intractable. The concurrent perturbation method (CPM) is deployed in this study to obtain the required number of cycles in a reasonable time frame. Lastly, previous numerical and experimental analyses of the TCC engine have shown that the flow dynamics at the time of ignition is correlated with the cycle-to-cycle variability. Hence, similar analysis is performed on the current simulation results to determine if this correlation effect is well-captured by the current modeling approach.


2003 ◽  
Vol 19 (5) ◽  
pp. 782-794 ◽  
Author(s):  
Ying Huang ◽  
Hong-Gye Sung ◽  
Shih-Yang Hsieh ◽  
Vigor Yang

2013 ◽  
Vol 136 (6) ◽  
Author(s):  
James Tyacke ◽  
Paul Tucker ◽  
Richard Jefferson-Loveday ◽  
Nagabushana Rao Vadlamani ◽  
Robert Watson ◽  
...  

Flows throughout different zones of turbines have been investigated using large eddy simulation (LES) and hybrid Reynolds-averaged Navier–Stokes-LES (RANS-LES) methods and contrasted with RANS modeling, which is more typically used in the design environment. The studied cases include low and high-pressure turbine cascades, real surface roughness effects, internal cooling ducts, trailing edge cut-backs, and labyrinth and rim seals. Evidence is presented that shows that LES and hybrid RANS-LES produces higher quality data than RANS/URANS for a wide range of flows. The higher level of physics that is resolved allows for greater flow physics insight, which is valuable for improving designs and refining lower order models. Turbine zones are categorized by flow type to assist in choosing the appropriate eddy resolving method and to estimate the computational cost.


2010 ◽  
Vol 86 (3-4) ◽  
pp. 533-561 ◽  
Author(s):  
Ville Anton Vuorinen ◽  
Harri Hillamo ◽  
Ossi Kaario ◽  
Mika Nuutinen ◽  
Martti Larmi ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document