数值传热学习题集

合集下载

传热学习题(附参考答案)

传热学习题(附参考答案)

传热学习题(附参考答案)一、单选题(共56题,每题1分,共56分)1.下列不能提高对流传热膜系数的是( )。

A、利用多管程结构;B、增大管径;C、在壳程内装折流挡板;D、冷凝时在管壁上开一些纵槽。

正确答案:B2.传热过程中当两侧流体的对流传热系数都较大时,影响传热过程的将是( )A、管壁热阻;B、管内对流传热热阻;C、污垢热阻;D、管外对流传热热阻;正确答案:C3.在一输送系统中,改变离心泵的出口阀门开度,不会影响( )A、泵的工作点B、管路特性曲线C、管路所需压头D、泵的特性曲线正确答案:D4.用l20℃的饱和蒸汽加热原油,换热后蒸汽冷凝成同温度的冷凝水,此时两流体的平均温度差之间的关系为(Δtm)并流( )(Δtm)逆流A、大于B、小于C、不定D、等于正确答案:D5.橡胶与塑料和纤维比较,正确的是( )。

A、模量最大B、强度最大C、结晶度最大D、Tg最低正确答案:D6.某单程列管式换热器,水走管程呈湍流流动,为满足扩大生产需要,保持水的进口温度不变的条件下,将用水量增大一倍,则水的对流传热膜系数为改变前的( )A、2倍B、1.149倍C、1.74倍D、不变正确答案:C7.按照离心泵完好标准,轻石脑油返输用离心泵机械密封允许泄漏量( )。

A、允许每分钟5滴B、比液化气的允许泄漏量多5滴C、允许每分钟15滴D、允许每分钟10滴正确答案:D8.对于列管式换热器,当壳体与换热管温度差( )时,产生的温度差应力具有破坏性,因此需要进行热补偿A、大于50℃B、大于60℃C、大于45℃D、大于55℃正确答案:A9.下列哪个选项是离心泵的主要部件( )A、电机B、轴封装置和轴向力平衡装置C、密封环D、叶轮和泵壳正确答案:D10.蒸汽中不凝性气体的存在,会使它的对流传热系数α值( )A、升高B、降低C、都可能D、不变正确答案:B11.以乙烯为原料经催化剂催化聚合而得的一种热聚性化合物是( )A、PB、PC、PVCD、PP正确答案:B12.关闭出口阀启动离心泵的原因是( )A、轴功率最大B、能量损失最小C、处于高效区D、启动电流最小正确答案:D13.若被输送液体的黏度增大时,离心泵的效率( )。

传热学习题(附参考答案)

传热学习题(附参考答案)

传热学习题(附参考答案)一、单选题(共56题,每题1分,共56分)1.安装在管路中的阀门( )A、需考虑流体方向B、不必考虑流体方向C、不必考虑操作时的方便D、不必考虑维修时的方便正确答案:A2.用水蒸气在列管换热器中加热某盐溶液,水蒸气走壳程。

为强化传热,下列措施中最为经济有效的是( )A、增大换热器尺寸以增大传热面积B、减少传热壁面厚度C、改单管程为双管程D、在壳程设置折流挡板正确答案:C3.下列不能提高对流传热膜系数的是( )。

A、增大管径;B、利用多管程结构;C、冷凝时在管壁上开一些纵槽。

D、在壳程内装折流挡板;正确答案:A4.离心泵是依靠离心力对流体作功,其作功的部件是( )。

A、电机B、泵轴C、叶轮D、泵壳正确答案:C5.对于活化能越大的反应,速率常数随温度变化越 ( )A、不确定B、大C、小D、无关正确答案:B6.有机玻璃是指( )。

A、聚甲基丙烯酸甲酯B、聚苯乙烯C、聚乙烯D、聚氯乙烯正确答案:A7.水在无相变时在圆形管内强制湍流,对流传热系数αi为1000W /(m2·℃)若将水的流量增加1倍,而其他条件不变,则αi为( )A、不变B、500C、2000D、1741正确答案:D8.不属于换热器检修内容的是( )A、清扫管束和壳体B、管束焊口、胀口处理及单管更换C、检查修复管箱、前后盖、大小浮头、接管及其密封面,更换垫片D、检查校验安全附件正确答案:D9.间歇反应器的一个生产周期不包括( )A、出料时间B、反应时间C、加料时间D、设备维修时间正确答案:D10.可在器内设置搅拌器的是( )换热器A、套管B、釜式C、热管D、夹套正确答案:D11.流体流量突然减少,会导致传热温差( )。

A、始终不变B、下降C、升高D、变化无规律正确答案:B12.下列哪个选项不是列管换热器的主要构成部件。

( )A、封头B、管束C、外壳D、蛇管正确答案:D13.列管式换热器一般不采用多壳程结构,而采用( )以强化传热效果A、翅片板B、折流挡板C、隔板D、波纹板正确答案:B14.裂解气深冷分离过程中采用的主要方法是( )A、精馏法B、特殊精馏法C、萃取法D、吸附法正确答案:A15.下面关于裂解气分离流程说法正确的是( )A、一套乙烯装置采用哪种流程,主要取决于流程对所需处理裂解气的适应性、能量消耗、运转周期及稳定性、装置投资等几个方面。

(完整版)传热学题库

(完整版)传热学题库

5. 温度边界层:对流换热时,在传热壁面附近形成的一层温度有很大变化(或温
度变化率很大)的薄层
四、简答题
1.导热问题常见的三类边界条件(文字叙述或用数学式表达)。 答:第一类边界条件:τ>0,tw=fw(x, τ)
第二类边界条件:τ>0, t n w
fw (x, )
第三类边界条件:τ>0,
对流换热表面传热系数 h=40 W/m2·K,铝的物性参数为ρ=2700kg/m3,c=0.9
kJ/kg·K,λ=240W/m·K。试用集总参数法确定该铝球由 400℃降至 100℃所需的
时间。(忽略辐射换热)
一、单项选择题
1. 温度梯度的方向是指向温度
A.增加方向 B.降低方向 C. 不变方向 D. 趋于零方向
q
1
1000 0.8
3.将初始温度为 400℃,重量为 40g 的铝球突然抛入 15℃的空气中。已知 对流换热表面传热系数 h=40 W/m2·K,铝的物性参数为ρ=2700kg/m3,c=0.9 kJ/kg·K,λ=240W/m·K。试用集总参数法确定该铝球由 400℃降至 100℃所需的 时间。(忽略辐射换热)
4. 实际物体的辐射力比同温度下黑体的辐射力 (黑度 <1)
A.大 B.小 C.一样 D.差不多
5.管内对流换热的流态判别是用
A. Gr B. Re C. Pe
D. Gr·Pr
6.导热系数 是物性参数下面说法正确的是 p6
A.与材料种类有关 B.与材料温度有关 C. 与材料种类、温度都
有关
7.灰体的吸收比与投射辐射的波长分布 p374
A.减少导热
B.减小对流换热
C.减少对流与辐射换热

传热学习题及参考答案

传热学习题及参考答案

《传热学》复习题一、判断题1.稳态导热没有初始条件。

()2.面积为A的平壁导热热阻是面积为1的平壁导热热阻的A倍。

()3.复合平壁各种不同材料的导热系数相差不是很大时可以当做一维导热问题来处理()4.肋片应该加在换热系数较小的那一端。

()5.当管道外径大于临界绝缘直径时,覆盖保温层才起到减少热损失的作用。

()6.所谓集总参数法就是忽略物体的内部热阻的近视处理方法。

()7.影响温度波衰减的主要因素有物体的热扩散系数,波动周期和深度。

()8.普朗特准则反映了流体物性对换热的影响。

()9. 傅里叶定律既适用于稳态导热过程,也适用于非稳态导热过程。

()10.相同的流动和换热壁面条件下,导热系数较大的流体,对流换热系数就较小。

()11、导热微分方程是导热普遍规律的数学描写,它对任意形状物体内部和边界都适用。

( )12、给出了边界面上的绝热条件相当于给出了第二类边界条件。

( )13、温度不高于350℃,导热系数不小于0.12w/(m.k)的材料称为保温材料。

( )14、在相同的进出口温度下,逆流比顺流的传热平均温差大。

( )15、接触面的粗糙度是影响接触热阻的主要因素。

( )16、非稳态导热温度对时间导数的向前差分叫做隐式格式,是无条件稳定的。

( )17、边界层理论中,主流区沿着垂直于流体流动的方向的速度梯度零。

( )18、无限大平壁冷却时,若Bi→∞,则可以采用集总参数法。

( )19、加速凝结液的排出有利于增强凝结换热。

( )20、普朗特准则反映了流体物性对换热的影响。

( )二、填空题1.流体横向冲刷n排外径为d的管束时,定性尺寸是。

2.热扩散率(导温系数)是材料指标,大小等于。

3.一个半径为R的半球形空腔,空腔表面对外界的辐射角系数为。

4.某表面的辐射特性,除了与方向无关外,还与波长无关,表面叫做表面。

5.物体表面的发射率是ε,面积是A,则表面的辐射表面热阻是。

6.影响膜状冷凝换热的热阻主要是。

传热部分习题集

传热部分习题集

第十章 传热部分习题集一、选择题1、单层平壁导热速率方程式为δλt A Q ∆=,方程式中Δt 是指:………[ ] 牛顿冷却定律t A Q ∆=α中的Δt 是指……………………………………[ ]A 、热流体之间的平均温度差B 、平壁两侧表面上的温度差C 、流体进出口温度差D 、固体壁面与流体主体之间的温度差2、在平壁的稳定热传导中,已知两层的导热发21λλ>,则两层平壁的导热速率Q 1和Q 2必满足如下关系: ……………………………………………………………[ ]A 、Q 1>Q 2B 、Q 1<Q 2C 、Q 1=Q 2D 、无法判断3、当流体的进出口温度一定时,在列管式换热器中,冷热液体分别采用逆流、折流和并流操作,比较其传热推动力的在大小,正确的顺序是…………………………[ ]A.Δt m 并>Δt m 逆>Δt m 折, B.Δt m 逆>Δt m 折>Δt m 并C. Δt m 逆>Δt m 并>Δt m 折D.Δt m 并>Δt m 折>Δt m 逆4、套管式换热器的内管介质为空气,套管环隙介质为饱和蒸汽,如果蒸汽压力和空气进口温度一定时,当空气流量增加时,则:(1)总传热系数K 将:……………………………………………………………[ ](2)空气的出口温度将:…………………………………………………………[ ]A 、增大B 、减少C 、基本不变D 、无法确定。

