数值模拟偏微分方程的三种方法

《数值模拟偏微分方程的三种方法》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、数值模拟偏微分方程的三种方法I.三者简介有限差分方法（FiniteDifferenceMethods）是数值模拟偏微分方程最早采用的方法，至今仍被广泛使用。该方法包括区域剖分和差商代替导数两个步骤。首先将求解区域划分为弟分网格，用有限个网格节点代替连续的求解区域。其次，利用Taylor级数展开等方法将偏微分方程屮的导数项在网格节点上用函数值的差商代替进行离散，从而建立以网格节点上的值为未知量的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且十分成熟的数值方法。茅商代替导数

2、后的格式称为有限茅分格式，从格式的精度来考虑，有一阶格式、二阶格式和高阶格式。从差分的空间离散形式来考虑，冇屮心格式和迎风格式。对于瞬态方程，考虑时间方向的离散，有显格式、隐格式、交替显隐格式等。口前常见的差分格式，主要是以上几种格式的组合，不同的组合构成不同的差分格式。差分方法主要适用于结构网格，网格的人小一•般根据问题模型和Courant稳定条件來决定。有限元方法（FiniteElementMethods）的基础是虚位移原理和分片多项式插值。该方法的构造过程包括以下三个步骤。首先，利用虚位移原理得到偏微分方程的弱形式，将

3、计算区域划分为有限个互不重叠的单元（三角形、以边形、四面体、六面体等），在每个单元上选择合适的节点作为求解函数的插值点，将偏微分方程屮的变量改写成由各变量或其导数的节点值与所选用的分片插值基函数组成的线性表达式，得到微分方程的离散形式。利用插值函数的局部支集性质及数值积分可以得到未知量的代数方程组。有限元方法有较完善的理论基础，具有求解区域灵活（复杂区域）、单元类型灵活（适于结构网格和非结构网格）、程序代码通用（数值模拟软件多数基于冇限元方法）等特点。有限元方法最早应用于结构力学，随着计算机的发展已经渗透到计算物理、流体力学

4、与电磁学等各个数值模拟领域。根据所采用的检验函数（虚位移函数）和插值函数的不同，有限元方法也分为多种计算格式。从检验函数的选择来说，有配置法、最小二乘法和伽辽金法，从计算单元网格的形状來划分，有三角形网格、四边形网格和多面体网格等，从插值函数的精度來划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于有限元方法，其基木思路和解题步骤可归纳为（1）建立积分方程，根据虚位移原理或方程余量，建立与微分方程初边值问题等价的积分表达式，这是有限元法的出发点。(2)区域单元剖分，根据求解区域的形状及实际

5、问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。区域单元划分釆用冇限元方法的前处理完成,并给出计算单元和节点编号相互之间的关系、节点的位置坐标，同时还需要列出问题的边界的节点号和相应的边值条件。(3)确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函数作为单元的形函数。有限兀方法屮的形函数是在单兀屮选取的，由于各单元具有规则的儿何形状，在选取形函数时可遵循一定的法则。(4)单元分析：将各个单元小的求解函数用单元形函数的线性组合表达式进行逼近;再将近似函数代入积分方程，并对单元区域进行

6、积分，可获得含冇待定系数(即单元中各节点的函数值)的单元矩阵与荷载。(5)总体合成：在得岀单元矩阵与荷载Z后，将区威中所冇单元矩阵与荷载按一定法则进行迭加，形成总体有限元方程。(6)边界条件的处理：一般边界条件有三种形式，分为木质边界条件(Dirichlet边界条件)、自然边界条件(Neumann边界条件)、混合边界条件(Cauchy边界条件)。对于自然边界条件，一般在积分表达式屮可自动得到满足。对于本质边界条件和混合边界条件，需按一定法则后对总休有限元方程进行修止。(7)解有限元方程：根据边界条件修正的总体有限元方程组，采

7、用适当的代数方程组求解器，求出各节点的函数值。有限体积法(FiniteVolumeMethod)又称为控制体积法。其基本思路是：将计算区域划分为一系列互不重叠的控制体，并使每个网格点周围冇一个控制体；将待求解的微分方程对每一个控制体积积分，便得岀一组离散方程。该方法的未知量为网格点上的函数值。为了求出控制体积的积分，须假定函数值在网格点控制体边界上的变化规律。从积分区域的选取方法來看，冇限体积法属于冇限元方法中检验函数取分片常数插值的了区域法;从未知量的近似方法看來，有限体积法屈于采用局部近似多项式插值逼近。有限体积法的基木

8、思路易于理解，能够保持物理量在控制体上的守恒性质，也即离散方程保持了微分方程物理量在控制体满足某种守恒原理的物理意义。这是有限体积法吸引人的优点。此外，在有限体积法中，插值函数只用于计算控制体积的积分，因此可以对微分方程小不同的项采取不同的插值函数。II.三者各有长短有限差分方法直观，经验