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

Zhaoli Guo

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

ResearcherUnclaimed source record
physics.flu-dynphysics.comp-phComputational Engineering, Finance, and Sciencecond-mat.softmath.NANumerical Analysisphysics.geo-ph

Catalog footprint

What is connected

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

19 published item(s)

preprint2022arXiv

Discrete unified gas kinetic scheme for the conservative Allen-Cahn equation

In this paper, the discrete unified gas kinetic scheme (DUGKS) with an improved microflux across the cell interface for the conservative Allen-Cahn equation (CACE) is proposed. In the context of DUGKS, the recovered kinetic equation from the flux evaluation with linear reconstruction in the previous DUGKS is analyzed. It is found that the calculated microflux across the cell interface is only the solution to the target kinetic equation with first order accuracy, which can result in an inaccurate CACE since the force term is involved or the first moment of the collision model has no conservation property. To correctly recover the kinetic equation up to the second order accuracy, the value of the distribution function that will propagate along the characteristic line with ending point at the cell interface is appropriated by the parabolic reconstruction instead of the linear reconstruction. To validate the accuracy of the present DUGKS for the CACE, several benchmark problems, including the diagonal translation of a circular interface, the rotation of a Zaleska disk and the deformation of a circular interface, have been simulated. Numerical results show that the present DUGKS scheme is able to capture the interface with improved accuracy when compared with the previous DUGKS.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-equilibrium transport of inhomogeneous shale gas under ultra-tight confinement

The non-equilibrium transport of inhomogeneous and dense gases highly confined by surface is encountered in many engineering applications. For example, in the shale gas production process, methane is extracted from ultra-tight pores under high pressure so the gas is inhomogeneous and dense. Currently, the complex non-equilibrium transport of inhomogeneous and dense gases where gas surface interactions play a key role is commonly investigated by molecular dynamics or on a continuum-assumption basis. Here, a tractable kinetic model based on the generalized Enskog equation and the mean-field theory is employed to couple the effects of the volume exclusion and the long-range intermolecular attraction forces. The interactions between gas molecules and confined surface are modelled by a 10-4-3 Lennard-Jones potential, which can capture gas surface adsorption. The cross-sectional density profiles of methane under different confinements are in good agreement with the molecular dynamics results reported in the literature, and the transport behaviors are validated by the non-equilibrium molecular dynamics. The velocity of methane flow in shale matrix is plug-like due to its dense characteristics in nanopores. The influence of pressure, temperature, pore size and shale composition on density and velocity profiles is analyzed quantitatively. Our results show that the Klinkenberg correction is not applicable to model shale gas flow in the production process; the Navier-Stokes model using the second-order slip boundary condition cannot produce the proper velocity profiles, and consequently fails to predict the accurate flow rate in nanopores. This study sheds new light on understanding the physics of non-equilibrium dense gas flows in shale strata.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-ideal gas dynamics under confinement: rarefaction effect, dense effect and molecular interaction

The effects of volume exclusion and long-range intermolecular attraction are investigated by the simplified kinetic model for surface-confined inhomogeneous fluids. Gas dynamics of the ideal gas, the hard-sphere fluid and the real gas are simulated by the Boltzmann equation, the Enskog equation and the simple kinetic equation, respectively. Only the Knudsen minimum appears for the ideal gas, while both the Knudsen minimum and the Knudsen maximum occur for the hard-sphere fluid and the real gas under certain confinements, beyond which the maximum and minimum may disappear. The Boltzmann equation and the Enskog equation overestimates and underestimates the mass flow rate of the real gas dynamics under confinement, respectively, where the volume exclusion and the long-range intermolecular attractive potential among molecules are not ignorable. With the increase of the channel width, gas dynamics of the hard-sphere fluid and the real gas tends to the Boltzmann prediction gradually. The density inhomogeneity, which hinders the flow under confinement, is more obvious when the solid fraction is larger. The anomalous slip occurs for real gas under constant confinement. The flow at a smaller Knudsen number (larger solid fraction or channel width) contributes more practical amount of mass transfer, although the rarefaction effects is more prominent at larger Knudsen numbers. The temperature has no effect on density and velocity profiles of the ideal gas and the hard-sphere fluid, but the energy parameter among the real gas molecules decreases with the increasing temperature and the real gas dynamics tends to the hard-sphere ones consequently.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A well-balanced gas kinetic scheme for Navier-Stokes equations with gravitational potential

