Author ' s personal copy A Cartesian grid embedded boundary method for solving the Poisson and heat equations with discontinuous coefficients in three dimensions

We present a method for solving Poisson and heat equations with discon- tinuous coefficients in two- and three-dimensions. It uses a Cartesian cut-cell/embedded boundary method to represent the interface between materi- als, as described in Johansen & Colella (1998). Matching conditions across the interface are enforced using an approximation to fluxes at the boundary. Overall second order accuracy is achieved, as indicated by an array of tests using non-trivial interface geometries. Both the elliptic and heat solvers are shown to remain stable and efficient for material coefficient contrasts up to 106, thanks in part to the use of geometric multigrid. A test of accuracy when adaptive mesh refinement capabilities are utilized is also performed. An example problem relevant to nuclear reactor core simulation is presented, demonstrating the ability of the method to solve problems with realistic physical parameters.

[1]  Zhilin Li A Fast Iterative Algorithm for Elliptic Interface Problems , 1998 .

[2]  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 .

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

[4]  Kazufumi Ito,et al.  Three-dimensional elliptic solvers for interface problems and applications , 2003 .

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

[6]  Tianbing Chen,et al.  Piecewise-polynomial discretization and Krylov-accelerated multigrid for elliptic interface problems , 2008, J. Comput. Phys..

[7]  Kazufumi Ito,et al.  Maximum Principle Preserving Schemes for Interface Problems with Discontinuous Coefficients , 2001, SIAM J. Sci. Comput..

[8]  Michael Oevermann,et al.  A Cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces , 2006, J. Comput. Phys..

[9]  Leslie Greengard,et al.  Electrostatics and heat conduction in high contrast composite materials , 2006 .

[10]  Franck Nicoud,et al.  Conservative High-Order Finite-Difference Schemes for Low-Mach Number Flows , 2000 .

[11]  William D. Henshaw,et al.  A composite grid solver for conjugate heat transfer in fluid-structure systems , 2009, J. Comput. Phys..

[12]  Xu-dong Liu,et al.  A numerical method for solving variable coefficient elliptic equation with interfaces , 2005 .

[13]  P. Colella,et al.  Embedded boundary grid generation using the divergence theorem, implicit functions, and constructive solid geometry , 2008 .

[14]  Michael Oevermann,et al.  A sharp interface finite volume method for elliptic equations on Cartesian grids , 2009, J. Comput. Phys..

[15]  C. Peskin Numerical analysis of blood flow in the heart , 1977 .

[16]  John B. Bell,et al.  An Adaptive Projection Method for Unsteady, Low-Mach Number Combustion , 1998 .

[17]  Anita Mayo,et al.  The rapid evaluation of volume integrals of potential theory on general regions , 1992 .

[18]  Peter McCorquodale,et al.  A Cartesian grid embedded boundary method for the heat equation on irregular domains , 2001 .

[19]  Shan Zhao,et al.  High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources , 2006, J. Comput. Phys..

[20]  Ray A. Berry,et al.  A CFD M&S PROCESS FOR FAST REACTOR FUEL ASSEMBLIES , 2008 .

[21]  M. Berger,et al.  Robust and efficient Cartesian mesh generation for component-based geometry , 1998 .

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

[23]  P. P. Starling The numerical solution of Laplace's equation , 1963 .

[24]  Li-Tien Cheng,et al.  A second-order-accurate symmetric discretization of the Poisson equation on irregular domains , 2002 .

[25]  L. Leibowitz,et al.  Thermodynamic and transport properties of sodium liquid and vapor , 1995 .

[26]  Edward H. Twizell,et al.  Second-order,L0-stable methods for the heat equation with time-dependent boundary conditions , 1996, Adv. Comput. Math..

[27]  R. Fedkiw,et al.  A Boundary Condition Capturing Method for Poisson's Equation on Irregular Domains , 2000 .

[28]  M S Day,et al.  Numerical simulation of laminar reacting flows with complex chemistry , 2000 .

[29]  David Trebotich,et al.  Performance of embedded boundary methods for CFD with complex geometry , 2008 .