Parallel modeling of three-dimensional variably saturated ground water flows with unstructured mesh using open source finite volume platform openfoam

Research output: Contribution to journalArticle

7 Citations (Scopus)

Abstract

We report the development of a three-dimensional open source model suGWFoam, based on the multiphysics platform OpenFOAM for variably saturated flow in porous media. The nonlinear Richards equation is spatially discretized on unstructured meshes using finite volume method in conjunction with the modified Picard iterative method, which inherently guarantees conservation of mass. A switching algorithm between different forms of the Richards equation is implemented and three choices of convergence criteria suitable for different simulation scenarios are provided. The linearized equation results in a symmetric matrix system which can be solved with one of the following Krylov subspace linear solvers: preconditioned conjugate gradient (PCG), incomplete Cholesky preconditioned conjugate gradient (ICCG) and geometric agglomerated algebraic multigrid (GAMG). Various preconditioners are runtime selectable. Parallel computation through domain decomposition shows good speedup for large scale simulations. Comparing to other available codes for ground water, the merits of this model reside in the architecture of OpenFOAM. Built upon the modularized and hierarchical design of the platform, this model has less than several hundreds of lines of code on its top level and yet has the full functionalities for three-dimensional unstructured mesh, parallel computing, and numerous choices of spatial and temporal discretization schemes. This liberates the researchers in the water resource community from the tedious and long process of coding. It is also extremely easy to couple groundwater flow with surface water hydrodynamics using existing flow solvers (such as channel flows and waves) in OpenFOAM. This model is fully functional and validations show excellent agreement with data in literature for a wide range of applications. It can be used for both fundamental researches and real world applications.

Original languageEnglish (US)
Pages (from-to)223-238
Number of pages16
JournalEngineering Applications of Computational Fluid Mechanics
Volume7
Issue number2
DOIs
StatePublished - Jan 1 2013

Fingerprint

Groundwater Flow
Groundwater flow
Unstructured Mesh
Finite Volume
Open Source
Richards Equation
Preconditioned Conjugate Gradient
Three-dimensional
Modeling
Algebraic multigrid
Cholesky
Convergence Criteria
Krylov Subspace
Flow in Porous Media
Multiphysics
Water Resources
Channel Flow
Discretization Scheme
Ground Water
Finite volume method

All Science Journal Classification (ASJC) codes

  • Computer Science(all)
  • Modeling and Simulation

Cite this

@article{8aa75898ceed4d71a17b7e3234a7a49f,
title = "Parallel modeling of three-dimensional variably saturated ground water flows with unstructured mesh using open source finite volume platform openfoam",
abstract = "We report the development of a three-dimensional open source model suGWFoam, based on the multiphysics platform OpenFOAM for variably saturated flow in porous media. The nonlinear Richards equation is spatially discretized on unstructured meshes using finite volume method in conjunction with the modified Picard iterative method, which inherently guarantees conservation of mass. A switching algorithm between different forms of the Richards equation is implemented and three choices of convergence criteria suitable for different simulation scenarios are provided. The linearized equation results in a symmetric matrix system which can be solved with one of the following Krylov subspace linear solvers: preconditioned conjugate gradient (PCG), incomplete Cholesky preconditioned conjugate gradient (ICCG) and geometric agglomerated algebraic multigrid (GAMG). Various preconditioners are runtime selectable. Parallel computation through domain decomposition shows good speedup for large scale simulations. Comparing to other available codes for ground water, the merits of this model reside in the architecture of OpenFOAM. Built upon the modularized and hierarchical design of the platform, this model has less than several hundreds of lines of code on its top level and yet has the full functionalities for three-dimensional unstructured mesh, parallel computing, and numerous choices of spatial and temporal discretization schemes. This liberates the researchers in the water resource community from the tedious and long process of coding. It is also extremely easy to couple groundwater flow with surface water hydrodynamics using existing flow solvers (such as channel flows and waves) in OpenFOAM. This model is fully functional and validations show excellent agreement with data in literature for a wide range of applications. It can be used for both fundamental researches and real world applications.",
author = "Xiaofeng Liu",
year = "2013",
month = "1",
day = "1",
doi = "10.1080/19942060.2013.11015466",
language = "English (US)",
volume = "7",
pages = "223--238",
journal = "Engineering Applications of Computational Fluid Mechanics",
issn = "1994-2060",
publisher = "Department of Civil and Environmental Engineering, The Hong Kong Polytechnic University",
number = "2",

}

