A similarity theory of approximate deconvolution models of turbulence

We apply the phenomenology of homogeneous, isotropic turbulence to the family of approximate deconvolution models proposed by Stolz and Adams. In particular, we establish that the models themselves have an energy cascade with two asymptotically different inertial ranges. Delineation of these gives insight into the resolution requirements of using approximate deconvolution models. The approximate deconvolution model's energy balance contains both an enhanced energy dissipation and a modification to the model's kinetic energy. The modification of the model's kinetic energy induces a secondary energy cascade which accelerates scale truncation. The enhanced energy dissipation completes the scale truncation by reducing the model's micro-scale from the Kolmogorov micro-scale.     

Existence of smooth attractors for the Navier-Stokes-omega model of turbulence

We consider the Navier-Stokes- model, given by

     

Finite element error analysis of a zeroth order approximate deconvolution model based on a mixed formulation

A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is sought. This is a low order model, but its importance lies in the insight that it provides for the analysis of the higher order models actually used in practice by the pioneers Stolz and Adams [N.A. Adams, S. Stolz, On the approximate deconvolution procedure for LES, Phys. Fluids 2 (1999) 1699-1701; N.A. Adams, S. Stolz, Deconvolution methods for subgrid-scale approximation in large eddy simulation, in: B.J. Geurts (Ed.), Modern Simul. Strategies for Turbulent Flow, Edwards, Philadelphia, 2001, pp. 21-44] and others. The higher order models have proven to be of high accuracy. However, stable discretizations of them have proven to be tricky and other stabilizations, such as time relaxation and eddy viscosity, are often added. We propose a discretization based on a mixed variational formulation that gives the correct energy balance. We show it to be unconditionally stable and prove convergence.     

Large eddy simulation for turbulent magnetohydrodynamic flows

We investigate the mathematical properties of a model for the simulation of large eddies in turbulent, electrically conducting, viscous, incompressible flows. We prove existence and uniqueness of solutions for the simplest (zeroth) closed MHD model (1.7), we show that its solutions converge to the solution of the MHD equations as the averaging radii converge to zero, and derive a bound on the modeling error. Furthermore, we show that the model preserves the properties of the 3D MHD equations: the kinetic energy and the magnetic helicity are conserved, while the cross helicity is approximately conserved and converges to the cross helicity of the MHD equations, and the model is proven to preserve the Alfvén waves, with the velocity converging to that of the MHD, as δ₁,δ₂ tend to zero. We perform computational tests that verify the accuracy of the method and compare the conserved quantities of the model to those of the averaged MHD.     

A mathematical and physical Study of multiscale deconvolution models of turbulence

This work presents a rigorous analysis of mathematical and physical properties for solutions of multiscale deconvolution turbulence models. We show that solutions of these models exactly conserve model quantities for the integral invariants of fundamental physical importance: kinetic energy, helicity, and (in two dimensions) enstrophy. The kinetic energy conservation is the key that allows us to next apply the phenomenology of homogeneous, isotropic turbulence to establish the existence of a model energy cascade and, in particular, that the cascade exhibits enhanced energy dissipation in a secondary accelerated cascade, which ends at the model's microscale (which we establish is larger than the Kolmogorov microscale). We also prove that the model dissipates energy at the same rate as true turbulent flow, ～ O(U³ ∕ L), independent of Reynolds number. Lastly, we prove the existence of global attractors for the model solutions; the proof of which also shows that solutions are actually one degree of regularity higher than previously known. Copyright © 2012 John Wiley & Sons, Ltd.     

On Taylor/Eddy solutions of approximate deconvolution models of turbulence

This article shows that so called general Green–Taylor solutions, also called Taylor solutions or eddy solutions, of the Navier–Stokes equations are also exact solutions to approximate deconvolution models of turbulence. Thus, these special structures in flows exist as exact features in the models studied and their persistence/transient behavior is exactly determined by their stability, not by the effects of modelling or truncation errors.     

On the well-posedness and conservation laws of a family of multiscale deconvolution models for the magnetohydrodynamics equations

