Source author record

Yanbiao Gan

Yanbiao Gan 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

12works
7topics
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

12 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

Effects of the initial perturbations on the Rayleigh-Taylor-Kelvin-Helmholtz instability system

In the paper, the effects of initial perturbations on the Rayleigh-Taylor instability (RTI), Kelvin-Helmholtz instability (KHI), and the coupled Rayleigh-Taylor-Kelvin-Helmholtz instability (RTKHI) systems are investigated using a multiple-relaxation-time discrete Boltzmann model. Six different perturbation interfaces are designed to study the effects of the initial perturbations on the instability systems. Based on the mean heat flux strength $D_{3,1}$, the effects of initial interfaces on the coupled RTKHI are examined in detail. The research is focused on two aspects: (i) the main mechanism in the early stage of the RTKHI, (ii) the transition point from KHI-like to RTI-like for the case where the KHI dominates at earlier time and the RTI dominates at later time. It is found that the early main mechanism is related to the shape of the initial interface, which is represented by both the bilateral contact angle $θ_{1}$ and the middle contact angle $θ_{2}$. The influence of inverted parabolic and inverted ellipse perturbations ($θ_{1}<90$) on the transition point of the RTKHI system is greater than that of other interfaces.

preprint2020arXiv

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). It is found that, in the dynamic KHI process, the mixing degree and fluid flow are similar for cases with various thermal conductivity and initial temperature configurations. Physically, both heat conduction and temperature exert slight influences on the formation and evolution of the KHI.

preprint2016arXiv

Thermo-hydrodynamic non-equilibrium effects on compressible Rayleigh-Taylor instability

The effects of compressibility on Rayleigh-Taylor instability (RTI) are investigated by inspecting the interplay between thermodynamic and hydrodynamic non-equilibrium phenomena (TNE, HNE, respectively) via a discrete Boltzmann model (DBM). Two effective approaches are presented, one tracking the evolution of the \emph{local} TNE effects and the other focussing on the evolution of the mean temperature of the fluid, to track the complex interfaces separating the bubble and the spike regions of the flow. It is found that, both the compressibility effects and the \emph{global} TNE intensity show opposite trends in the initial and the later stages of the RTI. Compressibility delays the initial stage of RTI and accelerates the later stage. Meanwhile, the TNE characteristics are generally enhanced by the compressibility, especially in the later stage. The global or mean thermodynamic non-equilibrium indicators provide physical criteria to discriminate between the two stages of the RTI.

preprint2014arXiv

Discrete Boltzmann modeling of liquid-vapor system

We further probe the Discrete Boltzmann Modeling(DBM) of the single-component two phase flows or the liquid-vapor system. There are two kinds of nonequilibrium effects in the system. The first is the Mechanical NonEquilibrium(MNE). The second is the Thermodynamic NonEquilibrium(TNE). The MNE is well described in the traditional fluid dynamic theory. The description of TNE resorts to the gas kinetic theory. Since based on the Boltzmann equation, the DBM makes possible to analyze both the MNE and TNE. The TNE is the main discussion of this work. A major purpose of this work is to show that the DBM results can be used to confirm and/or improve the macroscopic modeling of complex system.

preprint2014arXiv

Lattice Boltzmann kinetic modeling and simulation of thermal liquid-vapor system

We present a highly efficient lattice Boltzmann (LB) kinetic model for thermal liquid-vapor system. Three key components are as beow: (i) a discrete velocity model by Kataoka \emph{et al.} [Phys. Rev. E \textbf{69}, 035701(R)(2004)]; (ii) a forcing term $I_{i}$ aiming to describe the interfacial stress and recover the van der Waals equation of state by Gonnella \emph{et al.} [Phys. Rev. E \textbf{76}, 036703 (2007)]; and (iii) a Windowed Fast Fourier Transform (WFFT) scheme and its inverse by our group [Phys. Rev. E \textbf{84}, 046715 (2011)] for solving the spatial derivatives, together with a second-order Runge-Kutta (RK) finite difference scheme for solving the temporal derivative in the LB equation. The model is verified and validated by well-known benchmark tests. The results recovered from the present model are well consistent with previous ones[Phys. Rev. E \textbf{84}, 046715 (2011)] or theoretical analysis. The usage of less discrete velocities, high-order RK algorithm and WFFT scheme with 16th-order in precision makes the model more efficient by about $10$ times and more accurate than the original one.

preprint2013arXiv

Lattice BGK kinetic model for high speed compressible flows: hydrodynamic and nonequilibrium behaviors

We present a simple and general approach to formulate the lattice BGK model for high speed compressible flows. The main point consists of two parts: an appropriate discrete equilibrium distribution function (DEDF) $\mathbf{f}^{eq}$ and a discrete velocity model with flexible velocity size. The DEDF is obtained by $\mathbf{f}^{eq}=\mathbf{C}^{-1}\mathbf{M}$, where $\mathbf{M}$ is a set of moment of the Maxwellian distribution function, and $\mathbf{C}$ is the matrix connecting the DEDF and the moments. The numerical components of $\mathbf{C}$ are determined by the discrete velocity model. The calculation of $\mathbf{C}^{-1}$ is based on the analytic solution which is a function of the parameter controlling the sizes of discrete velocity. The choosing of discrete velocity model has a high flexibility. The specific heat ratio of the system can be flexible. The approach works for the one-, two- and three-dimensional model constructions. As an example, we compose a new lattice BGK kinetic model which works not only for recovering the Navier-Stokes equations in the continuum limit but also for measuring the departure of system from its thermodynamic equilibrium. Via adjusting the sizes of the discrete velocities the stably simulated Mach number can be significantly increased up to 30 or even higher. The model is verified and validated by well-known benchmark tests. Some macroscopic behaviors of the system due to deviating from thermodynamic equilibrium around the shock wave interfaces are shown.

preprint2012arXiv

Cellular automata model for elastic solid material

The Cellular Automaton (CA) modeling and simulation of solid dynamics is a long-standing difficult problem. In this paper we present a new two-dimensional CA model for solid dynamics. In this model the solid body is represented by a set of white and black particles alternatively positioned in the $x$- and $y$- directions. The force acting on each particle is represented by the linear summation of relative displacements of the nearest-neighboring particles. The key technique in this new model is the construction of eight coefficient matrices. Theoretical and numerical analyses show that the present model can be mathematically described by a conservative system. So, it works for elastic material. In the continuum limit the CA model recovers the well-known Navier equations. The coefficient matrices are related to the shear module and Poisson ratio of the material body. Compared with previous CA model for solid body, this model realizes the natural coupling of deformations in the $x$- and $y$- directions. Consequently, the wave phenomena related to the Poisson ratio effects are successfully recovered. This work advances significantly the CA modeling and simulation in the field of computational solid dynamics.

preprint2012arXiv

Lattice Boltzmann modeling and simulation of compressible flows

In this mini-review we summarize the progress of Lattice Boltzmann(LB) modeling and simulating compressible flows in our group in recent years. Main contents include (i) Single-Relaxation-Time(SRT) LB model supplemented by additional viscosity, (ii) Multiple-Relaxation-Time(MRT) LB model, and (iii) LB study on hydrodynamic instabilities. The former two belong to improvements of physical modeling and the third belongs to simulation or application. The SRT-LB model supplemented by additional viscosity keeps the original framework of Lattice Bhatnagar-Gross-Krook (LBGK). So, it is easier and more convenient for previous SRT-LB users. The MRT-LB is a completely new framework for physical modeling. It significantly extends the range of LB applications. The cost is longer computational time. The developed SRT-LB and MRT-LB are complementary from the sides of convenience and applicability.

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.