TY - JOUR

T1 - Parallel modeling of three-dimensional variably saturated ground water flows with unstructured mesh using open source finite volume platform openfoam

AU - Liu, Xiaofeng

PY - 2013/1/1

Y1 - 2013/1/1

N2 - We report the development of a three-dimensional open source model suGWFoam, based on the multiphysics platform OpenFOAM for variably saturated flow in porous media. The nonlinear Richards equation is spatially discretized on unstructured meshes using finite volume method in conjunction with the modified Picard iterative method, which inherently guarantees conservation of mass. A switching algorithm between different forms of the Richards equation is implemented and three choices of convergence criteria suitable for different simulation scenarios are provided. The linearized equation results in a symmetric matrix system which can be solved with one of the following Krylov subspace linear solvers: preconditioned conjugate gradient (PCG), incomplete Cholesky preconditioned conjugate gradient (ICCG) and geometric agglomerated algebraic multigrid (GAMG). Various preconditioners are runtime selectable. Parallel computation through domain decomposition shows good speedup for large scale simulations. Comparing to other available codes for ground water, the merits of this model reside in the architecture of OpenFOAM. Built upon the modularized and hierarchical design of the platform, this model has less than several hundreds of lines of code on its top level and yet has the full functionalities for three-dimensional unstructured mesh, parallel computing, and numerous choices of spatial and temporal discretization schemes. This liberates the researchers in the water resource community from the tedious and long process of coding. It is also extremely easy to couple groundwater flow with surface water hydrodynamics using existing flow solvers (such as channel flows and waves) in OpenFOAM. This model is fully functional and validations show excellent agreement with data in literature for a wide range of applications. It can be used for both fundamental researches and real world applications.

AB - We report the development of a three-dimensional open source model suGWFoam, based on the multiphysics platform OpenFOAM for variably saturated flow in porous media. The nonlinear Richards equation is spatially discretized on unstructured meshes using finite volume method in conjunction with the modified Picard iterative method, which inherently guarantees conservation of mass. A switching algorithm between different forms of the Richards equation is implemented and three choices of convergence criteria suitable for different simulation scenarios are provided. The linearized equation results in a symmetric matrix system which can be solved with one of the following Krylov subspace linear solvers: preconditioned conjugate gradient (PCG), incomplete Cholesky preconditioned conjugate gradient (ICCG) and geometric agglomerated algebraic multigrid (GAMG). Various preconditioners are runtime selectable. Parallel computation through domain decomposition shows good speedup for large scale simulations. Comparing to other available codes for ground water, the merits of this model reside in the architecture of OpenFOAM. Built upon the modularized and hierarchical design of the platform, this model has less than several hundreds of lines of code on its top level and yet has the full functionalities for three-dimensional unstructured mesh, parallel computing, and numerous choices of spatial and temporal discretization schemes. This liberates the researchers in the water resource community from the tedious and long process of coding. It is also extremely easy to couple groundwater flow with surface water hydrodynamics using existing flow solvers (such as channel flows and waves) in OpenFOAM. This model is fully functional and validations show excellent agreement with data in literature for a wide range of applications. It can be used for both fundamental researches and real world applications.

UR - http://www.scopus.com/inward/record.url?scp=84880160312&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84880160312&partnerID=8YFLogxK

U2 - 10.1080/19942060.2013.11015466

DO - 10.1080/19942060.2013.11015466

M3 - Article

AN - SCOPUS:84880160312

VL - 7

SP - 223

EP - 238

JO - Engineering Applications of Computational Fluid Mechanics

JF - Engineering Applications of Computational Fluid Mechanics

SN - 1994-2060

IS - 2

ER -