《热传导问题的数值.ppt》由会员分享，可在线阅读，更多相关《热传导问题的数值.ppt（30页珍藏版）》请在三一办公上搜索。
1、第三章 导热问题的数值解法,任好玲机电实验大楼B506,一、导热问题数值求解的基本思想及内部节点离散方程的建二、边界节点离散方程的建立及代数方程的求解三、非稳态导热问题的数值解法（自学）四、导热问题数值计算实例,主要内容,数值求解的基本思想及常用的数值求解方法有限差分法节点离散方程的建立泰勒级数展开法与热平衡法。节点离散方程(组)的求解 直接求解;简接求解高斯-赛德尔(Gauss-Seidel)迭代法非稳态导热问题数值求解的有关概念,重点：用热平衡法建立稳态导热问题的离散方程，数值求解的高斯-赛德尔(Gauss-Seidel)迭代法,主要内容,数值求解的基本方法及过程,求解导热问题的三种基本方
2、法：(1)理论分析法；(2)数值计算 法；(3)实验法,三种方法的基本求解过程 所谓理论分析方法，就是在理论分析的基础上，直接对微分方程在给定的定解条件下进行积分，这样获得的解称之为分析解，或叫理论解；数值计算法，把原来在时间和空间连续的物理量的场，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，从而获得离散点上被求物理量的值；并称之为数值解；实验法，就是在传热学基本理论的指导下，采用对所研究对象的传热过程所求量的方法,主要内容,数值求解方法的特点,分析法 能获得所研究问题的精确解，可以为实验和数值计算提供比较依据；局限性很大，对复杂的问题无法求解；分析解具
3、有普遍性，各种情况的影响清晰可见.,数值法：在很大程度上弥补了分析法的缺点，适应性强，特别对于复杂问题更显其优越性；与实验法相比成本低,实验法:是传热学的基本研究方法 适应性不好；费用昂贵.,有限差分法（finite-difference）有限元法（finite-element）边界元法（boundary-element）分子动力学模拟（MD）,导热问题数值求解的基本思想及内部节点离散方程的建立,物理问题的数值求解过程,建立控制方程及定解条件,确定节点（区域离散化）,建立节点物理量的代数方程,设立温度场的迭代初值,求解代数方程,是否收敛,解的分析,改进初场,是,否,导热问题数值求解的基本思想及
4、内部节点离散方程的建立,有限差分法,解方程,并用节点的解的集合(离散值)来代替原物体内的连续温度分布,离散：将连续体用网格分割成有限单元体,取节点：以单元体的中心点代表该单元体,建立节点离散方程：对每一单元体按一定方法，将针对微元体得出的导热微分方程简转化成针对有限单元体的节点离散方程(代数方程组),导热问题数值求解的基本思想及内部节点离散方程的建立,有限差分法,例题：二维矩形域内稳态无内热源，常物性的导热问题,控制方程：是指描写物理问题的微分方程 针对图示的导热问题，它的控制方程（即导热微分方程）为：,其四个边的边界条件为三个边界条件中的一种，三个边界条件为：,控制容积、网格线、节点、界面线
5、、步长,导热问题数值求解的基本思想及内部节点离散方程的建立,有限差分法,例题：二维矩形域内稳态无内热源，常物性的导热问题,建立节点物理量的代数方程（离散方程）,首先划分各节点的类型；其次，建立节点离散方程；最后，代数方程组的形成。,对节点(m,n)的代数方程，当 x=y 时，有,设立迭代初场,:直接解法与迭代解法，传热问题的有限差分法中主要采用迭代法。采用迭代法求解时，需对被求的温度场预先设定一个解，这个解称为初场，并在求解过程中不断改进,导热问题数值求解的基本思想及内部节点离散方程的建立,有限差分法,例题：二维矩形域内稳态无内热源，常物性的导热问题,求解代数方程组,求解时遇到的问题：线性；非
6、线性；收敛性等。,m=1 的左边界上各节点的温度已知外，其余(M-1)N 个节点均需建立离散方程，共有(M-1)N 个方程，则构成一个封闭的代数方程组。,线性代数方程组：代数方程一经建立，其各项系数在整个求解过程中不再变化；,非线性代数方程组：代数方程一经建立，其各项系数在整个求解过程中不断更新。,是否收敛判断：是指用迭代法求解代数方程是否收敛，即本次迭代计算所得之解与上一次迭代计算所得之解的偏差是否小于允许值。,解的分析：通过求解代数方程，获得物体中的温度分布，根据温度场应进一步计算通过的热流量，热应力及热变形等。因此，对于数值分析计算所得的温度场及其它物理量应作详细分析，以获得定性或定量上
7、的结论。,建立离散方程的常用方法,节点类型,这是导热问题数值计算的关键一步。要得出节点的离散方程，首先要了解该节点是哪种类型。,具有对流边界条件的外角顶；具有对流边界条件的平直边界节点；具有对流边界条件和对称绝热角顶；具有绝热边界条件的平直边界节点；具有对流边界条件的内角顶；内部节点。,建立离散方程的常用方法,泰勒级数展开法,对于节点（m+1,n）及（m-1,n）分别写出函数t对（m,n）点的泰勒级数,建立离散方程的常用方法,泰勒级数展开法,中心差分,二维稳态导热的微分方程式为,其在节点（m,n）处的中心差分方程式为,若 x=y 则有,建立离散方程的常用方法,控制容积平衡法(热平衡法),基本思
8、想：对每个有限大小的控制容积（元体）应用能量守恒，从而获得温度场的代数方程组，它从基本物理现象和基本定律出发，不必事先建立控制方程，依据能量守恒和Fourier导热定律即可。,能量守恒：导入元体的总热流量元体内热源生成热 导出元体的总热流量元体内能的增量,从节点(m-1,n)通过界面 w 传导到节点(m,n)的热流量,通过界面 e,n,s 传导给节点（m,n）的热流量也可求得,对元体(m,n).根据能量守恒定律可知：,建立离散方程的常用方法,控制容积平衡法(热平衡法),导入元体的热流量为正，导出元体的热量为负,说明：上述分析与推导是在笛卡儿坐标系中进行的；热平衡法概念清晰，过程简捷；热平衡法与
9、建立微分方程的思路与过程一致，但不同的是前者是有限大小的元体，后者是微元体；适用于非均匀网格；适用于导热系数为温度的函数或内热源分布不均匀的情形；适用于物体内的各个节点（内节点与边界节点）,物理概念清晰，推导过程简单；（2）,边界节点离散方程的建立及代数方程的求解,边界节点离散方程的建立,对于第一类边界条件的热传导问题，处理比较简单，因为已知边界的温度，可将其以数值的形式加入到内节点的离散方程中，组成封闭的代数方程组，直接求解。而对于第二类或第三类边界条件的导热问题，所有内节点的离散方程组成的代数方程组是不封闭的，因未知边界温度，因而应对位于该边界上的节点补充相应的代数方程，才能使方程组封闭，
10、以便求解。为了求解方便，这里我们将第二类边界条件及第三类边界条件合并起来考虑，用qw表示边界上的热流密度或热流密度表达式。为使结果更具一般性，假设物体具有内热源(不必均匀分布)。,边界节点离散方程的建立及代数方程的求解,边界节点离散方程的建立平直边界上的节点,为使讨论具有一般性，设物体具有内热源,平直边界上的节点,当,边界节点离散方程的建立及代数方程的求解,边界节点离散方程的建立外部角点,为使讨论具有一般性，设物体具有内热源,外部角点,当,边界节点离散方程的建立及代数方程的求解,边界节点离散方程的建立内部角点,为使讨论具有一般性，设物体具有内热源,内部角点,当,边界节点离散方程的建立及代数方程
11、的求解,边界节点离散方程的建立边界热流密度,绝热边界,平直边界,当,传入计算区域的热量为正,对流边界,外部角点,内部角点,当计算区域中出现曲线或倾斜边界时,常用阶梯型折线来模拟真实边界,边界节点离散方程的建立及代数方程的求解,边界节点离散方程的建立不规则边界处理,模拟或逼近精度与什么有关？,边界节点离散方程的建立及代数方程的求解,边界代数方程的求解方法,迭代法：先对要计算的场作出假设（设定初场），在迭代计算中不断予以改进，直到计算前的假定值与计算结果相差小于允许值为止的方法，称迭代计算收敛。,直接解法：通过有限次运算获得精确解的方法，如：矩阵求解，高斯消元法。,高斯赛德尔迭代法：每次迭代计算，
12、均是使用节点温度的最新值。雅可比迭代法：每次迭代计算，均用上一次迭代计算出的值。,边界节点离散方程的建立及代数方程的求解,边界代数方程的求解方法,设初场：,边界节点离散方程的建立及代数方程的求解,边界代数方程的求解方法收敛的准则,当有接近于零的t 时，第三个较好,=10-310-6,差分方程能收敛,边界节点离散方程的建立及代数方程的求解,边界代数方程的求解方法,例题1：高斯-赛德尔迭代法：,满足收敛条件,能否收敛？,设初值，一般设为零,边界节点离散方程的建立及代数方程的求解,边界代数方程的求解方法,例题1：高斯-赛德尔迭代法求1、2、3、4温度：,1,2,3,4,t=100,t=100,t=1
13、00,t=500,各节点均为内节，故：,小结,导热问题数值求解的基本思想数值计算是解决较复杂的导热问题的有效途径常用的数值求解方法：有限单元法和有限差分法*离散(分割)、取节点、列节点离散方程、解方程建立节点离散方程的两种方法:泰勒级数展开法与热平衡法。节点离散方程的建立是导热问题数值求解的重要环节(要求能根据不同情况，用热平衡法自行推导稳态导热问题的节点有限差分方程)。节点离散方程（组）的求解 用高斯-赛德尔(Gauss-Seidel)迭代法*求解二维稳态导热问题的有限差分方程组迭代原则、基本求解步骤、迭代方程选取对迭代结果的影响、获得收敛解的条件向前差分、向后差分与中心差分的概念,小结,例题解析,P191，题4-23(1)（要求0.02）(建议步长：x=y=0.25m),节点网格,小结,例题解析,P188，题4-9：,P188，题4-10：迭代方程:,小结,例题解析,
