导热微分方程ppt课件.ppt

《导热微分方程ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第三讲导热微分方程前言运用物理数学方法研究表示特定现象各物理量的关系。确定物理过程的量随空间时间的变化。在铸造凝固过程中我们要得到温度随时间空间位置变化的情况。数学物理方法：限制在一定时间间隔内（微元时间间隔dt），并在整个空间内只研究一个微元体（选取微元体的体积为dv）。从数学观点看，这些量是无穷小；物理观点看，这些量足够大，以致在所研究的范围内能忽略介质的不连续性，而作为连续介质处理。1.直角坐标系下的导热微分方程对所选取的微元体应用能量守恒定律：在dt时间内，由于导热从外部进入微元体的热量以

2、及微元体内热源所产生的热量等于微元体所包含的内能或焓的变化。焓（enthalpy）：H，系统热力学参数。定义：H=U+PV。是状态函数，即系统的状态一定，焓就是定值了。1.直角坐标系下的导热微分方程按照能量守恒定律，微元体的热平衡式可以表示为下列形式：导入微元体的总热流量+微元体内热源的生成热=微元体内能的增量+导出微元体的总热流量导入微元体的总热量：由Fourier导热定律，在dt时间内x方向通过x=x表面流入微元体的热流量为:同理：微元体内热源的生成热：若单位时间单位体积物体中内热源的发热率为

3、则在dt时间内，微元体的发热量为：导出微元体的总热量：在dt时间内，x方向通过x=x+dx表面流出微元体的热流量dQx+dx，可按Taylor级数展开如下：同理：微元体内能的增量：微元体内能的增量:微元体的热平衡式:+化简后的热平衡式：为拉普拉斯(laplance)算子2.其他坐标系的导热微分方程在工程实际中往往涉及柱面或球面对称的导热问题，边界条件给定在一个表面上，此表面具有坐标保持不变的性质。通过坐标变换，可以将直角坐标下的导热微分方程变换成相应坐标系的导热微分方程。圆柱坐标系：设圆柱坐标系下

4、任一点P（r，θ，z），该点在直角坐标系中的投影关系得出直角坐标与柱坐标的变换关系如下：x=rcosθr=(x2+y2)1/2y=rsinθθ=tan-1(y/x)z=z柱坐标系下导热微分方程的一般形式：球坐标系：3.直角坐标系下一般方程的特殊形式：原因：由于导热微分方程的一般形式存在非线性，往往难以直接进行数值求解。但在某些特殊条件下，可将上述问题简化为同问题所建立的条件相一致的最简单形式，以便进行数值求解。3.直角坐标系下一般方程的特殊形式：无内热源、常物性条件下导热微分方程的简化形式：稳态导

5、热：一维：二维：三维：3.直角坐标系下一般方程的特殊形式：无内热源、常物性条件下导热微分方程的简化形式：非稳态导热：一维：二维：三维：导温系数或热扩散系数（Thermaldiffusivity）:说明:（1）表征温度变化的速度，即表征物体内部温度趋于一致的能力。是物体热惯性的度量。（2）与物质的性质有关。例如：液体和气体具有较大的热惯性，它们的导温系数就小。金属的热扩散系数比型砂大几十倍，铸件在金属型中要比砂型中冷却的快。4.导热过程的定解条件定解条件：使微分方程得到特解的附加条件。对于导热微分方

6、程，通过数学方法原则上可以获得方程的通解。然而，对于实际工程问题而言，还要求得出既满足导热微分方程式，又满足根据具体问题规定的一些附加条件下的特解。非稳态导热问题的定解条件：初始条件：初始时刻温度分布。边界条件：物体边界上的温度或换热情况。稳态导热问题的定解条件：没有初始条件，仅有边界条件。边界条件：（1）第一类边界条件：给定了边界上的温度值。最简单的形式：给定边界温度保持不变，即Tw=常数。对于非稳态导热问题，这类边界条件给定温度为边界上空间位置与时间的函数，即：边界条件：（2）第二类边界条件：

7、给定边界上的热流密度值。最简单的形式：给定边界上的热流密度保持定值，即qw=常数。当qw=0，绝热边界条件。对于非稳态导热问题，这类边界条件给定温度为边界上空间位置与时间的函数，即：边界条件：（3）第三类边界条件：给定边界上物体与周围流体间的表面换热系数hc，及周围流体的温度Tf。此类边界条件可表示为：边界条件：（4）辐射换热边界条件：当考虑物体表面的热辐射时，辐射换热边界条件：Tf为已知环境温度Tw为物体表面温度边界条件：线性化处理，使得该类边界条件具有与对流换热边界条件类似的形式，即hr称为辐

8、射换热系数边界条件：在实际导热问题中，物体表面经常同时存在着对流和辐射两种换热，其边界条件为：边界条件：（5）金属/铸型界面换热条件：λ1：为铸件材料的导热系数w1，w2：分别代表铸件与铸型表面Tw1，Tw2：分别代表铸件和铸型的表面温度hi：为金属/铸型界面换热系数，可以通过经验数值或假设一定的分布函数形式处理最简单的处理方式为：假设金属/铸型界面处于理想接触状态，此时其边界条件的表达式为：同样适合于不同铸型材料间接触界面换热条件的处理。初始条件：非稳态问题的初始条件：（1）为t