Navier-Stokes-alpha model: LES equations with nonlinear dispersion

We present a framework for discussing LES equations with nonlinear dispersion. In this framework, we discuss the properties of the nonlinearly dispersive Navier-Stokes-alpha model of incompressible fluid turbulence --- also called the viscous Camassa-Holm equations and the LANS equations in the literature --- in comparison with the corresponding properties of large eddy simulation (LES) equations obtained via the approximate-inverse approach. In this comparison, we identify the spatially filtered NS-alpha equations with a class of generalized LES similarity models. Applying a certain approximate inverse to this class of LES models restores the Kelvin circulation theorem for the defiltered velocity and shows that the NS-alpha model describes the dynamics of the defiltered velocity for this class of generalized LES similarity models. We also show that the subgrid scale forces in the NS-alpha model transform covariantly under Galilean transformations and under a change to a uniformly rotating reference frame. Finally, we discuss in the spectral formulation how the NS-alpha model retains the local interactions among the large scales, retains the nonlocal sweeping effects of large scales on small scales, yet attenuates the local interactions of the small scales amongst themselves.Comment: 15 pages, no figures, Special LES volume of ERCOFTAC bulletin, to appear in 200

The Navier-Stokes-alpha model of fluid turbulence

We review the properties of the nonlinearly dispersive Navier-Stokes-alpha (NS-alpha) model of incompressible fluid turbulence -- also called the viscous Camassa-Holm equations and the LANS equations in the literature. We first re-derive the NS-alpha model by filtering the velocity of the fluid loop in Kelvin's circulation theorem for the Navier-Stokes equations. Then we show that this filtering causes the wavenumber spectrum of the translational kinetic energy for the NS-alpha model to roll off as k^{-3} for k \alpha > 1 in three dimensions, instead of continuing along the slower Kolmogorov scaling law, k^{-5/3}, that it follows for k \alpha < 1. This rolloff at higher wavenumbers shortens the inertial range for the NS-alpha model and thereby makes it more computable. We also explain how the NS-alpha model is related to large eddy simulation (LES) turbulence modeling and to the stress tensor for second-grade fluids. We close by surveying recent results in the literature for the NS-alpha model and its inviscid limit (the Euler-alpha model).Comment: 22 pages, 1 figure. Dedicated to V. E. Zakharov on the occasion of his 60th birthday. To appear in Physica

Consequences of Symmetries on the Analysis and Construction of Turbulence Models

Since they represent fundamental physical properties in turbulence (conservation laws, wall laws, Kolmogorov energy spectrum, ...), symmetries are used to analyse common turbulence models. A class of symmetry preserving turbulence models is proposed. This class is refined such that the models respect the second law of thermodynamics. Finally, an example of model belonging to the class is numerically tested.Comment: Published in SIGMA (Symmetry, Integrability and Geometry: Methods and Applications) at http://www.emis.de/journals/SIGMA

Three regularization models of the Navier-Stokes equations

We determine how the differences in the treatment of the subfilter-scale physics affect the properties of the flow for three closely related regularizations of Navier-Stokes. The consequences on the applicability of the regularizations as SGS models are also shown by examining their effects on superfilter-scale properties. Numerical solutions of the Clark-alpha model are compared to two previously employed regularizations, LANS-alpha and Leray-alpha (at Re ~ 3300, Taylor Re ~ 790) and to a DNS. We derive the Karman-Howarth equation for both the Clark-alpha and Leray-alpha models. We confirm one of two possible scalings resulting from this equation for Clark as well as its associated k^(-1) energy spectrum. At sub-filter scales, Clark-alpha possesses similar total dissipation and characteristic time to reach a statistical turbulent steady-state as Navier-Stokes, but exhibits greater intermittency. As a SGS model, Clark reproduces the energy spectrum and intermittency properties of the DNS. For the Leray model, increasing the filter width decreases the nonlinearity and the effective Re is substantially decreased. Even for the smallest value of alpha studied, Leray-alpha was inadequate as a SGS model. The LANS energy spectrum k^1, consistent with its so-called "rigid bodies," precludes a reproduction of the large-scale energy spectrum of the DNS at high Re while achieving a large reduction in resolution. However, that this same feature reduces its intermittency compared to Clark-alpha (which shares a similar Karman-Howarth equation). Clark is found to be the best approximation for reproducing the total dissipation rate and the energy spectrum at scales larger than alpha, whereas high-order intermittency properties for larger values of alpha are best reproduced by LANS-alpha.Comment: 21 pages, 8 figure

On the need for a nonlinear subscale turbulence term in POD models as exemplified for a high Reynolds number flow over an Ahmed body

We investigate a hierarchy of eddy-viscosity terms in POD Galerkin models to account for a large fraction of unresolved fluctuation energy. These Galerkin methods are applied to Large Eddy Simulation data for a flow around the vehicle-like bluff body call Ahmed body. This flow has three challenges for any reduced-order model: a high Reynolds number, coherent structures with broadband frequency dynamics, and meta-stable asymmetric base flow states. The Galerkin models are found to be most accurate with modal eddy viscosities as proposed by Rempfer & Fasel (1994). Robustness of the model solution with respect to initial conditions, eddy viscosity values and model order is only achieved for state-dependent eddy viscosities as proposed by Noack, Morzynski & Tadmor (2011). Only the POD system with state-dependent modal eddy viscosities can address all challenges of the flow characteristics. All parameters are analytically derived from the Navier-Stokes based balance equations with the available data. We arrive at simple general guidelines for robust and accurate POD models which can be expected to hold for a large class of turbulent flows.Comment: Submitted to the Journal of Fluid Mechanic

Stochastic representation of the Reynolds transport theorem: revisiting large-scale modeling

We explore the potential of a formulation of the Navier-Stokes equations incorporating a random description of the small-scale velocity component. This model, established from a version of the Reynolds transport theorem adapted to a stochastic representation of the flow, gives rise to a large-scale description of the flow dynamics in which emerges an anisotropic subgrid tensor, reminiscent to the Reynolds stress tensor, together with a drift correction due to an inhomogeneous turbulence. The corresponding subgrid model, which depends on the small scales velocity variance, generalizes the Boussinesq eddy viscosity assumption. However, it is not anymore obtained from an analogy with molecular dissipation but ensues rigorously from the random modeling of the flow. This principle allows us to propose several subgrid models defined directly on the resolved flow component. We assess and compare numerically those models on a standard Green-Taylor vortex flow at Reynolds 1600. The numerical simulations, carried out with an accurate divergence-free scheme, outperform classical large-eddies formulations and provides a simple demonstration of the pertinence of the proposed large-scale modeling

Inertial Range Scaling, Karman-Howarth Theorem and Intermittency for Forced and Decaying Lagrangian Averaged MHD in 2D

We present an extension of the Karman-Howarth theorem to the Lagrangian averaged magnetohydrodynamic (LAMHD-alpha) equations. The scaling laws resulting as a corollary of this theorem are studied in numerical simulations, as well as the scaling of the longitudinal structure function exponents indicative of intermittency. Numerical simulations for a magnetic Prandtl number equal to unity are presented both for freely decaying and for forced two dimensional MHD turbulence, solving directly the MHD equations, and employing the LAMHD-alpha equations at 1/2 and 1/4 resolution. Linear scaling of the third-order structure function with length is observed. The LAMHD-alpha equations also capture the anomalous scaling of the longitudinal structure function exponents up to order 8.Comment: 34 pages, 7 figures author institution addresses added magnetic Prandtl number stated clearl