基于WENO-分段线性格式的旋翼流场数值模拟

为了降低二阶JST(Jameson-schmidt-turkel)格式导致的数值耗散,发展了一套适合于格心格式的基于加权本质无振荡(Weighted essentially non-oscillatory,WENO)-分段线性格式的旋翼流场数值模拟方法。采用运动嵌套网格方法生成围绕旋翼的网格系统,主控方程选择Navier-Stokes(N-S)方程。为了更加有效地对桨尖涡等流动细节进行捕捉,在笛卡尔背景网格上采用七阶Roe-WENO格式计算对流通量;在桨叶贴体非结构网格上采用二阶精度的Roe-分段线性格式计算对流通量。时间离散采用了高效的双时间隐式LU-SGS(Lower upper symmetric Gauss-Seidel)方法进行时间推进。最后,应用上述方法对悬停状态的C-T(Caradonna-tung)旋翼和Helishape 7A旋翼进行了数值模拟,将数值计算结果与实验数据进行了对比,计算值与实验值吻合较好;并将桨尖涡模拟效果与二阶JST格式的模拟效果进行了对比。对比结果表明:本文方法能有效对旋翼流场进行计算,且在相同计算条件下,WENO-分段线性格式能够更有效地捕捉旋翼涡流场的流动特性,表明在计算旋翼涡流场时WENO-分段线性格式相比传统二阶JST格式具有更低的数值耗散。     

三角网格面自动铺丝定角度路径规划算法

研究三角网格曲面上的路径规划问题,以局部测地线为补充,提出了改进的定角度初始路径规划算法。该算法避免了传统参数化求解过程的无解和多解情况,尤其适用于曲率变化较大的复杂曲面和拼接曲面。采用等距覆盖算法来获得均匀的满铺路径,并分别提出基于路径样条线末端延伸和芯模曲面整体延拓的两种边界处理思路,以解决偏移路径无法到达曲面边界的问题。以某机翼模型为算例对上述算法进行验证,利用计算机辅助三维交互应用(Computer aided three-dimensional interactive application,CATIA)的逆向工程模块实现了完整的算法可视化,建立了三维建模软件与路径规划算法的系统关系。验证表明本文提出的算法可靠、系统、适应性强,具有工程应用价值。     

鸭式布局战斗机非常规机动的流场机理数值分析

以先进战斗机的非常规机动为对象,发展了一种适用于大幅度运动变化的非结构嵌套网格生成方法,建立了一整套非定常流场N-S方程数值求解方法。在对三角翼动态气动特性计算验证的基础上,模拟了飞机过失速机动条件下飞行姿态和来流速度的变化特征,对鸭式布局战斗机"眼镜蛇机动"的非定常涡结构、非定常气动力效应和气动特性进行了数值研究,揭示了鸭式布局战斗机"眼镜蛇机动"的非定常流场机理。     

多区域多重网格法在直升机粒子分离器中的应用

首次将多区域多重网格法用在一般曲线坐标系下交错网格布局的 SIMPLEC算法中,在直升机粒子分离器这类独特结构中的成功应用,突破多重网格方法在多连通域中运用的局限性,使复杂问题变得简单。结果表明,它不仅能加快流场求解的收敛速度,而且还能提高计算精度。     

两种多重网格法求解径流叶轮机械可压流动

为了加快数值求解叶轮机械内部流场的收敛速度,并改善计算结果,在数值求解径流叶轮机械可压流动的计算中,采用了由微元、二重和三重网格逐级转入和同时转入的两种多重网格法。对德国宇航局Ｋｒａｉｎ Ｈ设计的离心压气机内部流场的数值计算表明,两种多重网格法的计算结果差别不大,尽管在收敛速度方面存在差异,但均与测量数据吻合比较好。     

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

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

裂纹扩展的无网格数值模拟方法

疲劳断裂是航空材料的重要失效形式 ,由于裂纹尖端应力存在奇异性 ,传统有限元方法模拟裂纹沿任意路径扩展存在很多不足。作为一种新兴的数值模拟方法 ,无网格计算只需将求解问题离散为独立的节点 ,计算过程中可以实时跟踪裂纹尖端区域进行局布细化。将连续的裂纹扩展过程看作多个线性增量 ,每一个增量内裂纹扩展角根据应力强度因子确定 ,通过在裂纹尖端细化节点和引入外部基函数提高了计算精度。本文给出了应用无网格方法模拟裂纹扩展过程的关键技术和计算流程 ,通过对带有中心斜裂纹的 Ti-6 Al-4 V合金平板进行分析 ,预测得到的裂纹扩展路径与实验值吻合的较好。     

二维非结构网格的一个TVD型有限体积方法

考虑双曲型守恒律方程,对二维非结构三角形网格给出一种 TVD型有限体积方法,主要思想是在一阶单调格式的基础上,在每一个单元上对变量作单调线性重构函数,时间离散采用二阶 Runge-Kutta方法。通过计算分析了该方法的精度,对平面激波反射和空穴流动的计算结果表明该方法是成功的。     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

用隐式多重网格法计算三维粘性流动

求解高雷诺湍流流动时,边界层法向网格间距较流向和展向相比非常小,因此边界层中存在高度伸缩的网络,这将大大降低多重网格的求解效率。通过谐调的处理多重网格过程的各个细节,既提高了解的精度,又克服了网格展弦比的影响,使得计算效率提高了6～7倍。