Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
10357141 | Journal of Computational Physics | 2009 | 23 Pages |
Abstract
We present a second order sharp interface finite volume method for the solution of the three-dimensional elliptic equation â·(β(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 β-ratios and non-smooth interfaces - the method achieves second order of accuracy in the Lâ and L2 norm.
Keywords
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science Applications
Authors
M. Oevermann, C. Scharfenberg, R. Klein,