三维舵面绕流的N-S方程数值计算研究

通过自动生成复杂构型的非结构混合网格。采用格心格式的有限体积法以适应不同的网格单元，应用Spalart-Allmaras湍流模型进行绕三维舵面流动的N-S方程数值模拟，将部分计算铰链力矩结果与风洞试验结果进行了对比分析。说明了计算方法对粘性流动良好的适应性,分析了迎角、马赫数、雷诺数对舵面铰链力矩的影响以及计算方法的可靠性和计算结果的合理性。     

一种新的动网格方法及其应用

 介绍了一种新型的动网格方法。该方法以Delauney图为基础进行整个网格的移动变形：首先是生成满足Delauney特性的Delauney图;然后把对网格移动的特性映射到Delauney图的移动上;最后把网格点在Delauney图以映射的方式进行重新定位而得到新的网格点坐标,从而得到移动后的网格。该方法不需要迭代,因而效率高,并且能适用于任意拓扑结构的网格。以翼身组合体的机翼颤振、翼身融合体的变体过程的非定常流动等为例来说明该方法的高效性及适用面广等特点。     

三维动态非结构重叠网格Navier-Stokes方程并行算法

 提出了一种三维动态非结构重叠网格Navier-Stokes(N-S)方程的并行计算方法。N-S方程的空间离散采用格点有限体积方法,时间离散采用隐式的双时间步长方法。应用一方程Spalart-Allmaras(S-A)模型来计算湍流黏性。并行计算采用动态的区域分裂方法,在每一物理时间步利用METIS网格分区系统对网格进行分区。为了实现各CPU之间的负载均衡,每块网格都按CPU个数进行分区并对活动节点和非活动节点进行了加权处理。最后,通过对外挂物投放无黏流动的数值模拟和内埋武器弹舱开启黏性流动的数值模拟,验证了该并行程序的准确性、高性能并行计算以及处理复杂几何外形的能力。     

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

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

非结构嵌套网格的直升机旋翼/机身前飞流场数值模拟

建立了一个适用于前飞状态直升机单独旋翼、机身以及旋翼/机身组合体流场的数值计算方法及模型.在该方法中,使用非结构运动嵌套网格描述桨叶之间以及桨叶与机身之间的空间位置变化关系,控制方程采用惯性坐标系下的非定常Navier-Stokes(N-S)主控方程,空间方向上使用二阶迎风格式,时间方向上使用Lower-upper symmetric Gauss-Seidel(LU-SGS)格式.对于求解中遇到的旋翼/机身间距离过近所导致的贡献单元生成困难问题提出了一种新的解决方法.应用所建立的方法,分别对单独旋翼、机身及旋翼机身组合体的前飞绕流场进行了数值模拟与分析,计算结果与试验结果的对比表明本文建立的方法对直升机的前飞绕流场的计算和分析是很有效的.      

机体/动力装置一体化分析中的动力影响效应数值模拟

采用合适的进、排气边界条件, 对带动力的单独涡扇发动机模型和带动力的DLR-F6翼身架舱复杂组合体(WBPN)模型的流场进行了数值模拟和分析.运用Roe三阶迎风偏置通量差分裂方法和隐式近似因子分解方法求解雷诺平均Navier-Stokes方程, 基于点搭接多块网格技术生成高质量计算网格, 通过多重网格技术来加速收敛.着重对动力发动机进、排气边界条件的给定进行了推导, 以准确地模拟发动机动力效应.计算结果与实验值吻合良好, 验证了所采用计算方法、网格分块策略和边界条件处理方法的正确性.      

基于非局部理论的一维损伤模型的有限元算法

为了解决损伤局部化问题的网格依赖性,采用非局部理论,从损伤本构关系中引入特征长度入手,建立非局部损伤本构模型。提出了一种新的单元,将位移、非局部应变和非局部应变梯度作为基本的未知量,引入了C-1连续的形函数,并推导了非局部损伤本构模型的有限元格式。在ABAQUS中进行了二次开发,嵌入了本文提出的单元和本构模型,并用ABAQUS软件分析了一维损伤结构的受力情况。计算结果表明,数值模拟与理论推导结果的总体趋势吻合得较好,证明了本文方法的正确性与可行性。     

DPWⅢ机翼和翼身组合体构型数值模拟

采用亚跨超CFD软件平台(TRIP)数值模拟了DPW Ⅲ提供的DPW_W1/W2两种机翼构型和DPW-F6/F6_FX2B两种翼身组合体构型,主要目的是通过两种机翼构型和两种翼身组合体构型的数值模拟,研究网格密度对运输机构型气动特性计算结果的影响。数值模拟采用的多块对接网格来自AIAA CFD Drag Prediction Workshop Ⅲ(DPW Ⅲ),采用National Transonic Facility(NTF)的试验结果和CFL3D的计算结果作对比。详细研究了网格密度对两种机翼构型和翼身组合体的总体气动特性和压力分布的影响。采用SST两方程模型计算两种构型均得到了网格收敛结果,网格密度主要影响压差阻力而对摩擦阻力影响较小,计算结果较好地预测了机翼和翼身组合体外形优化前后总体气动特性的变化量。     

三维非结构自适应多重网格技术

发展了一种基于几何外形的非结构自适应多重网格技术。根据流场参数的变化梯度确定加密边,由加密准则对粗网格进行自适应剖分得到分布合理的较密网格。通过预先生成的初始极密表面网格来将边界网格的加密点投影到边界上,使得边界保持初始外形。为了提高流场求解速度,将自适应加密后的一系列网格作为多重网格,采用了多重网格算法进行计算。M6机翼的流场计算表明该方法可以准确捕捉激波位置,并且表面的压力系数分布与实验值相吻合。     

局部流动分离诱导结构振动问题的数值模拟

基于非结构网格技术和SST-DES模型,采用流体与结构紧耦合策略,对飞机机身突出物（如雷达罩）下游的蒙皮振动过程进行了数值模拟计算。在求解NS方程过程中,使用有限体积法,结合有限元法求解结构运动方程。在时间域内模拟了不同结构参数和飞机不同迎角下的飞机蒙皮的结构响应,建立了亚音速条件下壁面流固耦合分析方法。计算结果表明：机身突出物诱导出的分离涡的确可以激发下游壁面的振动,其壁面振动的幅度与结构本身参数密切相关,振动主要为发生在壁面的法向上,随着迎角增加,壁面的振动减弱。通过基本原理和算例可以看出,方法可以预测不同条件下壁面的动态响应过程,适用范围宽,可以用于壁面结构的设计。