A Finite Difference Method for Two-Phase Parabolic Obstacle-like Problem

• Main Author: Avetik Arakelyan
• Format: Article
• Language: English
• Published: Republic of Armenia National Academy of Sciences 2015-05-01
• Series: Armenian Journal of Mathematics
• Online Access: http://www.armjmath.sci.am/index.php/ajm/article/view/109

In this paper we treat the numerical approximation of the two-phase parabolic obstacle-like problem: \[\Delta u -u_t=\lambda^+\cdot\chi_{\{u>0\}}-\lambda^-\cdot\chi_{\{u<0\}},\quad (t,x)\in (0,T)\times\Omega,\] where $T < \infty, \lambda^+ ,\lambda^- > 0$ are Lipschitz continuous functions, and $\Omega\subset\mathbb{R}^n$ is a bounded domain. We introduce a certain variation form, which allows us to define a notion of viscosity solution. We use defined viscosity solutions framework to apply Barles-Souganidis theory. The numerical projected Gauss-Seidel method is constructed. Although the paper is devoted to the parabolic version of the two-phase obstacle-like problem, we prove convergence of the discretized scheme to the unique viscosity solution for both two-phase parabolic obstacle-like and standard two-phase membrane problem. Numerical simulations are also presented.