Parallel version of the spectral preconditioner for solving the Poisson equation
Authors
-
Aleksei A. Manaev
-
Vadim V. Lisitsa
Keywords:
Poisson equation
preconditioner
spectral decomposition
domain decomposition
Abstract
This paper presents a parallel version of the spectral preconditioner for the numerical solution of the Poisson equation in heterogeneous media. The action of the spectral preconditioner is based on the eigendecomposition of operators approximating derivatives and the corresponding boundary conditions along one spatial direction, followed by the solution of a series of independent one-dimensional problems along another spatial direction. The parallel version is based on the domain decomposition method using a polynomial approximation of the inverse matrix and a block-diagonal matrix, where each block is analogous in its action to the spectral preconditioner. It is shown that the proposed approach for constructing the preconditioner can significantly accelerate computations compared to standard preconditioning methods.
Section
Methods and algorithms of computational mathematics and their applications
References
- H. Andrä, N. Combaret, J. Dvorkin, et al., “Digital rock physics benchmarks — part II: Computing effective properties,” Computers and Geosciences. 50, 33–43 (2013).
doi 10.1016/j.cageo.2012.09.008
- Y. Bazaikin, B. Gurevich, S. Iglauer, et al., “Effect of CT image size and resolution on the accuracy of rock property estimates,” Journal of Geophysical Research: Solid Earth 122 (5), 3635–3647 (2017).
doi 10.1002/2016JB013575
- N. M. Evstigneev, O. I. Ryabkov, and K. M. Gerke, “Stationary Stokes solver for single-phase flow in porous media: A blastingly fast solution based on Algebraic Multigrid Method using GPU,” Advances in Water Resources 171, Article Number 104340 (2023).
doi 10.1016/j.advwatres.2022.104340
- T. S. Khachkova, V. V. Lisitsa, E. A. Gondul, et al., “Two-phase flow simulation algorithm for numerical estimation of relative phase permeability curves of porous materials,” Russian Journal of Numerical Analysis and Mathematical Modelling 39 (4), 209–221 (2024).
doi 10.1515/rnam-2024-0020
- M. Li, S. Foroughi, J. Zhao, et al., “Image-based pore-scale modelling of the effect of wettability on breakthrough capillary pressure in gas diffusion layers,” Journal of Power Sources 584, Article Number 233539 (2023).
doi 10.1016/j.jpowsour.2023.233539
- V. Shulakova, M. Pervukhina, T. M. Müller, et al., “Computational elastic up-scaling of sandstone on the basis of X-ray micro-tomographic images,” Geophysical Prospecting 61 (2), 287–301 (2012).
doi 10.1111/j.1365-2478.2012.01082.x
- G. V. Reshetova and T. S. Khachkova, “A numerical method to estimate the effective elastic moduli of rocks from two- and three-dimensional digital images of rock core samples,” Numerical Methods and Programming 18 (4), 416–433 (2017).
doi 10.26089/NumMet.v18r435
- T. S. Khachkova, V. V. Lisitsa, G. V. Reshetova, and V. A. Tcheverda, “Numerical estimation of electrical resistivity in digital rocks using GPUs,” Numerical Methods and Programming 21 (3), 306–318 (2020).
doi 10.26089/NumMet.v21r326
- C. Dorn and M. Schneider, “Lippmann-Schwinger solvers for the explicit jump discretization for thermal computational homogenization problems,” International Journal for Numerical Methods in Engineering 118 (11), 631–653 (2019).
doi 10.1002/nme.6030
- Y. Saad, Iterative Methods for Sparse Linear Systems(Society for Industrial and Applied Mathematics, Philadelphia, 2003).
- E. Haber, U. M. Ascher, D. A. Aruliah, and D. W. Oldenburg, “Fast Simulation of 3D Electromagnetic Problems Using Potentials,” Journal of Computational Physics 163 (1), 150–171 (2000).
doi 10.1006/jcph.2000.6545
- V. Kostin, S. Solovyev, A. Bakulin, and M. Dmitriev, “Direct frequency-domain 3D acoustic solver with intermediate data compression benchmarked against time-domain modeling for full-waveform inversion applications,” Geophysics 84 (4), 207–219 (2019).
doi 10.1190/geo2018-0465.1
- K. V. Voronin and S. A. Solovyev, “Solution of the Helmholtz problem using the preconditioned low-rank approximation technique,” Numerical Methods and Programming 16 (2), 268–280 (2015).
doi 10.26089/NumMet.v16r226
- H. Johansen and P. Colella, “A Cartesian Grid Embedded Boundary Method for Poisson’s Equation on Irregular Domains,” Journal of Computational Physics 147 (1), 60–85 (1998).
doi 10.1006/jcph.1998.5965
- K. Stüben, “A review of algebraic multigrid,” Journal of Computational and Applied Mathematics 128 (1–2), 281–309 (2001).
doi 10.1016/S0377-0427(00)00516-1
- C. Dick, M. Rogowsky, and R. Westermann, “Solving the Fluid Pressure Poisson Equation Using Multigrid—Evaluation and Improvements,” IEEE Transactions on Visualization and Computer Graphics 22 (11), 2480–2492 (2015).
doi 10.1109/TVCG.2015.2511734
- T. Khachkova, V. Lisitsa, G. Reshetova, and V. Tcheverda, “GPU-based algorithm for evaluating the electrical resistivity of digital rocks,” Computers and Mathematics with Applications 82 (9), 200–211 (2021).
doi 10.1016/j.camwa.2020.11.005
- A. A. Manaev and V. V. Lisitsa, “Spectral preconditioner for solving the Poisson equation,” Numerical Methods and Programming 26 (2), 111–128 (2025).
doi 10.26089/NumMet.v26r208
- A. A. Manaev, T. S. Khachkova, and V. V. Lisitsa, “Numerical algorithm to estimate formation factor of digital rocks,” Numerical Methods and Programming 26 (4), 479–501 (2025).
doi 10.26089/NumMet.v26r432
- J. J. Hasbestan, C. N. Xiao, and I. Senocak, “PittPack: An open-source Poisson’s equation solver for extreme-scale computing with accelerators,” Computer Physics Communications 254 (3), Article Number 107272 (2020).
doi 10.1016/j.cpc.2020.107272
- V. Lisitsa, A. Manaev, T. Khachkova, et al., “GPU-oriented numerical algorithm to estimate formation factor of porous materials,” Journal of Computational Science 96, Article Number 102829 (2026).
doi 10.1016/j.jocs.2026.102829
- A. A. Samarskii, The theory of difference schemes(CRC Press, 2001; Nauka, Moscow, 1983).
- E. Haber and U. M. Ascher, “Fast finite volume simulation of 3D electromagnetic problems with highly discontinuous coefficients,” SIAM Journal on Scientific Computing 22 (6), 1943–1961 (2001).
doi 10.1137/S1064827599360741
- S. Chandrasekaran, P. Dewilde, M. Gu, and N. Somasunderam, “On the Numeral Rank of the Off-Diagonal Blocks of Schur Complements of Discretized Elliptic PDEs,” SIAM Journal on Matrix Analysis and Applications 31 (5), 2261–2290 (2010).
doi 10.1137/090775932
- S. A. Solovyev, “Application of the low-rank approximation technique in the Gauss elimination method for sparse linear systems,” Numerical Methods and Programming 15 (3), 441–460 (2014).
https://www.mathnet.ru/eng/vmp263 Cited July 09, 2026.
- A. Pleshkevich, D. Vishnevskiy, and V. Lisitsa, “Sixth-order accurate pseudo-spectral method for solving one-way wave equation,” Applied Mathematics and Computation 359, 34–51 (2019).
doi 10.1016/j.amc.2019.04.029
- J. D. Hyman and C. L. Winter, “Stochastic generation of explicit pore structures by thresholding Gaussian random fields,” Journal of Computational Physics 277, 16–31 (2014).
doi 10.1016/j.jcp.2014.07.046