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.     

Simulation of the turbulent flow in a model hydraulic spool valve

The paper presents results of the simulation of the flow in a model hydraulic spool valve by means of Large Eddy Simulation (LES). The flow is characterized by a very low Reynolds number which is untypical for LES. The emphasis of the study is the analysis of the flow and the investigation of the influence of the subgrid-scale model. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Stochastic parameterization for large eddy simulation of geophysical flows

Recently, stochastic, as opposed to deterministic, parameterizations are being investigated to model the effects of unresolved subgrid scales (SGS) in large eddy simulations (LES) of geophysical flows. We analyse such a stochastic approach in the barotropic vorticity equation to show that (i) if the stochastic parameterization approximates the actual SGS stresses, then the solution of the stochastic LES approximates the ``true" solution at appropriate scale sizes; and that (ii) when the filter scale size approaches zero, the solution of the stochastic LES approaches the true solution.

     

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.     

Convergence of approximate deconvolution models to the mean Navier–Stokes equations

We consider a 3D Approximate Deconvolution Model ADM which belongs to the class of Large Eddy Simulation (LES) models. We aim at proving that the solution of the ADM converges towards a dissipative solution of the mean Navier–Stokes equations. The study holds for periodic boundary conditions. The convolution filter we first consider is the Helmholtz filter. We next consider generalized convolution filters for which the convergence property still holds.     

Large eddy simulation of an ethylene-air turbulent premixed V-flame

Large eddy simulation (LES) using a dynamic eddy viscosity subgrid scale stress model and a fast-chemistry combustion model without accounting for the finite-rate chemical kinetics is applied to study the ignition and propagation of a turbulent premixed V-flame. A progress variable c-equation is applied to describe the flame front propagation. The equations are solved two dimensionally by a projection-based fractional step method for low Mach number flows. The flow field with a stabilizing rod without reaction is first obtained as the initial field and ignition happens just upstream of the stabilizing rod. The shape of the flame is affected by the velocity field, and following the flame propagation, the vortices fade and move to locations along the flame front. The LES computed time-averaged velocity agrees well with data obtained from experiments.     

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.     

Numerical simulation of cavitation around a two-dimensional hydrofoil using VOF method and LES turbulence model

In this paper simulation of cavitating flow over the Clark-Y hydrofoil is reported using the large eddy simulation (LES) turbulence model and volume of fluid (VOF) technique. We applied an incompressible LES modelling approach based on an implicit method for the subgrid terms. To apply the cavitation model, the flow has been considered as a single fluid, two-phase mixture. A transport equation model for the local volume fraction of vapour is solved and a finite rate mass transfer model is used for the vapourization and condensation processes. A compressive volume of fluid (VOF) method is applied to track the interface of liquid and vapour phases. This simulation is performed using a finite volume, two phase solver available in the framework of the OpenFOAM (Open Field Operation and Manipulation) software package. Simulation is performed for the cloud and super-cavitation regimes, i.e., σ = 0.8, 0.4, 0.28. We compared the results of two different mass transfer models, namely Kunz and Sauer models. The results of our simulation are compared for cavitation dynamics, starting point of cavitation, cavity’s diameter and force coefficients with the experimental data, where available. For both of steady state and transient conditions, suitable accuracy has been observed for cavitation dynamics and force coefficients.     

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.     

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.     

On the structural stability of thermoelastic model of porous media

In the present paper we study the structural stability of the mathematical model of the linear thermoelastic materials with voids. We prove that the solutions of problems depend continuously on the constitutive quantities, which may be subjected to error or perturbations in the mathematical modelling process. Thus, we assume to have changes in the various coupling coefficients of the model and then we establish estimates of continuous dependence of solutions. We have to outline that such estimates play a central role in obtaining approximations to these kinds of problems. To derive a priori estimates for a solution we first establish appropriate bounds for the solutions of certain auxiliary problems. These are achieved by means of so‐called Rellich‐like identities. We also investigate how the solution in the coupled model behaves as some coupling coefficients tend to zero. Copyright © 2007 John Wiley & Sons, Ltd.     

