BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Guangcai Zhang

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

ResearcherUnclaimed source record
cond-mat.softcond-mat.mtrl-sciphysics.flu-dynphysics.comp-phnlin.CGcond-mat.stat-mechcond-mat.otherphysics.plasm-phComputationComputational Engineering, Finance, and ScienceComputational Geometrycond-mat.mes-hallNumerical Analysisphysics.class-ph

Catalog footprint

What is connected

28works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

28 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The different roles of Pu-oxide overlayers in the hydrogenation of Pu-metal: An ab initio molecular dynamics study based on vdW-DFT+U

Based on the van der Waals density functional theory (vdW-DFT)+U scheme, we carry out the ab initio molecular dynamics (AIMD) study of the interaction dynamics for H$_{2}$ impingement against the stoichiometric PuO$_{2}$(111), the reduced PuO$_{2}$(111), and the stoichiometric $α$-Pu$_{2}$O$_{3}$(111) surfaces. The hydrogen molecular physisorption states, which can not be captured by pure DFT+\textit{U} method, are obtained by employing the vdW-DFT+\textit{U} scheme. We show that except for the weak physisorption, PuO$_{2}$(111) surfaces are so difficult of access that almost all of the H$_{2}$ molecules will bounce back to the vacuum when their initial kinetic energies are not sufficient. Although the dissociative adsorption of H$_{2}$ on PuO$_{2}$(111) surfaces is found to be very exothermic, the collision-induced dissociation barriers of H$_{2}$ are calculated to be as high as $3.2$ eV and $2.0$ eV for stoichiometric and reduced PuO$_{2}$ surfaces, respectively. Unlike PuO$_{2}$, our AIMD study directly reveals that the hydrogen molecules can penetrate into $α$-Pu$_{2}$O$_{3}$(111) surface and diffuse easily due to the $25$\ native O vacancies located along the $\langle $111$\rangle $ diagonals of $α$-Pu$_{2}$O$_{3}$ matrix. By examining the temperature effect and the internal vibrational excitations of H$_{2}$, we provide a detailed insight into the interaction dynamics of H$_{2}$ in $α$-Pu$_{2}$O$_{3}$. The optimum pathways for hydrogen penetration and diffusion, the corresponding energy barriers ($1.0$ eV and $0.53$ eV, respectively) and rate constants are systematically calculated. Overall, our study fairly reveals the different interaction mechanisms between H$_{2}$ and Pu-oxide surfaces, which have strong implications to the interpretation of experimental observations.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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)].

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Molecular dynamics study of hcp/fcc nucleation and growth in bcc iron driven by uniaxial strain

Molecular dynamics simulations are performed to investigate the structural phase transition in body-centered cubic (bcc) single crystal iron under high strain rate loading. We study the nucleation and growth of the hexagonal-close-packed (hcp) and face-centered-cubic (fcc) phases, and their crystal orientation dependence. Results reveal that the transition pressures are less dependent on the crystal orientations ($\mathtt{\sim}$14 GPa for loading along [001], [011], and [111] directions). However, the pressure interval of mixed phase for [011] loading is much shorter than loading along other orientations. And the temperature increased amplitude for [001] loading is evidently lower than other orientations. The hcp/fcc nucleation process is presented by the topological medium-range-order analysis. For loading along [001] direction, we find that the hcp structure occurs firstly and grows into laminar morphology in the (011)$_{\text{bcc}}$ planes with a little fcc atoms as intermediate structure. For loading along [011] and [111] directions, both the hcp and fcc structures nucleation and growth along the \{110\}$_{\text{bcc}}$ planes are observed, whose morphology is also discussed.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access