网格自适应与并行计算在气动力计算中的应用

来源 :南京航空航天大学 | 被引量 : 6次 | 上传用户:Holden
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
计算的准确性与效率一直是计算流体力学研究的重点之一。本文从网格自适应与并行计算两方面对计算的准确性与效率进行了研究。针对自适应方法中的关键问题,提出了新的捕捉流场中关键区域的方法;针对现存串行程序与硬件条件,发展了一种非结构并行计算算法,并结合自适应技术发展了一套并行自适应方法。基于非结构有限体积法,发展了一套适用于不可压/可压、定常/非定常、无粘/有粘、二维/三维流场的数值模拟求解器,数值计算了二维翼型、三维机翼与翼身组合体的气动力。空间离散采用二阶精度的Jameson格心格式。时间离散,对定常流动采用五步Runge-Kutta与LU-SGS方法;对非定常流动采用双时间步长法。对流通量采用JST格式。湍流模型选取一方程SA模型。合理运用了一些加速收敛措施,并针对不可压流引入低速预处理加速收敛技术。最后,通过算例对准确性、收敛性进行了验证。为了提高计算准确性、减少计算资源,采用了h型网格自适应技术。首先给出二维四边形、三维六面体网格为初始稀疏网进行流场的计算。在初始网格流场解的基础上,本文提出速度紊乱度构造指示函数对粘性流动中的附面层、湍流区与尾迹区进行捕捉并局部加密网格。为了简洁高效地捕捉激波与涡结构区域,文中提出以单位体积总能为指示变量构造统一的指示函数捕捉上述区域。对于局部网格的加密,二维情况下,一个四边形网格单元将被切分为四个小的四边形网格单元;三维情况下,一个六面体网格单元将被切分为八个小的六面体网格单元。但是对二维/三维物面第一层网格加密时,将不对法向方向进行切分。自适应后的网格不同层级网格间会出现悬挂点。文中采用按边循环的非结构求解器将解决悬挂点问题,同时保证了边/面上的通量守恒。新产生的网格单元处的流场值从初始流场解插值获得。这样自适应后的流场计算可在初始流场解的基础上进行,节约了计算时间。针对三维数值模拟中网格数量巨大的特点与目前的计算机硬件条件,流场迭代计算中使用了并行计算方法。本文提出只在流场迭代部分采取并行措施,其余部分全部由0进程完成;初始网格与自适应后网格的划分采用在程序中调用METIS的方法;并行计算时,进程间信息的交换采用MPICH2函数接口。重点研究了分区边界处理策略,数据高效交换及快速的输入输出问题。METIS网格划分软件开源且易于获取,针对不同的网格数提供不同的调用接口,具有分区质量高、速度快、最小化子域边界数量的优点。MPICH2是目前主流的并行函数接口标准。因此上述方法对现有串行程序的改造简单快速,适用性广。由于自适应前后都保证了负载的平衡,此方法具有较高的并行效率。通过算例对并行自适应程序进行了验证。
其他文献
随着信息技术和知识经济的发展,现代服务业的内涵发生了变化,当今世界已进入服务经济时代,尤其是随着现代科学技术向服务业的加速渗透,使之在现代经济发展过程中占有越来越重
由于信息技术在教学中的重要作用,在学科教学中的普遍应用信息技术,实现信息技术与学科教学的融合已成为发展的必然趋势。因此,笔者对我校信息技术与学科教学整合进行了详细
劳动法学理论和司法实践中,集体谈判成为解决集体劳资纠纷的主要途径。随着富士康、本田事件的发生,劳资集体谈判越来越成为社会热点,政府和工会正在推动企业与工会之间的工
针对河北省中小学校突发公共事件展开研究,阐述中小学校突发公共事件的类型以及存在的问题。探究体育教育在中小学校突发事件应急管理中的作用:突发公共事件发生前的预防和导
针对吊管机运行中存在的问题,如履带式车辆行驶在纵向坡面,会出现倾覆危险、维修成本高等,进行全方位分析,并简单介绍了改进吊管机结构的必要性,提出吊管机结构改进要点与施
介绍酸溶分解法、焦硫酸钾分解法、碱性熔融法、微波消解法、离心浇铸等试样前处理技术在分析Cr系铁合金元素时的应用特点及效果。
对三相异步电动机进行运行维护,是保证其基本运行效率的重要环节。围绕三相异步电动机,重点分析了其运行维护的几点措施,以期为设备运行的有效性提供支持。
前几天发生的一起伤人事件,把人们的视线重新引到淡出公众视野很久的上海家化身上。9月18日傍晚,上海家化董事长兼总经理谢文坚在公司办公楼前遇袭,所幸伤势并不严重。长期在外
会议
介绍了目前用于我国森林火灾监测的美国NOAA系列气象卫星和EOS卫星.NOAA/AVHRR气象卫星信息源具有周期短、密度高、多时相动态遥感的特殊能力,但分辨率相对较低.EOS卫星在林
选用紫丁香(SyringAoblAtA)为试验材料,对紫丁香的实生苗进行4种处理:CK(土壤含水量为试验用土最大田间持水量的75%)、W1(土壤含水量为试验用土最大田间持水量的60%)、W2(土