基于Navier-Stokes方程的复杂流动数值模拟精度与并行计算研究

来源 :西安电子科技大学 | 被引量 : 0次 | 上传用户:trittt
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
随着计算机和计算技术的发展,计算流体力学(CFD)在工程应用中取得了长足的进步,在航空、航天、石油、能源等工业领域发挥着巨大的作用.如今,CFD技术已经成为航空、航天飞行器气动设计中的必要手段.但是CFD还存在一些具有挑战性的问题.当CFD用于解决复杂流动问题时,其精度和效率仍然是CFD研究人员面临的难题.本文主要围绕CFD的计算精度和效率问题开展了研究,并进行了实际外形的空-空导弹复杂流动数值模拟.本文主要工作概括如下:1.基于超立方体概念,提出了一种通用的绕翼身组合体外形的高质量多块拼接结构化网格构建方法.用这种方法生成了绕DLR-F4翼身组合体外形的三套不同密度的“超立方体”计算网格.通过求解雷诺平均Navier-Stokes (RANS)方程耦合Baldwin-Lomax (BL)代数模型和Spalart-Allmaras (SA)一方程模型,数值模拟了DLR-F4翼身组合体流场.在“超立方体”网格上获得了与实验数据吻合良好,且优于其他软件和其他网格的计算结果.研究表明:本文高质量网格构建方法可行,用此方法生成的网格质量高.这种网格构建方法适用于任意具有相似外形的翼身组合体网格生成.2.选用AIAA第二界阻力研讨会(DPWⅡ)提供的粗、中、细不同密度的多块结构对接网格,通过求解雷诺平均Navier-Stokes方程耦合Spalart-Allmaras(S A)、Wilcox’s k-ω(Wilcox)和Menter’s k-co SST (SST)三种湍流模型,数值模拟了DLR-F6翼/身(WB)和翼/身/挂/舱(WBNP)两种外形的流场.研究了阻力计算精度、考察了网格和湍流模型对组合体构型总体气动特性的影响.计算采用Roe格式以准确地捕捉激波.为提高格式的精度、消除激波附近的非物理振荡采用了MUSCL (Monotone Upstream centered Scheme for Conservation Laws)插值和三阶迎风偏置限制器.通过一系列精心设计的测试算例,基于SA湍流模型考察了CFD对网格尺度的敏感性,评估了三种不同湍流模型的性能.取得了与实验符合良好的计算结果.SA模型预测的翼身组合体的总体阻力优于其他软件的结果.与实验对比发现:使用SA湍流模型,网格密度对压差阻力影响较大,对机翼和吊舱表面的压力系数分布影响较小;湍流模型对机翼表面的激波位置有影响,对阻力尤其是摩擦阻力影响较大,对翼根处的流动分离有影响,SA模型预测的分离区最大,而Wilcox预测的最小.总体来说,SST模型的计算结果与实验数据最接近,预测的阻力一致小于SA和Wilcox模型的阻力;SA模型的升力与实验吻合较好;Wilcox模型高估了升力和阻力.计算还发现,与实验相同攻角的条件下运行CFD程序,获得的机翼表面压力系数分布与实验吻合较好,而与之相匹配升力条件下的结果不然.3.通过求解多个标准模型的流动问题,系统地分析研究了低耗散E-CUSP (LDE)迎风格式在可压缩流动中的性能.对比格式选用广泛使用的Roe格式,控制方程为Euler或雷诺平均Navier-Stokes方程耦合Spalart-Allmaras (SA)湍流模型,采用全隐式无分裂Gauss-Seidel线松驰时间推进方法.计算表明:LDE格式能够准确地捕捉陡峭激波,计算结果与实验数据吻合良好.相对于Roe格式,LDE格式由于不需要矩阵运算,所以计算量小、收敛速度快.进一步使用LDE格式耦合SA模型数值模拟了绕细长旋成体超声速、大攻角分离流动.准确地模拟了背风面的流动分离和主分离涡的发展过程;捕捉到了精细的二次涡和具有相当分辨率的横向分离.预测到的力和力矩与实验吻合良好.研究表明,LDE格式耦合SA模型具有高的模拟精度和计算效率、鲁棒性好,可用于模拟复杂的分离流动.4.通过求解雷诺平均Navier-Stokes方程耦合Spalart-Allmaras湍流模型数值模拟了具有复杂外形的多翼空-空导弹复杂流场.计算方法为有限体积多步Runge-Kutta时间推进方法.用求解椭圆型偏微分方程并结合代数方法生成了连续拼接多块结构化网格.对具有滚转角或无滚转角的多翼空-空导弹的气动力和力矩进行了计算,获得了与实验数据吻合良好且优于其它程序的计算结果.5.提出了拆分方式有限制的目标为极小化最大完工时间的平行机排序问题.与以往研究不同,对工件的拆分方式进行了限制,该问题是NP难的.发展了求解该问题的一个近似算法,并用此算法解决了异构并行机系统并行计算中的负载平衡问题.通过实验仿真,并与传统方法相比较,验证了这个算法的可行性和有效性.将上述排序问题的研究成果应用于基于分区策略的CFD并行计算,针对异构并行环境,给出了一种负载平衡算法,解决了异构并行机系统CFD并行计算中的数据划分和任务分配问题.对处理机之间的通信方式进行了优化,减少了处理机之间的通信次数和进程之间的相互等待时间,提高了计算/通信比.在工作站集群环境下实现了并行计算,并用等效率标准和等速度标准对并行算法的可扩展性进行了测评.结果表明,本文并行算法具有良好的可扩展性,发展的负载平衡算法性能良好,并行效率高.本文发展的负载平衡算法适用于MIMD结构计算机上基于分区策略并行计算中的负载平衡问题.
其他文献
研究职工心理健康状况的评价与分析问题。首先运用MATLAB编程统计各类因素下每个人为自己评5、4、3、2、1分数项目的个数。然后运用层次分析法建立评价模型,确定各个项目对心
《水浒》女性人物空间存在呈现出独特的间断展现程式,设置程式手法主要有人物排布均匀化和人物配置单元化。这种手法完善了豪杰人物形象的塑造,凸显了女性人物存在的社会价值
文章运用文献资料、问卷调查、数理分析等研究方法,对岭南地域非物质文化遗产传承与高校武术教学融合现状进行研究,研究发现:"非遗"知识教育宣传工作不到位,学生不了解"非遗"
目的:本课题主要研究女性肥胖者腹腔内脂肪与中医证素的相关性关系,探寻腹腔内脂肪面积差异的关键性相关证素,为肥胖的既病防变提供参考依据。方法:采集150例成年女性肥胖者
司法警察职业道德的特征:1.与人民警察职业道德具有职业性质上的共同性;2.与法官职业道德在职业实践活动中的相关性;3.在具体内容上具有法定性、操作性和实践性。司法警察职
从我国A股上市公司近年来股份回购实例的统计出发,利用事件研究法建立模型,重点对云天化、长春高新和九芝堂等公司股份回购后的市场反应和公司价值进行实证分析,总结中国上市
令日本人一直不断憧憬的对象,没有什么能比得上富士山了。对于“富士山乃日本的第一山”的说法,大概没什么人对此持有异议。但是日本的富士山形象到底是如何成立的?其形成要
文章采用文献资料法、逻辑分析法等方法,结合青少年的成长发育对体育游戏的本质及其价值进行了探索性的研究,指出了体育游戏的本质属性有:自由性、无效用性、娱乐性、趣味性
<正>韩国首尔有9条地铁,乘坐地铁成了市民的首选,因为坐地铁实在太方便了,只要"见色行事"。打开首尔地铁交通图,五颜六色,仔细一看,却十分有序。9种颜色分别代表9条地铁线路:
中国特色案例制度与中国法文化之契合性可以分别从以例辅律之传统、比附类推之方法、人为称首之取向、案例中国特色等四个层面予以阐释。在本体论上,基于以例辅律之传统例律