Effect of bulk viscosity in supersonic flow past spacecraft

In this paper, we consider the effect of bulk viscosity in various hydrodynamic problems. We numerically study this effect on the front structure of the one-dimensional stationary shock wave and on the flow past blunt body. We estimate the effect of the bulk viscosity coefficient (BVC) on the heat transfer and drag of a sphere in a supersonic flow, apparently for the first time, by the numerical solution of parabolized Navier–Stokes equations. The solution is obtained by an original fast convergent method of global iterations of the longitudinal pressure gradient. The directions of further investigations of bulk viscosity are suggested.     

Zero-viscosity-capillarity limit to the planar rarefaction wave for the 2D compressible Navier–Stokes–Korteweg equations

We shall consider the two-dimensional (2D) isentropic Navier–Stokes–Korteweg equations which are used to model compressible fluids with internal capillarity. Formally, the 2D isentropic Navier–Stokes–Korteweg equations converge, as the viscosity and the capillarity vanish, to the corresponding 2D inviscid Euler equations, and we do justify this for the case that the corresponding 2D inviscid Euler equations admit a planar rarefaction wave solution. More precisely, it is proved that there exists a family of smooth solutions for the 2D isentropic compressible Navier–Stokes–Korteweg equations converging to the planar rarefaction wave solution with arbitrary strength for the 2D Euler equations. A uniform convergence rate is obtained in terms of the viscosity coefficient and the capillarity away from the initial time. The key ingredients of our proof are the re-scaling technique and energy estimate, in which we also introduce the hyperbolic wave to recover the physical viscosities and capillarity of the inviscid rarefaction wave profile.     

An iterative solver for the Oseen problem and numerical solution of incompressible Navier–Stokes equations

Incompressible unsteady Navier–Stokes equations in pressure–velocity variables are considered. By use of the implicit and semi‐implicit schemes presented the resulting system of linear equations can be solved by a robust and efficient iterative method. This iterative solver is constructed for the system of linearized Navier–Stokes equations. The Schur complement technique is used. We present a new approach of building a non‐symmetric preconditioner to solve a non‐symmetric problem of convection–diffusion and saddle‐point type. It is shown that handling the differential equations properly results in constructing efficient solvers for the corresponding finite linear algebra systems. The method has good performance for various ranges of viscosity and can be used both for 2D and 3D problems. The analysis of the method is still partly heuristic, however, the mathematically rigorous results are proved for certain cases. The proof is based on energy estimates and basic properties of the underlying partial differential equations. Numerical results are provided. Additionally, a multigrid method for the auxiliary convection–diffusion problem is briefly discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical Simulation of the Flow over a Segment-Conical Body on the Basis of Reynolds Equations

Numerical simulation was used to study the 3D supersonic flow over a segment-conical body similar in shape to the ExoMars space vehicle. The nonmonotone behavior of the normal force acting on the body placed in a supersonic gas flow was analyzed depending on the angle of attack. The simulation was based on the numerical solution of the unsteady Reynolds-averaged Navier–Stokes equations with a two-parameter differential turbulence model. The solution of the problem was obtained using the in-house solver HSFlow with an efficient parallel algorithm intended for multiprocessor super computers.     

Comparison of computations of asymptotic flow models in a constricted channel

We aim at comparing computations with asymptotic models issued from incompressible Navier–Stokes at high Reynolds number: the Reduced Navier–Stokes/Prandtl (RNS/P) equations and the Double Deck (DD) equations. We treat the case of the steady two dimensional flow in a constricted pipe. In particular, finite differences and finite element solvers are compared for the RNS/P equations. It results from this study that the two codes compare well. Numerical examples also illustrate the interest of these asymptotic models as well as the flexibility of the finite element solver.     

Vanishing viscosity and Debye-length limit to rarefaction wave with vacuum for the 1D bipolar Navier–Stokes–Poisson equation

