Abstract:In contrast with the diffusion equation which smoothens the initial data to C ∞ C^\infty for t > 0 t>0 (away from the corners/edges of the domain), the subdiffusion equation only exhibits limited spatial regularity. As a result, one generally cannot expect high-order accuracy in space in solving the subdiffusion equation with nonsmooth initial data. In this paper, a new splitting of the solution is constructed for high-order finite element approximations to the subdiffusion equation with nonsmooth initial data. The method is constructed by splitting the solution into two parts, i.e., a time-dependent smooth part and a time-independent nonsmooth part, and then approximating the two parts via different strategies. The time-dependent smooth part is approximated by using high-order finite element method in space and convolution quadrature in time, while the steady nonsmooth part could be approximated by using smaller mesh size or other methods that could yield high-order accuracy. Several examples are presented to show how to accurately approximate the steady nonsmooth part, including piecewise smooth initial data, Dirac–Delta point initial data, and Dirac measure concentrated on an interface. The argument could be directly extended to subdiffusion equations with nonsmooth source data. Extensive numerical experiments are presented to support the theoretical analysis and to illustrate the performance of the proposed high-order splitting finite element methods.

An adaptive low-rank splitting approach for the extended Fisher–Kolmogorov equation

Full-rank and low-rank splitting methods for the Swift–Hohenberg equation

Finite element analysis of extended Fisher-Kolmogorov equation with Neumann boundary conditions

A low-rank Lie-Trotter splitting approach for nonlinear fractional complex Ginzburg-Landau equations

A stable multiplicative dynamical low-rank discretization for the linear Boltzmann-BGK equation

Krylov-based Adaptive-Rank Implicit Time Integrators for Stiff Problems with Application to Nonlinear Fokker-Planck Kinetic Models

A multi-fidelity adaptive dynamical low-rank based optimization algorithm for fission criticality problems

Projected exponential methods for stiff dynamical low-rank approximation problems

High-order Adaptive Rank Integrators for Multi-scale Linear Kinetic Transport Equations in the Hierarchical Tucker Format

Asymptotic-Preserving and Energy Stable Dynamical Low-Rank Approximation

High-order splitting finite element methods for the subdiffusion equation with limited smoothing property

A low-rank algorithm for strongly damped wave equations with visco-elastic damping and mass terms

Accelerating the simulation of kinetic shear Alfvén waves with a dynamical low-rank approximation

Partially explicit splitting scheme with explicit-implicit-null method for nonlinear multiscale flow problems

Interpolatory dynamical low-rank approximation for the 3+3d Boltzmann-BGK equation

Adaptive Low-Rank Methods: Problems on Sobolev Spaces

A space-time adaptive low-rank method for high-dimensional parabolic partial differential equations

A rank-adaptive robust integrator for dynamical low-rank approximation

A new 2-level implicit high accuracy compact exponential approximation for the numerical solution of nonlinear fourth order Kuramoto–Sivashinsky and Fisher–Kolmogorov equations