Source author record

Yingjun Li

Yingjun Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Discrete Boltzmann modeling of high-speed compressible flows with various depths of non-equilibrium

The non-equilibrium high-speed compressible flows present wealthy applications in engineering and science. With the deepening of Thermodynamic Non-Equilibrium (TNE), higher-order non-conserved kinetic moments of the distribution function are needed to capture the main feature of the flow state and evolution process. Based on the ellipsoidal statistical Bhatnagar-Gross-Krook model, Discrete Boltzmann Models (DBMs) that consider various orders (from the first up to the sixth order) of TNE effects are developed to study flows in various depths of TNE. Specifically, at first, two types of one-dimensional Riemann problems and a Couette flow are used to show the model's capability to capture large flow structures with zero-order and first-order TNE effects, respectively. Then, a shock wave structure given by Direct simulation Monte Carlo is used to verify the model's capability to capture fine structures at the level of mean free path of molecules. Further, we focus on the TNE degree of two colliding fluids. A five-component vector $\mathbf{S}_{TNE} = (τ, Δ\mathbf{u}, ΔT, \bm{Δ_{2}^{*}},\bm{Δ_{3,1}^{*}})$ is introduced to roughly characterize the TNE degree. It is found that the TNE strengths obtained from various perspectives are different. These findings demonstrate that the inadequacy of focusing only on the few kinetic moments appearing in Navier-Stokes increases with the degree of discreteness and deviation from thermodynamic equilibrium. Finally, a two-dimensional free jet is simulated to indicate that, to obtain satisfying hydrodynamic quantities, the DBM should include at least up to the third-order TNE effects.

preprint2021arXiv

Delineation of the flow and mixing induced by Rayleigh-Taylor instability through tracers

Rayleigh-Taylor-instability(RTI) induced flow and mixing are of great importance in both nature and engineering scenarios. To capture the underpinning physics, tracers are introduced to make a supplement to discrete Boltzmann simulation of RTI in compressible flows. Via marking two types of tracers with different colors, the tracer distribution provides a clear boundary of two fluids during the RTI evolution. Fine structures of the flow and thermodynamic nonequilibrium behavior around the interface in a miscible two-fluid system are delineated. Distribution of tracers in its velocity phase space makes a charming pattern showing quite dense information on the flow behavior, which opens a new perspective for analyzing and accessing significantly deep insights into the flow system. RTI mixing is further investigated via tracer defined local mixedness. The appearance of Kelvin-Helmholtz instability is quantitatively captured by mixedness averaged align the direction of the pressure gradient. The role of compressibility and viscosity on mixing are investigated separately, both of which show two-stage effect. The underlying mechanism of the two-stage effect is interpreted as the development of large structures at the initial stage and the generation of small structures at the late stage. At the late stage, for a fixed time, a saturation phenomenon of viscosity is found that further increase of viscosity cannot see an evident decline in mixedness. The mixing statues of heavy and light fluids are not synchronous and the mixing of a RTI system is heterogenous. The results are helpful for understanding the mechanism of flow and mixing induced by RTI.

preprint2020arXiv

Varied fusion reaction probability induced by ion stopping modification in laser-driven plasma with different temperature

The dynamics of nuclear reaction in plasma is a fundamental issue in many high energy density researches, such as the astrophysical reactions and the inertial confinement fusion. The effective reaction cross-sections and ion stopping power in plasma need to be taken into account to analyze the reactivity. In this research, we have experimentally investigated the from D-D reactions from interactions between deuteron beams and deuterated polystyrene (CD) plasma, driven by two laser pulses respectively. The neutron yields, plasma density and deuteron energy loss in plasma have been measured, and the plasma temperature and deuteron stopping power have been analyzed from simulations. It is shown that, compared with a cold target, the reaction probability in plasma conditions can be enhanced or suppressed, which is ascribed to the deuteron stopping power modifications in plasma. In hotter CD plasma, the energy loss of moderate energetic deuterons reduces, which leads to higher D-D reaction probability, while the contrary happens in colder plasma. This work provides new understanding of fusion reactions in plasma environment.

preprint2016arXiv

Numerical simulation and experimental study of PbWO4/EPDM and Bi2WO6/EPDM for the shielding of γrays

The MCNP5 code was employed to simulate the γ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 γ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 γray shielding systems were studied. It was found that the γ-ray shielding performance can be influenced not only by the material thickness ratio but also by the arrangement of the composites.

