The Iterative Method for the Poisson Equation Based on Domain Decomposition

Jie Sun,Rui Zhang

DOI: https://doi.org/10.54097/vd9zpw84

2023-01-01

Abstract:This paper is concerned with 2D Poisson equation. A region-segmentation-based iteration algorithm for the numerical solution of the Poisson equation is presented. The problem can be solved in parallel, where numerical results test our method is effective.

What problem does this paper attempt to address?

An Approach to Numerically Solving the Poisson Equation

Zhichen Feng,Zheng-Mao Sheng

DOI: https://doi.org/10.1088/0031-8949/90/6/065603

2015-01-01

Physica Scripta

Abstract:We introduce an approach for numerically solving the Poisson equation by using a physical model, which is a way to solve a partial differential equation without the finite difference method. This method is especially useful for obtaining the solutions in very many free-charge neutral systems with open boundary conditions. It can be used for arbitrary geometry and mesh style and is more efficient comparing with the widely-used iterative algorithm with multigrid methods. It is especially suitable for parallel computing. This method can also be applied to numerically solving other partial differential equations whose Green functions exist in analytic expression.
Iterative Methods for the Poisson Equation in L-shape Domain Based on Domain Decomposition Method

Xiao-liang Cheng,Abdul Wasim Shaikh

DOI: https://doi.org/10.1016/j.amc.2005.12.023

IF: 4.397

2006-01-01

Applied Mathematics and Computation

Abstract:In this paper we propose the finite difference method for the second-order elliptic equation in the L-shape domain. The scheme is based on the non-overlap domain decomposition method. We use the coarse mesh difference scheme at the interface grids and fine mesh difference scheme at the interior mesh points. We derive the error estimate O(h2+H3) with mesh size h and coarse mesh size H. Then we discuss some iterative methods for our scheme and derive the upper bound of the convergent rate 1−H. Finally we present some numerical experiments.
On an interior-exterior nonoverlapping domain decomposition method for the Poisson--Boltzmann equation

Xuanyu Liu,Yvon Maday,Chaoyu Quan,Hui Zhang

DOI: https://doi.org/10.48550/arxiv.2204.07958

2022-01-01

Abstract: A nonoverlapping domain decomposition method is studied for the linearized Poisson--Boltzmann equation, which is essentially an interior-exterior transmission problem with bounded interior and unbounded exterior. This problem is different from the classical Schwarz alternating method for bounded nonoverlapping subdomains well studied by Lions in 1990, and is challenging due to the existence of unbounded subdomain. To obtain the convergence, a new concept of interior-exterior Sobolev constant is introduced and a spectral equivalence of related Dirichlet-to-Neumann operators is established afterwards. We prove rigorously that the spectral equivalence results in the convergence of interior-exterior iteration. Some numerical simulations are provided to investigate the optimal stepping parameter of iteration and to verify our convergence analysis.
Fast Algorithm for the Three-Dimensional Poisson Equation in Infinite Domains

Chunxiong Zheng,Xiang Ma

DOI: https://doi.org/10.1093/imanum/draa051

2020-01-01

IMA Journal of Numerical Analysis

Abstract:This paper is concerned with a fast finite element method for the three-dimensional Poisson equation in infinite domains. Both the exterior problem and the strip-tail problem are considered. Exact Dirichlet-to-Neumann (DtN)-type artificial boundary conditions (ABCs) are derived to reduce the original infinite-domain problems to suitable truncated-domain problems. Based on the best relative Chebyshev approximation for the square-root function, a fast algorithm is developed to approximate exact ABCs. One remarkable advantage is that one need not compute the full eigensystem associated with the surface Laplacian operator on artificial boundaries. In addition, compared with the modal expansion method and the method based on Pad$\acute{\textrm{e}}$ approximation for the square-root function, the computational cost of the DtN mapping is further reduced. An error analysis is performed and numerical examples are presented to demonstrate the efficiency of the proposed method.
An Additive Overlapping Domain Decomposition Method for the Helmholtz Equation

Wei Leng,Lili Ju

DOI: https://doi.org/10.1137/18M1196170

2018-07-11

Abstract:In this paper, we propose and analyze an additive domain decomposition method (DDM) for solving the high-frequency Helmholtz equation with the Sommerfeld radiation condition. In the proposed method, the computational domain is partitioned into structured subdomains along all spatial directions, and each subdomain contains an overlapping region for source transferring. At each iteration all subdomain PML problems are solved completely in parallel, then all horizontal, vertical and corner directional residuals on each subdomain are passed to its corresponding neighbor subdomains as the source for the next iteration. This DDM method is highly scalable in nature and theoretically shown to produce the exact solution for the PML problem defined in ${\mathbb{R}}^2$ in the constant medium case. A slightly modified version of the method for bounded truncated domains is also developed for its use in practice and an error estimate is rigorously proved. Various numerical experiments in two and three dimensions are conducted on the supercomputer "Tianhe-2 Cluster" to verify the theoretical results and demonstrate excellent performance of the proposed method as an iterative solver or a preconditioner.

