传热学上机指导书.doc

《传热学上机指导书.doc》由会员分享，可在线阅读，更多相关《传热学上机指导书.doc（8页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、五、计算机实习指导书本指导书是为配合本科生传热学课中计算机应用方面的教学而编写的。应用计算机解决工程实际问题，是现代工程技术人员所必备的技能。在传热学课程中引入计算机实习的目的，是使学生初步掌握用计算机求解传热问题的技能，从而提高学生应用计算机解决工程实际问题的能力，同时也加深对所学习的传热学内容的理解。大量的传热问题能够用计算机求解。研究如何用计算机求解传热问题的专门知识数值传热学（或称计算传热学）已经发展成了传热学的一个分支学科。传热学课中所涉及的只是数值传热学的初步知识。因此，本次计算机实习也仅仅是作为数值传热学的入门。本指导书给出了三个练习题及相应的算法。这三个练习题分别涉及了一维稳态

2、导热、二维稳态导热和一维非稳态导热。要求学生在掌握问题的数值计算方法的基础上，独立编写计算机程序并用所编的程序计算出这三个练习题的数值结果。1 练习题一：一维稳态导热的数值计算11 物理问题图1示出了一个等截面直肋，处于温度t=80的流体中。肋表面与流休之间的对流换热系数为h=45W/m2.，肋基处温度tw=300，肋端绝热。肋片由铝合金制成，其导热系数为=110W/m，肋片厚度为=0.01m，高度为H=0.1m。试计算肋内的温度分布及肋的总换热量。1 2数学描述及其解析解引入无量纲过余温度，则以无量纲温度描述的肋片导热微分方程及其 (1-1) (1-2) (1-3)其中(其中符号含义与教科书

3、杨世铭陶文铨编著传热学相同，以下同)。上述数学模型的解析解为: (1-4)按式(1-4)计算得到的在肋内各点的温度由表1给出。表1 等截面直肋内各点的温度坐标xm00.010.020.030.040.050.060.070.080.090. 1温度t300.00286.56274.84264.76256.13248.97243.19238.11235.58233.70233.091 3数值离散131 区域离散在对方程(1-1)(1-3)进行数值离散之前,应首先进行计算区域的离散。计算区域的离散如图1所示，总节点数取N。132 微分方程的离散由于方程(1-1)在计算区域内部处处成立，因而对图1所

4、示的各离散点亦成立。对任一节点i有：用在节点i的二阶差分代替在节点i的二阶导数,得:整理上式成迭代形式: (i=2,3,N-1) (1-5)133 边界条件离散上面得到的离散方程式(1-5)，对所有内部节点都成立，因此每个内部节点都可得出一个类似的方程。事实上，式(1-5)表达的是一个代数方程组。但这个方程组的个数少于未知数 (i=1,2, ,N)的个数。因此，还需要根据边界条件补充进两个方程后代数方程组才封闭。左边界(x=0)为第一类边界条件，温度为已知，因此可以根据式(1-2)直接补充一个方程为：右边界为第二类边界条件，由图1中边界节点N的向后差分来代替式(1-3)中的导数，得：将此式整理

5、为迭代形式，得：134 最终的离散格式 (i=2,3,N-1) （1-6）135 代数方程组的求解及其程序代数方程组有各种求解方法，较为有效而简便的方法是高斯-赛德尔迭代方法。式（1-6）已给出了代数方程组的迭代形式。在实际计算中，应首先假定一个温度场的初始分布，即给出各节点的温度初值：将这些初值代入方程组（1-6）中进行迭代计算，直至收敛。假如第K步迭代已完成，即为已知，则K+1次迭代的计算式为： (i=2,3,N-1) （1-7）根据式（1-7）编写程序的工作由学生自行完成。计算结果可与解析解比较。 通过数值计算可以发现，由于对肋端的绝热边界条件离散时采用了一阶精度的向后差分，因此当网格数

6、较少时，数值计算结果与解析解在肋端附近有较大的差别。对肋端的绝热边界条件采用具有二阶精度的元体平衡法进行数值离散时，得到以下离散格式：计算可以发现，采用这个格式当网格数较少时，亦能得到精度较高的数值计算结果。2 练习题二：二维稳态导热的数值计算2 1物理问题图2示出了一矩形区域，其边长L=W=1，假设区域内无内热源，导热系数为常数，三个边温度为T1=0，一个边温度为T2=1，求该矩形区域内的温度分布。2 2 数学描述对上述问题的微分方程及其边界条件为：(2-1) (2-2) 作为参考，以下给出该问题的解析解: (2-3)表2列出了由式(2-3)计算得到的在平面区域内各不同位置的温度值。表2 平

7、面区域内温度分布(用式2-3计算得到)坐标温度T坐标温度Txyxy0.250.250.0680.250.750.4320.500.250.0950.500.750.5400.750.250.0680.750.750.4321.000.250.41.000.7530.250.500.1820.251.000.9990.500.500.250.501.001.0240.750.500.1820.751.000.9991.000.5021.001.006523 数值离散2.3.1 区域离散区域离散如图2所示，x方向总节点数为N，y方向总节点数为M，区域内任一节点用i,j表示。2.3.2 方程的离散对

