3.3 微分方程的有限容积法离散解析.pptVIP

3.3 微分方程的有限容积法离散 胡 茂 彬 /~humaobin/ humaobin@ 3.3.1 控制容积积分法离散 此方法基于对控制方程的积分 控制容积积分法实施步骤 积分：将守恒型控制微分方程在控制容积及其时间间隔内作空间、时间积分，能够解除某些微分号的先积出来。 2) 作型线假设：根据离散精度要求，选择函数及其导数对时间、空间的分布曲线(型线) 。 3) 进一步积分：按照型线对积分式各项作定积分，并整理成节点上的离散代数方程。 常用型线选择 – 空间 分段线性分布 常用型线选择 – 时间 分段线性分布 例：一维有源对流扩散方程的离散 1) 积分： 2) 型线假设 非稳态项： 需对空间做假设，取阶梯分布 对 流 项：需对时间分布取假设 a) 阶梯显式： c) 取随 t 分段线性分布算术平均格式 Crank-Nicolson 格式 扩散项：需对时间分布取假设 a) 阶梯显式： c) 取随 t 分段线性分布算术平均格式 Crank-Nicolson 格式 源项：需对时间、空间分布取假设 a) 阶梯显式： c) 取随 t 分段线性分布算术平均格式 Crank-Nicolson 格式 界面上的函数值和导数值 分段线性分布 界面上的函数值和导数值 离散格式：各种情况的组合 显式 全隐式 Crank-Nicolson 格式 均匀网格，不可压 整理出最终的离散表达式 关于型线假设的讨论 只有适当选择型线，才能导出离散方程；一旦离散方程建立，型线使命完成 ? 这与有限元方法不同 型线可变化多样：即同一方程不同物理量可选不同的型线；同一物理量不同坐标也可选不同的型线；同一物理量同一坐标在方程不同项中也可选不同型线 型线选择多样 ? 离散方程的多样性 3.3.2 控制容积平衡法离散 此法不需要控制方程，直接通过 物理量平衡条件得到离散格式 等式右端各项的时间层次型线假设可采取三种不同方案：阶梯显式、阶梯全隐式、C-N格式 3.3.3 控制容积法离散方程的四条基本规则 离散结果通用型式 四条基本规则 控制容积界面上的相容性——控制体任意界面上各类通量（质量流量、动量流量、能量流量）都只能有一个值，这个值由界面的控制体共同确定 离散方程求解函数各系数都是正值——只有这样，才能保证解的物理真实性 3. 源项线化时需取负斜率—只有这样，才能保证正系数原则（上一条）不致遭破坏 4.源项为常数的稳态问题满足邻近系数求和 * 阶梯分布 阶梯分布 b) 阶梯隐式： b) 阶梯全隐式： b) 阶梯全隐式： 空间取阶梯分布，而时间可取： 阶梯分布