A higher-order finite-volume discretization method for Poisson's equation in cut cell geometries

We present a method for generating higher-order finite volume discretizations for Poisson's equation on Cartesian cut cell grids in two and three dimensions. The discretization is in flux-divergence form, and stencils for the flux are computed by solving small weighted least-squares linear systems. Weights are the key in generating a stable discretization. We apply the method to solve Poisson's equation on a variety of geometries, and we demonstrate that the method can achieve second and fourth order accuracy in both truncation and solution error for these examples. We also show that the Laplacian operator has only stable eigenvalues for each of these examples.

[1]  Michael J. Aftosmis,et al.  Robust and efficient Cartesian mesh generation for component-based geometry , 1997 .

[2]  L. Greengard,et al.  A Fast Poisson Solver for Complex Geometries , 1995 .

[3]  Mark F. Adams,et al.  Algebraic multigrid methods for constrained linear systems with applications to contact problems in solid mechanics , 2003, Numer. Linear Algebra Appl..

[4]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[5]  Eloy Romero,et al.  A parallel implementation of Davidson methods for large-scale eigenvalue problems in SLEPc , 2014, TOMS.

[6]  Phillip Colella,et al.  HIGH-ACCURACY EMBEDDED BOUNDARY GRID GENERATION USING THE DIVERGENCE THEOREM , 2015 .

[7]  Phillip Colella,et al.  A cartesian grid embedded boundary method for the heat equation and poisson's equation in three dimensions , 2004 .

[8]  Vicente Hernández,et al.  SLEPc: Scalable Library for Eigenvalue Problem Computations , 2002, VECPAR.

[9]  William Gropp,et al.  Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries , 1997, SciTools.

[10]  L. Greengard,et al.  Regular Article: A Fast Adaptive Multipole Algorithm in Three Dimensions , 1999 .

[11]  S. Pirzadeh.  Advanced Unstructured Grid Generation for Complex Aerodynamic Applications , 2013 .

[12]  R. Fedkiw,et al.  A fourth order accurate discretization for the Laplace and heat equations on arbitrary domains, with applications to the Stefan problem , 2005 .

[13]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[14]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[15]  L. Greengard,et al.  A Direct Adaptive Poisson Solver of Arbitrary Order Accuracy , 1996 .

[16]  P. Colella,et al.  A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains , 1998 .