A fast riemann solver with constant covolume applied to the random choice method

The Riemann problem for the unsteady one-dimensional Euler equations together with the constant-covolume equation of state is solved exactly. The solution is then applied to the random choice method to solve the general initial-boundary value problem for the Euler equations. The iterative procedure to find p*, the pressure between the acoustic waves, involves a single algebraic (non-linear) equation, all other quantities follow directly throughout the x–t plane, except within rarefaction fans where an extra iterative procedure is required. The solution is validated against existing exact results both directly and in conjunction with the random choice method.     

Asymptotic Stability of Rarefaction Wave for the Navier–Stokes Equations for a Compressible Fluid in the Half Space

This paper is concerned with the asymptotic stability towards a rarefaction wave of the solution to an outflow problem for the Navier–Stokes equations in a compressible fluid in the Eulerian coordinate in the half space. This is the second one of our series of papers on this subject. In this paper, firstly we classify completely the time-asymptotic states, according to some parameters, that is the spatial-asymptotic states and boundary conditions, for this initial boundary value problem, and some pictures for the classification of time-asymptotic states are drawn in the state space. In order to prove the stability of the rarefaction wave, we use the solution to Burgers’ equation to construct a suitably smooth approximation of the rarefaction wave and establish some time-decay estimates in L ^p -norm for the smoothed rarefaction wave. We then employ the L ²-energy method to prove that the rarefaction wave is non-linearly stable under a small perturbation, as time goes to infinity. P. Zhu was supported by JSPS postdoctoral fellowship under P99217.     

Stochastic differential equations for the kinetic and hydrodynamic equations

A finite system of stochastic interacting particles is considered. The system approximates the solutions of the kinetic equations (the Boltzmann equation, the Boltzmann-Enskog equation) as well as the solutions describing the macroscopic evolution of fluids: the Euler and the Navier-Stokes hydrodynamic equations.     

Interference between a centered rarefaction wave and an inclined shock wave

The gas flow in the zone of interaction between an oblique shock and a centered isentropic rarefaction wave is studied using the direct statistical simulation method for solving the Boltzmann equation. The data of calculations of the shock and rarefaction wave structures, flow fields, and streamlines are given for the free-stream Mach number M_∞ = 6, 4 and 2. The formation of the interaction zone is simulated by a gas flow past a double-plane wedge in which the break of the generating line leads to formation of the centered isentropic rarefaction wave. The results of calculations of this flow in solving the Boltzmann equation are given in the Euler approximation.     

Invariant Manifolds for Steady Boltzmann Flows and Applications

We develop a theory of invariant manifolds for the steady Boltzmann equation and apply it to the study of boundary layers and nonlinear waves. The steady Boltzmann equation is an infinite dimensional differential equation, so the standard center manifold theory for differential equations based on spectral information does not apply here. Instead, we employ a time-asymptotic approach using the pointwise information of Green’s function for the construction of the linear invariant manifolds. At the resonance cases when the Mach number at the far field is around one of the critical values of ?1, 0 or 1, the truly nonlinear theory arises. In such a case, there are wave patterns combining the fast decaying Knudsen-type and slow varying fluid-like waves. The key Knudsen manifolds consisting of only Knudsentype layers are constructed through delicate analysis of identifying the singular behavior around the critical Mach numbers. Around Mach number ± 1, the fluidlike waves are compressive and expansive waves; and around the Mach number 0, they are linear thermal layers. The quantitative analysis of the fluid-like waves is done using the reduction of dimensions to the center manifolds.Two-scale nonlinear dynamics based on those on the Knudsen and center manifolds are formulated for the study of the global dynamics of the combined wave patterns. There are striking bifurcations in the transition of evaporation to condensation and in the transition of the Milne’s problem with a subsonic far field to one with a supersonic far field. The analysis of these wave patterns allows us to understand the Sone Diagram for the study of the complete condensation boundary value problem. The monotonicity of the Boltzmann shock profiles, a problem that initially motivated the present study, is shown as a consequence of the quantitative analysis of the nonlinear fluid-like waves.     

The Acoustic Limit for the Boltzmann Equation

The acoustic equations are the linearization of the compressible Euler equations about a spatially homogeneous fluid state. We first derive them directly from the Boltzmann equation as the formal limit of moment equations for an appropriately scaled family of Boltzmann solutions. We then establish this limit for the Boltzmann equation considered over a periodic spatial domain for bounded collision kernels. Appropriately scaled families of DiPerna-Lions renormalized solutions are shown to have fluctuations that converge entropically (and hence strongly in L ¹) to a unique limit governed by a solution of the acoustic equations for all time, provided that its initial fluctuations converge entropically to an appropriate limit associated to any given L ² initial data of the acoustic equations. The associated local conservation laws are recovered in the limit. Accepted: October 22, 1999     

