High-order approximations for an incompressible viscous flow on a rough boundary

In this paper, we propose approximations of fluid flow that could be used for obtaining wall laws of higher order. We consider the two-dimensional laminar fluid flow, modeled by the incompressible Stokes system in a straight channel with a rough side. The roughness is periodic and the ratio of the amplitude of the rough part and the size of the flow domain is denoted by ?, being a small number. We impose periodic boundary conditions on the flow. We generalize the boundary layers needed for the construction of flow approximations of higher order with respect to ?. The existence of the layers and their features are discussed. Finally we give the error estimates for the approximations and establish an explicit wall law.     

Numerical simulation of axisymmetric afterbody flows with jet exhaust

A fifth-order accurate compact difference scheme was used to compute the flow over an axisymmetric afterbody with jet exhaust. The solution was based on the mass-averaged Navier-Stokes equations combined with a two-parameter differential model of turbulence. The computations were performed on a specially generated mesh such that the flow in the exterior and interior of the nozzle could be described simultaneously. Numerical results are presented for various external flow conditions and various chamber pressures.     

Implicit finite element schemes for stationary compressible particle-laden gas flows

The derivation of macroscopic models for particle-laden gas flows is reviewed. Semi-implicit and Newton-like finite element methods are developed for the stationary two-fluid model governing compressible particle-laden gas flows. The Galerkin discretization of the inviscid fluxes is potentially oscillatory and unstable. To suppress numerical oscillations, the spatial discretization is performed by a high-resolution finite element scheme based on algebraic flux correction. A multidimensional limiter of TVD type is employed. An important goal is the efficient computation of stationary solutions in a wide range of Mach numbers. This is a challenging task due to oscillatory correction factors associated with TVD-type flux limiters and the additional strong nonlinearity caused by interfacial coupling terms. A semi-implicit scheme is derived by a time-lagged linearization of the nonlinear residual, and a Newton-like method is obtained in the limit of infinite CFL numbers. The original Jacobian is replaced by a low-order approximation. Special emphasis is laid on the numerical treatment of weakly imposed boundary conditions. It is shown that the proposed approach offers unconditional stability and faster convergence rates for increasing CFL numbers. The strongly coupled solver is compared to operator splitting techniques, which are shown to be less robust.     

求解一维定常对流扩散反应方程的混合型紧致差分格式

借助显式紧致格式和隐式紧致格式的思想,基于截断误差余项修正,并结合原方程本身,构造出了一种求解一维定常对流扩散反应方程的高精度混合型紧致差分格式.格式仅用到三个点上的未知函数值及一阶导数值,而一阶导数值利用四阶Pade格式进行计算,格式整体具有四阶精度.数值实验结果验证了格式的精确性和可靠性.     

四维抛物型方程的高精度分支稳定的显式差分格式

对四维抛物型方程构造了一个高精度显格式,格式的稳定性条件为r=Δt/Δx2=△t/Δy2=△t/△z2=Δt/Δw2<1/2,截断误差阶达到O(Δt2 Δx4),通过数值实验,表明本文理论分析的正确性和文中格式较同类格式的优越性.     

Finite element simulation of compressible particle-laden gas flows

A macroscopic two-fluid model of compressible particle-laden gas flows is considered. The governing equations are discretized by a high-resolution finite element method based on algebraic flux correction. A multidimensional limiter of TVD type is employed to constrain the local characteristic variables for the continuous gas phase and conservative fluxes for a suspension of solid particles. Special emphasis is laid on the efficient computation of steady state solutions at arbitrary Mach numbers. To avoid stability restrictions and convergence problems, the characteristic boundary conditions are imposed weakly and treated in a fully implicit manner. A two-way coupling via the interphase drag force is implemented using operator splitting. The Douglas-Rachford scheme is found to provide a robust treatment of the interphase exchange terms within the framework of a fractional-step solution strategy. Two-dimensional simulation results are presented for a moving shock wave and for a steady nozzle flow.     

Ergodic measures of geodesic flows on compact Lie groups

Let Ψ be the geodesic flow associated with a two-sided invariant metric on a compact Lie group. In this paper, we prove that every ergodic measure μ of Ψ is supported on the unit tangent bundle of a flat torus. As an application, all Lyapunov exponents of μ are zero hence μ is not hyperbolic. Our underlying manifolds have nonnegative curvature (possibly strictly positive on some sections), whereas in contrast, all geodesic flows related to negative curvature are Anosov hence every ergodic measure is hyperbolic.     

Fourth-order alternating direction implicit compact finite difference schemes for two-dimensional Schrödinger equations

In this paper, alternating direction implicit compact finite difference schemes are devised for the numerical solution of two-dimensional Schrödinger equations. The convergence rates of the present schemes are of order O(h⁴+τ²). Numerical experiments show that these schemes preserve the conservation laws of charge and energy and achieve the expected convergence rates. Representative simulations show that the proposed schemes are applicable to problems of engineering interest and competitive when compared to other existing procedures.     

On the structure of internal dissipation of composite compact schemes for gasdynamic simulation

The structure of the internal dissipation terms in composite compact schemes intended for gasdynamic simulation is considered. The main cause of the insufficient stability of high-order accurate schemes is indicated. A method for controlling the dissipative properties of schemes is proposed that makes it possible to compute compressible gas flows with strong shock waves. The supersonic turbulent unsteady flow past a two-dimensional cavity directed toward the stream is computed.     

A large deviations principle for stochastic flows of viscous fluids

We study the well-posedness of a stochastic differential equation on the two dimensional torus ${\mathbb{T}}^2$, driven by an infinite dimensional Wiener process with drift in the Sobolev space $L^2(0,T;H^1\left({\mathbb{T}}^2\right))$. The solution corresponds to a stochastic Lagrangian flow in the sense of DiPerna Lions. By taking into account that the motion of a viscous incompressible fluid on the torus can be described through a suitable stochastic differential equation of the previous type, we study the inviscid limit. By establishing a large deviations principle, we show that, as the viscosity goes to zero, the Lagrangian stochastic Navier–Stokes flow approaches the Euler deterministic Lagrangian flow with an exponential rate function.     

Direct numerical simulation of turbulent flows in eccentric pipes

A numerical algorithm was developed for solving the incompressible Navier-Stokes equations in curvilinear orthogonal coordinates. The algorithm is based on a central-difference discretization in space and on a third-order accurate semi-implicit Runge-Kutta scheme for time integration. The discrete equations inherit some properties of the original differential equations, in particular, the neutrality of the convective terms and the pressure gradient in the kinetic energy production. The method was applied to the direct numerical simulation of turbulent flows between two eccentric cylinders. Numerical computations were performed at Re = 4000 (where the Reynolds number Re was defined in terms of the mean velocity and the hydraulic diameter). It was found that two types of flow develop depending on the geometric parameters. In the flow of one type, turbulent fluctuations were observed over the entire cross section of the pipe, including the narrowest gap, where the local Reynolds number was only about 500. The flow of the other type was divided into turbulent and laminar regions (in the wide and narrow parts of the gap, respectively).     

The Neumann stability of high-order symmetric schemes for convection-diffusion problems

We find a series of sufficient conditions for the Neumann stability of the leapfrog-Euler scheme of arbitrary even order in the spatial coordinates for the three-dimensional convection-diffusion equation.     

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.     

Exponential stability for a nonlinear one-dimensional heat-conductive viscous real gas

In this paper, we prove the existence and the exponential stability in H₊ⁱ (i=1,2), an incomplete metric subspace of Hⁱ×Hⁱ×Hⁱ (i=1,2), of a nonlinear C₀-semigroup S(t) for a nonlinear one-dimensional heat-conductive viscous real gas in bounded domain Ω=(0,1). The results in this paper improve those previously related results.     

Software package for 3D viscous gas flow simulation on multiprocessor computer systems

A parallel software package designed for the numerical simulation of three-dimensional viscous gas flows is presented. The numerical algorithm is based on kinetically consistent difference schemes used on locally refined grids. The software package has been tested in various super-and subsonic flow problems. It provides an opportunity for the direct simulation of turbulent flows. The efficiency of parallelization is analyzed depending on the problem size and the number of processors.     

Global spherically symmetric flows for a viscous radiative and reactive gas in an exterior domain

We consider the equations for a viscous, compressible, radiative and reactive gas (pressure $P=R\rho \theta +a{\theta }^4/3$, internal energy $e=c_v\theta +a{\theta }^4/\rho$) over an unbounded exterior domain in ${\mathbb{R}}^n$, where $n\ge 2$ is the space dimension. The existence, uniqueness, and large-time behavior of global spherically symmetric solutions are established for large initial data. The key point in the analysis is to deduce certain uniform a priori estimates on the solutions, especially on lower and upper bounds of the specific volume and temperature.     

Application of compact and multioperator schemes to the numerical simulation of acoustic fields generated by instability waves in supersonic jets

Acoustic fields generated by instability waves in supersonic jets were numerically simulated. A seventh-order multioperator scheme was used to solve the Euler equations linearized about the mean flow field in an axisymmetric turbulent jet. The mean field was computed using fifth-order compact approximations of the convective terms under conditions similar to experimental data. The numerical results were found to agree well with the experiment.     

Application of compound finite difference schemes for the calculation of time-dependent flows with narrow thermal layers

A procedure for the calculation of nonsteady hydrodynamic flows with narrow thermal layers is described. It is used for the simulation of high speed heating by friction under the shock compression of a plate on a wedge.