Error analysis of a weighted least-squares finite element method for 2-D incompressible flows in velocity–stress–pressure formulation

In this paper we are concerned with a weighted least-squares finite element method for approximating the solution of boundary value problems for 2-D viscous incompressible flows. We consider the generalized Stokes equations with velocity boundary conditions. Introducing the auxiliary variables (stresses) of the velocity gradients and combining the divergence free condition with some compatibility conditions, we can recast the original second-order problem as a Petrovski-type first-order elliptic system (called velocity–stress–pressure formulation) in six equations and six unknowns together with Riemann–Hilbert-type boundary conditions. A weighted least-squares finite element method is proposed for solving this extended first-order problem. The finite element approximations are defined to be the minimizers of a weighted least-squares functional over the finite element subspaces of the H¹ product space. With many advantageous features, the analysis also shows that, under suitable assumptions, the method achieves optimal order of convergence both in the L²-norm and in the H¹-norm. © 1998 B. G. Teubner Stuttgart—John Wiley & Sons, Ltd.     

Efficient Solution of the Compressible Euler Equations for Low Mach Numbers

This paper investigates the benefits of local preconditioning for the compressible Euler equations to predict nearly incompressible fluid flow. The AUSMDV(P) upwind method by Edwards and Liou is employed to maintain the spatial accuracy of the method for low Mach numbers. The results indicate excellent solution quality and fast convergence to steady state for compressible as well as nearly incompressible fluid flow. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non‐homogeneous Navier–Stokes systems with order‐parameter‐dependent stresses

We consider the Navier–Stokes system with variable density and variable viscosity coupled to a transport equation for an order‐parameter c. Moreover, an extra stress depending on c and ?c, which describes surface tension like effects, is included in the Navier–Stokes system. Such a system arises, e.g. for certain models of granular flows and as a diffuse interface model for a two‐phase flow of viscous incompressible fluids. The so‐called density‐dependent Navier–Stokes system is also a special case of our system. We prove short‐time existence of strong solution in L^q‐Sobolev spaces with q>d. We consider the case of a bounded domain and an asymptotically flat layer with a combination of a Dirichlet boundary condition and a free surface boundary condition. The result is based on a maximal regularity result for the linearized system. Copyright © 2010 John Wiley & Sons, Ltd.     

On singular limit for the compressible Navier–Stokes system with non-monotone pressure law

We consider a singular limit for the compressible Navier–Stokes system with general non-monotone pressure law in the asymptotic regime of low Mach number and large Reynolds numbers. We show that any dissipative weak solution approaches the solution of incompressible Euler equation both for well-prepared initial data and ill-prepared initial data.     

On compactness of the velocity field in the incompressible limit of the full Navier–Stokes–Fourier system on large domains

The incompressible limit for the full Navier–Stokes–Fourier system is studied on a family of domains containing balls of the radius growing with a speed that dominates the inverse of the Mach number. It is shown that the velocity field converges strongly to its limit locally in space, in particular, the effect of the sound waves is eliminated by means of the local decay estimates for the acoustic wave equation. Copyright © 2008 John Wiley & Sons, Ltd.     

Life span of 2‐D irrotational compressible fluids in the halfplane

We consider the Euler equations of barotropic inviscid compressible fluids in the half plane. It is well known that, as the Mach number goes to zero, the compressible flows approximate the solution of the equations of motion of inviscid, incompressible fluids. In dimension two such limit solution exists on any arbitrary time interval, with no restriction on the size of the initial velocity. It is then natural to expect the same for the compressible solution, if the Mach number is sufficiently small. We consider smooth irrotational solutions. First, we study the life span, i.e. the largest time interval T(ε) of existence of classical solutions, when the initial data are a small perturbation of size εfrom a constant state. For the proof of this result we use a combination of energy and decay estimates. Then, the estimate of the life span allows to show, by a suitable scaling of variables, the existence of irrotational solutions on any arbitrary time interval, for any small enough Mach number. Copyright © 2002 John Wiley & Sons, Ltd.     

On the spectrum for the artificial compressible system

Stability of stationary solutions of the incompressible Navier–Stokes system and the corresponding artificial compressible system is considered. Both systems have the same sets of stationary solutions and the incompressible system is obtained from the artificial compressible one in the zero limit of the artificial Mach number ? which is a singular limit. It is proved that if a stationary solution of the incompressible system is asymptotically stable and the velocity field of the stationary solution satisfies an energy-type stability criterion by variational method with admissible functions being only potential flow parts of velocity fields, then it is also stable as a solution of the artificial compressible one for sufficiently small ?. The result is applied to the Taylor problem.     

Numerical Modelling of Newtonian and Non-Newtonian Fluids Flow

The work deals with numerical modelling of 2D/3D laminar incompressible viscous flows for Newtonian and non–Newtonian fluids. The unsteady system of Navier–Stokes equations with steady boundary conditions in the form of an artificial compressibility method is solved by multistage Runge–Kutta finite volume method. Steady state solution is achieved for t→∞. Convergence is followed by steady residual behaviour. For unsteady solution high compressibility coefficient β² is considered. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Incompressible Limit and the Initial Layer of the Compressible Euler Equation in ℝn+

We consider the incompressible limit of the compressible Euler equation in the half-space ℝⁿ₊. It is proved that the solutions of the non-dimensionalized compressible Euler equation converge to the solution of the incompressible Euler equation when the Mach number tends to zero. If the initial data v^∞₀ do not satisfy the condition ‘∇⋅v^∞₀=0’, then the initial layer will appear. © 1997 B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Positivity-Preserving Runge-Kutta Discontinuous Galerkin Method on Adaptive Cartesian Grid for Strong Moving Shock

In order to suppress the failure of preserving positivity of density or pressure, a positivity-preserving limiter technique coupled with $h$-adaptive Runge-Kutta discontinuous Galerkin (RKDG) method is developed in this paper. Such a method is implemented to simulate flows with the large Mach number, strong shock/obstacle interactions and shock diffractions. The Cartesian grid with ghost cell immersed boundary method for arbitrarily complex geometries is also presented. This approach directly uses the cell solution polynomial of DG finite element space as the interpolation formula. The method is validated by the well documented test examples involving unsteady compressible flows through complex bodies over a large Mach numbers. The numerical results demonstrate the robustness and the versatility of the proposed approach.     

Large‐eddy simulation of compressible rectangular duct flow

Turbulent flow inside ducts differs significantly from the statistically one‐dimensional turbulent channel flow due to the additional confining lateral walls. Secondary flows in the corners are found which are known to affect the turbulence structure close to the walls. Most computations have been performed for incompressible or only slightly compressible flows and with square duct geometries (e.g. [1,2]). In the present study we address the effect of the aspect ratio on the flow structure in high subsonic rectangular duct flow using LES. Mean‐flow characteristics are analyzed and compared to square duct and channel flow data. We find good agreement with available incompressible DNS square duct data for low Mach number computations. For aspect ratios as high as five, the mean‐flow profile in the center of the longer side of the duct along the wall normal direction differs distinctively from channel flow for the present Mach 0.7 computations.     

Coupling of compressible and incompressible flow regions using the multiple pressure variables approach

In many cases, multiphase flows are simulated on the basis of the incompressible Navier–Stokes equations. This assumption is valid as long as the density changes in the gas phase can be neglected. Yet, for certain technical applications such as fuel injection, this is no longer the case, and at least the gaseous phase has to be treated as a compressible fluid. In this paper, we consider the coupling of a compressible flow region to an incompressible one based on a splitting of the pressure into a thermodynamic and a hydrodynamic part. The compressible Euler equations are then connected to the Mach number zero limit equations in the other region. These limit equations can be solved analytically in one space dimension that allows to couple them to the solution of a half‐Riemann problem on the compressible side with the help of velocity and pressure jump conditions across the interface. At the interface location, the flux terms for the compressible flow solver are provided by the coupling algorithms. The coupling is demonstrated in a one‐dimensional framework by use of a discontinuous Galerkin scheme for compressible two‐phase flow with a sharp interface tracking via a ghost‐fluid type method. The coupling schemes are applied to two generic test cases. The computational results are compared with those obtained with the fully compressible two‐phase flow solver, where the Mach number zero limit is approached by a weakly compressible fluid. For all cases, we obtain a very good agreement between the coupling approaches and the fully compressible solver. Copyright © 2014 John Wiley & Sons, Ltd.     

Non‐linear stability and convection for laminar flows in a porous medium with Brinkman law

The non‐linear stability of plane parallel shear flows in an incompressible homogeneous fluid heated from below and saturating a porous medium is studied by the Lyapunov direct method.In the Oberbeck–Boussinesq–Brinkman (OBB) scheme, if the inertial terms are negligible, as it is widely assumed in literature, we find global non‐linear exponential stability (GNES) independent of the Reynolds number R. However, if these terms are retained, we find a restriction on R (depending on the inertial convective coefficient) both for a homogeneous fluid and a mixture heated and salted from below. In the case of a mixture, when the normalized porosity ε is equal to one, the laminar flows are GNES for small R and for heat Rayleigh numbers less than the critical Rayleigh numbers obtained for the motionless state. Copyright © 2003 John Wiley & Sons, Ltd.     

