Acta mathematica scientia,Series B ›› 2012, Vol. 32 ›› Issue (1): 237-258.doi: 10.1016/S0252-9602(12)60015-0

• Articles • Previous Articles     Next Articles

A MATHEMATICAL THEORY FOR LES CONVERGENCE

H. Lim T. Kaman1, Y. Yu, V. Mahadeo1, Y. Xu, H. Zhang,1 J. Glimm1, S. Dutta2 D. H. Sharp3, B. Plohr3   

  1. 1.Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY 11794, USA|2.College of St. Catherine, Madison, NJ, USA|3.Los Alamos National Laboratory, Los Alamos, NM 87544, USA
  • Received:2011-11-02 Online:2012-01-20 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.

Abstract:

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

CLC Number: 

  • 76F25
Trendmd