Numerical Analysis
Domain Decomposition Method for Poisson--Boltzmann Equations based on Solvent Excluded Surface

Abhinav Jha,Benjamin Stamm

2024-04-05

Abstract:In this paper, we develop a domain decomposition method for the nonlinear Poisson-Boltzmann equation based on a solvent-excluded surface widely used in computational chemistry. The model relies on a nonlinear equation defined in $\mathbb{R}^3$ with a space-dependent dielectric permittivity and an ion-exclusion function that accounts for steric effects. Potential theory arguments transform the nonlinear equation into two coupled equations defined in a bounded domain. Then, the Schwarz decomposition method is used to formulate local problems by decomposing the cavity into overlapping balls and only solving a set of coupled sub-equations in each ball. The main novelty of the proposed method is the introduction of a hybrid linear-nonlinear solver used to solve the equation. A series of numerical experiments are presented to test the method and show the importance of the nonlinear model.

Numerical Analysis,Number Theory
A Multiscale Domain Decomposition Algorithm For Boundary Value Problems For Eikonal Equations

Lindsay Martin,Richard Tsai

DOI: https://doi.org/10.48550/arXiv.1805.04198

2018-06-12

Abstract:In this paper, we present a new multiscale domain decomposition algorithm for computing solutions of static Eikonal equations. The new method is an iterative two-scale method that uses a parareal-like update scheme in combination with standard Eikonal solvers. The purpose of the two scales is to accelerate convergence and maintain accuracy. We adapt a weighted version of the parareal method for stability, and the optimal weights are studied via a model problem. Numerical examples are given to demonstrate the method.

Numerical Analysis
Block iterative methods for linear algebraic equation and domain decomposition method for incompressible viscous flow

Anlu Ren,Jianfeng Zou,Yongxia Zhou,Guang Wang Li

2002-01-01

Journal of Hydrodynamics

Abstract:It was proved numerically that the Domain Decomposition Method CDDM) with one layer overlapping grids is identical to the block iterative method of linear algebra equations. The results obtained by using DDM could be in reasonable agreement with the results of full-domain simulation. With the three dimensional solver developed by the authors, the flow field in a pipe was simulated by using the full-domain DDM with one layer overlapping grids and with patched grids respectively. Both of the two cases led to the convergent solution. Further research shows the superiority of the DDM with one layer overlapping grids to the DDM with patched grids. A comparison between the numerical results obtained by the authors and the experimental results given by Enayet[3] shows that the numerical results are reasonable.
A Numerical Study of the Poisson Equation with Arbitrary Domains Based on the Extrapolation Technique

Hui Zhang,Pan Yu,Xueying Shi,Lin Li

DOI: https://doi.org/10.1088/1742-6596/2455/1/012017

2023-01-01

Journal of Physics Conference Series

Abstract:Abstract The Poisson equation with arbitrary domains is considered in this paper using an extrapolation technique. The two-dimensional problem is primarily examined, and the extrapolation formulas O2 and O4 are thoroughly deduced. Not only the effect of different extrapolations on boundary point errors is analyzed, but a strategy for selecting extrapolation points is also provided. Finally, numerical accuracy results are presented.
Domain Decomposition Methods for the Monge-Ampère equation

Yassine Boubendir,Jake Brusca,Brittany Froese Hamfeldt,Tadanaga Takahashi

2023-06-03

Abstract:We introduce a new overlapping Domain Decomposition Method (DDM) to solve the fully nonlinear Monge-Ampère equation. While DDMs have been extensively studied for linear problems, their application to fully nonlinear partial differential equations (PDE) remains limited in the literature. To address this gap, we establish a proof of global convergence of these new iterative algorithms using a discrete comparison principle argument. Several numerical tests are performed to validate the convergence theorem. These numerical experiments involve examples of varying regularity. Computational experiments show that method is efficient, robust, and requires relatively few iterations to converge. The results reveal great potential for DDM methods to lead to highly efficient and parallelizable solvers for large-scale problems that are computationally intractable using existing solution methods.

Numerical Analysis
A Novel Domain Decomposition Method for Highly Oscillating Partial Differential Equations

Y. Duan,P. F. Tang,T. Z. Huang

DOI: https://doi.org/10.1016/j.enganabound.2009.05.002

IF: 3.25

2009-01-01

Engineering Analysis with Boundary Elements

