Entropy considerations in numerical simulations of non-equilibrium rarefied flows

Abstract Non-equilibrium rarefied flows are encountered frequently in supersonic flight at high altitudes, vacuum technology and in microscale devices. Prediction of the onset of non-equilibrium is important for accurate numerical simulation of such flows. We formulate and apply the discrete version of Boltzmann’s H -theorem for analysis of non-equilibrium onset and accuracy of numerical modeling of rarefied gas flows. The numerical modeling approach is based on the deterministic solution of kinetic model equations. The numerical solution approach comprises the discrete velocity method in the velocity space and the finite volume method in the physical space with different numerical flux schemes: the first-order, the second-order minmod flux limiter and a third-order WENO schemes. The use of entropy considerations in rarefied flow simulations is illustrated for the normal shock, the Riemann and the two-dimensional shock tube problems. The entropy generation rate based on kinetic theory is shown to be a powerful indicator of the onset of non-equilibrium, accuracy of numerical solution as well as the compatibility of boundary conditions for both steady and unsteady problems.

[1]  Ludwig Boltzmann,et al.  Lectures on Gas Theory , 1964 .

[2]  B. Perthame,et al.  The Gaussian-BGK model of Boltzmann equation with small Prandtl number , 2000 .

[3]  C. Schrock,et al.  Continuum Onset Parameter Based on Entropy Gradients Using Boltzmann's H-Theorem , 2005 .

[4]  Mikhail S. Ivanov,et al.  Numerical simulation of shock wave propagation in microchannels using continuum and kinetic approaches , 2009 .

[5]  L. Mieussens.  Discrete-Velocity Models and Numerical Schemes for the Boltzmann-BGK Equation in Plane and Axisymmetric Geometries , 2000 .

[6]  Taku Ohwada,et al.  Structure of normal shock waves: Direct numerical analysis of the Boltzmann equation for hard-sphere molecules , 1993 .

[7]  Alina A. Alexeenko,et al.  Modeling of Viscous Shock Tube Using ES-BGK Model Kinetic Equations , 2009 .

[8]  P. Bhatnagar,et al.  A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems , 1954 .

[9]  V. V. Aristov,et al.  Unified solver for rarefied and continuum flows with adaptive mesh and algorithm refinement , 2007, J. Comput. Phys..

[10]  Lowell H. Holway,et al.  New Statistical Models for Kinetic Theory: Methods of Construction , 1966 .

[11]  Florin Iancu,et al.  Efficiency of shock wave compression in a microchannel , 2006 .

[12]  G. Bird.  Molecular Gas Dynamics and the Direct Simulation of Gas Flows , 1994 .

[13]  David B. Goldstein,et al.  Hybrid Euler/Particle Approach for Continuum/Rarefied Flows , 1998 .

[14]  Luc Mieussens,et al.  Numerical comparison of Bhatnagar–Gross–Krook models with proper Prandtl number , 2004 .

[15]  R. C. TOMTaNSON.  COMPRESSIBLE FLOW , 2002 .

[16]  G. F. Naterer,et al.  Entropy and the Second Law Fluid Flow and Heat Transfer Simulation , 2003 .

[17]  Alina A. Alexeenko,et al.  Non-Equilibrium Flow Modeling Using High-Order Schemes for the Boltzmann Model Equations , 2008 .

[18]  M. Krook.  Continuum equations in the dynamics of rarefied gases , 1959, Journal of Fluid Mechanics.

[19]  Luc Mieussens,et al.  Macroscopic Fluid Models with Localized Kinetic Upscaling Effects , 2006, Multiscale Model. Simul..

[20]  Christopher R. Schrock,et al.  Calculation of Entropy Generation Rates via DSMC with Application to Continuum/Equilibrium Onset , 2005 .

[21]  Shee‐Mang Yen,et al.  Temperature Overshoot in Shock Waves , 1966 .

[22]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[23]  Jonathan M. Burt,et al.  A hybrid particle approach for continuum and rarefied flow simulation , 2009, J. Comput. Phys..

[24]  Chi-Wang Shu,et al.  Numerical Comparison of WENO Finite Volume and Runge–Kutta Discontinuous Galerkin Methods , 2001, J. Sci. Comput..

[25]  Martin Brouillette,et al.  Shock waves at microscales , 2003 .

[26]  S. Osher,et al.  Uniformly High-Order Accurate Nonoscillatory Schemes. I , 1987 .