变形飞机动态气动特性数值模拟研究

根据近年来国内外提出的折叠翼变形飞机方案,设计了一种折叠翼变形飞机的3D数模,利用作者发展的动态混合网格技术和非定常计算方法,对该折叠翼变形飞机机翼折叠和展开过程中的动态气动特性进行了数值模拟,并与准定常计算结果进行了对比分析.计算结果表明:机翼在展开/折叠的过程中,飞机的气动特性会发生剧烈变化,在变形飞行控制率设计时需要引起高度的重视.     

HyperFLOW 软件数值模拟TrapWing 高升力外形北大核心CSCD

采用High-Lift研讨会提供的梯形翼外形(TrapWing),利用自主研发的基于结构、非结构网格的通用CFD软件HyperFLOW进行了数值计算,以评估其对复杂外形低速流的模拟能力。分别采用了三套不同拓扑结构的计算网格,包括两套非结构/混合网格和一套多块结构网格,每套网格又分为粗、中、细三种密度不同的网格数量以考察其网格收敛性。利用Richardson插值法,对计算结果开展了可信度分析。结果表明,不管是结构网格还是非结构网格,HyperFLOW均建立了可接受的可信度;对于高升力外形数值模拟,SA湍流模型要比SST湍流模型模拟的更准确;在失速迎角附近,现有的二阶精度解算器仍需持续改进。     

任意平面域的三角形网络和混合网络生成

本文以阵面推进法为基础,给出了一种可对任意平面域进行三角形网络和混合网络自动生成的方法。该方法包括了易于生成与调控的矩形背景网络技术,简便的外形输入和初始陈面划分技术,加速查寻的数据结构,以及节点松驰和基于Ｄｅｌａｕｎａｙ准则的对角线交换技术。网格生成的实例表明,方法能快速可靠地生成任意平面域的高质量用角形网格和适于粘性计算的三角形／四边形混合网络,具有较好的通用性与实用性。     

非结构网格及混合网格复杂无粘流场并行计算方法研究

在分布式存储网络异构编程PVM环境下,本文将原串行计算程序改造为并行计算程序,提高了计算效率.利用非结构网格,对球体及航天飞机外形绕流流场进行了计算,取得了较好的性能结果.最后将该方法应用于混合网格的计算,求解了多种复杂外形的无粘绕流,如翼身组合体和DDF1C导弹,得到了满意的结果.     

E 级计算给 CFD 带来的机遇与挑战

E 级(Exascale)计算机有望在2020年前后投入使用,E 级计算将给计算科学和科学研究带来革命性的变化。计算流体力学(CFD)作为超大规模计算机应用的重要领域之一,将迎来前所未有的发展机遇,同时也将面临极其严峻的技术挑战。本文对当前国内外超大规模 CFD 计算的现状进行了概述,探讨了未来 CFD 的发展趋势,并对 E 级计算给 CFD 带来的机遇与挑战进行了分析,最后提出了适应未来 E 级计算的 CFD 发展思路与建议。为了实现 E 级计算的宏伟目标,CFD 与计算机科学、应用数学等学科的“协同设计”势在必行。我们期待本文的分析能对我国的高性能计算应用、CFD 发展有一定的启示作用。     

基于动态混合网格的多体分离数值模拟方法

在以往发展的动态混合网格技术的基础上,耦合六自由度运动方程计算和非定常流场计算,建立了一种多体分离问题的非定常数值模拟方法。首先采用刚体运动的六自由度运动方程,建立刚体运动轨迹计算模块,然后以该模块为纽带将以往建立的非定常计算方法与改进的基于Delaunay背景网格插值方法和局部网格重构方法相结合的动态混合网格生成方法耦合起来,形成了精度较高、效率较快的多体分离数值模拟方法。利用该方法对典型多体分离问题进行了数值模拟,计算结果与风洞实验结果吻合较好。     

基于网格缩小的非结构网格梯度重构及制造解精度测试与验证

网格收敛性研究是验证与确认中的重要手段,本文针对一般非结构网格难以自相似加密的问题,实现了基于网格缩小的精度测试方法。在典型的各向同性和各向异性拉伸网格上,利用网格缩小精度测试方法分别考察了梯度重构精度以及制造解流动模拟精度,将网格缩小精度测试结果与传统的网格加密精度测试结果以及理论分析结果进行对比,验证了网格缩小精度测试方法与网格加密精度测试方法及理论分析的一致性。最后将网格缩小精度测试方法应用到各向异性拉伸网格粘性制造解精度测试中,得到了预期中的结果,并初步考察了梯度重构方法、网格类型对数值模拟精度的影响,显示出网格缩小精度测试方法在一般非结构网格精度测试上的优势,具有一定应用前景。     

海鸥翼折转运动的数值模拟及分析

扑翼的非定常流动及气动特性分析对微型飞行器设计具有重要意义。本文采用动态混合网格技术以及非定常数值计算方法,对海鸥翼的拍动过程进行了数值研究。采用基于径向基函数的插值技术实现网格的变形,为了提高插值效率,发展了基于最大物面误差的参考点选择算法。基于文献观测数据建立了海鸥翼的三维模型,并设计了相应的拍动以及变形规律,对拍动角、折转角的影响进行了分析。分析结果表明,折转角可以减小上拍过程的不利影响,对提高整个扑翼周期的时均升力、减少时均阻力和能耗是有益的。进一步通过调整折转的相位来增加翼的折转时间,并减少其展开的时间,对提高扑翼过程的平均气动力特性是有益的。     

鱼类“钻石”状群游流动干扰及其减阻机理的数值研究

观察研究发现,在鱼类的群游中,鱼体可以有效利用相互之间有益的干扰,减少游动阻力,增强自己的推力和游动效率.本文基于作者以往建立的动态混合网格生成方法以及不可压缩流非定常计算方法,对四鱼呈"钻石"阵型排列的群游进行了数值模拟,通过不同横向间距、不同流向间距、不同相位差等游动状态下流场结构、力学特性的比较,分析了鱼体之间的横向干扰和流向干扰,对其降低阻力、减少功耗的流动机理进行了探讨.     

基于Delaunay背景网格插值和局部网格重构的变形体动态混合网格生成技术

建立了一种基于Delaunay背景网格插值方法和局部网格重构方法相结合的变形体动态混合网格生成方法.首先利用作者发展的定态混合网格生成技术生成变形体初始网格,即变形体附近利用结构化的四边形网格,外场采用自适应Cartesian网格,中间由三角形网格过渡.当物体运动或变形时, 首先利用Liu和Qin提出的基于Delaunay背景网格的插值方法对变形体附近的网格进行变形;当物体运动或变形位移很大导致局部网格质量急剧下降、甚至网格相交时,则在局部重新生成网格.利用该方法生成了多种变形体外形的动态混合网格,如单个鱼体、双鱼串列巡游,昆虫单翼、双翼扑动,多段翼型变形等,并在鱼体巡游的非定常数值模拟中得到了应用.