非线性非一致介质二维Maxwell方程leap-frog Crank-Nicolson多区域Legendre-tau配置谱方法

该文研究了具有非线性电导率的非一致介质二维Maxwell方程的数值求解方法,提出了多区域Legendre-tau配置谱方法.该数值方法空间上达到谱精度,时间上是二阶精度.时间方向采用leap-frog Crank-Nicolson三层格式进行离散.非线性项放在中间已知层采用谱配置法显式处理,线性项采用Legendre谱方法隐式处理.利用显隐数值格式,既有较好的稳定性又方便算法实施.基于合理的弱形式,不需要使用额外附加的连接条件,以自然边界条件的方式处理交界面条件.定义不同次数多项式逼近空间,构建一致的数值格式.详细证明了半离散和全离散数值格式的稳定性和收敛性,并得到L2-范数的最优误差估计.算例中,利用快速Legendre变换在Chebyshev点上计算非线性项,提高算法效率.数值结果证实了该数值方法求解此类非线性问题的有效性,并且没有因为解的间断而损失谱精度.