Extension of a simplified Baldwin–Lomax model to nonequilibrium turbulence: Model,analysis and algorithms

Complex turbulence not at statistical equilibrium is impossible to simulate using eddy viscosity models due to a backscatter. This research presents the way to correct the Baldwin–Lomax model for nonequilibrium effects and gives an analysis of the energy evolution in the corrected model. Furthermore, a finite element approximation of the corrected eddy model with first‐order and second‐order time discretization is also presented. A numerical test is given to support the theory.     

Energy dissipation in the Smagorinsky model of turbulence

The Smagorinsky model often severely over-dissipates flows and, consistently, previous estimates of its energy dissipation rate blow up as $Re\to \infty$. This report estimates time averaged model dissipation, $〈{\epsilon }_S〉$, under periodic boundary conditions as$〈{\epsilon }_S〉\le 2\frac{U^3}{L}+{\mathrm{Re}}^{-1}\frac{U^3}{L}+\frac{32}{27}{C}_S^2{\left(\frac{\delta }{L}\right)}^2\frac{U^3}{L}\text{,}$ where $U,L$ are global velocity and length scales and $C_S\simeq 0.1,\delta <1$ are model parameters. Thus, in the absence of boundary layers, the Smagorinsky model does not over dissipate.     

Damping Functions correct over‐dissipation of the Smagorinsky Model

This paper studies the time‐averaged energy dissipation rate ?ε _{S M D} (u )? for the combination of the Smagorinsky model and damping function. The Smagorinsky model is well known to over‐damp. One common correction is to include damping functions that reduce the effects of model viscosity near walls. Mathematical analysis is given here that allows evaluation of ?ε _{S M D} (u )? for any damping function. Moreover, the analysis motivates a modified van Driest damping. It is proven that the combination of the Smagorinsky with this modified damping function does not over‐dissipate and is also consistent with Kolmogorov phenomenology. Copyright © 2017 John Wiley & Sons, Ltd.     

Numerical approximation of the Smagorinsky model by a two-level method

A two-level method for discretizing the Smagorinsky model for the numerical simulation of turbulent flows is proposed. In the two-level algorithm, the solution to the fully nonlinear coarse mesh problem is utilized in a single-step linear fine mesh problem. When modeling parameters are chosen appropriately, the error in the two-level algorithm is comparable to the error in solving the fully nonlinear problem on the fine mesh. We provide an a priori error estimate for the two-level method, which yields appropriate scalings between the coarse and fine mesh-sizes (H and h, respectively), and the radius of the spatial filter used in the Smagorinsky model (δ). In addition, we provide an algorithm in which a coarse mesh correction is performed to further enhance the accuracy. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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     

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.     

反映强流动曲率效应的非线性湍流模型

首先定性地分析了流线曲率效应对流场湍流结构的影响,然后以U型槽道流为典型算例,对多种湍流模型进行了评估．评估的模型包括:线性涡粘性模型,二阶和三阶非线性涡粘性模型,二阶显式代数应力模型和Reynolds应力模型．评估结果表明,性能良好的三阶非线性涡粘性模型,如黄于宁等人发展的HM模型以及CLS模型,可以较好地描述流线的曲率效应对湍流结构的影响,如凸曲率作用下内壁附近湍流强度的衰减和凹曲率作用下外壁附近湍流的增强,并且较好地确定了管道下游的分离点位置和分离泡长度,其预测的结果和实验符合较好,与Reynolds力模型的结果十分接近,因此可以较好地应用于具有曲率效应的工程湍流的计算．     

Numerical analysis of two ensemble eddy viscosity numerical regularizations of fluid motion

This report analyzes an efficient ensemble regularization algorithm for under‐resolved and convection dominated flows (including ones at higher Reynolds numbers). Computing an ensemble simultaneously allows each realization to access ensemble data. This allows use of means and fluctuations in regularizations used for each realization. The combined approach of ensemble time stepping and ensemble regularizations allows direct calculation of the turbulent viscosity coefficient and gives an unconditionally stable algorithm. It also suggests reconsidering an old but not as well‐developed definition of the mixing length. This mixing length vanishes at solid walls without van Driest damping, increases stability, and improves flow predictions in our preliminary tests. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 630–651, 2015     

Numerical analysis of a transient eddy current axisymmetric problem involving velocity terms

The aim of this article is to analyze a transient axisymmetric electromagnetic model involving velocity terms in the Ohm's law. To this end, we introduce a time‐dependent weak formulation leading to a degenerate parabolic problem and establish its well posedness.We propose a finite‐element method for space discretization and prove well posedness and error estimates. Then, we combine it with a backward Euler time discretization and prove stability and error estimates. Finally, numerical results assessing the performance of the method are reported. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

湍流不同尺度间的近程和共振作用规律及应用

