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.
References
- Konovalov AN. Fictitious domain method in filtration problems of a two-phase incompressible fluid, taking into account capillary forces. Chislennye metody mekhaniki sploshnoi sredy. 1972;3(5):52–68. Russian.
- Konovalov AN, Konuh GV, Tsurikov NV. About principles of building iterative processes in the fictitious domain method. In: Variatsionnye metody v zadachakh chislennogo analiza: sbornik nauchnykh trudov. Novosibirsk: Siberian Branch of the Russian Academy of Sciences; 1986. p. 37–52. Russian.
- Konovalov АN. Zadachi fil’tratsii mnogofaznoi neszhimaemoi zhidkosti [Filtration problems of multi-phase incompressible fluid]. Novosibirsk: Nauka; 1988. Russian.
- Vabishevich PN, Gassiev RV, Pulatov PA. Computational realization of fictitious domain method for elliptic equations on the basis of changing-triangular method. Zhurnal vychislitel’noi matematiki i matematicheskoi fiziki. 1987;27(9):1381–1387. Russian.
- Vabishevich PN. Metod fiktivnykh oblastei v zadachakh matematicheskoi fiziki [Fictitious domain method in the problems of mathematical physics]. Moscow: URSS; 2016. Russian.
- Samarsky AA. Teoriya raznostnykh skhem [The theory of finite-difference schemes]. Moscow: Nauka; 1989. Russian.
- Turovets S, Volkov V, Zherdetsky A, Prakonina A, Malony AD. А 3D finite-difference BiCG iterative solver with the Fourier – Jacobi preconditioner for the anisotropic EIT/EEG forward problem. Computational and Mathematical Methods in Medicine. 2014;2014:12. DOI: 10.1155/2014/426902.
- Samarskii AA, Mazhukin VI, Matus PP, Shishkin GI. Monotone difference schemes for equations with mixed derivatives. Matematicheskoe modelirovanie. 2001;13(2):17–26. Russian.
- Rybak IV. Monotone and conservative difference schemes for elliptic equations with mixed derivatives. Mathematical Modelling and Analysis. 2004;9(2):169 –178.
- Volkov VM, Prakonina AU. Finite-difference schemes and iterative methods for multidimensional elliptic equations with mixed derivatives. Vesci Nacyjanal’naj akadjemii navuk Belarusi. Seryja fizikamatjematychnyh navuk. 2018;54(4):454 – 459. Russian. DOI: 10.29235/1561-2430-2018-54-4-454-459.
- Barrett R, Berry M, Chan TF, Demmel J, Donato J, Dongarra J, et al. Templates for the solution of linear systems: building blocks for iterative methods. Philadelphia: SIAM; 1994. 143 p.
- Martynenko SI. Universal multigrid technology for the numerical solution of partial differential equations on structured grids. Vychislitel’nye metody i programmirovanie. 2000;1(1):83–102. Russian.
Copyright (c) 2019 Journal of the Belarusian State University. Mathematics and Informatics
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
The authors who are published in this journal agree to the following:
- The authors retain copyright on the work and provide the journal with the right of first publication of the work on condition of license Creative Commons Attribution-NonCommercial. 4.0 International (CC BY-NC 4.0).
- The authors retain the right to enter into certain contractual agreements relating to the non-exclusive distribution of the published version of the work (e.g. post it on the institutional repository, publication in the book), with the reference to its original publication in this journal.
- The authors have the right to post their work on the Internet (e.g. on the institutional store or personal website) prior to and during the review process, conducted by the journal, as this may lead to a productive discussion and a large number of references to this work. (See The Effect of Open Access.)