4-1_导热问题的数值解法.pptVIP

4-3 非稳态导热问题的数值解法 (2) 热平衡法 三、建立节点物理量的代数方程（离散方程） 右边第三类边界 显示格式存在稳定性问题：如果节点 tn(i) 前面的系数小于零，则数值解出现不稳定的震荡结果。 显示格式 显示格式：格式右边全部为第 i 时间层的温度值，只要 i 时间层温度已知，即可计算得到 i+1 时间层的温度。 即：空间步长?x和时间步长?τ的选取有限制 显示格式的稳定性条件： 4-3 非稳态导热问题的数值解法 四、一维导热问题显式差分格式的不稳定性讨论 隐式格式 隐式格式：空间离散采用（i+1）时层的值。 隐式格式不存在稳定性问题，对时间步长和空间步长没有限制，但是计算量较大。 非稳态导热节点离散方程的两种格式： 4-3 非稳态导热问题的数值解法 2．讨论一维导热问题显式差分格式稳定性限制的物理意义 由节点显式差分方程可以看出， 前的系数可能是负值。当系数为负时， 时刻节点温度越高 ， 时刻的温度却越低，违背了热力学第二定律。 4-3 非稳态导热问题的数值解法 表面温度在不同 Δτ 时的计算值 4-3 非稳态导热问题的数值解法 4-3 非稳态导热问题的数值解法 1．首先选择空间坐标间隔 Δx ，即距离步长。对二维问题一般使 Δy=Δx　。 2．对显式格式差分方程，根据方程的稳定性条件选择允许的最大时间步长 Δτ ； 在稳定性