On the benefits of ODT-based stochastic turbulence modeling

We summarize the group's progress in applying, analyzing, and improving ODT and ODT-based stochastic turbulence models like ODTLES. Compared to DNS these models span a wider range of scales while compared to RANS/LES (i) the molecular effects are retained and (ii) no assumption of scale separation is made. In this regard ODTLES has more properties of DNS than of standard LES. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulations de l'écoulement turbulent marin avec un modèle de déconvolution

We display a continous equation for the deconvolution process that generalizes the Van Cittert algorithm in the case of oceanic boundary conditions for a given fixed wind. We deduce a LES model for which we have existence and uniqueness of a strong solution. Finally, we display several numerical simulations showing the practical interest of the model. To cite this article: A.-C. Bennis et al., C. R. Acad. Sci. Paris, Ser. I 347 (2009).     

Global existence and stability of solution of a reaction-diffusion model for cancer invasion

In this paper, we study a mathematical model of cancer invasion proposed by Gatenby and Gawlinski. The model is a strongly coupled degenerate reaction-diffusion system. Very few mathematical results are known for this system. We investigate the global existence of classical solutions for the system by using energy estimates and the bootstrap arguments, and global asymptotic stability of equilibrium points of the system by Lyapunov functions.     

Evaluating the effect of gap size in a single function mathematical programming model for the three-group classification problem

This study examines the impact that the size of the classification gap can have on the classificatory performance of a mathematical programming based discriminant model. In mathematical programming based models that project the discriminant scores onto a line, the discriminant score of an observation may fall into the gap between adjacent group intervals; thus there is no clear cut way to determine the group in which the observation should be classified. We examine a procedure that we refer to as the split gap approach. The split gap approach is defined as a strategy of estimating the performance of a mathematical programming based model using a nonzero gap size to separate group intervals and then splitting the gap between adjacent group intervals to classify future observations. Studies that propose models with a classification gap generally do not assess the effect of the gap on the performance of the model. This paper investigates this effect. A theoretical assessment and a Monte Carlo simulation are used to determine the impact of different gap sizes on a mixed integer programming model using a single function classification model for the three-group case.     

Mathematical model of the Galathea-belt toroidal magnetic trap

We consider a mathematical model of equilibrium configurations of plasma, magnetic field, and electric field in a toroidal trap with two ring conductors with current loaded into plasma. We present the mathematical apparatus of the model based on the numerical solution of boundary value problems for the Grad–Shafranov equation (a differential equation of elliptic type for the magnetic flux function), solution methods for these problems, and numerically obtained properties of equilibrium configurations. We indicate the differences in configurations in the toroidal trap and in its analog straightened into a cylinder.     

Existence of solutions to the biaxial nematic liquid crystals with two order parameter tensors

The recent paper considers some mathematical problems on a tensor model with two order parameters to the biaxial nematic liquid crystals. The energy in this model contains a cubic term and is unbounded from below. We study the gradient flow generated by this energy in two dimensions. We focus on establishing global existence of classical and weak solutions.     

Estimates for the probability of ruin with special emphasis on the possibility of large claims

The present paper investigates, for the general Andersen model, the asymptotic behaviour of the probability of ruin function when the initial risk reserve tends to infinity. Whereas the exponential (Cramér) case is well understood, in the past, less attention has been paid to a systematic study of a model taking big claim sizes into account. We give a thorough treatment of the latter and also review previously known but mostly scattered results to show how they all follow from essentially one mathematical model.     

On the global asymptotic stability of the equilibrium of the Lotka-Volterra equations in a varying environment

We consider a mathematical model describing the evolution of prey and predator populations in a varying environment. The model is a system of ordinary differential equations with a unique nontrivial equilibrium. We derive sufficient conditions for the global asymptotic stability of the equilibrium.