有限元法有限差分法有限体积法的区别

有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。

该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。

有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。

对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。

从差分的空间形式来考虑,可分为中心格式和逆风格式。

考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。

目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。

差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。

其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。

通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。

有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。

采用不同的权函数和插值函数形式,便构成不同的有限元方法。

有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。

在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。

在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。

根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。

从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。

不同的组合同样构成不同的有限元计算格式。

对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的平方误差最小;在配置法中,先在计算域内选取N个配置点。

令近似解在选定的N个配置点上严格满足微分方程,即在配置点上令方程余量为0。

插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。

有限元插值函数分为两大类,一类只要求插值多项式本身在插值
点取已知值,称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身,还要求它的导数值在插值点取已知值,称为哈密特(Hermite)多项式插值。

单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等。

常采用的无因次坐标是一种局部坐标系,它的定义取决于单元的几何形状,一维看作长度比,二维看作面积比,三维看作体积比。

在二维有限元中,三角形单元应用的最早,近来四边形等参元的应用也越来越广。

对于二维三角形和四边形电源单元,常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。

对于有限元方法,其基本思路和解题步骤可归纳为
(1)建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的出发点。

(2)区域单元剖分,根据求解区域的形状及实际问题的物理特点,将区域剖分为若干相互连接、不重叠的单元。

区域单元划分是采用有限元方法的前期准备工作,这部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标,同时还需要列出自然边界和本质边界的节点序号和相应的边界值。

(3)确定单元基函数,根据单元中节点数目及对近似解精度的要求,选择满足一定插值条件的插值函数作为单元基函数。

有限元方法中的基函数是在单元中选取的,由于各单元具有规则的几何形状,在选取基函数时可遵循一定的法则。

(4)单元分析:将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近;再将近似函数代入积分方程,并对单元区域进行积分,可获得含有待定系数(即单元中各节点的参数值)的代数方程组,称为单元有限元方程。

(5)总体合成:在得出单元有限元方程之后,将区域中所有单元有限元方程按一定法则进行累加,形成总体有限元方程。

(6)边界条件的处理:一般边界条件有三种形式,分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。

对于自然边界条件,一般在积分表达式中可自动得到满足。

对于本质边界条件和混合边界条件,需按一定法则对总体有限元方程进行修正满足。

(7)解有限元方程:根据边界条件修正的总体有限元方程组,是含所有待定未知量的封闭方程组,采用适当的数值计算方法求解,可求得各节点的函数值。

有限体积法(Finite Volume Method)又称为控制体积法。

其基本思路是:将计算区域划分为一系列不重复的控制体积,并使每个网格点周围有一个控制体积;将待解的微分方程对每一个控制体积积分,便得
出一组离散方程。

其中的未知数是网格点上的因变量的数值。

为了求出控制体积的积分,必须假定值在网格点之间的变化规律,即假设值的分段的分布的分布剖面。

从积分区域的选取方法看来,有限体积法属于加权剩余法中的子区域法;从未知解的近似方法看来,有限体积法属于采用局部近似的离散方法。

简言之,子区域法属于有限体积发的基本方法。

有限体积法的基本思路易于理解,并能得出直接的物理解释。

离散方程的物理意义,就是因变量在有限大小的控制体积中的守恒原理,如同微分方程表示因变量在无限小的控制体积中的守恒原理一样。

限体积法得出的离散方程,要求因变量的积分守恒对任意一组控制体积都得到满足,对整个计算区域,自然也得到满足。

这是有限体积法吸引人的优点。

有一些离散方法,例如有限差分法,仅当网格极其细密时,离散方程才满足积分守恒;而有限体积法即使在粗网格情况下,也显示出准确的积分守恒。

就离散方法而言,有限体积法可视作有限单元法和有限差分法的中间物。

有限单元法必须假定值在网格点之间的变化规律(既插值函数),并将其作为近似解。

有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。