In this paper, we consider the one-dimensional (1D) compressible bipolar Navier–Stokes–Poisson equations. We know that when the viscosity coefficient and Debye length are zero in the compressible bipolar Navier–Stokes–Poisson equations, we have the compressible Euler equations. Under the case that the compressible Euler equations have a rarefaction wave with one-side vacuum state, we can construct a sequence of the approximation solution to the one-dimensional bipolar Navier–Stokes–Poisson equations with well-prepared initial data, which converges to the above rarefaction wave with vacuum as the viscosity and the Debye length tend to zero. Moreover, we also obtain the uniform convergence rate. The results are proved by a scaling argument and elaborate energy estimate.     

On stationary two‐dimensional flows around a fast rotating disk

We study the two‐dimensional stationary Navier–Stokes equations describing flows around a rotating disk. The existence of unique solutions is established for any rotating speed, and qualitative effects of a large rotation are described precisely by exhibiting a boundary layer structure and an axisymmetrization of the flow.     

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.     

Regularity of the 3D Navier–Stokes equations with viewpoint of 2D flow

The regularity of 2D Navier–Stokes flow is well known. In this article we study the relationship of 3D and 2D flow, and the regularity of the 3D Naiver–Stokes equations with viewpoint of 2D equations. We consider the problem in the Cartesian and in the cylindrical coordinates.     

Retracted: Unsteady flow and heat transfer on a rotating disk in a porous medium

This article has been retracted. See retraction notice DOI: 10.1002/mma.850 . An unsteady flow and heat transfer in a porous medium of a viscous incompressible fluid over a rotating disk in an otherwise ambient fluid are studied. The unsteadiness in the flow field is caused by the angular velocity of the disk which varies with time. The new self‐similar solution of the Navier–Stokes and energy equations is obtained numerically. The solution obtained here is not only the solution of the Navier–Stokes equations, but also of the boundary layer equations. Also, for a simple scaling factor, it represents the solution of the flow and heat transfer in the forward stagnation‐point region of a rotating sphere or over a rotating cone. The asymptotic behaviour of the solution for a large porosity or for a large independent variable is also examined. The surface shear stresses in the radial and tangential directions and the surface heat transfer increase as the acceleration parameter increases. Also, the surface shear stress in the radial direction and the surface heat transfer decrease with increasing porosity, but the surface shear stress in the tangential direction increases. Copyright © 2006 John Wiley & Sons, Ltd.     

On the regularity of the solutions to the 3D Navier–Stokes equations: a remark on the role of the helicity

We show that if velocity and vorticity are orthogonal at each point (and they become orthogonal fast enough) then solutions of the 3D Navier–Stokes equations are smooth. This condition implies that the helicity is identically zero and, in a certain sense, the flow resembles the 2D geometric situation. To cite this article: L.C. Berselli, D. Córdoba, C. R. Acad. Sci. Paris, Ser. I 347 (2009).     

The vanishing viscosity as a selection principle for the Euler equations: The case of 3D shear flow

We show that for a certain family of initial data, there exist non-unique weak solutions to the 3D incompressible Euler equations satisfying the weak energy inequality, whereas the weak limit of every sequence of Leray–Hopf weak solutions for the Navier–Stokes equations, with the same initial data, and as the viscosity tends to zero, is uniquely determined and equals the shear flow solution of the Euler equations corresponding to this initial data. This simple example suggests that, also in more general situations, the vanishing viscosity limit of the Navier–Stokes equations could serve as a uniqueness criterion for weak solutions of the Euler equations.     

Diffusive wave in the low Mach limit for non-viscous and heat-conductive gas

The low Mach number limit for one-dimensional non-isentropic compressible Navier–Stokes system without viscosity is investigated, where the density and temperature have different asymptotic states at far fields. It is proved that the solution of the system converges to a nonlinear diffusion wave globally in time as Mach number goes to zero. It is remarked that the velocity of diffusion wave is proportional with the variation of temperature. Furthermore, it is shown that the solution of compressible Navier–Stokes system also has the same phenomenon when Mach number is suitably small.     

A discontinuous finite element solution of the Boltzmann kinetic equation in collisionless and BGK forms for macroscopic gas flows

