Iterative realization of finite difference schemes in the fictitious domain method for elliptic problems with mixed derivatives

Abstract

Development of efficient finite difference schemes and iterative methods for solving anisotropic diffusion problems in an arbitrary geometry domain is considered. To simplify the formulation of the Neumann boundary conditions, the method of fictitious domains is used. On the example of a two-dimensional model problem of potential distribution in an isolated anisotropic ring conductor a comparative efficiency analysis of some promising finite-difference schemes and iterative methods in terms of their compatibility with the fictitious domain method is carried out. On the basis of numerical experiments empirical estimates of the asymptotic dependence of the convergence rate of the biconjugate gradient method with Fourier – Jacobi and incomplete LU factorization preconditioners on the step size and the value of the small parameter determining the continuation of the conductivity coefficient in the fictitious domain method are obtained. It is shown, that for one of the considered schemes the Fourier – Jacobi preconditioner is spectrally optimal and allows to eliminate the asymptotical dependence of the iterations number to achieve a given accuracy both on the value of the step size and the value of the small parameter in the fictitious domain method.