收藏本站
《中国空气动力研究与发展中心》 2013年
收藏 | 手机打开
二维码
手机客户端打开本文

基于混合网格的高阶精度DG/FV混合算法研究

李明  
【摘要】:非结构/混合网格具有适于离散复杂几何外形和易于网格自适应从而更精细和高效地捕捉流场细微结构等优点,因此近年来其网格生成方法和求解Euler、N-S方程的计算方法得到了迅速发展。目前几乎所有的CFD商业软件及绝大部分专业流场解算器,无论其使用有限体积方法(FVM)、有限差分方法(FDM)、有限元方法(FEM)等等都是基于二阶精度方法。这些方法在实际工程应用中已取得很大的成功,然而对于很多问题,如计算气动声学(CAA)、旋涡主导的流动问题、湍流的大涡模拟(LES)和直接数值模拟(DNS)等等,需要采用高阶精度计算方法(三阶或以上精度)。二阶精度方法在求解上述问题时存在对旋涡模拟的过度耗散和色散等缺点,需要使用低数值耗散和色散的高阶精度方法来提高旋涡和其它复杂的分离、非定常流动问题的预测精度。因此基于非结构/混合网格的高阶精度方法在近一二十年受到广泛的关注。 目前,大多数非结构/混合网格上的数值方法都源于有限体积方法或者有限元方法。在众多的高阶精度计算方法中,间断Galerkin有限元方法(DGM)因其计算精度高、提高精度阶无需拓展网格模板、并行算法容易实现、且具有很多良好的数学性质等优点,已成为CFD领域卓有成效的数值方法之一。然而间断Galerkin方法也存在一些缺点,主要是内存需求量和计算CPU资源耗费巨大。相对于传统的二阶精度DG方法,二阶精度有限体积方法因为其不需要求解额外的高阶基函数的自由度(DOFs)和数值积分,其在使用相同网格时的计算机内存和CPU资源需求量都要小很多。但是对于高阶有限体积方法来说,多维的高阶重构需要的模板很大。因而不具有DG方法的紧致特性,这在大规模并行计算时会带来大量的信息交换。另外,高阶重构模板的选择具有不确定性和任意性,从而使得重构非常复杂,且重构矩阵的计算也要耗费大量的CPU资源,特别是在三维情况下。 为了解决上述问题,张涵信提出了一个结合有限体积方法和有限元方法优点的具有三阶精度的混合方法。张来平、刘伟等提出了“静态重构”和“动态重构”的概念用于构造高阶精度数值方法,并基于静动态“混合重构”的思想,构造了一类针对三角形和三角形/笛卡儿混合网格的高阶精度方法用于求解守恒律方程,命名为DG/FV混合方法。在这类混合方法中,每个单元中的分布多项式的低阶系数使用通常的DG方法进行求解(称为动态重构);对于高阶系数,类似有限体积方法,利用本单元和邻近单元上的低阶系数信息进行重构得到(称为静态重构)。此混合方法构造简单、适于复杂外形、方便推广获得更高阶格式。另外,此方法模板紧致,只和相邻单元进行数据交换,因而通讯量小,非常适合于并行计算。 本文理论分析了上述高阶精度DG/FV混合方法的耗散、色散特性及稳定性条件;并将其推广应用于二维混合网格下的N-S方程求解;定性理论分析和数值比较了其数值精度和计算效率。同时,为了进一步加速定常流的收敛速度,发展了基于Newton/Gauss-Seidel迭代的隐式计算方法。计算结果表明,此DG/FV混合方法达到了设计的4~5阶精度,且隐式方法相对于显式Range-Kutta方法,收敛速度提高约1~2个量级。 本文共分为六章。 第一章中,首先简要回顾了计算流体力学中基于非结构网格/混合网格下的计算方法及其高阶精度、高分辨率格式最新研究进展;然后讨论了几种常见的基于非结构/混合网格的高阶精度数值方法的优缺点,指出了可以进行改进的地方;最后对本文工作做一简要介绍。 第二章详细介绍了DG/FV混合方法,包括方法的基本思想、具体构造方法、空间重构算法、显隐式时间离散方法、数值通量格式、高斯数值积分方法、限制器、边界条件及曲边界处理方法等。 在第三章中讨论了DG/FV混合方法的数值特性,包括方法的数值色散和耗散特性;给出了方法的稳定性条件;并利用线性和非线性标量方程对DG/FV混合方法的计算精度进行了验证;定性理论分析和数值验证了此方法的计算量和内存消耗方面的优势。 第四章和第五章测试了一些数值算例来验证本文DG/FV混合方法的精度及计算效率。在第四章中,对一些无粘算例进行了数值模拟,包括等熵涡流动、亚声速圆柱绕流、Bump流动、一维Sod激波管问题、Shu-Osher问题及前台阶的超声速流动问题等。第五章中给出了一些粘性流动算例,如平板间的Couette流动、平板层流流动、粘性圆柱绕流、方腔流动及NACA0012翼型绕流等。上述数值计算中,将计算精度及效率与同阶DG方法结果进行了对比分析,并且了选择典型算例对比分析了显隐式时间离散方法的计算效率和隐式方法参数对计算效率的影响。 第六章为结束语,给出本文工作总结及未来可能的研究方向。
【学位授予单位】:中国空气动力研究与发展中心
【学位级别】:博士
【学位授予年份】:2013
【分类号】:O35

