第一章_常用数值分析方法§4_有限差分法与有限单元法
- 格式:ppt
- 大小:720.00 KB
- 文档页数:15
有限差分方法(Finite Differential Method)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以泰勒级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
从差分的空间形式来考虑,可分为中心格式和逆风格式。
考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
有限元法(Finite Element Method)的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
有限差分法原理有限差分法(Finite Difference Method)是一种常见的数值分析方法,广泛应用于工程、物理、经济等领域的数值模拟和计算中。
它的基本原理是将微分方程转化为差分方程,通过在空间和时间上进行离散,将连续的问题转化为离散的问题,从而用计算机进行求解。
有限差分法在实际工程中具有重要的应用价值,本文将对有限差分法的原理进行详细介绍。
有限差分法的基本思想是将求解的区域进行网格划分,然后利用差分近似代替微分运算,通过有限差分近似的方式将微分方程转化为代数方程组,进而求解出数值解。
有限差分法的核心在于如何进行差分近似,以及如何选择合适的差分格式。
一般来说,差分格式可以分为前向差分、后向差分、中心差分等不同类型,根据不同问题的特点和求解精度的要求,选择合适的差分格式对问题进行离散化处理。
在空间上进行离散化时,通常采用均匀网格划分的方法,将求解区域划分为若干个小区间,每个小区间内的差分近似都可以通过相似的方式进行处理。
而在时间上进行离散化时,则需要根据具体问题选择合适的时间步长,通过逐步迭代的方式求解出时间上的数值解。
有限差分法的原理可以用一个简单的一维热传导方程来进行说明。
假设有一根长度为L的杆,其温度分布满足一维热传导方程,即∂u/∂t = α∂^2u/∂x^2,其中u(x,t)表示杆上某一点的温度分布,α为热传导系数。
我们可以将空间上的区域进行均匀网格划分,时间上进行等间隔的离散化,然后利用差分近似代替微分运算,最终得到一个关于时间和空间上温度分布的差分方程组,通过迭代计算得到数值解。
有限差分法作为一种数值计算方法,其精度和稳定性受到网格划分和时间步长的影响。
通常来说,网格划分越精细,时间步长越小,数值解的精度越高,但计算量也会相应增加。
因此,在实际应用中需要根据具体问题的要求和计算资源的限制进行合理的选择。
总之,有限差分法是一种重要的数值计算方法,通过将微分方程转化为差分方程,利用计算机进行求解,可以有效地解决实际工程中的复杂问题。
有限差分法和有限元法
有限差分法(Finite Difference Method)和有限元法(Finite Element Method)是两种常用的数值计算方法,用于求解偏微分方程的数值解。
有限差分法是通过将求解区域离散化为网格,然后在各个网格节点处用差分逼近偏微分方程中的导数项,将偏微分方程转化为代数方程组。
通过求解这个方程组,可以得到离散节点上的数值解。
有限差分法适用于一维、二维或三维的问题,可用来处理线性或非线性、稳定或非稳定的偏微分方程。
有限差分法的优点是简单易实现,容易理解和计算,但是对于复杂的几何形状和边界条件,离散网格的选择可能会对精度和计算结果产生较大的影响。
有限元法则是通过将求解区域划分为互不重叠的有限元,每个有限元内部采用局部函数近似原方程,然后将所有有限元的近似解拼接在一起,形成整个求解区域上的近似解。
有限元法通常在每个有限元上构造基函数,通过求解代数方程组确定基函数的系数,从而得到整个求解区域上的数值解。
有限元法适用于一维、二维或三维的问题,能够处理各种几何形状和边界条件,适用范围更广。
有限元法的优点是对复杂几何形状的适应性好,精度高,但是相对于有限差分法而言,复杂度较高,需要更多的计算量和计算时间。
总体来说,有限差分法更适用于简单的几何形状和边界条件,而有限元法更适用于复杂的几何形状和边界条件。
两种方法在
实际的工程和科学计算中都有广泛的应用,选择哪种方法取决于具体问题的性质和求解的要求。
常用数值分析方法1.插值方法插值是通过已知数据点的近似值,获得未知位置上的函数值。
常用的插值方法包括拉格朗日插值、牛顿插值和分段线性插值等。
插值方法通常用于数据的光滑处理、曲线拟合和函数逼近等问题。
2.数值微分与积分方法数值微分是通过有限差分等方法,对实际问题的函数进行求导。
数值积分则是通过数值方法求解复杂函数的积分。
常用的数值微分与积分方法包括欧拉法、龙格-库塔法和辛算法等。
3.非线性方程求解非线性方程求解是求解形如f(x)=0的方程,其中f(x)是一个非线性函数。
常用的非线性方程求解方法包括二分法、牛顿法和割线法等。
这些方法基于不同的数学原理来逼近方程的根。
4.线性方程组求解线性方程组求解是求解形如Ax=b的方程组,其中A是一个矩阵,b 是一个向量。
常用的线性方程组求解方法包括高斯消元法、LU分解和迭代法等。
这些方法可以高效地求解大规模的线性方程组。
5.最小二乘法最小二乘法是一种用于拟合实验或观测数据的方法。
它通过最小化观测数据与理论模型之间的残差平方和,得到最佳的参数估计。
最小二乘法广泛应用于曲线拟合、回归分析和信号处理等领域。
6.数值优化数值优化是在约束条件下求解最优化问题的方法。
常用的数值优化方法包括梯度下降法、共轭梯度法和拟牛顿法等。
这些方法可以在函数复杂或维度高的情况下,有效地寻找最优解。
7.偏微分方程数值解法偏微分方程数值解法是用数值方法解决偏微分方程的方法。
常用的数值解法包括有限差分法、有限元法和谱方法等。
这些方法广泛应用于物理学、工程学和金融学等领域,可以模拟和预测复杂现象。
总之,数值分析方法在科学和工程领域中起着重要的作用。
通过数学和计算机的结合,数值分析使得复杂计算变得简单,从而有效解决各种实际问题。
有限差分,有限元,有限体积等等离散方法的区别介绍一、区域离散化所谓区域离散化,实质上就是用一组有限个离散的点来代替原来连续的空间。
实施过程是;把所计算的区域划分成许多互不重迭的子区域,确定每个子区域的节点位置及该节点所代表的控制容积。
节点:需要求解的未知物理量的几何位置;控制容积:应用控制方程或守恒定律的最小几何单位。
一般把节点看成是控制容积的代表。
控制容积和子区域并不总是重合的。
在区域离散化过程开始时,由一系列与坐标轴相应的直线或曲线簇所划分出来的小区域称为子区域。
网格是离散的基础,网格节点是离散化物理量的存储位置。
大家都知道,常用的离散化方法有:有限差分法,有限元法,有限体积法。
1. 有限差分法是数值解法中最经典的方法。
它是将求解区域划分为差分网格,用有限个网格节点代替连续的求解域,然后将偏微分方程(控制方程)的导数用差商代替,推导出含有离散点上有限个未知数的差分方程组。
这种方法发展比较早,比较成熟,较多用于求解双曲线和抛物线型问题。
用它求解边界条件复杂、尤其是椭圆型问题不如有限元法或有限体积法方便。
2. 有限元法是将一个连续的求解域任意分成适当形状的许多微小单元,并于各小单元分片构造插值函数,然后根据极值原理(变分或加权余量法),将问题的控制方程转化为所有单元上的有限元方程,把总体的极值作为各单元极值之和,即将局部单元总体合成,形成嵌入了指定边界条件的代数方程组,求解该方程组就得到各节点上待求的函数值。
对椭圆型问题有更好的适应性。
有限元法求解的速度较有限差分法和有限体积法慢,在商用CFD软件中应用并不广泛。
目前的商用CFD软件中,FIDAP采用的是有限元法。
3. 有限体积法又称为控制体积法,是将计算区域划分为网格,并使每个网格点周围有一个互不重复的控制体积,将待解的微分方程对每个控制体积积分,从而得到一组离散方程。
其中的未知数十网格节点上的因变量。
子域法加离散,就是有限体积法的基本方法。
就离散方法而言,有限体积法可视作有限元法和有限差分法的中间产物。
有限差分,有限元,有限体积等离散方法的区别介绍1 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
从差分的空间形式来考虑,可分为中心格式和逆风格式。
考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
2 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。
三者各有所长:有限差分法:直观,理论成熟,精度可选。
但是不规则区域处理繁琐,虽然网格生成可以使FDM应用于不规则区域,但是对区域的连续性等要求较严。
使用FDM的好处在于易于编程,易于并行。
有限元方法:适合处理复杂区域,精度可选。
缺憾在于内存和计算量巨大。
并行不如FDM和FVM直观。
不过FEM的并行是当前和将来应用的一个不错的方向。
有限容积法:适于流体计算,可以应用于不规则网格,适于并行。
但是精度基本上只能是二阶了。
FVM的优势正逐渐显现出来,FVM在应力应变,高频电磁场方面的特殊的优点正在被人重视。
比较一下:有限容积法和有限差分法:一个区别就是有限容积法的截差是不定的(跟取的相邻点有关,积分方法离散方程),而有限差分就可以直接知道截差(微分方法离散方程)。
有限容积法和有限差分法最本质的区别是,前者是根据积分方程推导出来的(即对每个控制体积分),后者直接根据微分方程推导出来,所以前者的精度不但取决于积分时的精度,还取决与对导数处理的精度,一般有限容积法总体的精度为二阶,因为积分的精度限制,当然有限容积法对于守恒型方程导出的离散方程可以保持守恒型;而后者直接由微分方程导出,不涉及积分过程,各种导数的微分借助Taylor展开,直接写出离散方程,当然不一定有守恒性,精度也和有限容积法不一样,一般有限差分法可以使精度更高一些。
当然二者有联系,有时导出的形式一样,但是概念上是不一样的。
至于有限容积法和有限元相比,有限元在复杂区域的适应性对有限容积是毫无优势可言的,至于有限容积的守恒性,物理概念明显的这些特点,有限元是没有的。
目前有限容积在精度方面与有限元法有些差距。
有限元方法比有限差分优越的方面主要在能适应不规则区域,但是这只是指的是传统意义上的有限差分,现在发展的一些有限差分已经能适应不规则区域。
对于椭圆型方程,如果区域规则,传统有限差分和有限元都能解,在求解效率,这里主要指编程负责度和收敛快慢、内存需要,肯定有限差分有优势。
传热学_燕山大学中国大学mooc课后章节答案期末考试题库2023年1.求解导热问题时,常用的数值计算法有:()参考答案:有限单元法_有限差分法2.导热问题数值求解的基本步骤包括:()参考答案:代数方程组的求解_问题的数学描述_区域离散_节点离散方程的建立3.求解传热问题的基本方法有哪些?()参考答案:理论分析_实验测量_数值计算4.解析法的特点:()参考答案:能够得到精确解_解有普遍性5.采用有限差分法计算导热问题,可以得到节点处温度值的精确解。
参考答案:错误6.有限差分法较有限元法更适合求解不规则物体边界的导热问题。
参考答案:错误7.高斯-赛德尔迭代法比简单迭代法计算精度高。
参考答案:错误8.下面()种说法是不正确的。
参考答案:同类物理现象就是彼此相似的现象9.流体横掠单管时,相似准则中的定型尺寸取为()。
参考答案:管外径10.流体在圆管内流动时,相似准则中的定型尺寸取为()。
参考答案:管内径11.对流换热系数不是流体的物性参数,是一个反应对流换热复杂过程的多变量的函数。
参考答案:正确12.自然对流换热与强制对流换热都是对流换热,属于是同类物理现象。
参考答案:错误13.不符合肋片描述的是()参考答案:减小换热量作用14.对温度梯度描述错误的是()参考答案:沿温度增加方向为负15.黑体是指能够吸收投射到其表面上的全部辐射能的物体。
参考答案:正确16.导热是固体内部或固体之间的传热过程。
参考答案:错误17.工业上,多数工程材料作为灰体处理不会引起较大的误差。
参考答案:正确18.导热问题数值求解的基本思想是:()参考答案:用求解区域上离散点的温度集合代替连续的温度场,计算难度小_用代数方程代替微分方程19.一般而言,与非金属材料相比,金属材料的导热系数()参考答案:较大20.五种具有实际意义的换热过程为:导热、对流换热、复合换热、传热过程和()参考答案:辐射换热21.热量传递的三种基本方式是()参考答案:热对流、导热、热辐射22.由表面 1 和表面 2 组成的封闭系统中,【图片】()【图片】参考答案:可能大于、等于或小于23.流体纵掠平板时,对流动边界层的描述,下列第()种说法是不准确的。
有限元法与有限差分法的主要区别有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式与高阶格式。
从差分的空间形式来考虑,可分为中心格式与逆风格式。
考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况与柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分与二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。
通过对时间与空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
有限元方法的基础是变分原理与加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数与插值函数形式,便构成不同的有限元方法。
有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。