3.2.3基于Taylor展开的有限差分离散试卷.pptVIP

3.2.3 基于Taylor展开的有限差分离散 胡 茂 彬 /~humaobin/ humaobin@ 1．基于 Taylor 展开，应用差分算子符号直接构成差分格式 例 3.1 具有均匀速度的一维对流方程 1) 时间向前，空间两步中心 (FTCS) 整 理 Courant 数： 对流方程离散的重要参数 2) 空间向前时间向前 (FTFS) 3) 时间向前、空间向后(FTBS) 例 3.2 常扩散系数的一维扩散方程 (1) 时间向前、空间二阶中心 (FTCS) 2) 时间两步中心、空间二阶中心(CTCS) 3) 时间向后、空间二阶中心(BTCS) 4) Crank-Nicolson格式 (CN) 2．基于Taylor展开构成差分格式的待定系数方法 例3.3 对匀流方程设计时间向前，空间二阶精度的迎风格式 迎风格式：从来流上游方向寻找依赖区 a0 时，应该选上游节点 j-2, j-1, j a0 时，上游节点是 j+2, j+1, j a 0 情况 代入之后 系数： 结果 a 0 的情况： 时间导数项 时间向前一阶差分 结果： a0 a0 3.2.4 其它构成导数有限差分离散的方法 1.多项式拟合求解函数构造函数导数的差分格式 用关于求解函数的 拟合多项式来逼近函数 A) 线性拟合 令 j 点坐标为 0 B) 二次拟合 令 j 点坐标为 0 多项式拟合方法在处理边界条件中的应用 例3.4 二维导热区域固壁边界问题 已知：热流 (定义输入壁面热流为正)，内部温度 未知：壁面温度 壁面附近温度二次拟合 壁面坐标为 0，向上为正方向 由热流确定壁面温度 另一个问题：由壁面温度求热流 三次多项式拟合温度 热流表达式 2 预估－校正多步法构成差分格式 第一步称为预估步 其后各步称为校正步 A 交替方向隐式离散格式(ADI-Alternative Direction Implicit) 例：2维无源非稳态导热方程 Peaceman-Rachford： 分两步进行，每步推进 1) 对 x 方向取隐式，对 y 方向取显式 2) 对 y 方向取隐式，对 x 方向取显式 1). 2). 3维的情况- 须做调整以使计算稳定 Brain Douglas 提出的 3D 恒稳格式 B MacCormack 格式这是一类格式，基于预估-校正思想 例：匀流方程 预估步：FTFS - 体现迎风思想 校正步：FTBS，并与原始值取平均 评述 该方法构成格式在预估步必须考虑迎风性 尽管每一步的精度都是一阶，两步合起来的最终精度为全二阶 该格式广泛应用于流动和传热问题的对流项离散 3. 分裂差分离散方法 与预估校正方法的区别是： 预估校正方法基于原方程离散 而分裂差分方法改造原方程再进行离散 例：一维线性对流扩散方程 二者等价 例：二维非定常线性热传导方程 4. 半离散化方法 适用于高超音速飞行等对时间导数离散精度要求非常高的问题 Motivation 前差、后差格式精度太低 中心差分格式，需要三个甚至更多个时层的离散方程，计算起步困难 方法：对时间导数项不离散，而只离散空间导数项，最终解常微分方程 例：一维匀流方程 和分裂 离散化 * 扩散数 精度每升高一阶，涉及的节点增加 1 个，展开项数也要增加 1 项。 j j-1 j-2 令： 加上原来值，取平均。 原始值 *