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

Yonghao Zhang

Yonghao Zhang 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 Sciencephysics.geo-ph

Catalog footprint

What is connected

11works
4topics
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

11 published item(s)

preprint2020arXiv

General synthetic iteration scheme for non-linear gas kinetic simulation of multi-scale rarefied gas flows

The general synthetic iteration scheme (GSIS) is extended to find the steady-state solution of nonlinear gas kinetic equation, removing the long-standing problems of slow convergence and requirement of ultra-fine grids in near-continuum flows. The key ingredients of GSIS are that the gas kinetic equation and macroscopic synthetic equations are tightly coupled, and the constitutive relations in macroscopic synthetic equations explicitly contain Newton's law of shear stress and Fourier's law of heat conduction. The higher-order constitutive relations describing rarefaction effects are calculated from the velocity distribution function, however, their constructions are simpler than our previous work (Su et al. Journal of Computational Physics 407 (2020) 109245) for linearized gas kinetic equations. On the other hand, solutions of macroscopic synthetic equations are used to inform the evolution of gas kinetic equation at the next iteration step. A rigorous linear Fourier stability analysis in periodic system shows that the error decay rate of GSIS can be smaller than 0.5, which means that the deviation to steady-state solution can be reduced by 3 orders of magnitude in 10 iterations. Other important advantages of the GSIS are (i) it does not rely on the specific form of Boltzmann collision operator and (ii) it can be solved by sophisticated techniques in computational fluid dynamics, making it amenable to large scale engineering applications. In this paper, the efficiency and accuracy of GSIS is demonstrated by a number of canonical test cases in rarefied gas dynamics.

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
preprint2019arXiv

Can we find steady-state solutions to multiscale rarefied gas flows within dozens of iterations?

One of the central problems in the study of rarefied gas dynamics is to find the steady-state solution of the Boltzmann equation quickly. When the Knudsen number is large, i.e. the system is highly rarefied, the conventional iteration scheme can lead to convergence within a few iterations. However, when the Knudsen number is small, i.e. the flow falls in the near-continuum regime, hundreds of thousands iterations are needed, and yet the "converged" solutions are prone to be contaminated by accumulated error and large numerical dissipation. Recently, based on the gas kinetic models, the implicit unified gas kinetic scheme (UGKS) and its variants have significantly reduced the iterations in the near-continuum flow regime, but still much higher than that of the highly rarefied gas flows. In this paper, we put forward a general synthetic iteration scheme (GSIS) to find the steady-state solutions of general rarefied gas flows within dozens of iterations at any Knudsen number. As the GSIS does not rely on the specific kinetic model/collision operator, it can be naturally extended to quickly find converged solutions for mixture flows and even flows involving chemical reactions. These two superior advantages are also expected to accelerate the slow convergence in simulation of near-continuum flows via the direct simulation Monte Carlo method and its low-variance version.

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
preprint2018arXiv

High-Order Implicit Hybridizable Discontinuous Galerkin Method for the Boltzmann Equation

The high-order hybridizable discontinuous Galerkin (HDG) method combining with an implicit iterative scheme is used to find the steady-state solution of the Boltzmann equation with full collision integral on two-dimensional triangular meshes. The velocity distribution function and its trace are approximated in the piecewise polynomial space of degree up to 4. The fast spectral method (FSM) is incorporated into the DG discretization to evaluate the collision operator. Specific polynomial approximation is proposed for the collision term to reduce the computational cost. The proposed scheme is proved to be accurate and efficient.

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

Lattice Boltzmann modelling of intrinsic permeability

Lattice Boltzmann method (LBM) has been applied to predict flow properties of porous media including intrinsic permeability, where it is implicitly assumed that the LBM is equivalent to the incompressible (or near incompressible) Navier-Stokes equation. However, in LBM simulations, high-order moments, which are completely neglected in the Navier-Stokes equation, are still available through particle distribution functions. To ensure that the LBM simulation is correctly working at the Navier-Stokes hydrodynamic level, the high-order moments have to be negligible. This requires that the Knudsen number (Kn) is small so that rarefaction effect can be ignored. In this technical note, we elaborate this issue in LBM modelling of porous media flows, which is particularly important for gas flows in ultra-tight media.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Beyond the Knudsen number: assessing thermodynamic non-equilibrium in gas flows

For more than 150 years the Navier-Stokes equations for thermodynamically quasi-equilibrium flows have been the cornerstone of modern computational fluid dynamics that underpins new fluid technologies. However, the applicable regime of the Navier-Stokes model in terms of the level of thermodynamic non-equilibrium in the local flowfield is not clear especially for hypersonic and low-speed micro/nano flows. Here, we re-visit the Navier-Stokes model in the framework of Boltzmann statistics, and propose a new and more appropriate way of assessing non-equilibrium in the local flowfield, and the corresponding appropriateness of the Navier-Stokes model. Our theoretical analysis and numerical simulations confirm our proposed method. Through molecular dynamics simulations we reveal that the commonly-used Knudsen number, or a parametric combination of Knudsen and Mach numbers, may not be sufficient to accurately assess the departure of flowfields from equilibrium, and the applicability of the Navier-Stokes model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kinetic diffuse boundary condition for high-order lattice Boltzmann model with streaming-collision mechanism

The implementation of the kinetic diffuse boundary condition with the characteristic streaming-collision mechanism is studied for the high-order lattice Boltzmann (LB) models. The obtained formulation is also tested and validated numerically for three high-order LB models for both isothermal and thermal Couette flows. The streaming-collision mechanism ensures that high-order LB models can retain particle feature while go beyond the Navier-Stokes hydrodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lattice ellipsoidal statistical BGK model for thermal non-equilibrium flows

A thermal lattice Boltzmann model is constructed on the basis of the ellipsoidal statistical Bhatnagar-Gross-Krook (ES-BGK) collision operator via the Hermite moment representation. The resulting lattice ES-BGK model uses a single distribution function and features an adjustable Prandtl number. Numerical simulations show that using a moderate discrete velocity set, this model can accurately recover steady and transient solutions of the ES-BGK equation in the slip-flow and early transition regimes in the small Mach number limit that is typical of microscale problems of practical interest. In the transition regime in particular, comparisons with numerical solutions of the ES-BGK model, direct Monte Carlo and low-variance deviational Monte Carlo simulations show good accuracy for values of the Knudsen number up to approximately 0.5. On the other hand, highly non-equilibrium phenomena characterized by high Mach numbers, such as viscous heating and force-driven Poiseuille flow for large values of the driving force, are more difficult to capture quantitatively in the transition regime using discretizations chosen with computational efficiency in mind such as the one used here, although improved accuracy is observed as the number of discrete velocities is increased.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multiscale lattice Boltzmann approach to modeling gas flows

For multiscale gas flows, kinetic-continuum hybrid method is usually used to balance the computational accuracy and efficiency. However, the kinetic-continuum coupling is not straightforward since the coupled methods are based on different theoretical frameworks. In particular, it is not easy to recover the non-equilibrium information required by the kinetic method which is lost by the continuum model at the coupling interface. Therefore, we present a multiscale lattice Boltzmann (LB) method which deploys high-order LB models in highly rarefied flow regions and low-order ones in less rarefied regions. Since this multiscale approach is based on the same theoretical framework, the coupling precess becomes simple. The non-equilibrium information will not be lost at the interface as low-order LB models can also retain this information. The simulation results confirm that the present method can achieve model accuracy with reduced computational cost.

0 citations0 reviewsNo code linkedOpen access