Acta mathematica scientia, Series B >
A MATHEMATICAL THEORY FOR LES CONVERGENCE
Received date: 2011-11-02
Online published: 2012-01-20
Supported by
This work is supported in part by the Nuclear Energy University Program of
the Department of Energy, project NEUP-09-349, Battelle Energy Alliance LLC 00088495 (subaward with DOE
as prime sponsor), Leland Stanford Junior University 2175022040367A (subaward with DOE as prime sponsor),
Army Research O?ce W911NF0910306. The work of D.H. Sharp was supported by the US Department of
Energy under Contract DEAC 5206NA25396.
Practical simulations of turbulent processes are generally cutoff, with a grid
resolution that stops within the inertial range, meaning that multiple active regions and
length scales occur below the grid level and are not resolved. This is the regime of large
eddy simulations (LES), in which the larger but not the smaller of the turbulent length
scales are resolved. Solutions of the fluid Navier-Stokes equations, when considered in the
inertial regime, are conventionally regarded as solutions of the Euler equations. In other
words, the viscous and di?usive transport terms in the Navier-Stokes equations can be
neglected in the inertial regime and in LES simulations, while the Euler equation becomes
fundamental.
For such simulations, significant new solution details emerge as the grid is refined. It
follows that conventional notions of grid convergence are at risk of failure, and that a new,
and weaker notion of convergence may be appropriate. It is generally understood that the
LES or inertial regime is inherently fluctuating and its description must be statistical in
nature. Here we develop such a point of view systematically, based on Young measures,
which are measures depending on or indexed by space time points. In the Young measure
d (ξ)ν ξ x,t, the random variable ξ of the measure is a solution state variable, i.e., a solution
dependent variable, representing momentum, density, energy and species concentrations,
while the space time coordinates,x,t , serve to index the measure.
Theoretical evidence suggests that convergence via Young measures is suffciently weak to
encompass the LES/inertial regime; numerical and theoretical evidence suggests that this
notion may be required for passive scalar concentration and thermal degrees of freedom.
Our objective in this research is twofold: turbulent simulations without recourse to ad-
justable parameters (calibration) and extension to more complex physics, without use of
additional models or parameters, in both cases with validation through comparison to
experimental data.
Key words： numerical methods; turbuent mixing
H. Lim , T. Kaman , Y. Yu , V. Mahadeo , Y. Xu , H. Zhang , J. Glimm , S. Dutta , D. H. Sharp , B. Plohr . A MATHEMATICAL THEORY FOR LES CONVERGENCE[J]. Acta mathematica scientia, Series B, 2012 , 32(1) : 237 -258 . DOI: 10.1016/S0252-9602(12)60015-0
[1] Aluie H. Compressible turbulence: The cascade and its locality. Submitted, 2011
[2] Aluie H. Scale decomposition in compressible turbulence. J Fluid Mechanics, Submitted, 2011
[3] Aluie H. Scale locality and the inertial range in compressible turbulence. J Fluid Mechanics, Submitted,
2011
[4] Batchelor G. The Theory of Homogeneous Turbulence. Cambridge: Cambridge University Press, 1955
[5] Chen G -Q, Glimm J. Kolmogorov’s theory of turbulence and inviscid limit of the Navier-Stokes equations
in R . Commun Math Phys, 2010, In Press
[6] Dimonte G, Youngs D L, Dimits A, et al. A comparative study of the turbulent Rayleigh-Taylor instabil-
ity using high-resolution three-dimensional numerical simulations: The alpha-group collaboration. Phys
Fluids, 2004, 36: 1668–1693
[7] Ding X, Chen G -Q, Luo P. Convergence of the Lax-Friedrichs scheme for isentropic gas dynamics I and
II. Acta Mathematica Scientia, 1985, 5: 415–432, 433–472
[8] Donzis D A, Yeung P K. Resolution e?ects and scalilng in numerical simu- lations of turbulence. Phys D,
2010, 239: 1278–1287
[9] Evans L C. Weak Convergence Methods for Nonlinear Partial Di?erential Equations. Providence RI:
American Mathematical Society, 1990
[10] Gamba M, Miller V A, Mungal M G, Hanson R K. Ignition and flame structure in a compact inlet/scramjet
combustor model//AIAA-2011-2366, 2011. 17th AIAAInternational Space Planesand HypersonicSystems
and Technologies Conference
[11] Gamba M, Mungal M G, Hanson R K. Ignition and near-wall burning in transverse hydrogen jets in
supersonic crossflow//AIAA-2011-0319, 2011. 49th AIAA Aerospace Sciences Meeting and Aerospace
Exposition
[12] Gangbo W, Westerdickenberg M. Optimal transport for the system of isentropic Euler equations. Comm
Partial Di?er Equ, 2009, 34: 1041–1073
[13] Glimm J, Sharp D H, Kaman T, Lim H. New directions for Rayleigh-Taylor mixing. Philosophical Transac-
tions of The Royal Society A: Turbulent mixing and beyond. 2011, Submitted for publication; Los Alamos
National Laboratory National Laboratory preprint LA UR 11-00423. Stony Brook University preprint
number SUNYSB-AMS-11-01
[14] Hahn M, Drikakis D, Youngs D L, Williams R J R. Richtmyer-Meskov turbulent mixing arising from an
inclined material interface with realistic surface perturbations and reshocked flow. Phys Fluids, 2011, 23:
0461011-11
[15] Kaman T, Lim H, Yu Y, et al. A numerical method for the simulation of turbulent mixing and its basis
in mathematical theory//Lecture Notes on Numerical Methods for Hyperbolic Equations: Theory and
Applications: Short Course Book. London: CRC/Balkema, 2011: 105–129
[16] Kolmogorov A N. Doklady Akad. Nauk SSSR, 1941, 32:141
[17] Leliss C D, Szekelyhidi L. The euler equations as a di?erential inclusion. Ann Math, 2009, 170: 1471–1436
[18] Leliss C D, Szekelyhidi L. On admissibility criteria for weak solutions of the Euler equations. Arch Rational
Mech Appl, 2010, 195: 225–260
[19] Lim H, Iwerks J, Glimm J, Sharp D H. Nonideal Rayleigh-Taylor mixing. Proc Nation Acad Sci, 2010,
107(29): 12786–12792
[20] Lim H, Iwerks J, Yu Y, Glimm J, Sharp D H. Verification and validation of a method for the simulation
of turbulent mixing. Physica Scripta, 2010, T142: 014014
[21] Lim H, Yu Y, Glimm J, Li X -L, Sharp D H. Chaos, transport, and mesh convergence for fluid mixing.
Acta Mathematicae Applicatae Sinica, 2008, 24: 355–368
[22] Lim H, Yu Y, Glimm J, Li X L, Sharp D H. Subgrid models in turbulent mixing. Astronomical Society of
the Pacific Conference Series, 2008, 406:42
[23] Lim H, Yu Y, Glimm J, Li X L, Sharp D H. Subgrid models for mass and thermal di?usion in turbulent
mixing. Physica Scripta, 2010, T142: 014062
[24] Lim H, Yu Y, Glimm J, Sharp D H. Mathematical, physical and numerical principles essential for models
of turbulent mixing. IMA Volumes in Mathematics and its Applications. Nonlinear Conservation Laws
and Applications, 2010, 153: 405–414
[25] Lim H, Yu Y, Glimm J, Sharp D H. Nearly discontinuous chaotic mixing. High Energy Density Physics,
2010, 6: 223–226
[26] Mal’ek J, NeˇcasJ,RokytaM,Rouˇziˇcka M. Weak and Measure-valued Solutions to Evolutionary PDEs.
London: Chapman & Hall, 1996
[27] Margolin L G, Rider W J, Grinstein F F. Modelilng turbulent flow through implicit les. J Turbulence,
2006, 7:1–27
[28] Margolin L G, Smolarkiewcz P K, Wyszogrodzki A A. Implicit turbulence modeling for high Reynolds
number flows. Trans ASME, 2002, 124: 862–867
[29] Masser T O. Breaking Temperature Equilibrium in Mixed Cell Hydrodynamics[D]. State University of
New York at Stony Brook, 2007
[30] McComb W D. The Physics of Fluid Turbulence. Oxford: Oxford University Press, 1990
[31] Moin P, Squires K, Cabot W, Lee S. A dynamic subgrid-scale model for compressible turbulence and scalar
transport. Phys Fluids, 1991, A3: 2746–2757
[32] Monin A S, Yaglom A M. Statistical Fluid Mechanics: Mechanics of Turbulence. Cambridge, MA: MIT
Press, 1971
[33] Mueschke N, Schilling O. Investigation of Rayleigh-Taylor turbulence and mixing using direct numerical
simulation with experimentally measured initial conditions. i. Comparison to experimental data. Phys
Fluids, 2009, 21: 014106 1–19
[34] Mueschke N, Schilling O. Investigation of Rayleigh-Taylor turbulence and mixing using direct numerical
simulation with experimentally measured initial conditions. ii. Dynamics of transitional flow and mixing
statistics. Phys Fluids, 2009, 21: 014107 1–16
[35] Perna R D. Global existence of solutions to nonlinear hyperbolic systems of conservation laws. J Di?er
Equ, 1967, 20: 187–212
[36] Sche?er V. An inviscid flow with compact support in space-time. J Geom Anal, 1993, 3: 343–401
[37] Shnirelman A. On the nonuniqueness of weak solutions of the Euler equations. Comm Pure Appl Math,
1997, 50: 1261–1286
[38] Smeeton V S, Youngs D L. Experimental investigation of turbulent mixing by Rayleigh-Taylor instability
(part 3). AWE Report Number 0 35/87, 1987
[39] Thorner B, Drikakis D, Youngs D L, Williams R J R. The influence of initial conditions on turbulent
mixing due to Richtmyer-Meshkov insstability. J Fluid Mech, 2010, 654: 99–139
[40] Yakhot A, Orszag S A, Yakhot V, Israeli M. Renormalization group analysis of formulation of large-eddy
simulation. J Sci Comp, 1989, 4: 139–157
[41] Yakhot V, Orszag S A. Renormalization group analysis of turbulence. Phys Rev Lett, 1986, 57: 1722–1724
[42] Yakhot V, Orszag S A. Renormalization group analysis of turbulence I: Basic theory. J Sci Comp, 1986,
1:3–52
/
| 〈 |
|
〉 |