We present a mathematical study of a large eddy simulation (LES) model for the incompressible magnetohydrodynamics equations. The classical closure problem arising for LES models is solved with the multiscale deconvolution technique developed by Dunca in [11]. We prove the model admits unique, regular weak solutions and provide a mathematical study of the modeling error.     

On the convergence of an approximate deconvolution model to the 3D mean Boussinesq equations

In this paper, we study an LES model for the approximation of large scales of the 3D Boussinesq equations. This model is obtained using the approach first described by Stolz and Adams, based on the Van Cittern approximate deconvolution operators, and applied to the filtered Boussinesq equations. Existence and uniqueness of a regular weak solution are provided. Our main objective is to prove that this solution converges towards a solution of the filtered Boussinesq equations, as the deconvolution parameter goes to zero. Copyright © 2014 John Wiley & Sons, Ltd.     

Conservation laws of turbulence models

The conservation of mass, momentum, energy, helicity, and enstrophy in fluid flow are important because these quantities organize a flow, and characterize change in the flow's structure over time. In turbulent flow, conservation laws remain important in the inertial range of wave numbers, where viscous effects are negligible. It is in the inertial range where energy, helicity (3d), and enstrophy (2d) must be accurately cascaded for a turbulence model to be qualitatively correct. A first and necessary step for an accurate cascade is conservation; however, many turbulent flow simulations are based on turbulence models whose conservation properties are little explored and might be very different from those of the Navier-Stokes equations.We explore conservation laws and approximate conservation laws satisfied by LES turbulence models. For the Leray, Leray deconvolution, Bardina, and Nth order deconvolution models, we give exact or approximate laws for a model mass, momentum, energy, enstrophy and helicity. The possibility of cascades for model quantities is also discussed.     

Numerical simulation of UV disinfection reactors: Evaluation of alternative turbulence models

Six turbulence models, including standard k–ε, k–ε RNG, k–ω (88), revised k–ω (98), Reynolds stress transport model (RSTM), and two-fluid model (TFM), were applied to the simulation of a closed conduit polychromatic UV reactor. Predicted flow field and turbulent kinetic energy were compared with the experimental data from a digital particle image velocimetry (DPIV). All of the predicted flow fields were combined with a multiple segment source summation (MSSS) fluence rate model and three different microbial response kinetic models to simulate the disinfection process at two UV lamp power conditions. Microbial transport was simulated using the Lagrangian particle tracking method. The results show that the fluence distributions and the effluent inactivation levels were sensitive to the turbulence model selection. The level of sensitivity was a function of the operating conditions and the UV response kinetics of the microorganisms. Simulations with operating conditions that produced higher log inactivation or utilized microorganisms with higher UV sensitivity showed greater sensitivity to the turbulence model selection. In addition, a broader fluence distribution was found with turbulence models that predicted a larger wake region behind the lamps.     

Energy analysis and improved regularity estimates for multiscale deconvolution models of incompressible flows

This paper presents new analytical results and the first numerical results for a recently proposed multiscale deconvolution model (MDM) recently proposed. The model involves a large‐eddy simulation closure that uses a novel deconvolution approach based on the introduction of two distinct filtering length scales. We establish connections between the MDM and two other models, and, on the basis of one of these connections, we establish an improved regularity estimate for MDM solutions. We also prove that the MDM preserves Taylor‐eddy solutions of the Navier–Stokes equations and therefore does not distort this particular vortex structure. Simulations of the MDM are performed to examine the accuracy of the MDM and the effect of the filtering length scales on energy spectra for three‐dimensional homogeneous and isotropic flows. Numerical evidence for all tests clearly indicates that the MDM gives very accurate coarse‐mesh solutions and that this multiscale approach to deconvolution is effective. Copyright © 2015 John Wiley & Sons, Ltd.     

A simple and stable scale-similarity model for large Eddy simulation: Energy balance and existence of weak solutions

