多种微分方程数值计算的方法分析.docVIP

多种微分方程数值计算的方法分析 　　 摘要：本文首先从理论上分析了两种主要用于解偏微分方程的数值方法――有限元方法和有限差分法的基本思想和主要的步骤；其次，分别用工程上有很大应用的两类方程进行两种数值方法的数值算例分析；最后，总结两种方法的主要相同点和离散时的重要区别。 　　 关键词：偏微分方程；数值计算；有限元方法；有限差分法 　　 　　 一、引言 　　 在工程应用、科学研究中，所建立的数学模型很大程度上都是偏微分方程。但是，除了极少数特殊类型的偏微分方程能用解析的方法求得其精确解外，大多数情况下要得出解的解析表达式是非常困难的，因此，就需要应用数值方法来近似逼近精确解。 　　 大规模的计算中，往往要用到计算机，然而，电子计算机只能储存有限个数据，做有限次运算，所以任何一种用计算机解题的方法，都必须要把连续问题离散化，最终化成有限形式的线性代数方程组。 　　 二、有限元方法的基本理论 　　 1.基本步骤及理论。有限元法是求解边值问题的数值方法。有限元法求解偏微分方程的基本思想就是传统的Ritz-Galerkin法，但是它运用样条函数提供了一种选取“局部基函数”或“分片多项式空间”的技巧，克服了Ritz-Galerkin法选取基函数的固有困难，它已成为求解偏微分方程，特别是线性椭圆型偏微分方程的一种有效的数值方法。 　　 有限元方法的基本步骤可以归纳如下： 　　 （1）把原问题转化为变分形式。 　　 （2）选定单元形状，对区域进行剖分。 　　 （3）构造节点基函数，形成有限元空间。 　　 （4）以某种方法给出单元各状态变量的离散关系，形成单元刚度矩阵，并且组装成总刚矩阵，有限元法最终导致联立方程组，求解各节点处的函数值。 　　 （5）收敛性及误差估计，对于计算所得的结果，将通过与设计准则提供的准许值比较来评价并确定是否需要重复计算。 　　 2.数值算例。下面以一类二维带复合边界条件的偏微分方程 进行有限元解法的说明。 　　 -αΔu+βu=s，x∈Ωu=d，x∈Dα■+?酌u=r，x∈?祝 （1） 　　 其中，■为u沿单位外法向的方向导数，Ω为有界区域，边界D∪?祝，?鄣Ω，α，β，?酌，r，s，d为给定常数。 　　 本文对区域采取三角剖分，在每个单元（e）上，在该单元上的有限元解可以表示为：u=uiφi+ujφj+ukφk。其中，φi为节点（x，y）的基函数，选取如下： 　　 φi=ai+bix+ciy （2） 　　 满足：φi（xi，yi）=1，φi（xj，yj）=0，φi（xk，yk）=0，将上述条件代入（2），可具体求得基函数的表达式。同理，其他节点处的基函数也是如此形式及求法。 　　 下面将是在变分形式下，推导出方程（1）的有限元解法的代数方程组形式，即方程（1）的弱解。 　　 设u是方程（1）的解，则u在任意单元（e）上也满足方程（1），记单元（e）为Ω，则变分形式如下： 　　 ■（-αΔu+βu）vdσ=■svdσ，?坌v∈H10（Ω） （3） 　　运用Green公式，（3）可以化为： 　　 ■α?塄u??塄v+βuvdσ=■svdσ （4） 　　令u=■uiφi，v=φj，代入上式，整理得： 　　 ■ui■[α（■■+■■）+βφiφj]dσ= 　　 ■sφjdσ （5） 　　其中，对j=1，2，3均成立。写成矩阵的形式即： 　　 （αA+βB）U=b （6） 　　其中，U=（u1，u2，u3）T为待求矩阵， 　　 Aij=■■■+■■dσ （7） 　　 Bij=■φiφjdσ，bj=■sφjdσ （8） 　　 根据上述，对每一个单元（e）都可以求出如式（7），（8）形式的矩阵，即为单刚矩阵，然后将各单元的单刚矩阵组装成总刚矩阵，并根据边界条件删除相关的节点，即可得到最终的线性方程组。 　　 三、有限差分法基本理论 　　 1.基本步骤及理论。有限差分法是微分方程和积分方程数值解的一种重要方法，基本思想是把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组。 　　 有限差分法的基本步骤可以归纳如下： 　　 （1）对求解区域进行网格剖分。 　　 （2）构造逼近微分方程定解问题的差分格式（差分格式有多种，可以根据方程的类型选择适当的差分格式）。 　　 （3）差分方程的解法，差分解的存在唯一性，收敛性及稳定性的研究分析。 　　 2.数值算例。下面以一类一维四阶半线性抛物型方程进行有限差分法的说明。 　　 u1+Δ2u=|u