非线性抛物型方程组的交替偏向有限元格式及其收敛性分析.pdf

山东大学硕士学位论文 厂 中文摘要 f在微分方程的数值计算中，问题的维数往往是制约运算速 度1}寺关键因素，许多常规算法下O(M2d)(其中M为每一维上剖 分的最大点数、d为问题的维数)的计算工作量是不可忍受的， 甚至有些问题有限时间内无法计算．Douglas教授等提出的交 替方向算法大大改进了微分方程的数值计算方法【“，O(M州) 自々计算工作量也实现了许多复杂问题的单机可计算性，而且这 种计算工作量6々改进随着维数越高越明显。 在所有微分方程问题中，线性问题的分析与计算是最简单 的，而且也最有代表性，对非线性问题往往可以在线性问题的 基础上运用一定的技巧及非线性项的特定性质(例如Lipschitz 一层已知量代替的办法来处理，从而使得方程组可裂变为L个 单个自々方程分别求解． 本文共分为五章： 在第一章中，我们提出了如下问题，并叙述了非线性系数 所满足的条件，给出了问题的弱形式． Z X[0，TJ (z，t)∈n ◆ u)m一[V。，A(x，u)V。u]_B(x，u)V。u+f(x，u) Z0)=uo(x) o∈Q 似此“ Zt、=0 (。，t)∈aQ×【0，T】 在第二章中给出了该方程组运用交替方向算法所得到的数 2 值计算格式o j d—l项 d一1项 d一2项 ’(a固旆+揣椭痢(，。Q。商+ 丽2AAt，。啪而前)砌。Q+揣面前吼) c-。G。…。㈨(2订一2矿1)+淼妒 其中： 倦．示方程组中的第f个方程． (砂)≈=((B”(u)V。U“，K。。，…，蛆)l+(，“(∽，K。，氓．，{d)f ，V’、一 1≤七≤叫 一／n【A“(u)Vz旷，VzK”。，’id】f出) 在第三章中对该格式的进行了误差估计，得到了如下误差 估计定理一 ／定理3．1： 当前言中所有的条件都满足，A2mnzh，t)， △；sCh2，定义：en=：Un一札n，并且如果有： IIuo一Ⅳ。IIxo-≤C(h’+1+(△t)2)， l|u1—1彤11I础≤C(h’+1+(A02)， 那么，对任意h存在一个正常数e(u)使e“有如下估计： ∑ 吣+IleⅣlI‰+ 3A至(II篆署晤+¨篙哥峪+ n21 一‘ 山东大学硕士学位论文 ==E=自====!==自==自========一! 接下来在第四章中针对计算格式，从矩阵的角度仔细分析 了算法计算工作量上的优势． 在第五章中针对一个具体的计算实例，运用Fortran语言编 制了计算程序，并列举了部分计算数据． 最后在附录中给出了核心算法的Fortran语言描述． 关键字： ， 非线性抛物型方程组V