Multigrid Poisson Solver for Complex Geometries Using Finite Difference Method
Abstract
We present an efficient numerical method, inspired by transformation optics, for solving the Poisson equation in complex and arbitrarily shaped geometries. The approach operates by mapping the physical domain to a uniform computational domain through coordinate transformations, which can be applied either to the entire domain or selectively to specific boundaries inside the domain. This flexibility allows both homogeneous (Laplace equation) and inhomogeneous (Poisson equation) problems to be solved efficiently using iterative or fast direct solvers, with only the material parameters and source terms modified according to the transformation. The method is formulated within a finite difference framework, where the modified material properties are computed from the coordinate transformation equations, either analytically or numerically. This enables accurate treatment of arbitrary geometric shapes while retaining the simplicity of a uniform grid solver. Numerical experiments confirm that the method achieves second-order accuracy , and offers a straightforward pathway to integrate fast solvers such as multigrid methods on the uniform computational grid.
Cite
@article{arxiv.2602.01888,
title = {Multigrid Poisson Solver for Complex Geometries Using Finite Difference Method},
author = {Deepak Gautam and Bhooshan Paradkar},
journal= {arXiv preprint arXiv:2602.01888},
year = {2026}
}