The Incompressible Limit of the Non-Isentropic Euler Equations

We study the Euler equations for slightly compressible fluids, that is, after rescaling, the limits of the Euler equations of fluid dynamics as the Mach number tends to zero. In this paper, we consider the general non-isentropic equations and general data. We first prove the existence of classical solutions for a time independent of the small parameter. Then, on the whole space ℝ ^d , we prove that the solution converges to the solution of the incompressible Euler equations. Accepted December 1, 2000?Published online April 23, 2001     

Singular Limits of the Equations of Magnetohydrodynamics

This paper studies the asymptotic limit for solutions to the equations of magnetohydrodynamics, specifically, the Navier–Stokes–Fourier system describing the evolution of a compressible, viscous, and heat conducting fluid coupled with the Maxwell equations governing the behavior of the magnetic field, when Mach number and Alfvén number tends to zero. The introduced system is considered on a bounded spatial domain in \mathbbR³{\mathbb{R}^{3}}, supplemented with conservative boundary conditions. Convergence towards the incompressible system of the equations of magnetohydrodynamics is shown.     

The Fully Compressible Semi-Geostrophic System from Meteorology

The semi-geostrophic equations are an approximation to the 3-D Euler equations for an atmosphere where the effects of rotation dominate. They are used by meteorologists to model the formation of fronts. Mathematically rigorous results were obtained by Benamou and Brenier and by Cullen and Gangbo in the incompressible case. In this paper we extend the results of Benamou and Brenier to the fully compressible case and show the existence of weak solutions to a reformulation of these equations in so-called dual variables. The Monge-Kantorovich theory appears as a material tool in our study. (Accepted September 23, 2002) Published online March 6, 2003 Communicated by Y. Brenier     

Fourth‐order compact formulation of Navier–Stokes equations and driven cavity flow at high Reynolds numbers

A new fourth‐order compact formulation for the steady 2‐D incompressible Navier–Stokes equations is presented. The formulation is in the same form of the Navier–Stokes equations such that any numerical method that solve the Navier–Stokes equations can easily be applied to this fourth‐order compact formulation. In particular, in this work the formulation is solved with an efficient numerical method that requires the solution of tridiagonal systems using a fine grid mesh of 601 × 601. Using this formulation, the steady 2‐D incompressible flow in a driven cavity is solved up to Reynolds number with Re = 20 000 fourth‐order spatial accuracy. Detailed solutions are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

A manufactured solution for a two-dimensional steady wall-bounded incompressible turbulent flow

This paper presents a manufactured solution (MS), resembling a two-dimensional, steady, wall-bounded, incompressible, turbulent flow for RANS codes verification. The specified flow field satisfies mass conservation, but requires additional source terms in the momentum equations. To also allow verification of the correct implementation of the turbulence models transport equations, the proposed MS exhibits most features of a true near-wall turbulent flow. The model is suited for testing six eddy-viscosity turbulence models: the one-equation models of Spalart and Allmaras and Menter; the standard two-equation k–ε model and the low-Reynolds version proposed by Chien; the TNT and BSL versions of the k–ω model.     

Time‐accurate intermediate boundary conditions for large eddy simulations based on projection methods

Projection methods are among the most adopted procedures for solving the Navier–Stokes equations system for incompressible flows. In order to simplify the numerical procedures, the pressure–velocity de‐coupling is often obtained by adopting a fractional time‐step method. In a specific formulation, suitable for the incompressible flows equations, it is based on a formal decomposition of the momentum equation, which is related to the Helmholtz–Hodge Decomposition theorem of a vector field in a finite domain. Owing to the continuity constraint also in large eddy simulation of turbulence, as happens for laminar solutions, the filtered pressure characterizes itself only as a Lagrange multiplier, not a thermodynamic state variable. The paper illustrates the implications of adopting such procedures when the decoupling is performed onto the filtered equations system. This task is particularly complicated by the discretization of the time integral of the sub‐grid scale tensor. A new proposal for developing time‐accurate and congruent intermediate boundary conditions is addressed. Several tests for periodic and non‐periodic channel flows are presented. This study follows and completes the previous ones reported in (Int. J. Numer. Methods Fluids 2003; 42, 43 ). Copyright © 2005 John Wiley & Sons, Ltd.     

Local pressure preconditioning method for steady incompressible flows

The convergence and accuracy characteristics of the preconditioned incompressible Euler and Navier–Stokes equations are studied. An object-oriented C++ numerical code has been developed for solving the inviscid and viscous, steady, incompressible flows problems. The code is based on the cell-centred finite volume method. In this scheme, two-dimensional incompressible Euler and Navier–Stokes equations are modified by a robust artificial compressibility (AC) and a local preconditioning matrix of pressure-sensor type. The preconditioned equations are solved with the Jameson's numerical approach, i.e. artificial dissipation and artificial viscosity terms under the form of a fourth- and second-order derivative, respectively. An explicit four-stage Runge–Kutta integration algorithm is applied to obtain the steady-state condition. The computed results include the steady-state solution of flow past the NACA-hydrofoils and a circular cylinder in free stream, for which the numerical results are compared with numerical works of other researchers. Good agreement is observed. The effects of AC parameter, artificial viscosity and dissipation factor, and local preconditioning coefficient on convergence rate and solution accuracy are tested by computing flow over the NACA0012 hydrofoil. In addition, some important design criteria of a preconditioner, such as stiffness reduction, hyperbolicity, symmetrisability, accuracy preservation for M → 0, and M-property have been examined analytically.     