5、有四种液体,其中适宜在列管式壳程中流动………………………………[ ]A 、不清洁的和易结垢的流体B 、腐蚀性的流体C 、被冷却的流体D 、压力高的流体6、若传热系数K 满足公式:o i K αα111+=,当0α>>i α时,为了提高K 值,最有效的途径是……………………………………………………………………………[ ]A 、0α提高B 、i α提高C 、0α,i α同时提高D 、不能确定8、由两层不同材料构成的等厚平壁,已知λ1 >λ2 , Δt 1与Δt 2 R 1与R 2,Q 1与Q 2之间的关系为 :…………………………………………………………………………[ ]A、Q1<Q2 R 1>R2 Δt 1=Δt 2B、Q1=Q2 R 1>R2 Δt 1>Δt 2C、Q1>Q2 R 1<R2 Δt 1<Δt 2D、Q1=Q2 R1<R2Δt1<Δt29、在间壁式换热器中,用100℃的饱和水蒸汽将空气由20℃加热到60℃,换热器的管壁温度接近于:……………………………………………………………………[ ]A、40℃B、60℃C、100℃D、80℃10、为了减少室外设备的热损失,保温层外包有一层金属皮应该是:………[ ]A、表面光滑、色泽较浅。

4传热(习题及解答2007版)

4传热(习题及解答2007版)

四传热单层平壁导热4.1 红砖平壁墙,厚度为500 mm,一侧温度为200 ℃,另一侧温度为30 ℃,设红砖的平均导热系数可取0.57 W/(m·℃),试求:(1)单位时间,单位面积导过的热量q为多少W/m2,并换算成工程单位制kCal/ (m2·h) 。

(2)距离高温侧350 mm处的温度为多少?(3)如红砖导热系数与温度关系可表示为:λ=0.51+5×10 -4t,则q又为多少?多层平壁导热4.2 某燃烧炉的炉墙由三种砖依次砌成:第一层为耐火砖,厚b1=0.23 m,λ1=1.05 W/(m·℃);第二层为绝热砖,λ2 =0.151 W/(m·℃);第三层为普通砖,b3 =0.24 m,λ3 =0 .93 W/(m·℃)。

若已知耐火砖内侧温度为1000 ℃,耐火砖与绝热砖接触处温度为940 ℃,绝热砖与普通砖接触面温度不得超过138℃;试求:(1)绝热砖层的厚度;(2)普通砖外侧温度。

4.3 平壁炉壁由三种材料组成,其厚度和导热系数如下:序号材料厚度b(mm)导热系数λ(W/m·℃)1(内层)耐火砖200 1.072 绝热砖100 0.143 钢板 6 45若耐火砖层内表面温度t1=1150 ℃,钢板外表面温度t4 =30 ℃,试计算导热的热通量。

又实测通过炉壁的热损失为300 W/m2,如计算值与实测不符,试分析原因并计算附加热阻。

多层圆筒导热4.4 一外径为100 mm的蒸汽管,外包一层50 mm绝热材料A,λA =0.07 W/(m·℃),其外再包一层25 mm的绝热材料B,λB =0.087 W/(m·℃)。

设A的内侧温度为170 ℃,B外侧温度为38 ℃。

试求每米管上的热损失及A、B界面的温度。

4.5 Φ60×3 mm铝合金管(其导热系数可取为45 W/(m·℃)),外包一层厚30 mm石棉,之外再包一层厚30 mm的软木,石棉和软木的导热系数分别为0.16 W/(m·℃)和0.04 W/(m·℃) 。

数值传热学 习题答案

数值传热学 习题答案

数值传热学习题答案数值传热学习题答案数值传热学是热力学的一个重要分支,主要研究热量在物质中传递的机理和规律。

在实际工程中,我们经常会遇到各种与传热有关的问题,通过数值计算可以得到准确的答案。

下面我将为大家提供一些数值传热学习题的答案,希望能够帮助大家更好地理解和应用这门学科。

1. 一个铝制热交换器的表面积为10平方米,其表面温度为100摄氏度,环境温度为20摄氏度。

已知铝的导热系数为200 W/(m·K),求热交换器的传热速率。

答:根据传热定律,传热速率与传热面积、传热系数和温度差之间成正比。

传热速率 = 传热系数× 传热面积× 温度差。

将已知数据代入公式中,可得传热速率= 200 × 10 × (100 - 20) = 160,000 W。

2. 一个房间的尺寸为5米× 5米× 3米,墙壁和天花板的厚度为0.2米,墙壁和天花板的导热系数为0.5 W/(m·K),室内温度为25摄氏度,室外温度为10摄氏度。

求房间的传热损失。

答:房间的传热损失可以通过计算墙壁和天花板的传热速率来得到。

墙壁和天花板的传热速率 = 传热系数× 传热面积× 温度差。

墙壁和天花板的传热面积 = 2 × (5 × 5) + 2 × (5 × 3) = 70平方米。

将已知数据代入公式中,可得墙壁和天花板的传热速率= 0.5 × 70 × (25 - 10) = 525 W。

因此,房间的传热损失为525瓦特。

3. 一个水箱的体积为1立方米,初始温度为20摄氏度,水的密度为1000千克/立方米,比热容为4186 J/(千克·摄氏度),水箱的表面积为2平方米,表面温度为100摄氏度。

已知水的传热系数为0.6 W/(m^2·K),求水箱内水的温度随时间的变化。

数值传热学习题答案(汇总版)

数值传热学习题答案(汇总版)

2-4-9
= rP rS
式(2-4-9)也可以写成 a PTP = a E TE + aW TW + b 的形式。而且两种结果是一致的。
2—6:
n n TE −TW dT P , n = 解:将 , dx 2x n n TE −2TPn + TW d 2T P , n = , dx2 x 2
dk = f (x ) 代入原方程,得: dx

2-4-4
rk rk a E = , aW = , a P = a E + aW , b x w x e
= SrP r ,
式(2-4-4)可以写成 a PTP = a E TE + aW TW + b 的形式。 2. 再用 Taylor 展开法导出 k
2 2 uE + uP u = , 2 2 e
2 2 uW + uP u = 2 2 w
t u ut N − uP y = (y ) , n n
t
t ut u p − uS y = (y ) 。 s s
t
(y ) n = (y ) s = y
n n n n TE −TW TE −2TPn + TW k + f (x ) +S=0 整理得: 2x x 2
4kT P= 2k + xf ( x)T E+2k − xf ( x)T W +2x 2 S
− 2k 时, a E 会成为负值, x 2k 当 f(x)> 时, aW 会成为负值。 x
rk dr = rk r r dr dr dr
w
e
1 d

数值传热学部分习题答案

数值传热学部分习题答案

