Runge-Kutta型波形松弛方法的A-稳定

摘要：

描述芯片或电力系统运行规律的常用数学模型是高维微分代数方程组,其中的微分方程组太大,诸如线性多步法和Runge-Kutta(RK)法等经典数值方法均不能有效求解.为求解这些微分方程组,学者们提出了波形松弛(WR)方法.多数情况下,这些微分方程组是刚性的,求解他们需要稳定性好的隐式方法,尤其需要A-稳定的方法.此外,RK方法是使用最广泛的常微分方程的数值方法.然而,迄今为止尚未发现RK型WR方法A-稳定的研究.本文研究了 RK型WR方法的A-稳定性,获得了方法A-稳定的充分条件.常见A-稳定的RK方法有Gauss-Legendre方法、Radau ⅠA 方法、Radau ⅡA 方法、Lobatto ⅢA 方法、Lobatto ⅢB 方法、Lobatto ⅢC方法,而且并非RK方法A-稳定,相应的RK型WR方法也A-稳定.在一个假设下,本文所得结果说明当选择低阶Radau ⅠA方法,或Radau ⅡA方法,或Lobatto ⅢC方法为底层方法时,存在分裂方式使得RK型WR方法是A-稳定的.

关键词：A-稳定、波形松弛方法、Runge-Kutta方法、刚性方程

所属期刊栏目：44

分类号：O241.81;O175.1;TM13

资助基金：福建省自然科学基金项目;福建省发树慈善基金

在线出版日期：2023-10-07（万方平台首次上网日期，不代表论文的发表时间）

页数：共10页

页码：327-336