preprint2015arXiv

Multiple-relaxation-time lattice Boltzmann kinetic model for combustion

To probe both the Hydrodynamic Non-Equilibrium (HNE) and Thermodynamic Non-Equilibrium (TNE) in the combustion process, a two-dimensional Multiple-Relaxation-Time (MRT) version of Lattice Boltzmann Kinetic Model(LBKM) for combustion phenomena is presented. The chemical energy released in the progress of combustion is dynamically coupled into the system by adding a chemical term to the LB kinetic equation. Beside describing the evolutions of the conserved quantities, the density, momentum and energy, which are what the Navier-Stokes model describes, the MRT-LBKM presents also a coarse-grained description on the evolutions of some non-conserved quantities. The current model works for both subsonic and supersonic flows with or without chemical reaction. In this model both the specific-heat ratio and the Prandtl number are flexible, the TNE effects are naturally presented in each simulation step. The model is verified and validated via well-known benchmark tests. As an initial application, various non-equilibrium behaviours, including the complex interplays between various HNEs, between various TNEs and between the HNE and TNE, around the detonation wave in the unsteady and steady one-dimensional detonation processes are preliminarily probed. It is found that the system viscosity (or heat conductivity) decreases the local TNE, but increase the global TNE around the detonation wave, that even locally, the system viscosity (or heat conductivity) results in two kinds of competing trends, to increase and to decrease the TNE effects. The physical reason is that the viscosity (or heat conductivity) takes part in both the thermodynamic and hydrodynamic responses.

preprint2014arXiv

Polar coordinate lattice Boltzmann kinetic modeling of detonation phenomena

A novel polar coordinate lattice Boltzmann kinetic model for detonation phenomena is presented and applied to investigate typical implosion and explosion processes. In this model, the change of discrete distribution function due to local chemical reaction is dynamically coupled into in the modified lattice Boltzmann equation, which could recovery the Navier-Stokes equations, including contribution of chemical reaction, via the Chapman-Enskog expansion. For the numerical investigations, the main focuses are the nonequilibrium behaviors in these processes. The system at the disc center is always in its thermodynamic equilibrium. The internal kinetic energies in different degrees freedoms around the detonation front do not coincide due to the fluid viscosity. They show the maximum difference at the inflexion point where the pressure has the largest spatial derivative. The dependence of the reaction rate on the pressure, influences of the shock strength and reaction rate on the departure amplitude of the system from its local thermodynamic equilibrium are probed.

preprint2011arXiv

FFT-LB modeling of thermal liquid-vapor systems

We further develop a thermal LB model for multiphase flows. In the improved model, we propose to use the FFT scheme to calculate both the convection term and external force term. The usage of FFT scheme is detailed and analyzed. By using the FFT algorithm spatiotemporal discretization errors are decreased dramatically and the conservation of total energy is much better preserved. A direct consequence of the improvement is that the unphysical spurious velocities at the interfacial regions can be damped to neglectable scale. Together with the better conservation of total energy, the more accurate flow velocities lead to the more accurate temperature field which determines the dynamical and final states of the system. With the new model, the phase diagram of the liquid-vapor system obtained from simulation is more consistent with that from theoretical calculation. Very sharp interfaces can be achieved. The accuracy of simulation results are also verified by the Laplace law. The FFT scheme can be easily applied to other models for multiphase flows.

preprint2011arXiv

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 $γ$ 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_ρ}$. After the initial transient period and before the vortex has been well formed, the linear growth rates, $γ_v$ and $γ_ρ$, vary with ${D_{v}}$ and ${D_ρ}$ approximately in the following way, $\lnγ_{v}=a-bD_{v}$ and $γ_ρ=c+e\ln D_ρ ({D_ρ}<{D_ρ^{E}})$, where $a$, $b$, $c$ and $e$ are fitting parameters and ${D_ρ^{E}}$ is the effective interaction width of density transition layer. When ${D_ρ}>{D_ρ^{E}}$ the linear growth rate $γ_ρ$ 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)].

preprint2011arXiv

Phase separation in thermal systems: LB study and morphological characterization

We investigate thermal and isothermal symmetric liquid-vapor separations via a FFT-Thermal Lattice Boltzmann (FFT-TLB) model. Structure factor, domain size and Minkowski functionals are employed to characterize the density and velocity fields as well as to understand the configurations and the kinetic processes. Compared with the isothermal phase separation, the freedom in temperature prolongs the Spinodal Decomposition (SD) stage and induces different rheological and morphological behaviors in the thermal system. After the transient procedure, both the thermal and isothermal separations show power-law scalings in domain growth; while the exponent for thermal system is lower than that for isothermal system. With respect to the density of field, the isothermal system presents more likely bicontinuous configurations with narrower interfaces, while the thermal system presents more likely configurations with scattered bubbles. Heat creation, conduction and lower interfacial stresses are main reasons for the differences in thermal system. Different from the case with isothermal phase separation, the release of latent heat causes the changing of local temperature which results in new local mechanical balance. When the Prandtl number becomes smaller, the system approaches thermodynamical equilibrium more quickly. The increasing of mean temperature makes lower the interfacial stress in the following way: $σ=σ_{0}[(T_{c}-T)/(T_{c}-T_{0})]^{3/2}$, where $T_{c}$ is the critical temperature and $σ_{0}$ is the interfacial stress at a reference temperature $T_{0}$, which is the main reason for lower growth exponent in thermal case.

preprint2011arXiv

Prandtl number effects in MRT Lattice Boltzmann models for shocked and unshocked compressible fluids

For compressible fluids under shock wave reaction, we have proposed two Multiple-Relaxation-Time (MRT) Lattice Boltzmann (LB) models [F. Chen, et al, EPL \textbf{90} (2010) 54003; Phys. Lett. A \textbf{375} (2011) 2129.]. In this paper, we construct a new MRT Lattice Boltzmann model which is not only for the shocked compressible fluids, but also for the unshocked compressible fluids. To make the model work for unshocked compressible fluids, a key step is to modify the collision operators of energy flux so that the viscous coefficient in momentum equation is consistent with that in energy equation even in the unshocked system. The unnecessity of the modification for systems under strong shock is analyzed. The model is validated by some well-known benchmark tests, including (i) thermal Couette flow, (ii) Riemann problem, (iii) Richtmyer-Meshkov instability. The first system is unshocked and the latter two are shocked. In all the three systems, the Prandtl numbers effects are checked. Satisfying agreements are obtained between new model results and analytical ones or other numerical results.

preprint2010arXiv

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$.

preprint2010arXiv

Multiple-relaxation-time lattice Boltzmann model for compressible fluids

We present an energy-conserving multiple-relaxation-time finite difference lattice Boltzmann model for compressible flows. This model is based on a 16-discrete-velocity model. The collision step is first calculated in the moment space and then mapped back to the velocity space. The moment space and corresponding transformation matrix are constructed according to the group representation theory. Equilibria of the nonconserved moments are chosen according to the need of recovering compressible Navier-Stokes equations through the Chapman-Enskog expansion. Numerical experiments showed that compressible flows with strong shocks can be well simulated by the present model. The used benchmark tests include (i) shock tubes, such as the Sod, Lax, Sjogreen, Colella explosion wave and collision of two strong shocks, (ii) regular and Mach shock reflections, and (iii) shock wave reaction on cylindrical bubble problems. The new model works for both low and high speeds compressible flows. It contains more physical information and has better numerical stability and accuracy than its single-relaxation-time version.

preprint2010arXiv

Three-Dimensional Lattice Boltzmann Model for High-Speed Compressible Flows

A highly efficient three-dimensional (3D) Lattice Boltzmann (LB) model for high speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara[Phys. Rev. E 69, 056702 (2004)]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative (NND) scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, (i) Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; (ii) reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.

preprint2009arXiv

Simulating liquid-vapor phase separation under shear with lattice Boltzmann method

We study liquid-vapor phase separation under shear via the Shan-Chen lattice Boltzmann model. Besides the rheological characteristics, we analyze the Kelvin-Helmholtz(K-H) instability resulting from the tangential velocity difference of the fluids on two sides of the interface. We discuss also the growth behavior of droplets. The domains being close to the walls are lamellar-ordered, where the hydrodynamic effects dominate. The patterns in the bulk of the system are nearly isotropic, where the domain growth results mainly from the diffusion mechanism. Both the interfacial tension and the K-H instability make the liquid-bands near the walls tend to rupture. When the shear rate increases, the inequivalence of evaporation in the upstream and coagulation in the downstream of the flow as well as the role of surface tension makes the droplets elongate obliquely. Stronger convection makes easier the transferring of material particles so that droplets become larger.