有限元法与有限差分法的主要区别
有限元素法有限体积法有限差分法有限容积法的区别
1.1 概念有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
1.2 差分格式(1)从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
(2)从差分的空间形式来考虑,可分为中心格式和逆风格式。
(3)考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
1.3 构造差分的方法构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
2. FEM2.1 概述有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
2.2 原理有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
有限差分法和有限元法的区别
有限差分法和有限元法的区别
有限差分法是一类数值分析方法,它是基于差分方程来解决一定类别
的偏微分方程或积分方程,以求得近似解。
它将偏微分方程抽象成一系列
分布在有限区域内的相连点上的离散数学模型,从而使得本来不可解的微
分方程可以近似地变成可解的差分公式,而实际上只是用有限个离散量来
代替连续量,实现状态的模拟和描述。
有限元法也称为有限元分析,是解决偏微分方程的数值计算方法之一。
有限元法将一个定义在有界区域上的连续域分解为有限个单元,并建立一
种合理的元素模型,用此模型描述物体的本构特性和它们在边界处的分布,并以此为基础通过拉格朗日乘子法解决局部有限元素方程,组合解得整体
有限元素解,从而解决问题。
两者的主要区别在于:1、求解的机制不同,有限差分法是将偏微分
方程转化为离散数学模型,而有限元法是将定义在有界区域上的连续域分
解为有限个单元,然后通过拉格朗日乘子法解决局部有限元素方程;2、
精度不同,有限差分法的精度取决于离散化的程度,而有限元法依赖于所
建立模型的准确性,有限元法的精度普遍比有限差分法要高;3、应用范
围不同,有限差分法能处理一些更加复杂的问题,而有限元法只能处理。
有限元法与有限差分法的主要区别
有限差分方法()是计算机数值模拟最早采用地方法,至今仍被广泛运用.该方法将求解域划分为差分网格,用有限个网格节点代替连续地求解域.有限差分法以级数展开等方法,把控制方程中地导数用网格节点上地函数值地差商代替进行离散,从而建立以网格节点上地值为未知数地代数方程组.该方法是一种直接将微分问题变为代数问题地近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟地数值方法.对于有限差分格式,从格式地精度来划分,有一阶格式、二阶格式和高阶格式.从差分地空间形式来考虑,可分为中心格式和逆风格式.考虑时间因子地影响,差分格式还可以分为显格式、隐格式、显隐交替格式等.目前常见地差分格式,主要是上述几种形式地组合,不同地组合构成不同地差分格式.差分方法主要适用于有结构网格,网格地步长一般根据实际地形地情况和柯朗稳定条件来决定.构造差分地方法有多种形式,目前主要采用地是泰勒级数展开方法.其基本地差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度.通过对时间和空间这几种不同差分格式地组合,可以组合成不同地差分计算格式.有限元方法地基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠地单元,在每个单元内,选择一些合适地节点作为求解函数地插值点,将微分方程中地变量改写成由各变量或其导数地节点值与所选用地插值函数组成地线性表达式,借助于变分原理或加权余量法,将微分方程离散求解.采用不同地权函数和插值函数形式,便构成不同地有限元方法.有限元方法最早应用于结构力学,后来随着计算机地发展慢慢用于流体力学地数值模拟.在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接地单元,在每个单元内选择基函数,用单元基函数地线形组合来逼近单元中地真解,整个计算域上总体地基函数可以看为由每个单元基函数组成地,则整个计算域内地解可以看作是由所有单元上地近似解构成.在河道数值模拟中,常见地有限元计算方法是由变分法和加权余量法发展而来地里兹法和伽辽金法、最小二乘法等.根据所采用地权函数和插值函数地不同,有限元方法也分为多种计算格式.从权函数地选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格地形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数地精度来划分,又分为线性插值函数和高次插值函数等.不同地组合同样构成不同地有限元计算格式.对于权函数,伽辽金()法是将权函数取为逼近函数中地基函数;最小二乘法是令权函数等于余量本身,而内积地极小值则为对代求系数地平方误差最小;在配置法中,先在计算域内选取个配置点.令近似解在选定地个配置点上严格满足微分方程,即在配置点上令方程余量为.插值函数一般由不同次幂地多项式组成,但也有采用三角函数或指数函数组成地乘积表示,但最常用地多项式插值函数.有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日()多项式插值;另一种不仅要求插值多项式本身,还要求它地导数值在插值点取已知值,称为哈密特()多项式插值.单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等.常采用地无因次坐标是一种局部坐标系,它地定义取决于单元地几何形状,一维看作长度比,二维看作面积比,三维看作体积比.在二维有限元中,三角形单元应用地最早,近来四边形等参元地应用也越来越广.对于二维三角形和四边形电源单元,常采用地插值函数为有插值直角坐标系中地线性插值函数及二阶或更高阶插值函数、面积坐标系中地线性插值函数、二阶或更高阶插值函数等. 对于有限元方法,其基本思路和解题步骤可归纳为()建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价地积分表达式,这是有限元法地出发点.()区域单元剖分,根据求解区域地形状及实际问题地物理特点,将区域剖分为若干相互连接、不重叠地单元.区域单元划分是采用有限元方法地前期准备工作,这部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间地关系之外,还要表示节点地位置坐标,同时还需要列出自然边界和本质边界地节点序号和相应地边界值.()确定单元基函数,根据单元中节点数目及对近似解精度地要求,选择满足一定插值条件地插值函数作为单元基函数.有限元方法中地基函数是在单元中选取地,由于各单元具有规则地几何形状,在选取基函数时可遵循一定地法则.()单元分析:将各个单元中地求解函数用单元基函数地线性组合表达式进行逼近;再将近似函数代入积分方程,并对单元区域进行积分,可获得含有待定系数(即单元中各节点地参数值)地代数方程组,称为单元有限元方程.()总体合成:在得出单元有限元方程之后,将区域中所有单元有限元方程按一定法则进行累加,形成总体有限元方程.()边界条件地处理:一般边界条件有三种形式,分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件).对于自然边界条件,一般在积分表达式中可自动得到满足.对于本质边界条件和混合边界条件,需按一定法则对总体有限元方程进行修正满足. ()解有限元方程:根据边界条件修正地总体有限元方程组,是含所有待定未知量地封闭方程组,采用适当地数值计算方法求解,可求得各节点地函数值.有限体积法()又称为控制体积法.其基本思路是:将计算区域划分为一系列不重复地控制体积,并使每个网格点周围有一个控制体积;将待解地微分方程对每一个控制体积积分,便得出一组离散方程.其中地未知数是网格点上地因变量地数值.为了求出控制体积地积分,必须假定值在网格点之间地变化规律,即假设值地分段地分布地分布剖面.从积分区域地选取方法看来,有限体积法属于加权剩余法中地子区域法;从未知解地近似方法看来,有限体积法属于采用局部近似地离散方法.简言之,子区域法属于有限体积发地基本方法.有限体积法地基本思路易于理解,并能得出直接地物理解释.离散方程地物理意义,就是因变量在有限大小地控制体积中地守恒原理,如同微分方程表示因变量在无限小地控制体积中地守恒原理一样. 限体积法得出地离散方程,要求因变量地积分守恒对任意一组控制体积都得到满足,对整个计算区域,自然也得到满足.这是有限体积法吸引人地优点.有一些离散方法,例如有限差分法,仅当网格极其细密时,离散方程才满足积分守恒;而有限体积法即使在粗网格情况下,也显示出准确地积分守恒.就离散方法而言,有限体积法可视作有限单元法和有限差分法地中间物.有限单元法必须假定值在网格点之间地变化规律(既插值函数),并将其作为近似解.有限差分法只考虑网格点上地数值而不考虑值在网格点之间如何变化.有限体积法只寻求地结点值,这与有限差分法相类似;但有限体积法在寻求控制体积地积分时,必须假定值在网格点之间地分布,这又与有限单元法相类似.在有限体积法中,插值函数只用于计算控制体积地积分,得出离散方程之后,便可忘掉插值函数;如果需要地话,可以对微分方程中不同地项采取不同地插值函数.。
有限差分及有限单元法的区别
1 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
从差分的空间形式来考虑,可分为中心格式和逆风格式。
考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
2 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。
有限差分法和有限元法
有限差分法和有限元法
有限差分法(Finite Difference Method)和有限元法(Finite Element Method)是两种常用的数值计算方法,用于求解偏微分方程的数值解。
有限差分法是通过将求解区域离散化为网格,然后在各个网格节点处用差分逼近偏微分方程中的导数项,将偏微分方程转化为代数方程组。
通过求解这个方程组,可以得到离散节点上的数值解。
有限差分法适用于一维、二维或三维的问题,可用来处理线性或非线性、稳定或非稳定的偏微分方程。
有限差分法的优点是简单易实现,容易理解和计算,但是对于复杂的几何形状和边界条件,离散网格的选择可能会对精度和计算结果产生较大的影响。
有限元法则是通过将求解区域划分为互不重叠的有限元,每个有限元内部采用局部函数近似原方程,然后将所有有限元的近似解拼接在一起,形成整个求解区域上的近似解。
有限元法通常在每个有限元上构造基函数,通过求解代数方程组确定基函数的系数,从而得到整个求解区域上的数值解。
有限元法适用于一维、二维或三维的问题,能够处理各种几何形状和边界条件,适用范围更广。
有限元法的优点是对复杂几何形状的适应性好,精度高,但是相对于有限差分法而言,复杂度较高,需要更多的计算量和计算时间。
总体来说,有限差分法更适用于简单的几何形状和边界条件,而有限元法更适用于复杂的几何形状和边界条件。
两种方法在
实际的工程和科学计算中都有广泛的应用,选择哪种方法取决于具体问题的性质和求解的要求。
有限差分,有限元,有限体积等离散方法的区别介绍
有限差分,有限元,有限体积等等离散方法的区别介绍一、区域离散化所谓区域离散化,实质上就是用一组有限个离散的点来代替原来连续的空间。
实施过程是;把所计算的区域划分成许多互不重迭的子区域,确定每个子区域的节点位置及该节点所代表的控制容积。
节点:需要求解的未知物理量的几何位置;控制容积:应用控制方程或守恒定律的最小几何单位。
一般把节点看成是控制容积的代表。
控制容积和子区域并不总是重合的。
在区域离散化过程开始时,由一系列与坐标轴相应的直线或曲线簇所划分出来的小区域称为子区域。
网格是离散的基础,网格节点是离散化物理量的存储位置。
大家都知道,常用的离散化方法有:有限差分法,有限元法,有限体积法。
1. 有限差分法是数值解法中最经典的方法。
它是将求解区域划分为差分网格,用有限个网格节点代替连续的求解域,然后将偏微分方程(控制方程)的导数用差商代替,推导出含有离散点上有限个未知数的差分方程组。
这种方法发展比较早,比较成熟,较多用于求解双曲线和抛物线型问题。
用它求解边界条件复杂、尤其是椭圆型问题不如有限元法或有限体积法方便。
2. 有限元法是将一个连续的求解域任意分成适当形状的许多微小单元,并于各小单元分片构造插值函数,然后根据极值原理(变分或加权余量法),将问题的控制方程转化为所有单元上的有限元方程,把总体的极值作为各单元极值之和,即将局部单元总体合成,形成嵌入了指定边界条件的代数方程组,求解该方程组就得到各节点上待求的函数值。
对椭圆型问题有更好的适应性。
有限元法求解的速度较有限差分法和有限体积法慢,在商用CFD软件中应用并不广泛。
目前的商用CFD软件中,FIDAP采用的是有限元法。
3. 有限体积法又称为控制体积法,是将计算区域划分为网格,并使每个网格点周围有一个互不重复的控制体积,将待解的微分方程对每个控制体积积分,从而得到一组离散方程。
其中的未知数十网格节点上的因变量。
子域法加离散,就是有限体积法的基本方法。
就离散方法而言,有限体积法可视作有限元法和有限差分法的中间产物。
计算流体力学中有限差分法、有限体积法和有限元法的区别
有限元法,有限差分法和有限体积法的区别1. FDM1.1 概念有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
1.2 差分格式(1)从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
(2)从差分的空间形式来考虑,可分为中心格式和逆风格式。
(3)考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
1.3 构造差分的方法构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
2. FEM2.1 概述有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
2.2 原理有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域的解可以看作是由所有单元上的近似解构成。
计算电磁场理论中的有限差分法与有限元法
计算电磁场理论中的有限差分法与有限元法电磁场理论是电磁学的重要组成部分,研究电磁场的分布和变化规律对于解决实际问题具有重要意义。
在计算电磁场中,有限差分法和有限元法是两种常用的数值计算方法。
本文将从理论原理、应用范围和优缺点等方面对这两种方法进行探讨。
有限差分法是一种将连续问题离散化的方法,通过将连续的电磁场分割成网格,然后在每个网格上进行离散计算。
这种方法的基本思想是将微分方程转化为差分方程,然后利用差分方程进行求解。
有限差分法的优点是简单易懂,计算过程直观,适用于各种电磁场问题的求解。
然而,由于差分法中的网格离散化会引入一定的误差,所以在计算精度上存在一定的限制。
与有限差分法相比,有限元法是一种更加精确的数值计算方法。
有限元法将电磁场问题的求解区域划分为有限个小单元,然后在每个小单元上建立适当的插值函数,通过求解代数方程组得到电磁场的近似解。
有限元法的优点是可以处理复杂的几何形状和材料特性,适用于各种边界条件和非线性问题。
然而,有限元法的计算过程相对较为复杂,需要对问题进行合理的离散化和网格划分,同时对于大规模问题,计算量也较大。
在实际应用中,根据具体问题的特点和求解要求,选择合适的数值计算方法是十分重要的。
对于简单的电磁场问题,如一维导线的电流分布,可以选择有限差分法进行求解。
而对于复杂的电磁场问题,如三维空间中的电磁波传播,有限元法更适合。
此外,有限差分法和有限元法还可以结合使用,通过将两种方法的优点相结合,提高计算精度和效率。
除了理论原理和应用范围,有限差分法和有限元法的优缺点也值得关注。
有限差分法的优点是简单易懂,计算过程直观,而且对于一些简单问题可以得到较为准确的结果。
然而,由于差分法中的网格离散化会引入一定的误差,对于复杂问题的求解精度有限。
相比之下,有限元法可以处理复杂的几何形状和材料特性,适用于各种边界条件和非线性问题,计算精度较高。
然而,有限元法的计算过程相对复杂,需要对问题进行合理的离散化和网格划分,同时对于大规模问题计算量较大。
有限元法 有限差分法 有限体积法的区别
三者各有所长:有限差分法:直观,理论成熟,精度可选。
但是不规则区域处理繁琐,虽然网格生成可以使FDM应用于不规则区域,但是对区域的连续性等要求较严。
使用FDM的好处在于易于编程,易于并行。
有限元方法:适合处理复杂区域,精度可选。
缺憾在于内存和计算量巨大。
并行不如FDM和FVM直观。
不过FEM的并行是当前和将来应用的一个不错的方向。
有限容积法:适于流体计算,可以应用于不规则网格,适于并行。
但是精度基本上只能是二阶了。
FVM的优势正逐渐显现出来,FVM在应力应变,高频电磁场方面的特殊的优点正在被人重视。
比较一下:有限容积法和有限差分法:一个区别就是有限容积法的截差是不定的(跟取的相邻点有关,积分方法离散方程),而有限差分就可以直接知道截差(微分方法离散方程)。
有限容积法和有限差分法最本质的区别是,前者是根据积分方程推导出来的(即对每个控制体积分),后者直接根据微分方程推导出来,所以前者的精度不但取决于积分时的精度,还取决与对导数处理的精度,一般有限容积法总体的精度为二阶,因为积分的精度限制,当然有限容积法对于守恒型方程导出的离散方程可以保持守恒型;而后者直接由微分方程导出,不涉及积分过程,各种导数的微分借助Taylor展开,直接写出离散方程,当然不一定有守恒性,精度也和有限容积法不一样,一般有限差分法可以使精度更高一些。
当然二者有联系,有时导出的形式一样,但是概念上是不一样的。
至于有限容积法和有限元相比,有限元在复杂区域的适应性对有限容积是毫无优势可言的,至于有限容积的守恒性,物理概念明显的这些特点,有限元是没有的。
目前有限容积在精度方面与有限元法有些差距。
有限元方法比有限差分优越的方面主要在能适应不规则区域,但是这只是指的是传统意义上的有限差分,现在发展的一些有限差分已经能适应不规则区域。
对于椭圆型方程,如果区域规则,传统有限差分和有限元都能解,在求解效率,这里主要指编程负责度和收敛快慢、内存需要,肯定有限差分有优势。
计算流体力学中有限差分法、有限体积法和有限元法的区别
有限元法,有限差分法和有限体积法的区别1. FDM1.1 概念有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
1.2 差分格式(1)从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
(2)从差分的空间形式来考虑,可分为中心格式和逆风格式。
(3)考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
1.3 构造差分的方法构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
2. FEM2.1 概述有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
2.2 原理有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域的解可以看作是由所有单元上的近似解构成。
有限差分与有限元方法的比较
有限差分与有限元方法的比较有限差分与有限元方法是数值计算中常见的两种方法。
它们都是将复杂的微分方程转化为离散的差分或元素方程,并通过计算机进行求解。
本文就来进行比较分析有限差分和有限元方法的优劣。
一、有限差分方法有限差分方法是在微分方程在特定点求导近似值的简单方法。
方法就是把空间划分成一些简单形状的网格,再在每个小网格中取一个点,代表整个区域内的数值。
利用这些节点上的信息,可以构造微分方程在这些点上的近似值,然后将这些方程组成一个巨大的线性方程组,并用高斯消元法求解。
由于有限差分方法简单快捷,所以它广泛应用于求解流体力学、电磁场理论、热传导等问题。
然而,有限差分方法的缺点也同样明显。
首先,它要求网格非常细密才能得到准确的解,这会导致计算量极大。
其次,有限差分方法只能用于求解正则区域的微分方程,无法处理分界面上出现的奇异性。
最后,它只能求解边界定值问题,不适合处理较为复杂的自适应边值问题。
二、有限元方法有限元方法通过将一个区域分割成许多小被称为有限元的部分,在每个元素内都有一个形函数来描述解,通过将每个元素内的形函数拼合成整个区域内的解。
有限元方法的最大优点在于它具有极高的适应性。
它可以处理非正则区域、自适应边值问题和分界面上出现的奇异性。
此外,它可以使用自适应网格来自适应选择精度以平衡精度和计算成本。
然而,有限元方法在训练和解决问题时的算法比有限差分更复杂,在面对复杂系统时会导致鲁棒性差,精度不够高。
从构建模型以及计算量来看,有限元方法通常比有限差分昂贵。
三、比较与评价综上所述,有限差分和有限元方法都有其优点和缺点,需要根据具体问题选择最合适的方法。
有限差分方法简单快捷,适合边界定值问题,但需要一个非常细密的网格,且只能求解正则区域的偏微分方程。
而有限元方法适用范围广,具有很强的适应性,但训练和解决问题的算法比有限差分更复杂,对处理复杂系统时,鲁棒性不足,精度不够高。
因此,在实际应用中,应根据具体问题的特点选择最合适的方法来获得更准确的结果。
数值计算中的有限元和有限差分方法
数值计算中的有限元和有限差分方法数值计算是一种利用数字来求解数学问题的技术。
在各个领域中,数值计算都被广泛应用,尤其是在工程计算中具有重要的地位。
有限元和有限差分方法是数值计算的两个重要工具,本文将介绍它们的原理、优缺点以及应用。
一、有限元方法有限元方法(Finite Element Method,简称FEM)是一种适用于工程力学、流体力学、热传导等问题的数值计算方法。
首先将问题区域离散化成若干个小区域,每个小区域称为有限元;然后通过对每个有限元的变形、应力和应变的计算,得到整个问题的解。
有限元方法的基本原理是建立一个局部变形和应力的数学模型,借助于位移和应力的离散函数来代表局部信息,并将不连续的位移和应力函数在结点处相互连接,形成一个连续作用的整体模型,从而求解整个问题的解。
通过该方法可以精确地求解各种材料构件的形变、应变以及应力分布等问题,并且具有灵活性和广泛性。
有限元方法的优点是求解精度较高,分析结果可靠。
可以分析复杂的问题以及非线性问题,并可进行多物理场耦合分析。
此外,还可以基于现有的有限元软件进行建模分析,避免重复造轮子。
然而,它也存在限制,例如建模时需要对问题进行适当的假设,并且需要对材料力学性质等信息有一定的了解。
此外,考虑更复杂的物理现象时,需要使用更高阶的元来表示求解方程,这会导致计算量增加,计算时间增长。
二、有限差分法有限差分方法(Finite Difference Method,简称FDM)是一种常用的求解微分方程的数值计算方法。
该方法将微分方程中的导数用有限差分的形式表示出来,从而将连续问题离散化成为一个离散点问题,并通过计算在各个离散点上函数值的差分,从而得到微分方程的数值解。
有限差分方法的基本思想是将连续函数转化为离散函数,然后在离散点上近似求解微分方程。
该方法简单易懂,计算量小,代码实现相对容易。
因此,将微分方程离散化是数值计算中经常采用的方法。
与有限元方法相比,有限差分方法在处理一些简单问题的时候表现更好,计算速度快,精度也有保障。
有限元法有限差分法有限体积法
有限元法有限差分法有限体积法
有限元法 (Finite Element Method)、有限差分法 (Finite Difference Method)、有限体积法 (Finite Volume Method) 都是
常见的数值方法,用于求解各种科学和工程问题的偏微分方程。
有限元法是一种离散化技术,将待求解的问题分解成多个简化部分,并分别通过逼近方法近似求解。
这种技术能提供问题在离散域的数值解,并且能够适应复杂的几何形状和物理特性。
有限差分法是一种通过近似求解微分方程的方法,通过将求解区域离散化成网格,并在网格交点处近似微分方程来进行计算。
这种方法对解析完全或者解析复杂的问题很有效,但是对于复杂的几何形状和物理特性有一定的限制。
有限体积法是一种通过求解离散的控制体积内的物理量平衡方程来求解偏微分方程的方法。
这种方法适用于处理包含物质交换或流量的宏观物理问题,并且能够直接处理不规则网格。
总的来说,这三种数值方法各有适用范围,需要根据实际问题的特点进行选择。
有限元法与有限差分法的主要区别
有限元法与有限差分法的主要区别有限元法与有限差分法的主要区别有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
从差分的空间形式来考虑,可分为中心格式和逆风格式。
考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
有限差分,有限元,有限体积等离散方法的区别介绍
有限差分,有限元,有限体积等离散方法的区别介绍1 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
从差分的空间形式来考虑,可分为中心格式和逆风格式。
考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
2 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
有限差分-有限元-有限体积等的区别介绍
有限差分,有限元,有限体积等离散方法的区别介绍1 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
从差分的空间形式来考虑,可分为中心格式和逆风格式。
考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
2 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
有限元法,有限差分法和有限体积法的区别
有限元法,有限差分法和有限体积法的区别有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
从差分的空间形式来考虑,可分为中心格式和逆风格式。
考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
有限差分,有限元,有限体积等离散方法的区别介绍
有限差分,有限元,有限体积等等离散方法的区别介绍一、区域离散化所谓区域离散化,实质上就是用一组有限个离散的点来代替原来连续的空间。
实施过程是;把所计算的区域划分成许多互不重迭的子区域,确定每个子区域的节点位置及该节点所代表的控制容积。
节点:需要求解的未知物理量的几何位置;控制容积:应用控制方程或守恒定律的最小几何单位。
一般把节点看成是控制容积的代表。
控制容积和子区域并不总是重合的。
在区域离散化过程开始时,由一系列与坐标轴相应的直线或曲线簇所划分出来的小区域称为子区域。
网格是离散的基础,网格节点是离散化物理量的存储位置。
大家都知道,常用的离散化方法有:有限差分法,有限元法,有限体积法。
1. 有限差分法是数值解法中最经典的方法。
它是将求解区域划分为差分网格,用有限个网格节点代替连续的求解域,然后将偏微分方程(控制方程)的导数用差商代替,推导出含有离散点上有限个未知数的差分方程组。
这种方法发展比较早,比较成熟,较多用于求解双曲线和抛物线型问题。
用它求解边界条件复杂、尤其是椭圆型问题不如有限元法或有限体积法方便。
2. 有限元法是将一个连续的求解域任意分成适当形状的许多微小单元,并于各小单元分片构造插值函数,然后根据极值原理(变分或加权余量法),将问题的控制方程转化为所有单元上的有限元方程,把总体的极值作为各单元极值之和,即将局部单元总体合成,形成嵌入了指定边界条件的代数方程组,求解该方程组就得到各节点上待求的函数值。
对椭圆型问题有更好的适应性。
有限元法求解的速度较有限差分法和有限体积法慢,在商用CFD软件中应用并不广泛。
目前的商用CFD软件中,FIDAP采用的是有限元法。
3. 有限体积法又称为控制体积法,是将计算区域划分为网格,并使每个网格点周围有一个互不重复的控制体积,将待解的微分方程对每个控制体积积分,从而得到一组离散方程。
其中的未知数十网格节点上的因变量。
子域法加离散,就是有限体积法的基本方法。
就离散方法而言,有限体积法可视作有限元法和有限差分法的中间产物。
有限差分,有限元,有限体积等的区别介绍
有限差分,有限元,有限体积等离散方法的区别介绍1 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
从差分的空间形式来考虑,可分为中心格式和逆风格式。
考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
2 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
有限元法与有限差分法的主要区别
有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
从差分的空间形式来考虑,可分为中心格式和逆风格式。
考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。
根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。
从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。
不同的组合同样构成不同的有限元计算格式。
对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的平方误差最小;在配置法中,先在计算域内选取N个配置点。
令近似解在选定的N个配置点上严格满足微分方程,即在配置点上令方程余量为0。
插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。
有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身,还要求它的导数值在插值点取已知值,称为哈密特(Hermite)多项式插值。
单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等。
常采用的无因次坐标是一种局部坐标系,它的定义取决于单元的几何形状,一维看作长度比,二维看作面积比,三维看作体积比。
在二维有限元中,三角形单元应用的最早,近来四边形等参元的应用也越来越广。
对于二维三角形和四边形电源单元,常采用的插值函数为有La grange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。
对于有限元方法,其基本思路和解题步骤可归纳为(1)建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的出发点。
(2)区域单元剖分,根据求解区域的形状及实际问题的物理特点,将区域剖分为若干相互连接、不重叠的单元。
区域单元划分是采用有限元方法的前期准备工作,这部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标,同时还需要列出自然边界和本质边界的节点序号和相应的边界值。
(3)确定单元基函数,根据单元中节点数目及对近似解精度的要求,选择满足一定插
值条件的插值函数作为单元基函数。
有限元方法中的基函数是在单元中选取的,由于各单元具有规则的几何形状,在选取基函数时可遵循一定的法则。
(4)单元分析:将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近;再将近似函数代入积分方程,并对单元区域进行积分,可获得含有待定系数(即单元中各节点的参数值)的代数方程组,称为单元有限元方程。
(5)总体合成:在得出单元有限元方程之后,将区域中所有单元有限元方程按一定法则进行累加,形成总体有限元方程。
(6)边界条件的处理:一般边界条件有三种形式,分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。
对于自然边界条件,一般在积分表达式中可自动得到满足。
对于本质边界条件和混合边界条件,需按一定法则对总体有限元方程进行修正满足。
(7)解有限元方程:根据边界条件修正的总体有限元方程组,是含所有待定未知量的封闭方程组,采用适当的数值计算方法求解,可求得各节点的函数值。
有限体积法(Finite Volume Method)又称为控制体积法。
其基本思路是:将计算区域划分为一系列不重复的控制体积,并使每个网格点周围有一个控制体积;将待解的微分方程对每一个控制体积积分,便得出一组离散方程。
其中的未知数是网格点上的因变量的数值。
为了求出控制体积的积分,必须假定值在网格点之间的变化规律,即假设值的分段的分布的分布剖面。
从积分区域的选取方法看来,有限体积法属于加权剩余法中的子区域法;从未知解的近似方法看来,有限体积法属于采用局部近似的离散方法。
简言之,子区域法属于有限体积发的基本方法。
有限体积法的基本思路易于理解,并能得出直接的物理解释。
离散方程的物理意义,就是因变量在有限大小的控制体积中的守恒原理,如同微分方程表示因变量在无限小的控制体积中的守恒原理一样。
限体积法得出的离散方程,要求因变量的积分守恒对任意一组控制体积都得到满足,对整个计算区域,自然也得到满足。
这是有限体积法吸引人的优点。
有一些离散方法,例如有限差分法,仅当网格极其细密时,离散方程才满足积分守恒;而有限体积法即使在粗网格情况下,也显示出准确的积分守恒。
就离散方法而言,有限体积法可视作有限单元法和有限差分法的中间物。
有限单元法必须假定值在网格点之间的变化规律(既插值函数),并将其作为近似解。
有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。
有限体积法只寻求的结点值,这与有限差分法相类似;但有限体积法在寻求控制体积的积分时,必须假定值在网格点之间的分布,这又与有限单元法相类似。
在有限体积法中,插值函数只用于计算控制体积的积分,得出离散方程之后,便可忘掉插值函数;如果需要的话,可以对微分方程中不同的项采取不同的插值函数。
如有侵权请联系告知删除,感谢你们的配合!。