飞机/外挂物组合体绕流的数值计算

采用嵌合体埋入方法计算带外挂翼-身组合体复杂外形的流动,将飞机/外挂物外形的流动区域分为机翼-机身和外挂物两个子区域,在机翼-机身子域上求解欧拉方程,在外挂物子域上求解全速势方程,通过两个流场的耦合迭代计算全流场的解。对C-5A飞机及某战斗机的计算表明,数值模拟与实验吻合。实现了不同子域流场使用不同方程的计算方法。     

外挂物投放技术的数值模拟

应用Icemcfd和Cart3d软件对飞机投放外挂物的流场进行了数值模拟：对欧拉方程用有限体积法进行离散，运用自动及人工干预网格细化、多重网格技术的加速收敛、自动位置和网格更新的六自由度投掷分析，得到外挂物的分离轨迹。计算结果与AEDC的试验结果作对比，两者吻合较好。     

翼下外挂物对机翼振动及颤振特性的影响研究

应用MSC／NASTRAN软件计算分析了某型飞机翼下不同位置挂装外挂物后，对机翼振动特性和机翼／外挂物系统颤振特性的影响，计算结果表明：翼下外挂物对机翼二弯和机翼一扭模态的影响较大。     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

高速风洞品字形低阻弹外挂测力试验研究

在ＦＬ－２４风洞中，采用８台天平分别对母机，六枚“品”字形低阻弹及其外挂组合本同时进行了外挂测力试验，结果表明，各外挂物受机翼下洗流影响，其法向力，俯仰力矩特性与其单个外挂物在自由流中呈现的特性大不相同，翼下挂低阻弹及其挂架受侧洗流影响较大，零侧滑角时产季了横向气动特性，处于机翼复杂流场下的外挂物，除滚转力矩量值不大外侧向力和偏航力矩量值分别与法向力和俯仰力矩量值相当。     

外挂物在飞机上悬挂的气动相关性研究

悬挂外挂物的气动相关性问题包括确定外挂物对飞机气动特性的影响;确定处于主机干扰流场中外挂物的气动特性;确定带外挂物飞机的颤振特性及确定外挂物安全投放条件等。本文试图通过歼×飞机2000多次风洞试验及其它试验和必要的计算分析对外挂物悬挂的气动相关性问题进行综合性研究和评述并给出一些结果,对一些研究方法进行了一般性讨论     

飞机外挂物离机过渡过程的有限元算法.

研究飞机外挂物在按给定轨迹下沉旋转的离机过渡过程中，流场非结构网格的调整和再生成方法以及过渡过程欧拉方程有限元数值解法。中指出外挂物运动时流场网格调整或再生成仅需在围绕在外挂物的一小区内进行以减少计算量，而且固定网格的欧拉方程的有限元数值解法也可适用于运动网格情况，只要将固定网格的守恒通量项Fk用运动网格的Ek来代替，并建立了Fk和Ek之间的关系。末给出了带外挂物飞机的实例计算结果。     

飞机外挂物投放安全性的参数影响分析

研究了飞机在不同飞行高度、速度和姿态下,不同外挂物质量对投放安全性的影响。采用CFD方法计算外挂物所受气动力,通过六自由度刚体动力学方程计算外挂物的投放轨迹和姿态变化。准定常方法在计算外挂物下落轨迹的每个步长内,认为外挂物和载机都处在定常流场中,载机在投放过程中保持原有运动状态。通过算例给出了典型外挂物在不同速度、高度、姿态下投放的运动轨迹,分析了对安全投放有重要影响的因素,提出了减少外挂物投放风险的建议。     

外挂物投放问题的二维非结构动态网格技术研究

以动态网格守恒型方程为控制方程,在二维非结构网格基础上发展了弹簧式光顺法和局部网格重构法两种动态网格变形技术。结合运动轨迹计算方程,将动态网格法应用于飞机外挂物投放问题中,成功模拟了导弹自由投放的整个过程,获得了详细的包括弹体运动速度、运动轨迹和受力情况等在内的丰富的流场信息。成功的应用不仅验证了动态网格法的可行性;而且也表明,动态网格法的进一步完善与发展将成为解决外挂物投放问题的有力工具。     

计算机模拟铰接式直升机的外挂施工

为了给直升机施工的安全分析提供理论基础和求解非线性外挂直升机运动方程,本文提出了用计算机模拟直升机外挂飞行的方法。所用直升机运动方程模拟程序采用美国宇航局兰利研究中心发布的铰接式单旋翼直升机模拟程序。期望航迹及外挂载荷的模拟,根据航迹及载荷反求操纵输入及配平参数的方法和运动方程的积分等将由本文给出。     

