计算流体力学讲义第3讲-差分方法1.ppt

计算流体力学讲义 第三讲 有限差分法（1） 知识点： 差分方法的理论基础 （相容、收敛、稳定性；Lax等价定理；精度、修正方程; 守恒性） 差分格式的构造 差分格式的Fourier分析 传统计算方法： 有限差分法， 有限体积法 ， 有限元法， 谱方法（谱元法）等； 最近发展的方法: 基于粒子的算法（格子-Boltzmann, BGK），无网格 优点 缺点 适用范围 有限差分法 简单成熟，可构造高精度格式 处理复杂网格不够灵活 相对简单外形的高精度计算 有限体积法 守恒性好，可处理复杂网格 不易提高精度（二阶以上方法复杂） 复杂外形的工程计算 有限元法 基于变分原理，守恒性好 对于复杂方程处理困难 多用于固体力学等 谱方法 精度高 外形、边界条件简单 简单外形的高精度计算 粒子类方法 算法简单，可处理复杂外形 精度不易提高 复杂外形的工程计算 第三讲 有限差分法（1） （教材第3.1、3.2节及第4章） * Copyright by Li Xinliang Copyright by Li Xinliang * 流体力学 理论研究 实验研究 数值研究 计算流体力学 （数值计算技术、 计算方法研究） 理论研究： 格式推导、 稳定性分析，精度/误差分析，…… 实验研究： 数值实验， 采用实际问题考核方法的正确性 数值研究： 采用数值计算推导格式、考察精度/稳定性/分辨率…… “计算流体力学”作为一个学科， 其研究手段依然包括理论、实验及数值模拟。 与 的依赖关系 Copyright by Li Xinliang * 举例说明： 研究“计算流体力学” 学科的理论手段、实验手段及计算手段 研究CFD的理论手段 例：Fourier分析 线性系统： 线性方程+ 线性格式 任意函数都可分解为三角函数的叠加 差分系统 (解差分方程) 初始值 数值解（特定时刻离散的函数值） 记为： 是差分算子，把离散函数（有限点列） 映射为另一个离散函数 {vi} 与{ui}的依赖关系 线性系统，可大为简化 波数空间单一的依赖关系： 原理： 线性系统，输入一个波，只能输出一个波（且波数不变）。 非线性系统会产生多个谐波 线性差分系统： 针对一个单波， 研究经过差分系统后的变化就可以了解该系统。 Fourier误差分析； Fourier稳定性分析 理论分析的局限性： 对于复杂系统（非线性方程、非线性格式）非常困难 Copyright by Li Xinliang * 研究CFD的实验手段 例： 精度分析 思想： 通过具体算例来研究（考核，分析…）差分方法 典型的文章： 提出方法+理论分析 + 算例验证 差分离散 理论方法，Taylor展开，求余项 。 对于复杂（如非线性）格式，难度大。 实验方法，通过算例考核精度 精确解 为该离散函数的模 常用的模： 1 模： 2 模： 无穷模： 计算误差 分析误差对网格步长的依赖关系 斜率为精度的阶数 （通常用最小二乘法计算） 斜率为精度的阶数n Copyright by Li Xinliang * 常用的验证算例（“实验验证”） 考核方法通常找一些难度大的（条件苛刻、极端）的算例。否则，无法突出方法的优越性。 1维算例： Sod 激波管， Shu-Osher, 方波/尖波…… 2维算例： 前/后台阶、双马赫反射、二维Riemann问题、漩涡-激波干扰、翼型扰流、圆柱绕流 3维复杂算例： 各向同性湍流的DNS, 槽道湍流的DNS, 激波-边界层干扰的DNS Shu-Osher问题的计算结果 （Li et al. Init. J. Num. Fluid. 2005) 航空领域权威的考核算例—— DPW标准计算模型 常用一、二维算例整理后已发到流体中文网 Copyright by Li Xinliang * 研究CFD的计算手段 例： 差分格式构造 理论方法： 手工推导系数（工作量大） 数值方法： 通过数值手段推导系数 数值求解，获得系数 格式优化； 通过数值计算手段进行 Fourier分析; …… § 3.1 有限差分法基本原理 1. 差分方法的基本原理 离散点上利用Taylor展开，把微分转化成差分 … j-2 j-1 j j+1 … （等距网格） 多维问题，各方向独自离散；（时间同样考虑） 比有限体积法计算量小；