Search

A high-performance gas kinetic solver using multi-level parallelization is developed to enable pore-scale simulations of rarefied flows in porous media. The Boltzmann model equation is solved by the discrete velocity method with an iterative scheme. The multi-level MPI/OpenMP parallelization is implemented with the aim to efficiently utilise the computational resources to allow direct simulation of rarefied gas flows in porous media based on digital rock images for the first time. The multi-level parallel approach is analyzed in details confirming its better performance than the commonly-used MPI processing alone for an iterative scheme. With high communication efficiency and appropriate load balancing among CPU processes, parallel efficiency of 94% is achieved for 1536 cores in the 2D simulations, and 81% for 12288 cores in the 3D simulations. While decomposition in the spatial space does not affect the simulation results, one additional benefit of this approach is that the number of subdomains can be kept minimal to avoid deterioration of the convergence rate of the iteration process. This multi-level parallel approach can be readily extended to solve other Boltzmann model equations., Comment: 23 pages, 9 figures; Submitted to Computer Physics Communication; Uploaded to ResearchGate

In a recent paper by Lasseux, Vald\'{e}s-Parada and Porter (J.~Fluid~Mech. \textbf{805} (2016) 118-146), it is found that the apparent gas permeability of the porous media is a nonlinear function of the Knudsen number. However, this result is highly questionable, because the adopted Navier-Stokes equations and the first-order velocity-slip boundary condition are first-order (in terms of the Knudsen number) approximations of the Boltzmann equation and the kinetic boundary condition for rarefied gas flows. Our numerical simulations based on the Bhatnagar-Gross-Krook kinetic equation and regularized 20-moment equations prove that the Navier-Stokes equations with the first-order velocity-slip boundary condition are only accurate at a very small Knudsen number limit, where the apparent gas permeability is a linear function of the Knudsen number., Comment: a comment paper regarding the application range of the Navier-Stokes equation with the first-order velocity-slip boundary condition in porous media

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., Comment: 25 pages,31 figures

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.