Collar网格的生成及其在嵌合体技术中的应用

本文应用双曲型方程生成表面网格和空间网格的方法成生Ｃｏｌｌａｒ网格,并将Ｃｏｌｌａｒ网格应用到嵌合体技术中,在保证计算网格的生成方便快捷而且网格质量高的前提下,有效地解决了为物体结合部的内外边界点寻找插值单元的问题。经数值计算验证,用本文方法求解有结合部的复杂组合体绕流可以得到精确可靠的数值结果。     

A new grid deformation technology with high quality and robustness based on quaternion

Quality and robustness of grid deformation is of the most importance in the field of aircraft design, and grid in high quality is essential for improving the precision of numerical simulation. In order to maintain the orthogonality of deformed grid, the displacement of grid points is divided into rotational and translational parts in this paper, and inverse distance weighted interpolation is used to transfer the changing location from boundary grid to the spatial grid. Moreover, the deformation of rotational part is implemented in combination with the exponential space mapping that improves the certainty and stability of quaternion interpolation. Furthermore, the new grid deformation technique named ‘‘layering blend deformation' is built based on the basic quaternion technique, which combines the layering arithmetic with transfinite interpolation(TFI) technique. Then the proposed technique is applied in the movement of airfoil, parametric modeling, and the deformation of complex configuration, in which the robustness of grid quality is tested. The results show that the new method has the capacity to deal with the problems with large deformation, and the ‘‘layering blend deformation' improves the efficiency and quality of the basic quaternion deformation method significantly.     

燃气轮机燃烧室中贴体网格的生成与应用

运用贴体坐标系统分别对燃气轮机燃烧室的冷态和热态流场进行了数值模拟。燃烧室的曲线壁面边界用Ｔｈｏｍｐｓｏｎ的非正交贴体坐标网格处理，并把ＳＩＭＰＬＥ算法应用到曲线坐标系下求解各守恒方程。采用双方程ｋ－ε模型来描述紊流特性；采用按简单的单步化学反应假设的Ａｒｒｈｅｎｉｕｓ－ＥＢＵ模型表述紊流燃烧；采用热通量法辐射模型来估算辐射通量和壁面温度。另外，还研究了冷态时不同旋流数和紊流强度对气流结构的影响。     

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

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

一种包含运动边界的高精度流场数值计算方法

为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。     

二维轴流压气机转—静非定常Euler方程数值计算

将Euler方程写在包括转子和静子的绝对坐标系下,通过在转—静之间设计一个可以随时间扭动的平行四边形网格计算了二维轴流压气机转—静非定常流场。由于整个方程统一写在动网格下,转—静之间没有使用插值。      

基于自由型曲面的贴体网格生成

本文建立了基于自由型曲面的曲面网格生成的数学模型,生成了曲面网格,把它作为三维Ｄｉｒｉｃｈｌｅｔ问题的边界条件,求解三维Ｐｏｉｓｓｏｎ方程,获得了满足数值模拟要求的三维贴体网格。采用上述方法,完成了基于自由型曲面的某直升机调座舱三维贴体网格的生成,结果令人满意。     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

Numerical simulation of aerodynamic interaction for a tilt rotor aircraft in helicopter mode

A rotor CFD solver is developed for simulating the aerodynamic interaction phenomenon among rotor, wing and fuselage of a tilt rotor aircraft in its helicopter mode. The unsteady Navier–Stokes equations are discretized in inertial frame and embedded grid system is adopted for describing the relative motion among blades and nacelle/wing/fuselage. A combination of multi-layer embedded grid and ‘‘extended hole fringe" technique is complemented in original grid system to tackle grid assembly difficulties arising from the narrow space among different aerodynamic components, and to improve the interpolation precision by decreasing the cell volume discrepancy among different grid blocks. An overall donor cell searching and automatic hole cutting technique is used for grid assembly, and the solution processes are speeded up by introduction of Open MP parallel method. Based on this solver, flow fields and aerodynamics of a tilt rotor aircraft in hover are simulated with several rotor collective angles, and the corresponding states of an isolated rotor and rotor/wing/fuselage model are also computed to obtain reference solution.Aerodynamic interference influences among the rotor and wing/fuselage/nacelle are analyzed,and some meaningful conclusions are drawn.