On Admissibility Criteria for Weak Solutions of the Euler Equations

We consider solutions to the Cauchy problem for the incompressible Euler equations satisfying several additional requirements, like the global and local energy inequalities. Using some techniques introduced in an earlier paper, we show that, for some bounded compactly supported initial data, none of these admissibility criteria singles out a unique weak solution. As a byproduct, in more than one space dimension, we show bounded initial data for which admissible solutions to the p-system of isentropic gas dynamics in Eulerian coordinates are not unique.     

Higher order spectral/hp finite element models of the Navier–Stokes equations

In this paper, we present higher order least-squares finite element formulations for viscous, incompressible, isothermal Navier–Stokes equations using spectral/hp basis functions. The second-order Navier–Stokes equations are recast as first-order system of equations using stresses as auxiliary variables. Both steady-state and transient problems are considered. For a better coupling of pressure and velocity, especially in transient flows, an iterative penalisation strategy is employed. The outflow-type boundary conditions are applied in a weak sense through the least-squares functional. The formulation is verified by solving various benchmark problems like the lid-driven cavity, backward-facing step and flow over cylinder problems using direct serial solver UMFPACK.     

Applications of the artificial compressibility method for turbulent open channel flows

A three‐dimensional (3‐D) numerical method for solving the Navier–Stokes equations with a standard k–ε turbulence model is presented. In order to couple pressure with velocity directly, the pressure is divided into hydrostatic and hydrodynamic parts and the artificial compressibility method (ACM) is employed for the hydrodynamic pressure. By introducing a pseudo‐time derivative of the hydrodynamic pressure into the continuity equation, the incompressible Navier–Stokes equations are changed from elliptic‐parabolic to hyperbolic‐parabolic equations. In this paper, a third‐order monotone upstream‐centred scheme for conservation laws (MUSCL) method is used for the hyperbolic equations. A system of discrete equations is solved implicitly using the lower–upper symmetric Gauss–Seidel (LU‐SGS) method. This newly developed numerical method is validated against experimental data with good agreement. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical computation of an impinging jet with uniform wall suction

The paper presents numerical predictions of a turbulent axisymmetric jet impinging onto a porous plate, based on a finite volume method of solving the Navier-Stokes equations for an incompressible air jet with the K–ε turbulence model. The velocity and pressure terms of the momentum equations are solved by the SIMPLE (semi-implicit method for pressure-linked equation) method. In this study, non-uniform staggered grids are used. The parameters of interest include the nozzle-to-wall distance and the suction velocity. The results of the present calculations are compared with available data reported in the literature. It is found that suction effects reduce the boundary layer thickness and increase the velocity gradient near the wall.     

Explicit formulations for evaluation of velocity gradients using boundary‐domain integral equations in 2D and 3D viscous flows

In this paper, explicit boundary‐domain integral equations for evaluating velocity gradients are derived from the basic velocity integral equations. A free term is produced in the new strongly singular integral equation, which is not included in recent formulations using the complex variable differentiation method (CVDM) to compute velocity gradients (Int. J. Numer. Meth. Fluids 2004; 45 :463–484; Int. J. Numer. Meth. Fluids 2005; 47 :19–43). The strongly singular domain integrals involved in the new integral equations are accurately evaluated using the radial integration method (RIM). Considerable computational time for evaluating integrals of velocity gradients can be saved by using present formulation than using CVDM. The formulation derived in this paper together with those presented in reference (Int. J. Numer. Meth. Fluids 2004; 45 :463–484) for 2D and in (Int. J. Numer. Meth. Fluids 2005; 47 :19–43) for 3D problems constitutes a complete boundary‐domain integral equation system for solving full Navier–Stokes equations using primitive variables. Three numerical examples for steady incompressible viscous flow are given to validate the derived formulations. Copyright © 2007 John Wiley & Sons, Ltd.     

Incompressible SPH simulation of wave breaking and overtopping with turbulence modelling

In this paper a truly incompressible version of the smoothed particle hydrodynamics (SPH) method is presented to investigate the surface wave overtopping. SPH is a pure Lagrangian approach which can handle large deformations of the free surface with high accuracy. The governing equations are solved based on the SPH particle interaction models and the incompressible algorithm of pressure projection is implemented by enforcing the constant particle density. The two‐equation k–ε model is an effective way of dealing with the turbulence and vortices during wave breaking and overtopping and it is coupled with the incompressible SPH numerical scheme. The SPH model is employed to reproduce the experiment and computations of wave overtopping of a sloping sea wall. The computations are validated against the experimental and numerical data found in the literatures and good agreement is observed. Besides, the convergence behaviour of the numerical scheme and the effects of particle spacing refinement and turbulence modelling on the simulation results are also investigated in further detail. The sensitivity of the computed wave breaking and overtopping on these issues is discussed and clarified. Copyright © 2005 John Wiley & Sons, Ltd.     

