A high-order hybridizable discontinuous Galerkin method with fast convergence to steady-state solutions of the gas kinetic equation

Research output: Contribution to journalArticle

  • 1 Citations

Abstract

The mass flow rate of Poiseuille flow of rarefied gas through long ducts of two-dimensional cross-sections with arbitrary shape is critical in the pore-network modeling of gas transport in porous media. Here, for the first time, the high-order hybridizable discontinuous Galerkin (HDG) method is used to find the steady-state solution of the linearized Bhatnagar-Gross-Krook equation on two-dimensional triangular meshes. The velocity distribution function and its traces are approximated in piecewise polynomial spaces (of degree up to 4) on the triangular meshes and mesh skeletons, respectively. By employing a numerical flux that is derived from the first-order upwind scheme and imposing its continuity weakly on the mesh skeletons, global systems for unknown traces are obtained with fewer coupled degrees of freedom when compared to the original discontinuous Galerkin formulation. To achieve fast convergence to the steady-state solution, a diffusion-like equation for flow velocity, which is asymptotic-preserving into the fluid dynamic limit, is solved by the HDG simultaneously on the same meshes. The proposed HDG-synthetic iterative scheme is proved to be accurate and efficient. Specifically, for flows in the near-continuum regime, numerical simulations have shown that, to achieve the same level of accuracy, our scheme could be faster than the conventional iterative scheme by two orders of magnitude, also it is faster than the synthetic iterative scheme based on the finite difference discretization in the spatial space by one order of magnitude. In addition, the implicit HDG method is more efficient than an explicit discontinuous Galerkin gas kinetic solver, as well as the implicit discontinuous Galerkin scheme when the degree of approximating polynomial is larger than 2. The HDG-synthetic iterative scheme is ready to be extended to simulate rarefied gas mixtures and the Boltzmann collision operator.
LanguageEnglish
Pages973-991
Number of pages19
JournalJournal of Computational Physics
Volume376
Early online date1 Sep 2018
DOIs
StatePublished - 1 Jan 2019

Fingerprint

Kinetic theory of gases
Galerkin method
Galerkin methods
kinetic equations
mesh
Polynomials
Fluid dynamics
Velocity distribution
gases
Gas mixtures
Flow velocity
rarefied gases
Ducts
Distribution functions
Flow of gases
Porous materials
musculoskeletal system
Flow rate
Fluxes
polynomials

Keywords

  • hybridizable discontinuous Galerkin
  • gas kinetic equation
  • synthetic iterative scheme
  • fast convergence

Cite this

@article{643e0829095a40148d0d352e15dec460,
title = "A high-order hybridizable discontinuous Galerkin method with fast convergence to steady-state solutions of the gas kinetic equation",
abstract = "The mass flow rate of Poiseuille flow of rarefied gas through long ducts of two-dimensional cross-sections with arbitrary shape is critical in the pore-network modeling of gas transport in porous media. Here, for the first time, the high-order hybridizable discontinuous Galerkin (HDG) method is used to find the steady-state solution of the linearized Bhatnagar-Gross-Krook equation on two-dimensional triangular meshes. The velocity distribution function and its traces are approximated in piecewise polynomial spaces (of degree up to 4) on the triangular meshes and mesh skeletons, respectively. By employing a numerical flux that is derived from the first-order upwind scheme and imposing its continuity weakly on the mesh skeletons, global systems for unknown traces are obtained with fewer coupled degrees of freedom when compared to the original discontinuous Galerkin formulation. To achieve fast convergence to the steady-state solution, a diffusion-like equation for flow velocity, which is asymptotic-preserving into the fluid dynamic limit, is solved by the HDG simultaneously on the same meshes. The proposed HDG-synthetic iterative scheme is proved to be accurate and efficient. Specifically, for flows in the near-continuum regime, numerical simulations have shown that, to achieve the same level of accuracy, our scheme could be faster than the conventional iterative scheme by two orders of magnitude, also it is faster than the synthetic iterative scheme based on the finite difference discretization in the spatial space by one order of magnitude. In addition, the implicit HDG method is more efficient than an explicit discontinuous Galerkin gas kinetic solver, as well as the implicit discontinuous Galerkin scheme when the degree of approximating polynomial is larger than 2. The HDG-synthetic iterative scheme is ready to be extended to simulate rarefied gas mixtures and the Boltzmann collision operator.",
keywords = "hybridizable discontinuous Galerkin, gas kinetic equation, synthetic iterative scheme, fast convergence",
author = "Wei Su and Peng Wang and Yonghao Zhang and Lei Wu",
year = "2019",
month = "1",
day = "1",
doi = "10.1016/j.jcp.2018.08.050",
language = "English",
volume = "376",
pages = "973--991",
journal = "Journal of Computational Physics",
issn = "0021-9991",

}

