还剩16页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
第四章复习题
1、试简要说明对导热问题进行有限差分数值计算的基本思想与步骤
1、试说明用热平衡法建立节点温度离散方程的基本思想
1、推导导热微分方程的步骤和过程与用热平衡法建立节点温度离散方程的过程十分相似,___前者得到的是精确描述,而后者解出的确实近似解
1、第三类边界条件边界节点的离散那方程,也可用将第三类边界条件表达式中的一阶导数用差分公式表示来建立试比较这样建立起来的离散方程与用热平衡建立起来的离散方程的异同与优劣5.对绝热边界条件的数值处理本章采用了哪些方法?试分析比较之.6.什么是非稳态导热问题的显示格式?什么是显示格式计算中的稳定性问题?7.用高斯-塞德尔迭代法求解代数方程时是否一定可以得到收敛德解?不能得出收敛的解时是否因为初场的假设不合适而造成?8.有人对一阶导数你能否判断这一表达式是否正确,___?一般性数值计算4-
1、采用计算机进行数值计算不仅是求解偏微分方程的有力工具,而且对一些复杂的经验公式及用无穷级数表示的分析解,也常用计算机来获得数值结果试用数值方法对Bi=
0.1110的三种情况计算下列特征方程的根并用计算机查明,当时用式(3-19)表示的级数的第一项代替整个级数(计算中用前六项之和来替代)可能引起的误差解,不同Bi下前六个根如下表所示Biμ1μ2μ3μ4μ5μ
60.
10.
31113.
17316.
29919.
435412.
574315.
71431.
00.
86033.
42566.
43739.
529312.
645315.
7713101.42__
4.
30587.
228110.
200313.
214216.2594Fo=
0.2及
0.24时计算结果的对比列于下表Fo=
0.2Bi=
0.1Bi=1Bi=10第一项的值
0.
948790.
629450.11866前六和的值
0.
951420.
643390.12248比值
0.
997240.
978330.96881Fo=
0.2Bi=
0.1Bi=1Bi=10第一项的值
0.
996620.
965140.838__前六项和的值
0.
9940.
950640.82925比值
1.
0021.
015251.01163Fo=
0.24Bi=
0.1Bi=1Bi=10第一项的值
0.
945130.
611080.10935前六项的值
0.
946880.
61980.11117比值
0.
998140.
986940.98364Fo=
0.24Bi=
0.1Bi=1Bi=10第一项的值
0.
992770.
936980.77311前六项和的值
0.
991010.
927910.76851比值
1.
001771.
009781.005984-
2、试用数值计算证实,对方程组用高斯-赛德尔迭代法求解,其结果是发散的,并分析其原因解将上式写成下列迭代形式假设初值为0,迭代结果如下迭代次数
0123402.
52.
6252.
093752.63281250-
0.
750.4375-
1.
1718751.
2617182501.25-
0.
06252.078125-
0.__453125显然,方程迭代过程发散因为迭代公式的选择应使每一个迭代变量的系数总大于或等于式中其他变量的系数绝对值代数和4-
3、试对附图所示的常物性,无内热源的二维稳态导热问题用高斯-赛德尔迭代法计算之值解温度关系式为开始时假设取℃;℃得迭代值汇总于表迭代次数
020201515126.
2522.
812521.
562514.
84375228.
5937523.
35937522.
10937515.
1171875328.
867187523.
4960937522.
2460756515.
18554258428.
9355425823.
5302712922.
2802712915.
20263565528.
9526356523.
5388178222.
2888178215.20690__
1628.95690__
23.
5409544622.
29095544515..20797723其中第五次与第六次相对偏差已小于迭代终止4-
4、试对附图所示的等截面直肋的稳态导热问题用数值方法求解节点2,3的温度图中.肋高H=4cm纵剖面__导热系数解对于2点可以列出节点2节点3由此得,,,于是有,,代入得,,,,,离散方程的建立4-
5、试将直角坐标中的常物性无内热源的二维稳态导热微分方程化为显式差分格式,并指出其稳定性条件(解常物性无内热源二维非稳态方程微分方程为扩散项取中心差分,非稳态项取向前差分所以有稳定性条件4-
6、极坐标中常物性无内热源的非稳态导热微分方程为试利用本题附图中的符号,列出节点(ij)的差分方程式解将控制方程中的各阶导数用相应的差分表示式代替,可得也可采用热平衡法对于图中打阴影线的控制容积写出热平衡式得对等式两边同除以并简化,可以得出与上式完全一样相同的结果4-
7、一金属短圆柱在炉内受热厚被竖直地移植到空气中冷却,底面可以认为是绝热的为用数值法确定冷却过程中柱体温度的变化,取中心角为1rad的区域来研究(如本题附图所示)已知柱体表面发射率,自然对流表面传热系数,环境温度,金属的热扩散率,试列出图中节点(1,1),(M1)Mn及(MN)的离散方程式在r及z方向上网格是各自均分的解应用热平衡法来建立四个节点点离散方程节点(1,1)节点(m,1)节点(m,n)4-
8、一个二维物体的竖直表面收液体自然对流冷却,为考虑局部表面传热系数的影响,表面传热系数采用来表示试列出附图所示的稳态无内热源物体边界节点(Mn)的温度方程,并对如何求解这一方程提出你的看法设网格均分解利用热平衡法,将h写为,其中为上一次迭代值,则方程即可线性化4-
9、在附图所示的有内热源的二维导热区域中,一个界面绝热,一个界面等温(包括节点4),其余两个界面与温度为的流体对流换热,h均匀,内热源强度为试列出节点1,2,5,6,9,10的离散方程式解节点1;节点2;节点5;节点6;节点9;节点10当以上诸式可简化为节点1;节点2;节点5节点6;节点9;节点10一维稳态导热计算4-
10、一等截面直肋,高H厚,肋根温度为,流体温度为,表面传热系数为h肋片导热系数为将它均分成4个节点(见附图),并对肋端为绝热及为对流边界条件(h同侧面)的两种情况列出节点2,3,4的离散方程式设H=45cm=50W/m.K℃,℃,计算节点2,3,4的温度(对于肋端的两种边界条件)解采用热平衡法可列出节点
2、
3、4的离散方程为节点2;节点3;节点4肋端绝热,肋端对流其中将已知条件代入可得下列两方程组肋端绝热肋端对流由此解得肋端绝热,,;肋端对流,,肋端对流换热的条件使肋端温度更接近于流体温度4-
11、复合材料在____及化工等工业中日益得到广泛的应用附图所示为双层圆筒壁,假设层间接触紧密,无接触热阻存在已知W/m.K,℃,℃,试用数值方法确定稳态时双层圆筒壁截面上的温度分布解采用计算机求解,答案从略采用热平衡法对两层管子的各离散区域写出能量方程,进行求解;如果采用Taylor展开法列出方程,则需对两层管子单独进行,并引入界面上温度连续及热流密度连续的条件,数值计算也需分两区进行,界面耦合截面的温度分布定性地示于上图中4-
12、有一水平放置的等截面直杆,根部温度℃,其表面上有自然对流散热,,其中,d为杆直径,杆高H=10cm,直径d=1cm=50W/m.K,℃不计辐射换热试用数值方法确定长杆的散热量(需得出与网格无关的解杆的两端可认为是绝热的解数值求解过程略,Q=
2.234W4-13在上题中考虑长杆与周围环境的辐射换热,其表面发射率为
0.8,环境可作为温度为的大空间,试重新计算其导热量解数值求解过程略,Q=
3.320W4-
14、有如附图所示的一抛物线肋片,表面形线方程为肋根温度及内热源恒定,流体表面传热系数h流体温度为常数定义试
(1)建立无量纲温度的控制方程;
(2)在无量纲参数下对上述控制方程进行数量计算确定无量纲温度的分布解无量纲温度方程为数值计算结果示于下图中,无量纲温度从肋根的1变化到肋端的
0.852一维非稳态导热计算4-
15、一直径为1cm长4cm的钢制圆柱形肋片,初始温度为25℃,其后,肋基温度突然升高到200℃,同时温度为25℃的气流横向掠过该肋片,肋端及两侧的表面传热系数均为100试将该肋片等分成两段(见附图),并用有限差分法显式格式计算从开始加热时刻起相邻4个时刻上的温度分布(以稳定性条件所允许的时间间隔计算依据)已知=43W/m.K,(提示节点4的离散方程可按端面的对流散热与从节点3到节点4的导热相平衡这一条件列出)解三个节点的离散方程为节点2节点3节点4以上三式可化简为稳定性要求,即,代入得,如取此值为计算步长,则,于是以上三式化成为时间点12340200252525△
200128.8125252△
200128.
8155.
8055.093△
200137.
9573.
6472.544△
200143.
0486.
7085.30在上述计算中,由于之值正好使,因而对节点2出现了在及2时刻温度相等这一情况如取为上值之半,则,,,于是有对于相邻四个时层的计算结果如下表所示时间点12340200252525△
20076.9125252△
200102.
8632.
7032.533△
200116.
9842.
6342.234△
200125.
5152.
5751.944-
16、一厚为
2.54cm的钢板,初始温度为650℃,后置于水中淬火,其表面温度突然下降为
93.5℃并保持不变试用数值方法计算中心温度下降到450℃所需的时间已知建议将平板8等分,取9个节点,并把数值计算的结果与按海斯勒计算的结果作比较解数值求解结果示于下图中随着时间步长的缩小,计算结果逐渐趋向于一个恒定值,当=
0.00001s时,得所需时间为
3.92s如图所示,横轴表示时间步长从1秒,
0.1秒,
0.01秒,
0.001秒,
0.0001秒,
0.00001秒的变化;纵轴表示所需的冷却时间(用对数坐标表示)4-
17、一火箭燃烧器,壳体内径为400mm厚10mm壳体内壁上涂了一层厚为2mm的包裹层火箭发动时,推进剂燃烧生成的温度为3000℃的烟气,经燃烧器端部的喷管喷住大气大气温度为30℃设包裹层内壁与燃气间的表面传热系数为2500W/m.K,外壳表面与大气间的表面传热系数为350,外壳材料的最高允许温度为1500℃试用数值法确定为使外壳免受损坏,燃烧过程应在多长时间内完成包裹材料的=
0.3W/m.K,a=解采用数值方法解得4-
18、锅炉汽包从冷态开始启动时,汽包壁温随时间变化为控制热应力,需要计算汽包内壁的温度场试用数值方法计算当汽包内的饱和水温度上升的速率为1℃/min3℃/min时,启动后10min20min及30min时汽包内壁截面中的温度分布及截面中的最大温差启动前,汽包处于100℃的均匀温度汽包可视为一无限长的圆柱体,外表面绝热,内表面与水之间的对流换热十分强烈汽包的内径外半径热扩散率解数值方法解得部分结果如下表所示汽包壁中的最大温差,K启动后时间,min温升速率,K/min
13107.
13621.
41209.
46328.
393010.
1930.574-
19、有一砖墙厚为,=
0.85W/m.K,室内温度为℃,h=6起初该墙处于稳定状态,且内表面温度为15℃后寒潮入侵,室外温度下降为℃,外墙表面传热系数如果认为内墙温度下降
0.1℃是可感到外界温度起变化的一个定量判据,问寒潮入侵后多少时间内墙才感知到?解采用数值解法得t=7900s4-
20、一冷柜,起初处于均匀的温度(20℃)后开启压缩机,冷冻室及冷柜门的内表面温度以均匀速度18℃/h下降柜门尺寸为保温材料厚8cm,=
0.02W/m.K冰箱外表面包裹层很薄,热阻可忽略而不计柜门外受空气自然对流及与环境之间辐射的加热自然对流可按下式计算其中H为门高表面发射率通过柜门的导热可看作为一维问题处理试计算压缩机起动后2h内的冷量损失解取保温材料的,用数值计算方法得冷量损失为4-
21、一砖砌墙壁,厚度为240mm,=
0.81W/m.K设冬天室外温度为24h内变化如下表所示室内空气温度℃且保持不变;外墙表面传热系数为10,内墙为6试用数值方法确定一天之内外墙,内墙及墙壁中心处温度随时间的变化取设上述温度工况以24h为周期进行变化时刻/h00010020030040050060070080090010001100温度/-
5.9-
6.2-
6.6-
6.7-
6.8-
6.9-
7.2-
7.7-
7.6-
7.0-
4.9-
2.3时刻/h120013001400150016001700180019002000210022002300温度/-
1.
02.
41.
81.
81.
60.5-
1.6-
2.8-
3.5-
4.3-
4.8-
5.3解采用数值解法得出的结果如下表所示时刻/h012345678环境温度/-
5.9-
6.2-
6.6-
6.7-
6.8-
6.9-
7.2-
7.7-
7.6外墙温度/-
1.70-
2.19-
2.44-
2.76-
2.85-
2.93-
3.01-
3.26-
3.67墙壁中心温度/
3.
653.
323.
152.
922.
872.
812.
752.
592.31内墙温度/
8.
998.
828.
738.
618.
588.
558.
528.
438.28时刻/h91011121314151617环境温度/-7-
4.9-
2.3-
12.
41.
81.
81.
60.5外墙温度/-
3.58-
3.07-
1.
340.
781.
874.
634.
154.
143.97墙壁中心温度/
2.
362.
703.
875.
326.
057.
957.
627.
627.51内墙温度/
8.
318.
499.
119.
8710.
2611.
2611.
1011.
1011.10时刻/h181920212223环境温度/-
1.6-
2.8-
3.5-
4.3-
4.8-
5.3外墙温度/
3.
061.
340.36-
0.22-
0.87-
1.29墙壁中心温度/
6.
095.
735.
054.
664.
213.93内墙温度/
10.
7110.
109.
739.
539.
309.14__稳态导热问题4-
22、如附图所示,一矩形截面的空心电流母线的内外表面分别与温度为的流体发生对流换热,表面传热系数分别为,且各自沿周界是均匀的,电流通过壁内产生均匀热源今欲对母线中温度分布进行数值计算,试
(1)划出计算区域
(1)对该区域内的温度分布列出微分方程式及边界条件;
(1)对于图中内角顶外角顶及任一内部节点列出离散方程式(),设母线的导热系数为常数4-
23、一个长方形截面的冷空气通道的尺寸如附图所示假设在垂直于纸面的方向上冷空气及通道墙壁的温度变化很小,可以忽略试用数值方法计算下列两种情况下通道壁面的温度分布及每米长度上通过壁面的冷量损失
(1)内外壁分别维持在10℃及30℃
(1)内外壁与流体发生对流换热,且有℃,,℃,解此题应采用计算机求解如有墙角导热的热点模拟实验设备,则计算参数(如h,及网格等)可以取得与实验设备的参数相一致,以把计算结果与实测值作比较根据对称性,取1/4区域为计算区域数值计算解出,对于给定壁温的情形,每米长通道的冷损失为
39.84W,对于第三类边界条件为
30.97W(取壁面导热系数)内外表面为给定壁温时等温线分布如下图所示第三类边界条件的结果定性上类似4-
24、为了提高现代燃气透平的进口燃气温度以提高热效率,在燃气透平的叶片内部开设有冷却通道以使叶片金属材料的温度不超过允许值,为对叶片中的温度分布情况作一估算,把附图a所示的截片形状简化成为附图b所示的情形已知,试计算
(1)截面中最高温度及其位置;
(2)单位长度通道上的热量解根据对称性选择1/4区域为计算区域,采用网格,取壁面时得单位长度的传热量为
987.8W,等温线分布如图所示截面中最高温度发生在左上角,该处温度为
1419.9综合分析与分析、论述题4-
25、工业炉的炉墙以往常用红砖和耐火砖组成由于该两种材料的导热系数较大,散热损失较严重,为了节省能量,近年来国内广泛采用在耐火砖上贴一层硅酸纤维毡,如附图所示今用以下的非稳态导热简化模型来评价黏贴硅酸纤维毡的收益设炉墙原来处于与环境平衡的状态,s时内壁表面突然上升到550℃并保持不变这一非稳态导热过程一直进行到炉墙外表面的对流,辐射热损失与通过墙壁的导热量相等为止在炉墙升温过程中外表面的总表面传热系数由两部分组成,即自然对流引起的部分及辐射部分其中为外表面温度,为内表面温度,为简化计算,设三种材料的导热系数分别为W/m.K,W/m.K,W/m.K试计算每平方炉墙每平方__上由于粘贴了硅酸纤维毡而在炉子升温过程中节省的能量解采用数值计算方法,详细过程从略4-
26、空气在附图所示的一长方形截面的送风管道中作充分发展的层流流动,其z方向的动量方程简化为而且上式可看成是源项为的一常物性导热方程试用数值方法求解这一方程并计算fRe之值f为阻力系数,Re为特征长度为当量直径计算时可任取一个值,并按a/b=
0.5及1两种情形计算解假设壁温为常数,则不同a/b下换热充分发展时的fRe及Nu数的分析解为a/bNufRe
12.
98570.
53.39624-
27、一家用烤箱处于稳定运行状态,箱内空气平均温度℃,气体与内壁间的表面传热系数外壁面与20℃的周围环境间的表面传热系数烤箱保温层厚30mmW/m.K保温层两侧的护板用金属制成且很薄,分析中可不予考虑,然后,突然将烤箱调节器开大,风扇加速,内壁温度突然上升到185℃,设升温过程中烤箱外壁面与环境间的表面传热系数可用计算,环境温度仍保持为20℃,为烤箱外壁面温度,c之值与运行时一样试确定烤箱内壁温度跃升后到达新的稳定状态所需时间解需采用数值方法求解,过程从略小论文题目4-
28、一厚为
2.54cm的钢管,初始温度为16℃其后,温度为572℃的液态金属突然流过管内,并经历了10s液态金属与内壁面间的表面传热系数h=
2.84钢管可以按平壁处理,其外表面的散热由对流及辐射两条路径,并分别可按及计算,,周围环境温度=20℃试用有限差分法确定在液态金属开始流入后的18s时截面上的温度分布已知钢管的41W/m.K,,c=536J/kg.K解在钢管壁厚方向上取27个点,以内壁为坐标原点,沿着壁厚方向为x正方向,数值计算结果如下位置/cm
00.
10.
20.
30.
40.
50.
60.
70.
80.9温度/
216.
0215.
6214.
6213.
0210.
7207.
9204.
6200.
8196.
6192.1位置/cm
11.
11.
21.
31.
41.
51.
61.
71.
81.9温度/
187.
3182.
3177.
2172.
0166.
9161.
8157.
0152.
5148.
2144.4位置/cm
22.
12.
22.
32.
42.
52.54温度/
141.
0138.
1135.
7133.
9132.
6132.
0131.9用图形表示如下4-
29、为对两块平板的对接焊过程(见附图a)进行计算,对其物理过程作以下简化处理钢板中的温度场仅是x及时间的函数;焊枪的热源作用在钢板上时钢板吸收的热流密度,为电弧有效加热半径,为最大热流密度;平板上下表面的散热可用计算,侧面绝热;平板的物性为常数,熔池液态金属的物性与固体相同;固体熔化时吸收的潜热折算成当量的温升值,即如设熔化潜热为L,固体比热容为c则当固体达到熔点后要继续吸收相当于使温度升高(L/c)的热量,但在这一吸热过程中该温度不变这样,附图a所示问题就简化为附图b所示的一维稳态导热问题试
(1)列出该问题的数学描写;
(2)计算过程开始后
3.4s内钢板中的温度场,设在开始的
0.1s内有电弧的加热作用已知,h=
12.6=
41.9W/m.KL=255kJ/kg℃,H=12cm解取初始温度与环境温度均为该问题的数学描写为0xHt0;;0;0为了更好分辨热源附近的温度场宜采用非均分网格计算得出开始加热后的
3.4s内钢板中的温度分布如下图所示4-
30、在壁厚为7cm的铸铁模型中铸造14cm厚的黄铜板设此问题可按一维问题处理,试确定达到铜版完全凝固所需的时间计算时作以下简化处理液体铜在瞬间内充满形腔;液体铜及铸型的初始温度各自均匀;液体铜内无自然对流,固液体铜内均为导热;液体铜与固体铜的物性相同且为常数;铸件与铸型之间接触良好,不存在空气隙;铸型外两表面与周围环境间的散热可用表示;液体铜在固定的凝固点下凝固,凝固过程中释放出的熔化潜热可折算成相当于使物体温度升高(L/c)的热量,但在潜热释放过程中该温度应一直保持为经过这样一番简化后所计算的问题变为如附图所示的双层平板的一维导热问题试
(1)列出该问题的数学描写;
(2)在下列条件下计算使钢板完全凝固所需的时间已知铸型初温℃,液体铜初温为1100℃,℃h=4,W/m.K,W/m.K,,L=
167.5kJ/kg℃解设铸型厚为,铸件半厚为,则有0x,t0,,0,,,,数值计算结果得出所需时间为
304.9s4-
31、建筑物采暖的一种方式是在__地板下设置热空气通道,如附图所示设地板下的水泥混凝土层的一侧绝热,地面温度℃热空气通道截面尺寸为150mm并在混凝土层中对称布置,通道壁温保持为80℃试计算单位长度热空气通道的传热量,并从计算结果中整理出此种情形下形状因子S之值解单位长度传热量为
122.5W,形状因子为S=
3.101m4-
32、试用数值方法确定如附图所示圆管外正方形翅片的肋效率已知=40mm翅片厚W/m.K据文献〔10〕分析,此时肋效率可以画成的曲线形成,并以为参数这里是一假想半径,以为半径的圆的__等于所研究翅片的__在h=10~100的范围内进行计算,并把结果表示成的曲线解计算结果如下图所示4-33.有一块印制电路板如附图a所示.中间为
0.8mm厚的铜板,导热系数为,其两侧为玻璃纤维环氧树脂板层,铜板底端被冷却到40℃,其他三个侧面可以认为绝热,金属板上__的发热元件及其功耗如图所示.假定通过玻璃纤维环氧树脂板层的散热可以不计,试用数值计算确定铜板中的温度分布.根据元件确定的网格划分示于附图b中.。