Modeling groundwater flow by lattice Boltzmann method in curvilinear coordinates

Abstract

In order to promote the use of the lattice Boltzmann method (LBM) for the simulation of isotropic groundwater°ow in a con¯ned aquifer with arbitrary geometry, Poisson's equation was transformed into a curvilinear coordinate system. With the metric function between the physical and the computational domain established, Poisson's equation written in Cartesian coordinates was transformed in curvilinear coordinates. Following, the appropriate equilibrium function for the D2Q9 square lattice has been de¯ned. The resulting curvilinear formulation of the LBM for groundwater°ow is capable of modeling°ow in domains of complex geometry with the opportunity of local re¯ning/coarsening of the computational mesh corresponding to the complexity of the°ow pattern and the required accuracy. Since the proposed form of the LBM uses the transformed equation of°ow implemented in the equilibrium function,¯nding a solution does not require supplementary procedures along the curvilinear boundaries, nor in the zones requiring mesh density adjustments. Thus, the basic concept of the LBM is completely maintained. The improvement of the proposed LBM over the previously published classical methods is completely veri¯ed by three examples with analytical solutions. The results demonstrate the advantages of the proposed curvilinear LBM in modeling groundwater°ow in complex°ow domains.