Abstract
Mathematically, the typical difference of discrete Boltzmann model (DBM) from the traditional hydrodynamic one is that the Navier-Stokes (NS) equations are replaced by a discrete Boltzmann equation. But physically, this replacement has a significant gain: a DBM is roughly equivalent to a hydrodynamic model supplemented by a coarse-grained model of the thermodynamic non-equilibrium (TNE) effects, where the hydrodynamic model can be and can also be beyond the NS. Via the DBM, it is convenient to perform simulations on systems with flexible Knudsen number. The observations on TNE are being obtaining more applications with time.
Keywords
- Boltzmann equation
- discrete Boltzmann modeling
- Navier-Stokes
- compressible flows
- non-equilibrium effects
1. Introduction
Generally, compressible flow is frequently referred also to as gas dynamics which is the branch of fluid mechanics that deals with flows having significant changes in fluid density. That is because gases, mostly, display such behavior. In fact, all materials are compressible. Besides liquids, gases and plasmas, to some extent, the plastic solids under strong shock can also be modeled as compressible flows. Because, in the last case, that the strength of material is negligibly smaller than that of the shock. Flows with a Mach number less than 0.3 are usually treated as being incompressible for that the variation of density due to velocity is less than 5% in that case. The study on compressible flow is relevant to various fields, such as high-speed aircraft, rocket motors, jet engines, high-speed entry into a planetary atmosphere, gas pipelines, etc. where shock wave and/or detonation play a significant role. In this chapter, the following several kinds of flows including high Mach number flows with combustion, multiphase flows with phase separation and complex flows with hydrodynamic instability1 are taken as examples, and all flows are treated uniformly as being compressible. It is hopeful that the methods and ideas developed in this chapter may be adapted, more or less, to other kinds of complex flows.
Some common and typical features of these flows are as below: each of them possesses multi-scale structures and/or kinetic modes. There exist plenty of interfaces inside the system. The interfaces include material interfaces and mechanical interfaces. Each of them experiences complicated competitions between behaviors in various spatial-temporal scales. The forcing and responsive processes inside the system are very complicated. Such a system generally shows pronounced non-equilibrium behaviors.
1.1. Traditional models
The traditional macroscopic models of compressible flows are generally based on Navier-Stokes (NS) equations, even Euler equations. The model of Euler equations assumes that the system is always at its local thermodynamic equilibrium (LTE). The NS model considers the thermodynamic non-equilibrium (TNE) via the viscous stress and heat flux. The viscous stress and heat flux compose a set of convenient and effective description of the TNE. But the description is also quite dense and coarse-grained. Many specific TNE behaviors are invisible under NS description, even though they are helpful for understanding the specific TNE status. Besides, since it includes only the first-order term of Knudsen number in the Chapman-Enskog expansion [1], it is reasonable only when the Knudsen number is very small. It cannot be used to access deeper TNE behaviors. To access the complicated non-equilibrium behaviors, one possible solution is to use the molecular dynamics (MD) [2] or direct simulation Monte Carlo (DSMC) [3]. The MD simulation can help to understand some fundamental mechanisms from the atomic level. But the spatial and temporal scales it can access are too small to be comparable with experiments. The DSMC simulation has a similar constraint of computational cost.
The NS model is not enough to capture so complex non-equilibrium behaviors while the MD and DSMC cannot access spatial-temporal scales that are large enough. Under such conditions, a kinetic approach based on the non-equilibrium statistical mechanics (NESM) [1, 4, 5] is preferable.
1.2. Non-equilibrium statistical mechanics
Statistical mechanics is a branch of theoretical physics. It was developed to study the average behavior of a mechanical system with uncertain state by using probability theory. Microscopic mechanical laws do not contain concepts such as temperature, heat, or entropy; however, statistical mechanics shows how these concepts arise from the natural uncertainty about the state of a system when that system is prepared in practice. Statistical mechanics provides exact methods to connect thermodynamic quantities to microscopic behavior.
The NESM is based on mechanics and some necessary assumptions. The concept of macroscopic observation and assumption of coarse-grained density are cornerstones. The Liouville equation [4, 5] is the most fundamental governing equation when without quantum fluctuations. It describes the
where
where
Here
and
When
We need to simplify the model via considering some simpler cases. If considering only the case where correlations among three and more particles are negligible, the two-particle interaction is relevant to their distance, and the two-particle distribution function can be written as the product of two single particle distribution functions, specifically,
we obtain the Boltzmann equation [1, 4],
with
Here is the single particle distribution function. Compared with the MD and Liouville equation, the Boltzmann equation is a much coarse-grained model. The purposes of establishing Boltzmann equation are to define and calculate the entropy in non-equilibrium state, to derive, prove, even modify the fundamental hydrodynamic differential equation [5].
1.3. Non-equilibrium statistical mechanics and macroscopic description
All hydrodynamic quantities are some kinds of kinetic moments of the distribution function and can be expressed as:
where the particle mass has been assumed to be 1,
and
respectively. Compared with Boltzmann equation, in Navier-Stokes equations,
where
From molecular dynamics to Boltzmann equation, to Navier-Stokes, and further to Euler equations, in each step, the description becomes coarse-grained, and the contained physical information becomes less [6]. The switching of model in each step corresponds to the state of system under consideration gets closer to its thermodynamic equilibrium, the behavior is simpler, consequently the system can be described by fewer physical variables. For a given non-equilibrium system, the switching of model in each step corresponds to the spatial-temporal scale that we use to observe the system becomes larger, consequently more smaller structures and quicker kinetic modes are invisible. What obtained are the remaining larger structures and slower kinetic modes. Based on Boltzmann equation, the most relevant TNE effects accompanying the hydrodynamic behaviors can be studied, in addition to the general hydrodynamic behaviors described by the hydrodynamic model.
2. Discrete Boltzmann theory
2.1. Discrete Boltzmann modeling
From Boltzmann equation to DBM, two steps of coarse-grained physical modelings are needed. The principle for coarse-grained modeling is that the physical quantities used to measure the system must keep the same values after simplification.
Step 1: Linearization of the collision term
Even though, compared with MD or Liouville equation, Boltzmann equation is a much coarse-grained model, its collision term is still too complicated to be solved in most practical cases. The simplest way to simplify is to introduce a local equilibrium velocity distribution function,
The physical meaning of the linearized collision model is thus: collisions of molecules result in that
where
It should be note that the Single-Relaxation-Time model works when all the kinetic modes approaching to thermodynamic equilibrium share more or less the same relaxation time. For more complicated cases where the relaxation times of different kinetic modes approaching to thermodynamic equilibrium are significantly different, the multiple-relaxation-time (MRT) collision model is needed [15].
Step 2. Discretization of the particle velocity space
We first consider the case without the force term. The discrete Boltzmann equation reads,
where
where the left hand side gives the kinetic moments of
where the left hand side gives the kinetic moments of
where
Via the same idea, it is straight forward to formulate MRT-DBM. To ensure the relaxation times have clear physical meanings, in the MRT-DBM, the collision term is first calculated in the kinetic moment space, and then transformed back to the discrete velocity space. To ensure DBM describe reasonable flow behaviors, a correction term is needed [15].
According to the Chapman-Enskog analysis, to access system which is deeper into thermodynamic non-equilibrium, higher order terms in Knudsen number should be considered. As a result, the requirement of
2.2. Non-equilibrium: definition and measuring
Once the concept of equilibrium is defined, the concept of non-equilibrium is clear. In classical mechanics, a particle is in mechanical equilibrium if the net force on that particle is zero. By extension, a physical system made up of many parts is in mechanical equilibrium if the net force on each of its individual parts is zero. In addition to defining mechanical equilibrium in terms of force, there are many alternative definitions for mechanical equilibrium which are all mathematically equivalent. In terms of momentum, a system is in equilibrium if the momentum of its parts is all constant. In terms of velocity, the system is in equilibrium if velocity is constant. In a rotational mechanical equilibrium, the angular momentum of the object is conserved and the net torque is zero. More generally in conservative systems, equilibrium is established at a point in configuration space where the gradient of the potential energy with respect to the generalized coordinates is zero. Similarly, a fluid system is in fluid mechanical equilibrium or hydrodynamic equilibrium if the net force on each of its ‘fluid particles’ (small fluid elements) is zero and without temperature gradient around the ‘fluid particle’.
For ideal gas system, in thermodynamic equilibrium there are no net macroscopic flows of matter or energy either within a system or between systems. In non-equilibrium systems, by contrast, there are net flows of matter or energy. Global thermodynamic equilibrium means that the relevant intensive parameters are homogeneous throughout the whole system, while local thermodynamic equilibrium means that those intensive parameters are varying in space and time, but are varying so slowly that, for any point, one can assume thermodynamic equilibrium in some neighborhood about that point. Rarefied gases at ordinary temperatures behave very nearly like ideal gas and the Maxwell speed distribution is an excellent approximation for such gases. Thus, it forms the basis of the kinetic theory of gases.
It is evident that Euler equations are used to investigate fluid flows which are at local thermodynamic equilibrium but mechanical non-equilibrium, while NS, Burnett and Super-Burnett equations are used to investigate fluid flows which are at mechanical and thermodynamic non-equilibrium. Only when all kinds of kinetic moments of (
In many practical cases, it is neither necessary to know all the details of
A centrally important motivation of DBM is to check, measure and analyze the non-equilibrium state and effects [6, 16, 17, 18]. The DBM presents two sets of measures for the TNE. One set is dynamically from the difference of the kinetic moments of
where
Entropy production is a highly concerned quantity in both physics and engineering studies. From the physics side, it is helpful for understanding the complex non-equilibrium behaviors. From the engineering side, a process with lower entropy production may have higher energy transformation efficiency. Following the way of defining entropy equilibrium equation in the non-equilibrium thermodynamics, a new entropy equilibrium equation can be obtained as follows [19],
where
and
are the entropy flux and entropy production rate, respectively.
The TNE behaviors are very complex and difficult to quantitatively investigate. Finding a convenient and efficient method to characterize the TNE status and effects is the corner stone, DBM presents such a potential approach [16, 17, 18, 19, 20, 21, 22, 23, 24, 25].
2.3. DBM versus CFD
The traditional CFD needs to first know the exact form of the hydrodynamic equations, then design numerical algorithm according to the properties of those equations. The DBM is a coarse-grained model derived from the Boltzmann equation. In principle, it can be formulated and applied to simulate flows without knowing the exact form of the hydrodynamic equations, only if necessary kinetic moments of
3. Applications
3.1. Combustion system
Combustion has long been playing a dominant role in the transportation and power generation. To improve combustion efficiency and decrease pollution, in recent years, some new combustion concepts have been proposed. For example, pulsed detonation engine, spinning detonation engine, microscale combustion, nanopropellent, partially premixed and stratified combustion, plasma-assistant combustion, cool flames, etc. All these new combustion concepts involve complicated non-equilibrium chemical and transport processes [26].
The chemical reaction process is very complex and may include varieties of reaction mechanism. So far, most of the chemical reaction kinetic models are phenomenological. As the first step, we consider only the simplified form of Lee-Tarver chemical reaction rate law [21]. Considering the thermal initiation, the reaction kinetic is described by
where
where
Non-equilibrium quantities defined in Eq. (22) are used to study a simple case of detonation [21]. For the case of CJ detonation shown in Figure 1(a), the corresponding non-equilibrium quantities,
Figure 2 gives the pressure P and corresponding non-equilibrium quantity
Figure 3 shows some numerical results aiming to investigate the main mechanisms resulting in entropy increase and their relative importance in the combustion system [19]. It is clear that, in the checked cases, the most pronounced contribution to entropy increase is from the chemical reaction,
3.2. Multiphase flow with phase separation
Phase separation is an important branch in the field of multiphase flows. It is also a kind of non-equilibrium phase transition. The key step for modeling phase separation is to incorporate the non-ideal gas effects into the discrete Boltzmann equation. Enskog equation can be regarded as an extension of Boltzmann equation under the hard-ball molecule model [1]. Although the specific treatments may be different, the aims are the same. Those are all to replace the equation of state of ideal gas with a more practical one.
In 2007 Gonnella, Lamura, and Sofonea (GLS) [27] introduced an appropriate inter-particle interaction to the Watari’s model [28] to describe van der Waals fluids. The evolution equation of GLS model reads:
where the external force term
In a recent work [20], the GLS model was further developed to be a kinetic model which can be used to access both the hydrodynamic non-equilibrium and the thermodynamic non-equilibrium. To roughly and averagely estimate the derivation amplitude from the thermodynamic equilibrium, a TNE strength can be defined as
Figure 4 shows that the maximum value point can work as a physical criterion to discriminate the two stages, spinodal decomposition and domain growth, of phase separation. The TNE strength increases with time in the first stage while decreases with time in the second stage. More details are referred to the original publication [20].
3.3. Rayleigh-Taylor interfacial instability
Rayleigh-Taylor instability (RTI) occurs at the interface between two fluids with different densities. The compressible RTI system can be described by [23, 24]
where
Generally, the depth of the mixing layer is an important parameter to measure the evolution of RTI. For incompressible RTI, the measurement is readily performed by tracing the constant density. However, for the compressible case, how to measure the mixing layer remains a thorny problem. Here we present two independent interface-tracking methods. One is by tracking the mean temperature of the upper and bottom fluids while the other is by tracking the maximum values of TNE characteristic quantities, such as
3.4. Compressible flows under shock
Figure 6 shows the profiles of physical quantities in the process where the shock wave passes outwards from the heavy medium to the light one. Figure 6(a) and (b) are for the case without and with initial perturbations at the interface, respectively. From left to right, one can find three kind of interfaces, the rarefaction wave, material interface and shock wave, which are indicated by dashed lines [25].
According to the TNE information, the main feature of actual particle velocity distribution function can be qualitatively recovered. Figure 7 shows an example, where the interface is not perturbed initially. The details are referred to Ref. [25]. DBM simulations [25] show that the shear stress exists only for the oblique shocking. As a consistent correspondence, MD results [29] show that fluctuating shear stresses exist if observed in a scale with a few angstroms, while their mean value becomes negligibly small when being averaged over a scale with several hundred angstroms.
A comparison of the DBM results and those of the MD is shown in Figure 8. Figure 8(a) and (b) shows the DBM results for the cases where the interface is not and is perturbed initially, respectively. Figure 8(c) shows the shear stresses from MD simulation. Since the MD uses particle description, the existence of locally fluctuating shear stress corresponds to the observation in the case with perturbed interface; the observation that mean value of shear stresses becomes negligibly small in a much larger scale roughly correspond to the case with non-perturbed interface for the DBM simulations.
3.5. Shock wave in plasma
Figure 9 shows an example for that the TNE effects can be used to physically discriminate shock wave in plasma from those in common fluid. From the first two rows, the two TNE quantities,
4. Summary
Understanding compressible flows need more time. DBM presents a convenient way to model and simulate systems with trans-scale Knudsen number. Mathematically, the only difference of discrete Boltzmann from the traditional hydrodynamic modeling is that the NS equations are replaced by a discrete Boltzmann equation. But physically, this replacement has a significant gain: a DBM is roughly equivalent to a hydrodynamic model supplemented by a coarse-grained model of the TNE, where the hydrodynamic model can be and can also beyond the NS. The TNE provided by DBM has been used to investigate non-equilibrium effect during detonation process, to discriminate different stages of phase separation, to recover actual particle velocity distribution function qualitatively, to track the interfaces of different fluid, and to discriminate shock wave in plasma from those in common fluid. More use of those TNE quantities are further being discovered with the deeper investigation of the compressible and complex flows.
Acknowledgments
We warmly thank Profs. Wei Kang, Zhihua Chen, Yanbiao Gan, Feng Chen, Chuandong Lin, Huilin Lai, Zhipeng Liu, Bo Yan, et al. for helpful discussions. We acknowledge support of the National Natural Science Foundation of China (under Grant Nos. 11475028 and 11772064), Science Challenge Project (under Grant No. JCKY2016212A501 and TZ2016002), and Science Foundation of Laboratory of Computational Physics.
References
- 1.
Chapman S, Cowling TG. The Mathematical Theory of Non-uniform Gases. third ed. New York: Cambridge: University Press; 1970 - 2.
Rapaport DC. The Art of Molecular Dynamics Simulation. 2nd ed. New York: Cambridge: University Press; 2004 - 3.
Bird GA. Molecular Gas Dynamics and the Direct Simulation of Gas Flows. New York: Clarendon Press; 2003 - 4.
Chen S. Non-equilibrium Statistical Mechanics. Beijing: Science Press; 2010 - 5.
Shen H. Statistical Mechanics. Hefei: China University of Sci. & Tech. Press; 2011 - 6.
Xu A, Zhang G, Ying Y. Progress of discrete Boltzmann modeling and simulation of combustion system. Acta Physica Sinica. 2015; 64 (18):184701 - 7.
Bhatnagar PL, Gross EP, Krook M. A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component systems. Physical Review. 1954; 94 (3):511 - 8.
Succi S. The Lattice Boltzmann Equation for Fluid Dynamics and beyond. New York: Clarendon Press; 2001 - 9.
Qian YH, d'Humières D, Lallemand P. Lattice BGK models for Navier-stokes equation. EPL (Europhysics Letters). 1992; 17 (6):479 - 10.
Chen H, Chen S, Matthaeus WH. Recovery of the Navier-stokes equations using a lattice-gas Boltzmann method. Physical Review A. 1992; 45 (8):R5339 - 11.
Holway LH Jr. New statistical models for kinetic theory: Methods of construction. The. Physics of Fluids. 1966; 9 (9):1658-1673 - 12.
Shakhov EM. Generalization of the Krook kinetic relaxation equation. Fluid. Dynamics (Pembroke, Ont.). 1968; 3 (5):95-96 - 13.
Larina IN, Rykov VA. Kinetic model of the Boltzmann equation for a diatomic gas with rotational degrees of freedom. Computational Mathematics & Mathematical Physics. 2010; 50 (12):2118-2130 - 14.
Liu G. A method for constructing a model form for the Boltzmann equation. Physics of Fluids. 1990; 2 (2):277-280 - 15.
Chen F, Xu A, Zhang G, et al. Two-dimensional MRT LB model for compressible and incompressible flows. Frontiers of Physics. 2014; 9 (2):246-254 - 16.
Xu A, Zhang G, Gan Y, et al. Lattice Boltzmann modeling and simulation of compressible flows. Frontiers of Physics. 2012; 7 (5):582-600 - 17.
Xu A, Zhang G, Li Y, et al. Modeling and simulation of nonequilibrium and multiphase complex systems: Lattice Boltzmann kinetic theory and application. Progress in Physics. 2014; 34 (3):136-167 - 18.
Xu A, Zhang G, Gan Y. Progress in studies on discrete Boltzmann modeling of phase separationprocess. Mechanics in Engineering. 2016; 38 (4):361-374 - 19.
Zhang Y, Xu A, Zhang G, et al. Kinetic modeling of detonation and effects of negative temperature coefficient. Combustion & Flame. 2016; 173 :483-492 - 20.
Gan Y, Xu A, Zhang G, et al. Discrete Boltzmann modeling of multiphase flows: Hydrodynamic and thermodynamic non-equilibrium effects. Soft Matter. 2015; 11 (26):5336 - 21.
Yan B, Xu A, Zhang G, et al. Lattice Boltzmann model for combustion and detonation. Frontiers of Physics. 2013; 8 (1):94-110 - 22.
Xu A, Lin C, Zhang G, et al. Multiple-relaxation-time lattice Boltzmann kinetic model for combustion. Physical Review E Statistical Nonlinear & Soft Matter Physics. 2015; 91 (4):043306 - 23.
Lai H, Xu A, Zhang G, et al. Nonequilibriumthermohydrodynamic effects on the Rayleigh-Taylor instability in compressible flows. Physical Review E. 2016; 94 (1):023106 - 24.
Chen F, Xu A, Zhang G. Viscosity, heat conductivity, and Prandtl number effects in the Rayleigh–Taylor instability. Frontiers of Physics. 2016; 11 (6):114703 - 25.
Lin C, Xu A, Zhang G, et al. Polar-coordinate lattice Boltzmann modeling of compressible flows. Physical Review E Statistical Nonlinear & Soft Matter Physics. 2014; 89 (1):013307 - 26.
Ju Y. Recent progress and challenges in fundamental combustion research. Advances in Mechanics. 2014; 44 (1):26-97 - 27.
Gonnella G, Lamura A, Sofonea V. Lattice Boltzmann simulation of thermal nonideal fluids. Physical Review E. 2007; 76 (3):036703 - 28.
Watari M, Tsutahara M. Two-dimensional thermal model of the finite-difference lattice Boltzmann method with high spatial isotropy. Physical Review E Statistical Nonlinear & Soft Matter Physics. 2003; 67 (2):036306 - 29.
Liu H. Shock waves and hydrodynamic instabilities under the condition of inertial confinement fusion. [Ph.D thesis]. Beijing: Peking University; 2017 - 30.
Liu Z. Non-equilibrium Characteristic of Shock Waves in Plasma, Post-Doctoral Research Report. Beijing: Institute of Applied Physics and Computational Mathematics; 2017
Notes
- It is clear that there exist significant overlaps among those classifications.