An efficient method for solving highly anisotropic elliptic equations

DOI10.1016/J.JCP.2011.06.022MaRDI QIDQ414000zbMATH OpenOpenAlexFDO

Authors Edward Santilli, Alberto Scotti

Publication date 8 May 2012

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1103.0298

zbMATH Keywords

boundary value problems for second-order elliptic systemscondition numberconvergenceeigenfunction expansionfinite difference methoditerative methods for linear systemsmulti-scale methodsNavier-Stokes equationNeumann boundary conditionsnumerical examplesPoisson equationpreconditioners for iterative methodssolution of discretized equations

Mathematics Subject Classification ID

Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Boundary value problems for second-order elliptic equations (35J25) Navier-Stokes equations (35Q30) Preconditioners for iterative methods (65F08) Iterative numerical methods for linear systems (65F10) Finite difference methods for boundary value problems involving PDEs (65N06)

Abstract: Solving elliptic PDEs in more than one dimension can be a computationally expensive task. For some applications characterised by a high degree of anisotropy in the coefficients of the elliptic operator, such that the term with the highest derivative in one direction is much larger than the terms in the remaining directions, the discretized elliptic operator often has a very large condition number - taking the solution even further out of reach using traditional methods. This paper will demonstrate a solution method for such ill-behaved problems. The high condition number of the D-dimensional discretized elliptic operator will be exploited to split the problem into a series of well-behaved one and (D-1)-dimensional elliptic problems. This solution technique can be used alone on sufficiently coarse grids, or in conjunction with standard iterative methods, such as Conjugate Gradient, to substantially reduce the number of iterations needed to solve the problem to a specified accuracy. The solution is formulated analytically for a generic anisotropic problem using arbitrary coordinates, hopefully bringing this method into the scope of a wide variety of applications.

Recommendations

• Résolution numérique d'un problème elliptique fortement anisotrope en deux dimensions par une méthode de paramétrisation. (Numerical resolution of a highly anisotropic two-dimensional elliptic problem by a parameterization method).
• scientific article; zbMATH DE number 4213374
• A hybrid method for anisotropic elliptic problems based on the coupling of an asymptotic-preserving method with the asymptotic limit model
• Numerical convergence of a parameterisation method for the solution of a highly anisotropic two-dimensional elliptic problem
• Solution of high order compact discretized 3D elliptic partial differential equations by an accelerated multigrid method

Cites work

• A numerical study of optimized sparse preconditioners
• Accurate projection methods for the incompressible Navier-Stokes equations
• An asymptotic preserving scheme for strongly anisotropic elliptic problems
• An asymptotic preserving scheme for the Euler equations in a strong magnetic field
• An asymptotic-preserving method for highly anisotropic elliptic equations based on a micro-macro decomposition
• Fourier Analysis of Iterative Methods for Elliptic pr
• scientific article; zbMATH DE number 991835 (Why is no real title available?)
• scientific article; zbMATH DE number 1049347 (Why is no real title available?)
• scientific article; zbMATH DE number 5245629 (Why is no real title available?)
• Iterative Krylov Methods for Large Linear Systems
• Monotone finite volume schemes for diffusion equations on unstructured triangular and shape-regular polygonal meshes
• Preserving monotonicity in anisotropic diffusion

Cited in

(6)

• High-order mimetic finite differences for anisotropic elliptic equations
• Deflation-accelerated preconditioning of the Poisson-Neumann Schur problem on long domains with a high-order discontinuous element-based collocation method
• An efficient high-order numerical solver for diffusion equations with strong anisotropy
• The stratified ocean model with adaptive refinement (SOMAR)
• Axial Green's function method for multi-dimensional elliptic boundary value problems
• Some Efficient Algorithms for Anisotropic Laplace-like Equation over an Exterior Two-Dimensional Domain

Describes a project that uses

Uses Software

• PETSc

This page was built for publication: An efficient method for solving highly anisotropic elliptic equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q414000)