习题4-2一维稳态导热问题的控制方程:022=+∂∂S xTλ 依据本题给定条件,对节点2节点3采用第三类边界条件具有二阶精度的差分格式,最后得到各节点的离散方程: 节点1: 1001=T节点2: 1505105321-=+-T T T 节点3:75432=+-T T求解结果:852=T ,403=T对整个控制容积作能量平衡,有:02150)4020(15)(3=⨯--⨯=∆+-=∆+x S T T h x S q f f B即:计算区域总体守恒要求满足习题4-5在4-2习题中,如果25.03)(10f T T h -⨯=,则各节点离散方程如下:节点1: 1001=T节点2: 1505105321-=+-T T T节点3:25.03325.032)20(4015])20(21[-⨯+=-⨯++-T T T T对于节点3中的相关项作局部线性化处理,然后迭代计算; 求解结果:818.822=T ,635.353=T (迭代精度为10-4)迭代计算的Matlab 程序如下: x=30; x1=20;while abs(x1-x)>0.0001a=[1 0 0;5 -10 5;0 -1 1+2*(x-20)^(0.25)]; b=[100;-150; 15+40*(x-20)^(0.25)]; t=a^(-1)*b; x1=x; x=t(3,1);endtcal=t习题4-12的Matlab程序%代数方程形式A i T i=C i T i+1+B i T i-1+D imdim=10;%计算的节点数x=linspace(1,3,mdim);%生成A、C、B、T数据的基数;A=cos(x);%TDMA的主对角元素B=sin(x);%TDMA的下对角线元素C=cos(x)+exp(x); %TDMA的上对角线元素T=exp(x).*cos(x); %温度数据%由A、B、C构成TDMAcoematrix=eye(mdim,mdim);for n=1:mdimcoematrix(n,n)=A(1,n);if n>=2coematrix(n,n-1)=-1*B(1,n);endif n<mdimcoematrix(n,n+1)=-1*C(1,n);endend%计算D矢量D=(coematrix*T')';%由已知的A、B、C、D用TDMA方法求解T%消元P(1,1)=C(1,1)/A(1,1);Q(1,1)=D(1,1)/A(1,1);for n=2:mdimP(1,n)=C(1,n)/(A(1,n)-B(1,n)*P(1,n-1));Q(1,n)=(D(1,n)+B(1,n)*Q(1,n-1))/(A(1,n)-B(1,n)*P(1,n-1));end%回迭Tcal(1,mdim)=Q(1,mdim);for n=(mdim-1):-1:1Tcal(1,n)=P(1,n)*Tcal(1,n+1)+Q(1,n);endTcom=[T;Tcal];%绘图比较给定T值和计算T值plot(Tcal,'r*')hold onplot(T)结果比较如下,由比较可知两者值非常切合(在小数点后8位之后才有区别):习题4-14充分发展区的温度控制方程如下:)(1rTr r r x T uc p ∂∂∂∂=∂∂λρ 对于三种无量纲定义w b w T T T T --=Θ、∞∞--=ΘT T T T w 、w w T T T T --=Θ∞进行分析如下1)由wb wT T T T --=Θ得:w w b T T T T +Θ-=)(由T 可得:x T x T x T T T x T w b w w b ∂∂Θ-+∂∂Θ=∂+Θ-∂=∂∂)1(])[(rT r T T r T T T r T w w b w w b ∂∂Θ-+∂Θ∂-=∂+Θ-∂=∂∂)1()(])[( 由b T 与r 无关、Θ与x 无关以及x T ∂∂、rT∂∂的表达式可知,除了w T 均匀的情况外,该无量纲温度定义在一般情况下是不能用分离变量法的; 2)由∞∞--=ΘT T T T w 得: ∞∞+Θ-=T T T T w )(由T 可得:xT x T T T x T w w ∂∂Θ=∂+Θ-∂=∂∂∞∞])[(rT r T T r T T T r T w w w ∂∂Θ+∂Θ∂-=∂+Θ-∂=∂∂∞∞∞)(])[( 由b T 与r 无关、Θ与x 无关以及x T ∂∂、rT∂∂的表达式可知,在常见的四种边界条件中除了轴向及周向均匀热流const q w =的情况外,有0=∂∂rT w,则该无量纲温度定义是可以用分离变量法的; 3)由wwT T T T --=Θ∞得: w w T T T T +Θ-=∞)(由T 可得:xT x T T T x T w w w ∂∂Θ-=∂+Θ-∂=∂∂∞)1(])[(rT r T T r T T T r T w w w w ∂∂Θ-+∂Θ∂-=∂+Θ-∂=∂∂∞∞)1()(])[( 同2)分析可知,除了轴向及周向均匀热流const q w =的情况外,有0=∂∂rT w,该无量纲温度定义是可以用分离变量法的;习题4-181)采用柱坐标分析,写出统一的稳态柱坐标形式动量方程:r r x x w r v r r r u x ∂∂+∂∂∂∂=∂∂+∂∂+∂∂1)()(1)(1)(φλφρθφρφρx 、r 和θ分别是圆柱坐标的3个坐标轴,u 、v 和w 管内的流动方向;对于管内的层流充分发展有:0=v 、0=w ,0=∂∂xu; 并且x 方向的源项:x pS ∂∂-=r 方向的源项:r pS ∂∂-= θ方向的源项:θ∂∂-=pr S 1 由以上分析可得到圆柱坐标下的动量方程: x 方向: 0)(1)(1=∂∂-∂∂∂∂+∂∂∂∂x pu r r r u r r r θλθλ r 方向:0=∂∂r pθ方向:0=∂∂θp边界条件: R r =,0=u0=r ,0=∂∂r u ;对称线上,0=∂∂θu 不考虑液体的轴向导热,并简化分析可以得到充分发展的能量方程为:)(1)(1θλθλρ∂∂∂∂+∂∂∂∂=∂∂Tr r r T r r r x T uc p 边界条件: R r =,w q r T =∂∂λ;0=r ,0=∂∂rTπθ/0=,0=∂∂-θλT2)定义无量纲流速:dxdp R uU 2-=λ并定义无量纲半径:R r /=η;将无量纲流速和无量纲半径代入x 方向的动量方程得:0))1((1))1((122=∂∂-∂-∂∂∂+∂-∂∂∂xp U dx dp R R R R U dx dp R RR R θληλθηηλληηη 上式化简得:01)1(1)(1=+∂∂∂∂+∂∂∂∂θηθηηηηηU U 边界条件:1=η,0=U0=η,0=∂∂ηU ;对称线上,0=∂∂θU定义无量纲温度:λ/0R q T T b-=Θ其中,0q 是折算到管壁表面上的平均热流密度,即:Rq q wπ=0; 由无量纲温度定义可得:b T Rq T +Θ=λ将T 表达式和无量纲半径η代入能量方程得:)(1)(100θληλθηηλληηηρ∂Θ∂∂∂+∂Θ∂∂∂=∂∂R q R R R R q R R R x T uc b p 化简得:)1(1)(10θηθηηηηηρ∂Θ∂∂∂+∂Θ∂∂∂=∂∂x T u c q R b p (1)由热平衡条件关系可以得:mm m b m p b p p RU U q R u u R q A u u dx dT A u c x T u c x T uc 020221221)(===∂∂=∂∂ππρρρ 将上式代入式(1)可得:)1(1)(12θηθηηηηη∂Θ∂∂∂+∂Θ∂∂∂=m U U 边界条件:0=η,0=∂Θ∂η;1=η,R q q w πη10==∂Θ∂0=θ,0=∂Θ∂θ;πθ=,0=∂Θ∂θ单值条件: 由定义可知:0/0=-=ΘλR q T T b b b 且: ⎰⎰Θ=ΘAAb UdAUdA即得单值性条件:0=Θ⎰⎰AA UdAUdA 3)由阻力系数f 及Re 定义有:228)(21/Re ⎪⎭⎫ ⎝⎛=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=D D U D u u dx dp D f e m e m me νρ 且:m W b m W b m W R q T T D T T q Nu ,0,,0~2)/(2Θ=-=-=λλ5-21.一维稳态无源项的对流-扩散方程如下所示:xx u 22∂∂Γ=∂∂φφρ (取常物性)边界条件如下:L L x x φφφφ====,;,00上述方程的精确解如下:11)/(00--=--⋅PeL x Pe L e e φφφφ Γ=/uL Pe ρ 2.将L 分成20等份,所以有:∆=P Pe 201 2 3 4 5 6 ………… …………… 17 18 19 20 21 对于中心差分、一阶迎风、混合格式和QUICK 格式分别分析如下: 1) 中心差分中间节点: 2)5.01()5.01(11-∆+∆++-=i i i P P φφφ 20,2 =i2) 一阶迎风中间节点: ∆-∆++++=P P i i i 2)1(11φφφ 20,2 =i3) 混合格式当1=∆P 时,中间节点:2)5.01()5.01(11-∆+∆++-=i i i P P φφφ20,2 =i当10,5=∆P 时,中间节点: 1-=i i φφ 20,2 =i 4) QUICK 格式*12111)35(8122121⎥⎦⎤⎢⎣⎡---++++++=+--∆∆-∆∆+∆i i i i i i i P P P P P φφφφφφφ 2≠i *1111)336(8122121⎥⎦⎤⎢⎣⎡--++++++=+-∆∆-∆∆+∆i i i i i i P P P P P φφφφφφ 2=i数值计算结果与精确解的计算程序如下:%except for HS, any other scheme doesnt take Pe<0 into consideration %expression of exact solutiony=dsolve('a*b*Dy=c*D2y','y(0)=y0,y(L)=yL','x')y=subs(y,'L*a*b/c','t')y=simple(subs(y,'a*b/c*x','t*X'));ysim=simple(sym(strcat('(',char(y),'-y0)','/(yL-y0)')))y=sym(strcat('(',char(ysim),')*(yL-y0)','+y0'))% in the case of Pe=0y1=dsolve('D2y=0','y(0)=y0,y(L)=yL','x')y1=subs(y1,'-(y0-yL)/L*x','(-y0+yL)*X')%grid Pe numbertt=[1 5 10];%dimensionless lengthm=20;%mdim is the number of inner nodemdim=m-1;X=linspace(0,1,m+1);%initial value of variable during calculationy0=1;yL=2;%cal exact solutionfor n=1:size(tt,2)t=m*tt(1,n);if t==0yval1(n,:)=eval(y1);elseyval1(n,:)=eval(y);endend%extra treatment because max number in MATLAB is 10^308if max(isnan(yval1(:)))yval1=yval1';yval1=yval1(:);indexf=find(isnan(yval1));for n=1:size(indexf,1)if rem(indexf(n,1),size(X,2))==0yval1(indexf(n),1)=yL;elseyval1(indexf(n),1)=y0;endendyval1=reshape(yval1,size(X,2),size(yval1,1)/size(X,2));yval1=yval1';end%CD solutiond=zeros(size(tt,2),mdim);a=repmat([1],size(tt,2),mdim);for n=1:size(tt,2)t=tt(1,n);b(n,:)=repmat([0.5*(1-0.5*t)],1,mdim);c(n,:)=repmat([0.5*(1+0.5*t)],1,mdim);d(n,1)=0.5*(1+0.5*tt(1,n))*y0;d(n,mdim)=0.5*(1-0.5*tt(1,n))*yL;endc(:,1)=0;b(:,mdim)=0;%numerical cal by using TDMA subfuctionyval2=TDMA(a,b,c,d,mdim);yval2=[repmat([1],size(tt,2),1),yval2,repmat([2],size(tt,2),1)]; Fig(1,X,yval1,yval2,tt);title('CD Vs. Exact Solution')% FUS solutiond=zeros(size(tt,2),mdim);a=repmat([1],size(tt,2),mdim);for n=1:size(tt,2)t=tt(1,n);b(n,:)=repmat([1/(2+t)],1,mdim);c(n,:)=repmat([(1+t)/(2+t)],1,mdim);d(n,1)=(1+tt(1,n))/(2+tt(1,n))*y0;d(n,mdim)=1/(2+tt(1,n))*yL;endc(:,1)=0;b(:,mdim)=0;%numerical cal by using TDMA subfuctionyval3=TDMA(a,b,c,d,mdim);yval3=[repmat([1],size(tt,2),1),yval3,repmat([2],size(tt,2),1)]; Fig(2,X,yval1,yval3,tt);title('FUS Vs. Exact Solution')% HS solutiond=zeros(size(tt,2),mdim);a=repmat([1],size(tt,2),mdim);for n=1:size(tt,2)t=tt(1,n);if t>2b(n,:)=repmat([0],1,mdim);c(n,:)=repmat([1],1,mdim);d(n,1)=y0;elseif t<-2b(n,:)=repmat([1],1,mdim);c(n,:)=repmat([0],1,mdim);d(n,mdim)=yL;elseb(n,:)=repmat([0.5*(1-0.5*t)],1,mdim);c(n,:)=repmat([0.5*(1+0.5*t)],1,mdim);d(n,1)=0.5*(1+0.5*t)*y0;d(n,mdim)=0.5*(1-0.5*t)*yL;endendc(:,1)=0;b(:,mdim)=0;% numerical cal by using TDMA subfuctionyval4=TDMA(a,b,c,d,mdim);yval4=[repmat([1],size(tt,2),1),yval4,repmat([2],size(tt,2),1)]; Fig(3,X,yval1,yval4,tt);title('HS Vs. Exact Solution')%QUICK Solutiond=zeros(size(tt,2),mdim);a=repmat([1],size(tt,2),mdim);for n=1:size(tt,2)t=tt(1,n);b(n,:)=repmat([1/(2+t)],1,mdim);c(n,:)=repmat([(1+t)/(2+t)],1,mdim);d(n,1)=(1+tt(1,n))/(2+tt(1,n))*y0;d(n,mdim)=1/(2+tt(1,n))*yL;endc(:,1)=0;b(:,mdim)=0;%numerical cal by using TDMA subfuctionyval5=zeros(size(tt,2),mdim);yval5com=yval5+1;counter=1;%iterativewhile max(max(abs(yval5-yval5com)))>10^-10if counter==1yval5com=TDMA(a,b,c,d,mdim);endfor nn=1:size(tt,2)for nnn=1:mdimif nnn==1d(nn,nnn)=((6*yval5com(nn,nnn)-3*y0-3*yval5com(nn,nnn+1))*tt(1,nn))/(8*(2+tt(1, nn)))+((1+tt(1,nn))/(2+tt(1,nn))*y0);elseif nnn==2d(nn,nnn)=((5*yval5com(nn,nnn)-3*yval5com(nn,nnn+1)-yval5com(nn,nnn-1)-y0)*tt (1,nn))/(8*(2+tt(1,nn)));elseif nnn==mdimd(nn,nnn)=((5*yval5com(nn,nnn)-3*yL-yval5com(nn,nnn-1)-yval5com(nn,nnn-2))*tt (1,nn))/(8*(2+tt(1,nn)))+(1/(2+tt(1,nn))*yL);elsed(nn,nnn)=((5*yval5com(nn,nnn)-3*yval5com(nn,nnn+1)-yval5com(nn,nnn-1)-yval5 com(nn,nnn-2))*tt(1,nn))/(8*(2+tt(1,nn)));endendendyval5=TDMA(a,b,c,d,mdim);temp=yval5;yval5=yval5com;yval5com=temp;counter=counter+1;endyval5=yval5com;yval5=[repmat([1],size(tt,2),1),yval5,repmat([2],size(tt,2),1)];Fig(4,X,yval1,yval5,tt);title('QUICK Vs. Exact Solution')%-------------TDMA SubFunction------------------function y=TDMA(a,b,c,d,mdim)%form a b c d resolve yval2 by using TDMA%eliminationp(:,1)=b(:,1)./a(:,1);q(:,1)=d(:,1)./a(:,1);for n=2:mdimp(:,n)=b(:,n)./(a(:,n)-c(:,n).*p(:,n-1));q(:,n)=(d(:,n)+c(:,n).*q(:,n-1))./(a(:,n)-c(:,n).*p(:,n-1));end%iterativey(:,mdim)=q(:,mdim);for n=(mdim-1):-1:1y(:,n)=p(:,n).*y(:,n+1)+q(:,n);end%-------------ResultCom SubFunction------------------ function y=ResultCom (a,b,c)for n=1:max(size(c,2))y(2*n-1,:)=a(n,:);y(2*n,:)=b(n,:);end%-------------Fig SubFunction------------------ function y=Fig(n,a,b,c,d)figure(n);plot(a,b);hold onplot(a,c,'*');str='''legend(';for n=1:size(d,2)if n==size(d,2)str=strcat(str,'''''Pe=',num2str(d(1,n)),''''')''');elsestr=strcat(str,'''''Pe=',num2str(d(1,n)),''''',');endendeval(eval(str));精确解与数值解的对比图,其中边界条件给定10=φ,2=L φ。

数值传热学第二章部分习题参考答案

数值传热学第二章部分习题参考答案

习题2-4 [解]1.先用控制容积积分法得出离散方程: 以r 乘式01=+⎪⎭⎫⎝⎛S dr dT rk dr d r ,并对图2-2所示的控制容积P 作积分: wewe dr dT rk dr dT rk dr dr dT rk dr d r r⎪⎭⎫⎝⎛-⎪⎭⎫ ⎝⎛=⎪⎭⎫ ⎝⎛⎰1 2-4-1 ()E P e eT T dT dr r δ-⎛⎫=⎪⎝⎭ 2-4-1-1 ()P W w wT T dT dr r δ-⎛⎫= ⎪⎝⎭2-4-1-2将式(2-4-1-1)、式(2-4-1-2)代入式(2-4-1)可以得到:()()W P wP E e ewT T x rk T T x rk dr dr dT rk dr d r r-⎪⎭⎫⎝⎛--⎪⎭⎫ ⎝⎛=⎪⎭⎫ ⎝⎛⎰δδ1 2-4-2 222ee wP w r r rSdr S Sr r -==∆⎰2-4-3根据式(2-4-2)、式(2-4-3)可以得到:P E W P e w e w rk rk rk rk T T T Sr r x x x x δδδδ⎡⎤⎛⎫⎛⎫⎛⎫⎛⎫+=++∆ ⎪ ⎪ ⎪ ⎪⎢⎥⎝⎭⎝⎭⎝⎭⎝⎭⎣⎦2-4-4令e E x rk a ⎪⎭⎫ ⎝⎛=δ,wW x rk a ⎪⎭⎫⎝⎛=δ,W E P a a a +=,P b Sr r =∆,式(2-4-4)可以写成b T a T a T a W W E E P P ++=的形式。

2. 再用Taylor 展开法导出022=++S drdTr k dr T d k的离散方程。

将点E T 对点P T 作Taylor 展开,有:()() +++=!2222e Pe P P E x dr Td x drdTT T δδ 2-4-5再将点W T 对点P T 作Taylor 展开,有:()() ++-=!2222wPw P P W x dr Td x drdTT T δδ 2-4-6根据式(2-4-5)、式(2-4-6)可以计算出dr dT ,22drTd ()[]()()[]()[]()()[]()()[]222222w e e w We P w e E w x x x x T x T x x T x dr dT δδδδδδδδ+---= 2-4-7 ()()()[]()()[]()()[]()222222we e w We P w e E w x x x x T x T x x T x dr T d δδδδδδδδ+++-= 2-4-8 将式(2-4-7)、式(2-4-8)代入上面的非守恒型方程,整理成(并考虑到常物性、均分网格):222P P P P E W P kr kr kr k k T T T r rS r r r ⎡⎤⎡⎤=++-+∆⎢⎥⎢⎥∆∆∆⎣⎦⎣⎦2-4-9 令12e P E kr ra k r r ⎡⎤=+=⎢⎥∆∆⎣⎦,12w P W kr r a k r r⎡⎤=-=⎢⎥∆∆⎣⎦,W E P a a a +=,P b r rS =∆式(2-4-9)也可以写成b T a T a T a W W E E P P ++=的形式。

数值传热学习题集

数值传热学习题集

简答题集锦1.流动与传热数值模拟的基本任务是什么?(把原来在时间域及空间域上连续的物理量的场,如速度场和压力场,用一系列有限个离散点上的变量值的集合来代替,通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组,然后求解代数方程组获得场变量的近似值CFD可以看做是在流动基本方程(质量守恒方程飞动量守恒方程、能量守恒方程)控制下对流动的数值模拟。

通过这种数值模拟,我们可以得到极其复杂问题的流场内各个位置上的基本物理量(如速度、压力、温度、浓度等)的分布,以及这些物理量随时间的变化情况,确定旋涡分布特性、空化特性及脱流区等。

)2.数值模拟过程如何实现,主要步骤是那些?(建模、网格划分、坐标系、数学方程、求解、后处理)a.建立反映工程问题或物理过程本质的数学模型;b.选择与计算区域的边界相适应的坐标系;c.建立网格;d.建立离散方程;e.求解代数方程组;f.后处理,显示计算结果3.建立离散方程有哪些主要方法?比较说明各种方法的优缺点?(有限差分、有限体积、有限元、有限分析等)4什么叫控制方程?常见的控制方程有哪几个?各用在什么场合?5试写出控制方程的通用形式,并说明通用形式中各项的意义?(写明通式,以及各个方程中通式的表达形式)6推导x 方向的动量控制方程中的源项u S 的表达式。

由此证明当密度和黏度为常数时,u S 变为0。

X 方向N-S 方程:Mx S xw z u z x v y u y divu x u x x p Dt Du +∂∂+∂∂∂∂+∂∂+∂∂∂∂++∂∂∂∂+∂∂-=)][()]([)2(μμλμρ)()())()())())()()()()()][()]([)2(gradu div divu xz w y v x u x gradu div S divu xz w y v x u x S S divu xz w y v x u x gradu div S xw z x v y x u x z u z y u y x u x S xw z u z x v y u y divu x u x Mx u Mx Mx Mx μλμμλμλμμμμμμμμμμλμ+∂∂+∂∂+∂∂+∂∂∂∂=++∂∂+∂∂+∂∂+∂∂∂∂=+∂∂+∂∂+∂∂+∂∂∂∂+=+∂∂∂∂+∂∂∂∂+∂∂∂∂+∂∂∂∂+∂∂∂∂+∂∂∂∂=+∂∂+∂∂∂∂+∂∂+∂∂∂∂++∂∂∂∂((()()( 因为密度为常数,所以连续性方程:0=∂∂+∂∂+∂∂z w y v x u ρρρ 推 得:0=∂∂+∂∂+∂∂z w y v x u 所以:Su= 0)()=∂∂+∂∂+∂∂+∂∂∂∂divu x z w y v x u x λμ(7区域离散为分几种,说明各自的特点。

传热计算习题附详细答案

传热计算习题附详细答案

传热计算题1.在一内径为0.25cm的管轴心位置上,穿一直径为 0.005cm的细导线,用以测定气体的导热系数。

当导线以0.5A 的电流时,产生的电压降为0.12V/cm,测得导线温度为167℃,空心管内壁温度为150℃。

试求充入管内的气体的导热系数试分析仪器精度以外造成结果误差的客观原因。

2.有两个铜质薄球壳,内球壳外径为0。

015m,外球壳内径为 0.1m,两球壳间装入一种其导热系数待测的粉粒料。

内球用电加热,输入功率为 50w,热量稳定地传向外球,然后散发到周围大气中。

两球壁上都装有热电偶,侧得内球壳的平均温度为120℃,外求壳的平均温度为50℃,周围大气环境温度为20℃;设粉粒料与球壁贴合,试求:(1)待测材料的导热系数(2)外球壁对周围大气的传热系数3.有一面积为10cm2带有保护套的热电偶插入一输送空气的长管内,用来测量空气的温度。

已知热电偶的温度读数为300℃,输气管的壁温为 200℃,空气对保护套的对流传热系数为60w/m2.k,该保护套的黑度为 0.8,试估算由于辐射造成的气体温度测量误差。

并叙述减小测量误差的途径。

已知 Stefan-Bohzman常数σ=5.67×10-9w/m2k 。

4.用两个结构尺寸相同的列管换热器按并联方式加热某中料液。

换热器的管束由32根长 3m 的Ф25×3mm 的钢管组成。

壳程为120℃的饱和蒸汽。

料液总流量为20m3/h,按相等流量分配到两个换热器中作湍流流动,由 25℃加热到 80℃。

蒸汽冷凝对流传热系数为8Kw/m2.℃,管壁及污垢热阻可不记,热损失为零,料液比热为 4.1KJ/kg.℃,密度为 1000kg/m3。

试求:(1)管壁对料液的对流传热系数(2)料液总流量不变,将两个换热器串联,料液加热程度有何变化?(3)此时蒸汽用量有无变化?若有变化为原来的多少倍?(两者情况下蒸汽侧对流传热系数和料液物性不变)5.某厂现有两台单壳程单管程的列管式空气加热器,每台传热面积为A0=20m2(管外面积),均由128根Ф25×2.5mm的钢管组成。

数值传热学习题解答(汇总版)

数值传热学习题解答(汇总版)

习题1-7解:由于对称性,取半个通道作为求解区域。

常物性不可压缩流体,二维层流、稳态对流换热的控制方程组为: 质量守恒方程0=∂∂+∂∂yv x u 动量守恒方程 ()()⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂−=∂∂+∂∂22221y u x u x py vu x uu νρ ()()⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂−=∂∂+∂∂22221y v x v y p y vv x uv νρ 能量守恒方程 ()()⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=∂∂+∂∂2222y T xT a y vT x uT 边界条件:进口截面 ()0,,===v c T y u u in ; 平板通道上(下)壁面 0,0=∂∂==yTv u ; 中心线上对称条件: 0,0u T v y y∂∂===∂∂; 出口截面0,0,0=∂∂=∂∂=∂∂xT x v x u ; 或者写:采用数值传热学的处理方法。

图1-10 习题1-7的图示本题如果采用整个通道作为计算区域,应该扣除0.5 分2-3.解:由u x u ∂∂=()xuu ∂∂21=η22y u ∂∂得: 其守恒形式为:()xuu ∂∂=2η22y u ∂∂ 对方程两端在t ∆时间间隔内对其控制容积积分得:()dxdydt x uu t t t nsew ⎰⎰⎰∆+∂∂=⎰⎰⎰∆+∂∂t t t e w n s dydxdt y u 222η()()[]dxdt y u y u dydt uu uu s n t t t ewtt t w e n s ][2⎪⎪⎭⎫ ⎝⎛∂∂−⎪⎪⎭⎫ ⎝⎛∂∂=−⎰⎰⎰⎰∆+∆+η 将()()2)(PE e uu uu uu +=, ()()()2P W w uu uu uu +=,()n PN n y u u y u δ−=⎪⎪⎭⎫ ⎝⎛∂∂,()sSP s y u u y u δ−=⎪⎪⎭⎫ ⎝⎛∂∂。

y y y s n ∆==)()(δδ 带入,得:xdt y u u u ydt uu uu t t t S P N tt tW E ∆∆+−=∆⎥⎦⎤⎢⎣⎡−⎰⎰∆+∆+]2[22)()(η t x yu u u t y uu uu tSt P t N t W t E ∆∆∆+−=∆∆−222)()(η整理得离散方程为:()()0242=∆−+−∆−yu u u xuu uu t P t S t N tWt E η2—3:解:由2221()u 2u u ux x y η∂∂∂===∂∂∂得:原方程的守恒形式为: 222()2u ux yη∂∂=∂∂ 对方程两端在t ∆时间间隔内对其控制容积积分,把可积的部分积出后得:22()t tsne wtu u dtdy +∆−⎰⎰= 2t te wtn s u u dtdx y y η+∆⎡⎤⎛⎫⎛⎫∂∂−⎢⎥ ⎪ ⎪∂∂⎝⎭⎝⎭⎣⎦⎰⎰选定2u 随y 而变化的型线,这里取为阶梯式,即在控制容积内沿y 方向不变,则2222()=y ()t tt ts ne we w ttu u dtdy u u dt +∆+∆−∆−⎰⎰⎰选定2u 随t 而变化的规律,这里采用阶梯式显式,则22()t tewty u u dt +∆∆−⎰= ()()22t t e w u u t y ⎡⎤−∆∆⎢⎥⎣⎦选定uy∂∂随x 而变化的型线,这里取为阶梯式,即在控制容积内沿x 方向不变,则22t tt t e wtt n s n s u u u u dtdx x dt y y y y ηη+∆+∆⎡⎤⎡⎤⎛⎫⎛⎫⎛⎫⎛⎫∂∂∂∂−=∆−⎢⎥⎢⎥ ⎪ ⎪ ⎪ ⎪∂∂∂∂⎝⎭⎝⎭⎝⎭⎝⎭⎣⎦⎣⎦⎰⎰⎰ 选定uy∂∂随t 而变化的规律,这里采用阶梯显式,则 2t ttn s u u x dt y y η+∆⎡⎤⎛⎫⎛⎫∂∂∆−⎢⎥ ⎪ ⎪∂∂⎝⎭⎝⎭⎣⎦⎰= 2t t n s u u t x y y η⎡⎤⎛⎫⎛⎫∂∂−∆∆⎢⎥ ⎪ ⎪∂∂⎢⎥⎝⎭⎝⎭⎣⎦进一步选取u 随x,y 分段线性变化,则2222E Pe u u u += , 222w 2W P u u u +=()nt PtN ty uu y u δ−=⎪⎪⎭⎫ ⎝⎛∂∂n , ()stSt p ts y u u y u δ−=⎪⎪⎭⎫ ⎝⎛∂∂。

传热习题集

传热习题集

一、填空题1、按热量传递的途径不同,我们一般把传热分为 、 、 三种方式。

2、按冷热流体接触的方式不同,我们一般把传热分为 、 、 三种方式。

3、常见的加热介质有 、 、 、 、 等。

4、按传热管的结构形式可分为 、 、 、 换热器等。

5、列管换热器是由 、 、 、 和封头等部分组成。

6、设置列管式换热器折流挡板的目的是 。

7、列管式换热器热补偿的方式有 、 、 。

8、在间壁式换热器中,总传热过程由下列步骤所组成:首先是热流体和管外壁间的_____传热,将热量传给管外壁面;然后,热量由管的外壁面以____________方式传给管的内壁面最后,热量由管的内壁面和冷流体间进行_________________传热。

9、采用饱和蒸汽加热某一种冷流体,若保持冷流体的进口温度T 1/、加热蒸汽压力P 不变,现冷体流量增加,则冷流体出口温度T 2/_______,传热速率Φ______,总传热系数K _______ ,传热平均温度差△T m ________。

10、某圆形管道外有两层厚度相等的保温材料A 和B ,温度分布线如右上图(b )中所示,则λA ______λB (填“﹥”或 “﹤”),将______层材料放在里层时保温效果更好。

11、若间壁侧流体的传热过程α1,α2相差较大(α1<<α2),K 值接近_____________侧的值。

12、金属固体的导热系数是随温度的升高而 ,非金属固体的导热系数是随温度的升高而 。

水和甘油的导热系数是随着温度的升高而 。

13、在导热速率方程式为δλt S Q ∆=,中λ称为: ,单位是 ;14、两流体在列管换热器中并流传热时,冷流体的最高极限出口温度为 ;在逆流传热时,冷流体的最高极限出口温度为 因此,在 流传热时,载热体的用量少 。

15、传热的基本方式可分为 、 、 三种。

间壁式换热器传热过程强化的途径主要有 、 、三种。

16. 在平壁稳定热传导过程中,通过三层厚度相同的平壁,每层间温度变化如图所示。

数值传热学习题答案

数值传热学习题答案

数值传热学习题答案【篇一:习题解答】s=txt>1.1 什么是汇编语言?汇编语言的特点是什么?答:为了克服机器语言难以记忆、表达和阅读的缺点,人们采用具有一定含义的符号作为助忆符,用指令助忆符、符号地址等组成的符号指令称为汇编格式指令(或汇编指令)。

汇编语言是汇编指令集、伪指令集和使用它们规则的统称。

① 55h,70h,70h,65h,72h② 53h,6ch,6fh,77h③ 43h,6fh,6dh,70h,75h,74h,65h,72h ④ 57h,68h,61h,74h 1.7 求下列带符号十进制数的8位基2码补码。

① +127② ?2③ ?128④ +2 答:汇编语言的特点是:(1)执行速度快。

(2)程序短小。

(3)可以直接控制硬件。

(4)可以方便地编译。

(5)辅助计算机工作者掌握计算机体系结构。

(6)程序编制耗时,可读性差。

(7)程序可移植性差。

1.2 把下列十进制数转换成二进制数、八进制数、十六进制数。

①127② 1021 ③ 0.875④ 6.25 答:① 1111111b;177q;7fh ② 1111111101;1775q;3fdh ③0.111 b;0.7q;0.eh ④ 110.01b;6.2q;6.4h1.3 把下列二进制数转换成十进制数。

① 1001.11 ② 101011.10011 ③ 111.011④1011.1答:① 9.75d ② 43.59375d③ 7.375d④ 11.5d1.4 把下列八进制数转换成十进制数。

① 573.06 ② 75.23 ③ 431.7④ 123.45 答:① 379.09375d ② 61.296875d③ 281.875④ 83.5781251.5 把下列十六进制数转换成十进制数。

① 0d5.f4 ② 8ba.7c ③ 0b2e.3a④ 6ec.2d 答:① 213.953125d ② 2234.484375③ 2862.2265625 ④1772.175781251.6 把下列英文单词转换成ascii编码的字符串。

数值传热学习题答案权威版

数值传热学习题答案权威版

数值传热学习题答案权威版习题4-2⼀维稳态导热问题的控制⽅程:022=+??S xTλ依据本题给定条件,对节点2节点3采⽤第三类边界条件具有⼆阶精度的差分格式,最后得到各节点的离散⽅程:节点1: 1001=T 节点2: 1505105321-=+-T T T 节点3:75432=+-T T 求解结果:852=T ,403=T对整个控制容积作能量平衡,有:02150)4020(15)(3=?+-?=?+-=?+x S T T h x S q f f B即:计算区域总体守恒要求满⾜习题4-5在4-2习题中,如果25.03)(10f T T h -?=,则各节点离散⽅程如下:节点1: 1001=T节点2: 1505105321-=+-T T T节点3:25.03325.032)20(4015])20(21[-?+=-?++-T T T T对于节点3中的相关项作局部线性化处理,然后迭代计算;求解结果:818.822=T ,635.353=T (迭代精度为10-4)迭代计算的Matlab 程序如下: x=30; x1=20;while abs(x1-x)>0.0001a=[1 0 0;5 -10 5;0 -1 1+2*(x-20)^(0.25)]; b=[100;-150; 15+40*(x-20)^(0.25)]; t=a^(-1)*b; x1=x; x=t(3,1); endtcal=t习题4-12的Matlab程序%代数⽅程形式A i T i=C i T i+1+B i T i-1+D imdim=10;%计算的节点数A=cos(x);%TDMA的主对⾓元素B=sin(x);%TDMA的下对⾓线元素C=cos(x)+exp(x); %TDMA的上对⾓线元素T=exp(x).*cos(x); %温度数据%由A、B、C构成TDMAcoematrix=eye(mdim,mdim);for n=1:mdimcoematrix(n,n)=A(1,n);if n>=2coematrix(n,n-1)=-1*B(1,n);endif ncoematrix(n,n+1)=-1*C(1,n);endend%计算D⽮量D=(coematrix*T')';%由已知的A、B、C、D⽤TDMA⽅法求解T%消元P(1,1)=C(1,1)/A(1,1);Q(1,1)=D(1,1)/A(1,1);for n=2:mdimP(1,n)=C(1,n)/(A(1,n)-B(1,n)*P(1,n-1));Q(1,n)=(D(1,n)+B(1,n)*Q(1,n-1))/(A(1,n)-B(1,n)*P(1,n-1)); end%回迭Tcal(1,mdim)=Q(1,mdim);for n=(mdim-1):-1:1Tcal(1,n)=P(1,n)*Tcal(1,n+1)+Q(1,n);endTcom=[T;Tcal];%绘图⽐较给定T值和计算T值plot(Tcal,'r*')hold on结果⽐较如下,由⽐较可知两者值⾮常切合(在⼩数点后8位之后才有区别):习题4-14充分发展区的温度控制⽅程如下:)(1rTr r r x T uc p =??λρ对于三种⽆量纲定义w b w T T T T --=Θ、∞∞--=ΘT T T T w 、ww T T T T --=Θ∞进⾏分析如下1)由wb wT T T T --=Θ得:w w b T T T T +Θ-=)(由T 可得:x T x T x T T T x T w b w w b ??Θ-+??Θ=?+Θ-?=??)1(])[(rT r T T r T T T r T w w b w w b ??Θ-+?Θ?-=?+Θ-?=??)1()(])[( 由b T 与r ⽆关、Θ与x ⽆关以及x T ??、rT的表达式可知,除了w T 均匀的情况外,该⽆量纲温度定义在⼀般情况下是不能⽤分离变量法的; 2)由∞∞--=ΘT T T T w 得: ∞∞+Θ-=T T T T w )(由T 可得:xT x T T T x T w w ??Θ=?+Θ-?=??∞∞])[(rx T ??、rT的表达式可知,在常见的四种边界条件中除了轴向及周向均匀热流const q w =的情况外,有0=??rT w,则该⽆量纲温度定义是可以⽤分离变量法的; 3)由wwT T T T --=Θ∞得: w w T T T T +Θ-=∞)(由T 可得:xT x T T T x T w w w ??Θ-=?+Θ-?=??∞)1(])[(rT r T T r T T T r T w w w w ??Θ-+?Θ-=+Θ-?=??∞∞)1()(])[( 同2)分析可知,除了轴向及周向均匀热流const q w =温度定义是可以⽤分离变量法的;习题4-181)采⽤柱坐标分析,写出统⼀的稳态柱坐标形式动量⽅程:S r r r r r r x x w r v r r r u x +++=??+??+??)(1)(1)()(1)(1)(θφλθφλφλφρθφρφρ x 、r 和θ分别是圆柱坐标的3个坐标轴,u 、v 和w 分别是其对应的速度分量,其中x 是管内的流动⽅向;对于管内的层流充分发展有:0=v 、0=w ,0=??xu;并且x ⽅向的源项:x pS ??-=r ⽅向的源项:r pS ??-=r S 1由以上分析可得到圆柱坐标下的动量⽅程: x ⽅向:0)(1)(1=??-+x pu r r r u r r r θλθλ r ⽅向:0=??r pθ⽅向:0=??θp边界条件: R r =,0=u0=r ,0=??r u ;对称线上,0=??θu不考虑液体的轴向导热,并简化分析可以得到充分发展的能量⽅程为:)(1)(1θλθλρ+=??Tr r r T r r r x T uc p 边界条件: R r =,w q r T =??λ;0=r ,0=??rTπθ/0=,0=??-θλT2)定义⽆量纲流速:dxdp R uU 2-=λ并定义⽆量纲半径:R r /=η;将⽆量纲流速和⽆量纲半径代⼊x ⽅向的动量⽅程得:0))1((1))1((122=??-?-+?-xp U dx dp R R R R U dx dp R RR R θληλθηηλληηη上式化简得:ηθηηηηηU U 边界条件:1=η,0=U0=η,0=??ηU;对称线上,0=??θU 定义⽆量纲温度:λ/0R q T T b -=Θ其中,0q 是折算到管壁表⾯上的平均热流密度,即:Rq q wπ=0;由⽆量纲温度定义可得:b T Rq T +Θ=λ将T 表达式和⽆量纲半径η代⼊能量⽅程得:)(1)(100θληλθηηλληηηρ?Θ+?Θ=??R q R R R R q R R R x T uc b p 化简得:)1(1)(10θηθηηηηηρ?Θ+?Θ=??x T u c q R b p (1)由热平衡条件关系可以得:mm m b m p b p p RU U q R u u R q A u u dx dT A u c x T u c x T u c 020221221)(===??=??ππρρρ将上式代⼊式(1)可得:)1(1)(12θηθηηηηη?Θ+?Θ=m U U 边界条件:0=η,0=?Θ?η;1=η,R q q w πη10==?Θ?0=?Θ?θ;πθ=,0=?Θθ单值条件:由定义可知:0/0=-=ΘλR q T T b b b 且: ??Θ=ΘAAb U d AU d A 即得单值性条件:0=Θ??AA UdAUdA3)由阻⼒系数f 及Re 定义有:228)(2/Re ??? ??=-=D D U D u u dx dp D f e m e m me νρ且:m W b m W b m W R q T T D T T q Nu ,0,,0~2)/(2Θ=-=-=λλ5-21.⼀维稳态⽆源项的对流-扩散⽅程如下所⽰:xx u 22??Γ=??φφρ(取常物性)边界条件如下:L L x x φφφφ====,;,00上述⽅程的精确解如下:11)/(00--=--?PeL x Pe L e e φφφφΓ=/uL Pe ρ 2.将L 分成20等份,所以有:=P Pe 201 2 3 4 5 6 ………… …………… 17 18 19 20 21 对于中⼼差分、⼀阶迎风、混合格式和QUICK 格式分别分析如下: 1)中⼼)5.01()5.01(11-?+?++-=i i i P P φφφ 20,2 =i2)⼀阶迎风中间节点: ?-?++++=P P i i i 2)1(11φφφ 20,2 =i3)混合格式当1=?P 时,中间节点:2)5.01()5.01(11-?+?++-=i i i P P φφφ20,2 =i当10,5=?P 时,中间节点: 1-=i i φφ 20,2 =i 4) QUICK 格式*12111)35(8122121---++++++=+--??-??+?i i i i i i i P P P P P φφφφφφφ 2≠i *1111)336(8122121--++++++=+-?-??+?i i i i i i P P P P P φφφφφφ 2=i数值计算结果与精确解的计算程序如下:%except for HS, any other scheme doesnt take Pe<0 into consideration %expression of exact solution y=dsolve('a*b*Dy=c*D2y','y(0)=y0,y(L)=yL','x')y=subs(y,'L*a*b/c','t')y=simple(subs(y,'a*b/c*x','t*X'));ysim=simple(sym(strcat('(',char(y),'-y0)','/(yL-y0)')))y=sym(strcat('(',char(ysim),')*(yL-y0)','+y0'))% in the case of Pe=0y1=dsolve('D2y=0','y(0)=y0,y(L)=yL','x')%grid Pe numbertt=[1 5 10];%dimensionless lengthm=20;%mdim is the number of inner nodemdim=m-1;X=linspace(0,1,m+1);%initial value of variable during calculationy0=1;yL=2;%cal exact solutionfor n=1:size(tt,2)t=m*tt(1,n);if t==0yval1(n,:)=eval(y1);elseyval1(n,:)=eval(y);endend%extra treatment because max number in MATLAB is 10^308 if max(isnan(yval1(:)))yval1=yval1';yval1=yval1(:);indexf=find(isnan(yval1));for n=1:size(indexf,1)if rem(indexf(n,1),size(X,2))==0yval1(indexf(n),1)=yL;elseyval1(indexf(n),1)=y0;endendyval1=reshape(yval1,size(X,2),size(yval1,1)/size(X,2));yval1=yval1';endd=zeros(size(tt,2),mdim);a=repmat([1],size(tt,2),mdim);for n=1:size(tt,2)t=tt(1,n);b(n,:)=repmat([0.5*(1-0.5*t)],1,mdim);c(n,:)=repmat([0.5*(1+0.5*t)],1,mdim);d(n,1)=0.5*(1+0.5*tt(1,n))*y0;d(n,mdim)=0.5*(1-0.5*tt(1,n))*yL;endc(:,1)=0;b(:,mdim)=0;%numerical cal by using TDMA subfuctionyval2=TDMA(a,b,c,d,mdim);yval2=[repmat([1],size(tt,2),1),yval2,repmat([2],size(tt,2),1)]; Fig(1,X,yval1,yval2,tt); title('CD Vs. Exact Solution')% FUS solutiond=zeros(size(tt,2),mdim);a=repmat([1],size(tt,2),mdim);for n=1:size(tt,2)t=tt(1,n);b(n,:)=repmat([1/(2+t)],1,mdim);c(n,:)=repmat([(1+t)/(2+t)],1,mdim);d(n,1)=(1+tt(1,n))/(2+tt(1,n))*y0;d(n,mdim)=1/(2+tt(1,n))*yL;endc(:,1)=0;b(:,mdim)=0;%numerical cal by using TDMA subfuctionyval3=TDMA(a,b,c,d,mdim);yval3=[repmat([1],size(tt,2),1),yval3,repmat([2],size(tt,2),1)]; Fig(2,X,yval1,yval3,tt); title('FUS Vs. Exact Solution')% HS solutiond=zeros(size(tt,2),mdim);a=repmat([1],size(tt,2),mdim);t=tt(1,n);if t>2b(n,:)=repmat([0],1,mdim);c(n,:)=repmat([1],1,mdim);d(n,1)=y0;elseif t<-2b(n,:)=repmat([1],1,mdim);c(n,:)=repmat([0],1,mdim);d(n,mdim)=yL;elseb(n,:)=repmat([0.5*(1-0.5*t)],1,mdim);c(n,:)=repmat([0.5*(1+0.5*t)],1,mdim);d(n,1)=0.5*(1+0.5*t)*y0;d(n,mdim)=0.5*(1-0.5*t)*yL;endendc(:,1)=0;b(:,mdim)=0;% numerical cal by using TDMA subfuctionyval4=TDMA(a,b,c,d,mdim);yval4=[repmat([1],size(tt,2),1),yval4,repmat([2],size(tt,2),1)]; Fig(3,X,yval1,yval4,tt); title('HS Vs. Exact Solution')%QUICK Solutiond=zeros(size(tt,2),mdim);a=repmat([1],size(tt,2),mdim);for n=1:size(tt,2)t=tt(1,n);b(n,:)=repmat([1/(2+t)],1,mdim);c(n,:)=repmat([(1+t)/(2+t)],1,mdim);d(n,1)=(1+tt(1,n))/(2+tt(1,n))*y0;d(n,mdim)=1/(2+tt(1,n))*yL;endc(:,1)=0;b(:,mdim)=0;%numerical cal by using TDMA subfuctionyval5=zeros(size(tt,2),mdim);yval5com=yval5+1;counter=1;%iterativewhile max(max(abs(yval5-yval5com)))>10^-10if counter==1yval5com=TDMA(a,b,c,d,mdim);endfor nn=1:size(tt,2)for nnn=1:mdimif nnn==1d(nn,nnn)=((6*yval5com(nn,nnn)-3*y0-3*yval5com(nn,nnn+1))*tt(1,nn))/(8*(2+tt(1, nn)))+((1+tt(1,nn))/(2+tt(1,nn))*y0);elseif nnn==2d(nn,nnn)=((5*yval5com(nn,nnn)-3*yval5com(nn,nnn+1)-yval5com(nn,nnn-1)-y0)*tt (1,nn))/(8*(2+tt(1,nn)));elseif nnn==mdimd(nn,nnn)=((5*yval5com(nn,nnn)-3*yL-yval5com(nn,nnn-1)-yval5com(nn,nnn-2))*tt (1,nn))/(8*(2+tt(1,nn)))+(1/(2+tt(1,nn))*yL); elsed(nn,nnn)=((5*yval5com(nn,nnn)-3*yval5com(nn,nnn+1)-yval5com(nn,nnn-1)-yval5 com(nn,nnn-2))*tt(1,nn))/(8*(2+tt(1,nn))); endendendyval5=TDMA(a,b,c,d,mdim);temp=yval5;yval5=yval5com;yval5com=temp;counter=counter+1;endyval5=yval5com;yval5=[repmat([1],size(tt,2),1),yval5,repmat([2],size(tt,2),1)];Fig(4,X,yval1,yval5,tt);title('QUICK Vs. Exact Solution')%-------------TDMA SubFunction------------------function y=TDMA(a,b,c,d,mdim)%form a b c d resolve yval2 by using TDMA%eliminationp(:,1)=b(:,1)./a(:,1);q(:,1)=d(:,1)./a(:,1);for n=2:mdimp(:,n)=b(:,n)./(a(:,n)-c(:,n).*p(:,n-1));q(:,n)=(d(:,n)+c(:,n).*q(:,n-1))./(a(:,n)-c(:,n).*p(:,n-1));end%iterativey(:,mdim)=q(:,mdim);for n=(mdim-1):-1:1y(:,n)=p(:,n).*y(:,n+1)+q(:,n);end%-------------ResultCom SubFunction------------------ function y=ResultCom (a,b,c)for n=1:max(size(c,2))y(2*n-1,:)=a(n,:);y(2*n,:)=b(n,:);end%-------------Fig SubFunction------------------ function y=Fig(n,a,b,c,d)figure(n);plot(a,b);hold onplot(a,c,'*');str='''legend(';for n=1:size(d,2)if n==size(d,2)str=strcat(str,'''''Pe=',num2str(d(1,n)),''''')''');elsestr=strcat(str,'''''Pe=',num2str(d(1,n)),''''',');endendeval(eval(str));精确解与数值解的对⽐图,其中边界条件给定10=φ,2=L φ。

传热习题集

传热习题集

本章总结
Biot准数 定义: Bi= hδ /λ= (δ /λ)/(1/h) δ 特征尺度 厚度、半径 物理意义: 内部导热热阻与表面对流热阻之比。 表征换热过程中各点温度趋于一致的能力
平壁一维稳态导热
t1
t2
q
模型
∂ 2t ∂ 2t ∂ 2t Φ ∂t = a + 2+ 2 + 2 ∂z ρc ∂τ ∂y ∂x
δ
x
d t =0 2 dx
2
边界条件:
温度分布
x=0
t= t 2 − t1
t = t1 ;
x + t1
x =δ
t = t2
t1 − t n +1
0.9619
qL = t1 − t 3 1 2π λ 1 = 1 ( ( 1 ln d2 d1 ln + 1
λ2
ln
d3 d2
)
57874.2 3.54 1200 − ln 753 ℃ = = 2π × 3.8296 2.94
1200 − 150 1 3.54 + 1 0.9619 ln 3.77 3.54 )
δ
多层平壁
dt t1 − t 2 =λ 热流量计算 q = −λ dx δ
q=
∑δ
i =1
n
i
/ λi
dt h1 (tf1 − t1 ) = −λ dx dt ( ) h2 t 2 − tf2 = −λ dx
d 2t =0 2 dx
x =0
x =δ
dt ∆t q = −λ = dx 1 + δ + 1 h1 λ h2
解: 由 Bi = h ⋅ r = 116 × 0.050 = 0.036 ≤ 0.1× 1

传热习题课计算题

传热习题课计算题

1、现测定一传热面积为2m2的列管式换热器的总传热系数K值。

已知热水走管程,测得其流量为1500kg/h,进口温度为80℃,出口温度为50℃;冷水走壳程,测得进口温度为15℃,出口温度为30℃,逆流流动。

(取水的比热c p=4.18×103J/kg·K)解:换热器的传热量:Q =q m c p (T 2-T 1)=1500/3600×4.18×103×(80-50)=52.25kW传热温度差△t m :热流体 80 → 50 冷流体 30 ← 15△t 1=50, △t 2=352355021<=∆∆t t 传热温度差△t m 可用算数平均值:5.4223550221=+=∆+∆=∆t t t m ℃⋅=⨯⨯=∆=23/6155.4221025.52m W t A Q K m ℃2、一列管换热器,由φ25×2mm 的126根不锈钢管组成。

平均比热为4187J/kg·℃的某溶液在管内作湍流流动,其流量为15000kg/h ,并由20℃加热到80℃,温度为110℃的饱和水蒸汽在壳方冷凝。

已知单管程时管壁对溶液的传热系数αi 为520W/m 2·℃,蒸汽对管壁的传热系数α0为1.16×104W/m 2·℃,不锈钢管的导热系数λ=17W/m·℃,忽略垢层热阻和热损失。

试求:管程为单程时的列管长度(有效长度)(总传热系数:以管平均面积为基准,00111d d b d d K m i m i ⋅++⋅=αλα)解:传热量:Q =q m c p (t 2-t 1) =15000/3600×4187×(80-20) ≈ 1.05×106W总传热系数:(以管平均面积为基准)1111152023210002171116102325004Kd d b d d K i m i m =⋅++⋅=⋅++⨯⋅αλα .. 解得: K =434.19W/m 2·℃ 对数平均温差: 1101102080△t 1=90 △t 2=30∆∆∆∆∆t t t t t m =-=-=12129030905461ln ln .℃ 传热面积: Q KA t m m =∆ A Q K t m m m==⨯⨯=∆10510434195461442862....A n d L m m =π; 列管长度:L A n d m m m ==⨯⨯≈π44281263140023487....3、有一列管式换热器,装有φ25×2.5mm钢管320根,其管长为2m,要求将质量流量为8000kg/h的常压空气于管程由20℃加热到85℃,选用108℃饱和蒸汽于壳程冷凝加热之。

传热习题

传热习题

传热习题一. 填空与选择1.①一列管换热器,列管规格为φ38×3,管长4m, 管数127根,则外表面积F1=_________________,而以内表面积计的传热面积F2=_____________________。

②两流体的间壁换热过程中 Q= m Cp△t计算式中, △t 表示_________________。

答:①F1=127×4π×0.038 ,F2=127×4π×0.032;②T1-T2 或t2– t1 。

2.①稳定热传导是指______________________________不随时间而改变。

②套管换热器,内管外径d1、外管内径d 2,则环隙的当量直径de=______________。

答:①传热系统中各点的温度仅随位置变;②d2– d1 。

3.①物体黑度是指在__________温度下,物体的__________与__________之比,在数值上它与同一温度下物体的__________相等。

②两固体间辐射传热速率公式为____________________________________________。

t04a05030答:①同一,辐射能力,黑体辐射能力,吸收率;②Q1——2=C1——2φA[(T1/100)4 - (T2/100)4 ] 。

4.①将单程列管式换热器改为双程的作用是__________,但这将使________减小,______ 增大。

②对流过程是__________和__________之间的传热过程。

答:①提高管内α,△t m,阻力;②流体,壁面5.①为了减少室外设备的热损失,保温层外所包的一层金属皮应该是________A.表面光滑,颜色较浅B.表面粗糙,颜色较深C.表面粗糙,颜色较浅②某一套管换热器用管间饱和蒸汽加热管内空气,设饱和蒸汽温度为100℃ , 空气进口温度为20℃,出口温度为80℃,问此套管换热器内管壁温应是___________ A.接近空气平均温度B.接近饱和蒸汽和空气的平均温度C.接近饱和蒸汽温度答:① A;② C 。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

简答题集锦
1.流动与传热数值模拟的基本任务是什么?
(把原来在时间域及空间域上连续的物理量的场,如速度场和压力场,用一系列有限个离散点上的变量值的集合来代替,通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组,然后求解代数方程组获得场变量的近似值CFD可以看做是在流动基本方程(质量守恒方程飞动量守恒方程、能量守恒方程)控制下对流动的数值模拟。

通过这种数值模拟,我们可以得到极其复杂问题的流场内各个位置上的基本物理量(如速度、压力、温度、浓度等)的分布,以及这些物理量随时间的变化情况,确定旋涡分布特性、空化特性及脱流区等。


2.数值模拟过程如何实现,主要步骤是那些?
(建模、网格划分、坐标系、数学方程、求解、后处理)
a.建立反映工程问题或物理过程本质的数学模型;
b.选择与计算区域的边界相适应的坐标系;
c.建立网格;
d.建立离散方程;
e.求解代数方程组;
f.后处理,显示计算结果
3.建立离散方程有哪些主要方法?比较说明各种方法的优缺点?(有限差分、有限体积、有限元、有限分析等)
4什么叫控制方程?常见的控制方程有哪几个?各用在什么场合?
5试写出控制方程的通用形式,并说明通用形式中各项的意义?(写明通式,以及各个方程中通式的表达形式)
6推导x 方向的动量控制方程中的源项u S 的表达式。

由此证明当密度
和黏度为常数时,u S 变为0。

X 方向N-S 方程:
Mx S x
w z u z x v y u y divu x u x x p Dt Du +∂∂+∂∂∂∂+∂∂+∂∂∂∂++∂∂∂∂+∂∂-=)][()]([)2(μμλμρ)()())()())())()()()()()][()]([)2(gradu div divu x
z w y v x u x gradu div S divu x
z w y v x u x S S divu x
z w y v x u x gradu div S x
w z x v y x u x z u z y u y x u x S x
w z u z x v y u y divu x u x Mx u Mx Mx Mx μλμμλμλμμμμμμμμμμλμ+∂∂+∂∂+∂∂+∂∂∂∂=++∂∂+∂∂+∂∂+∂∂∂∂=+∂∂+∂∂+∂∂+∂∂∂∂+=+∂∂∂∂+∂∂∂∂+∂∂∂∂+∂∂∂∂+∂∂∂∂+∂∂∂∂=+∂∂+∂∂∂∂+∂∂+∂∂∂∂++∂∂∂∂((()()( 因为密度为常数,所以连续性方程:
0=∂∂+∂∂+∂∂z w y v x u ρρρ 推 得:0=∂∂+∂∂+∂∂z w y v x u 所以:Su= 0)()=∂∂+∂∂+∂∂+∂∂∂∂divu x z w y v x u x λμ(
7区域离散为分几种,说明各自的特点。

(内节点法、外节点法)
先节点后界面
*在边界处控制0个容积,源项或不稳态项表现得较好
8.什么叫控制方程离散化?意义是什么?
答:控制方程的离散化是指:在将求解区域进行离散的前提下,建立关于各离散点上变量的值之间的制约关系,将连续的偏微分方程组及其定解条件按照某种方法比如Taylor级数展开法或是控制容积法等方法遵循特定的规则在计算区域的离散网格上转化为代数方程组,以得到连续系统的离散数值逼近解。

意义:相比于用解析法求得控制方程的解,将控制方程离散化更容易分析物理现象的具体规律,求得离散函数的解。

9常用的控制方程离散化方法有哪些?各有何特点?
常用的控制方程离散化方法有:
Taylor级数法、控制容积法以及多项式拟合法1.Taylor级数法:借助Taylor级数展开给出各阶导数的差商
表达式,将方程中的各阶层数用相应的差商表达式代替导数。

特点:a.着眼于节点上的微分平衡,不需要分布函数;b.在均匀网格,采用外节点法以及源项采用线性分布的情况下,用Taylor级
数法得到较为精确的结果;C.主要的差分格式几乎都可以用
Taylor级数展开法得到。

2.多项式拟合法:其思想是假设计算区域中解变量呈局部的一
维分布,用一个多项式拟合控制方程,然后求导代到微分方程。

3.控制容积法:将控制方程在控制容积上积分从而得到离散化方程的
离散化方法。

其特点:a.需要假定待求变量的分布函数; b.着眼于控制容积上流的平衡,在尺寸的选择上较为灵活; c.采用均匀网格至关重要,采用内节点法划分网格无法实现均匀网格划分,在采用均匀网格的前提条件下,控制容积法能得到更为精确的结果。

d.控制容积法简捷、物理概念清楚、以及得到广泛应用.
10简述有限体积法的基本思想。

(这种方法从描述流动与传热问题的守恒型控制方程出发,对它在控制容积上积分,在积分过程中需要对界面上被求函数本身(对流通量)及其一阶导数(扩散通量)的构成方式作出假设,从而形成不同的离散方程。


11简述有限差分法的基本思想。

(有限差分法的第一步是离散求解域,其次对每个节点所描述的流动与传热问题的偏微分方程中的导数项用相应的差分表达式来代替,从而在每个节点上形成一
个代数方程,其中包含了本节点及其邻点上所求量的未知值。


12简述瞬态问题与稳态问题之控制方程的区别,说明在时间域上离散控制方程的基本思想及方法,对比显式及全隐式时间积分方案的异同,给出这两种方案下所生成的瞬态对流-扩散问题的离散方程。

答:1)瞬态的控制方程比稳态的控制方程多了一个时间导数项(瞬态项)。

2)显式及全隐式时间积分方案的异同:
显式:按每一层的初始时刻之值来计算,所形成的离散方程
隐式:按每一层的终了时刻之值来计算。

3)显式:对流项
扩散项:
全隐式:
13为什么计算流体力学只能得到真实流场的近似数值解?
主要因素有三点:数学模型的简化、离散格式的近似、迭代过程的近似
14有哪些生成二阶导数的差分格式,简述其方法。

答:Taylor级数法:借助Taylor级数展开给出各阶导数的差商表达式,将方程中的各阶层数用相应的差商表达式代替导数。

多项式拟合法:由一阶导数的差分格式生成二阶导数,其思想是假设计算区域中解变量呈局部的一维分布,用一个多项式
拟合控制方程,然后求导代到微分方程。

15有限体积法在离散积分形式方程时有那些步骤?
(划分控制体,用数值积分方法将积分方程离散为代数方程,将控制体表面的物理量近似成控制点上物理量的插值形式)
16采用Simple方法求解NS方程时采用的压力校正方法有那些步骤,为什么SIMPLE方法要采用交错网格?
(假定压力场,求解动量方程得到速度场近似解,利用动量校正方程和连续性方程得到压力校正方程,求解后得到压力场的校正量校正压力场和速度场,利用新的压力场和速度场重新进行迭代,直至收敛。

采用交错网格是为了避免出现计算结果的压力交错现象)
17利用Talyer级数法生成等距网格上的一阶偏导数的三阶精度中心差分格式(提示:以i点为中心节点,利用i-1,i+1,i+2点的展开式)答案:略
补充:
18.试写出圆管内流体流动的控制方程及其边界条件(假定没有热交换),并写出CFD分析时的求解步骤
答案在讲义:《无限大平壁的瞬态非稳态导热》第81—84页
19.假设S 为常数,k 与x 有关,在均匀网格计算域中,w e x x x )()(δδ-=∆,推导
方程022=++S dr
dT r k dr T d k 的离散方程,假定: 222)()2(x T T T k dr T d k P W E p ∆-+=以及x
T T dr dT W E ∆-=2)( 把dk/dr 看成一给定的量,并分析在什么情况下系数a E 和a W 会变成负值,以及a E 和a W 会变成负值后对求解分布有什么影响?
20.求解离散方程过程中产生误差的原因?
答:模型误差、截断误差、离散误差以及传入误差,最主要是由离散误差所造成的,因划分的网格数目及步长有限,导致只能用有限个控制容积进行离散,因此产生误差。

21.什么是离散方程的相容性、收敛性和稳定性?。

相关文档
最新文档