Source author record

Aiguo Xu

Aiguo Xu 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

35works
16topics
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

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

preprint2022arXiv

Discrete Boltzmann modeling of Rayleigh-Taylor instability: effects of interfacial tension, viscosity and heat conductivity

The Rayleigh-Taylor Instability (RTI) in compressible flow with inter-molecular interactions is probed via the Discrete Boltzmann Method (DBM). The effects of interfacial tension, viscosity and heat conduction are investigated. It is found that the influences of interfacial tension on the perturbation amplitude, bubble velocity, and two kinds of entropy production rates all show differences at different stages of RTI evolution. It inhibits the RTI evolution at the bubble acceleration stage, while at the asymptotic velocity stage, it first promotes and then inhibits the RTI evolution. Viscosity and heat conduction inhibit the RTI evolution. Viscosity shows a suppressive effect on entropy generation rate related to heat flow at the early stage but a first promotive and then suppressive effect on entropy generation rate related to heat flow at a later stage. Heat conduction shows a promotive effect on entropy generation rate related to heat flow at an early stage. Still, it offers a first promotive and then suppressive effect on entropy generation rate related to heat flow at a later stage. By introducing the morphological boundary length, we found that the stage of exponential growth of interface length with time corresponds to the bubble acceleration stage. The first maximum point of interface length change rate and the first maximum point of the change rate of entropy generation rate related to viscous stress can be used as a new criterion for RTI to enter the asymptotic velocity stage.

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.

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

Morphological and non-equilibrium analysis of coupled Rayleigh-Taylor-Kelvin-Helmholtz instability

In this paper, the coupled Rayleigh-Taylor-Kelvin-Helmholtz instability(RTI, KHI and RTKHI, respectively) system is investigated using a multiple-relaxation-time discrete Boltzmann model. Both the morphological boundary length and thermodynamic nonequilibrium (TNE) strength are introduced to probe the complex configurations and kinetic processes. In the simulations, RTI always plays a major role in the later stage, while the main mechanism in the early stage depends on the comparison of buoyancy and shear strength. It is found that, both the total boundary length $L$ of the condensed temperature field and the mean heat flux strength $D_{3,1}$ can be used to measure the ratio of buoyancy to shear strength, and to quantitatively judge the main mechanism in the early stage of the RTKHI system. Specifically, when KHI (RTI) dominates, $L^{KHI} > L^{RTI}$ ($L^{KHI} < L^{RTI}$), $D_{3,1}^{KHI} > D_{3,1}^{RTI}$ ($D_{3,1}^{KHI} < D_{3,1}^{RTI}$); when KHI and RTI are balanced, $L^{KHI} = L^{RTI}$, $D_{3,1}^{KHI} = D_{3,1}^{RTI}$. A second sets of findings are as below: For the case where the KHI dominates at earlier time and the RTI dominates at later time, the evolution process can be roughly divided into two stages. Before the transition point of the two stages, $L^{RTKHI}$ initially increases exponentially, and then increases linearly. Hence, the ending point of linear increasing $L^{RTKHI}$ can work as a geometric criterion for discriminating the two stages. The TNE quantity, heat flux strength $D_{3,1}^{RTKHI}$, shows similar behavior. Therefore, the ending point of linear increasing $D_{3,1}^{RTKHI}$ can work as a physical criterion for discriminating the two stages.

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.

preprint2018arXiv

Comparative study of discrete Boltzmann model and Navier-Stokes

Discrete Boltzmann model (DBM) is a type of coarse-grained mesoscale kinetic model derived from the Boltzmann equation. Physically, it is roughly equivalent to a hydrodynamic model supplemented by a coarse-grained model for the relevant thermodynamic non-equilibrium (TNE) behaviours. The Navier-Stokes (NS) model is a traditional macroscopic hydrodynamic model based on continuity hypothesis and conservation laws. In this study, the two models are compared from two aspects, physical capability and computational cost, by simulating two kinds of flow problems including the thermal Couette flow and a Mach 3 step problem. In the cases where the TNE effects are weak, both the two models give accurate results for the hydrodynamic behaviour. Besides, DBM can provide more detailed non-equilibrium information, while the NS is more efficient if concern only the density, momentum, energy and their derived quantities. It is concluded that, if the TNE effects are strong or are to be investigated, the NS is insufficient while DBM is a good choice. While in the cases where the TNE effects are weak and only the macro flow fields are to be studied, the NS is more preferable.

