有限差分基础白.pptVIP

对于一维热流公式为 第三十页，共七十四页，2022年，8月28日 将对应节点（i, j）的微分方程写成如下形式 式中，θ为加权系数，取值范围为0≤θ≤1。 加权差分格式 加权差分格式为 第三十一页，共七十四页，2022年，8月28日 简化为 当θ=0时，显式格式。 当θ=1时，完全隐式格式。 当θ=1/2时，C-N格式。 当θ=2/3时，加辽金格式。 θ从0到1变化时，可得到不同的差分格式。 对于给定的Δt和Δ，随着θ的增长，计算精度下降，稳定性却越能得到保证。 隐式格式 第三十二页，共七十四页，2022年，8月28日 对于二、三类边界条件，在边界外设立虚节点，使边界节点变换为内节点： 用中心差商近似一阶微商； 边界节点取内节点差分方程。 2. 非稳态问题的边界条件 在开始进行计算的一瞬间（t=0），边界温度突然由T0变为Tw不大合理。因此，实际计算时，应作适当处理。 （1）给定温度Tw 第一步计算(t=0时) 边界节点温度为Tw/2； 完成第一步计算之后 固定温度边界节点保持Tw温度 第三十三页，共七十四页，2022年，8月28日 为提高整个差分格式的计算精度，常用中心差商来代替边界上的一阶微商。为此，在边界外设虚假节点。 （2）给定换热边界条件 在具有边长Δ为正方形网格的二维矩形区域中，有两类节点：一类是边上（如节点1），另一类在角上（