A lattice Boltzmann model for the generalized Burgers-Huxley equation

In this paper we develop a lattice Boltzmann model for the generalized Burgers-Huxley equation (GBHE). By choosing the proper time and space scales and applying the Chapman-Enskog expansion, the governing equation is recovered correctly from the lattice Boltzmann equation, and the local equilibrium distribution functions are obtained. Excellent agreement with the exact solution is observed, and better numerical accuracy is obtained than the available numerical result. The results indicate the present model is satisfactory and efficient. The method can also be applied to the generalized Burgers-Fisher equation and be extended to multidimensional cases.     

Lattice Boltzmann method for the generalized Kuramoto-Sivashinsky equation

In this paper, a lattice Boltzmann model with an amending function is proposed for the generalized Kuramoto-Sivashinsky equation that has the form u_t+uu_x+αu_xx+βu_xxx+γu_xxxx=0. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. It is found that the numerical results agree well with the analytical solutions.     

Hydrodynamic limits for the Boltzmann process

We study the behavior of the nonlinear Markov process associated to the Boltzmann equation under both hyperbolic and parabolic space-time scalings. In the first case the limit of the process is the solution of an o.d.e. with vector field given by a solution of the Euler equation, while in the second case the limit of the process, in the incompressible case, turns out to be a diffusion process whose drift is a solution of the incompressible Navier-Stokes equation.     

On the Convergence of the Boltzmann Equation for Semiconductors Toward the Energy Transport Model

The diffusion limit of the Boltzmann equation of semiconductors is analyzed. The dominant collisions are the elastic collisions on one hand and the electron–electron collisions with the Pauli exclusion terms on the other hand. Under a nondegeneracy hypothesis on the distribution function, a lower bound of the entropy dissipation rate of the leading term of the Boltzmann kernel for semiconductors in terms of a distance to the space of Fermi–Dirac functions is proved. This estimate and a mean compactness lemma are used to prove the convergence of the solution of the Boltzmann equation to a solution of the energy transport model.     

The hilbert expansion to the Boltzmann equation for steady flow

The Hilbert expansion to the Boltzmann equation is carried out for steady flow. It is shown that the first term in the Hubert series for the distribution function is a local Maxwellian leading to the steady Euler equations. The steady field equations that follow from the solution of the second term in the series are derived. The formulas for thermal conductivity and for viscosity of Hilbert that appear in the steady field equations of the second approximation are shown to be precisely the same as those obtained by Chapman and Enskog. The procedure to obtain higher approximations by Hubert's method is summarized.     

Derivation of Slip boundary conditions for the Navier-Stokes system from the Boltzmann equation

Rarefied gas flow behavior is usually described by the Boltzmann equation, the Navier-Stokes system being valid when the gas is less rarefied. Slip boundary conditions for the Navier-Stokes equations are derived in a rigorous and systematic way from the boundary condition at the kinetic level (Boltzmann equation). These slip conditions are explicitly written in terms of asymptotic behavior of some linear half-space problems. The validity of this analysis is established in the simple case of the Couette flow, for which it is proved that the right boundary conditions are obtained.     

Multiple-relaxation-time lattice Boltzmann model for the convection and anisotropic diffusion equation

A lattice Boltzmann model with a multiple-relaxation-time (MRT) collision operator for the convection–diffusion equation is presented. The model uses seven discrete velocities in three dimensions (D3Q7 model). The off-diagonal components of the relaxation-time matrix, which originate from the rotation of the principal axes, enable us to take into account full anisotropy of diffusion. An asymptotic analysis of the model equation with boundary rules for the Dirichlet and Neumann-type (specified flux) conditions is carried out to show that the model is first- and second-order accurate in time and space, respectively. The results of the analysis are verified by several numerical examples. It is also shown numerically that the error of the MRT model is less sensitive to the variation of the relaxation-time coefficients than that of the classical BGK model. In addition, an alternative treatment for the Neumann-type boundary condition that improves the accuracy on a curved boundary is presented along with a numerical example of a spherical boundary.     

Solution of the BGK model kinetic equation for very hard particle interaction

We study the Bhatnagar-Gross-Krook model kinetic equation with a velocity-dependent collision frequency. We derive the conditions that must be verified in order to keep the main physical properties of the Boltzmann equation, i.e.,H-theorem and conservation laws. The particular case of the so-called VHP interaction is considered, and the resulting kinetic equation is solved for a homogeneous and isotropic gas. Overpopulation phenomena are observed and analyzed for some kinds of initial conditions. The results are compared, where possible, with the exact solution of the Boltzmann equation.     

Study on the unified algorithm for three-dimensional complex problems covering various flow regimes using Boltzmann model equation

