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.
- Boltzmann equation
- discrete Boltzmann modeling
- compressible flows
- non-equilibrium effects
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 , 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)  or direct simulation Monte Carlo (DSMC) . 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
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,
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 .
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,
respectively. Compared with Boltzmann equation, in Navier-Stokes equations,
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 . 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
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 .
Step 2. Discretization of the particle velocity space
We first consider the case without the force term. The discrete Boltzmann equation reads,
where the left hand side gives the kinetic moments of
where the left hand side gives the kinetic moments of
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 .
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
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 ,
are the entropy flux and entropy production rate, respectively. is a rate function describing the process of chemical reaction occurring in the system, where
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.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 .
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 . Considering the thermal initiation, the reaction kinetic is described by
where is the local equilibrium distribution function taking into account chemical reaction effect.
Non-equilibrium quantities defined in Eq. (22) are used to study a simple case of detonation . For the case of CJ detonation shown in Figure 1(a), the corresponding non-equilibrium quantities, and , are shown in Figure 1(b) and (c), respectively. Interestingly, at the position of von Neumann pressure peak, the system is not far from but is nearly in its thermodynamic equilibrium. The internal energies in two degrees of freedom deviate from their mean value in opposite direction with the same amplitude. The internal energy in each degree of freedom deviates from the mean value in opposite direction before and after the von Neumann pressure peak. The amplitude in front of the von Neumann pressure peak is larger. The physical reasons are as follows. The mechanical non-equilibrium is the driving force for TNE. When a strong shock comes, the density, temperature and flow velocity increases abruptly so that the system does not have enough time to relax to its thermodynamic equilibrium. With decrease of the gradients of density, temperature and flow velocity, the system gets relatively more time for thermodynamic relaxation. At the position of von Neumann pressure peak, the system has been close to its thermodynamic equilibrium. After the von Neumann pressure peak, the density and flow velocity decrease quickly so that the system does not have enough time for thermodynamic relaxation again. The total deriving force for TNE in front of the von Neumann pressure peak is larger than behind.
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 . 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 . 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)  introduced an appropriate inter-particle interaction to the Watari’s model  to describe van der Waals fluids. The evolution equation of GLS model reads:
where the external force term
In a recent work , 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 .
3.3. Rayleigh-Taylor interfacial instability
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 . The latter method is based on the fact that takes its maximum value at the position of the interface along the
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 .
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. . DBM simulations  show that the shear stress exists only for the oblique shocking. As a consistent correspondence, MD results  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, and , show quite similar behaviors around shock wave and/or detonation wave in common fluid, even though they two have different physical meanings. However, the two quantities show qualitative difference around shock wave in plasma .
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.
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.
- It is clear that there exist significant overlaps among those classifications.