关于非结构网格生成中的背景函数

背景信息对于用阵面推进法生成非结构网格是一个关键的环节.文章基于众多文献中提到的背景信息提供方法,提出一种更为实用的背景信息提供技术.此方法不仅避免了不必要的机时,内存的消耗,又使网格具有光滑性,网格的生成具有灵活性,方便性的特点.基于生成的网格采用有限体积法研究无粘高超音速化学非平衡流的流场分布,生成的网格能够成功的运用于高超音速化学非平衡流的数值计算。     

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

复杂外形网格生成技术及应用

以某型双立尾战斗机为例,给出了实际复杂型号飞行器结构网格生成的一种分块方案及实现过程,用一种矢性三次多项式插值的方法生成相邻块公共交界面的网格,用Hilgenstock方法修正源项,求解椭圆型方程生成块内空间网格,使网格线与边界面夹角、与边界面间距控制为期望的值.在所生成的双立尾战斗机分块结构网格中求解Euler方程,获得了合理的全机气动力系数、翼面压力分布和绕流图谱结果.     

2D约束Delaunay剖分生成表面模型的表面网格

提出一种生成表面模型的表面非结构化网格的方法.该方法仅假定表面模型的每个面片在取定适当的投影平面时可看成是一个单值函数.对表面模型的每个面片,首先在其相应的投影面上进行二维约束Delaunay剖分,然后对网格顶点进行插值.由于采用联动剖分的策略使得单独生成的各个表面网格在共同边界处匹配.与通常的生成表面网格的网格前沿法相比,本方法无需预先对域的边界进行离散,边界的离散体现在约束Delaunay剖分中恢复限定线段的边界细分过程中,减少了用户干预.通过合理地指定2D约束Delaunay三角化时网格单元尺寸分布函数,可以有效控制表面网格对表面模型的逼近精度和实现自适应的表面网格.实验结果证明了本算法的有效性.     

非结构重叠网格显式装配算法

针对重叠网格中洞映射法占用过多物理内存的问题,发展了一种改进型洞映射法;基于相邻单元搜索法,发展了一种基于相邻阵面的贡献单元搜索法;通过将割补法与隐式切割技术相结合提出了一种非结构重叠网格显式装配算法。该算法首先生成一套包围物面的笛卡儿网格,其次存储所有与物面边界相交的笛卡儿网格信息,最后根据所存储笛卡儿网格与所需判断的网格单元的相对位置来判断其是否为洞内单元。在成功判断出所有洞内单元后,以当前洞边界为初始阵面推进,同时以各个网格单元的物面距离为判别标准对重叠区域进行优化,生成最终插值边界。所提算法优化了传统非结构重叠网格装配过程,具有物理内存占用低,贡献单元搜索次数少以及计算效率高等特点。通过2个典型复杂流动算例验证了所提算法的准确性与适用性。      

3点逐步r型自适应网格算法

为了根据某个流场参数梯度较合理地分布网格同时又不增加网格点的数量,提出一种基以总变差为基础的r型3点逐步自适应网格算法,通过3点逐步调整网格点的位置以减小网格线的扭曲,生成质量较好的计算网格.以绕二维圆柱和三维双椭球的超声速流动为算例,用Euler方程模拟,以流场的压力总变差来进行自适应,通过分别用所提出的自适应方法生成的自适应网格与用原始网格计算所得的流场结果对比,计算结果表明,用自适应网格计算所得的激波比用初始网格计算所得的激波薄,采用自适应方法所计算出的流场数值解具有更高的分辨率.      

结构/非结构网格中的数值优化计算

非结构网格因其在处理外形方面所具有的特殊优越性而倍受关注.在结构/非结构网格中做了粘性流动的数值优化计算,优化计算包括单目标和双目标气动优化,以及反设计计算,即从初始翼型不断修型而实现给定翼型的外形及压力分布的目标.算例及其结果分析证明了该方法用于翼型的修型和设计计算的有效性.      

燃烧室三维贴体网格的生成方法研究

