Discrete Boltzmann trans-scale modeling of high-speed compressible flows

Lattice Boltzmann Simulation of Spatial Fractional Convection-Diffusion Equation

The space fractional advection-diffusion equation is a crucial type of fractional partial differential equation, widely used for its ability to more accurately describe natural phenomena. Due to the complexity of analytical approaches, this paper focuses on its numerical investigation. A lattice Boltzmann model for the spatial fractional convection-diffusion equation is developed, and an error analysis is carried out. The spatial fractional convection-diffusion equation is solved for several examples. The validity of the model is confirmed by comparing its numerical solutions with those obtained from other methods The results demonstrate that the lattice Boltzmann method is an effective tool for solving the space fractional convection-diffusion equation.

Effects of Nonextensive Electrons on Dust-Ion Acoustic Waves in a Collisional Dusty Plasma with Negative Ions

The effects of nonextensive electrons on nonlinear ion acoustic waves in dusty negative ion plasmas with ion-dust collisions are investigated. Analytical results show that both solitary and shock waves are supported in this system. The wave propagation is governed by a Korteweg-de Vries Burgers-type equation. The coefficients of this equation are modified by the nonextensive parameter q. Numerical calculations indicate that the amplitude of solitary wave and oscillatory shock can be obviously modified by the nonextensive electrons, but the monotonic shock is little affected.

Thermodynamic nonequilibrium effects in bubble coalescence: A discrete Boltzmann study

The thermodynamic nonequilibrium (TNE) effects in a coalescence process of two initially static bubbles under thermal conditions are investigated by a discrete Boltzmann model. The spatial distributions of the typical nonequilibrium quantity, i.e., nonorganized momentum fluxes (NOMFs), during evolutions are investigated in detail. The density-weighted statistical method is used to highlight the relationship between the TNE effects and the morphological and kinetics characteristics of bubble coalescence. The results show that the xx component and yy component of NOMFs are antisymmetrical; the xy component changes from an antisymmetric internal and external double quadrupole structure to an outer octupole structure during the coalescence process. Moreover, the evolution of the averaged xx component of NOMFs provides two characteristic instants, which divide the nonequilibrium process into three stages. The first instant, when the averaged xx component of the NOMFs reaches its first local minimum, corresponds to the moment when the mean coalescence speed gets the maximum, and at this time the ratio of minor and major axes is about 1/2. The second instant, when the averaged xx component of the NOMFs gets its second local maximum, corresponds to the moment when the ratio of minor and major axes becomes 1 for the first time. It is interesting to find that the three quantities, TNE intensity, acceleration of coalescence, and the slope of boundary length, show a high degree of correlation and attain their maxima simultaneously. The surface tension and the heat conduction accelerate the process of bubble coalescence, while the viscosity delays it. Both the surface tension and the viscosity enhance the global nonequilibrium intensity, whereas the heat conduction restrains it. These TNE features and findings present some insights into the kinetics of bubble coalescence.

Particles on demand for flows with strong discontinuities

Particles-on-demand formulation of kinetic theory [B. Dorschner, F. Bösch and I. V. Karlin, Phys. Rev. Lett. 121, 130602 (2018)0031-900710.1103/PhysRevLett.121.130602] is used to simulate a variety of compressible flows with strong discontinuities in density, pressure, and velocity. Two modifications are applied to the original formulation of the particles-on-demand method. First, a regularization by Grad's projection of particles populations is combined with the reference frame transformations in order to enhance stability and accuracy. Second, a finite-volume scheme is implemented which allows tight control of mass, momentum, and energy conservation. The proposed model is validated with an array of challenging one- and two-dimensional benchmarks of compressible flows, including hypersonic and near-vacuum situations, Richtmyer-Meshkov instability, double Mach reflection, and astrophysical jet. Excellent performance of the modified particles-on-demand method is demonstrated beyond the limitations of other lattice Boltzmann-like approaches to compressible flows.

Mesoscopic kinetic approach for studying nonequilibrium hydrodynamic and thermodynamic effects of shock wave, contact discontinuity, and rarefaction wave in the unsteady shock tube