Some remarks about the resolution of high velocity flows near low densities

High velocity flows which are exposed to strong rarefaction waves and creating low densities regions in it present difficulties and inaccuracies for numerical resolution. In particular, variables related to the internal energy are wrongly evaluated. Use of classical schemes solving the Euler equations in conservative variables introduces significant errors in the determination of temperature. We recommend to employ a non-conservative formulation of the energy equation. Results found to be more accurate in using the present internal energy formulation. In order to have the formulation available for both shock and strong rarefaction waves, we propose a hybrid formulation of conservative and non-conservative ones, depending on a shock indicator. The results are compared with exact solutions and show a significant improvement of the accuracy. The method is then extended to two-dimensional cases. Received 28 March 1997 / Accepted 18 June 1997     

A new flux-limiting approach–based kinetic scheme for the Euler equations of gas dynamics

This paper proposes a new kinetic-theory-based high-resolution scheme for the Euler equations of gas dynamics. The scheme uses the well-known connection that the Euler equations are suitable moments of the collisionless Boltzmann equation of kinetic theory. The collisionless Boltzmann equation is discretized using Sweby's flux-limited method and the moment of this Boltzmann level formulation gives a Euler level scheme. It is demonstrated how conventional limiters and an extremum-preserving limiter can be adapted for use in the scheme to achieve a desired effect. A simple total variation diminishing criteria relaxing parameter results in improving the resolution of the discontinuities in a significant way. A 1D scheme is formulated first and an extension to 2D on Cartesian meshes is carried out next. Accuracy analysis suggests that the scheme achieves between first- and second-order accuracy as is expected for any second-order flux-limited method. The simplicity and the explicit form of the conservative numerical fluxes add to the efficiency of the scheme. Several standard 1D and 2D test problems are solved to demonstrate the robustness and accuracy.     

旋转爆轰的三维结构和侧向稀疏波的影响

基于带化学反应的三维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对圆环形燃 烧室内的旋转爆轰进行了数值模拟,讨论了旋转爆轰波的三维结构及侧向稀疏波对旋转爆轰波阵面的影响。 数值结果表明,爆轰波能够以旋转方式沿预混气层稳定传播。在侧向稀疏波作用下,爆轰波阵面发生变形。 与理想的C-J爆轰相比,爆轰波强度和爆轰参数都有所下降。     

Interaction of Rarefaction Waves of the Two-Dimensional Self-Similar Euler Equations

We construct classical self-similar solutions to the interaction of two arbitrary planar rarefaction waves for the polytropic Euler equations in two space dimensions. The binary interaction represents a major type of interaction in the two-dimensional Riemann problems, and includes in particular the classical problem of the expansion of a wedge of gas into vacuum. Based on the hodograph transformation, the method employed here involves the phase space analysis of a second-order equation and the inversion back to (or development onto) the physical space.     

Kinetic theory based multi-level adaptive finite difference WENO schemes for compressible Euler equations

In this paper we proposed the kinetic framework based fifth-order adaptive finite difference WENO schemes abbreviated as WENO-AO-K schemes to solve the compressible Euler equations, which are quasi-linear hyperbolic equations that can admit discontinuous solutions like shock and contact waves. The formulation of the proposed schemes is based on the kinetic theory where one can recover the Euler equations by applying a suitable moment method strategy to the Boltzmann equation. The kinetic flux vector splitting strategy is used in WENO-AO framework, which produces the computationally expensive error and exponential functions. Thus, to reduce the computational cost, a physically more relevant peculiar velocity based splitting strategy is used, which is more efficient than the kinetic flux vector splitting. High order of accuracy in time is achieved using the third-order total variation diminishing Runge–Kutta (TVD-RK) scheme. Several one- and two-dimensional test cases are solved for the compressible Euler equations using the proposed fifth-order WENO-AO-K schemes and the results are compared with conventional WENO-AO scheme. Proposed schemes capture the complex flow features in a smooth region accurately, and discontinuity is also well resolved. Error analysis of the proposed schemes shows optimal convergence rates in various norms.     

Investigation of the hysteresis phenomena in steady shock reflection using kinetic and continuum methods

The problem of transition of planar shock waves over straight wedges in steady flows from regular to Mach reflection and back was numerically studied by the DSMC method for solving the Boltzmann equation and finite difference method with FCT algorithm for solving the Euler equations. It is shown that the transition from regular to Mach reflection takes place in accordance with detachment criterion while the opposite transition occurs at smaller angles. The hysteresis effect was observed at increasing and decreasing shock wave angle. Received September 1, 1995 / Accepted November 20, 1995     

On the role of turbulence in detonation induced by Mach stem reflection