Abstract:This paper is devoted to designing a novel domain decomposition method (DDM) for highly oscillating partial differential equations (PDE), especially, where the asymmetric meshless collocation method using radial basis functions (RBF), also Kansa's method is applied for a numerical solutions. It is found that the numerical error become worse when the original solution become more oscillating. To conquer this defect, we use a novel domain decomposition method which is motivated by time parallel algorithm. This DDM is based on a decomposition of computational domain by a coarse centers and a finer distribution of distinct centers. A corrector is designed to obtain better numerical solution after several iteration. Theoretical analysis and numerical examples are given to demonstrate the accuracy and efficiency of the proposed algorithm.
Convergence of the Point Integral Method for the Poisson Equation with Dirichlet Boundary on Point Cloud

Zuoqiang Shi,Jian Sun

DOI: https://doi.org/10.48550/arxiv.1312.4424

2013-01-01

Abstract:The Poisson equation on manifolds plays an fundamental role in many applications. Recently, we proposed a novel numerical method called the Point Integral method (PIM) to solve the Poisson equations on manifolds from point clouds. In this paper, we prove the convergence of the point integral method for solving the Poisson equation with the Dirichlet boundary condition.
Numerical study of a domain decomposition method for a two-scale linear transport equation

Xu Yang,François Golse,Zhongyi Huang,Shi Jin

DOI: https://doi.org/10.3934/nhm.2006.1.143

2006-01-01

Abstract:We perform a numerical study on a domain decomposition method proposed in [13] for the linear transport equation between a diffusive and a non-diffusive region. This method avoids iterating the diffusion and transport solutions as in a typical domain decomposition method. Our numerical results, in both one and two space dimensions, confirm the theoretical analysis of [13]. We also provide an improved second order method that provides a more accurate numerical solution than that proposed in [13].
Convergence of Finite Integral Method for Poisson Equation and Corresponding Eigenproblems with Dirichlet Boundary

Zuoqiang Shi,Jian Sun

2014-01-01

Abstract:The Laplace-Beltrami operator (LBO) is a fundamental object associated to Riemannian manifolds, which encodes all intrinsic geometry of the manifolds and has many desirable properties. Recently, we proposed a novel numerical method, Point Integral method (PIM), to discretize the Laplace-Beltrami operator on point clouds \cite{LSS}. In this paper, we analyze the convergence of Point Integral method (PIM) for Poisson equation with Neumann boundary condition on submanifolds isometrically embedded in Euclidean spaces.
An overlapping domain decomposition splitting algorithm for stochastic nonlinear Schroedinger equation

Lihai Ji

2023-09-07

Abstract:A novel overlapping domain decomposition splitting algorithm based on a Crank-Nisolson method is developed for the stochastic nonlinear Schroedinger equation driven by a multiplicative noise with non-periodic boundary conditions. The proposed algorithm can significantly reduce the computational cost while maintaining the similar conservation laws. Numerical experiments are dedicated to illustrating the capability of the algorithm for different spatial dimensions, as well as the various initial conditions. In particular, we compare the performance of the overlapping domain decomposition splitting algorithm with the stochastic multi-symplectic method in [S. Jiang, L. Wang and J. Hong, Commun. Comput. Phys., 2013] and the finite difference splitting scheme in [J. Cui, J. Hong, Z. Liu and W. Zhou, J. Differ. Equ., 2019]. We observe that our proposed algorithm has excellent computational efficiency and is highly competitive. It provides a useful tool for solving stochastic partial differential equations.

Numerical Analysis
Solution of 1D Poisson Equation with Neumann-Dirichlet and Dirichlet-Neumann Boundary Conditions, Using the Finite Difference Method

K. Talla,Serigne Gueye,C. Mbow

DOI: https://doi.org/10.4236/JEMAA.2014.610031

2014-09-10

Abstract:An innovative, extremely fast and accurate method is presented for Neumann-Dirichlet and Dirichlet-Neumann boundary problems for the Poisson equation, and the diffusion and wave equation in quasi-stationary regime; using the finite difference method, in one dimensional case. Two novels matrices are determined allowing a direct and exact formulation of the solution of the Poisson equation. Verification is also done considering an interesting potential problem and the sensibility is determined. This new method has an algorithm complexity of O(N), its truncation error goes like O(h2), and it is more precise and faster than the Thomas algorithm.

Mathematics,Physics
A parallel domain decomposition method for solving elliptic equations on manifolds

Lizhen Qin,Feng Wang,Yun Wang

2024-02-01

Abstract:We propose a new numerical domain decomposition method for solving elliptic equations on compact Riemannian manifolds. One advantage of this method is its ability to bypass the need for global triangulations or grids on the manifolds. Additionally, it features a highly parallel iterative scheme. To verify its efficacy, we conduct numerical experiments on some $4$-dimensional manifolds without and with boundary.

Numerical Analysis
A modified domain decomposition spectral collocation method for parabolic partial differential equations

