热流问题的数值计算.pptxVIP

热流问题的数值计算;4.1 一维导热问题;4.1.2 用控制容积积分法导出通用方程的 离散形式 假定源项S在任一控制容积内中的值可以表 示温度的线性函数： （2） Sc为常数，Sp为S = f(T) 的曲线在P点的斜 率，且规定恒为负值。;对方程（1）在控制容积内作积分，得 经整理，得 (3) ;令 则式（3）变为： 系数aE, aW分别代表了节点P, E之间以及W, P之 间导热阻力的倒数，反映了节点E，W处的温度对P点温度的影响程度，具有影响系数的意义。;4.1.3 界面上当量导热系数的确定方法;2.调和平均（harmonic mean） 利用界面上热流连续的原理，由Fourier 定律： 又 ; 因此 3.两种方法的比较 若 ，按算术平均法，当网格均分时， 表明P，E间的导热热阻主要由导热系数大的物体所决定，不符合传热学的基本原理。;按调和平均法： P, E间的热阻与调和平均方法的结果一致。? 调和平均法在多数场合都优于算术平均法。? 由此，系数aE, aW可以表示为： ;4.导热系数发生阶跃性变化时，阶跃面的两种处理方式 把物性阶跃面作为控制容积的分界面； 把物性的阶跃面设置成一个节点的位置，计算阶跃面上的热流密度时，精度高于第一种方法。;4.1.4 一维非稳态导热方程及其离散化 ;常用的T随 时间变化的型线可以表示为： F为0与1之间的加权因子。将T表达式代 入（2），得： ;化简整理后，得 其中 ; 取f = 0, 1,1/2可依次得显式，隐式和C-N格式。在直角坐标系中 当网格划分均匀时，对无内热源、 常物性的 导热问题三种格式可表示为： 显式： 隐式： ; C-N格式： 通过Von Neumann分析法可证明，对源项不随时间而变的问题，当 ???，格式（5）是绝对稳定的；当时，稳定的条件为 ;4.1.5 稳定的格式导致有物理意义的解;由式（5），可得 为网格Fourier数。 ; 从物理概念上，TP与前一时层的温度之比应永远为正值，并随时间的增长而趋近于0。在不同的f下， 随 的变化如图所示：图4-3。 只有全隐格式才能满足要求，任何f1的格式， 当 大于一定值后都会出现物理 上不真实的解。 ;4.3 多维非稳态导热方程的全隐格式 ;直角坐标的网格系统;非稳态项的积分： 扩散项： ;源项： 整理，得: ;其中 界面上的当量导热系数按调和平均法计算。 ;2.圆柱轴对称坐标系 ;其中;3、极坐标系;其中 都是相邻两节点间导热热阻的倒数， 具有热惯性意义， 热惯性越大，上一时层的温度对下一时层的影响越大。;4.3 源项及边界条件的处理 ;关于源项线性化的说明;因为离散方程都可表示为 线性代数方程迭代求解收敛的一个充分必 要条件是对角占优，即： 要求 ;4. 由代数方程迭代求解的公式 SP绝对值的大小影响到迭代过程中温度的变化速度，SP的绝对值越大，系统的惯性越大，相邻两次迭代之间TP的变化越小，收敛速度下降，但有利于克服迭代过程的发散；SP的绝对值小，可使变化率加快，但容易引起发散。;例1： ，取SC = 3，SP = -5。 例2：