计算气动声学高阶差分格式的GPU并行实现

以圆管构型的声传播为分析对象,研究了基于图形处理器GPU的计算气动声学(Computational Aeroacous-tics,CAA)高阶有限差分算法的并行实现,并与CPU串行及MPI并行实现作了对比分析。首先介绍了管道简化模型的2.5维线化欧拉方程和GPU的编程模式以及调优参考准则,然后给出了相关物理量的空间离散方法的GPU实现。数值实验的结果表明,与CPU串行及MPI并行程序的结果相比,使用GPU的程序实现在达到与MPI并行同样的计算效率时,可以使用更少的计算资源。较之cluster上串行算法,工作站上GPU并行算法在使用不同网格规模的情况下可达到的3倍多的加速比。     

某型飞机起落架结构件气动噪声仿真与试验研究

起落架噪声是飞机着陆阶段噪声的主要组成部分。以某型飞机前起落架为研究对象,通过分离涡模拟方法对其支柱及扭力臂结构件简化模型的周围流场进行非定常计算,利用Fw—H方程积分法对各部件表面产生的声场进行求解,分析缓冲支柱及扭力臂结构件气动噪声的产生机制、声源特性。对该飞机起落架支柱及扭力臂结构件进行声学风洞试验,通过麦克风对噪声的测量获得结构件噪声频谱特性。仿真及试验结果均表明：支柱及扭力臂结构件气动噪声包含支柱和扭力臂引起的钝体扰流噪声和两者相对位置引起的干扰噪声,支柱噪声对总噪声的贡献大于扭力臂噪声,噪声辐射特性具有偶极子声源的辐射特性。     

改进S-A湍流模型对横向射流的CFD模拟

通过在Spalart-Allmaras(S-A)湍流模型中加入应变率张量,来改善模型对横向射流(JICF)的计算精度。模拟了速度比分别为0.5和1.5工况下的横向射流。经实验值验证,结果表明:改进后的S-A模型预测精度有所提升;并对湍流模型改进部分对横向射流的作用机制做出了分析。将改进后S-A模型的计算结果与k-ε模型,雷诺应力模型(RSM)和大涡模拟(LES)的结果进行比较,发现S-A模型在速度比为0.5时精度最高,而在速度比为1.5工况下没有明显优势;并对这一现象做出了合理分析。      

研发铝合金的集成计算材料工程

用于铝合金的集成计算材料工程是将微观(10-10~10-8m)、细观(10-8~10-4m)、介观(10-4~10-2m)和宏观(10-2~10 m)等多尺度计算模拟和关键实验集成到铝合金设计开发的全过程中,通过成分-工艺-结构-性能的集成化,把铝合金的研发由传统经验式提升到以组织演化及其与性能相关性为基础的科学设计上,从而大大加快其研发速度,降低研发成本。本文详细阐述了原子尺度模拟、相图计算、相场、元胞自动机和有限元等计算模拟方法及微结构表征和性能测定的实验方法,论述了其在铝合金研发中所发挥的具体作用。基于集成计算材料工程,提出了从用户需要、设计制备和工业生产3个层面研发铝合金的具体框架。通过2个应用实例,展示了集成计算材料工程在铝合金研发中的强大功能,这也为新型铝合金及其它新材料的设计和开发提供了新模式。     

多块多重网格法及其跨声速转子内流并行模拟

为提高压气机内部复杂流动计算的收敛速度,发展了一种多块多重网格上的可压缩流动并行计算流体动力学(CFD)方法,数值模拟了跨声速转子NASA Rotor 35内部流动及气动性能,重点分析了对接分区方式和内界面通信模式对并行性能的影响。结果表明:对接分区方法和多块多重网格内界面数据处理方法能够保证串、并行CFD方法结果基本一致,且与实验数据吻合良好;并行精度基本不受分区数目和通信模式的影响,并行效率随着相对通信负载的增大而降低。所发展的并行计算方法对跨声速转子气动分析和设计具有一定的可靠性,为开发大型流体机械复杂流动的基础并行算法及软件提供了参考价值。     

