An improved dynamic subgrid-scale model and its application to large eddy simulation of stratified channel flows

In the present paper, a new dynamic subgrid-scale (SGS) model of turbulent stress and heat flux for stratified shear flow is proposed. Based on our calculated results of stratified channel flow, the dynamic subgrid-scale model developed in this paper is shown to be effective for large eddy simulation (LES) of stratified turbulent shear flows. The new SGS model is then applied to the LES of the stratified turbulent channel flow to investigate the coupled shear and buoyancy effects on the behavior of turbulent statistics, turbulent heat transfer and flow structures at different Richardson numbers.     

A formulation on large eddy simulation

This note presents a formulation on large eddy simulation with the idea motivated by Fourier series representation. A set of filtered equations are obtained.     

Discrete large eddy simulation

IntroductionIn previous worksll,2], we discussed several issues associated with the standard version oflarge eddy simulation (LES) such as filtering and averaging. By the standard version, we meanthe traditional practice of first constructing a set of field equations of motions for turbulentmotion and then discretizing the equations to suit computational simulations.In this paper, we revisit the issue of large eddy simulation with a view towards improvingthe filtering procedure that is used i…     

Regularity and uniqueness for the stationary large eddy simulation model

In the note we are concerned with higher regularity and uniqueness of solutions to the stationary problem arising from the large eddy simulation of turbulent ows. The system of equations contains a nonlocal nonlinear term, which prevents straightforward application of a difference quotients method. The existence of weak solutions was shown in A. Świerczewska: Large eddy simulation. Existence of stationary solutions to the dynamical model, ZAMM, Z. Angew. Math. Mech. 85 (2005), 593–604 and P. Gwiazda, A. Świerczewska: Large eddy simulation turbulence model with Young measures, Appl. Math. Lett. 18 (2005), 923–929.     

A dynamical approach to large eddy simulation of turbulent flows: existence of weak solutions

We consider a system of equations coming from turbulence models using a large eddy simulation (LES) technique. The idea of this approach bases on decomposing the velocity into a part containing large flow structures and a part consisting of small scales. The equations for large‐scale quantities are derived from the Navier–Stokes equations with an additional constitutive relation for the contribution of small eddies. The mathematical difficulties in this paper focus on the non‐linear and non‐local turbulent term. Copyright © 2005 John Wiley & Sons, Ltd.     

Energy dissipation bounds for hear flows for a model in large eddy simulation

This report gives an upper bound for the time average of the energy dissipation rate, including energy dissipation due to viscous dissipation and turbulent diffusion in a model for the motion of large eddies in a turbulent flow. The bound is only a function of the reference velocity U, the domain diameter L, the eddy size δ, and surprisingly, the Reynolds number.     

A large-eddy-based lattice Boltzmann model for turbulent flow simulation

In this paper a novel and simple large-eddy-based lattice Boltzmann model is proposed to simulate two-dimensional turbulence. Unlike existing lattice Boltzmann models for turbulent flow simulation, which were based on primitive-variables Navier–Stokes equations, the target macroscopic equations of the present model are vorticity-streamfunction equations. Thanks to the intrinsic features of vorticity-streamfunction equations, the present model is efficient, stable and simple for two-dimensional turbulence simulation. The advantages of the present model are validated by numerical experiments.     

Convergence of finite element approximations of large eddy motion

Fluid motion in many applications occurs at higher Reynolds numbers. In these applications dealing with turbulent flow is thus inescapable. One promising approach to the simulation of the motion of the large structures in turbulent flow is large eddy simulation in which equations describing the motion of local spatial averages of the fluid velocity are solved numerically. This report considers “numerical errors” in LES. Specifically, for one family of space filtered flow models, we show convergence of the finite element approximation of the model and give an estimate of the error. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 689–710, 2002; Published online in Wiley InterScience (www.interscience.wiley.com); DOI 10.1002/num.10027.     

Asymptotic behaviour of commutation errors and the divergence of the Reynolds stress tensor near the wall in the turbulent channel flow

The derivation of the space averaged Navier–Stokes equations for the large eddy simulation (LES) of turbulent incompressible flows introduces two groups of terms which do not depend only on the space averaged flow field variables: the divergence of the Reynolds stress tensor and commutation errors. Whereas the former is studied intensively in the literature, the latter terms are usually neglected. This note studies the asymptotic behaviour of these terms for the turbulent channel flow at a wall in the case that the commutation errors arise from the application of a non‐uniform box filter. To perform analytical calculations, the unknown flow field is modelled by a wall law (Reichardt law and 1/αth power law) for the mean velocity profile and highly oscillating functions model the turbulent fluctuations. The asymptotics show that near the wall, the commutation errors are at least as important as the divergence of the Reynolds stress tensor. Copyright © 2006 John Wiley & Sons, Ltd.     

The effect of modeling of velocity fluctuations on prediction of collection efficiency of cyclone separators

The effect of modeling of velocity fluctuations on the prediction of collection efficiency of cyclone separators has been numerically investigated using the Reynolds stress turbulence model (RSTM) and large eddy simulation (LES). The Eulerian–Lagrangian modeling approach of CFD code Fluent 6.3.26 has been employed to simulate the three dimensional, unsteady turbulent gas–solid flows in a Stairmand high efficiency cyclone. The simulated results have been compared with experimental observations available in the literature. The analysis of results shows that the RSTM and the LES have adequately predicted the mean flow field. Results of the present study demonstrate that the LES has good performance on prediction of fluctuating flow field and collection efficiency for each and every particle size. However, the performance of the RSTM is found poor in terms of prediction of velocity fluctuations and collection efficiency, especially for small particles. This relates to the precessing of the vortex core phenomenon, which is resolved more accurately by LES as compared to the RSTM simulation. The results suggest that the prediction of collection efficiency, especially for small particles is greatly influenced by the simulation of velocity fluctuations in cyclones.     

