数值计算方法介绍 现代数值计算方法

1、有限差分法(Finite difference method,FDM)

这是求得偏微分方程数值解的最古老方法,对简单的几何形状中的流动与传热问题也是一种最容易实施的方法、其基本的实施方法是,将求解区域用网格线的交点(节点)所组成的点的几何来代替,在每个节点上,描述所研究的流动与传热问题的偏微分方程中的每一个导数项用相应的差分表达式来代替,从而在每个节点上形成一个代数方程,其中包含了本节点及其附近一些节点上的所求量的未知值。求解这些代数方程组就获得了所需的数值解。在规则区域的结构化网格上,有限差分法是十分简单而有效的,而且很容易引入对流项的高阶格式。其不足的是离散方程的守恒特性难以保证,而最严重的缺点则是对不规则区域的适应性差。

2、有限容积法(Finite volume method,FVM)

有限容积法从描述流动与传热问题的守恒型控制方程出发,对它在控制容积上作积分,在积分过程中需要对界面上被求函数的本身(对流通量)及其一阶导数的(扩散通量)构成方式作出假设,这就形成了不同的格式。由于扩散项多是采用相当于二阶精度的线性插值,因而格式的区别主要表现在对流项上。用有限容积法导出的离散方程可以保证具有守恒性(只要界面上的插值方法对位于界面两侧的控制容积式一样的即可),对区域形状的适应性也比有限差分法要好。是目前应用最普遍的一种数值方法。

数值计算方法介绍 现代数值计算方法

3、有限元法(Finite element method,FEM)

有限元法中把计算区域划分成一组离散的容积或者叫元体(在二维情形下元件的形状常常是三角形或者四边形),然后通过对控制方程作积分来得出离散方程,它与有限容积法的主要区别在于:

(1)对每个元体要选定一个形状函数(最简单的为线性函数),通过元件中节点上的被求变量之值来表达该形状函数,并在积分之前把所有假设的形状函数带入到控制方程中去;

(2)控制方程在积分之前应乘上一个选定的权函数,并要求在整个区域上控制方程余量的加权平均值为零,从而导出一组关于节点上被求变量的代数方程。

有限元法的最大优点是,对于不规则几何区域的适应性好。有限元法在对流项的离散处理及不可压缩Naver0-Stokes方程的原始变量法求解方面不如有限容积法发展成熟。但随着有线容积法中非结构化网格的应用,有限容积法与有限差分法之间的差别正在缩小中。

4、有限分析法(Finite analytic method,FAM)

在有限分析法中也像有限差分法、有限容积法那样,用一系列网格线将计算区域进行离散,所不同的是在这里每个节点(网格的交点)与其相邻的四个网格(二维问题)组成一个计算单元,即每一个计算单元由一个内点及八个邻点所组成。在计算单元内把控制方程的非线性项(如N-S方程中的对流项)局部线性化(即认为流速已知),并对该单元边界上未知函数的变化型线作出假设,把所选型线表达式中的常数或者系数用单元边界节点的函数来表示。这样在该单元内的被求解问题转化成了第一类边界条件下的问题,设法找出其分析解,并利用这一分析解找出该单元的内节点及其八个邻点上未知函数值之间的代数表达式,这就是该内点的离散方程。逐一对求解区域的每一个节点建立离散方程,并对计算区域边界上不是第一类条件的每节点补充一个方程,就完成了整个计算区域的离散方程的建立过程。

有限分析法是80年代初发展起来的一种数值方法,它可以克服高Reynolds数下有限差分法及有限容积法的数值解容易发散或振荡的缺点,但其计算工作量较大,对计算区域几何形状的适应性也较差。

5、边界元法(Boudary element method,BEM)

上面四种方法都必须对整个区域作离散化处理,用分布在整个区域上的有限个节点上函数的近似值来代替连续问题的解。在边界元方法中应用格林函数公式,并通过选择适当的权函数把空间求解域上的偏微分方程转换成为其边界上的积分方程,它把求解区中任一点的求解变量(如温度)与边界条件联系了起来,通过离散化处理,由积分方程导出边界节点上未知值的代数方程。解出边界上的未知值后就可以利用边界积分方程来获得内部任一点的被求函数之值。边界元法的最大优点是,可以使求解问题的空间维数降低一阶,从而使计算工作量及其所需计算机容量大大减小,边界元法推广应用的一个最大限制是,需要已知所求解偏微分方程的格林函数基本解,虽然对不少偏微分方程这个基本解业已找出,但对Navier-Stokes方程这个非线性偏微分方程,至今尚未找到基本解。目前的一种处理方式是,把Navier-Stokes方程中的非线性项看作是扩散方程的源项并通过迭代的方式来求解,但一般只能获得Re较低情形的解。

6、谱分析法(Spectral method,SM)

