Abstract:This paper introduces a novel approach for the construction of bulk--surface splitting schemes for semi-linear parabolic partial differential equations with dynamic boundary conditions. The proposed construction is based on a reformulation of the system as a partial differential--algebraic equation and the inclusion of certain delay terms for the decoupling. To obtain a fully discrete scheme, the splitting approach is combined with finite elements in space and a BDF discretization in time. Within this paper, we focus on the second-order case, resulting in a $3$-step scheme. We prove second-order convergence under the assumption of a weak CFL-type condition and confirm the theoretical findings by numerical experiments. Moreover, we illustrate the potential for higher-order splitting schemes numerically.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to construct an efficient volume - surface splitting algorithm in semi - linear parabolic partial differential equations with dynamic boundary conditions. Specifically, the author proposes a new method to construct a second - order volume - surface splitting scheme for such problems. This method is based on reformulating the system as partial differential - algebraic equations and decoupling by introducing certain delay terms. To obtain a fully discrete scheme, the splitting method combines the finite element method in space and the BDF discretization in time. The focus of this paper is on the second - order case, resulting in a three - step scheme. The author proves that the scheme has second - order convergence under the assumption of a weak CFL condition and verifies the theoretical results through numerical experiments. In addition, the potential of higher - order splitting schemes is also demonstrated through numerical experiments. ### Key Points Summary: 1. **Problem Background**: - Dynamic boundary conditions have received increasing attention in recent years because they can reflect the effective properties of the domain surface. - In the parabolic setting, dynamic boundary conditions are particularly suitable for modeling heat sources or heat transfer on the boundary. - Mathematically, dynamic boundary conditions define differential equations on the boundary, so that the overall system can be interpreted as a coupled system. 2. **Research Objectives**: - Construct an efficient volume - surface splitting algorithm, especially for semi - linear parabolic partial differential equations with dynamic boundary conditions. - Propose a new method to achieve higher accuracy and efficiency by introducing delay terms for decoupling. 3. **Methods and Techniques**: - Reformulate the system as partial differential - algebraic equations (PDAE) and introduce delay terms for decoupling. - Combine the finite element method in space and the BDF discretization in time to obtain a fully discrete scheme. - Focus on the second - order case, resulting in a three - step scheme. 4. **Main Contributions**: - Prove that the proposed splitting scheme has second - order convergence under the assumption of a weak CFL condition. - Verify the theoretical results through numerical experiments and show the potential of higher - order splitting schemes. 5. **Conclusions**: - The proposed second - order volume - surface splitting scheme shows good convergence and stability when dealing with semi - linear parabolic partial differential equations with dynamic boundary conditions. - Numerical experiments indicate that strict CFL conditions may not be required in practical applications. ### Formula Summary: - **Partial Differential - Algebraic Equations (PDAE)**: \[ \begin{aligned} &\left[\begin{array}{cc} M_u & 0 \\ 0 & M_p \end{array}\right] \left[\begin{array}{c} \dot{u} \\ \dot{p} \end{array}\right] + \left[\begin{array}{cc} K_u & 0 \\ 0 & K_p \end{array}\right] \left[\begin{array}{c} u \\ p \end{array}\right] + \left[\begin{array}{cc} -B_u^T & 0 \\ 0 & B_p^T \end{array}\right] \lambda = \left[\begin{array}{c} f(u) \\ g(p) \end{array}\right], \\ &B_u u - B_p p = 0. \end{aligned} \] - **Three - Step Splitting Scheme**: \[ \begin{aligned} &M_{11} \text{DBDF}_\tau u_1^{n + 3}+K_{11} u_1^{n + 3}=f_1^{n + 3}-M_{12} w^{n + 3}-K_{12} u_2^{n + 3}, \\ &M_{21} \t

A second-order bulk--surface splitting for parabolic problems with dynamic boundary conditions

Splitting schemes for the semi-linear wave equation with dynamic boundary conditions

A Higher Order Unfitted Space-Time Finite Element Method for Coupled Surface-Bulk problems

An almost symmetric Strang splitting scheme for nonlinear evolution equations

Overcoming the order barrier two in splitting methods when applied to semilinear parabolic problems with non-periodic boundary conditions

On second-order and fourth-order elliptic systems consisting of bulk and surface PDEs: Well-posedness, regularity theory and eigenvalue problems

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

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

Compact finite-difference scheme for some Sobolev type equations with Dirichlet boundary conditions

Boundary corrections for splitting methods in the time integration of multidimensional parabolic problems

Parabolic equations with dynamical boundary conditions and source terms on interfaces

A two-stage fourth order time-accurate discretization for Lax--Wendroff type flow solvers II. High order numerical boundary conditions

A Cut Discontinuous Galerkin Method for Coupled Bulk-Surface Problems

Strang splitting structure-preserving high-order compact difference schemes for nonlinear convection diffusion equations

Second-Order Unconditionally Stable Direct Methods for Allen-Cahn and Conservative Allen-Cahn Equations on Surfaces

Additive schemes based domain decomposition algorithm for solving singularly perturbed parabolic reaction-diffusion systems

First- and Second-Order Unconditionally Stable Direct Discretization Methods for Multi-Component Cahn–Hilliard System on Surfaces

Decoupling multistep schemes for elliptic-parabolic problems

Nonlinearly stable flux reconstruction high-order methods in split form

A Fast Compact Exponential Time Differencing Method for Semilinear Parabolic Equations with Neumann Boundary Conditions.

Efficient Parallel Algorithms for Parabolic Problems.