Domain-decomposition method for parallel lattice Boltzmann simulation of incompressible flow in porous media

Junye Wang, Xiaoxian Zhang, A. Glyn Bengough, John W. Crawford

Research output: Contribution to journalArticle

  • 42 Citations

Abstract

The lattice Boltzmann method has proven to be a promising method to simulate flow in porous media. Its practical application often relies on parallel computation because of the demand for a large domain and fine grid resolution to adequately resolve pore heterogeneity. The existing domain-decomposition methods for parallel computation usually decompose a domain into a number of subdomains first and then recover the interfaces and perform the load balance. Normally, the interface recovery and the load balance have to be performed iteratively until an acceptable load balance is achieved; this costs time. In this paper we propose a cell-based domain-decomposition method for parallel lattice Boltzmann simulation of flow in porous media. Unlike the existing methods, the cell-based method performs the load balance first to divide the total number of fluid cells into a number of groups (or subdomains), in which the difference of fluid cells in each group is either 0 or 1, depending on if the total number of fluid cells is a multiple of the processor numbers; the interfaces between the subdomains are recovered at last. The cell-based method is to recover the interfaces rather than the load balance; it does not need iteration and gives an exact load balance. The performance of the proposed method is analyzed and compared using different computer systems; the results indicate that it reaches the theoretical parallel efficiency. The method is then applied to simulate flow in a three-dimensional porous medium obtained with microfocus x-ray computed tomography to calculate the permeability, and the result shows good agreement with the experimental data.
Original languageEnglish
Article number016706
JournalPhysical Review E
Volume72
Issue number1
DOIs
StatePublished - Jul 2005

Fingerprint

Porous materials
Domain decomposition methods
Fluids
Incompressible flow
Tomography
Computer systems
Recovery
X rays
Costs

Cite this

Wang, Junye; Zhang, Xiaoxian; Bengough, A. Glyn; Crawford, John W. / Domain-decomposition method for parallel lattice Boltzmann simulation of incompressible flow in porous media.

In: Physical Review E, Vol. 72, No. 1, 016706, 07.2005.

Research output: Contribution to journalArticle

@article{c5d963bf2e3945edab2e53581a96f6c5,
title = "Domain-decomposition method for parallel lattice Boltzmann simulation of incompressible flow in porous media",
abstract = "The lattice Boltzmann method has proven to be a promising method to simulate flow in porous media. Its practical application often relies on parallel computation because of the demand for a large domain and fine grid resolution to adequately resolve pore heterogeneity. The existing domain-decomposition methods for parallel computation usually decompose a domain into a number of subdomains first and then recover the interfaces and perform the load balance. Normally, the interface recovery and the load balance have to be performed iteratively until an acceptable load balance is achieved; this costs time. In this paper we propose a cell-based domain-decomposition method for parallel lattice Boltzmann simulation of flow in porous media. Unlike the existing methods, the cell-based method performs the load balance first to divide the total number of fluid cells into a number of groups (or subdomains), in which the difference of fluid cells in each group is either 0 or 1, depending on if the total number of fluid cells is a multiple of the processor numbers; the interfaces between the subdomains are recovered at last. The cell-based method is to recover the interfaces rather than the load balance; it does not need iteration and gives an exact load balance. The performance of the proposed method is analyzed and compared using different computer systems; the results indicate that it reaches the theoretical parallel efficiency. The method is then applied to simulate flow in a three-dimensional porous medium obtained with microfocus x-ray computed tomography to calculate the permeability, and the result shows good agreement with the experimental data.",
author = "Junye Wang and Xiaoxian Zhang and Bengough, {A. Glyn} and Crawford, {John W.}",
year = "2005",
month = "7",
doi = "10.1103/PhysRevE.72.016706",
volume = "72",
journal = "Physical Review E",
number = "1",

}

Domain-decomposition method for parallel lattice Boltzmann simulation of incompressible flow in porous media. / Wang, Junye; Zhang, Xiaoxian; Bengough, A. Glyn; Crawford, John W.

In: Physical Review E, Vol. 72, No. 1, 016706, 07.2005.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Domain-decomposition method for parallel lattice Boltzmann simulation of incompressible flow in porous media

AU - Wang,Junye

AU - Zhang,Xiaoxian

AU - Bengough,A. Glyn

AU - Crawford,John W.

PY - 2005/7

Y1 - 2005/7

N2 - The lattice Boltzmann method has proven to be a promising method to simulate flow in porous media. Its practical application often relies on parallel computation because of the demand for a large domain and fine grid resolution to adequately resolve pore heterogeneity. The existing domain-decomposition methods for parallel computation usually decompose a domain into a number of subdomains first and then recover the interfaces and perform the load balance. Normally, the interface recovery and the load balance have to be performed iteratively until an acceptable load balance is achieved; this costs time. In this paper we propose a cell-based domain-decomposition method for parallel lattice Boltzmann simulation of flow in porous media. Unlike the existing methods, the cell-based method performs the load balance first to divide the total number of fluid cells into a number of groups (or subdomains), in which the difference of fluid cells in each group is either 0 or 1, depending on if the total number of fluid cells is a multiple of the processor numbers; the interfaces between the subdomains are recovered at last. The cell-based method is to recover the interfaces rather than the load balance; it does not need iteration and gives an exact load balance. The performance of the proposed method is analyzed and compared using different computer systems; the results indicate that it reaches the theoretical parallel efficiency. The method is then applied to simulate flow in a three-dimensional porous medium obtained with microfocus x-ray computed tomography to calculate the permeability, and the result shows good agreement with the experimental data.

AB - The lattice Boltzmann method has proven to be a promising method to simulate flow in porous media. Its practical application often relies on parallel computation because of the demand for a large domain and fine grid resolution to adequately resolve pore heterogeneity. The existing domain-decomposition methods for parallel computation usually decompose a domain into a number of subdomains first and then recover the interfaces and perform the load balance. Normally, the interface recovery and the load balance have to be performed iteratively until an acceptable load balance is achieved; this costs time. In this paper we propose a cell-based domain-decomposition method for parallel lattice Boltzmann simulation of flow in porous media. Unlike the existing methods, the cell-based method performs the load balance first to divide the total number of fluid cells into a number of groups (or subdomains), in which the difference of fluid cells in each group is either 0 or 1, depending on if the total number of fluid cells is a multiple of the processor numbers; the interfaces between the subdomains are recovered at last. The cell-based method is to recover the interfaces rather than the load balance; it does not need iteration and gives an exact load balance. The performance of the proposed method is analyzed and compared using different computer systems; the results indicate that it reaches the theoretical parallel efficiency. The method is then applied to simulate flow in a three-dimensional porous medium obtained with microfocus x-ray computed tomography to calculate the permeability, and the result shows good agreement with the experimental data.

U2 - 10.1103/PhysRevE.72.016706

DO - 10.1103/PhysRevE.72.016706

M3 - Article

VL - 72

JO - Physical Review E

T2 - Physical Review E

JF - Physical Review E

IS - 1

M1 - 016706

ER -