Rotating compressible fluids under strong stratification

We consider the Navier–Stokes system written in the rotational frame describing the motion of a compressible viscous fluid under strong stratification. The asymptotic limit for low Mach and Rossby numbers and large Reynolds number is studied on condition that the Froude number characterizing the degree of stratification is proportional to the Mach number. We show that, at least for the well prepared data, the limit system is the same as for the problem without stratification—a variant of the incompressible planar Euler system.     

Numerical solution of transonic and subsonic flow of compressible inviscid and viscous fluid

The work presents two numerical solutions of compressible flows problems with high and very low Mach numbers. Both problems are numerically solved by finite volume method and the explicit MacCormack scheme using a grid of quadrilateral cells. Moved grid of quadrilateral cells is considered in the form of conservation laws using Arbitrary Lagrangian–Eulerian method. In the first case, inviscid transonic flow through cascade DCA 8% is presented and the numerical results are compared to experimental data. The second case, numerical solution of unsteady viscous flow in the channel for upstream Mach number M_∞=0.012 and frequency of the wall motions 100 Hz is presented. The unsteady case can represent a simplified model of airflow coming from the trachea, through the glottal region with periodically vibrating vocal folds to the human vocal tract. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A VOF Method for Isobaric Expandable Flows

In this paper a method is described to simulate isobaric multiphase flows at low Mach numbers with steep temperature gradients for fluids with non-negligible thermal expansivity. Governing equations and solution procedure are outlined. Further, a test case is shown in order to verify the model. Single phase natural convection flows with large temperature differences and either constant or temperature dependent transport properties were simulated to prove the solution of coupled Navier-Stokes and energy equations. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A direct Schur–Fourier decomposition for the efficient solution of high‐order Poisson equations on loosely coupled parallel computers

In this paper a parallel direct Schur–Fourier decomposition (DSFD) algorithm for the direct solution of arbitrary order discrete Poisson equations on parallel computers is proposed. It is based on a combination of a Direct Schur method and a Fourier decomposition and allows to solve each Poisson equation almost to machine accuracy using only one communication episode. Thus, it is well suited for loosely coupled parallel computers, that have a high network latency compared with the CPU performance. Several three‐dimensional direct numerical simulations (DNS) of wall‐bounded turbulent incompressible flows have been carried out using the DSFD algorithm. Numerical examples illustrating the robustness and scalability of the method on a PC cluster with a conventional 100 Mbits/s network are also presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Kinetic model of the Boltzmann equation for a diatomic gas with rotational degrees of freedom

A system of model kinetic equations is proposed to describe flows of a diatomic rarefied gas (nitrogen). A conservative numerical method is developed for its solution. A shock wave structure in nitrogen is computed, and the results are compared with experimental data in a wide range of Mach numbers. The system of model kinetic equations is intended to compute complex-geometry three-dimensional flows of a diatomic gas with rotational degrees of freedom.     

Zero Mach number limit in critical spaces for compressible Navier-Stokes equations

We are concerned with the existence and uniqueness of local or global solutions for slightly compressible viscous fluids in the whole space. In [6] and [7], we proved local and global well-posedness results for initial data in critical spaces very close to the one used by H. Fujita and T. Kato for incompressible flows (see [14]). In the present paper, we address the question of convergence to the incompressible model (for ill-prepared initial data) when the Mach number goes to zero. When the initial data are small in a critical space, we get global existence and convergence. For large initial data and a bit of additional regularity, the slightly compressible solution is shown to exist as long as the corresponding incompressible solution does. As a corollary, we get global existence (and uniqueness) for slightly compressible two-dimensional fluids.     

On Low Mach Number Preconditioning of Finite Volume Schemes

A finite volume method for inviscid unsteady flows at low Mach numbers is studied. The method uses a preconditioning of the dissipation term within the numerical flux function only. It can be observed by numerical experiments, as well as by analysis, that the preconditioned scheme yields a physically corrected pressure distribution and combined with an explicit time integrator it is stable if the time step Δt satisfies the requirement to be 𝒪(M ²) as the Mach number M tends to zero, whereas the corresponding standard method remains stable up to Δt = 𝒪(M ),M → 0, though producing unphysical results. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)