Double-distribution-function discrete Boltzmann model for combustion

A 2-dimensional discrete Boltzmann model for combustion is presented. Mathematically, the model is composed of two coupled discrete Boltzmann equations for two species and a phenomenological equation for chemical reaction process. Physically, the model is equivalent to a reactive Navier-Stokes model supplemented by a coarse-grained model for the thermodynamic nonequilibrium behaviours. This model adopts 16 discrete velocities. It works for both subsonic and supersonic combustion phenomena with flexible specific heat ratio. To discuss the physical accuracy of the coarse-grained model for nonequilibrium behaviours, three other discrete velocity models are used for comparisons. Numerical results are compared with analytical solutions based on both the first-order and second-order truncations of the distribution function. It is confirmed that the physical accuracy increases with the increasing moment relations needed by nonequlibrium manifestations. Furthermore, compared with the single distribution function model, this model can simulate more details of combustion.Comment: Accepted for publication in Combustion and Flam

Numerical simulation and experimental study of PbWO4/EPDM and Bi2WO6/EPDM for the shielding of {\gamma}rays

The MCNP5 code was employed to simulate the {\gamma}ray shielding capacity of tungstate composites. The experimental results were applied to verify the applicability of the Monte Carlo program. PbWO4 and Bi2WO6 were prepared and added into ethylene propylene diene monomer (EPDM) to obtain the composites, which were tested in the {\gamma}ray shielding. Both the theoretical simulation and experiments were carefully chosen and well designed. The results of the two methods were found to be highly consistent. In addition, the conditions during the numerical simulation were optimized and double-layer {\gamma}ray shielding systems were studied. It was found that the {\gamma}-ray shielding performance can be influenced not only by the material thickness ratio but also by the arrangement of the composites.Comment: 8 pages,7 figures,Submitted to Chin.Phy.

Polar coordinate lattice Boltzmann modeling of compressible flows

We present a polar coordinate lattice Boltzmann kinetic model for compressible flows. A method to recover the continuum distribution function from the discrete distribution function is indicated. Within the model, a hybrid scheme being similar to, but different from, the operator splitting is proposed. The temporal evolution is calculated analytically, and the convection term is solved via a modifiedWarming-Beam (MWB) scheme.Within theMWB scheme a suitable switch function is introduced. The current model works not only for subsonic flows but also for supersonic flows. It is validated and verified via the following well-known benchmark tests: (i) the rotational flow, (ii) the stable shock tube problem, (iii) the Richtmyer-Meshkov (RM) instability, and (iv) the Kelvin-Helmholtz instability. As an original application, we studied the nonequilibrium characteristics of the system around three kinds of interfaces, the shock wave, the rarefaction wave, and the material interface, for two specific cases. In one of the two cases, the material interface is initially perturbed, and consequently the RMinstability occurs. It is found that themacroscopic effects due to deviating from thermodynamic equilibrium around thematerial interface differ significantly from those around the mechanical interfaces. The initial perturbation at the material interface enhances the coupling of molecular motions in different degrees of freedom. The amplitude of deviation from thermodynamic equilibrium around the shock wave is much higher than those around the rarefaction wave and material interface. By comparing each component of the high-order moments and its value in equilibrium, we can draw qualitatively the main behavior of the actual distribution function

Multiple-Relaxation-Time Lattice Boltzmann Approach to Compressible Flows with Flexible Specific-Heat Ratio and Prandtl Number

A new multiple-relaxation-time lattice Boltzmann scheme for compressible flows with arbitrary specific heat ratio and Prandtl number is presented. In the new scheme, which is based on a two-dimensional 16-discrete-velocity model, the moment space and the corresponding transformation matrix are constructed according to the seven-moment relations associated with the local equilibrium distribution function. In the continuum limit, the model recovers the compressible Navier-Stokes equations with flexible specific-heat ratio and Prandtl number. Numerical experiments show that compressible flows with strong shocks can be simulated by the present model up to Mach numbers $Ma \sim 5$.Comment: Accepted for publication in EP

Lattice Boltzmann study on Kelvin-Helmholtz instability: the roles of velocity and density gradients

A two-dimensional lattice Boltzmann model with 19 discrete velocities for compressible Euler equations is proposed (D2V19-LBM). The fifth-order Weighted Essentially Non-Oscillatory (5th-WENO) finite difference scheme is employed to calculate the convection term of the lattice Boltzmann equation. The validity of the model is verified by comparing simulation results of the Sod shock tube with its corresponding analytical solutions. The velocity and density gradient effects on the Kelvin-Helmholtz instability (KHI) are investigated using the proposed model. Sharp density contours are obtained in our simulations. It is found that, the linear growth rate $\gamma$ for the KHI decreases with increasing the width of velocity transition layer ${D_{v}}$ but increases with increasing the width of density transition layer ${D_{\rho}}$. After the initial transient period and before the vortex has been well formed, the linear growth rates, $\gamma_v$ and $\gamma_{\rho}$, vary with ${D_{v}}$ and ${D_{\rho}}$ approximately in the following way, $\ln\gamma_{v}=a-bD_{v}$ and $\gamma_{\rho}=c+e\ln D_{\rho} ({D_{\rho}}<{D_{\rho}^{E}})$, where $a$, $b$, $c$ and $e$ are fitting parameters and ${D_{\rho}^{E}}$ is the effective interaction width of density transition layer. When ${D_{\rho}}>{D_{\rho}^{E}}$ the linear growth rate $\gamma_{\rho}$ does not vary significantly any more. One can use the hybrid effects of velocity and density transition layers to stabilize the KHI. Our numerical simulation results are in general agreement with the analytical results [L. F. Wang, \emph{et al.}, Phys. Plasma \textbf{17}, 042103 (2010)].Comment: Accepted for publication in PR