A series of experiments conducted by Chan has shown that while some shock waves may not be strong enough to induce detonation when they collide with an obstacle the resulting Mach stem will induce detonation if it collides with a subsequent obstruction. A series of numerical simulations, however, failed to demonstrate the expected results if either the Euler or laminar Navier-Stokes equations are solved. On the other hand, calculations using the Favre averaged Navier-Stokes equations with a k--F turbulence model are able to reproduce the experimental results, indicating that turbulent effects may play an important role in the ignition process. A detailed examination of the results shows that turbulence causes the formation of activated kernels in a similar process to that observed in deflagration-detonation transition. The simulations in this paper have been undertaken using a modern high resolution hydrocode and a reduced kinetics mechanism for hydrogen combustion. The paper describes the reduced mechanism, the solution methods employed in the hydrocode and discusses the results of the simulations and their implications. Received 28 October 1997 / Accepted 30 April 1998     

A multi-entropy-level lattice Boltzmann model for the one-dimensional compressible Euler equations

In this paper, we propose a new lattice Boltzmann model for the one-dimensional compressible Euler equations. The new model is based on a three-entropy-level and three-speed lattice Boltzmann equation by using a method of higher-order moments of the equilibrium distribution functions. In order to obtain the second-order accuracy model, we employ the ghost field distribution functions to remove the non-physical dissipation terms in the Euler equations. We also use the conditions of the higher-order moments of the ghost field equilibrium distribution functions to obtain the equilibrium distribution functions. The numerical examples show that the numerical results can be compared with those classical methods.     

From the Boltzmann Equations¶to the Equations of¶Incompressible Fluid Mechanics, I

We consider here the problem of deriving rigorously from Boltzmann's equation, globally in time and for general initial conditions, fluid mechanics equations such as the Navier-Stokes or Euler equations.     

A lattice Boltzmann model for the compressible Euler equations with second‐order accuracy

In this paper, we propose a new lattice Boltzmann model for the compressible Euler equations. The model is based on a three‐energy‐level and three‐speed lattice Boltzmann equation by using a method of higher moments of the equilibrium distribution functions. In order to obtain second‐order accuracy, we employ the ghost field distribution functions to remove the non‐physical viscous parts. We also use the conditions of the higher moment of the ghost field equilibrium distribution functions to obtain the equilibrium distribution functions. In the numerical examples, we compare the numerical results of this scheme with those obtained by other lattice Boltzmann models for the compressible Euler equations. Copyright © 2008 John Wiley & Sons, Ltd.     

Counter-examples to Concentration-cancellation

We study the existence and the asymptotic behavior of large amplitude high-frequency oscillating waves subjected to the two-dimensional Burger equation. This program is achieved by developing tools related to supercritical WKB analysis. By selecting solutions which are divergence free, we show that incompressible or compressible two-dimensional Euler equations are not locally closed for the weak L ² topology.     

Macroscopic description of steady and unsteady rarefaction effects in boundary value problems of gas dynamics

Four basic flow configurations are employed to investigate steady and unsteady rarefaction effects in monatomic ideal gas flows. Internal and external flows in planar geometry, namely, viscous slip (Kramer’s problem), thermal creep, oscillatory Couette, and pulsating Poiseuille flows are considered. A characteristic feature of the selected problems is the formation of the Knudsen boundary layers, where non-Newtonian stress and non-Fourier heat conduction exist. The linearized Navier–Stokes–Fourier and regularized 13-moment equations are utilized to analytically represent the rarefaction effects in these boundary-value problems. It is shown that the regularized 13-moment system correctly estimates the structure of Knudsen layers, compared to the linearized Boltzmann equation data.     

Vanishing Viscosity Limit of the Compressible Navier–Stokes Equations for Solutions to a Riemann Problem

We study the vanishing viscosity limit of the compressible Navier–Stokes equations to the Riemann solution of the Euler equations that consists of the superposition of a shock wave and a rarefaction wave. In particular, it is shown that there exists a family of smooth solutions to the compressible Navier–Stokes equations that converges to the Riemann solution away from the initial and shock layers at a rate in terms of the viscosity and the heat conductivity coefficients. This gives the first mathematical justification of this limit for the Navier–Stokes equations to the Riemann solution that contains these two typical nonlinear hyperbolic waves.     

Kinetic Models for a Gas Filled Porous Matrix. WAVE PROPAGATION

The propagation of small perturbation in a gas filled porous matrix is investigated. The skeleton is supposed rigid and governed by the energy balance equation, where the heat exchanged between the two phases is taken into account. The Boltzmann equation is written for the gas where the integrals of the collisions between gas and solid particles are evaluated as those for the particles of a mixture. Different choices of the time and space scales lead to models equations which hold for different rarefaction regimes. The wave propagation characteristics are then dealt with in various situations.