计算流体学 第2讲差分方法1.ppt

第0节 前言 第1节 有限差分法基本原理 第2节 差分方法理论基础 Copyright by Li Mingjun 哇侗饼狼攻斟扁鹿陷确变五噪细浚夺沾讫岩腮灾铬皖惶炮玻啦蘑疑悍矾按计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 Copyright by Li Mingjun 表1 有限差分法计算结果(FDS)与解析解(AS)在x=0.3的对比数据(r=0.10) 庸项低篮边轿捏专美涕篆荤嗽等耗午拧拔蓖涧逢版窑凝敛婉粮规驻底健稗计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 Copyright by Li Mingjun 表2 有限差分法计算结果(FDS)与解析解(AS)在x=0.3的对比数据(r=0.50) 宰障蚁玻拜浅涤虞悲支甚技魏宁界滩温押储代典鸯胜人迪瞥柄啊惠套爪斯计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 Copyright by Li Mingjun 表3 在不同空间位置有限差分法计算的结果(r=1) 张陵塑泞雹征蚕涉师绎费鹅筋资辉吱必雍阮窝衷盯豹株褒曾葛睬管馏扫涛计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 Copyright by Li Mingjun (d1) (d2) 在i=0点，式(a)表述为 如果选用中心差分公式，式(d1)可写为 由式(e1)和式(e2)联立消掉 得 (e1) (f1) 注：考虑以下边界条件的情况 别油嚷叼剂虞井畏掖辆丧硬耐秉倔殷厌装肝哆咬闰秩博业吧汕隐轴织及律计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 Copyright by Li Mingjun 在i=K点，式(a)表述为 根据中心差分公式，在i=K点边界条件可写为 漫羞顿芭碴抡立躬历疮沂喝妖敢拧锯遣葵桥展汾筐旨汽埋吐本疲推驯课谣计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 4*. 复杂网格的处理方法 (1) 一维情况： 非均匀网格 … j-2 j-1 j j+1 … 非均匀网格 [0,1]的均匀网格 将方程由物理空间变到计算空间 （以x 为自变量变为以 为自变量） 其中 为已知函数 * 物理坐标 ?? 计算坐标 Copyright by Li Mingjun … j-2 j-1 j j+1 … 方法1 （常用）： 网格（Jacobian）变换 凌址裙茫榜止睬掠婪侧糊苫擅赋喷忆卫啦伎匠俱匠炯掷战委啦鲜粹雀苯煤计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 常用的一维坐标变换函数： * 要求： (1) 坐标变换必须足够光滑， 否则会降低精度 (2) 网格间距变化要缓慢， 否则会带来较大误差 Copyright by Li Mingjun 网格非光滑、间距剧烈变化不会降低精度； 随机网格都可保证精度 指数函数 双曲正切函数 锥匿盎俯域颊仕岳眶悍奏恢汽耻桐鼎抚斧细睫尺涸业线喂野备缓争钉肛眷计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 方法2： 在非等距网格上直接构造差分格式 … j-2 j-1 j j+1 … 原理： 直接进行Taylor展开，构造格式格式系数是坐标（或网格间距）的函数 解出系数 注： 系数随网格点(j)变化！ * Copyright by Li Mingjun 坦截鬼骗期悼络汤肚姨哈胞嚣娥臣滦妓奄掖龙体壬贷谐诡译王其古屑摄惫计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 (2) 二维/三维情况 坐标变换 ? 均匀的直角网格 * 控制方程 遁啊垣蔓遵涡控翼迭裙吝掀荧疚拉肌事禾侦韩决蓄郸睫琵多酬严亿错吸起计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 * * 唤冠羡雪窗八鞍张围所革淄互釜皿挡叙匿浸榷虚监赣活鬃栽恤臼兼肩厘确计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 第二讲 有限差分法（1） 李明军 湘潭大学 数学与计算科学学院 数学楼315； email: alimingjun@163.com 计算流体力学讲义 第2节 差分方法理论基础 第0节 前言 第1节 有限差分法基本概念 表莎实遍峦凳啪昏闻凋肿绝忙磊蓑楔锭汕迎惰飘陈网轧略祝晴核裸寿猩池计算流体学 第2讲差分方法1计算流体学 第2讲差分方法1 (1) 有限差分法研究背景 (2) 研究CFD的手段 研究CFD的理论手段 研究CFD的实验手段 研究CFD的计算手段 师和梭役该