从Navier-Stokes方程出发,研究了湍流不同尺度间的相互作用规律,给出相近尺度间近程粘性应力的积分和微分表达式．引入极相近尺度之间共振相互作用的概念,得到共振粘性应力的微分表达式．利用共振粘性应力张量获得不含经验关系和常数、近似封闭的大涡模拟（LES）方程组．利用近程和共振粘性应力张量获得不含经验关系和常数、近似封闭的湍流多尺度方程组．讨论了湍流多尺度方程的性质及用于湍流计算的优点,尺度间相互作用的近程特性说明:多尺度模拟是湍流计算很有价值的方法,并列举了算例．     

Numerical analysis of a method for high Peclet number transport in porous media

A variationally consistent eddy viscosity discretization is presented in [W.J. Layton, A connection between subgrid scale eddy viscosity and mixed methods, Appl. Math. Comput. 133 (2002) 147-157] for the stationary convection diffusion problem. This discretization is extended to the evolutionary problem in [N. Heitmann, Subgridscale stabilization of time-dependent convection dominated diffusive transport, J. Math. Anal. Appl. 331 (2007) 38-50] with a near optimal error bound. In the following, we couple this discretization with the porous media problem. We present a comprehensive analysis of stability and error for the velocity field derived from the porous media problem. Next, using a backward Euler approximation for the time derivative we follow the inherited error in velocity through the coupling with the convection diffusion problem. The method is shown to be stable and the error near optimal and independent of the diffusion coefficient, ?.     

A dynamic subgrid-scale model for the large eddy simulation of stratified flow

A new dynamic subgrid-scale (SGS) model, including subgrid turbulent stress and heat flux models for stratified shear flow is proposed by using Yoshizawa’s eddy viscosity model as a base model. Based on our calculated results, the dynamic subgrid-scale model developed here is effective for the large eddy simulation (LES) of stratified turbulent channel flows. The new SGS model is then applied to the large eddy simulation of stratified turbulent channel flow under gravity to investigate the coupled shear and buoyancy effects on the near-wall turbulent statistics and the turbulent heat transfer at different Richardson numbers. The critical Richardson number predicted by the present calculation is in good agreement with the value of theoretical analysis     

Numerical analysis of a dual-phase-lag model involving two temperatures

In this paper, we numerically analyse a phase-lag model with two temperatures which arises in the heat conduction theory. The model is written as a linear partial differential equation of third order in time. The variational formulation, written in terms of the thermal acceleration, leads to a linear variational equation, for which we recall an existence and uniqueness result and an energy decay property. Then, using the finite element method to approximate the spatial variable and the implicit Euler scheme to discretize the time derivatives, fully discrete approximations are introduced. A discrete stability property is proved, and a priori error estimates are obtained, from which the linear convergence of the approximation is derived. Finally, some one-dimensional numerical simulations are described to demonstrate the accuracy of the approximation and the behaviour of the solution.     

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.     

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).     

Numerical analysis of a model of cure process for composites

We consider a composite material composed of fibres included in a resin which becomes solid when it is heated up (reaction of reticulation). The mathematical modelling of the cure process is given by a kinetic equation describing the evolution of the reaction of reticulation coupled with the heat equation. In this paper, we are interested in the computation of approximate solutions. We propose a family of discretized problems depending on two parameters (β₁, β₂) ε [0, 1]² which split the linear and non‐ linear terms in implicit and explicit parts. We prove the stability and convergence of the discretization for any (β₁, β₂) ε [½, 1 ] × [0, 1]. We present also some numerical results. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical analysis of a mean field model of superconducting vortices

A finite volume/element approximation of a mean field modelof superconducting vortices in one and two dimensions is presented.The solutions of these approximations are investigated. A finiteelement approximation of the free boundary problem which isa special case of the steady state solution of the model isalso studied. We present some computed results from these approximations. Received 3 December 1997. Accepted 17 May 2000.     

Global dynamical analysis of a heroin epidemic model on complex networks

In this paper, a heroin epidemic model on complex networks is proposed. By the next generation matrix, the basic reproduction number $R_0$ is obtained. If $R_0<1$, then the drug-free equilibrium is globally asymptotically stable. If $R_0>1$, there is an unique endemic equilibrium and it is also globally asymptotically stable. Our results show that if the degree of the network is large enough, the drug transmission always spreads. Sensitivity analysis of the basic reproduction number with the various parameters in the model are carried out to verify the important effects for control the drug transmission. Some simulations illustrate our theoretical results     

Numerical analysis of a model for phase separation of a multicomponent alloy

We consider a fully discrete implicit finite-element approximationof a model for the phase separation of a multi-component alloy.We prove existence, uniqueness and stability of the numericalsolution for a sufficiently small time step. We prove convergenceto the solution of the associated continuous problem. We performa linear stability analysis of the equation and describe somenumerical experiments.     

Numerical analysis of a mathematical model of the myocardial blood-supply system

A hydraulic model of the hemodynamics of the arterial part of the myocardium is considered, and a numerical analysis of the model is conducted. Computer experiments are used to investigate the dependence of blood flows on parametric and structural changes in the system. Translated from Prikladnaya Matematika i Informatika, No. 2, pp. 56–62, 1999.