preprint2016arXiv

Complex fields in heterogeneous materials under shock: modeling, simulation and analysis

In this mini-review we summarize the progress of modeling, simulation and analysis of shock responses of heterogeneous materials in our group in recent years. The basic methodology is as below. We first decompose the problem into different scales. Construct/Choose a model according to the scale and main mechanisms working at that scale. Perform numerical simulations using the relatively mature schemes. The physical information is transferred between neighboring scales in such a way: The statistical information of results in smaller scale contributes to establishing the constitutive equation in larger one. Except for the microscopic Molecular Dynamics (MD) model, both the mesoscopic and macroscopic models can be further classified into two categories, solidic and fluidic models, respectively. The basic ideas and key techniques of the MD, material point method and discrete Boltzmann method are briefly reviewed. Among various schemes used in analyzing the complex fields and structures, the morphological analysis and the home-built software, GISO, are briefly introduced. New observations are summarized for scales from the larger to the smaller.

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.

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

Discrete Boltzmann model of compressible flows with spherical or cylindrical symmetry

To study simultaneously the hydrodynamic and thermodynamic behaviors in compressible flow systems with spherical or cylindrical symmetry, we present a theoretical framework for constructing Discrete Boltzmann Model(DBM) with spherical or cylindrical symmetry in spherical or cylindrical coordinates. To this aim, a key technique is to use \emph{local} Cartesian coordinates to describe the particle velocity in the kinetic model. Thus, the geometric effects, like the divergence and convergence, are described as a \textquotedblleft force term\textquotedblright . Even though the hydrodynamic models are one- or two-dimensional, the DBM needs a Discrete Velocity Model(DVM) with 3 dimensions. We use a DVM with 26 velocities to formulate the DBM which recovers the Navier-Stokes equations with spherical or cylindrical symmetry in the hydrodynamic limit. For the system with \emph{global} cylindrical symmetry, we formulated also a DBM based on a DVM with 2 dimensions and 16 velocities. In terms of the nonconserved moments, we define two sets of measures for the deviations of the system from its thermodynamic equilibrium state. The extension of current model to the multiple-relaxation-time version is straightforward.

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.

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.

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.

preprint2013arXiv

Lattice Boltzmann model for combustion and detonation

In this paper we present a lattice Boltzmann model for combustion and detonation. In this model the fluid behavior is described by a finite-difference lattice Boltzmann model by Gan et al. [Physica A, 2008, 387: 1721]. The chemical reaction is described by the Lee-Tarver model [Phys. Fluids, 1980, 23: 2362]. The reaction heat is naturally coupled with the flow behavior. Due to the separation of time scales in the chemical and thermodynamic processes, a key technique for a successful simulation is to use the operator-splitting scheme. The new model is verified and validated by well-known benchmark tests. As a specific application of the new model, we studied the simple steady detonation phenomenon. To show the merit of LB model over the traditional ones, we focus on the reaction zone to study the non-equilibrium effects. It is interesting to find that, at the von Neumann peak, the system is nearly in its thermodynamic equilibrium. At the two sides of the von Neumann peak, the system deviates from its equilibrium in opposite directions. In the front of von Neumann peak, due to the strong compression from the reaction product behind the von Neumann peak, the system experiences a sudden deviation from thermodynamic equilibrium. Behind the von Neumann peak, the release of chemical energy results in thermal expansion of the matter within the reaction zone, which drives the system to deviate the thermodynamic equilibrium in the opposite direction. From the deviation from thermodynamic equilibrium, defined in this paper, one can understand more on the macroscopic effects of the system due to the deviation from its thermodynamic equilibrium.

preprint2013arXiv

Simulation study on cavity growth in ductile metal materials under dynamic loading