The hydrostatic equilibrium state is the consequence of the exact hydrostatic balance between hydrostatic pressure and external force. Standard finite volume or finite difference schemes cannot keep this balance exactly due to their unbalanced truncation errors. In this study, we introduce an auxiliary variable which becomes constant at isothermal hydrostatic equilibrium state and propose a well-balanced gas kinetic scheme for the Navier-Stokes equations with a global reconstruction. Through reformulating the convection term and the force term via the auxiliary variable, zero numerical flux and zero numerical source term are enforced at the hydrostatic equilibrium state instead of the balance between hydrostatic pressure and external force. Several problems are tested numerically to demonstrate the accuracy and the stability of the new scheme, and the results confirm that, the new scheme can preserve the exact hydrostatic solution. The small perturbation riding on hydrostatic equilibria can be calculated accurately. The viscous effect is also illustrated through the propagation of small perturbation and the Rayleigh-Taylor instability. More importantly, the new scheme is capable of simulating the process of converging towards hydrostatic equilibrium state from a highly non-balanced initial condition. The ultimate state of zero velocity and constant temperature is achieved up to machine accuracy. As demonstrated by the numerical experiments, the current scheme is very suitable for small amplitude perturbation and long time running under gravitational potential.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

GPU acceleration of an iterative scheme for gas-kinetic model equations with memory reduction techniques

This paper presents a Graphics Processing Units (GPUs) acceleration method of an iterative scheme for gas-kinetic model equations. Unlike the previous GPU parallelization of explicit kinetic schemes, this work features a fast converging iterative scheme. The memory reduction techniques in this method enable full three-dimensional (3D) solution of kinetic model equations in contemporary GPUs usually with a limited memory capacity that otherwise would need terabytes of memory. The GPU algorithm is validated against the DSMC simulation of the 3D lid-driven cavity flow and the supersonic rarefied gas flow past a cube with grids size up to 0.7 trillion points in the phase space. The performance of the GPU algorithm is assessed by comparing with the corresponding parallel CPU program using Message Passing Interface (MPI). The profiling on several models of GPUs shows that the algorithm has a medium to high level of utilization of the GPUs' computing and memory resources. A $190\times$ speedup can be achieved on the Tesla K40 GPUs against a single core of Intel Xeon-E5-2680v3 CPU for the 3D lid-driven cavity flow.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A comparative study of discrete velocity methods for rarefied gas flows

In the study of rarefied gas dynamics, the discrete velocity method (DVM) has been widely employed to solve the gas kinetic equations. Although various versions of DVM have been developed, their performance, in terms of accuracy and computational efficiency, is yet to be compreheively studied in the whole flow regime. Here, the traditional third-order time-implicit Godunov DVM (GDVM) and the recently developed discrete unified gas-kinetic scheme (DUGKS) are analysed in finding steady-state solutions of the force-driven Poiseuille and lid-driven cavity flows. With the molecular collision and free streaming being treated simultaneously, the DUGKS preserves the second-order accuracy in the spatial and temporal discretizations in all flow regimes. Towards the hydrodynamic flow regime, the DUGKS is not only faster than the GDVM when using the same spatial mesh, but also requires less spatial resolution than that of the GDVM to achieve the same numerical accuracy. From the slip to free molecular flow regimes, however, the DUGKS is slower than the GDVM, due to the complicated flux evaluation and the time step is less than the maximum effective time step of the GDVM. Therefore, the DUGKS is preferable for problems involving different flow regimes, particularly when the hydrodynamic flow regime is dominant. For highly rarefied gas flows, if the steady-state solution is concerned, the implicit DVM, which can boost the convergence significantly, is a better choice.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A simplification of the unified gas kinetic scheme

Unified gas kinetic scheme (UGKS) is an asymptotic preserving scheme for the kinetic equations. It is superior for transition flow simulations, and has been validated in the past years. However, compared to the well known discrete ordinate method (DOM) which is a classical numerical method solving the kinetic equations, the UGKS needs more computational resources. In this study, we propose a simplification of the unified gas kinetic scheme. It allows almost identical numerical cost as the DOM, but predicts numerical results as accurate as the UGKS. Based on the observation that the equilibrium part of the UGKS fluxes can be evaluated analytically, the equilibrium part in the UGKS flux is not necessary to be discretized in velocity space. In the simplified scheme, the numerical flux for the velocity distribution function and the numerical flux for the macroscopic conservative quantities are evaluated separately. The simplification is equivalent to a flux hybridization of the gas kinetic scheme for the Navier-Stokes (NS) equations and conventional discrete ordinate method. Several simplification strategies are tested, through which we can identify the key ingredient of the Navier-Stokes asymptotic preserving property. Numerical tests show that, as long as the collision effect is built into the macroscopic numerical flux, the numerical scheme is Navier-Stokes asymptotic preserving, regardless the accuracy of the microscopic numerical flux for the velocity distribution function.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comparison of the LBE and DUGKS methods for DNS of decaying homogeneous isotropic turbulence

