simpop
普通会员
流体力学数值方法有很多种,其数学原理各不相同,但有二点是所有方法都具备的,即离散化和代数化。总的来说其基本思想是:将原来连续的求解区域划分成网格或单元子区域,在其中设置有限个离散点(称为节点),将求解区域中的连续函数离散为这些节点上的函数值;通过某种数学原理,将作为控制方程的偏微分方程转化为联系节点上待求函数值之间关系的代数方程(离散方程),求解所建立起来的代表方程以获得求解函数的节点值。
有限差分法 +
这是最早采用的数值方法,它是将求解区域划分为矩形或正交曲线网格,在网格线交点(即节点)上,将控制方程中的每一个微商用差商来代替,从而将连续函数的微分方程离散为网格节点上定义的差分方程,每个方程中包含了本节点及其附近一些节点上的待求函数值,通过求解这些代数方程就可获得所需的数值解。
有限元法 +
有限元法的基本原理是把适定的微分问题的解域进行离散化,将其剖分成相连结又互不重叠的具有一定规则几何形状的有限个子区域(如:在二维问题中可以划分为三角形或四边形;在三维问题中可以划分为四面体或六面体等),这些子区域称之为单元,单元之间以节点相联结。函数值被定义在节点上,在单元中选择基函数(又称插值函数),以节点函数值与基函数的乘积的线性组合成单元的近似解来逼近单元中的真解。
边界元法 +
边界元法是在经典积分方程和有限元法基础上发展起来的求解微分方程的数值方法,其基本思想是:将微分方程相应的基本解作为权函数,应用加权余量法并应用格林函数导出联系解域中待求函数值与边界上的函数值与法向导数值之间关系的积分方程;令积分方程在边界上成立,获得边界积分方程,该方程表述了函数值和法向导数值在边界上的积分关系,而在这些边界值中,一部份是在边界条件中给定的,另一部份是待求的未知量,边界元法就是以边界积分方程作为求解的出发点,求出边界上的未知量;在所导出的边界积分方程基础上利用有限元的离散化思想,把边界离散化,建立边界元代数方程组,求解后可获得边界上全部节点的函数值和法向导数值;将全部边界值代入积分方程中,即可获得内点函数值的计算表达式,它可以表示成边界节点值的线性组合。
有限体积法 +
有限体积法又称为控制体积法,其导出离散方程的基本思路是:
有限分析法 +
有限分析法在某种意义上说是在有限元法基础上发展起来的一种数值方法,其基本思想是:将求解区域划分成矩形网格,网格线的交点为计算节点,每个节点与相邻的四个网格组成一个计算单元,即一个计算单元由一个中心节点与8个相邻节点组成;在每个单元中函数的近似解不是象有限元方法那样采用单元基函数的线性组合来表达,而是以单元中未知函数的分析解来表达;为了获得单元中的分析解,单元边界条件采用插值函数来逼近,在单元中把控制方程中非线性项局部线性化(如N-S方程中的对流项中认为其流速为已知,并对单元中待求函数的组合形式作出假设,找出其系数用单元边界节点上待求函数值表达的分析解;利用单元分析解确定单元中心节点与8个相邻节点间待求函数值之间关系的一个代数方程,称为单元有限分析方程;将所有内点上的单元有限分析方程联立,就构成总体有限分析方程,通过代数方程组求解,即可获得求解区域中全部离散点的函数值。
本文摘录自华南理工大学交通学院李志印、熊小辉、吴家鸣撰写的《计算流体力学常用数值方法简介》一文,该文发表于《广东造船》2004年第3期。