The Boltzmann simplified velocity distribution function equation describing the gas transfer phenomena from various flow regimes will be explored and solved numerically in this study. The discrete velocity ordinate method of the gas kinetic theory is studied and applied to simulate the complex multi-scale flows. Based on the uncoupling technique on molecular movement and colliding in the DSMC method, the gas-kinetic finite difference scheme is constructed to directly solve the discrete velocity distribution functions by extending and applying the unsteady time-splitting method from computational fluid dynamics. The Gauss-type discrete velocity numerical quadrature technique for different Mach number flows is developed to evaluate the macroscopic flow parameters in the physical space. As a result, the gas-kinetic numerical algorithm is established to study the three-dimensional complex flows from rarefied transition to continuum regimes. The parallel strategy adapted to the gas-kinetic numerical algorithm is investigated by analyzing the inner parallel degree of the algorithm, and then the HPF parallel processing program is developed. To test the reliability of the present gas-kinetic numerical method, the three-dimensional complex flows around sphere and spacecraft shape with various Knudsen numbers are simulated by HPF parallel computing. The computational results are found in high resolution of the flow fields and good agreement with the theoretical and experimental data. The computing practice has confirmed that the present gas-kinetic algorithm probably provides a promising approach to resolve the hypersonic aerothermodynamic problems with the complete spectrum of flow regimes from the gas-kinetic point of view of solving the Boltzmann model equation. Supported by the National Natural Science Foundation of China (Grant Nos. 90205009 and 10321002) and the National Parallel Computing Center     

A kinetic model for droplet growth in the transition regime

A variant of the moment expansion method, used in an earlier paper to describe the flow of a gas toward an absorbing sphere, is applied to a more realistic model of a droplet condensing from a supersaturated vapor. In the simplest version a spherical droplet absorbs all incoming vapor molecules, but spontaneously emits molecules with a Maxwellian distribution at the droplet temperature and with the corresponding saturated vapor density. From a solution of the stationary linearized Boltzmann equation with these boundary conditions we obtain expressions for the heat and mass currents toward the sphere as a function of the supersaturation and the temperature difference between the droplet and the vapor at infinity. For small droplet radii the known free flow limit is obtained in a natural way. From the calculated expressions for the heat and mass current we derive evolution equations for the radius and temperature of the droplet. The temperature evolves more rapidly and can thus be eliminated adiabatically; the resulting growth curve for the radius shows a sharp transition from a kinetically controlled regime for small radii to a regime dominated by heat conduction for large radii. The effect of incomplete absorption at the surface is also studied. The actual calculations are carried out for Maxwell molecules, with parameters corresponding to argon at 0.65T _c and 100% supersaturation.     

Numerical study on the gas-kinetic high-order schemes for solving Boltzmann model equation

The high-order compact finite difference technique is introduced to solve the Boltzmann model equation, and the gas-kinetic high-order schemes are developed to simulate the different kinetic model equations such as the BGK model, the Shakhov model and the Ellipsoidal Statistical (ES) model in this paper. The methods are tested for the one-dimensional unsteady shock-tube problems with various Knudsen numbers, the inner flows of normal shock wave for different Mach numbers, and the two-dimensional flows past ...     

Origin of the Landau-Lifshitz hydrodynamic fluctuations in nonequilibrium systems and a new method for reducing the Boltzmann equation

The Landau-Lifshitz fluctuating fluxes in fluctuating hydrodynamics are derived from the deterministic Boltzmann equation with the aid of a reduction method developed by Fujisaka and Mori. Thus it is shown that the hydrodynamic fluctuations innonequilibrium systems are generated by the reduction of variables from the-space distribution function to its five momentum moments, i.e., the hydrodynamic variables. This differs from the Bixon-Zwanzig and Fox-Uhlenbeck theories, in which the Landau-Lifshitz fluctuating fluxes are derived from the molecular fluctuating force in thestochastic Boltzmann-Langevin equation, which is, however, negligible in nonequilibrium systems. Thus the present method improves the Chapman-Enskog reduction method so as to include the hydrodynamic fluctuations generated by the reduction of variables.Supported in part by the Scientific Research Fund of the Ministry of Education.     

Application of microscopic transport model in the study of nuclear equation of state from heavy ion collisions at intermediate energies

The equation of state (EOS) of nuclear matter, i.e., the thermodynamic relationship between the binding energy per nucleon, temperature, density, as well as the isospin asymmetry, has been a hot topic in nuclear physics and astrophysics for a long time. The knowledge of the nuclear EOS is essential for studying the properties of nuclei, the structure of neutron stars, the dynamics of heavy ion collision (HIC), as well as neutron star mergers. HIC offers a unique way to create nuclear matter with high density and isospin asymmetry in terrestrial laboratory, but the formed dense nuclear matter exists only for a very short period, one cannot measure the nuclear EOS directly in experiments. Practically, transport models which often incorporate phenomenological potentials as an input are utilized to deduce the EOS from the comparison with the observables measured in laboratory. The ultrarelativistic quantum molecular dynamics (UrQMD) model has been widely employed for investigating HIC from the Fermi energy (40 MeV per nucleon) up to the CERN Large Hadron Collider energies (TeV). With further improvement in the nuclear mean-field potential term, the collision term, and the cluster recognition term of the UrQMD model, the newly measured collective flow and nuclear stopping data of light charged particles by the FOPI Collaboration can be reproduced. In this article we highlight our recent results on the studies of the nuclear EOS and the nuclear symmetry energy with the UrQMD model. New opportunities and challenges in the extraction of the nuclear EOS from transport models and HIC experiments are discussed.