- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
其中第五次与第六次相对偏差已小于104 迭代终止。
15.20263565 15.20690891 15..20797723
4-4、试对附图所示的等截面直肋的稳态导热问题用数值方 法求解节点 2,3 的温度。图中
t0 850C,t f 250C, h 30W /(m2.K) .肋高 H=4cm,纵 剖面面积 AL 4cm2 , 导热系数 20W /(m.K) 。
Fo=0.2 Bi=0.1 0.99662 0.994 1.002
百度文库
Fo=0.24 Bi=0.1
Bi=1 0.62945 0.64339 0.97833
Bi=1 0.96514 0.95064 1.01525
Bi=1
Bi=10 0.11866 0.12248 0.96881
Bi=10 0.83889 0.82925 1.01163
t4 t4
20 30
30
15
t4 1/ 4t2 t3 10 5
开始时假设取
t
0
1
t 20
20
℃;
t
0
3
t
0
4
15 ℃
得迭代值汇总于表
迭代次数
0
20
20
15
15
1
26.25
22.8125
21.5625
14.84375
2
28.59375
23.359375
22.109375
15.1171875
时用式(3-19)表示的级数的第一项代替整个级数(计
算中用前六项之和来替代)可能引起 的误差。
解: Bi 0.1 1.0 10
,不同 Bi 下前六个根如下表所示:
μ1
μ2
μ3
μ4
0.3111
3.1731
6.2991
9.4354
0.8603
3.4256
6.4373
9.5293
1.4289
4.3058
7.2281
3
28.8671875 23.49609375 22.24607565 15.18554258
4
28.93554258 23.53027129 22.28027129
5
28.95263565 23.53881782 22.28881782
6
28.9569089 23.54095446 22.290955445
因为迭代公式的选择应使每一个迭代变量的系数总大于或等于式中其他变量的系数绝对值
代数和。
4-3、试对附图所示的常物性,无内热源的二维稳态导热问题用高斯-赛德尔迭代法计算
t1,t2 ,t3 ,t4 之值。
解:温度关系式为:
t1 1/ 4t2 t3 40 30
t t
2 3
1/ 4t1 1/ 4t1
Bi=10
第一项的值 前六项的值
比值
第一项的值 前六项和的值
比值
0.94513 0.94688 0.99814
Fo=0.24 Bi=0.1 0.99277 0.99101 1.00177
0.61108 0.6198 0.98694
x0
Bi=1 0.93698 0.92791 1.00978
0.10935 0.11117 0.98364
2x3
x1
x3
3
x1
x2
假设 x2, x3 初值为 0,迭代结果如下:
迭代次数 0
1
2
3
4
x1
0
2.5
2.625
2.09375
2.6328125
x2
0
-0.75
0.4375 - 1.171875
1.26171825
x3
0
1.25
显然,方程迭代过程发散
-0.0625
2.078125
-0.89453125
10.2003
μ5 12.5743 12.6453 13.2142
μ6 15.7143 15.7713 16.2594
Fo=0.2 及 0.24 时计算结果的对比列于下表:
第一项的值 前六和的值
比值
Fo=0.2 Bi=0.1 0.94879 0.95142 0.99724
第一项的值 前六项和的值
比值
第四章
复习题 1、 试简要说明对导热问题进行有限差分数值计算的基本思想与步骤。 2、 试说明用热平衡法建立节点温度离散方程的基本思想。 3、 推导导热微分方程的步骤和过程与用热平衡法建立节点温度离散方程的过程十分相似,
为什么前者得到的是精确描述,而后者解出的确实近似解。 4、 第三类边界条件边界节点的离散那方程,也可用将第三类边界条件表达式中的一阶导数
解:对于 2 点可以列出:
节点
2:
t1 t 2 x
t3 t4 x
2hx(t1
t2)
0;
节点 3:
t2 t3 x2
h(t f
t1
)
2h
x 2
(t
f
t3) 0
。
由此得:
t1 t2
t3
t2
2hx 2
(t1
t2) 0 t2
,
t3
h
(t f
t3)
x2 h
2
(t
f
t3) 0
用差分公式表示来建立。试比较这样建立起来的离散方程与用热平衡建立起来的离散方 程的异同与优劣。 5.对绝热边界条件的数值处理本章采用了哪些方法?试分析比较之. 6.什么是非稳态导热问题的显示格式?什么是显示格式计算中的稳定性问题? 7.用高斯-塞德尔迭代法求解代数方程时是否一定可以得到收敛德解?不能得出收敛的解 时是否因为初场的假设不合适而造成?
,
t2
t1
t3
2hxH
2
t
f
2
2hxH
2
t3
t2
h tf
hx 2 2
t f
1
h
hx 2 2
hx 2
30 0.022 20 0.01
0.06 ,于是有: t2
t1 t2 0.12t f 2 0.12
,
t3
t2
30 / 20t f 0.03t f
1 30 / 20 0.03
t2 1.5t f 0.03tf =t2 1.53t f
2.53
2.53
,代入得:
2.12t2
t1 t2
1.53t f 2.53
0.12t f , 5.3636t2
Bi=10 0.77311 0.76851 1.00598
4-2、试用数值计算证实,对方程组
x1 2x2 2x3 1
x1
x2
x3
3
2x1 2x2 x3 5
用高斯-赛德尔迭代法求解,其结果是发散的,并分析其原因。
解:将上式写成下列迭代形式
x1 1/ 25 2x2 x3
x2
1/ 21
8.有人对一阶导数 你能否判断这一表达式是否正确,为什么? 一般性数值计算 4-1、采用计算机进行数值计算不仅是求解偏微分方程的有力工具,而且对一些复杂的经验 公式及用无穷级数表示的分析解,也常用计算机来获得数值结果。试用数值方法对
Bi=0.1,1,10 的三种情况计算下列特征方程的根
并用计算机查明,当