知网文化
【参考文献】
中国期刊全文数据库 前9条
1 蔡庆东,温功碧;二维可压无粘流的自适应流量修正有限元解[J];航空学报;1994年11期
2 段占元,童秉纲,姜贵庆;高分辨率有限差分-有限元混合方法及其在气动热计算中的应用[J];空气动力学学报;1997年04期
3 张来平;刘伟;贺立新;邓小刚;张涵信;;一种新的间断侦测器及其在DGM中的应用[J];空气动力学学报;2011年04期
4 黄兆林,李宏玮,毛国良;再入飞行器高超声速绕流的一步和二步有限元计算[J];空气动力学学报;1994年02期
5 朱刚,沈孟育;跨声速叶栅粘流计算的多级Taylor-Galerkin有限元法[J];空气动力学学报;1995年04期
6 张来平,张涵信,高树椿;矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用[J];空气动力学学报;1998年01期
7 贺立新;张来平;张涵信;;间断Galerkin有限元和有限体积混合计算方法研究[J];力学学报;2007年01期
8 张来平,张涵信;NND格式在非结构网格中的推广[J];力学学报;1996年02期
9 张来平;贺立新;刘伟;李明;张涵信;;基于非结构/混合网格的高阶精度格式研究进展[J];力学进展;2013年02期
【共引文献】
中国期刊全文数据库 前10条
1 彭轩;马驰;;折线布置采空区流场特征与三角煤柱防灭火技术研究[J];安徽建筑工业学院学报(自然科学版);2011年05期
2 王宗木;一个改进的非线性方程求解的迭代方法[J];安徽建筑工业学院学报(自然科学版);1996年04期
3 孙玉香;许勇;;简易性在计算方法设计中的应用研究[J];安徽师范大学学报(自然科学版);2008年05期
4 ;Numerical simulations of compressible mixing layers with a discontinuous Galerkin method[J];Acta Mechanica Sinica;2011年03期
5 郝庆一;;二维移动网格方法中边界网格点移动的数值研究[J];安庆师范学院学报(自然科学版);2008年02期
6 王婷婷;马庆元;郭继平;;平板层流边界层近似速度分布计算方法的改进[J];鞍山科技大学学报;2006年06期
7 牛晓奇;;理想洛伦兹曲线拟合及其方程特征[J];安阳师范学院学报;2009年02期
8 王训宏;王快社;徐可为;胡锴;;搅拌摩擦焊接过程塑性流动形态与试验分析[J];兵器材料科学与工程;2006年06期
9 崔嵬;非标准插值多项式的截断误差分析[J];保定师范专科学校学报;2004年02期
10 熊艳,毛军,薛琳;平板层流边界层定常流动存在相似性解区域的分析[J];北方交通大学学报;2003年04期
中国重要会议论文全文数据库 前10条
1 李亚军;姚军;;裂缝型介质等效连续模型流动模拟方法研究[A];渗流力学与工程的创新与实践——第十一届全国渗流力学学术大会论文集[C];2011年
2 郑耀;;计算技术在航空航天领域中的应用[A];首届全国航空航天领域中的力学问题学术研讨会论文集(下册)[C];2004年
3 杨炳渊;史哓鸣;梁强;;高超声速有翼导弹多场耦合动力学的研究和进展[A];第十届全国空气弹性学术交流会会议论文集[C];2007年
4 刘应华;刘冬欢;郑小平;;不连续温度场问题的间断Galerkin方法[A];损伤、断裂与微纳米力学进展:损伤、断裂与微纳米力学研讨会论文集[C];2009年
5 占欢;范育新;宫继双;齐海帆;程晓军;;扇环形截面管内爆震波特性数值模拟[A];第十四届全国激波与激波管学术会议论文集(上册)[C];2010年
6 季刚;钱战森;高振勋;;横向声速射流与超声速来流干扰流动声场数值模拟[A];第十四届全国激波与激波管学术会议论文集(下册)[C];2010年
7 徐驰;王友军;束忠明;;坨面形状对CVD合成石英玻璃沉积速率的影响[A];2010全国玻璃技术交流研讨会论文集[C];2010年
8 郑耀;;高端计算技术在深空探测科学研究中的作用[A];中国宇航学会深空探测技术专业委员会第一届学术会议论文集[C];2005年
9 孙岩;;模压成型过程中的复合材料在模腔内的力学分析[A];第九届中国摩擦密封材料技术交流暨产品展示会论文集[C];2007年
10 刘久玲;李顺洲;刘明华;何世堂;;声表面波气相色谱仪响应机理及检测限分析[A];泛在信息社会中的声学——中国声学学会2010年全国会员代表大会暨学术会议论文集[C];2010年
中国博士学位论文全文数据库 前10条
1 李佳;载人潜器阻力性能的数值和试验预报及外形优化研究[D];哈尔滨工程大学;2010年
2 刘志平;无陀螺捷联惯导系统若干关键技术研究[D];哈尔滨工程大学;2010年
3 邓锐;阻流板对双体船水动力性能影响的数值研究[D];哈尔滨工程大学;2010年
4 贺克羽;新型控制棒水力驱动系统研究[D];哈尔滨工程大学;2009年
5 顾璇;近地/海运动过程中的空气动力学分析[D];哈尔滨工程大学;2009年
6 黄燕;碟形越浪式波能发电装置的水动力性能研究[D];中国海洋大学;2010年
7 韩艳;海底管线的冲刷与导流防护技术研究[D];中国海洋大学;2010年
8 王春林;旋流自吸泵设计及内部流动研究[D];江苏大学;2010年
9 张华;冲击荷载作用下岩石动态损伤特性研究[D];昆明理工大学;2009年
10 闫洁;三维气固两相平板间射流中相间耦合作用的直接数值模拟[D];浙江大学;2009年
中国硕士学位论文全文数据库 前10条
1 赵文武;淮南矿区典型煤层水力冲孔技术参数研究及消突效果考察[D];河南理工大学;2010年
2 张彦奇;超细水泥渗透特性微观试验研究及理论分析[D];山东科技大学;2010年
3 刘博;叶片式流体设备数字化研发与优化研究[D];郑州大学;2010年
4 李琦;某16V型柴油机相继增压性能分析[D];哈尔滨工程大学;2010年
5 李铁军;多级除尘器气动性能分析[D];哈尔滨工程大学;2010年
6 侯鹏;基于双曲型方程RKDG法在流场及声场中的数值模拟研究[D];哈尔滨工程大学;2010年
7 王祎;载机舰在风、浪中的气流场数值模拟研究[D];哈尔滨工程大学;2010年
8 张吕;液体火箭发动机喷管螺旋槽数控加工技术[D];大连理工大学;2010年
9 吕秀丽;离心泵内部三维流场数值模拟研究[D];辽宁工程技术大学;2009年
10 朱夺宝;两类具有分布时滞的离散传染病模型的稳定性[D];兰州大学;2010年
【二级参考文献】
中国期刊全文数据库 前10条
1 蔡庆东,温功碧;二维可压无粘流的自适应流量修正有限元解[J];航空学报;1994年11期
2 张涵信;无波动、无自由参数的耗散差分格式[J];空气动力学学报;1988年02期
3 徐国群,张国富;不可压Navier-Stokes方程组的SUPG有限元数值解[J];空气动力学学报;1991年03期
4 张涵信,贺国宏,张雷;高精度差分求解气动方程的几个问题[J];空气动力学学报;1993年04期
5 段占元,童秉纲,姜贵庆;高分辨率有限差分-有限元混合方法及其在气动热计算中的应用[J];空气动力学学报;1997年04期
6 张来平;刘伟;贺立新;邓小刚;张涵信;;一种新的间断侦测器及其在DGM中的应用[J];空气动力学学报;2011年04期
7 王岩;计算流体力学中有限元的强间断处理[J];空气动力学学报;1994年01期
8 黄兆林,李宏玮,毛国良;再入飞行器高超声速绕流的一步和二步有限元计算[J];空气动力学学报;1994年02期
9 朱刚,沈孟育;跨声速叶栅粘流计算的多级Taylor-Galerkin有限元法[J];空气动力学学报;1995年04期
10 贺立新;张来平;张涵信;;间断Galerkin有限元和有限体积混合计算方法研究[J];力学学报;2007年01期
【相似文献】
中国期刊全文数据库 前10条
1 张逸强;定常不可压缩粘性流动的一种高阶精度解法[J];力学学报;1987年S1期
2 赵兴艳,苏莫明,苗永淼;基于压力的高阶精度隐式格式的研究及应用[J];数值计算与计算机应用;2002年01期
3 孙光甫;奇异摄动问题的一致收敛离散方法解的一致高阶精度外推[J];应用数学和力学;1990年03期
4 王刚;张德良;刘凯欣;王景焘;;二维CE/SE算法改进及其高阶精度格式[J];计算力学学报;2008年06期
5 赵兴艳,苏莫明,苗永淼;高分辨率激波捕捉格式及其CFD应用[J];计算物理;2001年02期
6 郑华盛;赵宁;;一类高精度TVD差分格式及其应用[J];应用力学学报;2005年04期
7 吴开腾;张莉;;一类三维基于通量分裂的高精度MmB差分格式[J];兵工学报;2010年S1期
8 苏莫明,赵兴艳,苗永淼;基于预处理方法的流动数值分析[J];西北工业大学学报;2001年03期
9 郑华盛;徐伟;李曦;胡结梅;;一类基于非均匀剖分的二维高精度差分格式[J];南昌航空大学学报(自然科学版);2008年02期
10 郑华盛;赵宁;朱君;;二维非结构网格上的高精度有限体积WENO格式[J];空气动力学学报;2010年04期
中国重要会议论文全文数据库 前10条
1 邓小刚;;高阶精度耗散加权紧致非线性格式[A];全国流体力学青年研讨会论文集[C];2001年
2 邓小刚;张毅锋;毛枚良;张涵信;;高阶精度格式加速收敛研究[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(上)[C];2007年
3 吴开腾;张莉;;一类三维基于通量分裂的高精度差分格式的应用研究[A];中国力学学会学术大会'2009论文摘要集[C];2009年
4 吴开腾;张莉;;一类三维基于通量分裂的高精度MmB差分格式[A];第九届全国冲击动力学学术会议论文集(下册)[C];2009年
5 何洋洋;高静怀;;任意高阶间断有限元法在弹性波模拟中的应用研究[A];中国地球物理学会第二十四届年会论文集[C];2008年
6 高智;向华;;摄动有限体积方法和重构近似高精度的意义[A];第四届《海峡两岸计算流体力学学术研讨会》论文集[C];2003年
7 成娟;;可压缩欧拉方程高精度ENO守恒拉氏格式[A];全国计算物理学会第六届年会和学术交流会论文摘要集[C];2007年
8 Abolghasem Pilechi;Peyman Badiei;杜晶;;管道周围液流的数值模拟[A];2011年非开挖技术会议论文集[C];2011年
9 王景焘;刘凯欣;;基于时空守恒元解元方法的流体弹塑性问题计算[A];中国力学学会学术大会'2005论文摘要集(上)[C];2005年
10 詹杰民;董志;李毓湘;周锦添;韦永康;;传统浅水型和Boussi nesq型方程模拟海啸传播过程准确性的研究[A];第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会论文集[C];2009年
中国博士学位论文全文数据库 前10条
1 李明;基于混合网格的高阶精度DG/FV混合算法研究[D];中国空气动力研究与发展中心;2013年
2 张毅锋;高阶精度格式(WCNS)加速收敛和复杂流动数值模拟的应用研究[D];中国空气动力研究与发展中心;2008年
3 张青洁;色散方程的一类高精度并行算法[D];山东大学;2009年
4 沈俊;生物型体竞争模型的高分辨率高精度方法[D];中国科学技术大学;2007年
5 郑华盛;流体力学高精度数值方法研究[D];南京航空航天大学;2005年
6 李大明;有限体积方法在守恒律中的应用[D];浙江大学;2002年
7 王昆;复杂水流的高分辨率数值模拟[D];大连理工大学;2009年
8 刘梅林;节点间断伽辽金有限元方法及其在计算电磁学中的应用研究[D];南京航空航天大学;2011年
9 李兴良;球面阴阳网格与新型多离散矩有限体积方法的应用研究[D];南京信息工程大学;2007年
10 杨旻;几类有限体积元及有限体积格式的数值分析[D];山东大学;2005年
中国硕士学位论文全文数据库 前10条
1 王洁;反应扩散方程组的单调紧ADI方法[D];华东师范大学;2009年
2 赵延生;非结构网格的ENO有限体积方法研究[D];国防科学技术大学;2004年
3 陈华;基于径向基函数的有限体积方法[D];南京航空航天大学;2005年
4 李松;高阶精度耗散紧致差分格式的研究与应用[D];中国空气动力研究与发展中心;2007年
5 李彬;双曲型守恒律方程的高分辨率有限体积方法[D];国防科学技术大学;2004年
6 钮海峰;二维非结构网格大粒子有限体积方法[D];中国人民解放军国防科学技术大学;2003年
7 周文文;Sobolev方程的两类数值解法[D];内蒙古大学;2013年
8 黄师化;基于复合数值方法的水流动画模拟[D];安徽大学;2005年
9 钱凌志;对流占优扩散问题的三种数值解法研究[D];新疆大学;2010年
10 李付鹏;基于有限差分法和有限体积法的水流动画模拟[D];安徽大学;2004年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026