导热问题数值解法.ppt

《导热问题数值解法.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第三节非稳态导热问题的数值解非稳态导热与稳态导热的主要区别在于控制方程中多了一个非稳态项，而扩散项的离散方法与稳态导热是一样的。因此，本节将重点讨论非稳态项的离散方法以及扩散项离散时所取时间层的不同对计算带来的影响。一维非稳态导热问题的离散x为空间坐标，将计算区域划分为（N-1）等份，得到N个空间节点；两节点之间的距离为称为空间步长；τ为时间坐标，将时间坐标上的计算区域划分为（i-1）等份，得到i个时间节点。从一个时间层到下一个时间层的间隔为Δτ，称为时间步长。空间网格与时间网格的交点，如（n,i），代表了时间—空间

2、区域中一个节点的位置，相应的温度记为。xτn,in+1,in-1,in,i-1n,i+1非稳态项的离散如果将函数t在节点（n,i+1）对点（n,i）作泰勒展开，可有于是有当足够小时，可略而不计，此时非稳态项的差分格式可表示为非稳态项的三种差分格式向前差分向后差分向中心差分一维非稳态导热内节点方程的建立控制方程差分方程化简结果一维非稳态导热内节点方程的建立控制方程差分方程两种差分格式的比较显式差分格式隐式差分格式两种差分格式的区别格式的形式不同计算工作量不同显式格式计算工作量小，隐式格式计算工作量大限制条件不同显式格式

3、对时间步长和空间步长有相互制约的要求，但隐式格式对时间步长及空间步长之间的关系没有任何要求。建立非稳态导热问题节点方程的热平衡法将研究区域离散化对各节点所代表的元体建立能量平衡关系式对非稳态导热问题该能量平衡关系式为从各个方向进入元体的热量之和等于该元体热力学能的变化量整理化简，得到各节点的差分方程无限大平板换热边界上节点方程的建立左图示出了一无限大平板的右侧面的一部分，其右侧面受到周围流体的冷却，表面传热系数为h，流体温度为边界节点为N节点N代表宽度为的元体NN-1N-2x无限大平板换热边界上节点方程的推导从左侧

4、面进入元体的热量从右侧面进入元体的热量元体自身热力学能的改变量根据能量守恒定律整理后得到差分方程的进一步演化考察方程中的定义称为网格傅立叶数定义称为网格毕渥数差分方程可演化为一维非问题导热问题的差分方程内节点差分方程换热边界上的差分方程上述方程都是用显式差分格式表示的数值解求解一维非稳态导热的实例物理模型：设有一块厚度为2的无限大平壁，初始温度为。在初始瞬间将它放置于温度为的流体中，流体与板面间的表面传热系数h为常数。试用数值解法确定在非稳态导热过程中板内的温度分布。x控制方程解：由于问题的对称性，只要研究一半即可，

5、此时，该问题的控制方程为区域离散化将所研究平板的一半N等分，共有N+1个节点，其中节点1在平板中心截面上，节点N在平板右侧面上，如图所示两个节点之间的距离为节点-1与节点2换热情况对称，固有相同的温度时间步长取xn=123N-1N-1差分方程上述问题的差分方程为方程组的求解利用上述方程组，从初始温度出发，即可依次求得第二时间层、第三时间层直到I时间层上的温度分布。至于空间步长及时间步长的选取，原则上步长越小，计算结果越接近于精确解，但是需要的计算机内存及计算时间则大大增加。此外，与的关系还受到显式差分格式稳定性的影响

6、。下面，我们从离散方程的结构来分析，说明稳定性限制的物理意义，再通过数值计算实例予以说明。显式差分格式稳定性分析由内部节点差分方程可见，在节点n上，i+1时刻的温度是在该点i时刻温度的基础上考虑了左右相邻两点温度的影响后得出的。现在，假设相邻两点的温度不变，那么合理的情况是:i时刻节点n的温度越高，则其相继时刻（i+1时刻）的温度也越高；反之，i时刻节点n的温度越低，则其相继时刻的温度也越低。所以，在差分方程中要满足这种合理性的条件，则差分方程中与前面的系数必须保持同方向变化。由于的系数大于零，因此前面的系数也必须大

7、于零。显式差分格式稳定性条件内节点差分方程稳定性条件一维非稳态导热，换热边界上节点差分方程稳定性条件具体计算实例题目：厚的无限大平板受对称的冷却，初始温度℃。在初始瞬间，平板突然被置于温度℃的流体中。已知平板的导热系数,。试用数值法求解其温度分布。取。解：区域离散化，取则采用如图所示的离散方法，计算结果列于下表计算区域离散图xn=12340.010.010.01差分方程上述问题的差分方程为计算结果1760.7273.150.667.2648010032200.991.463.284801001002-328.9220

8、.219.2104801001001001550-109.614860100100100100076543210itn对计算结果的说明从上表可以看出，从这一时刻起出现了这样的情况：各点温度随时间作忽高忽低的波动，并且波动幅度越来越大；某点温度越高反而使相继时刻的温度越低。这种现象是违背热力学第二定律的。因为这意味着，在该时间间隔中，从某一时刻起