一维稳态热传导方程的数值解法及其分享课件.pptVIP

§6 一维稳态热传导方程的数值解法及其应用;6.1 一维稳态导热的通用控制方程;对控制容积P做积分导出其方程的离散形式 对于源项S，常表示为温度的函数 S=SC+SPTP 式中： SC --- 常数， SP --- P点的斜率 TP --- P点的温度;对控制方程在控制容积P中进行积分得：;对扩散项T 随x 呈分段线性分布得： ;;系数aE，aW： 1）分别代表了节点P,E间及W,P间导热阻力的倒数（热导），其大小反映了节点E,W处的温度对P点温度的影响程度； 2）在aE，aW中出现了界面上当量导热系数λ e ，λw，在进行数值计算时物性数据及温度等变量都存放在节点上，必须找出如何由相邻两节点的导热系数来获得当量导热系数。(λ~t，如钢材与耐材结合时λ会有一突变，同一钢材加热时，铁素体向奥氏体转化);6.2 计算当量导热系数的两种方法;2）调和平均法 设在控制容积P，E的导热系数不等，则根据界面上热流密度连续的原则，由傅立叶定律有： ;由上两式有： ;6.3 两种方法的比较;若按调和平均法计算，由5式则导热热阻为 即温度TP将一直扩展到界面e，而温降TP- T E实际上发生在 内。 说明P，E间的导热热阻由导热系数小的决定，符合传热学原理。所以此种情况下，调和平均法符合。 对于表征输运特性的物性参数，如导热系数，动力黏度，调和平均法均优于算术平均法。;3）在导热系数突变时，使用该调和平均值不必采用极密的网格，且对有内热源且导热系数连续变化的场合也要比算术平均好得多。（把阶跃面作为控制容积的分界面） 4）把物性阶跃面设置成一个节点的位置比作为控制容积分界面，使计算结果会更加精确。（由于此种情况阶跃面两侧温度梯度不同，如按3处理，相当于用平均值来代替，采用此种方法处理时，物性阶跃面两侧温度梯度单独计算，提高了计算精度。） ;一维稳态导热方程的离散形式可表示成：;6.4 非线性问题的处理步骤; 其中 设初值为T*，迭代后新的温度分布为T， 下一次迭代部分改进 f-松弛系数，f=0~1为欠松弛，f〉1为超松弛。 f=0时重复计算，无更新， f=1时更新快，但易发散。 强烈非线性问题迭代需???用欠松弛方法促进收敛。 F=0.1时新值偶和进去很小，二次结果相差很小，误以为收敛，此时F调大些，如仍相差很小，则收敛。 ;6.5 源项处理方法;当S 与T相关时，需对其线性化。S可表示成 该未知量的线性函数。在控制容积P内可表示成 ;S形式可以不同， S c ，S p均可以是T的函数;1）当源项为未知量的函数时，线性化的处理比假定 源项为常数更合理。即; 为了保证代数方程迭代求解的收敛。Δν为控制容积的体积，线性代数方程迭代求解收敛的一个充分条件是对角占优，即 ;6.6 促进收敛的四项基本原则;3、源项的负斜率线性化 只有当 才能保证中心节点系数 4、相邻节点系数之和 因控制方程是微分方程，往往只包含有因变量Φ的导数项，若Φ满足，Φ+c也满足，微分方程这性质也反映在离散方程中。可以理解为中心节点值TP是各个相邻节点值Tnb的一个加权平均值。;课后问题：