A sharp interface finite volume method for elliptic equations on Cartesian grids

  • Authors:

  • M. Oevermann;C. Scharfenberg;R. Klein

  • Affiliations:

  • Technische Universität Berlin, Institut für Energietechnik, Fasanenstr. 89, 10623 Berlin, Germany;Freie Universität Berlin, Fachbereich Mathematik und Informatik, Arnimallee 14, 14195 Berlin, Germany;Freie Universität Berlin, Fachbereich Mathematik und Informatik, Arnimallee 14, 14195 Berlin, Germany

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2009

Quantified Score

Hi-index 31.50

Visualization

Abstract

We present a second order sharp interface finite volume method for the solution of the three-dimensional elliptic equation @?.(@b(x-)@?u(x-))=f(x-) with variable coefficients on Cartesian grids. In particular, we focus on interface problems with discontinuities in the coefficient, the source term, the solution, and the fluxes across the interface. The method uses standard piecewise trilinear finite elements for normal cells and a double piecewise trilinear ansatz for the solution on cells intersected by the interface resulting always in a compact 27-point stencil. Singularities associated with vanishing partial volumes of intersected grid cells are removed by a two-term asymptotic approach. In contrast to the 2D method presented by two of the authors in [M. Oevermann, R. Klein, A Cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces, Journal of Computational Physics 219 (2006) 749-769] we use a minimization technique to determine the unknown coefficients of the double trilinear ansatz. This simplifies the treatment of the different cut-cell types and avoids additional special operations for degenerated interface topologies. The resulting set of linear equations has been solved with a BiCGSTAB solver preconditioned with an algebraic multigrid. In various testcases - including large @b-ratios and non-smooth interfaces - the method achieves second order of accuracy in the L"~ and L"2 norm.