计算传热学
传热学计算总结
《传热学》计算题总结一、 题型导热换热、对流换热、辐射换热、换热器 二、 公式小结1、 平壁稳态导热 第一类边界条件： 1） 单层：xt t t t w w w δ121--=221/)(m W t t q w w -=δλ多层∑∑=+=+-=-=ni in n i i in R t t t t q 1,11111λλδ第三类边界条件：传热问题2112111h h t t q n i i f f ++-=∑=λδ单位W/m 22、 圆筒壁稳态导热第一类边界条件 单层：121121r r n r r nt t t t w w w =--()12212112212r r n l t t t t r r n lw w w w πλπλ-=-=Φ多层：∑=++-=Φn i ii i n w w r r n l t t 111,1121 λπ第三类边界条件：1211112121ln 2121+=+++-=∑n n i i if f l r h ri r r h t t q ππλπ单位：W/m3、 对流换热 牛顿冷却公式：[]W )(f w t t hA Φ-=吸放热热量（热对流）：tvc t t mc p f f p ∆=-=Φρ)(21平板对流换热表面换热系数h管内对流换热表面换热系数h ：n Nu Pr Re023.08.0=（紊流，流体被加热n=0.4,流体被冷却 n=0.3）对流换热解题步骤1）定性温度→查物性，下标f 由t f 确定，下标w 由t w 确定； 2）由Re 判断流态；3）据Re 选择准则关联式计算Nu f ； 4）计算h 。
注意：1）外掠平板定性温度tm=1/2(tw+tf);管内定性温度tf 2）外掠平板临界Re=5×105；管内临界Re=104 3）换热量据牛顿冷却公式计算。
4、辐射换热斯蒂芬－玻尔兹曼定律（四次方定律）: （黑体）两表面封闭体系的辐射换热量：（实际表面）几种特殊情况的简化式： （a ） X 1-2=1时：（其中一个表面为平面或凸表面）（b ）A 1=A 2 时：（两无限大平壁之间）(c) A 1/A 2≈0 时 （空腔与内包壁）遮热板：111)T T (21214241b 2,1-+-=εεσq5、换热器设计计算传热过程方程式m t kA ∆=Φ；minmax minmax t lnt t t t m ∆∆∆-∆=∆热平衡式)()(22221111t t c M t t c M '-''=''-'=Φ， 其中M 为质量流量kg/s,c 为定压比热，由对应算术平均温度确定。
计算重点 公式 传热学
w
T0
’
T0
绝热滞止过程：
气体在绝热流动过程中，因受到某种阻碍流速 降为零的过程。 在绝热滞止时的温度和压力称为滞止温度T0和 滞止压力p0。若过程为定熵滞止过程：
T0 T
c2 f 2c p
k k 1
T0 p0 p T
3、临界压力比
pcr cr p0
2 cr ( ) k 1
对气体性质及过程特性没有规定。
二、流量计算[校核]
根据连续方程，喷管各截面的质量流量相等。 但各种形式喷管的流量大小都受最小截面控制，因 而通常按最小截面（收缩喷管的出口截面、缩放喷 管的喉部截面）来计算流量，即：
渐缩喷管：
qm qm
A2 c f 2 v2 Acr c fcr vcr
缩放喷管：
§5.2 可逆循环分析及其热效率
三、逆卡诺循环卡诺制冷循环
T
T0
制冷
T0 q1 s2
T2 s1 T0 c T2 c
s
1 T0 1 T2
Rc
q2 T2
w
§5.2 可逆循环分析及其热效率
逆卡诺循环卡诺制热循环
T
T1
制热
T1
T0 s1 T1 s2 s
q1
’
1 T0 1 T1
Rc q2
5.2 可逆循环分析及其热效率
热效率：
T2 S ab q2 t 1 1 q1 T1S ab
T2 1 c T1 由于n可以为任何自然数，所以，在T1 和T2之间的可逆循环有无数个。
熵变计算公式：
T 2 Rg ln p 2 △S 1 2 △Sad cp ln p1 T1
计算传热学7-1例题编程
计算传热学7-1例题编程由于没有提供具体的例题,我这里以一个简单的例子进行编程：计算一个半径为5cm的球体在30秒内从100摄氏度降到50摄氏度的冷却过程中,球体内每一时刻的温度值,球体的导热系数为0.05。
首先,需要找到球体内每一时刻的温度变化率。
根据热传导方程,温度变化率$dT/dt$与导热系数$K$、球体的半径$r$、球体内每一时刻的温度$T$、球体的质量$m$、球体的比热容$c$、球体的表面积$A$有关,计算公式如下：$dT/dt = -(K/(mr*c))*(A/(4*pi*r^2))*(T-T0)$其中,$T0$为环境温度,本例中为50摄氏度。
可以将上述公式转化为程序：```pythonimport mathK = 0.05 #导热系数r = 0.05 #半径,单位为米m = (4/3)*math.pi*pow(r,3)*7850 #质量,假设密度为7850kg/m^3c = 480 #比热容,单位为J/(kg*K)A = 4*math.pi*pow(r,2) #表面积T0 = 50 #环境温度T = 100 #初始温度t = 0 #初始时间dt = 0.1 #时间间隔,单位为秒while t <= 30:dTdt = -(K/(m*c))*(A/(4*math.pi*pow(r,2)))*(T-T0)T = T + dTdt*dtt = t + dtprint("Time:", t, "Temperature:", round(T,2))```执行上述程序,可以得到每时刻的温度值输出结果,单位为摄氏度。
根据该程序,可以修改参数来计算不同条件下的传热过程。
传热学计算公式
Nu = 2+0.6(Re^1/2)(Pr^1/3) 。
F=Q/kK*△tm F 是换热器的有效换热面积。
Q 是总的换热量。
k 是污垢系数一般取0.8-0.9K。
是传热系数。
△tm 是对数平均温差。
传热学三种传热方式可以分开学。
传热学相较于理论力学，工程热力学，流体力学而言还是比较简单的，一般大学生掌握了高等数学完全可以自学的。
学习传热学必须有耐心，了解几种换热方式和常见的几个常数公式（努谢尔特数、格拉晓夫数、伯努利常数，傅里叶常数，而且常常推导下几个常用常数公式间的关系，你会惊奇地发现他们其实不少是远亲的），其实解决传热学问题绝大多数都是在和导热系数较劲，有时候是直接涉及。
扩展资料：
在热对流方面，英国科学家牛顿于1701年在估算烧红铁棒的温度时，提出了被后人称为牛顿冷却定律的数学表达式，不过它并没有揭示出对流换热的机理。
传热学作为学科形成于19世纪。
1804年，法国物理学家毕奥在热传导方面得出的平壁导热实验结果是导热定律的最早表述。
稍后，法国的傅里叶运用数理方法，更准确地把它表述为后来称为傅里叶定律的微分形式。
1860年，基尔霍夫通过人造空腔模拟绝对黑体，论证了在相同温度下以黑体的辐射率(黑度)为最大，并指出物体的辐射率与同温度下该物体的吸收率相等，被后人称为基尔霍夫定律。
传热学计算公式范文
传热学计算公式范文传热学是物理学的一个分支，研究能量在物体之间的传递过程。
在传热学中，有许多重要的计算公式可以用于解决热传导、对流和辐射等传热现象。
下面将介绍一些常见的传热学计算公式。
热传导是物质内部由高温区向低温区传递热量的过程。
热传导热量的大小与物体的温度差、物体的热导率以及物体的尺寸等因素有关。
下面是一些常用的热传导计算公式：1.热流密度公式:热流密度（q）是单位时间内通过单位面积的热量传递量，可以由下式计算：q = -k * (dT/dx)其中，k是物体的热导率，dT/dx是温度梯度。
2.热传导率（k）：物体的热传导率是描述物质导热能力的物理量，可以用以下公式计算：k=Q*L/(A*ΔT)其中，Q是通过物体的热量，L是物体的长度，A是传热的横截面积，ΔT是温度差。
3.热阻（R）：热阻是描述物质阻碍热传导的程度的物理量，可以用以下公式计算：R=L/(k*A)其中，L是物体的长度，k是物体的热导率，A是传热的横截面积。
对流是物体表面与流体之间的热传递方式，流体通过对流来接触物体表面并将热量带走。
对于对流传热的计算，常用的公式有：1.流体的对流换热公式:流体通过对流来接触物体表面并带走热量，可以由下式计算：q = h * A * (T - Tfluid)其中，h是对流换热系数，A是物体表面积，T是物体表面的温度，Tfluid是流体的温度。
2.对流换热系数（h）：对流换热系数描述了流体的传热能力，它可以由以下公式计算：h=(Nu*k__)/L其中，Nu是Nusselt数，k__是流体的导热系数，L是流体经过的长度。
3. Nusselt数（Nu）：Nusselt数描述了流动体系中传热性能的参数，可以通过以下公式计算：Nu=(h*L)/k__其中，h是对流换热系数，L是流体经过的长度，k__是流体的导热系数。
辐射传热是物体通过辐射来传递能量的过程，对于辐射传热的计算，常用的公式有：1.斯特藩-玻尔兹曼定律:斯特藩-玻尔兹曼定律描述了黑体辐射能量的传递率，可以用下式表示：q=σ*ε*A*(T1^4-T2^4)其中，σ是斯特藩-玻尔兹曼常数，ε是物体的辐射率，A是物体的面积，T1和T2是物体的温度。
计算传热学第一类边界条件离散方程
计算传热学第一类边界条件离散方程
传热学中的第一类边界条件，也称为热流量边界条件，是指在一个热传导过程中，表面的热流量（或热通量）是已知的。
对于此类边界条件，我们可以通过离散化的方法来建立传热方程，并求解出时间和空间上的温度分布。
以一维传热问题为例，考虑一个长L的直管，管的外部温度为T1，内部温度为T2，管壁的热流量为q，假设热传导系数为k，管的长度分为n个小段，依据传热学原理可以建立如下的离散方程：
$$\frac{T_{i-1}-2T_i+T_{i+1}}{\Delta x^2}=-\frac{q}{k},\quad i=1,2,\cdots,n-1, $$
其中$T_i$是第i个网格（或节点）的温度。
将第一类边界条件代入方程，得到：
$$\frac{T_1-2T_2+T_3}{\Delta x^2}=-\frac{q}{k}, $$
$$\frac{T_{n-2}-2T_{n-1}+T_n}{\Delta x^2}=-\frac{q}{k}, $$
同时，我们还需要加入一个边界条件，即：
$$T_1=T_2=T_1, $$
$$T_n=T_{n-1}=T_2. $$
通过将这些方程进行整理和线性求解，可以求得各个节点的温度。
具体的解法可以用三对角矩阵算法（Thomas Algorithm）进行求解，也可以使用迭代法（如高斯-赛德尔迭代法或松弛迭代法）进行求解。
总之，对于第一类边界条件，我们可以通过将离散化的方法来建立方程并求解，这是解决传热学问题中的一种常用方法。
计算重点公式传热学
计算重点公式传热学传热学是研究热能在物质之间传递的学科，涵盖了热传导、热对流和热辐射三种传热方式。
在工程和科学领域中，计算传热是非常重要的，可以用来优化和设计各种热能设备和系统。
下面将介绍一些重要的传热计算公式。
1.热传导计算公式热传导是通过分子间的相互作用传递热能的方式。
对于常见的一维热传导问题，可以使用傅里叶热传导定律进行计算：q = -kA(dT/dx)其中，q是单位时间内通过物体的热量流率，k是物质的热导率，A 是传热截面积，dT/dx是温度梯度。
如果传热是在不同的材料之间进行，还需要考虑热传导的界面热阻。
界面热阻的计算公式为：R=1/(hA)其中，R是界面热阻，h是对流传热系数。
2.热对流计算公式热对流是通过流体的对流传递热能的方式。
对于流体中的对流传热，可以使用牛顿冷却定律进行计算：q=hAΔT其中，q是单位时间内通过物体的热量流率，h是对流传热系数，A 是传热表面积，ΔT是流体和物体之间的温度差。
对流传热系数h可以通过实验测量或者经验公式进行估算，常用的计算公式有Nusselt数和普朗特数。
3.热辐射计算公式热辐射是通过物体表面的电磁辐射传递热能的方式。
对于黑体辐射，可以使用斯特藩—玻尔兹曼定律进行计算：q=σAε(T^4)其中，q是单位时间内通过物体的热量流率，σ是斯特藩—玻尔兹曼常数，A是物体的表面积，ε是物体的辐射率，T是物体的温度。
对于非黑体的辐射传热，还需要考虑辐射率和视觉系数等因素。
4.综合传热计算在实际问题中，常常会有多种传热方式同时存在。
此时，需要将不同传热方式的热流量进行累加，得到总的传热量。
根据能量守恒定律，可以得到以下综合传热公式：q_total = q_conduction + q_convection + q_radiation其中，q_total是总的热量流率，q_conduction是热传导的热量流率，q_convection是热对流的热量流率，q_radiation是热辐射的热量流率。
计算传热学在工程领域的应用研究
计算传热学在工程领域的应用研究引言计算传热学是研究热传输过程中热量、温度、流速等参数变化规律的一门学科。
在工程领域中，计算传热学在热能转换、能源利用与环境保护等方面发挥着重要的作用。
本文将从传热模型、传热计算方法和应用案例三个方面阐述计算传热学在工程领域的应用研究。
第一章传热模型传热模型是计算传热学中的基本概念，它描述了热量从高温区向低温区传递的过程。
传热模型可分为对流传热、导热传热和辐射传热三种。
对流传热指的是流体在物体表面与物体接触的同时将热量带走的现象。
导热传热则是指物体自身内部传递热能的过程。
辐射传热则是指物体表面和环境之间的热辐射现象。
第二章传热计算方法传热计算方法指的是以数学公式和模型来对传热过程进行计算和模拟的方法。
传热计算方法主要为传热系数法、有限元法、有限差分法以及计算流体力学等方法。
传热系数法可用于估算流体和固体间的传热系数，从而计算热传输过程。
有限元法和有限差分法则是通过分块分割的方法对物体进行离散化处理，并建立数学模型进行传热计算。
计算流体力学则是通过求解流体运动的数学方程来模拟流体在传热过程中的行为。
第三章应用案例在工程领域中，计算传热学的应用非常广泛。
以下列出几个应用案例。
1. 燃烧工程：在燃烧工程中，通过计算传热学模型和计算热转换效率，可以选择最适合的燃烧工艺和优化燃烧效率。
2. 电子工程：在电子设备中，通过计算传热学和优化散热系统的设计可以有效降低电子设备的温度，提高电子设备的可靠性。
3. 石化工程：在石化工程中，通过计算传热学和流体力学分析可以优化反应器的设计，提高生产效率。
4. 汽车工程：在汽车工程中，通过计算传热学和流体力学分析可以优化发动机散热系统的设计，提高汽车发动机的性能。
结论综上所述，计算传热学在工程领域中的应用研究发挥着重要的作用。
通过传热模型的建立和传热计算方法的运用，可以解决一些复杂的传热问题，并优化工程设计，提高工程效率和质量。
传热学热传导公式
传热学热传导公式
热传导的公式是：ut=ku。
热传导是介质内无宏观运动时的传热现象，其在固体、液体和气体中均可发生，但严格而言，只有在固体中才是纯粹的热传导，而流体即使处于静止状态，其中也会由于温度梯度所造成的密度差而产生自然对流，因此，在流体中热对流与热传导同时发生。
通常使用傅里叶定律来计算：Q = -kA(dT/dx)，其中，Q为单位时间内通
过某一面积的热量流(单位为瓦特W)、k为物质的热传导系数(单位为瓦特/米·开尔文W/(m·K))、A为热源和热汇之间的接触面积(单位为平方米m²)、dT/dx为温度梯度(单位为开尔文/K)，表示在长度为x的方向上，温度变化的速率。
以上内容仅供参考，建议查阅传热学书籍或咨询专业人士获取更准确的信息。
数值传热学
数值传热学数值传热又称计算传热，是传热学与数值方法相结合的一门交叉学科，它采用数值方法描述流动和传热问题的控制方程，并用计算机求解。
数值换热，其基本思想是将原始坐标在空间和时间上连续的物理量场（如速度场、温度场和浓度场等），用一系列有限个离散点上的数值来代替，通过一定的原理建立离散点变量值之间的关系代数方程（称为离散方程）。
通过求解所建立的代数方程组，得到求解变量的近似值。
1简介数值传热学（numerical heat transfer）数值传热学，又称计算传热学，是指对描写流动与传热问题的控制方程采用数值方法，通过计算机求解的一门传热学与数值方法相结合的交叉学科。
数值传热学的基本思想是把原来在空间与时间坐标中连续的物理量的场（如速度场，温度场，浓度场等），用一系列有限个离散点上的值的集合来代替，通过一定的原则建立起这些离散点变量值之间关系的代数方程（称为离散方程）。
求解所建立起来的代数方程已获得求解变量的近似值。
2发展简史数值传热学，主要由20世纪中叶，S.V. Patankar和D.B.Spalding 等人在总结前人的研究基础上所提出。
E.M.Sparrow对数值传热学的发展也起到了一定的促进作用。
国内比较知名的学者是陶文铨教授。
陶文铨3研究方法数值传热学常用的数值方法1.有限差分法历史上最早采用的数值方法，对简单几何形状中的流动与换热问题最容易实施的数值方法。
其基本点是：将求解区域中用于坐标轴平行的一系列网格的交点所组成的点的集合来代替，在每个节点上，将控制方程中每一个导数用相应的差分表达式来代替，从而在每个节点上，形成一个代数方程，每个方程中包括了本节点及其附近一些节点上的未知值，求解这些代数方程就获得了所需的数值解。
2.有限容积法将所计算的区域划分成一系列控制容积划分为一系列控制容积，每个控制容积都有一个节点做代表。
通过将守恒型的控制方程对控制容积坐积分导出离散方程。
在导出过程中，需要对界面上的被求函数本身及其一阶导数的构成做出假定，是目前流动与换热问题的数值计算中应用最广的一种方法。
传热学公式总结
传热学公式总结在物理学中，传热学是一个重要的分支领域，研究物质之间热量的传递方式和规律。
在实际应用中，我们常常需要利用传热学公式来计算热传导、对流和辐射等过程中的热量变化。
本文将对传热学中常用的公式进行总结和归纳，帮助读者更好地理解和应用相关知识。
1. 热传导方程热传导是物质内部由于温度差异而引起的热量传递过程。
热传导的速率可以根据傅里叶定律描述：q = -kA(dT/dx)其中，q表示单位时间内通过横截面A传导的热量，k为材料的热导率，dT/dx表示单位长度内温度的变化率。
这个公式说明了热量传导与温度梯度之间的关系，温度梯度越大，热传导速率就越大。
2. 热对流公式热对流是通过流体介质的热传递方式，常见于气体和液体中。
热对流可以根据牛顿冷却定律进行计算：q = hA(Ts - T∞)其中，q表示通过表面积A从物体表面传递的热量，h为热对流系数，Ts为表面温度，T∞为流体的远场/环境温度。
牛顿冷却定律的基本思想是热量传递与温度差和表面积之间成正比，而且逆向传热过程中的温度差往往比较小。
3. 辐射传热公式辐射传热是通过电磁波辐射的方式进行的，不需要物质介质。
具体的辐射传热公式可以根据斯特藩-玻尔兹曼定律给出：q = εσA(T⁴s - T⁴∞)其中，q为单位时间内通过表面积A传递的辐射热量，ε为发射率（表征表面辐射能力的一种无量纲值），σ为斯特藩-玻尔兹曼常数，Ts为表面温度，T∞为远场/环境温度。
斯特藩-玻尔兹曼定律说明了辐射热量与表面温度的四次方成正比，这意味着一个小的温度提高可以显著增加辐射传热率。
4. 复合热传递在实际情况中，热传递往往是多种传热方式的复合过程。
例如，一个物体既有热传导，又有对流和辐射。
在这种情况下，总的热传递可以通过下列公式求得：q = q₁ + q₂ + q₃其中，q₁、q₂和q₃分别表示通过热传导、热对流和辐射传递的热量。
根据具体情况，我们可以使用以上公式中的一个或多个来计算总的热传递。
传热学第十章传热过程和换热器计算
例题 2
某逆流套管式换热器，刚投入工作时的运行参数为：
t1 360C,t1 300C,t2 30C,t2 200C 已知 qm1cp1=2500 W/K, k = 800 W/(m2.K)。运行一年后发现， 在 qm1cp1，qm2cp2，及入口温度不变的情况下，由于积垢使 得冷流体只能加热到162℃. 确定此情况的
(d)由式 kAtm 求出换热量 ;
(e)比较 与 ，如果相差较大，再重新假设流体出口温度， 重复上述计算，直到满意为止。
10.5 传热的强化与削弱（自学）
传热工程技术是根据现代工业生产和科学实践的需要而发展 起来的科学与工程技术，其主要任务是按照工业生产和科学 实践的要求来控制和优化热量传递过程。
和换热量 。
计算步骤：
(a) 先假设一个流体的出口温度，热平 衡方程式求出换热量 和 另一个流体的 出口温度；
kAtm
qm1cp1 t1 t1
qm2cp2 t2 t2
(b) 根据流体的进、出口4个温度求平均温差 tm ；
(c) 计算换热面两侧的表面传热系数 h1, h，2 进而求得总传热系数k；
tm (tm )ctf
教程中图10-23~10-26分别给出了管壳式换热器和交叉流式 换热器的 。
值取决于无量纲参数 P和 R： P tc tc , th tc
R th th tc tc
式中：下标h、c分别表示两种流体，上角标 ` 表示进口，`` 表示出口，图表中均以P为横坐标，R为参量。
1. 通过平壁的传热
K
1
1
1
h1 h2
KAt
说明： (1) h1和h2的计算；（2）如果计及辐射，换热系 数应该采用等效换热系数(总表面传热系数)
传热学 数值计算
数值计算
4、一无限大平壁厚度为 0.3m，其导热系数为λ =36.4 W/ (m·K)。平壁两侧表面均给定 为第三类边界条件，即 h 1 =60W/ (m 2 ·K)，t f 1 =25℃；h 2 =300W/ (m 2 ·K)，t f 2 =215℃。 当平壁中具有均匀的内热源 q v =2 10 W / m 时，试计算沿平壁厚度的稳态温度分布。
k
0
k k 1
11
k 2Fo (t 10 t f Bi ) (1 2Fo 2Bi Fo) t11 t11
0
四、
计算过程
⑴ 设定初值：
t (1~11) 35 ℃;
36W / m k ;
Bi h / ;
根据不同的 Fo 计算Δ τ ：
qv x 2
)/2
h2 x
h2 x TRB
qv x 2
) /(1
)
|T[i]-t[i]|<=EPS
NO
IT=IT+1
YES
打印“t[i]” ， “IT”
YES
IT>K
NO
停止
⑷ 程序与计算结果 #include"iostream.h" #include"iomanip.h" #include"math.h" #include"stdio.h" #define N 16 void main() { int M,i,IT,flag; //定义节点个数
计算结果：
各节点温度： 节点 1 2 411.24 10 422.65 3 420.33 11 414.23 4 427.22 12 403.62 5
传热学计算题
传热学计算题传热学是研究热的传递过程以及与其他物质和能量之间的相互作用的学科。
在工程领域，传热学计算是非常重要的，可以帮助工程师设计和改进各种设备和系统。
本文将通过一个计算示例来探讨传热学的一些基本概念和计算方法。
假设我们有一个长方体的金属块，其尺寸为2m×1m×0.5m，温度为100℃。
我们想要计算该块体的传热过程，以确定其表面和内部的温度分布。
首先，我们需要确定金属块的热传导性能。
热导率（λ）是描述材料传热性能的物理量，单位为W/(m·K)。
根据金属的种类和温度，我们可以查找到对应的热导率数值。
假设该金属块是铜材料，其温度为100℃时的热导率为400 W/(m·K)。
接下来，我们将应用热传导定律来计算金属块内部的温度分布。
热传导定律可以表示为：q = -λ × A × (∂T/∂x)其中，q为单位时间内通过单位面积的热量传递率，λ为热导率，A 为面积，∂T/∂x为温度梯度。
我们可以将金属块分成若干个切片，以每个切片为单位进行计算。
假设我们将其分成10个切片，每个切片的厚度为0.05m。
首先，我们计算金属块表面的热量传递率。
由于金属块的表面积为2×(2×1+1×0.5+2×0.5) = 14m^2，我们可以将其代入热传导定律的公式中。
q_surface = -λ × A_surface × (∂T/∂x)其中，A_surface为表面积，(∂T/∂x)为温度梯度。
由于金属块的表面温度为100℃，环境温度为25℃，温度梯度为(100-25)/0.05 = 1500 K/m。
将数值代入公式，我们可以计算得到表面的热量传递率为：q_surface = -400 × 14 × 1500 = -8.4 × 10^6 W接下来，我们计算金属块内部的温度分布。
数值传热学
数值传热学数值传热学（numerical heat transfer）数值传热学，又称计算传热学，是指对描写流动与传热问题的控制方程采用数值方法，通过计算机求解的一门传热学与数值方法相结合的交叉学科。
数值传热学的基本思想是把原来在空间与时间坐标中连续的物理量的场（如速度场，温度场，浓度场等），用一系列有限个离散点上的值的集合来代替，通过一定的原则建立起这些离散点变量值之间关系的代数方程（称为离散方程）。
求解所建立起来的代数方程已获得求解变量的近似值。
数值传热学（numerical heat transfer）数值传热学，又称计算传热学，是指对描写流动与传热问题的控制方程采用数值方法，通过计算机求解的一门传热学与数值方法相结合的交叉学科。
数值传热学的基本思想是把原来在空间与时间坐标中连续的物理量的场（如速度场，温度场，浓度场等），用一系列有限个离散点上的值的集合来代替，通过一定的原则建立起这些离散点变量值之间关系的代数方程（称为离散方程）。
求解所建立起来的代数方程已获得求解变量的近似值。
数值传热学常用的数值方法1.有限差分法历史上最早采用的数值方法，对简单几何形状中的流动与换热问题最容易实施的数值方法。
其基本点是：将求解区域中用于坐标轴平行的一系列网格的交点所组成的点的集合来代替，在每个节点上，将控制方程中每一个导数用相应的差分表达式来代替，从而在每个节点上，形成一个代数方程，每个方程中包括了本节点及其附近一些节点上的未知值，求解这些代数方程就获得了所需的数值解。
2.有限容积法将所计算的区域划分成一系列控制容积划分为一系列控制容积，每个控制容积都有一个节点做代表。
通过将守恒型的控制方程对控制容积坐积分导出离散方程。
在导出过程中，需要对界面上的被求函数本身及其一阶导数的构成做出假定，是目前流动与换热问题的数值计算中应用最广的一种方法。
3.有限元法把计算区域划分为一系列原题（在二维情况下，元体多为三角形或四边形），由每个元体上去数个点作为节点，然后通过对控制方程做积分来获得离散方程。
传热学计算题
传热学计算题答案需要依据具体的题目和具体的数据来确定。
以下是一个简单的传热学计算题的示例，并提供了一些可能的解答方法。
题目：在一个直径为6英寸的铜制圆形散热器中，一个1kW的电热元件被放置在中心位置。
假设散热器的热传导系数为100W/m·K，环境温度为25℃，要求在1小时内将电热元件产生的热量散发到周围环境中。
解答：首先，我们需要根据传热方程来计算散热器的散热量。
假设散热器的厚度为d，那么根据传热方程：散热量= 热源功率+ 壁面对流散热量+ 热传导散热量，其中：散热量= 散热面积×散热系数×温差电热元件产生的热量为1kW，即1000W，环境温度为25℃，因此可以列出以下方程：1000 + 2πr2d ×k ×(T - 25) = Q其中，r为散热器半径，d为散热器厚度，k为铜的热传导系数，T为散热器表面温度。
根据题目中给出的铜制圆形散热器的直径和厚度，可以计算出散热面积和散热器表面温度。
假设散热器表面温度与环境温度相同，即T = 25℃。
代入数据后可得：Q = 2π×(6英寸/2)2 ×π×100W/m2 ×(25℃- 25℃) = 3600W即每小时散发出3600瓦的热量。
由于散热器直径为6英寸，厚度为未知数d，因此需要求解方程2πr2d ×k = Q。
解得：d = 0.78mm。
因此，散热器的厚度为0.78毫米。
考虑到实际情况中散热器的厚度通常要大于这个数值，所以这里给出的数值是一个近似值。
另一种可能的解答方法是利用热力学基本公式来进行计算。
假设散热器的表面传热系数为Nu，即努塞尔数，它可以由空气流动和温度分布等因素决定。
在某些情况下，可以使用Nu 数来简化传热计算。
在这种情况下，可以根据热力学基本公式：热量= 质量流量×比热×温度变化来计算散热量。
已知电热元件产生的热量为1kW，环境温度为25℃，铜制圆形散热器的直径为6英寸，即15.24厘米，厚度为未知数d。
计算传热学第9节-第3章 二维导热1
计算传热学
3 二维导热
二维稳态导热有限体积方程
0 非稳态项 （体平均量） =
((λ (әT/әx))e-(λ (әT/әx))w)∆y∆z +((λ (әT/әy))n-(λ (әT/әy))s)∆x∆z
+ SP ∆x∆y∆z
N n ∆y W w P s S y z x ∆x e E
源项 扩散项 （体平均量） （面平均量） (λ e(TE-TP)/δ xe-λ w(TP-TW)/δ xw)∆y∆z +(λ n(TN-TP)/δ yn-λ s(TP-TS)/δ ys)∆x∆z+ SP∆x∆y∆z = 0 (λ e/δ xe (TE-TP)-λ w/δ xw (TP-TW))∆y +(λ n/δ yn (TN-TP) –λ s/δ ys (TP-TS))∆x+ SP∆x∆y = 0 (λ e/δ xe)∆y TE- (λ e/δ xe)∆y) TP -(λ w/δ xw )∆y TP+ (λ w/δ xw )∆y TW +(λ n/δ yn)∆x TN- (λ n/δ yn)∆x TP -(λ s/δ ys)∆x TP+ (λ s/δ ys)∆x TS + SP∆x∆y = 0
10
计算传热学
3 二维导热
二维稳态导热有限体积方程
0 非稳态项 （体平均量） =
((λ (әT/әx))e-(λ (әT/әx))w)∆y∆z +((λ (әT/әy))n-(λ (әT/әy))s)∆x∆z
扩散项 （面平均量）
+ SP ∆x∆y∆z 源项 （体平均量）
N n ∆y W w P s S y z x ∆x e E
计算传热学概念
计算传热学（Computational Heat Transfer，CHT）又称数值传热学（Numerical Heat Transfer，NHT）:是指对描写流动与传热问题的控制方程采用数值解法通过计算机予以求解的一门学科。
数值解法是一种离散近似的计算方法。
它所能获得的解不象分析解那样是被研究区域中未知量的连续函数，而只是某些代表性的点（称为节点）上的近似值。
有限差分法：有限容积法：有限元法：流体热传导（heat transfer）：当流体中存在着温度差时，温度高的地方将向温度低的地方传送热量的现象。
扩散（diffusion）：当流体混合物存在着组元的浓度差时，浓度高的地方将向浓度低的地方输送该组元物质的现象。
守恒型方程:如果微分方程在任意有限区域V内积分所得表达式都表示该区域的守恒定律，则为守恒型方程。
否则为非守恒型方程.抛物型问题：特点：这类问题中因变量与时间有关，描写了物理上的非稳态导热问题。
离散方程：步进求解，即从已知的（或已求解出的）某一时层上的值出发，根据边界条件，将解一步一步向前推进。
椭圆型问题特点：描写了物理上的稳态问题，求解区域内各点之值是互相影响的。
如稳态导热过程，有回流的流动与换热。
离散方程：联立求解（整场求解）空间区域离散化（domain discretization）:实质：用有限个离散的点代替原来的连续空间。
实施：计算区域划分多个子区域（sub-domain）,定其节点位置及节点所代表的控制容积（control volume）。
4种几何要素：网格线、节点、控制容积、界面内接点法和外节点法的比较：(1) 边界节点所代表的控制容积不同(2)当网格不均分时，节点位置不同（3）当网格不均分时，界面位置不同控制容积积分法：三步：1、将守恒型控制方程在控制容积中及△t内对空间和时间积分；2、选择未知函数及其导数对空间及时间的分布曲线3、按选定的型线作出积分，并整理成关于节点上未知值的代数方程显式：如果在整个时间步长内均取初始时刻之值而仅在该步长的结束时刻取终了之值，为显式，反之为隐式。
传热学计算题
传热学计算题
当涉及到传热学的计算题时，通常会涉及热传导、对流和辐射等传热方式。
以下是一个简单的计算题示例：
问题：一块厚度为2 cm，面积为0.5 m²的铝板，其表面温度为100°C，远离其他热源的环境温度为20°C。
已知铝板的导热系数为 200 W/(m·K)，求铝板上的热流量。
解答：
首先，我们需要确定铝板两侧的温度差ΔT。
ΔT = 表面温度 - 环境温度 = (100°C - 20°C) = 80°C。
根据热传导的公式 Q = kAΔT/d，其中 Q 是热流量，k 是导热系数，A 是面积，ΔT 是温度差，d 是物体的厚度。
将已知值代入公式进行计算：
Q = (200 W/(m·K)) × (0.5 m²) × (80°C) / (0.02 m)
= 40000 W = 40 kW
因此，铝板上的热流量为 40 kW。
这个计算结果表示在给定的温度差
下，铝板单位时间内从高温端传输的能量量。
请注意，这只是一个基本的计算示例，实际的传热问题可能会更加复杂，涉及更多的因素和传热方式。
在解决具体的传热计算题时，需根据所给条件选择适当的公式，并注意单位的一致性。
传热学常用公式
传热学常用公式1、热传导热流量与热流密度的区别，前者是单位时间内通过传热面积的总热量，单位为W，后者是单位时间内通过单位面积的热量，单位为W/m2。
傅里叶公式：热阻：类似于电阻，可以把它看成是阻挡热量传递的阻力，热流量=温差（动力）/热阻。
热阻与导热系数成反比，热阻大，导热系数就小。
面积热阻：2、热对流（对流换热）热对流指的是流体层之间发生相对位移，冷热流体掺混产生热量传递。
而在生活及工程中，更常见的是对流换热，即流体与固体表面之间的热量传递，它们都包含有热传导和热对流两种传热方式。
同时需要注意的是，对流换热中的流体必须要处于流动状态，如果流体是静止的，那么它就变成了单纯的热传导了。
用于计算对流换热的为牛顿冷却公式：注意两个温度之间的温差单位可以为K或者摄氏度。
对流换热热阻为：3、热辐射（辐射换热）一切温度高于0K的物体都会以电磁波的形式向外发射出热量，物体在环境中不断的发射出电磁波，同时吸收其它物体发射过来的电磁波能量，这个综合过程称为辐射换热。
热辐射不需要介质，可以在真空中传播。
用于计算辐射换热量的公式为四次方定律公式，要注意这里的T 是大写的，温度单位只能是K。
4、传热过程传热过程专指热量从固体壁面一侧流体通过固体壁面传递到另一侧流体的过程。
它包括三个环节，分别属于对流换热、热传导、对流换热。
传热过程的总热阻R即为三个子环节的子热阻串联相加。
即为：在计算传热过程的热流量或热流密度时可以直接使用“动力/热阻”来计算。
我们只需记住上述这些面积热阻就行，面积热阻更为常用。
5、稳态过程在计算中如果已知传热过程为稳态过程，那么要知道这意味着传热过程的三个子环节的热流密度均相等。
如果不相等，某个环节的热流密度大，那么该处温度会逐渐增加，即为非稳态过程。
6、温度梯度：沿等温线法线方向的温度变化率，该方向的变化率最大。
温度梯度为矢量，其方向为沿等温线法线方向指向温度升高的方向。
可是用gradt来表示。
7、热扩散率a热扩散率a越大，温度变化传播越迅速，物体的温度能更快的趋于一致。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
计算传热学课程设计摘要：利用数值模拟技术构建应拟实脸是提高传热学教学及实脸效率、降低实脸成本的一条有效途径。
本学期基于Mtalab强大的数位计算功能、圈形处理能力,学习了传热学应用模拟实。
本文主要介绍了在本学期，老师在课堂上所讲的一些基本内容和我在本学期主要学习到了那些东西和课下作业在计算机上进行的数值模拟。
关键字：计算传热学Matlab 数值模拟一、介绍计算传热学又称数值传热学：是指对描写流动与传热问题的控制方程采用数值解法通过计算机予以求解的一门学科。
数值解法是一种离散近似的计算方法。
它所能获得的解不像分析解那样是被研究区域中未知的量的连续函数，而只是某些代表性的点（称为节点）上的近似解【1】。
MATLAB是在上个世纪80年代,美国几个主要从事数学研究的学者,针对矩阵运算的复杂性,开发的专门用于矩阵运算的程序,结果应用效果很好【2】。
二、内容本学期只要学习了五讲内容。
第一讲是微分方程近似解的概述，主要内容是一维翅片传热问题及其解析解、求连续的近似解析函数的方法包括了RITZ法，权余法，配置法，以及迦辽金法。
第二讲是导热的有限体积法，主要内容有一维导热问题的数值解和一维非稳态问题。
要求一维导热问题的数值解首先要将求解区域离散化，离散的方法有外节点法和内节点法；其次在把源项线性化控制方程的离散化；再次进行交界面参数的计算和跃界面的处理，交界面参数计算的方法有线性插值法、调和平均法、待求变量插值和Kirchhoff变换法；最后进行边界条件的处理和差分方程的求解。
一维非稳态问题主要需要进行非稳态项的处理、控制方程的离散化、边界条件的处理和求解。
第三讲是离散方程的求解，主要内容是一维离散方程的求解和对位离散方程的求解。
一维离散方程求解的方法有TDMA、Jacobi迭代、Gaus-Seidel迭代和SOR迭代；多维离散方程的求解方法有点迭代、线迭代和ADI。
第四讲对流-扩散的有限体积法，主要内容是稳定的一维对流-扩散、非稳定的一维对流扩散和多维的对流扩散。
解决稳定的和非稳定的一维对流-扩散问题主要运用的格式有中心差分格式、一阶逆风格式、混合格式、乘方格式和指数格式。
这几种格式各有优缺点，在运用是要根据具体情况选用，以使模拟的精度更高，更接近具体的数据，以便更好的解决在工程运用中的具体问题。
第五讲是对流-扩散的有限体积法,只要内容是在假设流畅已知的条件下，解决稳定的一维对流-扩散、非稳定的一维对流扩散和多维的对流扩散问题。
他们是通过加错网格的压力修正方法和同位网格的压力修正方法。
第六讲是流场计算，这一讲是通过有限体积计算不可压缩流体的流场。
求解不可压缩流动流场数值解的方法有联立求解个变量代数方程组的方法和分离式求解代数方程组的方法。
而流场求解的关键问题是压力梯度的离散和压力的求解。
三、课后作业地源热泵柱热源模型，孔壁半径65mm，岩土导热系数1.8，体积比热3x106，孔壁Q=50W/m无穷远处温度18 给出导热微分方程，用一维数值模型求岩土径向温度数值解与解析解比较这个问题我们可以通过matlab得出各点温度分布图形。
在本题中节点的划分我使用的是外节点法，采用元体平衡法进行方程分析，分别对节点1、2、…、n-1、n进行分析，此处为计算方便直接将节点N的温度等于无穷远处温度。
通过matlab得出的图形如下所示：隐式格式显示格式通过matlab模拟的两个图形的对比可知显式格式和隐式格式的图形对比可知在两个温度模拟的差别并不大，说明两种格式都能反映岩土径向温度的变化情况。
参考文献「1〕陶文栓.传热学[M〕.北京:高等教育出版社,2005.「2〕彭芳辟.数学物理方程的Matlab解法与可视化〔M」.北京:清华大学出版社,2004.附录非稳态隐式clear;qb=50;k=1.8;ri=0.065;T8=18;ro=2.065;rc=3000000;Sc=0;Sp=0;NumNode=100;dx=(ro-ri)/(NumNode-1);Xo=ro;Xi=ri;X=[Xi:dx:Xo];tstep=1000;dtmax=rc*dx*dx/2/k;dt=60;Ar=ri;AWi(1)=0;AEi(1)=(k*Ar*2)/dx+k;APi0(1)=rc*((Ar+dx/2)*(Ar+dx/2)-Ar*Ar)/dt; %°ë¸ö¿ØÖÆÌåAPi(1)=APi0(1)+AWi(1)+AEi(1)-Sp*((Ar+dx/2)^2-Ar^2)*dx; %ÒþʽBi0(1)=((Ar+dx/2)*(Ar+dx/2)-Ar*Ar)*Sc+qb/3.14;%½Úµã2 ---- NumNode-1·½³ÌϵÊý ´Ë²¿·Ö²»Ëæʱ¼ä±ä»¯for i=2:NumNode-1Ar=ri+(i-1)*dx; %iµãÃæ»ýÒò×Ó% Ar+dx/2; %e½çÃæÃæ»ýÒò×Ó% Ar-dx/2; %w½çÃæÃæ»ýÒò×ÓAWi(i)=(2*k*Ar)/dx-k;AEi(i)=(2*k*Ar)/dx+k;APi0(i)=(rc*((Ar+dx/2)*(Ar+dx/2)-(Ar-dx/2)*(Ar-dx/2)))/dt;%AP(i)=AP0(i); %ÏÔʽAPi(i)=APi0(i)+AWi(i)+AEi(i)-Sp*Ar*dx; %ÒþʽBi0(i)=((Ar+dx/2)*(Ar+dx/2)-(Ar-dx/2)*(Ar-dx/2))*Sc;%BP(i)=B0(i)+(AP0(i)-AW(i)-AE(i)+dx*dx*Sp)*T0(i);endAr=ri+(NumNode-1)*dx;AWi(NumNode)=0;AEi(NumNode)=0;APi0(NumNode)=1;%AP(NumNode)=AP0(NumNode);APi(NumNode)=APi0(NumNode);Bi0(NumNode)=0;%B0=B0';TT0(1:NumNode)=18; %ÈÎÒâ³õʼֵfor t=1:tstepBB(1,1)=Bi0(1)+APi0(1)*TT0(1); %BBÿһ²½±ä»¯AA(1,1:2)=[-APi(1) AEi(1)]; %AE*Tw-AP*Tp+AE*TE=-BP %µÚ2~NumNode-1½Úµãfor i=2:NumNode-1AA(i,i-1:i+1)=[AWi(i) -APi(i) AEi(i)]; %AE*Tw-AP*Tp+AE*TE=-BP BB(i,1)=Bi0(i)+APi0(i)*TT0(i);end%AE*Tw-AP*Tp+AE*TE=-BPAA(NumNode,NumNode-1:NumNode)=[AWi(NumNode) -APi(NumNode)]; BB(NumNode,1)=Bi0(NumNode)+APi0(NumNode)*18;BB=-BB; %AE*Tw-AP*Tp+AE*TE=-BPT2=inv(AA)*BB; %Ö±½Ó½â AA*T=BB £¿£¿£¿½á¹û²»¶Ô£¿%A3=[AW' -APi' AE']%T2=tdma1(A3,BB);TT0=T2'; %¸üгõÖµTT(t,:)=T2'; %±£´æµ±Ç°Öµ Òþʽ½á¹ûendfigure;plot(X,TT(500,:),'-x');figure;plot([1:500]*dt/60,TT([1:500],1),'k-');非稳态显式clear;qb=50;k=1.8;ri=0.065;T8=18;ro=1.065;rc=3000000;NumNode=101;dx=(ro-ri)/(NumNode-1);Xo=ro;Xi=ri;X=[Xi:dx:Xo];tstep=1000;dtmax=rc*dx*dx/2/k;dt=60;Ar=ri+dx/2;AW(1)=0;AE(1)=(2*3.14*k*Ar)/dx;AP0(1)=(rc*(Ar*Ar-ri*ri)*3.14)/dt; AP(1)=AP0(1);B0(1)=qb;for i=2:NumNode-1Ar=ri+(i-1)*dx;AW(i)=(2*3.14*k*(Ar-dx/2))/dx;AE(i)=(2*3.14*k*(Ar+dx/2))/dx;AP0(i)=(rc*2*3.14*Ar*dx)/dt;AP(i)=AP0(i);B0(i)=0;endAr=ri+(NumNode-1)*dx;AW(NumNode)=(2*3.14*k*(Ar-dx/2))/dx;AE(NumNode)=(2*3.14*k*(Ar+dx/2))/dx;AP0(NumNode)=(rc*2*3.14*Ar*dx)/dt;AP(NumNode)=AP0(NumNode);B0(NumNode)=0;B0=B0';T0(1:NumNode)=18;for t=1:tstepAr=ri+dx/2;yy=AP0(1)-AW(1)-AE(1);B(1,1)=B0(1)+yy*T0(1);T1(1)=(AE(1)*T0(1+1)+B(1))/AP(1); for i=2:NumNode-1Ar=ri+(i-1)*dx;yy=AP0(i)-AW(i)-AE(i);B(i,1)=B0(i)+yy*T0(i);T1(i)=(AW(i)*T0(i-1)+AE(i)*T0(i+1)+B(i))/AP(i);endAr=ri+(NumNode-1)*dx;yy=AP0(NumNode)-AW(NumNode)-AE(NumNode);B(NumNode,1)=B0(NumNode)+yy*T0(NumNode);T1(NumNode)=(AW(NumNode)*T0(NumNode-1)+AE(NumNode)*T8+B(NumNode))/AP(NumNode); T0=T1;T(t,:)=T1;endfigure;hold onplot(X,T(1,:),'k-');plot(X,T(300,:),'m--*');plot(X,T(600,:),'r:p');plot(X,T(900,:),'g-v');figure;plot([1:500]*dt/60,T([1:500],1),'k-');。