Existence and uniqueness analysis of a detached shock problem for the potential flow

We study a problem for two-dimensional steady potential and isentropic Euler equations in a bounded domain, where an artificial detached shock interacts with a wedge. Using the stream function, we obtain a free boundary problem for the subsonic state and the detached artificial shock curve and we prove that such configuration admits a unique solution in certain weighted Hölder spaces. The proof is based on various Hölder and Schauder estimates for second-order elliptic equations and fixed point theorems. Moreover, we pose an energy principle and remark that the physical attached shock is the minimizer of the energy functional.     

Asymptotic stability of a composite wave of two viscous shock waves for the one-dimensional radiative Euler equations

This paper is devoted to the study of the wellposedness of the radiative Euler equations. By employing the anti-derivative method, we show the unique global-in-time existence and the asymptotic stability of the solutions of the radiative Euler equations for the composite wave of two viscous shock waves with small strength. This method developed here is also helpful to other related problems with similar analytical difficulties.     

Two-dimensional travelling waves over a moving bottom

Two-dimensional travelling waves on an ideal fluid with gravity and surface tension over a periodically moving bottom with a small amplitude are studied. The bottom and the wave travel with a same speed. The exact Euler equations are formulated as a spatial dynamic system by using the stream function. A manifold reduction technique is applied to reduce the system into one of ordinary differential equations with finite dimensions. A homoclinic solution to the normal form of this reduced system persists when higher-order terms are added, which gives a generalized solitary wave—the homoclinic solution connecting a periodic solution.     

Cauchy problem with general discontinuous initial data along a smooth curve for 2-d Euler system

We study the Cauchy problems for the isentropic 2-d Euler system with discontinuous initial data along a smooth curve. All three singularities are present in the solution: shock wave, rarefaction wave and contact discontinuity. We show that the usual restrictive high order compatibility conditions for the initial data are automatically satisfied. The local existence of piecewise smooth solution containing all three waves is established.     

Nanofluid flow and heat transfer in a Brinkman porous channel with variable porosity

The problem of forced convection in a channel filled with a nanofluidsaturated porous medium is investigated, numerically. A finite difference Computational Fluid Dynamics (CFD) model with structured uniform grid system is employed to solve the momentum and energy equations. In modeling flow in the channel, the effects of flow inertia, variable porosity and Brinkman friction are taken into account. Studies are carried out for different nanoparticles with different volume fractions in the range 0%-4% and different nanoparticle diameters. Comparison made between our numerical and semi analytical Differential Transform Method (DTM) results with those in previous published research is found to be appropriate. Results show that increasing either nanoparticls volume fraction or pressure gradient parameter improves heat transfer. Further, for large quantities of nanoparticle concentration and pressure gradient, the channeling phenomenon is intensified.     

Global attractors for singular perturbations of the Cahn-Hilliard equations

We consider the singular perturbations of two boundary value problems, concerning respectively the viscous and the nonviscous Cahn-Hilliard equations in one dimension of space. We show that the dynamical systems generated by these two problems admit global attractors in the phase space , and that these global attractors are at least upper-semicontinuous with respect to the vanishing of the perturbation parameter.     

Node-pair finite volume/finite element schemes for the Euler equation in cylindrical and spherical coordinates

A numerical scheme is presented for the solution of the compressible Euler equations in both cylindrical and spherical coordinates. The unstructured grid solver is based on a mixed finite volume/finite element approach. Equivalence conditions linking the node-centered finite volume and the linear Lagrangian finite element scheme over unstructured grids are reported and used to devise a common framework for solving the discrete Euler equations in both the cylindrical and the spherical reference systems. Numerical simulations are presented for the explosion and implosion problems with spherical symmetry, which are solved in both the axial–radial cylindrical coordinates and the radial–azimuthal spherical coordinates. Numerical results are found to be in good agreement with one-dimensional simulations over a fine mesh.     

The lifespan of spherically symmetric solutions of the compressible Euler equations outside an impermeable sphere