一种高效的组合动导数分量数值模拟算法

基于传统网格速度法,发展了在流场网格不变的情况下,数值模拟沉浮运动、纯俯仰运动和俯仰振荡运动的算法,提出了一套直接求解旋转阻尼导数和洗流时差导数的方法,该方法不需要实时更新网格,减少了计算时间和所需内存,避免了负体积的出现。首先通过NACA 0006的阵风响应和NACA 0012翼型的强迫俯仰振荡,验证了程序和方法的正确性;然后数值模拟了BFM(basic finner missile)标准模型的沉浮运动、纯俯仰运动和俯仰振荡运动,并计算了其旋转阻尼导数、洗流时差导数和组合动导数,所得结果与动网格方法和风洞试验值吻合,最大误差不超过4%,进一步验证了该方法的正确性。     

边界层抽吸效应对分布式推进系统气动性能影响数值研究

为探究边界层抽吸（BLI）效应对飞机的气动性能影响,采用基于沿流线体积力模型的飞机/发动机一体化数值模拟方法对某分布式推进系统进行了计算和分析。结果显示,BLI效应主要影响飞机中心体部分及发动机外整流罩的气体流动,对翼身融合体（BWB）的融合段及外翼段的升阻力影响较小。保持飞行条件和飞行攻角不变,飞机的升、阻力系数均随着无量纲化的发动机流量增加而增大,存在最佳无量纲化的发动机流量对应最大升阻比。发动机的安装位置直接影响机身表面的局部超声区和整流罩外的激波分布,布局在机身尾缘处会获得更好的升阻比,最大升阻比对应的无量纲化的发动机流量为0.65。      

基于RANS/LES方法的超声速底部流场数值模拟

分别采用基于两方程k-ω剪切应力输运(SST)湍流模型的延迟DES(DDES)、更改的DDES(MDDES)和改进的DDES(IDDES)方法,并引入可压缩修正,结合三阶MUSCL-Roe和五阶WENO-Roe两种空间离散格式,针对超声速底部的复杂流动现象,开展了数值模拟研究。计算结果表明本文方法能够捕捉到超声速底部流动中丰富的湍流结构,通过分析计算结果对超声速底部的流动机理有了进一步的认识,为下一步的超声速底部流动减阻改进和雷诺平均NavierStokes/大涡模拟(RANS/LES)方法在非定常高可压缩性流动中的应用提供了参考。通过对比分析不同空间离散格式的计算结果研究了数值耗散对计算的影响,五阶WENO-Roe格式的计算结果与实验结果吻合良好;对不同RANS/LES混合方法的计算结果进行了对比分析,结果表明IDDES方法在近壁区的表现优于DDES和MDDES方法。     

DCS5在计算气动声学中的应用研究

对流项采用五阶线性耗散紧致格式(DCS5)离散,粘性项采用6阶中心紧致格式离散,时间项采用5级5阶龙格-库塔时间推进格式求解流体动力学方程组,并由此建立近声场的直接模拟方法。为了应用结构网格在复杂几何形状条件下对流场进行计算,发展了高精度的多块特征对接边界技术。完成了低速无粘振荡圆柱和粘性圆柱绕流等标准计算声学算例,获得了令人满意的结果。     

重叠网格技术中割补法的研究与改进

针对计算流体力学中重叠网格生成困难的问题,分析了重叠网格自动生成技术中割补法的主要特点和基本流程,指出其在洞映射单元属性的判断方法和洞边界外推过程中出现的"孤点"现象等方面的问题,提出了合理的改进方案,并将改进后的方法应用到自适应嵌入网格和DRAGON网格的生成中,最后编制出通用程序并给出了三维算例.实际应用表明,改进后的方法具有更稳定、自动化程度更高、适用范围更广的特点.