基于θ方法的波形松弛方法的A稳定

摘要

描述芯片或电力系统运行规律的常用数学模型是高维微分代数方程组,其中的微分方程组太大,线性多步法和Runge-Kutta法等经典数值方法均不能有效求解。为求解这些微分方程组,借鉴常微分方程经典数值方法的A稳定定义,提出了波形松弛方法A稳定(强A稳定),给出了基于θ方法的波形松弛方法A稳定(强A稳定)和非A稳定的条件,以及几个支持理论结果的数值算例。研究结果表明WR方法并非天然继承底层方法的A稳定性,为使波形松弛方法A稳定,需要使用A稳定的底层方法和适当的分裂函数,这为刚性方程WR方法的构造奠定了理论基础。此外,借鉴经典数值方法的B稳定定义,提出了波形松弛方法的B稳定(强B稳定),给出了波形松弛方法强B稳定的条件。