We consider smooth three-dimensional spherically symmetric Eulerian flows of ideal polytropic gases outside an impermeable sphere, with initial data equal to the sum of a constant flow with zero velocity and a smooth perturbation with compact support. Under a natural assumption on the form of the perturbation, we obtain precise information on the asymptotic behavior of the lifespan as the size of the perturbation tends to 0. When there is no sphere, so that the flow is defined in all space, corresponding results have been obtained in [P. Godin, The lifespan of a class of smooth spherically symmetric solutions of the compressible Euler equations with variable entropy in three space dimensions, Arch. Ration. Mech. Anal. 177 (2005) 479–511].     

Convergence and stability of finite element nonlinear Galerkin method for the Navier-Stokes equations

In this article, we present a new fully discrete finite element nonlinear Galerkin method, which are well suited to the long time integration of the Navier-Stokes equations. Spatial discretization is based on two-grid finite element technique; time discretization is based on Euler explicit scheme with variable time step size. Moreover, we analyse the boundedness, convergence and stability condition of the finite element nonlinear Galerkin method. Our discussion shows that the time step constraints of the method depend only on the coarse grid parameter and the time step constraints of the finite element Galerkin method depend on the fine grid parameter under the same convergence accuracy. Received February 2, 1994 / Revised version received December 6, 1996     

Conceptual linearization of Euler governing equations to solve high speed compressible flow using a pressure‐based method

The main objective of the current work is to introduce a new conceptual linearization strategy to improve the performance of a primitive shock‐capturing pressure‐based finite‐volume method. To avoid a spurious oscillatory solution in the chosen collocated grids, both the primitive and extended methods utilize two convecting and convected momentum expressions at each cell face. The expressions are obtained via a physical‐based discretization of two inclusive statements, which are constructed via a novel incorporation of the continuity and momentum governing equations. These two expressions in turn provide a strong coupling among the Euler conservative statements. Contrary to the primitive work, the linearization in the current work respects the definitions and essence of physics behind deriving the Euler governing equations. The accuracy and efficiency of the new formulation are then investigated by solving the shock tube as a problem with moving normal and expansion waves and the converging‐diverging nozzle as a problem with strong stationary normal shock. The results show that there is good improvement in performance of the primitive pressure‐based shock‐capturing method while its superior accuracy is not deteriorated at all. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Stability of steady multi-wave configurations for the full Euler equations of compressible fluid flow

We are concerned with the stability of steady multi-wave configurations for the full Euler equations of compressible fluid flow. In this paper, we focus on the stability of steady four-wave configurations that are the solutions of the Riemann problem in the flow direction, consisting of two shocks, one vortex sheet, and one entropy wave, which is one of the core multi-wave configurations for the two-dimensional Euler equations. It is proved that such steady four-wave configurations in supersonic flow are stable in structure globally, even under the BV perturbation of the incoming flow in the flow direction. In order to achieve this, we first formulate the problem as the Cauchy problem (initial value problem) in the flow direction, and then develop a modified Glimm difference scheme and identify a Glimm-type functional to obtain the required BV estimates by tracing the interactions not only between the strong shocks and weak waves, but also between the strong vortex sheet/entropy wave and weak waves. The key feature of the Euler equations is that the reflection coefficient is always less than $1$, when a weak wave of different family interacts with the strong vortex sheet/entropy wave or the shock wave, which is crucial to guarantee that the Glimm functional is decreasing. Then these estimates are employed to establish the convergence of the approximate solutions to a global entropy solution, close to the background solution of steady four-wave configuration.     

Existence and continuous dependence of large solutions for the magnetohydrodynamic equations

Global solutions of the nonlinear magnetohydrodynamic (MHD) equations with general large initial data are investigated. First the existence and uniqueness of global solutions are established with large initial data in H ¹. It is shown that neither shock waves nor vacuum and concentration are developed in a finite time, although there is a complex interaction between the hydrodynamic and magnetodynamic effects. Then the continuous dependence of solutions upon the initial data is proved. The equivalence between the well-posedness problems of the system in Euler and Lagrangian coordinates is also showed.     

