收藏本站
《南京航空航天大学》 2011年
收藏 | 手机打开
二维码
手机客户端打开本文

高超声速多场耦合及其GPU计算加速技术研究

张兵  
【摘要】:高超声速飞行器,由于其飞行环境中存在激波边界层干扰、高温流动等复杂物理现象,给研制工作带来了一系列新的课题,尤其体现在气动热弹性等问题中,也是目前国际上的研究热点之一。解决这些问题的基础已经远远超出常规的气动弹性力学范畴,必须计入高温气体效应、边界层效应、结构热传导和热辐射等研究内容。它们具有流-固-热多场耦合和计算密集型特点,需要计算流体力学(CFD)、计算结构动力学(CSD)、计算传热学(CTD)等多个学科的相互配合,对理论分析和计算能力都是一个很大的挑战,目前已有的成功分析还较少。本文针对上述问题开展了以下研究工作: 针对化学反应非平衡流N-S方程,对隐式时间离散格式中的对流项、粘性项及源项的雅可比矩阵表达式进行了推导。采用Menter’SST两方程湍流模型、Roe格式、AUSM+-up格式及LUSGS方法编制了高超声速CFD计算程序,并通过典型算例验证其正确性,为后续研究工作提供基本计算手段。 从GPU架构特点出发,发展了基于数据并行的隐式CFD求解方法,适用于结构和非结构网格。采用CUDA技术,通过数据结构和算法优化,在GPU上成功实现了高超声速CFD计算。分别在Intel Core2Quad3.0GHz CPU和NVIDIA GTX280GPU上进行了数值仿真,结果表明隐式格式计算速度是显式格式6倍以上,采用显式格式的加速比达到28倍,采用隐式格式加速比达到了28.7倍,同时加速比随问题规模的增加而增加。仿真结果和实验值吻合良好。该工作可为后续研究提供快速计算手段。 依据几何关系及插值特点,提出了一种带旋转修正的TFI动网格方法,用来解决传统TFI方法在粘性网格和网格大变形时引起的网格正交性问题。以典型二维及三维粘性网格为例进行了方法的有效性研究。结果表明,在显著大变形情况下,引入旋转修正得到了正交性和光顺性良好的变形网格。方法的计算效率虽然较传统TFI有所降低,但相比弹簧方法能提高1到2个量级。 在位移和温度插值方法上,鉴于常规的基于薄板变形理论的样条方法(TPS/IPS)对于复杂结构外形存在插值困难,且对于温度插值没有明显的物理意义,提出了采用高阶等参单元形函数代替TPS/IPS的方法。算例表明,本文方法不仅成功解决了温度插值问题,而且位移插值精度也优于TPS/IPS方法,并更具普适性。此外,针对压力和热流插值需要保持守恒性的特点,从有限元法和有限体法的单元特性出发,提出一种具有局部守恒特性的界面载荷插值方法,并通过算例验证了方法的有效性。 采用共享内存技术开发了适用于通用有限元分析(FEA)和计算流体力学(CFD)软件的多场耦合计算平台,并基于分区耦合方式实现了流固耦合传热计算。作为算例,计算了外壁冷却的喷管和高超音速圆柱绕流的耦合传热问题,结果与实验值吻合良好。针对类X-34飞行器的头部热防护结构,考虑材料非线性和辐射效应,对高超音速巡航状态下驻点温度和结构冷却系统功率随热防护层厚度的变化规律进行了研究。计算结果表明,驻点温度随厚度的变化并不明显,而冷却系统功率随厚度增加急剧降低。此外材料发射率非线性对结果影响较大。 采用Euler方程和N-S方程研究了超声速和高超声速壁板颤振中的湍流边界层效应。低超声速条件下的计算结果表明,湍流边界层对颤振边界影响较为明显,在马赫数1.2左右达到最大值,随后这种差别随马赫数增大而逐渐减小,颤振边界计算结果和试验值吻合较好。在高超声速阶段,湍流边界层效应对颤振动压仍有较明显影响。在马赫数8时,N-S方程的结果高出Euler方程20%左右。这说明,对于高超声速壁板颤振,湍流边界层效应是不可忽略的影响因素。 采用三阶活塞理论、Euler方程和N-S方程三种气动力模型对二元双楔翼型的气动弹性问题进行了研究。结果表明采用Euler和三阶活塞理论的颤振速度非常接近,但和N-S方程结果差别较大。这可以归结为两方面原因:一是激波边界层干扰引起流场特征发生变化,二是边界层厚度间接改变了结构外形。因此,粘性效应对此类高超声速气弹稳定性有显著影响。 对典型吸气式高超声速飞行器机翼的三维气动热弹性问题进行了研究。将其分解为静气动热弹性配平和气动热弹性响应问题。其中,对于静气动热弹性配平问题,为避免计算过程不收敛,提出了采用动态过程代替静态过程的计算方法。在海拔10公里标准大气条件下,对一典型非对称高超声速机翼结构进行了气动热弹性计算,其中,气动力模型采用理想气体和化学反应非平衡气体。为进行比较,还采用活塞理论和Euler方程计算了相应的颤振速度。结果表明,采用理想气体时,用活塞理论和Euler方程得到的颤振速度比N-S方程结果分别高164.3和98.7%。而同样在N-S方程中,采用理想气体比化学反应非平衡气体的颤振速度低6.1%。造成上述结果的主要原因是高温引起的热应力及材料特性变化改变了结构的动力学特性。对于这类气动热弹性问题需要采用包含化学反应非平衡效应的N-S方程来求解。
【学位授予单位】:南京航空航天大学
【学位级别】:博士
【学位授予年份】:2011
【分类号】:TP391.41;O354