Cavity growth in ductile metal materials under dynamic loading is investigated via the material point method. Two typical cavity effects in the region subjected to rarefaction wave are identified: (i) part of material particles flow away from the cavity in comparison to the initial loading velocity, (ii) local regions show weaker negative or even positive pressures. Neighboring cavities interact via coalescence of isobaric contours. The growth of cavity under tension shows staged behaviors. After the initial slow stage, the volume and the dimensions in both the tensile and transverse directions show linear growth rate with time until the global tensile wave arrives at the upper free surface. It is interesting that the growth rate in the transverse direction is faster than that in the tensile direction. The volume growth rate linearly increases with the initial tensile velocity. After the global tensile wave passed the cavity, both the maximum particle velocity in the tensile direction and the maximum particle velocity in the opposite direction increase logarithmically with the initial tensile speed. The shock wave reflected back from the cavity and compression wave from the free surface induce the initial behavior of interfacial instabilities such as the Richtmyer-Meshkov instability, which is mainly responsible for the irregularity in the morphology of deformed cavity. The local temperatures and distribution of hot spots are determined by the plastic work. Compared with the dynamical process, the heat conduction is much slower.

preprint2013arXiv

Two dimensional MRT LB model for compressible and incompressible flows

In the paper we extend the Multiple-Relaxation-Time (MRT) Lattice Boltzmann (LB) model proposed in [Europhys. Lett. \textbf{90}, 54003 (2010)] so that it is suitable also for incompressible flows. To decrease the artificial oscillations, the convection term is discretized by the flux limiter scheme with splitting technique. New model is validated by some well-known benchmark tests, including Riemann problem and Couette flow, and satisfying agreements are obtained between the simulation results and analytical ones. In order to show the merit of LB model over traditional methods, the non-equilibrium characteristics of system are solved. The simulation results are consistent with physical analysis.

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.

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

Cluster Identification and Characterization of Physical Fields

The description of complex configuration is a difficult issue. We present a powerful technique for cluster identification and characterization. The scheme is designed to treat with and analyze the experimental and/or simulation data from various methods. Main steps are as follows. We first divide the space using face or volume elements from discrete points. Then, combine the elements with the same and/or similar properties to construct clusters with special physical characterizations. In the algorithm, we adopt administrative structure of hierarchy-tree for spatial bodies such as points, lines, faces, blocks, and clusters. Two fast search algorithms with the complexity are realized. The establishing of the hierarchy-tree and the fast searching of spatial bodies are general, which are independent of spatial dimensions. Therefore, it is easy to extend the skill to other fields. As a verification and validation, we treated with and analyzed some two-dimensional and three-dimensional random data.

preprint2010arXiv

Control spiral wave dynamics using feedback signals from line detectors

We numerically study trajectories of spiral-wave-cores in excitable systems modulated proportionally to the integral of the activity on the straight line, several or dozens of equi-spaced measuring points on the straight line, the double-line and the contour-line. We show the single-line feedback results in the drift of core center along a straight line being parallel to the detector. An interesting finding is that the drift location in $y$ is a piecewise linear-increasing function of both the feedback line location and time delay. Similar trajectory occurs when replacing the feedback line with several or dozens of equi-spaced measuring points on the straight line. This allows to move the spiral core to the desired location along a chosen direction by measuring several or dozens of points. Under the double-line feedback, the shape of the tip trajectory representing the competition between the first and second feedback lines is determined by the distance of two lines. Various drift attractors in spiral wave controlled by square-shaped contour-line feedback are also investigated. A brief explanation is presented.

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

Shock Wave Response of Porous Materials: From Plasticity to Elasticity