在谱方法中被求解的函数用有限项的级数展开来表示。例如,有限项的傅里叶展开,多项式展开等。显然,只要这种级数中每一项的系数一确定,这个被求解的函数的近似形式也就完全确定了。因而与前面五种离散方法不同的是,在谱方法中要建立的代数方程式关于这些系数的代数方程,而不是节点上被求解函数值的代数方程。显然,级数的项越多,其精度也越高。建立未知系数的代数方程的基本方法是加权余数法,即首先将近似解带入控制方程(设控制方程中所有的项均移到了等号左边),再乘以近似解级数中的一个项称为权函数(不包括该项中待定的系数),然后对整个求解区域作积分,并要求该积分式等于零,就得出一个关于待定系数的代数方程,这样以系数解中每一个含有待定系数的项作权函数,就可以得到总数与待定系数数量相等的代数方程组。求解方程组,就得出了被求解函数的近似解。

谱分析法用于偏微分方程的近似求解始于20世纪70年代末,其优点是可以获得高精度的解,但不适宜用于编制通用程序,目前只是在比较简单的流动与传热问题中已经应用的比较成功。

7、数值积分变换法(Integral transformation method,ITM)

我们知道在偏微分方程的分离变量法求解中,一个被求函数表示了两个特征值问题所规定的特征函数乘积的线性组合。一般地说,任意一个函数可用相应的特征值问题的特征函数来表示。对不具备分析解的非线性偏微分方程,设法把它的解表示成一个特征值问题的解及一个降维的定解问题解的组合:其中特征值问题具有分析解,而定解问题则应包含该问题的诸多非线性复杂因素,因而要采用数值解法,这就是数值积分变换法。由于该定解问题降低了维数(常常变为常微分方程的初值问题或者两点边值问题),因而数值求解比较方便,这一方法是分析解法与数值解法的混合方法。对于数值解问题而言,只需要得出某个降维了的空间内(对二维问题或者是时间坐标或者是某个空间坐标)的数值解,整个求解区域内的值则可以用分析解与数值解的组合来得到。在获得降维的定解问题过程中需要进行积分变换与反变换,而该降维问题的求解是采用数值方法的,因而称其为数值积分变换法。这一方法的优点是计算精度可以较高,降维问题是一个常微分方程,有成熟的数值方法可以采用,但这种方法不容易形成通用程序,特征值问题的选取有一定的任意性,对非线性强烈的问题,计算工作量比较大。

8、格子-Boltzmann方法(Lattice-Boltzmann method,LBM)

格子-Boltzmann方法是基于分子运动论的一种模拟流体流动的数值方法。在上述各种数值方法中,把本质上市离散的介质先假定是连续的,在此基础上建立起了N-S方法,然后又再把它离散化。在LMB中不再基于连续介质的假设,而是把流体看成是许多只有质量没有体积的微粒所组成,这些微粒可以向空间若干个方向任意运动。通过其质量、动量守恒的原理,建立起表征质点在给定的时刻位于空间某一位置附近的概率密度函数。再通过统计的方法来获得质点微粒的概率密度分布函数与宏观运动参数间的关系。这方法从提出至今只有10余年的历史,但已显示出其巨大的发展潜力。

此外应用于流动与换热计算的数值方法还有控制容积有限元法(CVFEM)及微分求积法(DQM)等。

-----------------摘自《计算传热学的近代发展》.陶文权

  

爱华网本文地址 » http://www.413yy.cn/a/25101015/247304.html

更多阅读

年底双薪计算方法 年底双薪什么时候发

年底双薪计算方法——简介年底双薪,简单的理解就是你到某一个公司(工厂)上班,到了年底的时候,给你发两个月的工资,即一年可以领13个月的工资,但要符合其发放条件,下面详细介绍发放条件及计算方法年底双薪计算方法——工具/原料电

宝宝身高体重标准表2014 孩子的身高计算方法

现将宝宝身高体重标准计算公式方法介绍如下:宝宝体重标准测量公式1~6个月时体重(公斤)=出生体重(或3公斤)+月龄×0.6(公斤);7~12月体重(公斤)=出生体重(或3公斤)+月龄×0.5(公斤);2~10岁体重(公斤)=年龄×2+7(或8)。宝宝身高标准的增

女性生理周期计算方法|女性生理周期计算器 女性生理计算器

女性生理周期计算方法是每一个女性朋友都想知道的,那你想知道女性生理周期计算方法到底是怎么计算的呢?请看下面的内容介绍。女性生理周期计算方法:女性的排卵日期一般在下次月经来潮前的14天左右。下次月经来潮的第1天算起,倒数14天或

声明:《数值计算方法介绍 现代数值计算方法》为网友灰色格调分享!如侵犯到您的合法权益请联系我们删除