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