The main objective of this work is to perform a detailed comparison of the lattice Boltzmann equation (LBE) and the recently developed discrete unified gas-kinetic scheme (DUGKS) methods for direct numerical simulation (DNS) of the decaying homogeneous isotropic turbulence (DHIT) in a periodic box. The flow fields and key statistical quantities computed by both methods are compared with those from pseudo-spectral (PS) method. The results show that the LBE and DUGKS have almost the same accuracy when the flow field is well-resolved, and that the LBE is less dissipative and is slightly more efficient than the DUGKS, but the latter has a superior numerical stability, particularly for high Reynolds number flows. Therefore, the DUGKS method can be viewed as a viable tool for DNS of turbulent flows. It should be emphasized that the main advantage of the DUGKS when compared with the LBE method is its feasibility in adopting nonuniform meshes, which is critical for wall-bounded turbulent flows. The present work provides a basis for further applications of DUGKS in studying the physics of the turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Lattice Boltzmann simulation of separation phenomenon in a binary gaseous flow through a microchannel

Gas separation of a binary gaseous mixture is one of characteristic phenomena in the micro-scale flows that differ from the conventional size flows. In this work, the separation in a binary gas mixture flows through a microchannel is investigated by the lattice Boltzmann method with a diffuse-bounce-back (DBB) boundary condition. The separation degree and rate are measured in the He--Ar and Ne--Ar systems for different mole fractions, pressure ratios, and Knudsen numbers. The results show that the separation phenomenon in the He--Ar mixture is more obvious than that in the Ne--Ar mixture at the same mole fraction owing to the larger molecular mass ratio. In addition, the increase in the pressure ratio reduces the difference in the molecular velocities between the two species, and the separation phenomenon becomes weaker. However, the gas separation is enhanced with an increase in the Knudsen number. This is because the resulting rarefaction effect reduces the interactions between the gas molecules of the two species, and thus increases the difference in the molecular velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phonon Boltzmann equation-based discrete unified gas kinetic scheme for multiscale heat transfer

Numerical prediction of multiscale heat transfer is a challenging problem due to the wide range of time and length scales involved. In this work a discrete unified gas kinetic scheme (DUGKS) is developed for heat transfer in materials with different acoustic thickness based on the phonon Boltzmann equation. With discrete phonon direction, the Boltzmann equation is discretized with a second-order finite-volume formulation, in which the time-step is fully determined by the Courant-Friedrichs-Lewy (CFL) condition. The scheme has the asymptotic preserving (AP) properties for both diffusive and ballistic regimes, and can present accurate solutions in the whole transition regime as well. The DUGKS is a self-adaptive multiscale method for the capturing of local transport process. Numerical tests for both heat transfers with different Knudsen numbers are presented to validate the current method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A general multiple-relaxation-time lattice Boltzmann model for nonlinear anisotropic convection-diffusion equations

In this paper, based on the previous work [B. Shi, Z. Guo, Lattice Boltzmann model for nonlinear convection-diffusion equations, Phys. Rev. E 79 (2009) 016701], we develop a general multiple-relaxation-time (MRT) lattice Boltzmann model for nonlinear anisotropic convection-diffusion equation (NACDE), and show that the NACDE can be recovered correctly from the present model through the Chapman-Enskog analysis. We then test the MRT model through some classic CDEs, and find that the numerical results are in good agreement with analytical solutions or some available results. Besides, the numerical results also show that similar to the single-relaxation-time (SRT) lattice Boltzmann model or so-called BGK model, the present MRT model also has a second-order convergence rate in space. Finally, we also perform a comparative study on the accuracy and stability of the MRT model and BGK model by using two examples. In terms of the accuracy, both the theoretical analysis and numerical results show that a \emph{numerical} slip on the boundary would be caused in the BGK model, and cannot be eliminated unless the relaxation parameter is fixed to be a special value, while the \emph{numerical} slip in the MRT model can be overcome once the relaxation parameters satisfy some constrains. The results in terms of stability also demonstrate that the MRT model could be more stable than the BGK model through tuning the free relaxation parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A lattice Boltzmann method for binary fluids based on mass-conserved quasi-incompressible phase-field theory

