非稳态、源项的离散和离散方程解法.pptxVIP

扩散方程离散方法一维稳态扩散源项：阶梯型假设；局部线性化扩散项：对梯度进行插值，用中心差分算法界面导热系数的确定方法1. 算术平均法假设扩散系数在P和E间线性变化界面导热系数的确定方法2. 调和平均法uch bette二维稳态扩散三维稳态扩散非稳态项的离散一维非稳态扩散UNSTEADY ONE-DIMENSIONAL CONDUCTION需要对温度随时间的变化进行离散时间上的离散显性格式显性格式假设T在整个时间步长内为旧的值，在时刻变成新的值。?所有系数应该为正。对于恒定的扩散系数k和相等的网格步长 δxPE = δxWP = ?x,缺点：限制了最大的时间步长。当减少空间步长时，要用更小的空间步长（二次方）。显示格式更容易计算，不需要解方程组。Crank-Nicolson格式Crank-Nicolson格式Tp随时间是线性变化的。Crank–Nicolson格式基于中心差分在时间上具有二阶精度。当时间步长够小，可以得到比显式更精确的结果。计算的总的准确度还取决于空间的离散格式。故Crank–Nicolson格式常与空间中心差分一起用。全隐性全隐形假设在t 时刻， Tp 突然变为新的值，并在整个时间步长内保持不变。因为全隐性在时间上只有一阶精度，为保证结果的准确性，应使用小时间步长。总结时间步长大，温度的下降常常一开始快速下降，然后平坦。全隐性比Crank-Nicolson的线性更符