CFD基本理论.ppt

计计 算算 流流 体体 力力 学学ComputationalFluidDynamics 第二章 计算流体力学的基本理论 授课人：钱 昆 船舶工程学院趴族帅哮扩诌瘦各媒砚嘶佬记客腿玄弊帚指恤凤嘱蔷屯樟令拔糠酝禄司俘CFD基本理论CFD基本理论1.数值求解方法 以差分或积分方法，用代数方程代替微分方程，以代数运算代替微分运算，最终得到微分方程在离散点上的数值解。第二章 计算流体力学的基本理论 透阔流沈谋亭残迫甄膘揖蹈蜒响怂汕不邢伞靳汛晌兵蔽铸蚊罢秩窗壶叼阀CFD基本理论CFD基本理论1.1 典型数值求解方法 第二章 计算流体力学的基本理论 有限差分法有限体积法有限元方法 加权余量法有限分析法边界元方法 畸嫂搜侄您鸟枷轧述排铸籽蜂衔弦敷牲绍宿仙康渐仔陨仆入响凹蜘渣屯咙CFD基本理论CFD基本理论1.2微分方程分方程分类分类 第二章 计算流体力学的基本理论 当微分方程转化为差分方程并用数值方当微分方程转化为差分方程并用数值方法求解时，不同类型的微分方程，其数法求解时，不同类型的微分方程，其数值处理方法各异，其中包括提法的适定值处理方法各异，其中包括提法的适定性、物理解的性质、差分格式的适用性性、物理解的性质、差分格式的适用性等等扳张纯缆艰冗滔市汇窍膏匈时苞贾胀择一迅埃和鼠担康室绵言帕永往继茬CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 流体力学典型一阶拟线性微分方程组的分类流体力学典型一阶拟线性微分方程组的分类对于一阶拟线性微分方程组的向量形式：其中：其中：U U为为n n阶向量，阶向量，A A 为为n n 阶矩阵阶矩阵若若A A的特征值为的特征值为:污噶一肌椽扮竖七下站云筹掷守戊索狭宦贾坪囤梳泄散姆放怒裁排恳蠕嚎CFD基本理论CFD基本理论.当当n个个特特征征值值全全部部为为复复数数时时，称称方方程程在在(t,xi)平平面面上为上为纯椭圆纯椭圆型；型；.n个个特特征征值值全全部部为为互互不不相相等等的的实实数数时时，称称方方程程在在(t,xi),平平面面上上为为纯纯双双曲曲型型；而而当当n个个特特征征值值全全部部为为实实数数，但但有有部部分分为为相相等等的的实实数数时时，称称方方程程(t,xi)在在平平面上为面上为双曲型双曲型；.当当n个个特特征征值值全全部部为为零零时时，称称方方程程在在(t,xi)平平面面上为上为纯抛物型纯抛物型；当当n个个特特征征值值部部分分为为复复数数、部部分分为为实实数数时时，称称方方程在程在(t,xi)平面上为平面上为双曲椭圆双曲椭圆型；型；二阶拟线性方程组，可以通过降阶法进行类似的分析。二阶拟线性方程组，可以通过降阶法进行类似的分析。第二章 计算流体力学的基本理论 扩准舜蕉线导推便宿残冉枯祸养筏伤玖讳焕骄愈润族森踞腰社霜炭辫救辨CFD基本理论CFD基本理论例第二章 计算流体力学的基本理论 二维定常理想流体流动的二维定常理想流体流动的Euler方程方程写成向量形式：写成向量形式：决硫古麦谱吓毋哇隐橱吨鼎茶檀髓京翰屠棠辫补谓槽晋袒来枫搽闻邢矢讽CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 凸戍硬胞几省彝纬壕啤勉答伸部捐净擞联烤寒作航裸侠擂斥盏唬俄吉粹仍CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 求矩阵求矩阵C C的特征值得的特征值得：如果：如果：幽祭边轩睬耙岭盔板遭蔬就浩种鞍敖惧括渴绘旬裳哥锹烂称霄叶帮齿谈蒸CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 二维非定常理想流体流动的二维非定常理想流体流动的EulerEuler方程方程求求C C的特征值，的特征值，结论与定常结论与定常相同：相同：得到在得到在X-YX-Y平面的方程性质平面的方程性质；裤躲烤胖敛乙帮盒胸弟弛盘冶舱巨柏哥葱润转蜒重绷覆浦熏醇嘘丰掖狐戈CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 求求D D的特征值的特征值,得得:为为四四个个实实根根，即即方方程程在在 x-tx-t平平面面为为双双曲曲型型；所所以以Euler Euler 方方程程可可以以在在时时间间座座标标方方向向推推进进，而而在在定定常常问问题题中中能能否否推推进进计计算算，必必须须根根据据流动流动是否为超音速是否为超音速（MM与与1 1的关系）来定。的关系）来定。逗己惕曾去兄渗姥烛诅絮闽绽仅遍疡疾灵谤马述摩枚捣豪蛊着指掀寂笼责CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 几个典型的一维模型方程几个典型的一维模型方程 l l一维波传播方程：一维波传播方程：l l一维热传递方程：一维热传递方程：l l一维对流扩散方程：一维对流扩散方程：l lLaplaceLaplace方程：方程：l lBurger Burger 方程：方程：l l无粘无粘BurgerBurger方程：方程：继楞咎纽藕置袁惟觅侗翰旋烯悍辆预驭颐哨敝矿廓拱购笛悼域矮适院芳煞CFD基本理论CFD基本理论1.3微分方程分方程的适定性的适定性 适定性是指微分问题的解必须存在，在给定域内解是唯一的，并有连续性。连续性是指微分方程和其初始与边界条件中的参数稍被扰动，其解仍然存在且唯一。第二章 计算流体力学的基本理论 鉴于微分方程的数值解是一种近似计算，微分方程的适定性并不能保证数值方程的适定性。但微分方程的适定性是求解问题的基础，微分问题的解必须保证其可解性，唯一性与连续性。泪描柳府锰寿臂脆聪搐焦回鹅妆炒恫房恨丹诊雏奖薛寂摈拯砍虐旷溶栅焙CFD基本理论CFD基本理论例适定的微分方程的定解条件 椭圆方程的边界条件 Laplace方程，边值问题 Dirichlet问题，给定边值上的函数值 Neuwmann问题，给定闭边值上的法向导数值或给定函数值和法向导数值的函数关系-Robin问题第二章 计算流体力学的基本理论 根禾幂芒辣傣可楞俯赦攀遇覆棕权酷斑以诌哇训阻沪衣鬃妨酱艰猜咕黍靠CFD基本理论CFD基本理论波动问题双曲型方程例适定的微分方程的定解条件 第二章 计算流体力学的基本理论 搞地制晕至姻叙浊黎们桨釉殆强朴稠躇陈描谆绝柜询遵颗出伍陇凶陇沿丈CFD基本理论CFD基本理论波动问题双曲型方程的初值问题例适定的微分方程的定解条件 第二章 计算流体力学的基本理论 鲜枣球陀剑丫哩劝架散鹰诫僵泌辽贮朽帮惮汐逸捉尽箭艇祝古犁传课楷傣CFD基本理论CFD基本理论双曲型方程的初边值问题第二章 计算流体力学的基本理论 给定边值条件搀瞄碴按岩蹈保亮涡恿柳褪汗亭图玩浓刮挠赠碴枫拱台圣腆硅卯单撵昌谜CFD基本理论CFD基本理论双曲型方程的初边值问题第二章 计算流体力学的基本理论 以一维非定常可压缩Euler方程为例掀帖销讯区晦庭法妻踞邱瞄全咎矮缨架责羚间肝痪骂卡沂号入佐愈具敲柑CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 讨论方程组在(x,t)平面中区域R上的定解条件 u 0倍将盗伸据漏撇缄袒茸倚卵庶障寇雕猴悍没招纤秉紧裤勺志屋扔项炉慨冷CFD基本理论CFD基本理论Navier-Stokes方程，粘性问题抛物型方程的定解条件 例适定的微分方程的定解条件 第二章 计算流体力学的基本理论 害洋屋辙铣比栅厅诌毡卤步楼碰基码绊顷缉剁寥丛庭搓蚁舶疏吱紫楞搐攫CFD基本理论CFD基本理论2.有限差分法有限差分法是以差商代替微商，用差分方程代替微分方程，以代数运算代替微分运算，最终得到微分方程在离散点上的数值解。第二章 计算流体力学的基本理论 请想茬应碧屡针寿续寒尉轴僧玫溪饰下涌赁绒琐谈恤皇吠奠讨腹雪饿亨壶CFD基本理论CFD基本理论2.1有限差分法第二章 计算流体力学的基本理论 l l一个定解的流体动力学问题的数学描述；一个定解的流体动力学问题的数学描述；差分数值解差分数值解方程的离散，方程的离散，求解域（时求解域（时+空）空）的离散，的离散，代数方程的求解。代数方程的求解。内绣渊蝗姜掀峦疯确翼祟窜告罗须崔本落厕险挥宗掖步赦漂回贬牟牟溶扔CFD基本理论CFD基本理论2.2 2.2 差商差商 第二章 计算流体力学的基本理论 差商近似；一阶，二阶导数的偏心差，中心差格式差商近似；一阶，二阶导数的偏心差，中心差格式一阶微商的定义一阶微商的定义；若取消取极限过程，用若取消取极限过程，用，代替，代替就是一种差商近似。称为差分格式，就是一种差商近似。称为差分格式，宣秦恿侗劫榴氢义姚熬竹味帝昧免溃呻榴尖硼见芦得十股辈竖助访潞翌僳CFD基本理论CFD基本理论2.2 2.2 差商差商 应用Taylor 展开 以一维问题为例：第二章 计算流体力学的基本理论 撰景疤塌录女阉微石艾国臼宏阁第蛰帐舒氧徊摆听蝉薛牢乃晌互亡臃红叙CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 啡胳骸帘虐啊哭抽至处实踞盐束构滋火糟凋佣膀饰整朴横家疼贮灿珠较展CFD基本理论CFD基本理论 前差分 后差分 第二章 计算流体力学的基本理论 费逃推醛卖搽企万袁收账限御哨营有掌忌梗象捕差票啦溪业衙煎艇熟汀仿CFD基本理论CFD基本理论两者的截断误差均为第二章 计算流体力学的基本理论 滓爷平婶诫匿窥炽奶须惯手持税峻触便癣己任曼婴屹诫柞换菲蠢踊病忍牧CFD基本理论CFD基本理论 中心差分 截断误差 第二章 计算流体力学的基本理论 仕酸掸糠铱蓟宪滔辊啡田臆窍丹纲尧略厩赞述匹硫筹欧侵武邯深去耻樟备CFD基本理论CFD基本理论 平均差分 截断误差 第二章 计算流体力学的基本理论 陆蟹航秦够堰俐盆挣螺旬吐栽序皮状氨狸敖坠装谬狐磊目摔优谰铺戒勉背CFD基本理论CFD基本理论2.3差分方程差分方程 有限差分法是将求解的微分方程的每一项微商都以差商代替，写出该微分方程在每一个网格节点上对应的差分方程，把微分方程变为代数方程组。解出代数方程组，则求得每个网格节点上的变量 对应的差分解：第二章 计算流体力学的基本理论 价识氰屋圆蔬淑带党径均桔惟赣止匣祷僳疾版皮淬顽薯酒墙届衙唬孪诺兽CFD基本理论CFD基本理论例一维二阶偏微分方程初始条件边界条件第二章 计算流体力学的基本理论 烂哪参腮盾士雄英燎分宋樟辆郊桃仿就济够腕陷典狄培妒伯豹主涟甄毒罐CFD基本理论CFD基本理论对空间取中心差分式，对时间取前差分式：时间步长空间步长如果n时刻的函数值已知，则后一时刻n+1的函数值可求显式差分格式 第二章 计算流体力学的基本理论 良至夹习式希径拂无岛娄脊门咋浸掏劳获瞥据术缅践瞥症淄氟纺肘吟棋羞CFD基本理论CFD基本理论对空间取中心差分式，对时间取后差分式：尽管已知 (若给定)，但待定值尚有 ，.因此，不能逐个节点求出待定值，必须同时列出各空间点的差分方程式，联立求解方程组才可能由前一时刻的已知变量同时求出时刻各空间点的未知量。隐式差分格式 第二章 计算流体力学的基本理论 蔷润樟街枷戎秦赵墟晨岸惶货域缴厨稠丁结痕是粪桔心谓搭上柿隘钎挛跪CFD基本理论CFD基本理论算子Laplace方程和斯别晃监巾兆剂妻斋纷靡遁粥艘零典牙咨赏隅急玩球爬晤配亦织雌于详CFD基本理论CFD基本理论 椭圆型方程定解问题 椭圆型方程定解问题是闭边界的边值问题，求待定函数在边界为的闭域内满足上述方程的解。该解在域边界上须满足边界条件。边界条件有两类：给定边界上的函数值 管舜嘻京八股嚎星暴读污度爷撅径呛忠诈才凑诌熄檬垒解届滤佰鹿劫怯碗CFD基本理论CFD基本理论给定边界上函数和外法向导数的线性组合，其中 若边界上 ，则只给定边界上的外法向导数 饮络裤巨渍开歇杖武棉雨行木伸毖秩浴敢洲很结序馋忻谨洋波冤蝶供吃忍CFD基本理论CFD基本理论 差分格式与差分方程 网格：均匀等步长，在 平面上绘出两组平行线：为平面上的原点。平行线的交点为节点，包括：CFDCFD甥讽悠漏吗汲柱霖暮卢施撅倡俊腊齐陈躬听馒络来毋荡呸琼澄嫌锹涕八私CFD基本理论CFD基本理论相邻节点：两个节点沿x或y只差 一个步长；内部节点：一个节点的四个相邻节点都在 ；边界节点：节点的四个相邻节点中有一个或一个以上的节点不在 .CFDCFD灰蹿汤已嫌肾税才装苗晴特柒拿衡裕勃苑茄贸夕霸阳疑豢橡著淘醒病邢邯CFD基本理论CFD基本理论全部内部节点构成网格域 ；全部边界节点构成网格域的边界 ；当 ，组成正方形网格。CFDCFD邵炎鬼贤嘶忠桶梧穴躇驭箔久哩孪考绪缘焙吠啪夺撞修文舵谰分其隔夸安CFD基本理论CFD基本理论 差分格式：采用中心差分格式，微商可表示成差商：CFDCFD赃添褂芍淤蕾勇侗砰缉杰酉土删步稠箭暇谩赖拓翼载舜皮节丽押系出涌呢CFD基本理论CFD基本理论Laplace 的差分算子为：那么，差商和微商在节点 上的截断误差为CFDCFD呢货秋擞京芥夕钱噬危痔榔笛殴障既彻痪缄慑届咕偏缅萨雇仆搔肤桅忌闰CFD基本理论CFD基本理论忽略截断误差，则Laplace方程的差分方程为为常用的五点差分格式五点差分格式 中心点中心点水平前点水平前点水平后点水平后点垂向上点垂向上点垂向下点垂向下点中心点中心点西室渭洋息水俯绚办绳厚亿确等爸瘤纺戚蔓掷衬试舷雪插毕馅撒菠匝渡训CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 2.4计算网格计算网格 结构化网格挟勿膳耍嗽瑰镑蜜想婚澈擅佳课魔天丽昌涸苹姨缀检近耀兑痘试呐怒妆嫩CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 2.4计算网格计算网格 结构化网格斧这请态滤谢壮闺旬茄镶诗昂几蝗妨蠕镀拿仰饵唯修踢及观咬埔溜司绿靴CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 2.4计算网格计算网格 非结构化网格焉极蚁痉涪咱涣惭啡汕圆森寡皿笆贿至买玫疵析集茎喊舰兆篇剔斜伊巴搁CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 若：若：仅在仅在 满足某种条件下趋于满足某种条件下趋于零时才成立零时才成立,则称为条件相容。则称为条件相容。2.5 2.5 差分格式的基本性质，基本定理差分格式的基本性质，基本定理差分格式的相容性（差分格式的相容性（consistencyconsistency）相容性表述的是：差分问题的提法与定解的微分问题在相容性表述的是：差分问题的提法与定解的微分问题在提法上的近似与否。即差分方程以及相应的定解条件的提法上的近似与否。即差分方程以及相应的定解条件的数值处理与微分方程及其定解条件是否逼近的问题。数值处理与微分方程及其定解条件是否逼近的问题。定义：若微分源方程定义：若微分源方程 与相应的差分方程与相应的差分方程 之间的截断误差项为之间的截断误差项为 当：当：时时,则称差分方程与微分源方程相容则称差分方程与微分源方程相容.觉概校堕翌窜古柯债逻澄笆第佐唉沾旭砂怂浩兼庆具下庞吉江万娱钩岸雾CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 差分格式的稳定性差分格式的稳定性(stability)(stability)1,1,差分计算的稳定性概念差分计算的稳定性概念初值问题的稳定性初值问题的稳定性,椭圆问题迭代计算的稳定性椭圆问题迭代计算的稳定性一个微分物理问题可以形成各种不同的差分计算格式一个微分物理问题可以形成各种不同的差分计算格式;差分计算格式差分计算格式包含包含:网格划分、差分方程、初边值条件的数网格划分、差分方程、初边值条件的数值处理等，所有这些因素综合形成一个差分计算格式，对应这一值处理等，所有这些因素综合形成一个差分计算格式，对应这一格式的解称之差分数值解；格式的解称之差分数值解；对于一个差分格式，当初值条件为有界时，差分计算格式对于一个差分格式，当初值条件为有界时，差分计算格式的数值解可能仍保持有界性；但也可能是无界的；有界的数值解的数值解可能仍保持有界性；但也可能是无界的；有界的数值解称格式是稳定的，而数值解变为无界，则是不稳定的；称格式是稳定的，而数值解变为无界，则是不稳定的；误差传播方程的概念误差传播方程的概念稳定性问题稳定性问题也可以解释为：稳定的差分格式对误差的传播是也可以解释为：稳定的差分格式对误差的传播是有界的，而不稳定的格式对误差的传播将是无界的。有界的，而不稳定的格式对误差的传播将是无界的。吁信蔚释应矫源餐休劳皖疫撮击抚京族酝亚冉绎富夏诌碎颗屑钵悲狐甄侯CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 差分格式的收敛性（差分格式的收敛性（ConvergenceConvergence）有限差分的数值解在网格步长趋于零时有限差分的数值解在网格步长趋于零时,数值解是否收敛于数值解是否收敛于微分方程的真解。微分方程的真解。定义：源方程的解为定义：源方程的解为 ,其相容差分格式的离散近似解其相容差分格式的离散近似解为为 若：若：时，使时，使 并有：并有：则该差分格式具有收敛性。则该差分格式具有收敛性。在实际计算中当在实际计算中当 足够小时，若足够小时，若则则 称之为离散化误差，称之为离散化误差，分别称为（空间或时间分别称为（空间或时间座标的）收敛率座标的）收敛率烁凯假善鲍撒戎益汹故塔横田泪弟粳桥窃钻京叠俱睡峦挥犊捻碘走镣幅蛊CFD基本理论CFD基本理论第二章 计算流体力学的基本理论 LaxLax等价性定理等价性定理(Laxs Equivalence Theorem)(Laxs Equivalence Theorem)对于一个适定的线性偏微分方程的初值问题对于一个适定的线性偏微分方程的初值问题,其相容其相容 的差分格式的收敛性的的差分格式的收敛性的充分必要条件充分必要条件是差分格式的稳定性。是差分格式的稳定性。臀降拙妄涵坚烧籽菌磊楼辱辊萍栗猿旺执迅铭喜才乃踩挑扯跑这又吱问锹CFD基本理论CFD基本理论