In this paper, a lattice Boltzmann equation (LBE) model is proposed for binary fluids based on a quasi-incompressible phase-field model [J. Shen et al, Comm. Comp. Phys. 13, 1045 (2013)]. Compared with the other incompressible LBE models based on the incompressible phase-field theory, the quasi-incompressible model conserves mass locally. A series of numerical simulations are performed to validate the proposed model, and comparisons with an incompressible LBE model [H. Liang et al, Phys. Rev. E 89, 053320 (2014)] are also carried out. It is shown that the proposed model can track the interface accurately, and the predictions by the quasi-incompressible and incompressible models agree qualitatively well as the distribution of chemical potential is uniform, otherwise differ significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A modified lattice Bhatnagar-Gross-Krook model for convection heat transfer in porous media

The lattice Bhatnagar-Gross-Krook (LBGK) model has become the most popular one in the lattice Boltzmann method for simulating the convection heat transfer in porous media. However, the LBGK model generally suffers from numerical instability at low fluid viscosities and effective thermal diffusivities. In this paper, a modified LBGK model is developed for incompressible thermal flows in porous media at the representative elementary volume scale, in which the shear rate and temperature gradient are incorporated into the equilibrium distribution functions. With two additional parameters, the relaxation times in the collision process can be fixed at a proper value invariable to the viscosity and the effective thermal diffusivity. In addition, by constructing a modified equilibrium distribution function and a source term in the evolution equation of temperature field, the present model can recover the macroscopic equations correctly through the Chapman-Enskog analysis, which is another key point different from previous LBGK models. Several benchmark problems are simulated to validate the present model with the proposed local computing scheme for the shear rate and temperature gradient, and the numerical results agree well with analytical solutions and/or those well-documented data in previous studies. It is also shown that the present model and the computational schemes for the gradient operators have a second-order accuracy in space, and better numerical stability of the present modified LBGK model than previous LBGK models is demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Discrete unified gas kinetic scheme on unstructured meshes

The recently proposed discrete unified gas kinetic scheme (DUGKS) is a finite volume method for deterministic solution of the Boltzmann model equation with asymptotic preserving property. In DUGKS, the numerical flux of the distribution function is determined from a local numerical solution of the Boltzmann model equation using an unsplitting approach. The time step and mesh resolution are not restricted by the molecular collision time and mean free path. To demonstrate the capacity of DUGKS in practical problems, this paper extends the DUGKS to arbitrary unstructured meshes. Several tests of both internal and external flows are performed, which include the cavity flow ranging from continuum to free molecular regimes, a multiscale flow between two connected cavities with a pressure ratio of 10000, and a high speed flow over a cylinder in slip and transitional regimes. The numerical results demonstrate the effectiveness of the DUGKS in simulating multiscale flow problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Performance evaluation of the general characteristics based off-lattice Boltzmann and DUGKS methods for low speed continuum flows: A comparative study

The general characteristics based off-lattice Boltzmann scheme (BKG) proposed by Bardow et~al.(2006), and the discrete unified gas kinetic scheme (DUGKS) are two methods that successfully overcome the time step restriction by the collision time, which is commonly seen in many other kinetic schemes. Basically, the BKG scheme is a time splitting scheme, while the DUGKS is an un-split finite volume scheme. In this work, we first perform a theoretical analysis of the two schemes in the finite volume framework by comparing their numerical flux evaluations. It is found that the effects of collision term are considered in the reconstructions of the cell-interface distribution function in both schemes, which explains why they can overcome the time step restriction and can give accurate results even as the time step is much larger than the collision time. The difference between the two schemes lies in the treatment of the integral of the collision term, in which the Bardow's scheme uses the rectangular rule while the DUGKS uses the trapezoidal rule. The performance of the two schemes, i.e., accuracy, stability, and efficiency are then compared by simulating several two dimensional flows, including the unsteady Taylor-Green vortex flow, the steady lid-driven cavity flow, and the laminar boundary layer problem. It is observed that, the DUGKS can give more accurate results than the BKG scheme. Furthermore, the numerical stability of the BKG scheme decreases as the Courant-Friedrichs-Lewy (CFL) number approaches to 1, while the stability of DUGKS is not affected by the CFL number apparently as long as CFL<1. It is also observed that the BKG scheme is about one time faster than the DUGKS scheme with the same computational mesh and time step.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A coupled discrete unified gas-kinetic scheme for Boussinesq flows

