Discontinuous bubble scheme for elliptic problems with jumps in the solution

We propose a new numerical method to solve an elliptic problem with jumps both in the solution and derivative along an interface. By considering a suitable function which has the same jumps as the solution, we transform the problem into one without jumps. Then we apply the immersed finite element method in which we allow uniform meshes so that the interface may cut through elements to discretize the problem as introduced in [1–3]. Some convenient way of approximating the jumps of the solution by piecewise linear functions is suggested. Our method can also handle the case when the interface passes through grid points. We believe this paper presents the first resolution of such cases. Numerical experiments for various problems show second-order convergence in L2 and first order in H1-norms. Moreover, the convergence order is very robust for all problems tested.