有限体积法只寻求的结点值,这与有限差分法相类似;但有限体积法在寻求控制体积的积分时,必须假定值在网格点之间的分布,这又与有限单元法相类似。

在有限体积法中,插值函数只用于计算控制体积的积分,得出离散方程之后,便可忘掉插值函数;如果需要的话,可以对微分方程中不同的项采取不同的插值函数。

合集下载

偏微分方程数值方法

偏微分方程数值方法

偏微分方程数值方法偏微分方程(Partial Differential Equation, PDE)是数学中的一种重要的方程类型,它描述了一个函数的多个变量的变化关系。

解决偏微分方程的数值方法在科学和工程领域有着广泛的应用。

本文将介绍几种常见的偏微分方程数值方法,并对其进行详细阐述。

1. 差分法(Finite Difference Method):差分法是最早也是最直接的一种数值方法,它基于连续函数在一些点的导数可以用它的前向、后向或中心的差商来近似的思想。

偏微分方程的差分格式包括向前差分法、向后差分法和中心差分法等。

对于二维的偏微分方程,可以采用网格化的方式将空间离散化,然后利用差分法进行近似求解。

2. 有限元法(Finite Element Method):有限元法是一种基于原始形式或变分形式对偏微分方程进行离散化的方法。

在有限元法中,将求解域分割成许多小的、简单的几何单元,然后在每个单元上构建近似解函数和试验函数。

通过构建弱形式并应用基本的变分原理,可以得到离散化的方程组,并通过求解这个方程组来得到数值解。

3. 有限差分法(Finite Difference Method):有限差分法是一种将连续的偏微分方程离散化成差分方程的方法。

它与差分法的主要区别在于有限差分法不需要对求解域进行网格化,而是直接在连续的求解域上进行离散化。

将偏微分方程中的导数通过差商来近似,然后通过求解离散化的差分方程来得到数值解。

4. 有限体积法(Finite Volume Method):有限体积法是一种将偏微分方程离散化为离散体积元的方法。

在有限体积法中,将求解域划分成离散的控制体积,然后通过对控制体积的积分运算,将偏微分方程转化为离散的代数方程组。

然后通过求解得到的代数方程组,可以得到数值解。

以上介绍的只是几种常见的偏微分方程数值方法,实际上还有很多其他的方法,如边界元法(Boundary Element Method)、谱方法(Spectral Method)、逆问题方法(Inverse Problem Method)等。

有限元素法有限体积法有限差分法有限容积法的区别

有限元素法有限体积法有限差分法有限容积法的区别

1.1 概念有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。

该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。

有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。

1.2 差分格式(1)从格式的精度来划分,有一阶格式、二阶格式和高阶格式。

(2)从差分的空间形式来考虑,可分为中心格式和逆风格式。

(3)考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。

目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。

差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

1.3 构造差分的方法构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。

其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。

通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。

2. FEM2.1 概述有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。

采用不同的权函数和插值函数形式,便构成不同的有限元方法。

2.2 原理有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。

在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。

有限元法,有限差分法和有限体积法的区别

有限元法,有限差分法和有限体积法的区别

有限元法,有限差分法和有限体积法的区别有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。

该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。

有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。

对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。

从差分的空间形式来考虑,可分为中心格式和逆风格式。

考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。

目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。

差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。

其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。

通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。

有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。

采用不同的权函数和插值函数形式,便构成不同的有限元方法。

有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。

在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。

有限元计算固体力学和有限体积法计算流体力学的本质区别

有限元计算固体力学和有限体积法计算流体力学的本质区别

有限元计算固体力学和有限体积法计算流体力学的本质区别有限元计算固体力学和有限体积法计算流体力学的本质区别2010-12-23有限元计算固体力学和有限体积法计算流体力学的本质区别是什么?如题。

回答汇总:可以说一个是基于变分原理,一个是基于守恒原理。

有限体积法其实就是有限元中讲的加权残数法中的子域法。