Wei-Hua Luo,Liang Yin,Jun Guo

DOI: https://doi.org/10.3934/nhm.2024041

2024-09-15

Networks & Heterogeneous Media

Abstract:In this paper, utilizing Legendre polynomials as the basis functions in both space and time, we present a modified domain decomposition spectral method for 2-dimensional parabolic partial differential equations. For solving the obtained linear/nonlinear algebraic equations, a dimension expanding preconditioner is applied employing the obtained saddle construction of the coefficient matrix. Numerical examples are given to show the performance of the presented method and the efficiency of the preconditioner.

mathematics, interdisciplinary applications
A rapidly converging domain decomposition algorithm for a time delayed parabolic problem with mixed type boundary conditions exhibiting boundary layers

Aakansha,Sunil Kumar,Higinio Ramos

DOI: https://doi.org/10.1007/s12190-024-01987-8

2024-02-04

Journal of Applied Mathematics and Computing

Abstract:A rapidly converging domain decomposition algorithm is introduced for a time delayed parabolic problem with mixed type boundary conditions exhibiting boundary layers. Firstly, a space-time decomposition of the original problem is considered. Subsequently, an iterative process is proposed, wherein the exchange of information to neighboring subdomains is accomplished through the utilization of piecewise-linear interpolation. It is shown that the algorithm provides uniformly convergent numerical approximations to the solution. Our analysis utilizes some novel auxiliary problems, barrier functions, and a new maximum principle result. More importantly, we showed that only one iteration is needed for small values of the perturbation parameter. Some numerical results supporting the theory and demonstrating the effectiveness of the algorithm are presented.

mathematics, applied
Solution of the mixed boundary problem for the Poisson equation on two-dimensional irregular domains

M. M. Chuiko,O. M. Korolyova

DOI: https://doi.org/10.37661/1816-0301-2023-20-2-111-120

2023-06-29

Informatics

Abstract:Objectives . A finite-difference computational algorithm is proposed for solving a mixed boundary-value problem for the Poisson equation given in two-dimensional irregular domains. Methods . To solve the problem, generalized curvilinear coordinates are used. The physical domain is mapped to the computational domain (unit square) in the space of generalized coordinates. The original problem is written in curvilinear coordinates and approximated on a uniform grid in the computational domain.The obtained results are mapped on non-uniform boundary-fitted difference grid in the physical domain. Results . The second order approximations of mixed Neumann-Dirichlet boundary conditions for the Poisson equation in the space of generalized curvilinear coordinate are constructed. To increase the order of Neumann condition approximations, an approximation of the Poisson equation on the boundary of the domain is used. Conclusions . To solve a mixed boundary value problem for the Poisson equation in two-dimensional irregular domains, the computational algorithm of second-order accuracy is constructed. The generalized curvilinear coordinates are used. The results of numerical experiments, which confirm the second order accuracy of the computational algorithm, are presented.

The Iterative Method for the Poisson Equation Based on Domain Decomposition

An Approach to Numerically Solving the Poisson Equation

Iterative Methods for the Poisson Equation in L-shape Domain Based on Domain Decomposition Method

On an interior-exterior nonoverlapping domain decomposition method for the Poisson--Boltzmann equation

Fast Algorithm for the Three-Dimensional Poisson Equation in Infinite Domains

An Additive Overlapping Domain Decomposition Method for the Helmholtz Equation

Domain Decomposition Method for Poisson--Boltzmann Equations based on Solvent Excluded Surface

A Multiscale Domain Decomposition Algorithm For Boundary Value Problems For Eikonal Equations

Block iterative methods for linear algebraic equation and domain decomposition method for incompressible viscous flow

A Numerical Study of the Poisson Equation with Arbitrary Domains Based on the Extrapolation Technique

Domain Decomposition Methods for the Monge-Ampère equation

A Novel Domain Decomposition Method for Highly Oscillating Partial Differential Equations

Convergence of the Point Integral Method for the Poisson Equation with Dirichlet Boundary on Point Cloud

Numerical study of a domain decomposition method for a two-scale linear transport equation

Convergence of Finite Integral Method for Poisson Equation and Corresponding Eigenproblems with Dirichlet Boundary

An overlapping domain decomposition splitting algorithm for stochastic nonlinear Schroedinger equation

Solution of 1D Poisson Equation with Neumann-Dirichlet and Dirichlet-Neumann Boundary Conditions, Using the Finite Difference Method

A parallel domain decomposition method for solving elliptic equations on manifolds

A modified domain decomposition spectral collocation method for parabolic partial differential equations

A rapidly converging domain decomposition algorithm for a time delayed parabolic problem with mixed type boundary conditions exhibiting boundary layers

Solution of the mixed boundary problem for the Poisson equation on two-dimensional irregular domains