TY - JOUR

T1 - A high-order hybridizable discontinuous Galerkin method with fast convergence to steady-state solutions of the gas kinetic equation

AU - Su,Wei

AU - Wang,Peng

AU - Zhang,Yonghao

AU - Wu,Lei

PY - 2019/1/1

Y1 - 2019/1/1

N2 - The mass flow rate of Poiseuille flow of rarefied gas through long ducts of two-dimensional cross-sections with arbitrary shape is critical in the pore-network modeling of gas transport in porous media. Here, for the first time, the high-order hybridizable discontinuous Galerkin (HDG) method is used to find the steady-state solution of the linearized Bhatnagar-Gross-Krook equation on two-dimensional triangular meshes. The velocity distribution function and its traces are approximated in piecewise polynomial spaces (of degree up to 4) on the triangular meshes and mesh skeletons, respectively. By employing a numerical flux that is derived from the first-order upwind scheme and imposing its continuity weakly on the mesh skeletons, global systems for unknown traces are obtained with fewer coupled degrees of freedom when compared to the original discontinuous Galerkin formulation. To achieve fast convergence to the steady-state solution, a diffusion-like equation for flow velocity, which is asymptotic-preserving into the fluid dynamic limit, is solved by the HDG simultaneously on the same meshes. The proposed HDG-synthetic iterative scheme is proved to be accurate and efficient. Specifically, for flows in the near-continuum regime, numerical simulations have shown that, to achieve the same level of accuracy, our scheme could be faster than the conventional iterative scheme by two orders of magnitude, also it is faster than the synthetic iterative scheme based on the finite difference discretization in the spatial space by one order of magnitude. In addition, the implicit HDG method is more efficient than an explicit discontinuous Galerkin gas kinetic solver, as well as the implicit discontinuous Galerkin scheme when the degree of approximating polynomial is larger than 2. The HDG-synthetic iterative scheme is ready to be extended to simulate rarefied gas mixtures and the Boltzmann collision operator.

AB - The mass flow rate of Poiseuille flow of rarefied gas through long ducts of two-dimensional cross-sections with arbitrary shape is critical in the pore-network modeling of gas transport in porous media. Here, for the first time, the high-order hybridizable discontinuous Galerkin (HDG) method is used to find the steady-state solution of the linearized Bhatnagar-Gross-Krook equation on two-dimensional triangular meshes. The velocity distribution function and its traces are approximated in piecewise polynomial spaces (of degree up to 4) on the triangular meshes and mesh skeletons, respectively. By employing a numerical flux that is derived from the first-order upwind scheme and imposing its continuity weakly on the mesh skeletons, global systems for unknown traces are obtained with fewer coupled degrees of freedom when compared to the original discontinuous Galerkin formulation. To achieve fast convergence to the steady-state solution, a diffusion-like equation for flow velocity, which is asymptotic-preserving into the fluid dynamic limit, is solved by the HDG simultaneously on the same meshes. The proposed HDG-synthetic iterative scheme is proved to be accurate and efficient. Specifically, for flows in the near-continuum regime, numerical simulations have shown that, to achieve the same level of accuracy, our scheme could be faster than the conventional iterative scheme by two orders of magnitude, also it is faster than the synthetic iterative scheme based on the finite difference discretization in the spatial space by one order of magnitude. In addition, the implicit HDG method is more efficient than an explicit discontinuous Galerkin gas kinetic solver, as well as the implicit discontinuous Galerkin scheme when the degree of approximating polynomial is larger than 2. The HDG-synthetic iterative scheme is ready to be extended to simulate rarefied gas mixtures and the Boltzmann collision operator.

KW - hybridizable discontinuous Galerkin

KW - gas kinetic equation

KW - synthetic iterative scheme

KW - fast convergence

UR - https://www.journals.elsevier.com/journal-of-computational-physics

U2 - 10.1016/j.jcp.2018.08.050

DO - 10.1016/j.jcp.2018.08.050

M3 - Article

VL - 376

SP - 973

EP - 991

JO - Journal of Computational Physics

T2 - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -