TECHNICAL PAPERS: Injection and Fuel Air Mixture Preparation

Effect of Subgrid Modeling on the In-Cylinder Unsteady Mixing Process in a Direct Injection Engine

[+] Author and Article Information
K. Sone, S. Menon

School of Aerospace Engineering, Georgia Institute of Technology, Atlanta, GA 30332

J. Eng. Gas Turbines Power 125(2), 435-443 (Apr 29, 2003) (9 pages) doi:10.1115/1.1501918 History: Received May 23, 2001; Revised December 01, 2001; Online April 29, 2003
Copyright © 2003 by ASME
Your Session has timed out. Please sign back in to continue.


Anderson, R. W., Yang, J., Brehob, D. D., Vallance, J. K., and Whiteaker, R. M., 1996, “Understanding the Thermodynamics of Direct Injection Spark Ignition (DISI) Combustion Systems: An Analytical and Experimental Investigation,” SAE Paper No. 962018.
Iwamoto, Y., Noma, K., Nakayama, O., Yamauchi, T., and Ando, H., 1997, “Development of Gasoline Direct Injection Engine,” SAE Paper No. 970541.
Galperin, B., and Orszag, S., eds., 1993, Large Eddy Simulation of Complex Engineering and Geophysical Flows, Cambridge University Press, New York.
Menon,  S., 2000, “Subgrid Combustion Modelling for Large-Eddy Simulations,” Int. J. Engine Res., 1(2), pp. 209–227.
Smagorinsky,  J., 1963, “General Circulation Experiments With the Primitive Equations: I. The Basic Equations,” Mon. Weather Rev., 91, pp. 99–164.
Mungal,  M. G., and Dimotakis,  P. E., 1984, “Mixing and Combustion With Low Heat Release in a Turbulent Shear Layer,” J. Fluid Mech., 148, pp. 349–382.
Menon, S., McMurtry, P. A., and Kerstein, A. R., 1993, “A Linear Eddy Mixing Model for Large Eddy Simulation of Turbulent Combustion,” Large Eddy Simulation of Complex Engineering and Geophysical Flows, B. Galperin and S. Orszag, eds., Cambridge University Press, New York, pp. 287–314.
Han,  Z., and Reitz,  R. D., 1995, “Turbulence Modeling of Internal Combustion Engines Using RNG k-ε Models,” Combust. Sci. Technol., 106, pp. 267–295.
Celik,  I., Yavuz,  I., Smirnov,  A., Smith,  J., Amin,  E., and Gel,  A., 2000, “Prediction of In-Cylinder Turbulence for IC Engines,” Combust. Sci. Technol., 153, pp. 339–368.
Menon,  S., Yeung,  P.-K., and Kim,  W.-W., 1996, “Effect of Subgrid Models on the Computed Interscale Energy Transfer in Isotropic Turbulence,” Comput. Fluids, 25(2), pp. 165–180.
Kim,  W.-W., Menon,  S., and Mongia,  H. C., 1999, “Numerical Simulations of Reacting Flows in a Gas Turbine Combustor,” Combust. Sci. Technol., 143, pp. 25–62.
Kim,  W.-W., and Menon,  S., 1999, “A New Incompressible Solver for Large-Eddy Simulations,” Int. J. Numer. Methods Fluids, 31, pp. 983–1017.
Chakravarthy,  V. K., and Menon,  S., 2000, “Modeling of Turbulent Premixed Flames in the Flamelet Regime,” Combust. Sci. Technol., 162, pp. 1–50.
Chakravarthy,  V. K., and Menon,  S., 2000, “Subgrid Modeling of Premixed Flames in the Flamelet Regime,” Flow, Turbul. Combust., 65, pp. 131–161.
Amsden, A. A., O’Rourke, P. J., and Butler, T. D., 1989, “KIVA-II: A Computer Program for Chemically Reactive Flows With Sprays,” Technical Report LA-11560-MS, Los Alamos National Laboratory, May.
Pannala, S., and Menon, S., 2000, “On Large Eddy Simulations of Reacting Two-Phase Flows,” Technical Report CCL-00-006, Georgia Institute of Technology, Atlanta, GA, May.
Menon,  S., and Calhoon,  W., 1996, “Subgrid Mixing and Molecular Transport Modeling for Large-Eddy Simulations of Turbulent Reacting Flows,” Proc. Combust. Inst., 26 , pp. 59–66.
Kerstein,  A. R., 1991, “Linear-Eddy Modelling of Turbulent Transport. Part 6. Microstructure of Diffusive Scalar Mixing Fields,” J. Fluid Mech., 231, pp. 361–394.
Chakravarthy,  V. K., and Menon,  S., 2000, “Linear-Eddy Simulations of Schmidt and Reynolds Number Scaling of Scalar Mixing,” Phys. Fluids, 13, pp. 488–499.
Sone, K., Patel, N. V., and Menon, S., 2001, “KIVALES—Large-Eddy Simulations of Internal Combustion Engines. Part I: Theory and Formulation,” Technical Report CCL-2001-008, Georgia Institute of Technology, Atlanta, GA.
Metcalfe,  R. W., Orszag,  S. A., Brachet,  M. E., Menon,  S., and Riley,  J. J., 1987, “Secondary Instability of a Temporally Growing Mixing Layer,” J. Fluid Mech., 184, pp. 207–243.
Riley, J. J., and Metcalfe, R. W., 1980, “Direct Numerical Simulation of a Perturbed, Turbulent Mixing Layer,” AIAA Paper No. 80-0274.
Michalke,  A., 1964, “On the Inviscid Instability of the Hyperbolic-Tangent Velocity Profile,” J. Fluid Mech., 19, pp. 543–556.
Morse,  A. P., Whitelaw,  J. H., and Yianneskis,  M., 1979, “Turbulent Flow Measurements by Laser-Doppler Anemometry in Motored Piston-Cylinder Assemblies,” J. Fluids Eng., 101, pp. 208–216.
Haworth,  D. C., and Jansen,  K., 2000, “Large-Eddy Simulation on Unstructured Deforming Meshes: Towards Reciprocating IC Engines,” Comput. Fluids, 29, pp. 493–524.
Weclas,  M., Melling,  A., and Durst,  F., 1998, “Flow Separation in the Inlet Valve Gap of Piston Engines,” Prog. Energy Combust. Sci., 24, pp. 165–195.
Nelson, C., and Menon, S., 1998, “Unsteady Simulations of Compressible Spatial Mixing Layers,” AIAA Paper No. 98-0786.
Han, Z., Fan, L., and Reitz, R. D., 1997, “Multidimensional Modeling of Spray Atomization and Air-Fuel Mixing in a Direct-Injection Spark-Ignition Engine,” SAE Paper No. 970884.


Grahic Jump Location
Azimuthal vorticity at crank angle=210 (during injection). Contour intervals are fixed with Δω=1000(sec−1). Solid contours show positive vorticity while dotted contours show negative vorticity. Even with fuel injection and vaporization, KIVALES still shows turbulent flow field. (a) KIVA-3V, start of injection=180, baseline grid. (b) KIVALES, start of injection=180, baseline grid. (c) KIVALES, start of injection=180, fine grid.
Grahic Jump Location
Fuel mass fraction (gaseous) at crank angle=210. Contour intervals are ΔY=0.03. A bold line denotes the stoichiometric surface. Figures on the left show contours in the symmetric plane and the ones on the right show in the plane perpendicular to the symmetric plane (the location is shown as a line in the left figures). (a) KIVA-3V, (b) KIVALES, (c) KIVALES-LEM.
Grahic Jump Location
Temporal evolution of scalar (fuel mass fraction) mean (μ) and the variance (σ2) standard deviation (σ) is shown for comparison between gradient diffusion closure KIVALES and KIVALES-LEM. Light lines denote mean while black lines denote standard deviation. Note that the linear eddy model (LEM) predicts larger scale variance between bottom dead center and top dead center.
Grahic Jump Location
PDF of fuel mass fraction normalized by the standard deviation at crank angle=360 where mean and variance are nearly same for all three cases while the shape of PDF is different in each case.
Grahic Jump Location
Turbulent kinetic energy evolution with crank angle. Light lines indicate KIVA-3V while dark lines indicate KIVALES. KIVA-3V models are turbulent scales except the largest scale motions, therefore, the turbulent kinetic energy in KIVA-3V is larger than the subgrid scale turbulent kinetic energy in KIVALES.
Grahic Jump Location
Turbulent kinetic energy at crank angle=210. Contour intervals are fixed with Δk=1.0×105 (cm2/sec2) for KIVALES and Δk=2.0×105 (cm2/sec2) for KIVA-3V, which is twice the magnitude obtained with KIVALES. (a) KIVA-3V, start of injection=180. (b) KIVALES, start of injection, start of injection=180.
Grahic Jump Location
Resolved scale fundamental mode energy growth. Growth rate of KIVALES is 0.19, which is also predicted by an earlier study by 23 whereas KIVA-3V rapidly dissipates.
Grahic Jump Location
Comparison of KIVA-3V and KIVALES predictions of turbulent mean and rms profiles in the experimental configuration. Symbols represent experimental data (24), solid lines and dot-dashed lines represent KIVALES and KIVA-3V, respectively. Scale of 1.0 on the x-axis corresponds to 400 cm/s in (a) and 1200 cm/s in (b), respectively. (a) Mean velocity profile at 36 deg crank angle. (b) rms velocity profile at 36 deg crank angle.
Grahic Jump Location
Geometry used for the internal combustion engine simulations with two valves, one is the intake and the other is the exhaust. The symmetric boundary condition is employed on the plane of symmetry. The pictures shown here are the baseline geometry (51,000 cells).




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In