计算流体力学第讲差分方法.pptVIP

* Copyright by Li Mingjun 算例1：有限差分法求解抛物型方程 一维非定常热传导方程 初始条件 t=0, T=T0(x) 边界条件 既可以采用显示法也可以采用隐式法。显示格式: (a) (b) (c1) (c2) 3. 数值算例 第二十九页，共四十七页，2022年，8月28日 * Copyright by Li Mingjun x=0 x=L 第三十页，共四十七页，2022年，8月28日 * Copyright by Li Mingjun 第三十一页，共四十七页，2022年，8月28日 * Copyright by Li Mingjun 表1 有限差分法计算结果(FDS)与解析解(AS)在x=0.3的对比数据(r=0.10) 第三十二页，共四十七页，2022年，8月28日 * Copyright by Li Mingjun 表2 有限差分法计算结果(FDS)与解析解(AS)在x=0.3的对比数据(r=0.50) 第三十三页，共四十七页，2022年，8月28日 计算流体力学第讲差分方法 * 第一页，共四十七页，2022年，8月28日 第0节 前言 (1) 有限差分法研究背景 (2) 研究CFD的手段 研究CFD的理论手段 研究CFD的实验手段 研究CFD的计算手段 * 第二页，共四十七页，2022年，8月28日 * 传统计算方法： 有限差分法， 有限体积法 ， 有限元法， 谱方法（谱元法）等； 最近发展的方法: 基于粒子的算法（格子-Boltzmann, BGK），无网格 优点 缺点 适用范围 有限差分法 简单成熟，可构造高精度格式 处理复杂网格不够灵活 相对简单外形的高精度计算 有限体积法 守恒性好，可处理复杂网格 不易提高精度（二阶以上方法复杂） 复杂外形的工程计算 有限元法 基于变分原理，守恒性好 对于复杂方程处理困难 多用于固体力学等 谱方法 精度高 外形、边界条件简单 简单外形的高精度计算 LBM方法 算法简单，可处理复杂外形 精度不易提高 复杂外形的工程计算 (1) 有限差分法研究背景 * Copyright by Li Mingjun 第三页，共四十七页，2022年，8月28日 * * 流体力学 理论研究 实验研究 数值研究 计算流体力学 （数值计算技术、 计算方法研究） 理论研究： 格式推导、 稳定性分析，精度及误差分析，…… 数值实验： 采用实际问题考核方法的正确性 数值研究： 采用数值计算推导格式、考察精度/稳定性/分辨率…… “计算流体力学”作为一个学科， 其研究手段依然包括理论、实验及数值模拟。 Copyright by Li Mingjun 第四页，共四十七页，2022年，8月28日 * 与 的依赖关系 * (2) 研究CFD的手段 例：Fourier分析 线性系统： 线性方程+ 线性格式 任意函数都可分解为三角函数的叠加 差分系统 (解差分方程) 初始值 数值解（特定时刻离散的函数值） 记为： 是差分算子，把离散函数（有限点列） 映射 为另一个离散函数 {vi} 与{ui}的依赖关系 线性系统，可大为简化 波数空间单一的依赖关系： 原理： 线性系统，输入一个波，只能输出一个波（且波数不变）。 非线性系统会产生多个谐波 线性差分系统： 针对一个单波， 研究经过差分系统后的变化就可以了解该系统。 Fourier误差分析； Fourier稳定性分析 理论分析的局限性： 对于复杂系统（非线性方程、非线性格式）非常困难 Copyright by Li Mingjun 研究CFD的理论手段 第五页，共四十七页，2022年，8月28日 * * 研究CFD的实验手段: 例： 精度分析 思想： 通过具体算例来研究（考核，分析…）差分方法 典型的文章： 提出方法+理论分析 + 算例验证 差分离散 理论方法，Taylor展开，求余项 。 对于复杂（如非线性）格式，难度大。 实验方法，通过算例考核精度 精确解: 为该离散函数的模 计算误差: 分析误差对网格步长的依赖关系 斜率为精度的阶数 （通常用最小二乘法计算） 斜率为精度的阶数n Copyright by Li Mingjun 常用的模