One major drawback of the eddy viscosity subgrid‐scale stress models used in large‐eddy simulations is their inability to represent correctly with a single universal constant different turbulent fields in rotating or sheared flows, near solid walls, or in transitional regimes. In the present work a new eddy viscosity model is presented which alleviates many of these drawbacks. The model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model is based on an algebraic identity between the subgrid‐scale stresses at two different filtered levels and the resolved turbulent stresses. The subgrid‐scale stresses obtained using the proposed model vanish in laminar flow and at a solid boundary, and have the correct asymptotic behavior in the near‐wall region of a turbulent boundary layer. The results of large‐eddy simulations of transitional and turbulent channel flow that use the proposed model are in good agreement with the direct simulation data.

A dynamic subgrid‐scale eddy viscosity model

We present an overview of the lattice Boltzmann method (LBM), a parallel and efficient algorithm for simulating single-phase and multiphase fluid flows and for incorporating additional physical complexities. The LBM is especially useful for modeling complicated boundary conditions and multiphase interfaces. Recent extensions of this method are described, including simulations of fluid turbulence, suspension flows, and reaction diffusion systems.

Lattice boltzmann method for fluid flows

A comprehensive review of spatiotemporal pattern formation in systems driven away from equilibrium is presented, with emphasis on comparisons between theory and quantitative experiments. Examples include patterns in hydrodynamic systems such as thermal convection in pure fluids and binary mixtures, Taylor-Couette flow, parametric-wave instabilities, as well as patterns in solidification fronts, nonlinear optics, oscillatory chemical reactions and excitable biological media. The theoretical starting point is usually a set of deterministic equations of motion, typically in the form of nonlinear partial differential equations. These are sometimes supplemented by stochastic terms representing thermal or instrumental noise, but for macroscopic systems and carefully designed experiments the stochastic forces are often negligible. An aim of theory is to describe solutions of the deterministic equations that are likely to be reached starting from typical initial conditions and to persist at long times. A unified description is developed, based on the linear instabilities of a homogeneous state, which leads naturally to a classification of patterns in terms of the characteristic wave vector q0 and frequency ω0 of the instability. Type Is systems (ω0=0, q0≠0) are stationary in time and periodic in space; type IIIo systems (ω0≠0, q0=0) are periodic in time and uniform in space; and type Io systems (ω0≠0, q0≠0) are periodic in both space and time. Near a continuous (or supercritical) instability, the dynamics may be accurately described via "amplitude equations," whose form is universal for each type of instability. The specifics of each system enter only through the nonuniversal coefficients. Far from the instability threshold a different universal description known as the "phase equation" may be derived, but it is restricted to slow distortions of an ideal pattern. For many systems appropriate starting equations are either not known or too complicated to analyze conveniently. It is thus useful to introduce phenomenological order-parameter models, which lead to the correct amplitude equations near threshold, and which may be solved analytically or numerically in the nonlinear regime away from the instability. The above theoretical methods are useful in analyzing "real pattern effects" such as the influence of external boundaries, or the formation and dynamics of defects in ideal structures. An important element in nonequilibrium systems is the appearance of deterministic chaos. A greal deal is known about systems with a small number of degrees of freedom displaying "temporal chaos," where the structure of the phase space can be analyzed in detail. For spatially extended systems with many degrees of freedom, on the other hand, one is dealing with spatiotemporal chaos and appropriate methods of analysis need to be developed. In addition to the general features of nonequilibrium pattern formation discussed above, detailed reviews of theoretical and experimental work on many specific systems are presented. These include Rayleigh-Benard convection in a pure fluid, convection in binary-fluid mixtures, electrohydrodynamic convection in nematic liquid crystals, Taylor-Couette flow between rotating cylinders, parametric surface waves, patterns in certain open flow systems, oscillatory chemical reactions, static and dynamic patterns in biological media, crystallization fronts, and patterns in nonlinear optics. A concluding section summarizes what has and has not been accomplished, and attempts to assess the prospects for the future.

/pdf/pattern-formation-outside-of-equilibrium-1601b9znum.pdf

Pattern formation outside of equilibrium

We review the field of cavity optomechanics, which explores the interaction between electromagnetic radiation and nano- or micromechanical motion This review covers the basics of optical cavities and mechanical resonators, their mutual optomechanical interaction mediated by the radiation pressure force, the large variety of experimental systems which exhibit this interaction, optical measurements of mechanical motion, dynamical backaction amplification and cooling, nonlinear dynamics, multimode optomechanics, and proposals for future cavity quantum optomechanics experiments In addition, we describe the perspectives for fundamental quantum physics and for possible applications of optomechanical devices

Cavity Optomechanics

