传热学数值计算方法课件.pptVIP

Thermal第三章离散化方法本章任务：前两章已阐述了获得物理现象的理论预测有一些明显的好处，且一些有趣的现象是受一些微分方程支配的，这些方程已用通用方程归纳了，接下来的任务就是推导求解这个通用方程的方法—即所谓的离散化方法。推导前的假设：设通用变量?仅仅是一个变量的函数（一维问题）*1/38太原理工大学

Thermal§3.1数值方法的本质§3.1-1任务1.什么是微分方程的数值解？它是由一组可以构成因变量?分布的数组成的集合，即用一组数字表示待定变量在定义域内的分布。类似于在实验室中进行的实验，仪器的读数构成了所研究区域内被测物理量的分布（有限个离散点的值的集合）。2.实验数据的处理—拟合关系式通常采用的方法是利用实验结果，拟合出遵循的方程，然后利用方程即可计算出任意位置点上因变量的值。*2/38太原理工大学

Thermaleg.设?随x的变化遵循高次多项式的形式利用实验数据，并采用数值方法求得有限数量的多项式中的各项系数值，进而得到相应的?与x遵循的方程式，由此即可求得任意点处的?值了。若最终的兴趣是得到不同位置上的?值，则此方法有些不便，因为各个系数a本身没有什么特别的意义，但要求得?还必须进行代入过程。3.数值方法及任务建立一个把一系列给定点上的?值作为原始未知量的方程，实际上求解微分方程的多数方法均属此类。该方法的任务是提供一组关于这些未知量的代数方程并规定求解这组方程的算法。*3/38太原理工大学

Thermal§3.1-2离散化的概念1.离散化方法及基本思想把注意力集中在网格节点处的值，用离散的值取代包含在微分方程精确解中的连续信息，这样就离散了?的分布，这类数值方法叫离散化方法。根据实际研究对象，把定义域分为若干个有限的区域，在定义域内连续变化的待求变量场，由有限区域上的若干个点的待求变量值来表示，这就是离散化的基本思想。*4/38太原理工大学

Thermal2.离散化方程：所取网格节点上未知因变量?值的代数方程，此方程由支配?的微分方程推导而得。在推导过程中，需对网格节点之间?如何变化作某种假设，变量?在节点间的分布形式不同，推导离散化方程的方法也就不同；另外可以选择在整个计算域内满足一个简单表达式的分布；更为实际的方法还是采用分段分布，即将计算区域分布一定数量的子域或单元，每个子域可以有一个独立的分布假设。*5/38太原理工大学

Thermal§3.1-3离散化方程的结构1.离散化方程的结构一个离散化方程是连接一组网格节点处?值的代数关系式，由支配?的微分方程推导而得，并表示与该微分方程相同的物理信息。当节点数很多时，离散方程的解接近于相应微分方程的精确解，相邻点之间?变化很小，有关?分段分布的细节就不那么重要了。相应于一个已知的微分方程，离散化方程的形式决不是唯一的，这起因于分布假设以及推导方法的不同。网格节点数非常多的极限条件下，所有可能类型的离散化方程将会给出相同的解。*6/38太原理工大学

Thermal2.离散化方法常见的方法主要有：有限差分法和有限元法。两种方法的区别来自于选择分布和推导离散化方程的方法不同。本书主要关注的方法具有有限差分的外形，但它采用了典型的有限元方法所具有的思想，把此方法叫有限差分法可能在于它坚持遵守习惯的有限差分法做法。*7/38太原理工大学

Thermal§3.2推导离散化方程的方法对于一个已知的微分方程，可以用许多方法推导出所要求的离散化方程。§3.2-1泰勒级数公式1.定义：在有限差分法中，通过把控制方程中的各阶导数用相应的差分表达式来代替而形成离散方程。各阶导数的差分表达式可由泰勒级数展开而得，把这种建立离散方程的方法称为泰勒级数展开法。2.差分方程式的建立：ii+1i-2i-1i+2节点i两侧分别有i-2,i-1,i+1,i+2，各节点间距都为h,用泰勒级数展开有：*xh8/38太原理工大学

Thermal取左端及右端的前三项，并进行相加或相减，便可得中心差分的近似式：剩余项的最低阶导数前系数的次数用同样的方法可以得到略去截断误差O(h)的差分计算式：*9/38太原理工大学

Thermal为了提高精度，可以得到截断误差更高阶的差分表达式。3.几点说明①.差分表达式分子项系数的代数和为零；②.各阶导数差分表达式的量纲必须与导数的量纲一致，因而，一阶导数各个差分表达式的分母为?x，二阶为(?x);③.给出一个差分表达式时，必须指明是对哪个点建立的，同样的节点数，不同的建格式的点，导致不同的截断误差，如2对i点只有一阶截差，但对i+1点则是二阶导数具有二阶截差的表达