Splitting Schemes for Coupled Differential Equations: Block Schur-Based Approaches and Partial Jacobi Approximation

12/21/2022
by   Roberto Nuca, et al.
0

Coupled multi-physics problems are encountered in countless applications and pose significant numerical challenges. Although monolithic approaches offer possibly the best solution strategy, they often require ad-hoc preconditioners and numerical implementations. Sequential (also known as splitted, partitioned or segregated) approaches are iterative methods for solving coupled problems where each equation is solved independently and the coupling is achieved through iterations. These methods offer the possibility to flexibly add or remove equations from a model and to rely on existing black-box solvers for every specific equation. Furthermore, when problems are non-linear, inner iterations need to be performed even in monolithic solvers, therefore making a sequential iterative approach a viable alternative. The cost of running inner iterations to achieve the coupling, however, could easily becomes prohibitive, or, in some cases the iterations might not converge. In this work we present a general formulation of splitting schemes for continuous operators, with arbitrary implicit/explicit splitting, like in standard iterative methods for linear systems. By introducing a generic relaxation operator we find the conditions for the convergence of the iterative schemes. We show how the relaxation operator can be thought as a preconditioner and constructed based on an approximate Schur-complement. We propose a Schur-based Partial Jacobi relaxation operator to stabilise the coupling and show its effectiveness. Although we mainly focus on scalar-scalar linear problems, most results are easily extended to non-linear and higher-dimensional problems. Numerical tests (1D and 2D) for two PDE systems, namely the Dual-Porosity model and a Quad-Laplacian operator, are carried out to confirm the theoretical results.

READ FULL TEXT

page 18

page 20

research
07/30/2019

Iterative and Non-iterative Splitting approach of a stochastic Burgers' equation

In this paper we present iterative and noniterative splitting methods, w...
research
09/03/2019

High-order partitioned spectral deferred correction solvers for multiphysics problems

We present an arbitrarily high-order, conditionally stable, partitioned ...
research
12/11/2022

DOSnet as a Non-Black-Box PDE Solver: When Deep Learning Meets Operator Splitting

Deep neural networks (DNNs) recently emerged as a promising tool for ana...
research
12/14/2019

Operator splitting schemes for American options under the two-asset Merton jump-diffusion model

This paper deals with the efficient numerical solution of the two-dimens...
research
08/30/2021

Partitioned Coupling vs. Monolithic Block-Preconditioning Approaches for Solving Stokes-Darcy Systems

We consider the time-dependent Stokes-Darcy problem as a model case for ...
research
06/29/2021

Numerical approximation of singular Forward-Backward SDEs

In this work, we study the numerical approximation of a class of singula...
research
03/14/2023

Symmetric integration of the 1+1 Teukolsky equation on hyperboloidal foliations of Kerr spacetimes

This work outlines a fast, high-precision time-domain solver for scalar,...

Please sign up or login with your details

Forgot password? Click here to reset