第四章复习题1、 试简要说明对导热问题进行有限差分数值计算的基本思想与步骤。
2、 试说明用热平衡法建立节点温度离散方程的基本思想。
3、 推导导热微分方程的步骤和过程与用热平衡法建立节点温度离散方程的过程十分相似,为什么前者得到的是精确描述,而后者解出的确实近似解。
4、 第三类边界条件边界节点的离散那方程,也可用将第三类边界条件表达式中的一阶导数用差分公式表示来建立。
试比较这样建立起来的离散方程与用热平衡建立起来的离散方程的异同与优劣。
5.对绝热边界条件的数值处理本章采用了哪些方法?试分析比较之.6.什么是非稳态导热问题的显示格式?什么是显示格式计算中的稳定性问题?7.用高斯-塞德尔迭代法求解代数方程时是否一定可以得到收敛德解?不能得出收敛的解时是否因为初场的假设不合适而造成?8.有人对一阶导数()()()221,253x t t t xti n i n i n in ∆-+-≈∂∂++你能否判断这一表达式是否正确,为什么? 一般性数值计算4-1、采用计算机进行数值计算不仅是求解偏微分方程的有力工具,而且对一些复杂的经验公式及用无穷级数表示的分析解,也常用计算机来获得数值结果。
试用数值方法对Bi=0.1,1,10的三种情况计算下列特征方程的根:)6,2,1( =n n μ3,2,1,tan ==n Binn μμ并用计算机查明,当2.02≥=δτa Fo 时用式(3-19)表示的级数的第一项代替整个级数(计算中用前六项之和来替代)可能引起的误差。
Bi n n =μμtanFo=0.2及0.24时计算结果的对比列于下表:4-2、试用数值计算证实,对方程组⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=++=++=-+5223122321321321x x x x x x x x x用高斯-赛德尔迭代法求解,其结果是发散的,并分析其原因。
解:将上式写成下列迭代形式()()⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧--=-+=--=2131323213212/1252/1x x x x x x x x x假设3,2xx 初值为0,迭代结果如下:迭代次数 0 1 2 3 41x 0 2.5 2.625 2.09375 2.63281252x 0 -0.75 0.4375 - 1.171875 1.261718253x 0 1.25 -0.0625 2.078125 -0.89453125显然,方程迭代过程发散因为迭代公式的选择应使每一个迭代变量的系数总大于或等于式中其他变量的系数绝对值代数和。
4-3、试对附图所示的常物性,无内热源的二维稳态导热问题用高斯-赛德尔迭代法计算4321,,,t t t t 之值。
解:温度关系式为:()()()()⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧+++=+++=+++=+++=5104/115304/130204/130404/1324413412321t t t t t t t t t t t t 开始时假设取()()200201==t t ℃;()()150403==t t ℃ 得迭代值汇总于表迭代次数0 20 20 15 151 26.25 22.8125 21.5625 14.843752 28.59375 23.359375 22.109375 15.11718753 28.8671875 23.49609375 22.24607565 15.185542584 28.93554258 23.53027129 22.28027129 15.20263565 5 28.95263565 23.53881782 22.28881782 15.206908916 28.9569089 23.54095446 22.290955445 15..20797723 其中第五次与第六次相对偏差已小于410-迭代终止。
4-4、试对附图所示的等截面直肋的稳态导热问题用数值方法求解节点2,3的温度。
图中)./(30,25,852000K m W h C t C t f ===.肋高H=4cm,纵剖面面积,42cm A L =导热系数)./(20K m W =λ。
解:对于2点可以列出:节点2:;0)(2214321=-∆+∆-+∆-t t x h x t t x t t λδλδ节点3:0)(22)(23132=-∆+-+∆-t t xh t t h xt t f f δλδ。
由此得:0)(22122321=-∆+-+-t t x h t t t t λδ,0)(2)(32332=-∆+-+-t t hx t t h t t f f λδλδδ,⎪⎭⎫ ⎝⎛∆+⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛∆++=λδλδ22222312xH h t xH h t t t f⎪⎪⎭⎫ ⎝⎛∆++⎥⎦⎤⎢⎣⎡⎪⎪⎭⎫ ⎝⎛∆++=λδλλδλ2122223x h h t x h t ht t f f06.001.02002.03022=⨯⨯=∆λδx h ,于是有:12.0212.0212+++=f t t t t , ()53.253.153.203.05.103.020/30103.020/302f 223ff f f t t t t t t t t t +++=++++==,代入得:f ft t t t t 12.053.253.112.2212++=,f 2123036.053.153.23636.5tt t t t f +++=, f t t t 8336.153.23636.412+=,3636.48336.153.22ff t t t +=,Ct 8.5979.593636.484.4505.2153636.4258336.18553.22≅=+=⨯+⨯=,Ct 8.3875.3853.22553.18.593≅=⨯+=。
离散方程的建立4-5、试将直角坐标中的常物性无内热源的二维稳态导热微分方程化为显式差分格式,并指出其稳定性条件()y x ∆≠∆。
解:常物性无内热源二维非稳态方程微分方程为⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=∂∂2222y t x t a tτ扩散项取中心差分,非稳态项取向前差分:()()⎪⎪⎭⎫ ⎝⎛∆+-+∆+-=∆--+-++211211122y t t t x t t t a t t in i n i n i n i n i n i n i n τ所以有()i nin i n i n t y x a t t y x a t ⎥⎦⎤⎢⎣⎡⎪⎪⎭⎫ ⎝⎛∆+∆∆+++⎪⎪⎭⎫ ⎝⎛∆+∆∆=-++2211221112111ττ 稳定性条件2/1≤+∆∆y x Fo Fo4-6、极坐标中常物性无内热源的非稳态导热微分方程为⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂=∂∂2222211ϕτt r r t r r t a t试利用本题附图中的符号,列出节点(i,j )的差分方程式。
解:将控制方程中的各阶导数用相应的差分表示式代替,可得:1111111122222112k k i k k k k k k k k j t jt j t j j t j j i j i j i j j j t t t t t t t t t r a r r r r τϕ++-+---⎛⎫--+--+=+⋅+⋅ ⎪ ⎪∆∆∆∆⎝⎭,,,,,,,,,。
也可采用热平衡法。
对于图中打阴影线的控制容积写出热平衡式得:()()111k k k k k k i j i ji j i ji j i jjj j t t t t t t r r c r r r r ϕρλλτϕϕ+-+---∆∆=∆+∆+∆∆∆,,,,,,11,22kkkki j i ji j i j j j ttt t r r r r r r λϕλϕ+---∆∆⎛⎫⎛⎫+∆+-∆ ⎪ ⎪∆∆⎝⎭⎝⎭,,,对等式两边同除以jr r ϕ∆∆并简化,可以得出与上式完全一样相同的结果。
4-7、一金属短圆柱在炉内受热厚被竖直地移植到空气中冷却,底面可以认为是绝热的。
为用数值法确定冷却过程中柱体温度的变化,取中心角为1rad 的区域来研究(如本题附图所示)。
已知柱体表面发射率,自然对流表面传热系数,环境温度,金属的热扩散率,试列出图中节点(1,1),(M,1)(M,n)及(M,N )的离散方程式。
在r 及z 方向上网格是各自均分的。
解:应用热平衡法来建立四个节点点离散方程。
节点(1,1):2121211211111111222282k k k k k k t t t t t t r r z r z c z r λλρτ+⎡⎤⎛⎫---∆∆∆∆∆⎛⎫⎛⎫+⋅=⋅⎢⎥⎪ ⎪ ⎪∆∆∆⎝⎭⎝⎭⎢⎥⎝⎭⎣⎦,,,,,,节点(m ,1):()11,11112111;22222k k k k k k kk m m m m m m m m m m m m t t t t t t t t z z z z z r r r r c r r r r z λλλρτ-+----∆∆∆∆∆⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫-+++⋅∆=⋅∆⋅ ⎪⎪ ⎪⎪ ⎪∆∆∆∆⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭,,,,,,,节点(m ,n ):()()1,4411111003332242242422k k k k k k m n m nm n m n m n m n m m m m m m m m m m n m n tt t t t t r r r r r r r z r z r r z r r h t t T T c r z rλλεσρ---------+⎡-⎤+∆∆∆∆∆∆∆⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎡⎤-+++⨯-+-= ⎪⎪ ⎪ ⎪ ⎪ ⎪⎪ ⎪ ⎪ ⎪⎢⎥⎣⎦∆∆∆⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎣⎦,,,,,,,。
4-8、一个二维物体的竖直表面收液体自然对流冷却,为考虑局部表面传热系数的影响,表面传热系数采用25.11)(t t c h -=来表示。
试列出附图所示的稳态无内热源物体边界节点(M,n )的温度方程,并对如何求解这一方程提出你的看法。
设网格均分。
解:利用热平衡法:()()0.25M n f M n fh c t t tt =--,,,将h 写为()()0.25M n f M n f h c t t t t =--,,,其中M n t ,为上一次迭代值,则方程即可线性化。
4-9、在附图所示的有内热源的二维导热区域中,一个界面绝热,一个界面等温(包括节点4),其余两个界面与温度为f t的流体对流换热,h 均匀,内热源强度为Φ 。
试列出节点1,2,5,6,9,10的离散方程式。
解:节点1:()512111102242f t t t t x y x y yh t t y x λλφ--∆∆⎛⎫⎛⎫++∆∆-∆-= ⎪ ⎪∆∆⎝⎭⎝⎭;节点2:()32621210222t t t t t t y y x x y x x y λλλφ---∆∆⎛⎫⎛⎫++∆+∆∆= ⎪ ⎪∆∆∆⎝⎭⎝⎭; 节点5:()()159565510222f t t t t t t y x y x y yh t t y y x λλλφ---∆∆⎛⎫⎛⎫++∆+∆∆-∆-= ⎪ ⎪∆∆∆⎝⎭⎝⎭; 节点6:()()()()2676105560t t t t t t t tx y x y x y y x y x λλλλφ----∆+∆+∆+∆+∆∆=∆∆∆∆;节点9:()5910991022422f t t t t x y x y x y h t t y x λλφ--∆∆∆∆⎛⎫⎛⎫⎛⎫++∆∆-+-= ⎪ ⎪ ⎪∆∆⎝⎭⎝⎭⎝⎭; 节点10:()()9101110610101222f t t t t t t y y x x y xh h t x x y λλλφ---∆∆⎛⎫⎛⎫++∆+∆∆-∆-= ⎪ ⎪∆∆∆⎝⎭⎝⎭。