传热学 热对流 计算 (1)
1、水以1.5m /s 的速度流过内径为25mm的加热管。管的内壁温度保持100℃，水的进口温度为15℃。若要使水的出口温度达到85℃，求单位管长换热量（不考虑修正）。已知50℃的水λf =0.648 W/(m.K)，νf =0.566×10-6m2/s，Pr =3.54。
2、取外掠平板边界层的流动由层流转化为湍流的临界雷诺数5×105,试计算25℃的空气和水达到临界雷诺数时所需要的平板长度,取u =1m/s，ν空气=15.53×10-6m2/s,ν水=0.905×10-6。
3、试推导努谢尔特关于层流膜状凝结的理论解
4、用实验测定一薄壁管流体平均对流换热系数。蒸汽在管外凝结并维持管内壁温度为100℃。水在管内流动流量为G=0.5Kg/s，水温从15℃升到45℃。管的内径d=50mm，长L=4.5m。试求管内流体与壁面间的平均换热系数。已知水在30℃时c p=4.174KJ/(Kg.K)
5、以0.8m/s 的流速在内径为2.5cm 的直管内流动，管子内表面温度为60℃，水的平均温度为30℃，管长2m ，试求水所吸收的热量。已知30℃时水的物性参数为：Pr =5.42，c p =4.17KJ/(Kg.K)，λ=61.8×10-2
W/(m.K)，ρ=995.7Kg/m 3，μ =80.15×10-6 Kg/(m.s)；水60℃ 时的ν=0.4699×10-6 m 2/s ，水在管内流动准则方程式为
4
.08.0Pr Re 027.0f f f Nu =，适用条件：Re f =104-1.2×105，Pr f
=0.6-120，水与壁面间的换热温差Δt ≤30℃。
6、计算一空气横掠管束换热的空气预热器的对流换热量。已知管束有25排，每排12根光管，管外径25mm ，管长 1.5m ，叉排形式，横向管间距S 1=50mm ，纵向管间距S 2=38mm ，管壁温度120℃，空气来流速度u f =4m/s ，空气进口温度20℃，出口温度40℃。已知空气物性：λf =0.0267W/(m.K)，νf =16.0×10-6m 2/s ，Pr f =0.701。最大流速u max = u f
S 1/(S 1-d)；推荐关联式：m
w
f
f n f f c Nu ???
?
??=Pr Pr Pr Re 36.0（公式适
用条件：N ≥20，光管管束，Pr f =0.7～500，除Pr w 的定性温度为壁温外，其余定性温度为流体在管束中的平均温度。指数m 对气体m =0，对液体m =0.25，
c、n值查下表。（10分）
计算传热学
1、已知：一块厚度为0.1mm 的无限大平板，具有均匀内热源，q ＝50×103W/m 3，，导热系数K ＝10W/m.℃，一侧边界给定温度为75℃，另一侧对流换热，T f ＝25℃，，h=50W/m 2.℃,求解稳态分布。（边界条件用差分代替微分和能量平衡法），画图。(内,外节点) 2、试以下述一维非稳态导热问题为模型，编写求解一维非稳态扩散型问题的通用程序： 00 00000()()()() L L f x x x x L fL L x x x x T T k s c x x T k h T T W x T k h T T W x T T x τρτ =====???+=????=-+??-=-+?= 其中，x 是空间坐标变量，τ是时间坐标变量，T 是温度（分布），k 是材料的导热系数，s 是内热源强度，ρ是材料的密度，c 是材料的比热，h 0和h L 分别是x 0和x L 处流体与固体壁面间的换热系数，而T f0和T fL 分别是固体壁两侧流体的温度，W 0和W L 是x 0和x L 处(非对流换热)热流密度，T 0(x )是固体壁内初始温度分布。注意k 、ρ、c 、s 、h 0 、h L 、W 0和W L 均可以是温度T 和/或空间坐标x 的函数。 具体要求： 1) 将数学模型无量纲化； 2) 考虑各种可能的边界条件和初始条件组合 3) 提供完整的程序设计说明，包括数学推导过程和程序使用说明 3、对于有源项的一维稳态方程， s dx d T dx d u dx d +=)()(φφρ 已知 x=0,φ=0,x=1, φ=1.源项S=0.5-X 利用迎风格式、混合格式、乘方格式求解φ的分布.
传热学数值计算大作业2014011673
数值计算大作业 一、用数值方法求解尺度为100mm×100mm 的二维矩形物体的稳态导热问题。物体的导热系数λ为1.0w/m·K。边界条件分别为： 1、上壁恒热流q=1000w/m2; 2、下壁温度t1=100℃； 3、右侧壁温度t2=0℃; 4、左侧壁与流体对流换热，流体温度tf=0℃，表面传热系数 h 分别为1w/m2·K、10 w/m2·K、100w/m2·K 和1000 w/m2·K; 要求： 1、写出问题的数学描述； 2、写出内部节点和边界节点的差分方程； 3、给出求解方法； 4、编写计算程序(自选程序语言)； 5、画出4个工况下的温度分布图及左、右、下三个边界的热流密度分布图； 6、就一个工况下（自选）对不同网格数下的计算结果进行讨论； 7、就一个工况下（自选）分别采用高斯迭代、高斯——赛德尔迭代及松弛法（亚松弛和超松弛）求解的收敛性（cpu 时间，迭代次数）进行讨论； 8、对4个不同表面传热系数的计算结果进行分析和讨论。 9、自选一种商业软件（fluent 、ansys 等）对问题进行分析，并与自己编程计算结果进行比较验证（一个工况）。（自选项） 1、写出问题的数学描述 设H=0.1m 微分方程 22220t t x y ??+=?? x=0，0 y=H ，0 取步长δx=0.02。已知x=0,Φ=0;x=1,Φ=1.令k=ρｕ/Γ计算结果图表： 程序及数据结果： 追赶法： #include a[i]=2+0.02*k; b[i]=4; c[i]=2-0.02*k; f[i]=0; } tdma(a,b,c,f,x); for(i=0;i 4-1 解：采用区域离散方法A 时；网格划分如右图。内点采用中心差分 23278.87769.9 T T T === 22d T T=0dx - 有 i+1i 12 2+T 0i i T T T x ---=? 将2点,3点带入 32122 2+T 0T T T x --=? 即321 209T T -+= 432322+T 0T T T x --=?4321322+T 0T T T x --=? 即4 321 209 T T T -+-= 边界点4 （1）一阶截差 由x=1 1dT dx =，得 431 3 T T -= （2）二阶截差 11B M M q x x x T T S δδλλ -=++ 所以 434111. 1. 36311 T T T =++ 即 431 22293 T T -= 采用区域离散方法B 22d T T=0dx - 由控制容积法 0w e dT dT T x dT dT ????--?= ? ????? 所以代入2点4点有 322121011336 T T T T T ----= 即 239 028T T -= 544431011363 T T T T T ----= 即 34599 02828T T T -+= 对3点采用中心差分有 432 32 2+T 013T T T --=?? ??? 即 23499 01919 T T T -+= 对于点5 由x=1 1dT dx =，得 541 6 T T -= （1）精确解求左端点的热流密度 由 ()2 1 x x e T e e e -= -+ 所以有 ()22 20.64806911x x x x dT e e q e e dx e e λ -====- +=-=++ （2）由A 的一阶截差公式 21 0.247730.743113 x T T dT q dx λ =-=-= =?= （3）由B 的一阶截差公式 0 0.21640 0.649213 x dT q dx λ =-=-= = （4）由区域离散方法B 中的一阶截差公式： 210.108460.6504()B B T T dT dx x δ-?? ==?= ? ?? 通过对上述计算结果进行比较可得：区域离散B 有控制容积平衡法建立的离散方程与区域离散方程A 中具有二阶精度的格式精确度相当！ 4-3 解：将平板沿厚度方向3等分，如图 传热学练习题 一、填空题 1、在范德瓦耳斯方程中， 是考虑分子之间的斥力而引进的改正项，V an 2 2 是考虑到分子之间的 而引进的改正项。 2、在等压过程中，引进一个函数H 名为焓则其定义为 ，在此过程中焓的变化为 ，这正是等压过程中系统从外界吸收的热量。 3、所在工作于一定温度之间的热机，以 的效率为最高，这是著名的 。 4、一个系统的初态A 和终态B 给定后，积分 与可逆过程的路径无关，克劳修斯根据这个性质引进一个态函数熵，它的定义是 ，其中A 和B 是系统的两个平衡态。 5、在热力学中引入了一个态函数TS U F -=有时把TS 叫做 ，由于F 是一个常用的函数，需要一个名词，可以把它叫做 。 6、锅炉按用途可分为电站锅炉、___________ 锅炉和生活锅炉。 7、锅炉按输出介质可分为、___________ 、__________ 和汽水两用锅炉。 8、锅炉水循环可分为___________ 循环和_________ 循环两类。 9、如果温度场随时间变化，则为__________。 10、一般来说，紊流时的对流换热强度要比层流时__________。 11、导热微分方程式的主要作用是确实__________。 12、一般来说，顺排管束的平均对流换热系数要比叉排时__________。 13、膜状凝结时对流换热系数__________珠状凝结。 二、判断题 1、系统的各宏观性质在长时间内不发生任何变化，这样的状态称为热力学平衡态。 （ ） 2、温度是表征物体的冷热程度的，温度的引入和测量都是以热力学定律为基础的。 （ ） 3、所谓第一类永动机，就是不需要能量而永远运动的机器。 （ ） 4、自然界中不可逆过程是相互关联的，我们可以通过某种方法把两个不可逆过程联系起来。 （ ） 5、对于处在非平衡的系统，可以根据熵的广延性质将整个系统的熵定义为处在局域平衡的各部分的熵之和。（ ） 6、 测量锅炉压力有两种标准方法，一种是绝对压力，一种是相对压力都称为表压力。（ ） 沈阳航空航天大学 预测燃气涡轮燃烧室出口温度场 沈阳航空航天大学 2013年6月28日 计算传热学 图1模型结构和尺寸图 1.传热过程简述 计算任务是用计算流体力学/计算传热学软件Fluent求解通有烟气的法兰弯管包括管内烟气流体和管壁固体在内的温度分布，其中管壁分别采用薄壁和实体壁两种方法处理。在进行分析时要同时考虑导热、对流、辐射三种传热方式。 (1) 直角弯管内外壁面间的热传导。注意：如果壁面按薄壁处理时，则不用考虑此项，因为此时管壁厚度忽略不计，内壁和外壁温度相差几乎为零。 (2) 管道外壁面与外界环境发生的自然对流换热。由于流体浮生力与粘性力对自然对流的影响，横管与竖管对流换热系数略有不同的。计算公式也不一样。同时，管道内壁面同烟气发生的强制对流换热。 (3) 管道外壁和大空间（环境）发生辐射换热 通过烟气温度和流量，我们可以推断出管道内烟气为湍流流动。这在随后的模 沈阳航空航天大学 拟计算中可以得到证实。 2.计算方案分析 2.1 控制方程及简化 2.1.1质量守恒方程： 任何流动问题都要满足质量守恒方程，即连续方程。其积分形式为： 0vol A dxdydz dA t ρρ?+=?????? 式中，vol 表示控制体；A 表示控制面。第一项表示控制体内部质量的增量，第二项表示通 过控制面的净通量。 直角坐标系中的微分形式如下： ()()()0u v w t x y z ρρρρ????+++=???? 上式表示单位时间内流体微元体中质量的增加，等于同一时间段内流入该微元体的净增量。 对于定常不可压缩流动，密度ρ为常数，该方程可简化为 0u v w x y z ???++=??? 2.1.2动量守恒方程： 动量守恒方程也是任何流动系数都必须满足的基本定律。数学式表示为： F m dv dt δδ= 流体的粘性本构方程得到直角坐标系下的动量守恒方程，即N-S 方程： ()()()u u p div Uu div gradu S t x ρρμ??+=+-?? ()()()v v p div Uv div gradv S t y ρρμ??+=+-?? ()()()w w p div Uw div gradw S t z ρρμ??+=+-?? 该方程是依据微元体中的流体的动量对时间的变化率等于外界作用在该微元体上的各种力之和。式中u S 、v S 、w S 是动量方程中的广义源项。和前面方程一样上式 1、用简捷方法确定附图中的角系数X 12。 2、一直径为4cm 的小铜球，初始温度为500℃，突然放置于10℃的空气中，假设铜球表面与周围环境的对流换热系数为30W/(m 2.K)，试计算铜球冷却到200℃所需要的时间。已知铜球的比热c p =0.377KJ/(Kg.K)， ρ=8440Kg/m 3，λ=109W/(m.K)。 3、水以1.5m /s 的速度流过内径为25mm 的加热管。管的内壁温度保持 100℃，水的进口温度为15℃。若要使水的出口温度达到85℃，求单位管长换热量（不考虑修正）。已知50℃的水λf =0.648 W/(m.K)，νf =0.566×10-6m 2/s ，Pr =3.54 解答 1、（1）11,21,222,11==X X A X A ，则5.02/422 1,2122 ,1===R R X A A X ππ （2）同上125.02 /44/221,2122,1===R R X A A X ππ 2、首先检验是否可以采用集总参数法。 () M A V h Bi v 1.000183.0109 3/02.030/<=?==λ（其中M=1/3） 可以采用集总参数法， )/1(1014.143/02.0437*******.04304 33 s V c hA P -?=?????=ππρ () τ??-=--=---∞∞401014.14exp 1050010200t t t t 可得：hour s 186.0670==τ 3、定性温度502 85152=+=+=out in f t t t （℃） 446101074.610 556.0025.05.1Re >?=??==-f f v ud 流动为紊流。 3.278Pr Re 023.0 4.08.0 ==f f Nu 则单位管长换热量 ()()28311501007213025.014.3=-???=-=f w l t t dh q π（W/m ） 《传热学》复习题一 一、名词解释（3分×5=15分） 1、导热： 2、对流换热： 3、肋片效率： 4﹑膜状凝结： 5﹑灰体； 二、填空题（2分×5=10分） 1.空间辐射热阻可表示为：（）。 2、（）是热量传递的动力。 3、二维、常物性、无内热源、直角坐标系中的稳态导热微分方程式为（）。 4、角系数的确定方法有（）两种。 5、温度场是指）。 三、判断题（2分×5=10分） 1、一个灰表面如果是漫反射的，则一定是漫辐射的。（） 2、热量传递的三种基本方式是：导热，对流换热，辐射换热。（） 3﹑当一个表面的吸收率α=1时，可当作一个绝热表面来处理。（） 4、管内强迫对流换热时，假定条件相同，弯管的比直管的换热系数大。（） 5、一个表面的有效辐射一定不大于它的本身辐射。（）。 四、简答题（6分×4=24分） 1、说明哪些因素影响了对流换热 2、简述热辐射的三个特点. 3、说明大气层的温室效应 4、写出努谢尔特准则Nu的表达式并说明努谢尔特准则Nu的物理含义 五、计算题（41分） 1、（14分）有一气体冷凝器，气侧对流换热表面传热系数h1=95W/(㎡.k),壁厚为㎜,λ=。水侧对流换热表面传热系数h2=5800W/(㎡.k)。计算每个环节的热阻以及传热热阻 2、（13分）加热炉置于25℃的厂房内，加热炉外形尺寸为高2.5m、宽3.5m、 长4m，加热外表面温度均匀且维持55℃。如果不考虑辐射作用，试计算加热 炉炉墙的散热量。已知Νu=C n,C=,n=1/3,空气的物性参数为：λ=,Pr=,粘度ν =㎡/S. 3、（14分）温度为99℃的热水进入一个逆流式换热器，并将4℃的冷 水加热到32℃，冷水的流量为1.3Kg/S，热水的流量为2.6Kg/S，总传热系 数为830W/（㎡.K ）。试计算换热器面积为多少 （水的比热为C P= ） 、室内一根水平放置的无限长的蒸汽管道，其保温层外径d=583 mm，外表面 实测平均温度及空气温度分别为，此时空气与管道外 表面间的自然对流换热的表面传热系数h=3.42 W /(m2 K),墙壁的温度近似取为 室内空气的温度，保温层外表面的发射率 问：(1)此管道外壁的换热必须考虑哪些热量传递方式； (2)计算每米长度管道外壁的总散热量。(12分) 解： (1)此管道外壁的换热有辐射换热和自然对流换热两种方式。 (2)把管道每米长度上的散热量记为qi 当仅考虑自然对流时，单位长度上的自然对流散热 q i,c =二d h t =二dh (j - t f ) = 3.14 0.583 3.42 (48 - 23 ) 二156 .5(W / m) 近似地取墙壁的温度为室内空气温度，于是每米长度管道外表面与室内物体及墙壁 之间的辐射为： q i厂d (T； -T；) = 3.14 0.583 5.67 10》0.9 [(48 273)4-(23 273)4] = 274.7(W /m) 总的散热量为q i = q i,c +q i,r = 156.5 +274.7 = 431.2(W/m) 2、如图所示的墙壁，其导热系数为50W/(m- K),厚度为50mm在稳态情况下的 墙壁内的一维温度分布为：t=200-2000x 2，式中t的单位为°C, x单位为m 试 求： t (1) 墙壁两侧表面的热流密度； (2) 墙壁内单位体积的内热源生成的热量 2 t =200 —2000x 解：(1)由傅立叶定律: ① dt W q ' (―4000x) = 4000二x A dx 所以墙壁两侧的热流密度： q x _. =4000 50 0.05 =10000 (1)由导热微分方程 茫?生=0得: dx 扎 3、一根直径为1mm 勺铜导线，每米的电阻为2.22 10 。导线外包有厚度为 0.5mm 导热系数为0.15W/(m ? K)的绝缘层。限定绝缘层的最高温度为 65°C,绝 缘层的外表面温度受环境影响，假设为40°C 。试确定该导线的最大允许电流为多 少？ 解：(1)以长度为L 的导线为例，导线通电后生成的热量为I 2RL ,其中的一部分 热量用于导线的升温，其热量为心务中：一部分热量通过绝热层的 导热传到大气中，其热量为：门二 1 , d In 2 L d 1 根据能量守恒定律知：l 2RL -门 述二厶E = I 2RL -门 即 E = — L dT m = I 2RL - t w1 _tw2 4 di 1 , d 2 In 2 L d 1 q v 、d 2t ——' 2 dx =-(7000)= 4000 50 二 200000 W/m 3 t w1 - t w2 。 2 q x 卫=4000.： 0 = 0 1、水以1.5m /s 的速度流过内径为25mm的加热管。管的内壁温度保持100℃，水的进口温度为15℃。若要使水的出口温度达到85℃，求单位管长换热量（不考虑修正）。已知50℃的水λf =0.648 W/(m.K)，νf =0.566×10-6m2/s，Pr =3.54。 2、取外掠平板边界层的流动由层流转化为湍流的临界雷诺数5×105,试计算25℃的空气和水达到临界雷诺数时所需要的平板长度,取u =1m/s，ν空气=15.53×10-6m2/s,ν水=0.905×10-6。 3、试推导努谢尔特关于层流膜状凝结的理论解 4、用实验测定一薄壁管流体平均对流换热系数。蒸汽在管外凝结并维持管内壁温度为100℃。水在管内流动流量为G=0.5Kg/s，水温从15℃升到45℃。管的内径d=50mm，长L=4.5m。试求管内流体与壁面间的平均换热系数。已知水在30℃时c p=4.174KJ/(Kg.K) 5、以0.8m/s 的流速在内径为2.5cm 的直管内流动，管子内表面温度为60℃，水的平均温度为30℃，管长2m ，试求水所吸收的热量。已知30℃时水的物性参数为：Pr =5.42，c p =4.17KJ/(Kg.K)，λ=61.8×10-2 W/(m.K)，ρ=995.7Kg/m 3，μ =80.15×10-6 Kg/(m.s)；水60℃ 时的ν=0.4699×10-6 m 2/s ，水在管内流动准则方程式为 4 .08.0Pr Re 027.0f f f Nu =，适用条件：Re f =104-1.2×105，Pr f =0.6-120，水与壁面间的换热温差Δt ≤30℃。 6、计算一空气横掠管束换热的空气预热器的对流换热量。已知管束有25排，每排12根光管，管外径25mm ，管长 1.5m ，叉排形式，横向管间距S 1=50mm ，纵向管间距S 2=38mm ，管壁温度120℃，空气来流速度u f =4m/s ，空气进口温度20℃，出口温度40℃。已知空气物性：λf =0.0267W/(m.K)，νf =16.0×10-6m 2/s ，Pr f =0.701。最大流速u max = u f S 1/(S 1-d)；推荐关联式：m w f f n f f c Nu ??? ? ??=Pr Pr Pr Re 36.0（公式适 用条件：N ≥20，光管管束，Pr f =0.7～500，除Pr w 的定性温度为壁温外，其余定性温度为流体在管束中的平均温度。指数m 对气体m =0，对液体m =0.25， 1. 传热学的发展概述 18世纪30年代首先从英国开始的工业革命促进了生产力的空前发展。生产力的发展为自然科学的发展成长开辟了广阔的道路。传热学这一门学科就是在这种大背景下发展成长起来的。导热和对流两种基本热量传递方式早为人们所认识，第三种热量传递方式则是在1803年发现了红外线才确认的，它就是热辐射方式。在批判“热素说”确认热是一种运动的过程中，科学史上的两个著名实验起着关键作用。其一是1798年伦福特(B ．T ．Rumford)钻炮筒大量发热的实验，其二是 1799年戴维(H ．Davy)两块冰块摩擦生热化为水的实验。确认热来源于物体本身内部的运动开辟了探求导热规律的途径。1804年毕渥根据实验提出了一个公式，认为每单位时间通过每单位面积的导热热量正比例于两侧表面温差，反比例于壁厚，比例系数是材料的物理性质。傅里叶于1822年发表了他的著名论著“热的解析理论”，成功地完成了创建导热理论的任务。他提出的导热定律正确概括了导热实验的结果，现称为傅里叶定律，奠定了导热理论的基础。他从傅里叶定律和能量守恒定律推出的导热微分方程是导热问题正确的数学描写，成为求解大多数工程导热问题的出发点。他所提出的采用无穷级数表示理论解的方法开辟了数学求解的新途径。傅里叶被公认为导热理论的奠基人。在傅里叶之后，导热理论求解的领域不断扩大。同样，自1823年M. Navier 提出流动方程以来，通过1845 年 G.G. Stokes 的改进，完成了流体流动基本方程的创建任务。流体流动理论是更加复杂的对流换热理论的必要前提，1909和1915年W. Nusselt 开辟了在无量纲数原则关系正确指导下，通过实验研究对流换热问题的一种基本方法。1904 年，L. Prandtl 提出的对流边界层理论使流动微分方程得到了简化，1921年 E. Pohlhausen 基于流动边界层理论引进了热边界层的概念，为对流传热微分方程的理论求解建立了基础。在辐射传热研究方面，19世纪J. Stefan 根据实验确定了黑体辐射力正比于它的绝对温度的四次方的规律，1900年M.Planck 提出的量子假说奠定了热辐射传热理论基础。上述传热理论为传热分析解析、数值以及实验研究奠定了理论基础。还要特别提到的是，由于计算机的迅速发展，用数值方法对传热问题的分析研究取得了重大进展，在20世纪70年代已经形成一个新兴分支—数值传热学。近年来，数值传热学得到了蓬勃的发展[2-4]。 2. 传热分析计算理论 热量传递主要有三种传递形式，分别是热传导、热对流和热辐射。热传导是指两个相互接触良好的物体之间的能量交换或一个物体由于其自身温度梯度而 引起的内部能量的传递。其遵循傅里叶定律[5]：dT q dx λ=-，其中λ是热导率， dT dx 是温度梯度，q 是热流密度。热对流是指在物体与其周围介质之间发生的热量交换。热对流分为自然对流和强制对流，用牛顿冷却方程描述为()w f q h t t =-，其中h 为表面传热系数，w t 为物体表面的温度，f t 为物体周围流体的温度。一个 物体或两个物体之间通过电磁波形式进行的能量传递交换称为热辐射，通常由斯 传热学经典计算题 热传导 1. 用热电偶测量气罐中气体的温度。热电偶的初始温度为20℃，与气体的表面传热系数为()210/W m K ?。热电偶近似为球形，直径为0.2mm 。试计算插入10s 后，热电偶的过余温度为初始过余温度的百分之几？要使温度计过余温度不大于初始过余温度的1％，至少需要多长时间？已知热电偶焊锡丝的()67/W m K λ=?，7310ρ= 3/kg m ，()228/c J kg K =?。 解： 先判断本题能否利用集总参数法。 3 5100.110 1.491067hR Bi λ--??===?＜0.1 可用集总参数法。 时间常数 3 73102280.110 5.563103c cV c R hA h ρρτ-??===?= s 则10 s 的相对过余温度 0θθ=exp c ττ??-= ???exp 1016.65.56??-= ???％ 热电偶过余温度不大于初始过余温度1％所需的时间，由题意 0θθ=exp c ττ??- ??? ≤0.01 exp 5.56τ?? - ???≤0.01 解得 τ≥25.6 s 1、空气以10m/s 速度外掠0.8m 长的平板，C t f 080=,C t w 030=，计算 该平板在临界雷诺数c e R 下的c h 、全板平均表面传热系数以及换热量。 （层流时平板表面局部努塞尔数 3/12/1332.0r e x P R Nu =，紊流时平板表面局部努塞尔数3/15/40296.0r e x P R Nu =，板宽为1m ，已知5105?=c e R ，定性 温度C t m 055=时的物性参数为： )/(1087.22K m W ??=-λ，s m /1046.1826-?=ν，697.0=r P ） 解：(1)根据临界雷诺数求解由层流转变到紊流时的临界长度 C t t t w f m 055)(21=+=,此时空气得物性参数为： )/(1087.22K m W ??=-λ，s m /1046.1826-?=ν，697.0=r P )(92.0101046.1810565m u R X ul R c c e c e =???==?=-ν ν 由于板长是0.8m ，所以，整个平板表面的边界层的流态皆为层流 ? ==3/12/1332.0r e x P R hl Nu λ)/(41.7697.0)105(8.01087.2332.0332.023/12/1523/12 /1C m W P R l h r e c c ?=????==-λ (2)板长为0.8m 时，整个平板表面的边界层的雷诺数为： 561033.41046.188.010?=??==-νul R e 全板平均表面传热系数: )/(9.13697.0)1033.4(8.01087.2664.0664.023/12/1523/12 /1C m W P R l h r e c ?=????==-λ 全板平均表面换热量W t t hA w f 9.557)3080(18.09.13)(=-???=-=Φ 关系式 返回到上一层以下汇总了工程中最常见的几类对流换热问题的对流换热计算关系式，适用边界条件，已定准则的适用范围，特征尺寸与定性温度的选取方法。 一、掠过平板的强迫对流换热 应注意区分层流和湍流两种流态 ( 一般忽略过渡流段 ) ，恒壁温与恒热流两种典型的边界条件，以及局部 Nu 数和平均 Nu 数。 沿平板强迫对流换热准则数关联式汇总 注意：定性温度为边界层的平均温度，即。 二、管内强迫对流换热 (1) 流动状况不同于外部流动的情形，无论层流或者湍流都存在流动入口段和充分发展段，两者的长度差别很大。计算管内流动和换热时，速度必须取为截面平均速度。 (2) 换热状况管内热边界层也同样存在入口段和充分发展段，只有在流体的 Pr 数大致等于 1 的时候，两个边界层的入口段才重合。理解并准确把握两种典型边界条件 ( 恒壁温与恒热流 ) 下流体截面平均温 度的沿程变化规律，对管内对流换热计算有着特殊重要的意义。 (3) 准则数方程式要注意区分不同关联式所针对的边界条件，因为层流对边界条件的敏感程度明显高于湍流时。还需要特别指出，绝大多数管内对流换热计算式 5f 对工程上的光滑管，如果遇到粗糙管，使用类比率关系式效果可能更好。下表汇总了不同流态和边界条件下管内强迫对流换热计算最常用的一些准则数关联式。 (4) 非圆截面管道仅湍流可以用当量直径的概念处理非圆截面管道的对流换热问题。层流时即使用当量直径的概念也无法将不同截面形状管道换热的计算式全部统一。 常热流 层流，充分发展段， 常壁温 层流，充分发展段， 充 - 充分发展段，气体， - 充分发展段，液体， ； 紊流，充分发展段， 《传热学》上机大作业 二维导热物体温度场的数值模拟 学校：西安交通大学 姓名：张晓璐 学号:10031133 班级：能动A06 一．问题（4-23） 有一个用砖砌成的长方形截面的冷空气通道，形状和截面尺寸如下图所示，假设在垂直纸面方向冷空气和砖墙的温度变化很小，差别可以近似的予以忽略。在下列两种情况下计算：砖墙横截面上的温度分布；垂直于纸面方向上的每米长度上通过墙砖上的导热量。 第一种情况：内外壁分别维持在10C ?和30C ? 第二种情况：内外壁与流体发生对流传热，且有C t f ?=101， )/(2021k m W h ?=，C t f ?=302，)/(422k m W h ?=，K m W ?=/53.0λ 二．问题分析 1.控制方程 02222=??+??y t x t 2.边界条件 所研究物体关于横轴和纵轴对称，所以只研究四分之一即可，如下图： 对上图所示各边界： 边界1：由对称性可知：此边界绝热，0=w q 。 边界2：情况一：第一类边界条件 C t w ?=10 情况二：第三类边界条件 )()( 11f w w w t t h n t q -=??-=λ 边界3：情况一：第一类边界条件 C t w ?=30 情况二：第三类边界条件 )()( 22f w w w t t h n t q -=??-=λ 三：区域离散化及公式推导 如下图所示，用一系列和坐标抽平行的相互间隔cm 10的网格线将所示区域离散化，每个交点可以看做节点，该节点的温度近似看做节点所在区域的平均温度。利用热平衡法列出各个节点温度的代数方程。 第一种情况： 内部角点： 储运与建筑工程学院能源与动力工程系 计算传热学课程大作业报告 作业题目：代数方程组的求解 学生姓名：田 学号： 专业班级：能动1 2017年9月23日 目录 一、计算题目 (3) 二、离散方程 (3) 三、程序设计 (4) 3.1 高斯赛德尔迭代法 (4) 3.2 TDMA法 (5) 四、程序及计算结果验证 (6) 五、网格独立性考核.................... 错误！未定义书签。 3.1 高斯赛德尔迭代法 (7) 3.2 TDMA法 (8) 六、结果分析与结论 (8) 3.1 高斯赛德尔迭代法 (9) 3.2 TDMA法 (10) 一、计算题目 分别用高斯赛德尔迭代和TDMA 方法求解方程 2 2dx d dx d u φφρΓ= （1） 在Γ u ρ=-5，-1,0,1,5情况下的解，并表示在图中。 其中，x =0，φ=0；x =1，φ=1. 二、离散方程 采用控制容积法： 即??Γ=e 22w e w dx d dx d u φφ ρ（2） ) )()(()2 2 ( w W P e P E p w p e x x u δφφδφφφφφφρ---Γ=+- +（3） 假设均分网格，则有x x x w e ?==)()(δδ 上式则变为： )2(2)(W P E W E u x φφφφφρ+-Γ=-?（4） 即11)2()2(4-+?+Γ+?-Γ=Γi i i u x u x φρφρφ（5） 11)421()421(-+Γ ?-+Γ?-=i i i u x u x φρφρφ（6） 三、程序设计 3.1 高斯赛德尔迭代法 由已知公式 11)421()421(-+Γ ?-+Γ?-=i i i u x u x φρφρφ可设计高斯赛德尔迭代C 语言程序如下： #include 对流换热公式汇总与分析 【摘要】流体与固体壁直接接触时所发生的热量传递过程，称为对流换热，它已不是基本传热方式。本文尝试对对流换热进行简单分类并对无相变对流换热公式简单汇总与分析。 【关键词】对流换热 类型 公式 适用范围 对流换热的基本计算形式——牛顿冷却公式： )(f w t t h q -= )/(2m W 或2Am 上热流量 )(f w t t h -=Φ )(W 上式中表面传热系数h 最为关键，表面传热系数是众多因素的函数，即 ),,,,,,,,(l c t t u f h p f w μαρλ= 综上所述，由于影响对流换热的因素很多，因此对流换热的分析与计算将分类进行，本文所涉及的典型换热类型如表1所示。 表1典型换热类型 1. 受迫对流换热 1.1 内部流动 对流换热 无相变换热 受迫对流换热 内部流动换热 圆管内受迫流动 非圆形管内受迫流动 外部流动 外掠平板 外掠单管 外掠管束（光管；翅片管） 自然对流换热 无限空间 竖壁；竖管 横管 水平壁（上表面与下表面） 有限空间 夹层空间 混合对流换热 — — — — 受迫对流与自然对流并存 相变换热 凝结换热 垂直壁凝结换热 水平单圆管及管束外凝结换热 管内凝结换热 沸腾换热 大空间沸腾换热 管内沸腾换热（横管、竖管等） 1.1.1 圆管内受迫对流换热 (1)层流换热公式 西德和塔特提出的常壁温层流换热关联式为 14 .03/13/13/1)()(Pr Re 86.1w f f f f l d Nu μμ= 或写成 14 .03/1)()(86.1w f f f l d Pe Nu μμ= 式中引用了几何参数准则 l d ，以考虑进口段的影响。 适用范围：16700Pr 48.0<<，75.9)(0044.0< 5-2 解：根据课本p158式（5—1a ）得一维稳态无源项的对流－扩散方程如下所示： 2 2x x u ??Γ =??φ φρ （取常物性） 边界条件如下： L L x x φφφφ====,; ,00 由（5—2）得方程的精确解为： 1 1)/(00--=--?Pe L x Pe L e e φφφφ Γ=/uL Pe ρ 将L 分成15等份，有：?=P Pe 15 对于中心差分、一阶迎风、混合格式和QUICK 格式分别分析如下： 1） (CD)中心差分 节点离散方程： 2 )5.01()5.01(1 1-?+?++-=i i i P P φφφ 10,2 =i 2） 一阶迎风 节点离散方程： ? -?++++=P P i i i 2)1(1 1φφφ 10,2 =i 3） 混合格式 当1=?P 时，节点离散方程：2 )5.01()5.01(1 1-?+?++-= i i i P P φφφ ，10,2 =i 当10,5=?P 时，节点离散方程： 1-=i i φφ ， 10,2 =i 4） QUICK 格式，节点离散方程： ??? ???--++++++= +-?? -??+?)336(8122121 1111i i i i i i P P P P P φφφφφφ， 2=i ?? ????---++++++= +--? ? -??+?)35(8122121 12111i i i i i i i P P P P P φφφφφφφ， 2≠i 用matlab 编程如下：（本程序在x/L=0-1范围内取16个节点进行离散计算，假设y(1)= 0φ=0,y(16)=L φ=1,程序中Pa 为?P ，x 为题中所提的x/L 。由于本程序假设 y(1)=0φ=0,y(16)=L φ=1，所以 y y y y y y L =--=--=--0 10 )1()16()1(00φφφφ） Pa=input('请输入Pa=') x=0:1/15:1 Pe=15*Pa; y=(exp(Pe*x)-1)/(exp(Pe)-1) plot(x,y,'-*k') %精确解 hold on y(1)=0,y(16)=1; for i=2:15 y(i)=((1+0.5*Pa)*y(i-1)+(1-0.5*Pa)*y(i+1))/2; end plot(x,y(1:16),'-or') %中心差分 hold on for i=2:15 y(i)=((1+Pa)*y(i-1)+y(i+1))/(2+Pa); end plot(x,y(1:16),'-.>g') %一阶迎风 hold on for i=2:15 if Pa==1 y(i)=((1+0.5*Pa)*y(i-1)+(1-0.5*Pa)*y(i+1))/2; else y(i)=y(i-1) end end plot(x,y(1:16),'-+y') %混合格式 hold on for i=2:15 if i==2 y(i)=y(i+1)/(2+Pa)+(1+Pa)*y(i-1)/(2+Pa)+(Pa/(2+Pa))*(6*y(i)-3*y(i-1)-3*y(i+1))/8 else y(i)=y(i+1)/(2+Pa)+(1+Pa)*y(i-1)/(2+Pa)+(Pa/(2+Pa))*(5*y(i)-y(i-1)-y(i-2)-3*y(i+1))/8 end end plot(x, y(1:16),'- 计算传热学课程设计（报告）题目:充满多孔介质的长方形截面通道内充分发展对流换热问题的数值研究 学生姓名：朱鹏齐尚超杨鹏来芦旭红 学号：10123106 10123107 10123108 10123103 专业班级：热能与动力工程10－1班 指导教师：黄善波巩亮 2013年 7 月 5 日 热工一班 组长：朱鹏组员：芦旭红，齐尚超，杨鹏来 目录 1.设计题目 (3) 1.1设计题目 (3) 1.2已知参数 (4) 2.物理与数学模型.. .................................... ..5 2.1物理模型 (5) 2.2数学模型 (5) 3.数值处理与程序设计 (6) 3.1数学模型无量纲化 (6) 3.2数值求解 (8) 3.3程序编写. (11) 4.程序的验证 (12) 5.计算结果与分析 (14) 6.结论 (21) 7.参考文献 (21) 8.附录 (22) 1 设计题目 （多孔介质，矩形a/b,单方程）水在一长方形截面的通道中进行充分发展的层流流动，该通道内充满多孔介质。多孔介质具有良好的强化换热能力，孔隙率ε是其基本结构参数，据此可以计算渗透率K ，惯性系数CF ，有效导热系数ke,具体表达式见[5]。其内部充满流体时的流动和换热通常采用体积平均法进行建模，即不考虑区域内孔的微结构而假定区域内任意一点处既有流体相又有固体相。由于金属泡沫的固体骨架导热系数较高，因此对于其内部的对流换热，通常采用局部非平衡模型，即考虑区域内流体温度和固体温度的差异。填充孔隙率为ε=0.6的多孔介质，渗透率表示为： 23 2 150(1)d K εε= - 惯性系数表示为： 23 1.75F C = 有效导热系数ke 表示为： (1)e f s k k k εε=+- 沿流动方向的速度方程可以简化为 222220f w w p w x y z K μμ ρε?????+---= ?????? （1） 截面上的平均流速为wm=0.1m/s,dp/dz 的值是恒定的，可以通过下式得到： 2d d m f m p w z K μ ρ=-- （2） 其中，w 为沿流动方向的速度。换热方程为： 22222()e w m T T w a b k q x y w ab ????++=?? ????? （3） 其中z 为轴向。假设流动和换热都达到充分发展，外壁面为恒热流边界条件（qw=1000W/2 m ）,请基于局部热平衡模型，选取a=0.04m,b=0.02m,以20℃作为水物性的参考温度，参考数值传热学（陶 文铨著）4.8节的内容 1.1 设计题目 1. 分析孔隙率对渗透率K ，惯性系数CF ，有效导热系数ke 的影响规律，计算 dp dz 随m w 的变化关系； 2. 通过能量守恒将方程（3）化为更简单的形式进而消去z （仅对恒热流条件下实施）， 对流动和换热方程进行无量纲化处理；计算传热学中国石油大学(华东)第四章大作业
数值传热学陶文铨第四章作业
传热学练习题
传热学作业
传热学计算题
《传热学》复习题
传热学计算例题
传热学 热对流 计算 (1)
计算传热学-传热基本原理及其有限元应用
传热学经典计算题
对流换热计算式
西安交通大学传热学大作业
计算传热作业1
对流换热公式汇总与分析..
数值传热学第五章作业
计算传热学课程设计报告..