Researcher profile

Yonghao Zhang

Yonghao Zhang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

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

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.

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.

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.

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.

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.

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.