具体请参考《传热与流体流动的数值计算》第33页。

固体采用拉格朗日方法,流体采用欧拉方法。

1.FDM 1.1概念有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。

该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。

有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。

1.2差分格式(1)从格式的精度来划分,有一阶格式、二阶格式和高阶格式。

(2)从差分的空间形式来考虑,可分为中心格式和逆风格式。

(3)考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。

目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。

差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

1.3构造差分的方法构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。

其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。

通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。

2..FEM 2.1概述有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。

有限差分,有限元,有限体积等离散方法的区别介绍

有限差分,有限元,有限体积等离散方法的区别介绍

有限差分,有限元,有限体积等等离散方法的区别介绍一、区域离散化所谓区域离散化,实质上就是用一组有限个离散的点来代替原来连续的空间。

实施过程是;把所计算的区域划分成许多互不重迭的子区域,确定每个子区域的节点位置及该节点所代表的控制容积。

节点:需要求解的未知物理量的几何位置;控制容积:应用控制方程或守恒定律的最小几何单位。

一般把节点看成是控制容积的代表。

控制容积和子区域并不总是重合的。

在区域离散化过程开始时,由一系列与坐标轴相应的直线或曲线簇所划分出来的小区域称为子区域。

网格是离散的基础,网格节点是离散化物理量的存储位置。

大家都知道,常用的离散化方法有:有限差分法,有限元法,有限体积法。

1. 有限差分法是数值解法中最经典的方法。

它是将求解区域划分为差分网格,用有限个网格节点代替连续的求解域,然后将偏微分方程(控制方程)的导数用差商代替,推导出含有离散点上有限个未知数的差分方程组。

这种方法发展比较早,比较成熟,较多用于求解双曲线和抛物线型问题。

用它求解边界条件复杂、尤其是椭圆型问题不如有限元法或有限体积法方便。

2. 有限元法是将一个连续的求解域任意分成适当形状的许多微小单元,并于各小单元分片构造插值函数,然后根据极值原理(变分或加权余量法),将问题的控制方程转化为所有单元上的有限元方程,把总体的极值作为各单元极值之和,即将局部单元总体合成,形成嵌入了指定边界条件的代数方程组,求解该方程组就得到各节点上待求的函数值。

对椭圆型问题有更好的适应性。

有限元法求解的速度较有限差分法和有限体积法慢,在商用CFD软件中应用并不广泛。

目前的商用CFD软件中,FIDAP采用的是有限元法。

3. 有限体积法又称为控制体积法,是将计算区域划分为网格,并使每个网格点周围有一个互不重复的控制体积,将待解的微分方程对每个控制体积积分,从而得到一组离散方程。

其中的未知数十网格节点上的因变量。

子域法加离散,就是有限体积法的基本方法。

就离散方法而言,有限体积法可视作有限元法和有限差分法的中间产物。

有限差分、有限元区别

有限差分、有限元区别

有限差分方法(Finite Differential Method)是计算机数值模拟最早采用的方法,至今仍被广泛运用。

该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。

有限差分法以泰勒级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。

对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。

从差分的空间形式来考虑,可分为中心格式和逆风格式。

考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。

目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。

差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。

其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。

通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。

有限元法(Finite Element Method)的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。

采用不同的权函数和插值函数形式,便构成不同的有限元方法。

有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。

在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。

有限体积法 有限差分法 有限元法

有限体积法有限差分法有限元法有限体积法、有限差分法、有限元法是三种数学方法,它们分别用于求解偏微分方程问题。

在工程、物理、气象、地质和生物等领域中都有广泛的应用。

它们之间的区别在于采用不同的逼近方法和离散化技术。

有限体积法是一种数值方法,通过离散化空间来对流体动力学等宏观定律进行描述。

通过建立小区域的质量平衡方程,计算该区域内的物理量积分,并通过解析物理方程,确定小区域物理量的变化率。

