Source author record

Lianhua Zhu

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

6 published item(s)

preprint2020arXiv

Fast convergence and asymptotic preserving of the General Synthetic Iterative Scheme

Recently the general synthetic iteration scheme (GSIS) is proposed to find the steady-state solution of the Boltzmann equation~\cite{SuArXiv2019}, where various numerical simulations have shown that (i) the steady-state solution can be found within dozens of iterations at any Knudsen number $K$, and (ii) the solution is accurate even when the spatial cell size in the bulk region is much larger than the molecular mean free path, i.e. Navier-Stokes solutions are recovered at coarse grids. The first property indicates that the error decay rate between two consecutive iterations decreases to zero with $K$, while the second one implies that the GSIS is asymptotically preserving the Navier-Stokes limit. This paper is dedicated to the rigorous proof of both properties.

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.

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.

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.

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.