收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

非结构网格FVM在复杂几何结构的湍流反应流计算中的应用研究

雷国东  
【摘要】: 非结构网格有限体积法适合应用于复杂几何结构的流场计算,而SIMPLE系列算法适合应用于工业上广泛存在的低速流、不可压流计算,欧拉-拉格朗日方法适合应用于含稀疏颗粒相的流场计算。本文正是以这些理论为基础展开工作,包括CGNS网格前后处理,压力修正算法的实现和拓展,化学组分源项模型,颗粒相与气相之间的相互作用等。分述如下: 本文提出了两种将CGNS非结构网格点数据结构转化为非结构网格面数据结构的快速算法,以作为本文操作非结构网格的基础:内部网格面的整数hash表算法和内部双面壁边界网格面的实数hash表算法。由于CGNS非结构网格存储的是网格结点数据,例如网格结点坐标、每一个网格由哪些结点构成等数据信息;而有限体积法离散过程中,需要用到大量的面积分,所以更多的需要网格面数据结构。内部网格面的构造根据网格面两侧的网格在这个公共网格面上的结点序号相同的特点设计算法,双面壁网格面则根据重合在一起的两个边界网格面中心距离接近于零,或者说小于某一个特征小量,比如该网格面面积的算术平方根来设计算法。两种hash表算法具有快速,内存消耗量小于后续有限体积求解器内存峰值的特点,使前处理算法不会成为整个计算的瓶颈。 在压力修正方程的构造过程中,考虑温度修正量对密度修正量的影响,通过静焓方程构造温度修正量与压力修正量的关系,而这种做法的实际效果相当于增强了压力修正量方程线性代数方程组的主对角线。对于可压缩流计算,在流量计算时,在线性插值和迎风格式的基础上,网格面上的密度值采用低阶和高阶的混合格式,动量方程中的压力项处理为滞后的压力梯度项,即线性插值计算方法,和网格面上的压力积分之和,网格面上的压力值则可以采用低阶和高阶混合格式。这些做法兼顾了分辨率和计算的稳定性。 在喷雾燃烧数值模拟中考虑了湍流反应组分源项模型如EBU模型、ED模型、Arrinius/EBU/ED模型和Presumed PDF模型。对于颗粒相和气相之间的相互作用,本文一是将颗粒相对气相的反作用阻力中的部分阻力系数剥离,由拉格朗日系统转入欧拉系统,加入到气相动量方程线性代数方程组的主对角线上,这一系数为正实数,增强了主对角线的系数,也增强了含颗粒流的气相迭代计算稳定性;二是在采用时间推进方法计算稳态流场过程中,对耦合的颗粒相与连续相之间计算采用显式和隐式的计算方法进行了研究,这样做的目的在于相对传统的显式时间计算方式,可以采用一个较大的时间步长值进行迭代,从而减少时间推进时的计算量,使得喷雾燃烧数值模拟的计算稳定性得以增强同时减少计算量。 在课题组研制的自主知识产权流体力学软件开发过程中,本文作者承担了部分功能模块的研发工作,取得了上述研究成果。并就一些算例与商用软件和经典计算结果予以对比,验证了本文算法的正确性。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 李春娜;叶正寅;王刚;;基于二维非结构网格的GMRES隐式算法[J];西北工业大学学报;2007年05期
2 赵庆贺;白俊强;许晓平;袁晓军;;二维非结构网格Euler方程高效解法[J];科学技术与工程;2008年01期
3 周涤宇;刘杰;;一种改进的求解输运方程的并行Sn扫描算法[J];计算机工程与科学;2008年04期
4 彭恒初;吴望一;;基函数法在黏性可压缩流动中的研究与应用[J];北京大学学报(自然科学版);2009年06期
5 黄士力;徐健;汪德爟;陆东燕;;风暴潮变分同化模型[J];水动力学研究与进展A辑;2010年04期
6 朱军政;曹颖;;FVCOM模型在象山港三维潮流盐度计算中的应用[J];海洋环境科学;2010年06期
7 赵春霞;王凯;马海峰;魏延福;程舜;;二维水流数学模型在涉水工程中数值模拟研究[J];水科学与工程技术;2011年03期
8 张雅,刘淑艳,王保国;使用RSM模型对旋风除尘器内湍流各向异性的数值模拟[J];工程热物理学报;2005年S1期
9 任登凤,谭俊杰,张军;GMRES方法在含动边界流场中的应用[J];力学与实践;2005年05期
10 王学德;伍贻兆;夏健;;三维非结构网格DSMC方法的实现及其应用[J];计算力学学报;2007年02期
11 韦喜忠;黄振宇;洪方文;;基于非结构网格的螺旋桨周围流场大涡模拟[J];水动力学研究与进展A辑;2008年04期
12 巨海涛;吴宏春;;非结构网格中子输运方程角度相关再平衡加速算法[J];核动力工程;2008年05期
13 雷国东;柳贡民;张文平;朱明刚;;考虑密度与温度耦合关系的非结构网格SIMPLEC算法[J];计算物理;2009年01期
14 吴望一;谢文俊;;一种新型的数值计算方法——基函数法[J];中国科学(G辑:物理学 力学 天文学);2009年08期
15 吴望一;林光;;多项式基函数法[J];应用数学和力学;2009年09期
16 卢康明;李光炽;;非结构网格浅水方程隐式解法[J];水动力学研究与进展A辑;2010年02期
17 张亮亮;王春武;;非结构网格上求解二维H-J方程的一种WENO格式[J];高校应用数学学报A辑;2010年04期
18 刘枫;冯定华;丁国昊;潘沙;李桦;;任意多边形的非结构网格快速生成算法研究[J];国防科技大学学报;2010年06期
19 于守兵;王万战;余欣;;基于非结构网格的三维浅水模型[J];河海大学学报(自然科学版);2011年02期
20 朱敏;杨旭东;;一种导管螺旋桨非定常干扰流CFD方法[J];航空计算技术;2011年03期
中国重要会议论文全文数据库 前10条
1 付峥;刘凯欣;;可应用于真正非结构网格的CE/SE方法[A];北京力学会第17届学术年会论文集[C];2011年
2 王瑞利;林忠;闫伟;温万治;;结构网格与非结构网格的联合使用技术[A];中国工程物理研究院科技年报(2002)[C];2002年
3 郭正;刘君;瞿章华;;动网格技术在二维非结构网格中的实现[A];第十届全国计算流体力学会议论文集[C];2000年
4 王振亚;叶友达;;基于非结构网格的隐式方法的比较[A];计算流体力学研究进展——第十一届全国计算流体力学会议论文集[C];2002年
5 王平;朱自强;;曲面非结构网格生成法研究[A];第十届全国计算流体力学会议论文集[C];2000年
6 王瑞利;;非结构网格物理量重映方法[A];中国工程物理研究院科技年报(2000)[C];2000年
7 王刚;叶正寅;陈迎春;;三维非结构网格上的N—S方程解法研究[A];计算流体力学研究进展——第十一届全国计算流体力学会议论文集[C];2002年
8 倪国喜;姚彦忠;;结构与非结构网格局部优化方法[A];中国工程物理研究院科技年报(2005)[C];2005年
9 孙树立;陈斌;刘剑飞;袁明武;;非结构动网格生成技术研究进展[A];中国力学学会学术大会'2009论文摘要集[C];2009年
10 孙树立;陈斌;刘剑飞;袁明武;;混合网格变形技术研究进展[A];中国力学学会学术大会'2009论文摘要集[C];2009年
中国博士学位论文全文数据库 前10条
1 李万爱;非结构网格高精度数值方法的若干问题研究[D];清华大学;2012年
2 王兰;超燃冲压发动机整机非结构网格并行数值模拟研究[D];中国空气动力研究与发展中心;2007年
3 耿艳芬;城市雨洪的水动力耦合模型研究[D];大连理工大学;2006年
4 王志力;基于Godunov和Semi-Lagrangian法的二、三维浅水方程的非结构化网格离散研究[D];大连理工大学;2006年
5 王学德;高超声速稀薄气流非结构网格DSMC及并行算法研究[D];南京航空航天大学;2006年
6 何冰;基于Boltzmann方程三维可压缩高速流动并行算法及其应用研究[D];上海大学;2008年
7 汪迎春;运动界面模拟技术及在环境分层流问题中的应用[D];河海大学;2007年
8 潘勇;高超声速流场磁场干扰效应数值模拟方法研究[D];南京航空航天大学;2007年
9 陈锐林;强风作用下列车脱轨分析[D];中南大学;2008年
10 刘晨;复杂燃烧流场数值模拟方法研究[D];南京航空航天大学;2009年
中国硕士学位论文全文数据库 前10条
1 赵滋阳;非结构网格下气动热力耦合数值方法研究[D];南京航空航天大学;2009年
2 兰大辉;非结构网格上求解流场的初步研究[D];哈尔滨工业大学;2010年
3 戴阳豪;基于非结构网格的二维水流数值模拟研究[D];哈尔滨工程大学;2012年
4 刘光生;非结构网格下河流对近岸水环境影响的数值模拟[D];大连理工大学;2010年
5 赵延生;非结构网格的ENO有限体积方法研究[D];国防科学技术大学;2004年
6 吴晶;无结构网格的生成及其在有限体积法中的应用[D];河海大学;2006年
7 王靖宇;涡轮内非结构网格气动热力耦合数值计算方法研究[D];南京航空航天大学;2012年
8 朱思美;非结构网格上解Hamilton-Jacobi方程的数值方法研究[D];国防科学技术大学;2004年
9 王永;非结构化网格生成技术及在SIMPLE算法中的应用研究[D];天津大学;2005年
10 王文全;节能离心泵全流道数值模拟及性能预测[D];昆明理工大学;2005年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978