The rate of convergence of conjugate gradients
Numerische Mathematik
SIAM Journal on Numerical Analysis
Balancing domain decomposition for problems with large jumps in coefficients
Mathematics of Computation
A boundary condition capturing method for Poisson's equation on irregular domains
Journal of Computational Physics
Dual-Primal FETI Methods for Three-Dimensional Elliptic Problems with Heterogeneous Coefficients
SIAM Journal on Numerical Analysis
A numerical method for solving variable coefficient elliptic equation with interfaces
Journal of Computational Physics
Analysis of a Stokes interface problem
Numerische Mathematik
Journal of Computational Physics
Journal of Computational Physics
A sharp interface finite volume method for elliptic equations on Cartesian grids
Journal of Computational Physics
Journal of Computational and Applied Mathematics
Hi-index | 31.45 |
When the coefficients of an elliptic problem have jumps of several orders of magnitude across an embedded interface, many iterative solvers exhibit deteriorated convergence properties or a loss of efficiency and it is difficult to achieve high solution accuracies in the whole domain. In this paper we present an asymptotic solution approach for the elliptic problem @?@?(@b(x)@?u(x))=f(x) on a domain @W=@W^+@?@W^- with piecewise constant coefficients @b^+, @b^- with @b^+@?@b^- and prescribed jump conditions at an embedded interface @C separating the domains @W^+ and @W^-. We are in particular focusing on a problem related to fluid mechanics, namely incompressible two-phase flow with a large density ratio across the phase boundary, where an accurate solution of the velocity depends on the accurate solution of a pressure Poisson equation with equal local relative errors in the whole domain. Instead of solving the equation in a single solution step we decompose the problem into two consecutive problems based on an asymptotic analysis of the physical problem where each problem is asymptotically independent of the ratio of coefficients @e=@b^-/@b^+. The proposed methods lead to a robust and accurate solution of the elliptic problem using standard black-box iterative solvers.