4-1 导热问题的数值解法(课件1).pptVIP

表面温度在不同 Δτ 时的计算值 4-3 非稳态导热问题的数值解法 4-3 非稳态导热问题的数值解法 1．首先选择空间坐标间隔 Δx ，即距离步长。对二维问题一般使 Δy=Δx　。 2．对显式格式差分方程，根据方程的稳定性条件选择允许的最大时间步长 Δτ ； 在稳定性条件允许范围内，Δτ 越大，计算工作量越小，但精度较差；对一维问题，一般取　　　　 ，即可满足工程精度要求；　　　 对于隐式差分方程，Δx、Δτ 可任意选取，不必进行稳定性条件校核； 4-3 非稳态导热问题的数值解法 五、节点方程组的求解 3. 按题意给定的初始温度分布，确定各节点上的温度初值； 4. 根椐建立的差分方程组，求 Δτ 时刻各节点的温度 tn(1) ； 5. 再由 tn(1) 为初值，换用τ=2Δτ 　　　　（即 k=2 ），重复 4 计算出 tn(2) ，如此反复，最后得到 j 时刻的 tn(j) 　。 4-3 非稳态导热问题的数值解法 * * 区域离散化：用一系列与坐标轴平行的网格线把求解区域划分成若干个子区域，以网格线的交点作为需要确定温度值的空间位置，称为节点 ( 结点 ) ，节点的位置用该节点在两个方向上的标号m，n表示。 以节点为中心，在两个节点的中心处划分界限，定出节点所代表的小区域，叫控制体。二维控制体是一个个的矩形面积。网格的步长在每一个方向上可以均匀划分