Bulk equations and Knudsen layers for the regularized 13 moment equations

The order of magnitude method offers an alternative to the Chapman-Enskog and Grad methods to derive macroscopic transport equations for rarefied gas flows. This method yields the regularized 13 moment equations (R13) and a generalization of Grad’s 13 moment equations for non-Maxwellian molecules. Both sets of equations are presented and discussed. Solutions of these systems of equations are considered for steady state Couette flow. The order of magnitude method is used to further reduce the generalized Grad equations to the non-linear bulk equations, which are of second order in the Knudsen number. Knudsen layers result from the linearized R13 equations, which are of the third order. Superpositions of bulk solutions and Knudsen layers show good agreement with DSMC calculations for Knudsen numbers up to 0.5.      

Boltzmann equation and moment equations in curvilinear coordinates

Various forms of writing the Boltzmann equation in an arbitrary orthogonal curvilinear coordinate system are discussed. The derivation is presented of a general transport equation and moment equations containing moments of the distribution function no higher than the fourth. For a gas of Maxwellian molecules it is shown that the system of moment equations for flows which differ little from equilibrium flows transforms into the system of hydrodynamic equations. The resulting equations may be useful in solving problems on motions of a rarefied gas by the moment methods. The results are valid for both the Boltzmann equation and model kinetic equations.The author wishes to thank A. A. Nikol'skii for discussions and helpful comments.     

Asymptotics, structure, and integration of sound-proof atmospheric flow equations

Relative to the full compressible flow equations, sound-proof models filter acoustic waves while maintaining advection and internal waves. Two well-known sound-proof models, an anelastic model by Bannon and Durran’s pseudo-incompressible model, are shown here to be structurally very close to the full compressible flow equations. Essentially, the anelastic model is obtained by suppressing ? _t ρ in the mass continuity equation and slightly modifying the gravity term, whereas the pseudo-incompressible model results from dropping ? _t p from the pressure equation. For length scales small compared to the density and pressure scale heights, the anelastic model reduces to the Boussinesq approximation, while the pseudo-incompressible model approaches the zero Mach number, variable density flow equations. Thus, for small scales, both models are asymptotically consistent with the full compressible flow equations, yet the pseudo-incompressible model is more general in that it remains valid in the presence of large density variations. For the relatively small density variations found in typical atmosphere–ocean flows, both models are found to yield very similar results, with deviations between models much smaller than deviations obtained when using different numerical schemes for the same model. This in agreement with Smolarkiewicz and Dörnbrack (Int J Numer Meth Fluids 56:1513–1519, 2007). Despite these useful properties, neither model can be derived by a low-Mach number asymptotic expansion for length scales comparable to the pressure scale height, i.e., for the regime they were originally designed for. Derivations of these models via scale analysis ignore an asymptotic time scale separation between advection and internal waves. In fact, only the classical Ogura and Phillips model, which assumes weak stratification of the order of the Mach number squared, can be obtained as a leading-order model from systematic low Mach number asymptotic analysis. Issues of formal asymptotics notwithstanding, the close structural similarity of the anelastic and pseudo-incompressible models to the full compressible flow equations makes them useful limit systems in building computational models for atmospheric flows. In the second part of the paper, we propose a second-order finite-volume projection method for the anelastic and pseudo-incompressible models that observes these structural similarities. The method is applied to test problems involving free convection in a neutral atmosphere, the breaking of orographic waves at high altitudes, and the descent of a cold air bubble in the small-scale limit. The scheme is meant to serve as a starting point for the development of a robust compressible atmospheric flow solver in future work.     

Multi-shock structure of roll wavesStructure multi-chocs d'ondes à rouleaux

The special class of periodic travelling waves which is known as roll waves is investigated for nonhomogeneous hyperbolic equations of gas dynamics type. In this Note these equations are applied to shallow water flows in inclined open channels, but the results obtained are more general and far-reaching. The necessary conditions for the existence of a roll wave are derived. It is shown that for a nonconvex pressure term, multi-shock configurations of roll waves of finite amplitude exist. A new type of periodic travelling wave, which corresponds to the slug flow regime in two-layer flows, is found. To cite this article: A. Boudlal, V.Yu. Liapidevskii, C. R. Mecanique 332 (2004).     