In this paper, an alternative approach to the traditional continuum analysis of flow problems is presented. The traditional methods, that have been popular with the CFD community in recent times, include potential flow, Euler and Navier–Stokes solvers. The method presented here involves solving the governing equation of the molecular gas dynamics that underlies the macroscopic behaviour described by the macroscopic governing equations. The equation solved is the Boltzmann kinetic equation in its simplified collisionless and BGK forms. The algorithm used is a discontinuous Taylor–Galerkin type and it is applied to the 2D problems of a highly rarefied gas expanding into a vacuum, flow over a vertical plate, rarefied hypersonic flow over a double ellipse, and subsonic and transonic flow over an aerofoil. The benefit of this type of solver is that it is not restricted to continuum regime (low Knudsen number) problems. However, it is a computationally expensive technique.     

On singularity formation of a 3D model for incompressible Navier–Stokes equations

We investigate the singularity formation of a 3D model that was recently proposed by Hou and Lei (2009) in [15] for axisymmetric 3D incompressible Navier–Stokes equations with swirl. The main difference between the 3D model of Hou and Lei and the reformulated 3D Navier–Stokes equations is that the convection term is neglected in the 3D model. This model shares many properties of the 3D incompressible Navier–Stokes equations. One of the main results of this paper is that we prove rigorously the finite time singularity formation of the 3D inviscid model for a class of initial boundary value problems with smooth initial data of finite energy. We also prove the global regularity of the 3D inviscid model for a class of small smooth initial data.     

Global well-posedness of the velocity–vorticity-Voigt model of the 3D Navier–Stokes equations

The velocity–vorticity formulation of the 3D Navier–Stokes equations was recently found to give excellent numerical results for flows with strong rotation. In this work, we propose a new regularization of the 3D Navier–Stokes equations, which we call the 3D velocity–vorticity-Voigt (VVV) model, with a Voigt regularization term added to momentum equation in velocity–vorticity form, but with no regularizing term in the vorticity equation. We prove global well-posedness and regularity of this model under periodic boundary conditions. We prove convergence of the model's velocity and vorticity to their counterparts in the 3D Navier–Stokes equations as the Voigt modeling parameter tends to zero. We prove that the curl of the model's velocity converges to the model vorticity (which is solved for directly), as the Voigt modeling parameter tends to zero. Finally, we provide a criterion for finite-time blow-up of the 3D Navier–Stokes equations based on this inviscid regularization.     

Pontryagin’s Principle for Optimal Control Problem Governed by 3D Navier–Stokes Equations

This paper deals with the Pontryagin maximum principle for optimal control problems governed by 3D Navier–Stokes equations with pointwise control constraint. The obtained result is proved by using some results on regularity of solutions of the Navier–Stokes equations and techniques of optimal control theory.     

A factorization method for numerical solution of the Navier–Stokes equations for a viscous incompressible liquid

Some implicit difference scheme of approximate factorization is proposed for numerical solution of the Navier–Stokes equations for an incompressible liquid in curvilinear coordinates. Testing of the algorithm is carried out on the solution of the problems concerning the Couette and Poiseuille flows; and the results are presented of numerical simulation of a flow between the rotating cylinders with covers.     

Navier–Stokes equations with periodic boundary conditions and pressure loss

We present in this note the existence and uniqueness results for the Stokes and Navier–Stokes equations which model the laminar flow of an incompressible fluid inside a two-dimensional channel of periodic sections. The data of the pressure loss coefficient enables us to establish a relation on the pressure and to thus formulate an equivalent problem.     

Stability of the superposition of rarefaction wave and contact discontinuity for the non‐isentropic Navier–Stokes–Poisson system

This paper is devoted to the study of the nonlinear stability of the composite wave consisting of a rarefaction wave and a viscous contact discontinuity wave of the non‐isentropic Navier–Stokes–Poisson system with free boundary. We first construct the composite wave through the quasineutral Euler equations and then prove that the composite wave is time asymptotically stable under small perturbations for the corresponding initial‐boundary value problem of the non‐isentropic Navier–Stokes–Poisson system. Only the strength of the viscous contact wave is required to be small. However, the strength of the rarefaction wave can be arbitrarily large. In our analysis, the domain decomposition plays an important role in obtaining the zero‐order energy estimates. By introducing this technique, we successfully overcome the difficulty caused by the critical terms involved with the linear term, which does not satisfy the quasineural assumption for the composite wave. Copyright © 2016 John Wiley & Sons, Ltd.