In averaging the Navier-Stokes equations, the problem of closure arises. Scale-similarity models address closure by (roughly speaking) extrapolation from the (known) resolved scales to the (unknown) unresolved scales. In a posteriori tests, scale-similarity models are often the most accurate but can prove to be unstable when used in a numerical simulation. In this report, we consider the scale-similarity model given by

. We prove it is stable (solutions satisfy an energy inequality) and deduce from that the existence of weak solutions of the model.     

Performance of turbulence models for flows through rough pipes

The Reynolds-averaged Navier–Stokes (RANS) equations were solved along with turbulence models, namely k–ε, k–ω, Reynolds stress models (RSM), and filtered Navier–Stokes equations along with Large Eddy Simulation (LES) to study the fully-developed turbulent flows in circular pipes roughened by repeated square ribs with various spacings. Solutions of these flows were obtained using the commercial computational fluid dynamics (CFD) software Fluent. The numerical results were validated against experimental measurements and other numerical data published in literature. The performance of the turbulence models was compared and discussed. All the RANS models and LES model were observed to perform equally well in predicting the time-averaged flow statistics. However no instantaneous information can be obtained from the RANS results. Therefore, when a rough overview of the flow process in a pipe roughened by repeated ribs is needed, any one of the RANS models can be of value. On the other hand, the instantaneous as well as time-averaged flows could be studied with more insight using LES, albeit at a cost of CPU effort at least one order higher.     

Numerical analysis and computational testing of a high accuracy Leray‐deconvolution model of turbulence

We study a computationally attractive algorithm (based on an extrapolated Crank‐Nicolson method) for a recently proposed family of high accuracy turbulence models, the Leray‐deconvolution family. First we prove convergence of the algorithm to the solution of the Navier‐Stokes equations and delineate its (optimal) accuracy. Numerical experiments are presented which confirm the convergence theory. Our 3d experiments also give a careful comparison of various related approaches. They show the combination of the Leray‐deconvolution regularization with the extrapolated Crank‐Nicolson method can be more accurate at higher Reynolds number that the classical extrapolated trapezoidal method of Baker (Report, Harvard University, 1976). We also show the higher order Leray‐deconvolution models (e.g. N = 1,2,3) have greater accuracy than the N = 0 case of the Leray‐α model. Numerical experiments for the 2d step problem are also successfully investigated. Around the critical Reynolds number, the low order models inhibit vortex shedding behind the step. The higher order models, correctly, do not. To estimate the complexity of using Leray‐deconvolution models for turbulent flow simulations we estimate the models' microscale.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Finite element analysis of laminar and turbulent flows using LES and subgrid-scale models

Numerical simulations of laminar and turbulent flows in a lid driven cavity and over a backward-facing step are presented in this work. The main objectives of this research are to know more about the structure of turbulent flows, to identify their three-dimensional characteristic and to study physical effects due to heat transfer. The filtered Navier–Stokes equations are used to simulate large scales, however they are supplemented by subgrid-scale (SGS) models to simulate the energy transfer from large scales toward subgrid-scales, where this energy will be dissipated by molecular viscosity. Two SGS models are applied: the classical Smagorinsky’s model and the Dynamic model for large eddy simulation (LES). Both models are implemented in a three-dimensional finite element code using linear tetrahedral elements. Qualitative and quantitative aspects of two and three-dimensional flows in a lid-driven cavity and over a backward-facing step, using LES, are analyzed comparing numerical and experimental results obtained by other authors.     

Effects of plane strain on inhomogeneous turbulence

The influence of the mean plane strain on the turbulence transportation is investigated by large eddy simulation (LES) in the shearless turbulence mixing layer. It is found that the mean strains enhance the turbulent fluctuations in the mixing region. Compression in the inhomogeneous direction can greatly increase the transport of turbulent kinetic energy by triple correlation terms, while stretching in the inhomogeneous direction decreases the turbulence transportation. The gradient diffusion models for turbulent transportation are evaluated and it is found that the intermittency consideration can improve the prediction ability of the gradient-type models for the triple correlation terms. Project supported by the Sino-French Laboratory in Beijing and the National Natural Science Foundation of China (Grant No. 19572041).