8、于图2中所有的内部节点方程(2-1)都适用，因此可写为：用i,j节点的二阶中心差分代替上式中的二阶导数，得：上式整理成迭代形式：(i=2,3,N-1),(j=2,3,M-1)补充四个边界上的第一类边界条件得：(j=1,2,M)(j=1,2,M)(i=1,2,N)(i=1,2,N)2.4 计算程序计算程序由学生自行完成。3 练习题三：一维非稳态导热的数值计算非稳态导热问题由于有时间变量，其数值计算出现了一些新特点。在非稳态导热微分方程中，与时间因素相关的非稳态项是温度对时间的一阶导数，这给差分离散带来了新的特点。由于这个特点，可以采用不同的方法构造差分方程，从而得到几种不同的差分格式，即所谓的显

9、式、隐式和半隐式。我们仍然从一个具体问题出发来研究非稳态导热问题的数值计算。31 问题一块无限大平板（如图3所示），其一半厚度为L=0.1m，初始温度T0=1000，突然将其插入温度T=20的流体介质中。平板的导热系数=34.89W/m，密度=7800kg/m3，比热c=0.712J/kg，平板与介质的对流换热系数为h=233W/m2.，求平板内各点的温度分布。32 数学描述由于平板换热关于中心线是对称的，仅对平板一半区域进行计算即可。坐标x的原点选在平板中心线上，因而一半区域的非稳态导热的数学描述为：(3-1)(3-2)(3-3)(3-4)该数学模型的解析解为: （3-5）其中,为方程的根，

10、。表3给出了在平板表面(x=L)处由式(3-5)计算得到的不同时刻的温度值。表3平板表面各不同时刻温度值。时 间（S）12345678910温 度()981.84974.47968.88964.20960.11956.14953.08949.97947.07944.3433 数值离散331 计算区域的离散一维非稳态导热指的是空间坐标是一维的。若考虑时间坐标，则所谓的一维非稳态导热实际上是二维问题(见图4)，即：有时间坐标和空间坐标X两个变量。但要注意，时间坐标是单向的，就是说，前一时刻的状态会对后一时刻的状态有影响，但后一时刻的状态却影响不到前一时刻，图4示出了以X和为坐标的计算区域的离散，时

11、间从=0开始，经过一个个时层增加到K时层和K+1时层。332 微分方程的离散对于i节点，在K和K+1时刻可将微分方程(3-1)写成下面式子：(3-6)(3-7)将式(3-6)(3-7)的左端温度对时间的偏导数进行差分离散为：(3-9)(3-8)观察式(3-8)和(3-9)，这两个式子的右端差分式完全相同，但在两个式子中却有不同含义。对式(3-8)，右端项相对i点在K时刻的导数是向前差分。而在式(3-9)中，右端项是I点在K+1时刻的导数的向后差分。将式(3-8)和(3-9)分别代入式(3-6)和(3-7)，并将式(3-6)和(3-7)右端关于x的二阶导数用相应的差分代替，则可得到下列显式和隐式

12、两种不同的差分格式：显式: (3-10) (K=0,1,2, , i=2,3,N-1)全隐式： (3-11) (K=0,1,2, i=2,3,N-1)以上两式中的。从式(3-10)可见，其右端只涉及K时刻的温度，当从K=0（即=0时刻）开始计算时，在K=0时等号右端都是已知值，因而直接可计算出K=1时刻各点的温度。由K=1时刻的各点的温度值，又可以直接利用式(3-10)计算K=2时刻的各点的温度，这样一个时层一个时层的往下推，各时层的温度都能用式(3-10)直接计算出来，不要求解代数方程组。而对于式(3-11) 等号右端包含了与等号左端同一时刻但不同节点的温度，因而必须通过求解代数方程组才能求

13、得这些节点的温度值。333 边界条件的离散对于式(3-3)和(3-4)所给出的边界条件，可以直接用差分代替微分，也可以用元体平衡法给出相应的边界条件，亦有显式和隐式之分。通常，当内部节点采用显式时，边界节点也用显式离散；内部节点用隐式时，边界节点亦用隐式。边界节点的差分格式是显示还是隐式，取决于如何与内部节点的差分方程组合。用K+1时刻相应节点的差分，代替式(3-3)和(3-4)中的微分，可得到边界节点的差分方程： （3-12）334 最终的离散格式显式： （i=1,2,3,N） （初始值） （3-13） (i=2,3,N-1) （3-14） （3-15） （3-16）其中K=0,1,2,。当采用二阶精度的元体平衡法离散时，与式（3-15）和（3-16）对应的离散格式为： 其中。隐式： （初始值） （3-17） （3-18） (i=2,3,N-1) （3-19） （3-20）其中K=0,1,2,在用隐式差分计算时，每个时层都需要迭代求解代数方程组（3-17）式（3-20）。在每个时层计算时，都要先假定一个温度场（一般取上一时层的温度场为本时层的初始场），然后迭代计算直至收敛。34 程序程序由学生自行完成。