Shock wave reaction results in various characteristic regimes in porous material. The geometrical and topological properties of these regimes are highly concerned in practical applications. Via the morphological analysis to characteristic regimes with high temperature, we investigate the thermodynamics of shocked porous materials whose mechanical properties cover a wide range from hyperplasticity to elasticity. It is found that, under fixed shock strength, the total fractional area $A$ of the high-temperature regimes with $T \geq T_{th}$ and its saturation value first increase, then decrease with the increasing of the initial yield $σ_{Y0}$, where $T_{th}$ is a given threshold value of temperature $T$. In the shock-loading procedure, the fractional area $A(t)$ may show the same behavior if $T_{th}$ and $σ_{Y0}$ are chosen appropriately. Under the same $A(t)$ behavior, $T_{th}$ first increases then decreases with $σ_{Y0}$. At the maximum point $σ_{Y0M}$, the shock wave contributes the maximum plastic work. Around $σ_{Y0M}$, two materials with different mechanical properties may share the same $A(t)$ behavior even for the same $T_{th}$. The characteristic regimes in the material with the larger $σ_{Y0}$ are more dispersed.

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.

preprint2009arXiv

Simulation study of shock reaction on porous material

Direct modeling of porous materials under shock is a complex issue. We investigate such a system via the newly developed material-point method. The effects of shock strength and porosity size are the main concerns. For the same porosity, the effects of mean-void-size are checked. It is found that, local turbulence mixing and volume dissipation are two important mechanisms for transformation of kinetic energy to heat. When the porosity is very small, the shocked portion may arrive at a dynamical steady state; the voids in the downstream portion reflect back rarefactive waves and result in slight oscillations of mean density and pressure; for the same value of porosity, a larger mean-void-size makes a higher mean temperature. When the porosity becomes large, hydrodynamic quantities vary with time during the whole shock-loading procedure: after the initial stage, the mean density and pressure decrease, but the temperature increases with a higher rate. The distributions of local density, pressure, temperature and particle-velocity are generally non-Gaussian and vary with time. The changing rates depend on the porosity value, mean-void-size and shock strength. The stronger the loaded shock, the stronger the porosity effects. This work provides a supplement to experiments for the very quick procedures and reveals more fundamental mechanisms in energy and momentum transportation.

preprint2007arXiv

Three-dimensional multi-mesh material point method for solving collision problems

Contact algorithm between different bodies plays an important role in solving collision problems. Usually it is not easy to be treated very well. Several ones for material point method were proposed by Bardenhangen, Brackbill, and Sulsky\cite{Bardenhagen2000,Bardenhagen2001}, Hu and Chen\cite{Hu_Chen2003}. An improved one for three-dimensional material point method is presented in this paper. The improved algorithm emphasizes the energy conservation of the system and faithfully recovers opposite acting forces between contacting bodies. Contrasted to the one by Bardenhagen, both the normal and tangential contacting forces are more appropriately applied to the contacting bodies via the contacting nodes of the background mesh; Contrasted to the one by Hu and Chen, not only the tangential velocities but also the normal ones are handled separately in respective individual mesh. This treatment ensures not only the contact/sliding/separation procedure but also the friction between contacting bodies are recovered. The presented contact algorithm is validated via numerical experiments including rolling simulation, impact of elastic spheres, impact of a Taylor bar and impact of plastic spheres. The numerical results show that the multi-mesh material point method with the improved contact algorithm is more suitable for solving collision problems.

preprint2005arXiv

Scaling and hydrodynamic effects in lamellar ordering

We study the kinetics of domain growth of fluid mixtures quenched from a disordered to a lamellar phase. At low viscosities, in two dimensions, when hydrodynamic modes become important, dynamical scaling is verified in the form $C(\vec k, t) \sim L^α f[(k-k_M)L]$ where $C$ is the structure factor with maximum at $k_M$ and $L$ is a typical length changing from power law to logarithmic growth at late times. The presence of extended defects can explain the behavior of $L$. Three-dimensional simulations confirm that diffuse grain boundaries inhibit complete ordering of lamellae. Applied shear flow alleviates frustration and gives power-law growth at all times.

preprint2003arXiv

Phase separation of incompressible binary fluids with Lattice Boltzmann Methods

We introduce new versions of lattice Boltzmann methods (LBM) for incompressible binary mixtures where fluctuations of total density are inhibited. As a test for the improved algorithms we consider the problem of phase separation of two-dimensional binary mixtures quenched from a disordered state into the coexistence region. We find that the stability properties of LBM are greatly improved. The control of density fluctuations and the possibility of running longer simulations allow a more precise evaluation of the growth exponent.