导热微分方程及稳态导热

工学院机电工程教研室传输原理工学院机电工程教研室传输原理热量的传递方式特定条件下热量传递、温度分布的有关规律工学院机电工程教研室传输原理锂电池被加热至250摄氏度（华氏482度）后体积迅速膨胀直至爆炸，研究员将此过程称之为“热失控”。结果显示，锂电池内部结构在“热失控”前并不会发生变化，但在爆炸发生瞬间，电池内部铜质材料融化将温度提升至1000摄氏度，产生的热量从内部向外扩散导致“热失控”。锂电池通常被安装在设备内部并与其他电子元件相连，“热失控”不仅会使电池核心崩溃，还增加了设备内部短路和破坏附近物体的风险。3D热成像技术工学院机电工程教研室传输原理温度变化可使改变汽车颜色变换的热敏涂料诞生热敏反应涂料涂在尼桑SkylineR33汽车表面，当倾倒冷水时会使汽车变成橙色。工学院机电工程教研室传输原理通用导热能量方程工学院机电工程教研室传输原理基本概念一、温度场（Temperaturefield）温度场是指在各个时刻物体内各点温度分布的总称。由傅立叶定律知，物体的温度分布是坐标和时间的函数：,,,zyxft其中为空间坐标，为时间坐标。,,xyz工学院机电工程教研室传输原理基本概念温度场分类稳态温度场（定常温度场）:是指在稳态条件下物体各点的温度分布不随时间的改变而变化的温度场称稳态温度场非稳态温度场（非定常温度场）:是指在变动工作条件下，物体中各点的温度分布随时间而变化的温度场0t0t工学院机电工程教研室传输原理基本概念等温面与等温线等温面：同一时刻、温度场中所有温度相同的点连接起来所构成的面等温线：用一个平面与各等温面相交，在这个平面上得到一个等温线簇工学院机电工程教研室传输原理基本概念工学院机电工程教研室传输原理工学院机电工程教研室传输原理温度梯度沿等温面法线方向上的温度增量与法向距离比值的极限。温度梯度和热流密度nntnttnΔΔlimgrad0Δ•温度梯度是向量，垂直于等温面，正向朝着温度增加的方向；•温度梯度的方向是温度变化率最大的方向。△m△nmtnt工学院机电工程教研室传输原理kjiztytxttgrad温度梯度的解析定义：温度场中点处的温度梯度：),,(zyxft),,(zyxkjizyx－hamilton算子，经此演算，标量场变成了矢量场。tttijkxyz—nablagrad—gradient工学院机电工程教研室传输原理（4）热流密度热流密度是指单位时间经过单位面积所传递的热量，用q表示，单位为。2/Wm导热的热流密度与温度梯度成正比，即：nqnttgrad—导热系数，物性值。单位为W/(m·K)。负号是因为热流密度与温度梯度的方向相反。热流密度为矢量，其在x、y、z轴上的投影用傅立叶定律表示为：xtqxytqyztqz工学院机电工程教研室传输原理基本概念t+Δttt-Δt温度梯度和热流密度的方向都是在等温面的法线方向。由于热流是从高温处流向低温处，因而温度梯度和热流密度的方向正好相反。工学院机电工程教研室传输原理热流量是指单位时间内通过面积F所传递的热量，用Q表示，单位为W。dFFQq热流量矢量，它们和温度梯度位于等温面的同一法线上，且沿温度降低方向为正。总热量是指在时间内通过面积F所传递的热量，用Qτ表示，单位为J或kJ。00dddFFqQQ基本概念工学院机电工程教研室传输原理工学院机电工程教研室传输原理导热微分方程•确定导热体内的温度分布是导热理论的首要任务•傅里叶定律：•确定热流密度的大小，应知道物体内的温度场:2-grad[Wm]qt(,,,)tfxyz目的：确定导热体内部温度的分布,从而进一步用傅里叶定律计算换热量、计算热应力。(,,,)tfxyz工学院机电工程教研室传输原理导热微分方程一、导热微分方程的推导理论基础：傅里叶定律+能量守恒方程导热微分方程:根据能量守恒定律与傅立叶定律，建立导热物体中的温度场应满足的数学表达式。导热体内取一微元体，根据能量守恒定律，单位时间净导入微元体的热量加上微元体内热源生成的热量应等于微元体焓的增加量dv＝dv工学院机电工程教研室传输原理导热微分方程物理问题描述三维的非稳态导热体，且物体内有内热源（导热以外其它形式的热量，如化学反应能、电能等）。假设条件(1)所研究的物体是各向同性的连续介质；(2)热导率、比热容和密度均为已知；(3)内热源均匀分布，强度为[W/m3]；(4)导热体与外界没有功的交换。Φ工学院机电工程教研室传输原理导热微分方程1、导入与导出微元体的净热量d时间内、沿x轴方向、经x表面导入的热量：xxdqdydzdd时间内、沿x轴方向、经x+dx表面导出的热量：xdxxdxdqdydzddxxqqqxxdxx工学院机电工程教研室传输原理导热微分方程d时间内、沿x轴方向导入与导出微元体净热量xxxdxqdddxdydzdxd时间内、沿y轴方向导入与导出微元体净热量yyydyqdddxdydzdyd时间内、沿z轴方向导入与导出微元体净热量zzzdzqdddxdydzdz工学院机电工程教研室传输原理导热微分方程[导入与导出净热量]:()yxzdqqqddxdydzdxyz傅里叶定律：xtqxytqyztqz[()()()]dtttdxdydzdxxyyzz工学院机电工程教研室传输原理导热微分方程2、d时间微元体内热源的发热量vvqdxdydzd3、微元体在d时间内焓的增加量dxdydzdtcdv＝工学院机电工程教研室传输原理导热微分方程()()()vttttcqxxyyzz非稳态项源项扩散项dv＝这是笛卡尔坐标系中三维非稳态导热微分方程的一般表达式。其物理意义：反映了物体的温度随时间和空间的变化关系。工学院机电工程教研室传输原理傅里叶定律由此可见ɑ物理意义：①ɑ值大，即值大或c值小，说明物体的某一部分一旦获得热量，该热量能在整个物体中很快扩散，其内部各点温度扯平的能力越大。②ɑ越大，表示物体中温度变化传播的越快。所以，ɑ也是材料传播温度变化能力大小的指标，亦称导温系数。热扩散率表征物体被加热或冷却时，物体内各部分温度趋向于均匀一致的能力，所以ɑ反映导热过程动态特性，研究不稳态导热重要物理量。工学院机电工程教研室传输原理几种特殊情况①若物性参数λ，，ρ均为常数②无内热源，常物性：③稳态，常物性：④稳态，常物性，无内热源：222222vpqttttaxyzc简写为：2vppqtatcac为热扩散系数(导温系数)㎡/spc2tat2vq0t20t工学院机电工程教研室传输原理导热微分方程单值性条件：确定唯一解的附加补充说明条件，包括四项：几何、物理、初始、边界1、几何条件：说明导热体的几何形状和大小，如：平壁或圆筒壁；厚度、直径等2、物理条件：说明导热体的物理特征如：物性参数、c和的数值，是否随温度变化；有无内热源、大小和分布；工学院机电工程教研室传输原理导热微分方程3、初始条件：又称时间条件，反映导热系统的初始状态)0,,,(zyxft４、边界条件:反映导热系统在界面上的特征，也可理解为系统与外界环境之间的关系。t=f(y,z,τ)0x1x0x1x),,(zyfxt0x1x)thttx（工学院机电工程教研室传输原理导热微分方程导热微分方程＋单值性条件＋求解方法温度场2、导热微分方程的适用范围1）适用于q不很高，而作用时间长。同时傅立叶定律也适用该条件。2）若时间极短，而且热流密度极大时，则不适用。3）若属极底温度（-273℃）时的导热不适用。工学院机电工程教研室传输原理导热微分方程对于圆柱坐标系cqzttrrtrrtatv)11(2222222工学院机电工程教研室传输原理导热微分方程对于球坐标系cqtrtrrtrrratv]sin1)(sinsin1)(1[2222222工学院机电工程教研室传输原理导热微分方程课上作业：列出下列问题的数学描述：1.一块厚度为d的平板，两侧的温度分别为tw1和tw2。导热系数为常数2.已知一单层圆筒壁的内、外半径分别为r1、r2，导热系数为常量，无内热源，内、外壁面维持均匀恒定的温度tw1，tw2。rtw2r1r2tw1工学院机电工程教研室传输原理