The Schur complement method as a fast parallel solver for elliptic partial differential equations in oceanography

DOI10.1002/(SICI)1099-1506(199909)6:6%3C497::AID-NLA171%3E3.0.CO;2-RMaRDI QIDQ2760367zbMATH OpenOpenAlexFDO

Authors Natalja Rakowsky

Publication date 19 December 2001

Published in Numerical Linear Algebra with Applications (Search for Journal in Brave)

Full work available at URL https://doi.org/10.1002/(sici)1099-1506(199909)6:6%3C497::aid-nla171%3E3.0.co;2-r

zbMATH Keywords

direct methoddomain decompositionexplicit time-stepping schemehydrostatic primitive equationsice-ocean model BRIOSoceanographyparallel elliptic solverregular rectangular gridSchur complement methodtwo-dimensional elliptic PDEvariable anisotropic coefficientsvertically integrated stream function

Mathematics Subject Classification ID

Direct numerical methods for linear systems and matrix inversion (65F05) Multigrid methods; domain decomposition for boundary value problems involving PDEs (65N55) Parallel numerical computation (65Y05) Basic methods in fluid mechanics (76M99) Computational methods for problems pertaining to geophysics (86-08) Hydrology, hydrography, oceanography (86A05) Glaciology (86A40)

Review text

The method is described in the context of ice-ocean model BRIOS (Bremerhaven Regional Ice-Ocean Simulation System). After an introduction, hydrostatic primitive equations (a set of three-dimensional time-dependent nonlinear coupled PDEs) are presented and further modified by the introduction of a vertically integrated stream function. The stream function solves the Laplace or, after a transformation of coordinate system, a two-dimensional elliptic PDE with highly variable anisotropic coefficients on a regular rectangular grid. The latter equation is a part of an explicit time-stepping scheme in BRIOS. The equation is solved in parallel by means of domain decomposition approach which gives rise to the application of Schur complement method. It uses a direct solver and an explicitly inverted Schur complement matrix, as the direct approach turns out to be advantageous in problems BRIOS solves. Two solvers are presented and, together with other BRIOS features, are tested and discussed.

Recommendations

• On the Schwarz alternating method for oceanic models on parallel computers
• Algebraic two-level preconditioners for the Schur complement method
• scientific article; zbMATH DE number 2247848
• A fast semi-direct method for the numerical solution of non-separable elliptic equations in irregular domains
• Parallel ocean general circulation modeling

Cites work

• A capacitance matrix method for Dirichlet problem on polygon region
• A numerical method for the Study of the circulation of the world ocean
• A semi-spectral primitive equation ocean circulation model using vertical sigma and orthogonal curvilinear horizontal coordinates
• An application of the capacitance matrix method to accommodate masked land areas and island circulations in a primitive equation ocean model
• Iterative Methods for the Solution of Elliptic Problems on Regions Partitioned into Substructures
• Matrix-dependent prolongations and restrictions in a blackbox multigrid solver
• The Interface Probing Technique in Domain Decomposition

Cited in

(3)

• A scalable parallel Poisson solver for three-dimensional problems with one periodic direction
• scientific article; zbMATH DE number 2052752 (Why is no real title available?)
• Parallel direct Poisson solver for discretisations with one Fourier diagonalisable direction

This page was built for publication: The Schur complement method as a fast parallel solver for elliptic partial differential equations in oceanography

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