Equivalence of u–p and ζ–ψ formulations of the time-dependent Navier–Stokes equations

This paper deals with the non-stationary incompressible Navier--Stokes equations for two-dimensional flows expressed in terms of the velocity and pressure and of the vorticity and streamfunction. The equivalence of the two formulations is demonstrated, both formally and rigorously, by virtue of a condition of compatibility between the boundary and initial values of the normal component of velocity. This condition is shown to be the only compatibility condition necessary to allow for solutions of a minimal regularity, namely H¹ for the velocity, as in most current numerical schemes relying on spatial discretizations of local type.     

A locally conservative,discontinuous least‐squares finite element method for the Stokes equations

Conventional least‐squares finite element methods (LSFEMs) for incompressible flows conserve mass only approximately. For some problems, mass loss levels are large and result in unphysical solutions. In this paper we formulate a new, locally conservative LSFEM for the Stokes equations wherein a discrete velocity field is computed that is point‐wise divergence free on each element. The central idea is to allow discontinuous velocity approximations and then to define the velocity field on each element using a local stream‐function. The effect of the new LSFEM approach on improved local and global mass conservation is compared with a conventional LSFEM for the Stokes equations employing standard C⁰ Lagrangian elements. Copyright © 2011 John Wiley & Sons, Ltd.     

Decay in Time of Incompressible Flows

In this paper we consider the Cauchy problem for incompressible flows governed by the Navier-Stokes or MHD equations. We give a new proof for the time decay of the spatial L₂ L_2 norm of the solution, under the assumption that the solution of the heat equation with the same initial data decays. By first showing decay of the first derivatives of the solution, we avoid some technical difficulties of earlier proofs based on Fourier splitting.     

Similarity in stationary motions of gas and two-phase medium with incompressible component

In this paper we suggest the transformation between the equations for a perfect gas and the equations describing in one-velocity approach the two-phase medium with any volume occupied by the incompressible phase. It is proved that the motion of a two-phase medium in the transformed coordinate system is similar with certain accuracy to that of a perfect gas. It means that the solutions obtained for perfect gas can be used to solve wave problems for media with incompressible component. There is no necessity directly to solve the problem for medium with incompressible component, and it is only sufficient to transform the known solution of the similar problem for a homogeneous medium. Thus, the solutions of many hydrodynamic problems for multi-component media with incompressible phase can be obtained without solving the original set of equations. The scope for the suggested transformation is demonstrated by reference to the strong explosion in a two-phase medium.     

A promising boundary element formulation for three‐dimensional viscous flow

In this paper, a new set of boundary‐domain integral equations is derived from the continuity and momentum equations for three‐dimensional viscous flows. The primary variables involved in these integral equations are velocity, traction, and pressure. The final system of equations entering the iteration procedure only involves velocities and tractions as unknowns. In the use of the continuity equation, a complex‐variable technique is used to compute the divergence of velocity for internal points, while the traction‐recovery method is adopted for boundary points. Although the derived equations are valid for steady, unsteady, compressible, and incompressible problems, the numerical implementation is only focused on steady incompressible flows. Two commonly cited numerical examples and one practical pipe flow problem are presented to validate the derived equations. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical studies of slow viscous rotating flow past a sphere. II

The Navier–Stokes equations, which are the governing equations for a steady, viscous, incompressible fluid rotating about the z-axis with angular velocity ω, are linearized using the Oseen approximation. Two parameters, namely the Reynolds number Re = Ua/v and Re_ω = 2ωa²/v (the Reynolds number w.r.t. rotation), enter the linearized equations. These equations are solved by the Peaceman–Rachford ADI method and the resulting algebraic equations are solved by the SOR method. Streamlines are plotted and compared with the Oseen solution for the non-rotating case. The magnitude of the vorticity vector with increasing θ is also plotted.     

A high-order compact scheme for solving the 2D steady incompressible Navier-Stokes equations in general curvilinear coordinates

In this paper, a high-order compact finite difference algorithm is established for the stream function-velocity formulation of the two-dimensional steady incompressible Navier-Stokes equations in general curvilinear coordinates. Different from the previous work, not only the stream function and its first-order partial derivatives but also the second-order mixed partial derivative is treated as unknown variable in this work. Numerical examples, including a test problem with an analytical solution, three types of lid-driven cavity flow problems with unusual shapes and steady flow past a circular cylinder as well as an elliptic cylinder with angle of attack, are solved numerically by the newly proposed scheme. For two types of the lid-driven trapezoidal cavity flow, we provide the detailed data using the fine grid sizes, which can be considered the benchmark solutions. The results obtained prove that the present numerical method has the ability to solve the incompressible flow for complex geometry in engineering applications, especially by using a nonorthogonal coordinate transformation, with high accuracy.