CFD程序验证的虚构解方法及其边界精度匹配问题

采用网格收敛分析和虚构解理论对程序验证进行了理论探索与数值实验研究,并发展完善了该方法.给出了虚构解程序验证中一种比较简单的源项处理方案;给出了虚构解程序验证中气动方程组一种通用的虚构解系数选取方法,最后用虚构解方法对二维结构网格Euler流动计算程序ARC2D进行了验证,并发现了该程序边界条件实现精度与流场内部精度不相匹配问题,并给予了修正,数值实践证明虚构解方法是一种严格稳健有效的程序验证方法,用它完成的程序验证具有较高的可信度.     

基于网格划分的ADS-B地面站信号覆盖及选址分析

针对通用航空飞行活动的安全监视问题,合理的广播式自动相关监视（Automatic dependent surveillance-broadcast,ADS-B）地面站选址规划能够提高空域的监视范围,基于网格划分分析ADS-B地面站对真实高度的监视信号覆盖范围并进行选址分析。考虑地球曲率和地形遮蔽对信号传播距离的影响,采用数字地理高程数据建立地理模型,将空域网格化处理,用Xdraw算法计算网格覆盖性,确定ADS-B地面站对真实高度的最大理论覆盖范围,构建地面站最大覆盖模型进行空间规划布局。最后,对ADS-B地面站信号覆盖情况进行仿真,并以湖北省为例,计算出满足该区域覆盖需求的地面站布局。为实际工作中的地面站布局规划提供了依据。     

阵风响应及减缓的非定常数值模拟

基于网格速度和动网格方法,通过求解非定常N-S方程计算翼型副翼偏转时的阵风减缓响应,研究开环控制律参数和副翼缝隙对阵风减缓效果的影响。对于单个1-cos阵风,随着偏转角度的增大,减缓效果变大;由于粘性的影响,CFD计算得到的升力峰值减小量明显呈非线性的趋势,减缓效果低于准定常线性理论。对于连续sin阵风,通过调整副翼偏转的频率和最大偏转角,升力能达到很好的减缓效果,此时,副翼偏转频率应接近阵风频率,最大偏转角存在一个最优值。在阵风环境下,副翼间隙对升力系数影响不大,而力矩系数增大一倍,间隙会降低副翼偏转对升力的减缓效果。     

并行多重网格算法求解跨声速Euler方程

在工作站网络分布式存贮ＰＶＭ环境下组织了有效求解三维跨声速Ｅｕｌｅｒ方程多重网格算法的并行计算,讨论了如何将串行算法转变为并行算法,并有效提高并行效率的一些基本方法,最后分析了影响并行效率的一些瓶颈问题。     

基于CFD数值模拟的拉瓦尔喷管流场分析

结合一维等熵管流理论,采用CFD方法对某型拉瓦尔喷管典型工作的流场进行整场数值模拟,阐述了网格生成,流场初始化,以及RNG湍流模型使用方法,并对拉瓦尔喷管工作特性进行了气动分析。研究表明,CFD计算结合一维管流理论能明显提高计算效率,数值计算结果与理论相符较好。     

临近空间桨梢小翼螺旋桨布局气动增效研究

基于雷诺平均Navier—Stokes控制方程和非结构网格旋彬静止滑移面技术,开展了临近空间螺旋桨桨梢小翼构型气动布局与增效作用机理研究,采用静止与旋转域的多块网格布局策略及转／静滑移面技术,发展了高效的螺旋桨流场与性能分析方法;通过开展上翘小翼、上下翘小翼等多种桨梢小翼螺旋桨布局研究,初步建立了基础参数的性能影响规律,研究表明：桨梢小翼螺旋桨布局是临近空间工况环境提高气动效率的一种合理可行的技术途径。     

并行环境下外挂物动态分离过程的数值模拟

 并行工程的发展提出了计算流体力学的发展新概念,在此基础上研究开发了一个以非定常流计算为核心的多体分离数值仿真集成系统。该系统采用运动嵌套网格方法,能够处理飞行器在空间中的大尺度运动,结合六自由度运动学方程,具有实时仿真运动部件在气动力、弹射力和重力等耦合作用下运动的能力。基于外挂物投放数值模拟具有网格运动、非定常的特点,从网格重建、多重网格方法等方面探索了提高计算效率的途径和实现方法,结合并行化计算使分析周期大大缩短。以20个中央处理器(CPU)为例,并行计算和3层多重网格使计算速度提高约40倍,同时分布式存储、并发计算的特点缓解了对单节点的内存需求,解决了工程应用中受计算周期、网格规模等限制的难题。为了验证本文方法的正确性,对外挂物投放标模实验进行了数值仿真,计算得到的外挂物运动轨迹在0.5 s时间范围内与实验值基本吻合;对分离过程后半段出现的偏差,从气动力变化的时间历程上进行了说明,提出偏航力矩的误差和累积效应导致了预测出现偏离,对关键的分离初始阶段,该方法可以取得令人满意的效率和精度。     

先进直升机旋翼悬停状态气动性能计算

 为更好地模拟先进直升机旋翼流场,同时准确计算其悬停状态的气动性能,建立了一套基于Navier-Stokes/Euler方程的混合CFD方法。该方法的求解域由两部分组成：一是围绕旋翼桨叶周围的黏性区域,采用可压缩Navier-Stokes方程来模拟旋翼附近的黏性流动和近场尾涡的捕捉;二是离桨叶较远、黏性可以忽略的远场区域,用Euler方程来描述其流动。在该方法中,将三阶逆风格式（MUSCL）与通量差分分裂方法相结合,无需添加人工黏性,因而可有效地减少旋翼尾迹数值耗散。?阌诹鞒》智蠼庖约爸芷谛员呓缣跫氖凳?采用了嵌套网格方法,并给出旋翼网格与背景网格交界面的信息传递方式。应用所建立的计算方法,首先对二维翼型、三维M6机翼的流场进行了数值模拟,以验证计算方法;然后,着重计算了有实验结果可供对比的具有先进气动外形的HELISHAPE 7A模型旋翼和UH-60A直升机旋翼,通过计算得到旋翼表面压力分布、桨叶展向拉力系数分布、桨叶表面细节流动以及气动性能等,进一步验证了该方法的有效性。     

用隐式WENO格式计算悬停旋翼跨声速流场

发展了一种基于高阶迎风格式的计算悬停旋翼跨声速流场的隐式有限体积法。对流项采用Roe通量差分分裂格式,使用五阶WENO格式进行左右状态重构,并与MUSCL插值进行比较;粘性项采用中心有限体积法。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。数值模拟采用了一种能够有效传递网格间流场信息的重叠网格,其中使用了三层内边界和贡献边界的方法以便插值的直接进行。用该方法对一跨声速悬停旋翼粘性流场进行了数值计算,数值结果表明:与MUSCL格式相比,WENO格式对激波位置捕捉得更准确,具有更强的涡捕捉能力,同时还表明了WENO格式在很大程度上能够克服涡耗散问题。     

基于积分法的热流后处理方法

傅立叶定律统计热流时的扰动放大机制和数值计算格式的数值耗散是壁面热流网格依赖性强的主要原因,扰动放大机制同时也是热流收敛速度慢的主要原因。对边界层特征明显的流动的N-S方程数值求解,提出了积分法热流计算方法,热流由傅立叶定律统计的热流项和粘性功两项组成。平板层流、钝锥绕流气动加热算例表明计算结果的网格依赖性明显降低,收敛速度加快,为提高热流计算的可靠性和计算效率提供了一条新的思路。