结构/非结构混合网格数值模拟栅格翼

为了详细了解栅格翼内部的流动特性,设计满足飞行器性能要求的栅格翼,本文开展了多种不同单栅格外形的气动特性数值计算,对单个栅格的外形进行了初步选型研究;然后采用结构网格、结构与非结构混合网格求解NS方程,对带有两片栅格翼的导弹进行了数值模拟,结构网格计算结果与结构非结构混合网格计算结果一致,并与国外文献的计算结果和风洞试验结果进行了对比,三者的结果吻合较好.采用混合网格比结构网格的网格规模小许多,可以节省大量计算机内存和计算时间,是一种比较好的栅格翼数值模拟方法.     

某栅格翼导弹减阻外形设计及其阻力特性数值研究

介绍和分析了国内外栅格翼的研究现状和技术难点。为了详细了解栅格翼气动布局栅格内部的流动特性,以及其与阻力设计满足飞行器性能要求的栅格翼,本文开展了几种不同单栅格形状外形的气动特性数值计算。针对传统数值计算中遇到的栅格翼计算网格生成困难、工作量大的问题,利用非结构空间直角网格对于复杂外形易于生成的特点对单个栅格的外形进行了初步选型研究。对于常规导弹和栅格翼导弹的气动性能进行了对比。     

基于结构网格的栅格翼绕流数值模拟

针对栅格翼这种复杂的构造形式,采取有效的多块结构网格生成策略,对栅格翼导弹生成高质量贴体的结构网格。在此基础上,通过求解N-S方程,对栅格翼导弹亚、跨、超声速绕流流场进行了数值模拟。计算结果与风洞试验结果吻合较好,能够准确地捕捉到栅格框内复杂流场的流动特性,从而结构网格生成策略与计算方法得到验证,为栅格翼导弹绕流流场的精确数值模拟提供了技术途径。     

带控制律导弹投放数值模拟

基于非结构重叠网格技术,耦合舵面控制律模块,采用非定常雷诺平均NS方程,建立了模拟控制律作用下导弹投放分离动态特性的CFD数值方法.非结构重叠网格生成技术以“物面相交”方法确定洞边界,在多套网格重叠部分采用算术平均进行物理量的计算传递.采用机翼/挂架/带舵外挂物模型,以及不同的舵面控制规律展开了对比计算,分析了结果.并比较了某导弹有无控制律作用时从某战斗机分离的动态特性,研究了增加控制律的必要性.     

带控制舵弹体气动特性分析

本文用数值模拟方法研究了带控制舵弹体的气动特性。该问题外形复杂，在控制舵和弹身之间会出现明显的缝隙。为了解决此类外形的气动力计算问题，我们发展了结构网格和非结构网格的混合网格技术。这样一方面充分发挥了贴体结构网格质量好、效率高的优点，另一方面充分发挥了非结构网格几何灵活性好、易于模拟复杂外形的优点。数值计算与风洞实验相当吻合，说明本方法具有良好的精度。最后对两种外形进行大量计算并进行气动特性分析，给出了静稳定特性和配平特性的研究结果。     

前缘后掠式栅格翼升阻特性研究

栅格翼阻力大的缺点一直是需要解决的问题,前期研究显示,栅格翼翼面后掠能有效减小超声速阶段的波阻。本文采用数值计算的方法,对栅格翼前缘不同后掠方式、后掠角度及前缘后掠贴体型栅格翼进行了分析研究,数值结果显示栅格翼前缘后掠能有效减小波阻,其中栅格格间交接点为尖点的后掠方式减阻效果更好;栅格翼前缘后掠角度对其升力特性的影响较大,主要体现于对栅格内部气流壅塞现象的改善方面;和背风面为平面的栅格翼相比较,前缘后掠贴体型栅格翼的升阻特性表现更优。     

栅格翼导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术,对栅格翼导弹黏性流场进行数值模拟,预测其气动特性.计算的马赫数为0.7～2.5.用有限体积法和LU-SGS算法求解N-S方程.计算与实验数据以及栅格翼附近流场结构进行了比较.     

混合网格方法研究栅格翼亚跨超声速气动特性