Recently, the discrete unified gas-kinetic scheme (DUGKS) [Z. L. Guo \emph{et al}., Phys. Rev. E ${\bf 88}$, 033305 (2013)] based on the Boltzmann equation is developed as a new multiscale kinetic method for isothermal flows. In this paper, a thermal and coupled discrete unified gas-kinetic scheme is derived for the Boussinesq flows, where the velocity and temperature fields are described independently. Kinetic boundary conditions for both velocity and temperature fields are also proposed. The proposed model is validated by simulating several canonical test cases, including the porous plate problem, the Rayleigh-bénard convection, and the natural convection with Rayleigh number up to $10^{10}$ in a square cavity. The results show that the coupled DUGKS is of second order accuracy in space and can well describe the convection phenomena from laminar to turbulent flows. Particularly, it is found that this new scheme has better numerical stability in simulating high Rayleigh number flows compared with the previous kinetic models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A semi-Lagrangian gas-kinetic scheme for smooth flows

In this paper, a semi-Lagrangian gas-kinetic scheme is developed for smooth flows based on the Bhatnagar-Gross-Krook (BGK) equation. As a finite-volume scheme, the evolution of the average flow variables in a control volume is under the Eulerian framework, whereas the construction of the numerical flux across the cell interface comes from the Lagrangian perspective. The adoption of the Lagrangian aspect makes the collision and the transport mechanisms intrinsically coupled together in the flux evaluation. As a result, the time step is independent of the particle collision time and solely determined by the Courant-Friedrichs-Lewy (CFL) conditions. A set of simulations are carried out to validate the performance of the new scheme. The results show that with second-order spatial accuracy, the scheme exhibits low numerical dissipation, and can accurately capture the Navier-Stokers solutions for the smooth flows with viscous heat dissipation from the low-speed incompressible to hypersonic compressible regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Discrete unified gas kinetic scheme for all Knudsen number flows: II. Compressible case

This paper is a continuation of our earlier work [Z.L. Guo {\it et al.}, Phys. Rev. E {\bf 88}, 033305 (2013)] where a multiscale numerical scheme based on kinetic model was developed for low speed isothermal flows with arbitrary Knudsen numbers. In this work, a discrete unified gas-kinetic scheme (DUGKS) for compressible flows with the consideration of heat transfer and shock discontinuity is developed based on the Shakhov model with an adjustable Prandtl number. The method is an explicit finite-volume scheme where the transport and collision processes are coupled in the evaluation of the fluxes at cell interfaces, so that the nice asymptotic preserving (AP) property is retained, such that the time step is limited only by the CFL number, the distribution function at cell interface recovers to the Chapman-Enskog one in the continuum limit while reduces to that of free-transport for free-molecular flow, and the time and spatial accuracy is of second-order accuracy in smooth region. These features make the DUGKS an ideal method for multiscale compressible flow simulations. A number of numerical tests, including the shock structure problem, the Sod tube problem with different degree of non-equilibrium, and the two-dimensional Riemann problem in continuum and rarefied regimes, are performed to validate the scheme. The comparisons with the results of DSMC and other benchmark data demonstrate that the DUGKS is a reliable and efficient method for multiscale compressible flow computation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Volume-averaged macroscopic equation for fluid flow in moving porous media

Darcy's law and the Brinkman equation are two main models used for creeping fluid flows inside moving permeable particles. For these two models, the time derivative and the nonlinear convective terms of fluid velocity are neglected in the momentum equation. In this paper, a new momentum equation including these two terms are rigorously derived from the pore-scale microscopic equations by the volume-averaging method, which can reduces to Darcy's law and the Brinkman equation under creeping flow conditions. Using the lattice Boltzmann equation method, the macroscopic equations are solved for the problem of a porous circular cylinder moving along the centerline of a channel. Galilean invariance of the equations are investigated both with the intrinsic phase averaged velocity and the phase averaged velocity. The results demonstrate that the commonly used phase averaged velocity cannot serve as the superficial velocity, while the intrinsic phase averaged velocity should be chosen for porous particulate systems.

0 citations0 reviewsNo code linkedOpen access