这种方法适用于偏微分方程的求解,同时可以避免非物理现象的出现,在计算过程中也不会涉及到边界值问题。

有限差分法是一种离散化的数学方法,可以将一个连续的函数微分方程转换成一个差分方程。

在计算差分方程时,需要将函数在有限点处进行展开,将其转化为有限项的多项式。

这个多项式可以用于近似函数,从而求解微分方程的数值解。

有限差分法可以应用于所有类型的偏微分方程,包括椭圆型、双曲型和抛物型方程。

有限元法是一种基于函数空间分析的数学方法,用于解决连续性和光滑性强的问题。

将连续问题转化为一组代数方程,通过将求解域分成无限多的小元素或区域,将标量或矢量场用有限个基函数来逼近。

将这些基函数带入微分方程中,并将未知系数替换为求解域中的节点上的未知量,就可以得到代数方程组。

最终,通过解决代数方程组来计算微分方程的数值解。

总之,有限体积法、有限差分法和有限元法是三种常用的数值方法。

它们在求解各种复杂偏微分方程方面都具有优越性。

但是它们在适用条件、误差分析、计算量等方面都有各自独特的特点和限制,因此需要根据不同的实际应用来选择和使用。

有限元法与有限差分法的主要区别

有限元法与有限差分法的主要区别有限元法与有限差分法的主要区别有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。

该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。

有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。

对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。

从差分的空间形式来考虑,可分为中心格式和逆风格式。

考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。

目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。

差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。

其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。

通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。

有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。

采用不同的权函数和插值函数形式,便构成不同的有限元方法。

有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。

在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。

时域有限fdtd

【转帖】有限差分,有限元,有限体积等离散方法的区别介绍作者: cool_smile (站内联系TA) 发布: 2009-10-23以下介绍是本人从网络上搜集的,供计算数学虫子参考。

也许小木虫论坛有,我没搜索到。

欢迎大家补充内容。

1 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。

该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。

有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。

对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。

从差分的空间形式来考虑,可分为中心格式和逆风格式。

考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。

目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。

差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。

其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。

通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。

2 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。

采用不同的权函数和插值函数形式,便构成不同的有限元方法。

在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。

有限差分 有限元 有限体积

有限差分有限元有限体积有限差分、有限元和有限体积是数值计算方法中常用的三种离散化方法。

它们的核心思想是将微分方程式转化为一系列有限的点上的代数方程式,将连续问题转化为离散问题。

一、有限差分法有限差分法是将微分方程的导数用差商来逼近的方法,用差商来代替微分运算。

用区间的两个端点上的函数值之差来代替区间内函数导数的平均值。

在连续的区间上进行近似,大大减小了计算量。

有限差分法是一种较为简单的数值解法,适用于规则网格的微分方程求解,被广泛应用在流体力学、结构力学、电场问题等领域中。

二、有限元法有限元法是将求解域分成若干个划分元,然后在每个单元内用多项式函数逼近问题的解,最终利用点、线、面元件的连接关系来求解整体问题的一种方法。

该方法可以处理复杂的几何形状和物理变化,适用于非常规的边界条件和材料特性,解决超过几百万自由度的三维大规模问题。

三、有限体积法有限体积法是将求解域分成若干个控制体,对质量、能量、动量等守恒量在各个控制体上进行积分,从而推导出控制体内分布的方程。

该方法以区域的体积分为基础,在各个控制体内求解守恒方程。

该方法适用于复杂的多组分、多相流动的领域以及非稳态或非线性问题。

无论是有限差分、有限元还是有限体积法,其核心思想都是通过把连续的微分方程式离散求解,从而转化为一系列有限的点上的代数方程式,解决了连续问题转化为离散问题的过程,从而通过离散求解代数方程式来得到问题的解。

这三种数值计算方法的应用使科学计算得以更加高效、精确地进行,对现代计算、科学技术的推进起到了巨大的贡献。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
相关文档
最新文档