English

Higher-order generalized-$\alpha$ methods for parabolic problems

Numerical Analysis 2021-02-12 v1 Numerical Analysis

Abstract

We propose a new class of high-order time-marching schemes with dissipation user-control and unconditional stability for parabolic equations. High-order time integrators can deliver the optimal performance of highly-accurate and robust spatial discretizations such as isogeometric analysis. The generalized-α\alpha method delivers unconditional stability and second-order accuracy in time and controls the numerical dissipation in the discrete spectrum's high-frequency region. Our goal is to extend the generalized-alphaalpha methodology to obtain a high-order time marching methods with high accuracy and dissipation in the discrete high-frequency range. Furthermore, we maintain the stability region of the original, second-order generalized-alphaalpha method foe the new higher-order methods. That is, we increase the accuracy of the generalized-α\alpha method while keeping the unconditional stability and user-control features on the high-frequency numerical dissipation. The methodology solve k>1,kNk>1, k\in \mathbb{N} matrix problems and updates the system unknowns, which correspond to higher-order terms in Taylor expansions to obtain (3/2k)th(3/2k)^{th}-order method for even kk and (3/2k+1/2)th(3/2k+1/2)^{th}-order for odd kk. A single parameter ρ\rho^\infty controls the dissipation, and the update procedure follows the formulation of the original second-order method. Additionally, we show that our method is A-stable and setting ρ=0\rho^\infty=0 allows us to obtain an L-stable method. Lastly, we extend this strategy to analyze the accuracy order of a generic method.

Keywords

Cite

@article{arxiv.2102.05910,
  title  = {Higher-order generalized-$\alpha$ methods for parabolic problems},
  author = {Pouria Behnoudfar and Quanling Deng and Victor M. Calo},
  journal= {arXiv preprint arXiv:2102.05910},
  year   = {2021}
}