This paper presents a detailed description of a molecular velocity distribution-based mesoscopic kinetic approach that enables a better understanding of various nonequilibrium hydrodynamic and thermodynamic effects in shock waves, contact discontinuities, and rarefaction waves. This builds on the mesoscopic kinetic approach in a previous investigation into regular reflection shocks by further addressing the mesoscopic physical meaning of kinetic moments from the view of kinetics and the implications of the magnitude and sign of nonequilibrium kinetic moments. To deepen understanding of nonequilibrium effects, this work focuses on the one-dimensional unsteady shock tube problem, which contains the typical and essential features of the discontinuous flows, and has no interference of two-dimensional flow direction. The approach uses a lattice Boltzmann method to solve the flow field, and describes nonequilibrium effects through the nonequilibrium kinetic moments of molecular velocity distribution functions. The mechanism of nonequilibrium effect in discontinuous flows is further probed. This work develops the mesoscopic kinetic approach and clarifies the mesoscopic physics of shock waves, contact discontinuities, and rarefaction waves.

Parallel finite-volume discrete Boltzmann method for inviscid compressible flows on unstructured grids

In this paper, a finite-volume discrete Boltzmann method based on a cell-centered scheme for inviscid compressible flows on unstructured grids is presented. In the new method, the equilibrium distribution functions are obtained from the circle function in two-dimensions (2D) and the spherical function in three-dimensions (3D). Moreover, the advective fluxes are evaluated by Roe's flux-difference splitting scheme, the gradients of the density and total energy distribution functions are computed with a least-squares method, and the Venkatakrishnan limiter is employed to prevent oscillations. To parallelize the method we use a graph-based partitioning approach that also guarantees the load balancing. The method is validated by seven benchmark problems: (a) a 2D flow pasting a bump, (b) a 2D Riemann problem, (c) a 2D flow passing the RAE2822 airfoil, (d) flows passing the NACA0012 airfoil, (e) 2D supersonic flows around a cylinder, (f) an explosion in a 3D box, and (g) a 3D flow around the ONERA M6 wing. The benchmark tests show that the results obtained by the proposed method match well with the published results, and the parallel numerical experiments show that the proposed parallel implementation has close to linear strong scalability, and parallel efficiencies of 95.31% and 94.56% are achieved for 2D and 3D problems on a supercomputer with up to 4800 processor cores, respectively.

Multiple-relaxation-time discrete Boltzmann modeling of multicomponent mixture with nonequilibrium effects

A multiple-relaxation-time discrete Boltzmann model (DBM) is proposed for multicomponent mixtures, where compressible, hydrodynamic, and thermodynamic nonequilibrium effects are taken into account. It allows the specific heat ratio and the Prandtl number to be adjustable, and is suitable for both low and high speed fluid flows. From the physical side, besides being consistent with the multicomponent Navier-Stokes equations, Fick's law, and Stefan-Maxwell diffusion equation in the hydrodynamic limit, the DBM provides more kinetic information about the nonequilibrium effects. The physical capability of DBM to describe the nonequilibrium flows, beyond the Navier-Stokes representation, enables the study of the entropy production mechanism in complex flows, especially in multicomponent mixtures. Moreover, the current kinetic model is employed to investigate nonequilibrium behaviors of the compressible Kelvin-Helmholtz instability (KHI). The entropy of mixing, the mixing area, the mixing width, the kinetic and internal energies, and the maximum and minimum temperatures are investigated during the dynamic KHI process. It is found that the mixing degree and fluid flow are similar in the KHI process for cases with various thermal conductivity and initial temperature configurations, while the maximum and minimum temperatures show different trends in cases with or without initial temperature gradients. Physically, both heat conduction and temperature exert slight influences on the formation and evolution of the KHI morphological structure.

Hydrodynamic and Thermodynamic Nonequilibrium Effects around Shock Waves: Based on a Discrete Boltzmann Method