研究了燃烧室三维贴体网格人机交互生成的方法.使用该方法可以生成包括扩压器、内外环通道和火焰筒在内的燃烧室流场计算所需的整个网格.可以对网格的疏密和正交性进行调整.可以处理火焰筒上的主燃孔,掺混孔和气膜孔.并能够输出边界曲面上网格点的类型,用于流场计算.网格边界线的输入和边界线上的网格点分布采用人机交互方式生成.使用VISUAL BASIC语言编写程序界面.把整个网格生成及图形显示功能集成在一个软件中.     

B样条应用于超限插值的网格生成方法

采用B样条曲线曲面生成技术与传统网格生成方法中超限插值技术,并将两者结合起,形成一种计算网格的生成方法.它发挥了两者的优点和长处:能充分反映复杂外形的完全信息,提高对外形的适应程度;可控制切矢的边界条件,提高生成网格的质量.算例表明了上述特点.      

合成射流控制NACA0015翼型大攻角流动分离

为了研究合成射流激励器处于NACA0015翼型回流区时对其分离流动的控制,采用商用计算流体力学软件Fluent 6.1求解Reynolds平均Navier-Stokes方程,通过对翼型气动力特性、脱落漩涡结构以及射流孔口附近流动结构的分析,揭示了合成射流处于分离区时对边界层控制的机理.结果表明:当合成射流孔口处于回流区时仍可有效推迟翼面边界层分离点,缩小回流区范围,从而有效提高翼型的升力.当射流方向垂直于壁面,无量纲频率以及吹气速度比都等于1时,翼型平均升力系数提高40%左右.      

超临界层流翼型优化设计策略

针对超临界层流翼型设计问题,提出一种两轮优化策略。采用γ-Reθt转捩模型耦合剪切应力输运(SST)模式的湍流模型对翼型边界层转捩进行预测。翼型几何参数化建模采用形状分类函数转换(CST)方法,设计变量为描述翼型几何特征的参数。第1轮优化的目的是尽量提高层流区域的比例,气动分析模型为基于Kriging模型的代理模型,优化算法为遗传算法,通过优化获得满足约束要求的层流翼型。第2轮优化目的是对第1轮优化获得的翼型进行微调,进一步提高翼型的升阻比,气动分析直接采用CFD程序,优化算法采用基于梯度的优化算法。算例表明,应用本文提出的两轮优化策略,可将超临界翼型NASA SC(2)0412优化设计成超临界层流翼型,翼型的上下表面层流区比例分别达到了55.5%和47.0%,升阻比提高了38.1%。     

NFFD控制点分布对气动外形优化的影响

基于非均匀有理B样条的自由型面变形(NFFD)技术具有对变形对象普适性和控制点影响区域局部性的特点,广泛应用于气动外形优化。本文通过扩展控制体和合理布置外侧控制点,实现了NFFD技术在参数化曲面并改变曲面形状时,同步变形控制体内的表面网格和空间网格,并从理论上保证了控制体边界内外的网格协调。基于离散伴随方法求取梯度,分别采用拟牛顿(QN)法和序列二次规划(SQP)优化方法,通过从初始翼型NACA0012到标准翼型EH1590的反设计,研究了NFFD控制点数量和分布对设计结果的影响。在某飞翼标模单点全机升阻比优化应用中,改进控制点分布后获得了更高的升阻比,收敛速度显著提高。     

两种手形对游泳推进效率影响的数值模拟

针对采用何种手形游泳能产生较高的推进效率,运用计算流体力学(CFD,Computational Fluid Dynamics)方法,建立了国内某女运动员手掌和前臂在五指并拢和五指分开两种手形的三维非结构网格模型;计算了90°攻角划水时两种不同手形的手掌和前臂在不同来流速度下的推进阻力、推进升力及其系数值,依据计算结果对手周围的流场进行了分析,研究了不同手形对游泳推进效率的影响.结果表明压差阻力是推进阻力的主要组成部分,在90°攻角划水时五指并拢具有更高推进效率.      

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.