数学物理中一些非线性波方程的谱方法.pdfVIP

数学物理中一些非线性波方程的谱方法 李会元马和平 (上海大学数学系，上海，200436) jl娶本文讨论了数学物理中一些非线性渡方程(组)的谱逼近求解方法，考虑了各种相 应的谱离散格式及其多区域并行算法． 关■谰 非线性波方程，谱方法，谱粘性方法，多区域并行算法 0引言 近年来，偏微分方程的谱方法蓬勃发展，并应用于数学物理中许多问题的数值求解， 受到关注和欢迎．一方面．谱方法具有所谓“无穷阶”的收敛速度，即若原方程的解是无 穷次可徽的，则恰当的半离散谱方法收敛速度比Ⅳ_1的任何次幂都快(这里Ⅳ是所取基函 数的个数)；另一方面，由于在实际计算中一些快速变换(如FFT，FlJ等)的引入大大地 减少了计算量，使得谱方法更具有实用性．本文讨论和推广谱方法在数学物理中一些非线 性波方程数值求解方面的最新发展． 考虑下面的非线性波方程(组) I“(x，o)=Ho(曲， 石EA． 式中： r为正整数；盘，卢≥o；a，=a／a。；a，=a／a：；“=(“l，“2，…，Ⅱ．)7； 性守恒型方程等，在数学物理中有着广泛的应用． 愿建立了F∞fier拟谱方法逼近，时间离散是修改的蛙跳格式：求“。(f)=∑I¨瓯o)c…， 使得 s丁， a：ko)邶a：r+1‰(r)=如)，rE 。斫(卅a，，Ⅳ帅Ⅳ(。)一d “Ⅳ(o)。PⅣ”o， (o-2) 2a；u(o))． 。Ⅳ(f)=PⅣ(。o+zaf”(0)+；r Fouri口拟谱投影算子(插值算子)．我们同时在本文中约定 卟O)=韭兰；≠止生，F(f)=业型；4盟．进而文献【l】证明了格式的广义稳定性，并利用它得到 了式(o．2)的收敛性估计．同时，文献【1】还把此格式推广应用到a2+卢2≠O的情况． 应当看到，Fourier谱和拟谱方法格式比较简单，计算实施方便，因而其研究工作也相 对较为成熟和完善．目前，方程(0．1)的初边值问题和无界区域上的问题收到更多的关注．对 于微分方程的边值问题，应采用代数多项式逼近求解．其关键之处在于如何构造出稳定的、 高精度的谱离散格式．使之实施方便、快捷．为此，一要恰当选取逼近空间及其基函数， 使得全谱离散后的线性方程组系数矩阵稀疏有规律，且条件数较好；二要对非线性项选择 合理的谱投影算子．使之在计算上能尽量采用相应的快速变换，减少计算量；三要保证初 始值的处理同整个谱逼近格式相容，提高方程求解的精度．本文还特别注重各种方法的可 并行性，并力求将所采用的格式推广到高阶的非线性波方程的求解．下文具体给出了类 Burger方程、类xdV方程和Burgers方程外部区域问题的谱方法及其离散解法，井讨论非线 性守恒型方程的谱粘性方法及其后处理过程． 1Bu曙e玮型方程的ch哪shev．L％聊Idm谱方法及其多区域并行算法 考虑下面的类Burgers方程的齐次边界问题 石∈A，o≤t≤r， IatⅣ(耳，f)+a，妒(Ⅱ(工，f))一国：H(x，f)=，(工，f)， {H(一l，f)=“(1，r)=O， 0≤f≤丁， (1．1) I“(J，O)=“o(工)， 膏∈A． 半谱离散逼近格式为：求“。O)∈¨，使得 (12) Vv∈y-． I(HⅣ(o)v)+(a，“Ⅳ(o)，a，y)=(“o，v)+(a：“o，a，V)， 基于引言中提出的几点思想，我们取机(而=c。(k(工)一k+：(算))， 子．若在时间方向采用类似式(0．2)的修改的蛙跳