A shock wave that is characterized by sharp physical gradients always draws the medium out of equilibrium. In this work, both hydrodynamic and thermodynamic nonequilibrium effects around the shock wave are investigated using a discrete Boltzmann model. Via Chapman–Enskog analysis, the local equilibrium and nonequilibrium velocity distribution functions in one-, two-, and three-dimensional velocity space are recovered across the shock wave. Besides, the absolute and relative deviation degrees are defined in order to describe the departure of the fluid system from the equilibrium state. The local and global nonequilibrium effects, nonorganized energy, and nonorganized energy flux are also investigated. Moreover, the impacts of the relaxation frequency, Mach number, thermal conductivity, viscosity, and the specific heat ratio on the nonequilibrium behaviours around shock waves are studied. This work is helpful for a deeper understanding of the fine structures of shock wave and nonequilibrium statistical mechanics.

Dynamics of bond breaking and formation in polyethylene near shock front

In a systematic study of shock wave propagating in crystalline polyethylenes using molecular dynamics method and the electron force field (eFF) potential, we show that microscopic structure of shock front is significantly affected by the anisotropy of long carbon chain and the bond breaking and recombination dynamics. However, macroscopic properties measured in Hugoniot experiments, such as compression ratio and shock velocity, are not sensitive to carbon chain anisotropy and bond dynamics. Our work also display that hydrogen molecules are formed when the piston speed is in the region between 10 km/s and 30 km/s. However, carbon-hydrogen pair distribution function does not display an indication of carbon-hydrogen phase segregation.

Mesoscopic Simulation of the Two-Component System of Coupled Sine-Gordon Equations with Lattice Boltzmann Method

In this paper, a new lattice Boltzmann model for the two-component system of coupled sine-Gordon equations is presented by using the coupled mesoscopic Boltzmann equations. Via the Chapman-Enskog multiscale expansion, the macroscopical governing evolution system can be recovered correctly by selecting suitable discrete equilibrium distribution functions and the amending functions. The mesoscopic model has been validated by several related issues where analytic solutions are available. The experimental results show that the numerical results are consistent with the analytic solutions. From the mesoscopic point of view, the present approach provides a new way for studying the complex nonlinear partial differential equations arising in natural nonlinear phenomena of engineering and science.

Mesoscopic Simulation of the (2 + 1)-Dimensional Wave Equation with Nonlinear Damping and Source Terms Using the Lattice Boltzmann BGK Model

In this work, we develop a mesoscopic lattice Boltzmann Bhatnagar-Gross-Krook (BGK) model to solve (2 + 1)-dimensional wave equation with the nonlinear damping and source terms. Through the Chapman-Enskog multiscale expansion, the macroscopic governing evolution equation can be obtained accurately by choosing appropriate local equilibrium distribution functions. We validate the present mesoscopic model by some related issues where the exact solution is known. It turned out that the numerical solution is in very good agreement with exact one, which shows that the present mesoscopic model is pretty valid, and can be used to solve more similar nonlinear wave equations with nonlinear damping and source terms, and predict and enrich the internal mechanism of nonlinearity and complexity in nonlinear dynamic phenomenon.

Discrete Boltzmann modeling of unsteady reactive flows with nonequilibrium effects

A multiple-relaxation-time discrete Boltzmann model (DBM) is developed for compressible thermal reactive flows. A unified Boltzmann equation set is solved for hydrodynamic and thermodynamic quantities as well as higher order kinetic moments. The collision, reaction, and force terms are uniformly calculated with a matrix inversion method, which is physically accurate, numerically efficient, and convenient for coding. Via the Chapman-Enskog analysis, the DBM is demonstrated to recover reactive Navier-Stokes (NS) equations in the hydrodynamic limit. Both specific heat ratio and Prandtl number are adjustable. Moreover, it provides quantification of hydrodynamic and thermodynamic nonequilibrium effects beyond the NS equations. The capability of the DBM is demonstrated through simulations of chemical reactions in the free falling process, sound wave, thermal Couette flow, and steady and unsteady detonation cases. Moreover, nonequilibrium effects on the predicted physical quantities in unsteady combustion are quantified via the DBM. It is demonstrated that nonequilibrium effects suppress detonation instability and dissipate small oscillations of fluid flows.