论文标题
对哈密顿PDES的任意高阶线性隐式传播算法
Arbitrary high-order linearly implicit energy-preserving algorithms for Hamiltonian PDEs
论文作者
论文摘要
在本文中,我们提出了一种新的策略,以系统地构建具有任意准确性的汉密尔顿PDE的线性隐含能量的方案。这种新颖的策略基于新开发的指数标量变量(ESAV)方法,该方法可以消除哈密顿功能中非线性项的边界限制,并提供对辅助变量的完全显式离散化,而无需计算额外的内部产品,这使得它比传统的标量auxil auxil auxil auxil auxil auxil auxil auxil auxil auxial auxial auxil auxil auxil auxil auxil auxial var(save sav)。为了实现任意的高阶精度和能量保存,我们为解决方案变量和辅助变量利用了符号runge-kutta方法,其中非线性项的内部阶段值是通过在先前计算中已经从数值的溶液中明确得出的。提出了一种预测校正策略,以进一步提高准确性。然后使用傅立叶伪谱法来获得完全离散的方案。与SAV方案相比,这些ESAV方案中的解决方案变量和辅助变量现在被解耦。此外,当线性项为恒定系数时,可以使用快速傅立叶变换来明确求解解决方案变量。对三个哈密顿PDE进行了数值实验,以证明ESAV方案的效率和保护。
In this paper, we present a novel strategy to systematically construct linearly implicit energy-preserving schemes with arbitrary order of accuracy for Hamiltonian PDEs. Such novel strategy is based on the newly developed exponential scalar variable (ESAV) approach that can remove the bounded-from-blew restriction of nonlinear terms in the Hamiltonian functional and provides a totally explicit discretization of the auxiliary variable without computing extra inner products, which make it more effective and applicable than the traditional scalar auxiliary variable (SAV) approach. To achieve arbitrary high-order accuracy and energy preservation, we utilize the symplectic Runge-Kutta method for both solution variables and the auxiliary variable, where the values of internal stages in nonlinear terms are explicitly derived via an extrapolation from numerical solutions already obtained in the preceding calculation. A prediction-correction strategy is proposed to further improve the accuracy. Fourier pseudo-spectral method is then employed to obtain fully discrete schemes. Compared with the SAV schemes, the solution variables and the auxiliary variable in these ESAV schemes are now decoupled. Moreover, when the linear terms are of constant coefficients, the solution variables can be explicitly solved by using the fast Fourier transform. Numerical experiments are carried out for three Hamiltonian PDEs to demonstrate the efficiency and conservation of the ESAV schemes.