In this article the principles of the field operation and manipulation (FOAM) C++ class library for continuum mechanics are outlined. Our intention is to make it as easy as possible to develop reliable and efficient computational continuum-mechanics codes: this is achieved by making the top-level syntax of the code as close as possible to conventional mathematical notation for tensors and partial differential equations. Object-orientation techniques enable the creation of data types that closely mimic those of continuum mechanics, and the operator overloading possible in C++ allows normal mathematical symbols to be used for the basic operations. As an example, the implementation of various types of turbulence modeling in a FOAM computational-fluid-dynamics code is discussed, and calculations performed on a standard test case, that of flow around a square prism, are presented. To demonstrate the flexibility of the FOAM library, codes for solving structures and magnetohydrodynamics are also presented with appropriate test case results given. © 1998 American Institute of Physics.

A tensorial approach to computational continuum mechanics using object-oriented techniques

We develop the dynamic renormalization group (RNG) method for hydrodynamic turbulence. This procedure, which uses dynamic scaling and invariance together with iterated perturbation methods, allows us to evaluate transport coefficients and transport equations for the large-scale (slow) modes. The RNG theory, which does not include any experimentally adjustable parameters, gives the following numerical values for important constants of turbulent flows: Kolmogorov constant for the inertial-range spectrumCK=1.617; turbulent Prandtl number for high-Reynolds-number heat transferPt=0.7179; Batchelor constantBa=1.161; and skewness factor¯S3=0.4878. A differentialK-\(\bar \varepsilon \) model is derived, which, in the high-Reynolds-number regions of the flow, gives the algebraic relationv=0.0837 K2/\(\bar \varepsilon \), decay of isotropic turbulence asK=O(t−1.3307), and the von Karman constantκ=0.372. A differential transport model, based on differential relations betweenK,\(\bar \varepsilon \), andν, is derived that is not divergent whenK→ 0 and\(\bar \varepsilon \) is finite. This latter model is particularly useful near walls.

Renormalization group analysis of turbulence I. Basic theory

Turbulence models are developed by supplementing the renormalization group (RNG) approach of Yakhot and Orszag [J. Sci. Comput. 1, 3 (1986)] with scale expansions for the Reynolds stress and production of dissipation terms. The additional expansion parameter (η≡SK/■) is the ratio of the turbulent to mean strain time scale. While low‐order expansions appear to provide an adequate description for the Reynolds stress, no finite truncation of the expansion for the production of dissipation term in powers of η suffices−terms of all orders must be retained. Based on these ideas, a new two‐equation model and Reynolds stress transport model are developed for turbulent shear flows. The models are tested for homogeneous shear flow and flow over a backward facing step. Comparisons between the model predictions and experimental data are excellent.

/pdf/development-of-turbulence-models-for-shear-flows-by-a-double-37jwuzybet.pdf

Development of turbulence models for shear flows by a double expansion technique

Using renormalization-group methods and the postulated equivalence between the inertial-range structures of turbulent flows satisfying initial and boundary conditions and of flows driven by a random force, we evaluate the Kolmogorov constant (1.617) and Batchelor constant (1.161), skewness factor (0.4878), power-law exponent (1.3307) for the decay of homogeneous turbulence, turbulent Prandtl number (0.7179), and von K\'arm\'an constant (0.372). This renormalization-group technique has also been used to derive turbulent transport models.

Renormalization-group analysis of turbulence.

Complex fluid physics can be modeled using an extended kinetic (Boltzmann) equation in a more efficient way than using the continuum Navier-Stokes equations. Here, we explain this method for modeling fluid turbulence and show its effectiveness with the use of a computationally efficient implementation in terms of a discrete or "lattice" Boltzmann equation.

http://science.sciencemag.org/content/sci/301/5633/633.full.pdf

Extended Boltzmann kinetic equation for turbulent flows.

We reformulate the renormalization group (RNG) and theɛ-expansion for derivation of turbulence models. The procedure is developed for the Navier-Stokes equations and the transport equations for the kinetic energyK and energy dissipation rate ℰ. The derivation draws on the works of Yakhot and Orszag (1986) and Smith and Reynolds (1992), and all results are found at low order in the underlying perturbation expansion in powers ofɛ. The sum of the source terms in the ℰ-equation is known to beO(1) due to the balance at leading order ofO(R

 T
 1/2
 ) terms. Smith and Reynolds (1992) showed the cancellation of some of theO(R

 T
 1/2
 ) terms generated by the RNG procedure. Here we show that including the random-force contribution to ℰ-production results in the cancellation ofall theO(R

 1/2
 ) terms. We find that two of theO(1) terms in the RNG equation for the mean dissipation rate ℰ have the same form as those in the widely used model
 -equation. The values of the coefficients of the familiar terms are close to those used in practice. An extra production term is predicted which is small for slow distortions, but important for rapid distortions. Hence, it may be a term that should be added to the
 model equation. We believe that the present derivation places the
 model equation on a more solid theoretical basis.

Victor Yakhot

Papers

Renormalization group analysis of turbulence I. Basic theory

Development of turbulence models for shear flows by a double expansion technique

Renormalization-group analysis of turbulence.

Extended Boltzmann kinetic equation for turbulent flows.

The renormalization group, the e-expansion and derivation of turbulence models