哈尔滨工业大学 计算传热学 第一章 绪论-2013.ppt

计算传热学Numerical Heat TransferComputational Heat Transfer（32学时） 教材及参考书 考试方式 背景 1. 大多数传热问题得不到解析解； 2. 实际工程问题需要的时间和空间有限个 点上的值，并容许有一定的误差。 因此希望通过某些近似方法来获得物理问题在时间和空间有限个点上的近似解。 计算传热问题的基本步骤 主要方法 教学目的 第一节 描写传热与流动问题的 控制方程 导热 对流 对流 辐射传递方程（RTE） 第二节 控制方程的守恒与非守恒形式及单值条件 方程的守恒性 初始条件与边界条件 第三节 控制方程的数学分类及其 对解的影响 Steady state temperature distributed of an insulated rod Transient distribution of temperature in an insulated rod 第四节 计算传热学求解问题的基本思想和总体步骤 Numerical Heat Transfer, Part A: Application Numerical Heat Transfer, Part B: Fundamentals Published by Taylor Francis Ltd 计算传热学的专门期刊 手册 W. J. Minkowycz, E. M. Sparrow, Handbook of Numerical Heat Transfer, 2nd Edition, Wiley, 2006 Mass, momentum and energy equations 1.??连续方程 2.??动量方程 3.??能量方程 状态方程： 组分守恒方程 第一项：物理量随时间的变化率： 非稳态项 第二项：因流体运动携带引起的变化率： 对流项 第三项：因梯度引起扩散的作用 扩散项 第四项： 源项 方程各项的物理意义： 导热、对流、扩散控制方程的通用形式 通用控制方程中，对流项采用散度的形式表示，在数值计算中称该方程为守恒型方程。 这样做可以保证通过某一公共界面的”流密度”不会因所考虑的单元或控制体不同而不同。 要使方程定解，必须给定： （1） 给定 值：第一类边界条件 （2） 给定 ：第二类边界条件 （3） 给定 与 的关系：第三类边界条件 1.? 几何条件 2.? 物理条件： …… 3.? 边界条件 初始条件的确定，因问题的区别可有很大不同。 （1）??非稳态问题，求某一时刻解，要给出准确的初始条件。 （2）?非稳态问题，求发展到下一稳态的解，初始条件有随意性和收敛性，主要考虑计算的稳定性和收敛性。 （3） 稳态问题，用非稳态方法求解，目的是改变方程属性，提高稳定性，收敛性。 4. 初始条件 例子: 突扩区域内的流动和传热 温度、速度边界条件 计算域所有的边界上都要给出吗??? 出口边界上怎么给!!!!!........? 第一项：物理量随时间的变化率： 非稳态项 第二项：因流体运动携带引起的变化率： 对流项 第三项：因梯度引起扩散的作用： 扩散项 第四项： 源项 方程各项的物理意义： 导热、对流、扩散控制方程的通用形式 二阶二元偏微分方程 Elliptic Parabolic hyperbolic 影响区和依赖区 椭圆型方程 边值问题 抛物型方程 步进法(marching method) 初值问题 如:非稳态导热、边界层型流动与传热问题 双曲型方程 Domains of dependence for the (a) hyperbolic,(b) parabolic and (c) elliptic problem 坐标的方向性（Spalding, D. B.） 单通道坐标： 坐标中给定位置的状态只受该位置一边的状 态变化影响。 双通道坐标： 受两边状态变化的影响。 抛物型： 至少有一个单通道坐标。 椭圆、双曲型：没有单通道坐标。 例：二维非稳态导热，从空间坐标看是双通道，但从时间坐标看是单通道的，故为抛物型问题。 应当注意，当一个方程是抛物型时，计算存储量和计算时间都将大大减少。因此在很多算法中，都通过一定的方法，将方程变换成抛物型。 One-way coordinate Two-way coordinate 下面的框图来表示。 基本思想 、 间的关系 * 上课时间：10-17周，周三7-8节，周四5-6节 上课地点