导热微分方程剖析.pptVIP

第三讲 导热微分方程 前言 运用物理数学方法研究表示特定现象各物理量的关系。确定物理过程的量随空间时间的变化。 在铸造凝固过程中我们要得到温度随时间空间位置变化的情况。 数学物理方法： 限制在一定时间间隔内（微元时间间隔dt ），并在整个空间内只研究一个微元体（选取微元体的体积为 dv）。 从数学观点看，这些量是无穷小；物理观点看，这些量足够大，以致在所研究的范围内能忽略介质的不连续性，而作为连续介质处理。 1. 直角坐标系下的导热微分方程 对所选取的微元体应用能量守恒定律： 在 dt 时间内，由于导热从外部进入微元体的热量以及微元体内热源所产生的热量等于微元体所包含的内能或焓的变化。 焓（enthalpy）：H，系统热力学参数。定义：H=U+PV。是状态函数，即系统的状态一定，焓就是定值了。 1. 直角坐标系下的导热微分方程 按照能量守恒定律，微元体的热平衡式可以表示为下列形式： 导入微元体的总热流量+微元体内热源的生成热=微元体内能的增量+导出微元体的总热流量 导入微元体的总热量： 由Fourier导热定律，在 dt 时间内 x 方向通过x=x表面流入微元体的热流量为: 微元体内热源的生成热： 若单位时间单位体积物体中内热源的发热率为 则在dt 时间内，微元体的发热量为： 导出微元体的总热量： 在dt 时间内，x方向通过x = x+dx 表面流出微元体的热流量dQx+dx，可按Taylor级数展开如下： 微元体内能的增量： 微元体内能的增量: 微元体的热平衡式: 化简后的热平衡式： 2. 其他坐标系的导热微分方程 在工程实际中往往涉及柱面或球面对称的导热问题，边界条件给定在一个表面上，此表面具有坐标保持不变的性质。 通过坐标变换，可以将直角坐标下的导热微分方程变换成相应坐标系的导热微分方程。 圆柱坐标系： 设圆柱坐标系下任一点P（r，θ，z），该点在直角坐标系中的投影关系得出直角坐标与柱坐标的变换关系如下： x= r cosθ r=(x2+y2)1/2 y= r sinθ θ= tan-1(y/x) z= z 柱坐标系下导热微分方程的一般形式： 3.直角坐标系下一般方程的特殊形式： 原因： 由于导热微分方程的一般形式存在非线性，往往难以直接进行数值求解。但在某些特殊条件下，可将上述问题简化为同问题所建立的条件相一致的最简单形式，以便进行数值求解。 3.直角坐标系下一般方程的特殊形式： 无内热源、常物性条件下导热微分方程的简化形式： 稳态导热： 一维： 二维： 三维： 3.直角坐标系下一般方程的特殊形式： 无内热源、常物性条件下导热微分方程的简化形式： 非稳态导热： 一维： 二维： 三维： 导温系数或热扩散系数（Thermal diffusivity）: 4.导热过程的定解条件 定解条件： 使微分方程得到特解的附加条件。 对于导热微分方程，通过数学方法原则上可以获得方程的通解。然而，对于实际工程问题而言，还要求得出既满足导热微分方程式，又满足根据具体问题规定的一些附加条件下的特解。 非稳态导热问题的定解条件： 初始条件： 初始时刻温度分布。 边界条件： 物体边界上的温度或换热情况。 稳态导热问题的定解条件： 没有初始条件，仅有边界条件。 边界条件： （1）第一类边界条件： 给定了边界上的温度值。 最简单的形式：给定边界温度保持不变，即 Tw = 常数。 对于非稳态导热问题，这类边界条件给定温度为边界上空间位置与时间的函数，即： 边界条件： （2）第二类边界条件： 给定边界上的热流密度值。 最简单的形式：给定边界上的热流密度保持定值，即 qw= 常数。当qw= 0，绝热边界条件。 对于非稳态导热问题，这类边界条件给定温度为边界上空间位置与时间的函数，即： 边界条件： （3）第三类边界条件： 给定边界上物体与周围流体间的表面换热系数 hc ，及周围流体的温度 Tf