Abstract:We describe novel numerical methods for solving a class of high-order time-dependent PDEs on general geometries, which involve second-order derivatives in time and up-to fourth-order derivatives in space. This type of PDEs are widely used in applications such as the Boussinesq equation and in modeling thin-walled elastic structures such as beams, plates and shells, etc. High-order spatial derivatives together with general geometries bring a number of challenges for many numerical methods. In this paper, we resolve these challenges by discretizing the spatial derivatives in domains with general geometries using the composite overlapping grid method. The discretization on overlapping grids requires numerical interpolations to couple solutions on different component grids. However, the presence of interpolation equations breaks the symmetry of the overall spatial discretization, causing numerical instability in time-stepping schemes. To address this, a fourth-order hyper-dissipation term is included for stabilization. Investigation of incorporating the hyper-dissipation term into several time-stepping schemes for solving the semi-discrete system leads to the development of a series of algorithms. Accurate and stable numerical boundary conditions for Dirichlet and Neumann type boundaries are also developed for general geometries. Quadratic eigenvalue problems for a simplified model problem on 1D overlapping grids are considered to reveal the weak instability caused by interpolation between component grids. This model problem is also investigated for the stabilization effects of the proposed algorithms. Carefully designed numerical experiments and two benchmark problems concerning the Kirchhoff–Love plate model are presented to demonstrate the accuracy and efficiency of our approaches. This work shows that finite difference methods on overlapping grids are well-suited for solving high-order PDEs in complex domains for realistic applications.

Fourth-order Compact Block-Centered Splitting Domain Decomposition Method for Parabolic Equations

Decomposition Methods for Time‐domain Maxwell's Equations

Efficient Parallel Algorithms for Parabolic Problems.

A Five-Step Block Method Coupled with Symmetric Compact Finite Difference Scheme for Solving Time-Dependent Partial Differential Equations

Optimized Operator-Splitting Methods in Numerical Integration of Maxwell's Equations

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

A Fourth-Order Accurate Finite-Volume Method with Structured Adaptive Mesh Refinement for Solving the Advection-Diffusion Equation

Numerical solution of 2-D incompressible viscous flow with high order compact scheme using domain decomposition and matched method

Splitting Algorithms for the High-Order Compact Finite-Difference Schemes in Wave-Equation Modeling

A modified domain decomposition spectral collocation method for parabolic partial differential equations

A fourth-order compact solver for fractional-in-time fourth-order diffusion equations

Numerical Methods for Fourth-Order PDEs on Overlapping Grids with Application to Kirchhoff–Love Plates

Parallel D–D Type Domain Decomposition Algorithm for Optimal Control Problem Governed by Parabolic Partial Differential Equation

A Compact Exponential Scheme for Solving 1D Unsteady Convection-Diffusion Equation with Neumann Boundary Conditions

Fourth-order compact finite difference schemes for solving biharmonic equations with Dirichlet boundary conditions

Numerical solution of time dependent nonlinear partial differential equations using a novel block method coupled with compact finite difference schemes

A parallel domain decomposition method for solving elliptic equations on manifolds

A domain decomposing reduced-order method for solving the multi-domain transient heat conduction problems

A Fourth-Order Compact Finite Volume Method on Unstructured Grids for Simulation of Two Dimensional Incompressible Flow

A compression scheme for domain decomposition method in solving electromagnetic problems

High order compact augmented methods for Stokes equations with different boundary conditions