采用基于结构/非结构混合网格的CFD方法对栅格翼的亚跨声速气动特性进行了研究,计算方法经过试验数据的验证,可达到工程精度。对单独栅格翼的研究表明,跨声速壅塞时,流动通过栅格前方亚声速气流减速从栅格外侧溢流而实现流量调节;超声速壅塞时,流动通过在栅格翼上产生网状脱体激波使来流减速至亚声速来进行流量调节,且计算结果与理论估算上下临界马赫数范围一致。对某栅格翼气动布局导弹的计算结果表明,在跨声速壅塞区间内,壅塞有自我调节作用,通过溢流使栅格通道内维持相似流动,从而使栅格气动性能保持基本平稳;同时由于溢流使流量下降,栅格翼升力效率下降,导弹静稳定性减弱。     

翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。     

栅格冀导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术，对栅格翼导弹黏性流场进行数值模拟，预测其气动特性。计算的马赫数为0．7～2．5。用有限体积法和LU-SGS算法求解N-S方程。计算与实验数据以及栅格翼附近流场结构进行了比较。     

复杂积冰翼型气动性能分析

为了克服复杂外形难以划分结构化网格的问题, 采用结构/非结构混合多块网格对复杂积冰翼型的气动性能-进行了分析计算.在靠近积冰边界的内层采用了非结构网格, 在非结构网格外采用结构化网格以节省存储空间和计算时间.不同类型网格之间采用了点点对接的方式以简化通量重构的处理方式.整个流动区域采用SST湍流模型, 求解了雷诺平均N-S方程.数值计算结果表明结构/非结构混合多块网格在复杂积冰翼型气动性能计算上的可行性, 积冰对翼型的气动性能造成了严重的影响.      

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

X—38飞行器无粘绕流的非结构网格模拟

本文采用非结构网格，利用DLR所发展的计算软件TAU对X－38飞行器高超声速无粘绕流进行数值计算。为了生成有效的计算网格及提高对流场细节的分辨率，在求解过程中应用了网格自适应,并研究了网格密度对气动力的影响程度；同时与现有的实验结果和用结构网格得到的结果进行了比较，符合程度令人满意。     

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难，将分块网格策略与嵌套网格方法融合，提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中，将整体网格按照倾转旋翼机结构特性进行了分块生成，在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格，首先，引入Inverse Map方法来辅助部件网格定位，采用新颖的多向投影方法保证了部件间洞边界的连续性；其次，基于独立挖洞嵌套策略完成机身与背景网格的组合，并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场，背景网格区域采用欧拉方程，并运用了SPMD(single program, multiple data)模式的并行加速技术，建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例，验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究，分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性，计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性，并能够有效地用于倾转旋翼机的气动特性分析，加速比能超过5.0。      

基于非结构网格的气动弹性数值方法研究

基于非结构网格和网格变形技术,研究了耦合求解流体和弹性结构体互相作用问题的数值方法.气动弹性模型采用了耦合求解可压缩雷诺平均的Navier-Stokes方程和线性结构动力学方程的方法.为了能处理在复杂几何域内的问题,采用了基于混合单元的非结构网格有限体积方法.开发的方法在二维翼型简谐振动进行了验证并应用于AGARD 445.6翼型颤振边界的数值模拟.      

曲面形栅格翼气动特性研究

与常规栅格翼布局相比，曲面形栅格翼布局减小了栅格翼翼元的等固壁通道的长度，从而有效降低了翼元内气流的壅塞，另外，这种栅格翼布局结构简单，易于折叠，减小了飞行器的轮廓尺寸，具有很好的工程应用前景。本文通过数值模拟方法，在亚、跨、超声速条件下，研究了曲面形栅格翼布局的减阻效果和减阻机理，对比了曲面形栅格翼和常规后掠形栅格翼布局的减阻效果，并对不同迎风方式的曲面形栅格翼气动特性进行了分析。     

二维结构/非结构混合网格的生成及其流场模拟

采用分区求解方法在二维结构/非结构混合网格中对粘性绕流作了数值模拟。利用层面推进法生成绕物体的O形结构网格,利用推进阵面法在结构网格外生成非结构网格。在结构网格区内,采用B-L湍流模型,求解雷诺平均N-S方程。在非结构网格区内求解Euler方程。两区的解通过在交界面上满足通量守恒条件进行耦合。数值算例表明了方法的有效性和实用性。