High-order partitioned spectral deferred correction solvers for multiphysics problems

by   Daniel Z. Huang, et al.

We present an arbitrarily high-order, conditionally stable, partitioned spectral deferred correction (SDC) method for solving multiphysics problems using a sequence of pre-existing single-physics solvers. This method extends the work in [1, 2], which used implicit-explicit Runge-Kutta methods (IMEX) to build high-order, partitioned multiphysics solvers. We consider a generic multiphysics problem modeled as a system of coupled ordinary differential equations (ODEs), coupled through coupling terms that can depend on the state of each subsystem; therefore the method applies to both a semi-discretized system of partial differential equations (PDEs) or problems naturally modeled as coupled systems of ODEs. The sufficient conditions to build arbitrarily high-order partitioned SDC schemes are derived. Based on these conditions, various of partitioned SDC schemes are designed. The stability of the first-order partitioned SDC scheme is analyzed in detail on a coupled, linear model problem. We show that the scheme is conditionally stable, and under conditions on the coupling strength, the scheme can be unconditionally stable. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems including a simple linear system of ODEs, advection-diffusion-reaction systems, and fluid-structure interaction problems with both incompressible and compressible flows, where we verify the design order of the SDC schemes and study various stability properties. We also directly compare the accuracy, stability, and cost of the proposed partitioned SDC solver with the partitioned IMEX method in [1, 2] on this suite of test problems. The results suggest that the high-order partitioned SDC solvers are more robust than the partitioned IMEX solvers for the numerical examples considered in this work, while the IMEX methods require fewer implicit solves.


page 11

page 14


High order semi-implicit multistep methods for time dependent partial differential equations

We consider the construction of semi-implicit linear multistep methods w...

Very High-Order A-stable Stiffly Accurate Diagonally Implicit Runge-Kutta Methods

A numerical search approach is used to design high-order diagonally impl...

Efficient implementation of partitioned stiff exponential Runge-Kutta methods

Multiphysics systems are driven by multiple processes acting simultaneou...

Recycling Krylov Subspaces for Efficient Partitioned Solution of Aerostructural Adjoint Systems

Robust and efficient solvers for coupled-adjoint linear systems are cruc...

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

Coupled multi-physics problems are encountered in countless applications...

Time-Accurate and highly-Stable Explicit operators for stiff differential equations

Unconditionally stable implicit time-marching methods are powerful in so...

Multirate Partitioned Runge-Kutta Methods for Coupled Navier-Stokes Equations

Earth system models are complex integrated models of atmosphere, ocean, ...

Please sign up or login with your details

Forgot password? Click here to reset