手机知网App
【参考文献】
中国期刊全文数据库 前10条
1 刘学强;李青;柴建忠;N.Qin;;一种新的动网格方法及其应用[J];航空学报;2008年04期
2 杨超;许赟;谢长川;;高超声速飞行器气动弹性力学研究综述[J];航空学报;2010年01期
3 陈飞国;葛蔚;李静海;;复杂多相流动分子动力学模拟在GPU上的实现[J];中国科学(B辑:化学);2008年12期
4 周季夫;钟诚文;尹世群;解建飞;张勇;;基于GPGPU的Lattice-Boltzmann数值模拟算法[J];计算机辅助设计与图形学学报;2008年07期
5 张涵信;差分计算中激波上、下游解出现波动的探讨[J];空气动力学学报;1984年01期
6 张涵信;无波动、无自由参数的耗散差分格式[J];空气动力学学报;1988年02期
7 段占元,童秉纲,姜贵庆;高分辨率有限差分-有限元混合方法及其在气动热计算中的应用[J];空气动力学学报;1997年04期
8 段占元,童秉纲,姜贵庆;三维物体表面摩阻的有限元算法[J];空气动力学学报;2000年01期
9 梁德旺,王可;AUSM+格式的改进[J];空气动力学学报;2004年04期
10 杨亚政;李松年;杨嘉陵;;高超音速飞行器及其关键技术简论[J];力学进展;2007年04期
中国博士学位论文全文数据库 前4条
1 柳军;热化学非平衡流及其辐射现象的实验和数值计算研究[D];国防科学技术大学;2004年
2 李海燕;高超声速高温气体流场的数值模拟[D];中国空气动力研究与发展中心;2007年
3 贺旭照;高超声速飞行器气动力气动热数值模拟和超声速流动的区域推进求解[D];中国空气动力研究与发展中心;2007年
4 潘沙;高超声速气动热数值模拟方法及大规模并行计算研究[D];国防科学技术大学;2010年
【共引文献】
中国期刊全文数据库 前10条
1 王保国;郭延虎;刘秋生;沈孟育;;HIGH-ORDER ACCURATE AND HIGH-RESOLUTION UPWIND FINITE VOLUME SCHEME FOR SOLVING EULER/REYNOLDS-AVERAGED NAVIER-STOKES EQUATIONS[J];Acta Mechanica Sinica;1998年01期
2 王汝权;申义庆;;SOME WEIGHT-TYPE HIGH-RESOLUTION DIFFERENCE SCHEMES AND THEIR APPLICATIONS[J];Acta Mechanica Sinica;1999年04期
3 姜宗林;;ON DISPERSION-CONTROLLED PRINCIPLES FOR NON-OSCILLATORY SHOCK-CAPTURING SCHEMES[J];Acta Mechanica Sinica;2004年01期
4 ;Non-linear characteristics of Rayleigh-Taylor instable perturbations[J];Acta Mechanica Sinica;2008年02期
5 ;Numerical investigation on the flowfield of “swallowtail” cavity for supersonic mixing enhancement[J];Acta Mechanica Sinica;2009年01期
6 吴开腾;张莉;;一类三维基于通量分裂的高精度MmB差分格式[J];兵工学报;2010年S1期
7 李椿萱,戴吾伦;绕AOTV高超音速流场的数值计算[J];北京航空航天大学学报;1993年01期
8 王浩,吴颂平,陆烨;高超音速化学非平衡流的隐式格式迭代解法[J];北京航空航天大学学报;2002年01期
9 戴梧叶,刘宇,马彬;塞式喷管三维流场的数值模拟[J];北京航空航天大学学报;2002年05期
10 童晓艳,梁国柱,张国舟,安联;气动谐振加热效应的数值模拟[J];北京航空航天大学学报;2004年12期
中国重要会议论文全文数据库 前10条
1 李微;李跃明;;温度梯度下典型热防护结构声-振耦合响应分析[A];全国结构振动与动力学学术研讨会论文集[C];2011年
2 刘超峰;李海东;杨炳渊;;身-翼-舵组合体高超声速变攻角颤振计算[A];第十届全国振动理论及应用学术会议论文集(2011)下册[C];2011年
3 吴德全;王健平;;有限谱NND格式计算可压缩流[A];计算力学研究与进展——中国力学学会青年工作委员会第三届学术年会论文集[C];1999年
4 杨科之;杨秀敏;董军;邓国强;;抗内爆炸结构的等效荷载设计法[A];第十二届全国结构工程学术会议论文集第Ⅱ册[C];2003年
5 宋书恒;张树海;朱国林;;翼型非定常运动绕流的数值模拟[A];首届全国航空航天领域中的力学问题学术研讨会论文集(上册)[C];2004年
6 杨炳渊;史哓鸣;梁强;;高超声速有翼导弹多场耦合动力学的研究和进展[A];第十届全国空气弹性学术交流会会议论文集[C];2007年
7 潘沙;丁国昊;王文龙;刘枫;李桦;;一维激波管问题差分格式计算对比分析[A];第十三届全国激波与激波管学术会议论文集[C];2008年
8 李旭东;姜宗林;;散心爆轰波在球形固壁上的反射[A];第十三届全国激波与激波管学术会议论文集[C];2008年
9 刘凯礼;李博;雷雨冰;张堃元;;宽高比对二元高超声速进气道性能的影响研究[A];第一届高超声速科技学术会议论文集[C];2008年
10 张兵;韩景龙;;一种GPU上的CFD并行计算方法[A];第十一届全国空气弹性学术交流会会议论文集[C];2009年
中国博士学位论文全文数据库 前10条
1 陈艺冰;多介质流体力学计算的GKS方法研究[D];中国工程物理研究院;2010年
2 闫洁;三维气固两相平板间射流中相间耦合作用的直接数值模拟[D];浙江大学;2009年
3 吴锋;基于GPU并行计算的数值模拟与燃煤锅炉系统的优化研究[D];浙江大学;2010年
4 冯定华;高速流动精细数值模拟、实验研究及其在气动光学中的应用[D];国防科学技术大学;2010年
5 李金梁;箔条干扰的特性与雷达抗箔条技术研究[D];国防科学技术大学;2010年
6 徐江斌;基于气象数据的云景真实感模拟技术研究[D];国防科学技术大学;2010年
7 徐晓亮;热防护机理与烧蚀钝体绕流的涡方法研究[D];北京交通大学;2011年
8 张磊;离心风机多目标优化及旋转失速机理研究[D];华北电力大学(北京);2011年
9 汪利;可压缩湍流边界层统计特性的直接数值模拟研究[D];中国科学技术大学;2011年
10 王筱蓉;冲压增程炮弹进气道型面气动优化方法研究[D];南京理工大学;2010年
中国硕士学位论文全文数据库 前10条
1 朱缓;基于分形算法的海面动态仿真[D];哈尔滨理工大学;2010年
2 赵金山;流场/结构耦合传热算法在高超声速风洞结构传热分析中的应用[D];中国空气动力研究与发展中心;2010年
3 刘伟;基于混合网格的高阶间断Galerkin/有限体积混合算法的研究[D];中国空气动力研究与发展中心;2010年
4 罗万清;气流焓值对非烧蚀材料表面氧化反应影响的数值研究[D];中国空气动力研究与发展中心;2010年
5 李刘苗;海天背景目标红外辐射特性研究[D];西安电子科技大学;2011年
6 雷军命;耐高过载、耐高弹速涡轮发电机研究[D];西安电子科技大学;2010年
7 刘沛;二维DFT并行算法在CUDA架构上的实现[D];暨南大学;2011年
8 滕人达;基于ABEEMσπ模型能量求解并行处理的研究[D];辽宁师范大学;2011年
9 吴颀峰;基于GPU的数字图像处理算法研究[D];江西理工大学;2011年
10 郭俊;新型低密度C/C隔热材料的制备研究[D];中南大学;2009年
【二级参考文献】
中国期刊全文数据库 前10条
1 范晓樯,李桦,丁猛;钝头双锥体有攻角表面热流计算[J];兵工学报;2002年01期
2 吕晓斌,朱自强;二维跨音速Euler方程分区并行算法[J];北京航空航天大学学报;2000年02期
3 李君哲,阎超,柯伦,张书庭;气动热CFD计算的格式效应研究[J];北京航空航天大学学报;2003年11期
4 马辉,赵烈,王发民;高超声速三维化学非平衡流动的数值模拟[J];北京航空航天大学学报;2004年02期
5 马辉,王发民;高温化学非平衡计算中源项处理方法[J];北京航空航天大学学报;2004年07期
6 吴志刚,惠俊鹏,杨超;高超声速下翼面的热颤振工程分析[J];北京航空航天大学学报;2005年03期
7 刘宇,张振鹏;三维抛物化 N-S 方程全隐格式矢通量分裂求解[J];北京航空航天大学学报;1998年06期
8 孟凡颢,钟腾育;用壁板颤振理论解决某系列飞机的方向舵蒙皮裂纹故障[J];飞机设计;2000年04期
9 宋健;在推力和阻力作用下飞行器横向振动分析[J];中国工程科学;2000年10期
10 黄玉盈,钟伟芳,朱达善,李四平,王和慧;各向异性中厚度板壳的弹塑性大变形分析[J];工程力学;1995年03期
中国重要会议论文全文数据库 前1条
1 曹平宽;;万亿次规模航空CFD并行算法技术初探[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
中国博士学位论文全文数据库 前9条
1 郑忠华;双模态冲压发动机燃烧室流场的大规模并行计算及试验验证[D];中国人民解放军国防科学技术大学;2003年
2 侯中喜;超声速复杂流场并行数值分析及高阶格式研究[D];中国人民解放军国防科学技术大学;2000年
3 柳军;热化学非平衡流及其辐射现象的实验和数值计算研究[D];国防科学技术大学;2004年
4 朱君;流体力学数值方法及并行策略研究[D];南京航空航天大学;2006年
5 李海燕;高超声速高温气体流场的数值模拟[D];中国空气动力研究与发展中心;2007年
6 王兰;超燃冲压发动机整机非结构网格并行数值模拟研究[D];中国空气动力研究与发展中心;2007年
7 贺旭照;高超声速飞行器气动力气动热数值模拟和超声速流动的区域推进求解[D];中国空气动力研究与发展中心;2007年
8 王学德;高超声速稀薄气流非结构网格DSMC及并行算法研究[D];南京航空航天大学;2006年
9 李晓宇;高超声速飞行器一体化布局气动外形设计[D];国防科学技术大学;2007年
中国硕士学位论文全文数据库 前6条
1 李凰立;再入弹头的气动加热及热响应分析[D];西北工业大学;2001年
2 王可;AUSM+格式的研究和应用[D];南京航空航天大学;2002年
3 张江涛;CFD显式差分计算程序的自动并行技术研究[D];国防科学技术大学;2002年
4 于战华;跨声速欧拉方程并行算法研究及应用[D];南京航空航天大学;2004年
5 吕丽丽;高超声速气动热工程算法研究[D];西北工业大学;2005年
6 余奇华;三维高超声速绕流流场并行算法研究[D];南京航空航天大学;2006年
【相似文献】
中国期刊全文数据库 前10条
1 沈赤;王沛明;许志宏;;优势区相图三维图形log P_i-log P_j-T的绘制[J];计算机与应用化学;1993年01期
2 吴艳花;陈洁群;;基于GIS的地价管理应用信息系统解决方案[J];沿海企业与科技;2009年07期
3 张栒;湛成伟;邓辉文;唐甘翌;;一种快速减法聚类算法[J];西南师范大学学报(自然科学版);2006年03期
4 来元,冯焕清,诸葛斌;一种用于表面重建的网格点生成算法[J];计算机工程与应用;2002年03期
5 王梅珍;王府梅;王善元;;基于网格的织物虚拟场景真实感模拟及计算机实现[J];东华大学学报(自然科学版);2006年06期
6 江志国;车竞;;基于网格面元法的高超声速飞行器性能分析[J];计算机仿真;2011年03期
7 欧莽;钦明皖;汪继文;;基于非结构三角网格的两种自适应技术研究[J];电脑知识与技术(学术交流);2006年35期
8 刘辉;彭滢;龚斌;代斌;魏代政;;网格环境下期权定价BSDE模型的并行实现[J];华中科技大学学报(自然科学版);2011年S1期
9 郭福成;;基于多级网格搜索的固定单站无源定位方法[J];信号处理;2008年06期
10 谢丹;;Illustrator CS3表现产品的质感[J];电脑知识与技术;2011年19期
中国重要会议论文全文数据库 前10条
1 杨新颖;龚光红;王行仁;;CTDB地形数据库的多维信息可视化[A];'2008系统仿真技术及其应用学术会议论文集[C];2008年
2 王元清;曾文平;朱小羽;贺东哲;;复杂空间管节点有限元网格的自动剖分算法[A];第三届全国现代结构工程学术研讨会论文集[C];2003年
3 李艳;李守先;余真;;Vicon程序模拟三维喷管的可行性探讨[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
4 梁树献;刘小虎;程兴无;杨亚群;徐珉;;距离加权平均法计算流域平均面雨量[A];新世纪气象科技创新与大气科学发展——中国气象学会2003年年会“03.7淮河大水的水文气象学问题”分会论文集[C];2003年
5 刘静;;用时间动态小网格订正法进行宁夏扬黄新灌区温度及热量资源网格点推算[A];西部大开发 科教先行与可持续发展——中国科协2000年学术年会文集[C];2000年
6 姚彦忠;袁光伟;倪国喜;曹轶;;一种变分方式的平面网格生成与网格重分方法[A];中国工程物理研究院科技年报(2005)[C];2005年
7 张雪芬;陈东;邹春辉;陈怀亮;;格点化冬小麦发育期推算方法[A];中国气象学会2005年年会论文集[C];2005年
8 邵爱梅;邱崇践;;4DVAR中插值计算引发的不稳定性研究[A];中国气象学会2005年年会论文集[C];2005年
9 宋国杰;杨顶辉;陈亚丽;;求解声波和弹性波方程的简化近似解析离散化方法[A];中国地球物理学会第二十四届年会论文集[C];2008年
10 何长江;刘军;梁仙红;;三维弹塑性流体力学自适应欧拉方法研究[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
中国重要报纸全文数据库 前10条
1 ;网格点亮澳网[N];网络世界;2005年
2 记者 陈湘萍;掌控地价走势 预测未来地价[N];南宁日报;2008年
3 本报记者 黄穗诚;广州公布土地使用权基准地价[N];广东建设报;2006年
4 易芳;规划总图PPS介绍[N];建筑报;2000年
5 张睿;穗首批限户型地块今日揭标[N];南方日报;2006年
6 张睿;先查基准地价 后计补交地价[N];南方日报;2007年
7 本报记者 黄穗诚;广州规范土地有偿使用[N];广东建设报;2005年
8 卢风华;广州土地有偿使用出新标准[N];中国国土资源报;2005年
9 田桂丹;二手非商品房交易补地价[N];房地产时报;2007年
10 穗房宣、边防军;广州基准地价实现“大一统”[N];中国房地产报;2004年
中国博士学位论文全文数据库 前10条
1 査文舒;基于PEBI网格的油藏数值计算及其实现[D];中国科学技术大学;2009年
2 张学莹;多介质流体运动界面高分辨率Level-Set方法研究[D];南京航空航天大学;2006年
3 徐小文;可扩展并行代数多重网格算法研究[D];中国工程物理研究院;2007年
4 黄玲;加权本质非振荡格式和快速扫描法和在行人流模型中的应用[D];中国科学技术大学;2008年
5 刘鑫;面向化学非平衡流的CFD并行计算技术和大规模并行计算平台研究[D];解放军信息工程大学;2006年
6 刘新儒;带边界约束的网格和曲面生成理论与方法研究[D];中南大学;2011年
7 王森铭;飞秒强激光场中双原子分子电离与解离动力学研究[D];大连理工大学;2006年
8 孙江龙;绕潜艇三维厚边界层计算和二维自由面波动研究[D];华中科技大学;2006年
9 姚彦忠;流体力学计算中网格优化方法研究[D];中国工程物理研究院;2007年
10 杨建中;复杂多曲面数控加工刀具轨迹生成方法研究[D];华中科技大学;2007年
中国硕士学位论文全文数据库 前10条
1 王伟;地学应用中的网格数据可视化[D];吉林大学;2013年
2 吕珍;基于OpenGL的流场数据可视化技术[D];武汉理工大学;2007年
3 杜宏云;测地自由曲线及其性质研究[D];南京航空航天大学;2008年
4 尹清波;基于激光网格标记的视觉坐标测量技术的研究[D];青岛大学;2004年
5 周义勇;脉搏图像化检测方法研究[D];兰州理工大学;2006年
6 左莹;基于高阶交错网格的有限差分地震波场数值模拟[D];长安大学;2009年
7 吴贺贺;基于非结构动网格技术及运动物面粘流场数值计算方法研究[D];江苏科技大学;2012年
8 党占威;基于空间剖分的网格特征恢复方法[D];浙江大学;2013年
9 张永华;叶轮机CFD网格生成[D];南京航空航天大学;2007年
10 石磊;基于Richardson外推法的CFD网格误差研究[D];华北电力大学(北京);2011年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026