Modifications of Gasdynamic Equations of Higher Approximations of the Chapman-Enskog Method

The non-Navier-Stokes continuum models proposed earlier on the basis of a modification of the gasdynamic equations of the higher (starting from the Burnett) approximations of the Chapman-Enskog method for shock wave flow are generalized to include the case of three-dimensional flows of a simple (monatomic) gas. The models are tested on the problems of shock wave structure and cylindrical Couette rarefied gas flow.     

Rarefied hypersonic flow in a plane channel with a surface glow discharge

The gasdynamic structure of a hypersonic molecular nitrogen flow in a plane channel whose opposite surfaces are segmented electrodes for generating a continuous surface glow discharge is investigated using a two-dimensional computational model. The electrodynamic structure of the surface glow discharge in the hypersonic rarefied gas flow (distributions of the charged particle concentrations, current density, and electric potential) is studied. A two-dimensional conjugate electrical-gasdynamic model consisting of the continuity, Navier-Stokes, and energy conservation equations and the chargedparticle continuity equations in the ambipolar approximation is proposed. The real thermophysical and transport properties of molecular nitrogen are taken into account. It is shown that using a surface glow discharge in a hypersonic rarefied gas flow makes it possible effectively to modify the shock-wave flow structure and hence to consider this type of discharge as additional tool for controlling rarefied gas flows.     

Exact solutions of the one-dimensional Russo-Smereka kinetic equation

We obtain new classes of invariant solutions of the integrodifferential equations describing the propagation of nonlinear concentration waves in a rarefied bubbly fluid. For all the solutions obtained, trajectories of particle motion in phase space are calculated. The stability of some flows is studied in a linear approximation. In several cases, the construction of solutions reduces to an integrodifferential equation of the second kind, which can be solved by the iteration method. Novosibirsk State University, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 4, pp. 21–32, July–August, 2000.     

Gas-dynamic analogy in the theory of stratified liquid flows with a free boundary

Novel approximate mathematical models of long-wave theory describing flows of a density-stratified liquid with a free boundary are proposed. It is shown that in certain cases the equations of the novel models coincide with either the equations of nonisentropic gas dynamics with a polytropic equation of state for γ = 2 or the equations describing the dynamics of a mixture of two perfect gases.     

Multicomponent fluid flow analysis using a new set of conservation equations

In this work hydrodynamics of multicomponent ideal gas mixtures have been studied. Starting from the kinetic equations, the Eulerian approach is used to derive a new set of conservation equations for the multicomponent system where each component may have different velocity and kinetic temperature. The equations are based on the Grad's method of moment derived from the kinetic model in a relaxation time approximation (RTA). Based on this model which contains separate equation sets for each component of the system, a computer code has been developed for numerical computation of compressible flows of binary gas mixture in generalized curvilinear boundary conforming coordinates. Since these equations are similar to the Navier-Stokes equations for the single fluid systems, the same numerical methods are applied to these new equations. The Roe's numerical scheme is used to discretize the convective terms of governing fluid flow equations. The prepared algorithm and the computer code are capable of computing and presenting flow fields of each component of the system separately as well as the average flow field of the multicomponent gas system as a whole. Comparison of the present code results with those of a more common algorithm based on the mixture theory in a supersonic converging-diverging nozzle provides the validation of the present formulation. Afterwards, a more involved nozzle cooling problem with a binary ideal gas (helium-xenon) is chosen to compare the present results with those of the ordinary mixture theory. The present model provides the details of the flow fields of each component separately which is not available otherwise. It is also shown that the separate fluids treatment, such as the present study, is crucial when considering time scales on the order of (or shorter than) the intercollisions relaxation times.     

Dissipative Boussinesq equations

The classical theory of water waves is based on the theory of inviscid flows. However it is important to include viscous effects in some applications. Two models are proposed to add dissipative effects in the context of the Boussinesq equations, which include the effects of weak dispersion and nonlinearity in a shallow water framework. The dissipative Boussinesq equations are then integrated numerically. To cite this article: D. Dutykh, F. Dias, C. R. Mecanique 335 (2007).     

