An asymptotic solution approach for elliptic equations with discontinuous coefficients
Journal article, 2014

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 ∇; {dot operator} (β(x)∇;u(x)) = f(x) on a domain Ω = Ω + ∪ Ω - with piecewise constant coefficients β + , β - with β + ≫ β - and prescribed jump conditions at an embedded interface Γ separating the domains Ω + and Ω - . 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 ε = β - / β + . The proposed methods lead to a robust and accurate solution of the elliptic problem using standard black-box iterative solvers. © 2013 Elsevier Inc.

Author

Michael Oevermann

Chalmers, Applied Mechanics, Combustion and Propulsion Systems

Rupert Klein

Freie Universität Berlin

Journal of Computational Physics

0021-9991 (ISSN) 1090-2716 (eISSN)

Vol. 261 230-243

Areas of Advance

Transport

Energy

Subject Categories

Computational Mathematics

Fluid Mechanics and Acoustics

Roots

Basic sciences

DOI

10.1016/j.jcp.2013.12.027

More information

Created

10/7/2017