On the existence of an inertial manifold for a deconvolution model of the 2D mean Boussinesq equations

We show the existence of an inertial manifold (ie, a globally invariant, exponentially attracting, finite‐dimensional manifold) for the approximate deconvolution model of the 2D mean Boussinesq equations. This model is obtained by means of the Van Cittern approximate deconvolution operators, which is applied to the 2D filtered Boussinesq equations.     

Large eddy simulation of a sheet/cloud cavitation on a NACA0015 hydrofoil

A single fluid model of sheet/cloud cavitation is developed and applied to a NACA0015 hydrofoil. First, a cavity formation model is set up, based on a three-dimensional (3D) non-cavitation model of Navier–Stokes equations with a large eddy simulation (LES) scheme for weakly compressible flows. A fifth-order polynomial curve is adopted to describe the relationship between density coefficient ratio and pressure coefficient when cavitation occurs. The Navier–Stokes equations including cavitation bubble clusters are solved using the finite-volume approach with time-marching scheme, and MacCormack’s explicit-corrector scheme is adopted. Simulations are carried out in a 3D field acting on a hydrofoil NACA0015 at angles of attack 4°, 8° and 20°, with cavitation numbers σ = 1.0, 1.5 and 2.0, Re = 10⁶, and a 360 × 63 × 29 meshing system. We study time-dependent sheet/cloud cavitation structures, caused by the interaction of viscous objects, such as vortices, and cavitation bubbles. At small angles of attack (4°), the sheet cavity is relatively stable just by oscillating in size at the accumulation stage; at 8° it has a tendency to break away from the upper foil section, with the cloud cavitation structure becoming apparent; at 20°, the flow separates fully from the leading edge of the hydrofoil, and the vortex cavitation occurs. Comparisons with other studies, carried out mainly in the context of flow patterns on which prior experiments and simulations were done, demonstrate the power of our model. Overall, it can snapshot the collapse of cloud cavitation, and allow a study of flow patterns and their instabilities, such as “crescent-shaped regions.”     

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.     

On large Eddy simulation and variational multiscale methods in the numerical simulation of turbulent incompressible flows

Numerical simulation of turbulent flows is one of the great challenges in Computational Fluid Dynamics (CFD). In general, Direct Numerical Simulation (DNS) is not feasible due to limited computer resources (performance and memory), and the use of a turbulence model becomes necessary. The paper will discuss several aspects of two approaches of turbulent modeling—Large Eddy Simulation (LES) and Variational Multiscale (VMS) models. Topics which will be addressed are the detailed derivation of these models, the analysis of commutation errors in LES models as well as other results from mathematical analysis.     

A stratified model of flow in a coastal trench

A stratified model of the circulation in a deep, narrow trench adjacent to a coast is described. The flow in the trench is driven by surface wind stress, coastal runoff and inflow at one end. The model is being developed to investigate the Norwegian coastal current flowing through the Norwegian Trench.     

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     

A dynamic model of power metal-oxide-semiconductor field-effect transistor half-bridges for the fast simulation of switching induced electromagnetic emissions

Hard switching of semiconductors is the main source of conducted electromagnetic emissions (EME) in pulse-width modulation (PWM) driven power inverters. The requirements on the electromagnetic compatibility grow with the increasing number of installed electric motor drives and inductive power converters. An accurate prediction of the conducted EME requires a model which considers the switching transition of the power semiconductors and the parasitic elements. This typically leads to complex SPICE models, which are hardly suitable for fast dynamic simulations and model-based controller design. This paper presents a compact mathematical model of a low voltage half-bridge inverter, which is based on large-signal models for the individual components and allows for the fast simulation of the conducted EME and switching losses. The high accuracy of the proposed mathematical model is demonstrated by measurement results. In particular, it is shown that the model is able to accurately predict the conducted electromagnetic emissions up to 100 MHz.     

A numerical study on statistical temporal scales in inertia particle dispersion

The statistical temporal scales involved in inertia particle dispersion are analyzed numerically. The numerical method of large eddy simulation, solving a filtered Navier-Stokes equation, is utilized to calculate fully developed turbulent channel flows with Reynolds numbers of 180 and 640, and the particle Lagrangian trajectory method is employed to track inertia particles released into the flow fields. The Lagrangian and Eulerian temporal scales are obtained statistically for fluid tracer particles and three different inertia particles with Stokes numbers of 1, 10 and 100. The Eulerian temporal scales, decreasing with the velocity of advection from the wall to the channel central plane, are smaller than the Lagrangian ones. The Lagrangian temporal scales of inertia particles increase with the particle Stokes number. The Lagrangian temporal scales of the fluid phase ‘seen’ by inertia particles are separate from those of the fluid phase, where inertia particles travel in turbulent vortices, due to the particle inertia and particle trajectory crossing effects. The effects of the Reynolds number on the integral temporal scales are also discussed. The results are worthy of use in examining and developing engineering prediction models of particle dispersion.     

Stochastic modeling of unresolved scales in complex systems

Model uncertainties or simulation uncertainties occur in mathematical modeling of multiscale complex systems, since some mechanisms or scales are not represented (i.e., ‘unresolved’) due to a lack in our understanding of these mechanisms or limitations in computational power. The impact of these unresolved scales on the resolved scales needs to be parameterized or taken into account. A stochastic scheme is devised to take the effects of unresolved scales into account, in the context of solving nonlinear partial differential equations. An example is presented to demonstrate this strategy. Dedicated to Professor Peter E. Kloeden on the occasion of his 60th birthday     

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.