Navier–Stokes Computations in Micro Shock Tubes

Micro shock tube flows were simulated using unsteady 2D Navier–Stokes equations combined with boundary slip velocities and temperature jumps conditions. These simulations were performed using the parallel version of a multi-block finite-volume home code. Different initial low pressures and shock tube diameters allow to have the scaling ratio ReD/4L vary. The numerical results show a strong attenuation of the shock wave strength with a decrease of the hot flow values along the tube. When the scaling ratio decreases the shock waves can transform into compression waves. Comparison to the existing 1D models also shows the limit of these models.     

A phenomenological and extended continuum approach for modelling non-equilibrium flows

This paper presents a new technique that combines Grad’s 13-moment equations (G13) with a phenomenological approach to rarefied gas flows. This combination and the proposed solution technique capture some important non-equilibrium phenomena that appear in the early continuum-transition flow regime. In contrast to the fully coupled 13-moment equation set, a significant advantage of the present solution technique is that it does not require extra boundary conditions explicitly; Grad’s equations for viscous stress and heat flux are used as constitutive relations for the conservation equations instead of being solved as equations of transport. The relative computational cost of this novel technique is low in comparison to other methods, such as fully coupled solutions involving many moments or discrete methods. In this study, the proposed numerical procedure is tested on a planar Couette flow case, and the results are compared to predictions obtained from the direct simulation Monte Carlo method. This test case highlights the presence of normal viscous stresses and tangential heat fluxes that arise from non-equilibrium phenomena, which cannot be captured by the Navier–Stokes–Fourier constitutive equations or phenomenological modifications.      

Long-time Stability in Systems of Conservation Laws,Using Relative Entropy/Energy

We study the long-time stability of shock-free solutions of hyperbolic systems of conservation laws, under an arbitrarily large initial disturbance in L ²∩ L ^∞. We use the relative entropy method, a robust tool which allows us to consider rough and large disturbances. We display practical examples in several space dimensions, for scalar equations as well as isentropic gas dynamics. For full gas dynamics, we use a trick from Chen [1], in which the estimate is made in terms of the relative mechanical energy instead of the relative mathematical entropy.     

Global Uniqueness of Transonic Shocks in Two-Dimensional Steady Compressible Euler Flows

We prove that for the two-dimensional steady complete compressible Euler system, with given uniform upcoming supersonic flows, the following three fundamental flow patterns (special solutions) in gas dynamics involving transonic shocks are all unique in the class of piecewise C ¹ smooth functions, under appropriate conditions on the downstream subsonic flows: (i) the normal transonic shocks in a straight duct with finite or infinite length, after fixing a point the shock-front passing through; (ii) the oblique transonic shocks attached to an infinite wedge; (iii) a flat Mach configuration containing one supersonic shock, two transonic shocks, and a contact discontinuity, after fixing a point where the four discontinuities intersect. These special solutions are constructed traditionally under the assumption that they are piecewise constant, and they have played important roles in the studies of mathematical gas dynamics. Our results show that the assumption of a piecewise constant can be replaced by some weaker assumptions on the downstream subsonic flows, which are sufficient to uniquely determine these special solutions. Mathematically, these are uniqueness results on solutions of free boundary problems of a quasi-linear system of elliptic-hyperbolic composite-mixed type in bounded or unbounded planar domains, without any assumptions on smallness. The proof relies on an elliptic system of pressure p and the tangent of the flow angle w = v/u obtained by decomposition of the Euler system in Lagrangian coordinates, and a newly developed method for the L ^∞ estimate that is independent of the free boundaries, by combining the maximum principles of elliptic equations, and careful analysis of the shock polar applied on the (maybe curved) shock-fronts.     

Coupled nonlinear constitutive models for rarefied and microscale gas flows: subtle interplay of kinematics and dissipation effects