Local well-posedness for the homogeneous Euler equations

We introduce Triebel-Lizorkin-Lorentz function spaces, based on the Lorentz L^p,q-spaces instead of the standard L^p-spaces, and prove a local-in-time unique existence and a blow-up criterion of solutions in those spaces for the Euler equations of inviscid incompressible fluid in Rⁿ,n≥2. As a corollary we obtain global existence of solutions to the 2D Euler equations in the Triebel-Lizorkin-Lorentz space. For the proof, we establish the Beale-Kato-Majda type logarithmic inequality and commutator estimates in our spaces. The key methods of proof used are the Littlewood-Paley decomposition and the paradifferential calculus by J.M. Bony.     

Strong trajectory attractors for dissipative Euler equations

The 2D Euler equations with periodic boundary conditions and extra linear dissipative term Ru, R>0 are considered and the existence of a strong trajectory attractor in the space is established under the assumption that the external forces have bounded vorticity. This result is obtained by proving that any solution belonging the proper weak trajectory attractor has a bounded vorticity which implies its uniqueness (due to the Yudovich theorem) and allows to verify the validity of the energy equality on the weak attractor. The convergence to the attractor in the strong topology is then proved via the energy method.     

The zero-electron-mass limit in the hydrodynamic model for plasmas

The limit of the vanishing ratio of the electron mass to the ion mass in the isentropic transient Euler-Poisson equations with periodic boundary conditions is proved. The equations consist of the balance laws for the electron density and current density for a given ion density, coupled to the Poisson equation for the electrostatic potential. The limit is related to the low-Mach-number limit of Klainerman and Majda. In particular, the limit velocity satisfies the incompressible Euler equations with damping. The difference to the zero-Mach-number limit comes from the electrostatic potential which needs to be controlled. This is done by a reformulation of the equations in terms of the enthalpy, higher-order energy estimates and a careful use of the Poisson equation.     

A finite volume scheme for a model coupling free surface and pressurised flows in pipes

A model is derived for the coupling of transient free surface and pressurized flows. The resulting system of equations is written under a conservative form with discontinuous gradient of pressure. We treat the transition point between the two types of flows as a free boundary associated to a discontinuity of the gradient of pressure. The numerical simulation is performed by making use of a Roe-like finite volume scheme that we adapted to such discontinuities in the flux. The validation is performed by comparison with experimental results.     

On blow-up criteria of smooth solutions to the 3-D Euler equations in a bounded domain

We show that a smooth solution of the 3-D Euler equations in a bounded domain breaks down, if and only if a certain norm of vorticity blows up at the same time. Here this norm is weaker than bmo-norm.     

The formulation of linear theory of a reflected shock in cylindrical geometry

In this paper we formulate the linear theory for compressible fluids in cylindrical geometry with small perturbation at the material interface. We derive the first order equations in the smooth regions, boundary conditions at the shock fronts and the contact interface by linearizing the Euler equations and Rankine-Hugoniot conditions. The small amplitude solution formulated in this paper will be important for calibration of results from full numerical simulation of compressible fluids in cylindrical geometry.     

Analysis on linear stability of oblique shock waves in steady supersonic flow

An attached oblique shock wave is generated when a sharp solid projectile flies supersonically in the air. We study the linear stability of oblique shock waves in steady supersonic flow under three dimensional perturbation in the incoming flow. Euler system of equations for isentropic gas model is used. The linear stability is established for shock front with supersonic downstream flow, in addition to the usual entropy condition.     

Global entropy solutions to the relativistic Euler equations for a class of large initial data

We are concerned with global entropy solutions to the relativistic Euler equations for a class of large initial data which involve the interaction of shock waves and rarefaction waves. We first carefully analyze the global behavior of the shock curves, the rarefaction wave curves, and their corresponding inverse curves in the phase plane. Based on these analyses, we use the Glimm scheme to construct global entropy solutions to the relativistic Euler equations for the class of large discontinuous initial data.Received: May 23, 2004