The constitutive relations of gases in a thermal nonequilibrium (rarefied and microscale) can be derived by applying the moment method to the Boltzmann equation. In this work, a model constitutive relation determined on the basis of the moment method is developed and applied to some challenging problems in which classical hydrodynamic theories including the Navier–Stokes–Fourier theory are shown to predict qualitatively wrong results. Analysis of coupled nonlinear constitutive models enables the fundamentals of gas flows in thermal nonequilibrium to be identified: namely, nonlinear, asymmetric, and coupled relations between stresses and the shear rate; and effect of the bulk viscosity. In addition, the new theory explains the central minimum of the temperature profile in a force-driven Poiseuille gas flow, which is a well-known problem that renders the classical hydrodynamic theory a global failure.     

Thrust enhancement due to flexible trailing-edge of plunging foils

Drag reduction for hydrofoils is studied through thrust generation on foils plunging at low Strouhal numbers in order to simulate the action of the ocean waves. Force, deformation and flow field measurements are presented for a partially flexible plunging foil in water tunnel experiments. The foil is predominantly rigid with a short flexible trailing-edge plate of length: L=0.1c, 0.2c, or 0.3c. Using flexible plates, whose natural structural frequency is much higher than the frequency of the plunge oscillations, increases thrust compared to the rigid case. Flexibility is generally more effective for larger lengths of the flexible plate and smaller plunge amplitudes. The maximum observed is therefore for the largest length and smallest amplitude studied: L=0.3c and a=0.1c and equates to 28% more thrust than the rigid case. Optima are observed in the non-dimensional rigidity (λ) versus flap angle amplitude (δ, which is a measure of the relative deformation) parameter space. These occur at λ≈2 and δ≈7–13° for a wide range of flexible plate length and plunge amplitude. Whilst a satisfactory explanation of why there is an optimal flap amplitude remains unavailable, the case of optimal flap angle amplitude results in increased trailing-edge vortex circulation, giving a stronger reverse Kármán vortex street and thus a stronger time-averaged jet.     

Effect of dissipation on the interaction between long and short waves in bubbly liquids

The interaction of long and short waves in a rarefied monodisperse mixture of a weakly compressible viscous liquid and gas bubbles is considered. For taking the dissipation effects into account an effective-viscosity scheme is used. Four cases of dissipation are distinguished, namely, strong, medium, weak, and very weak dissipation. In the cases of moderate, weak, and very weak dissipation the equations of resonance and non-resonance wave interaction are derived using the multiscale method. The effect of "degeneration" of the interaction is detected in certain of the models constructed. In the case of "degeneration" a class of new models of the dissipative resonance interaction is constructed and investigated numerically. Ufa. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 126–138, July–August, 2000.     

Slow gas motions and the negative drag of a strongly heated spherical particle

In the slow flows of a strongly and nonuniformly heated gas, in the continuum regime (Kn → 0) thermal stresses may be present. The theory of slow nonisothermal continuum gas flows with account for thermal stresses was developed in 1969–1974. The action of the thermal stresses on the gas results in certain paradoxical effects, including the reversal of the direction of the force exerted on a spherical particle in Stokes flow. The propulsion force effect is manifested at large but finite temperature differences between the particle and the gas. This study is devoted to the thermal-stress effect on the drag of a strongly heated spherical particle traveling slowly in a gas for small Knudsen numbers (M ～ Kn → 0), small but finite Reynolds numbers (Re ≤ 1), a linear temperature dependence of the transport coefficients µ ∝ T, and large but finite temperature differences ((T _w ? T _∞ )/T _M8 ～ 1). Two different systems of equations are solved numerically: the simplified Navier-Stokes equations and the modified Navier-Stokes equations with account for the thermal stresses.     

Existence of homoenergetic affine flows for the Boltzmann equation

An existence theorem is proved for homoenergetic affine flows described by the Boltzmann equation. The result complements the analysis of Truesdell and of Galkin on the moment equations for a gas of Maxwellian molecules. Existence of the distribution function is established here for a large class of molecular models (hard sphere and angular cut-off interactions). Some of the data lead to an implosion and infinite density in a finite time, in agreement with the physical picture of the associated flows; for the remaining set of data, global existence is shown to hold.