基于Delaunay图映射的动弹网格技术及其应用

在多块混合网格的基础上发展了一套基于Delaunay图映射方法的动弹网格技术,实现了针对气动弹性计算的网格变形运动方法,并使用基于多块混合网格的RANS方程的有限体积法耦合结构静平衡方程,针对大展弦比无人机在实际使用中可能出现的显著的静气动弹性问题进行了计算、分析和研究。研究结果表明,所发展的动弹网格方法为粘性混合网格的变形提供了较好的解决方案;RANS方程解算器具有较高的精度和较强的适用性,具备求解粘性流体与结构耦合等复杂航空气动弹性问题的能力。     

自适应笛卡尔网格Ghost Cell方法研究

在笛卡尔网格中,利用Ghost Cell方法处理浸入边界,模拟二维无粘可压缩流。针对静止物体,比较了各种不同的Ghost Cell边界条件下的熵误差,总压误差,以及阻力系数。此外为提高激波分辨率,将Ghost Cell方法与基于叉树结构的自适应笛卡尔网格算法相结合,数值结果显示本文的方法是切实可行的。     

基于CFD/CSD耦合方法的旋翼气动弹性载荷计算分析

为提高直升机前飞状态下旋翼非定常气动弹性载荷的预估精度,在旋翼气动弹性综合分析方法中引入旋翼CFD模块,建立了一套基于CFD/CSD松耦合分析的计算方法和程序。为高效解决流固耦合方法中由于桨叶挥舞、扭转等弹性变形带来的旋翼贴体网格变形问题,采用基于代数变换方法的网格变形技术,桨叶运动变形量和旋翼气动力信息通过流固交接面传递。旋翼流场分析方法中,主控方程采用耦合S-A湍流模型的Navier-Stokes方程,围绕旋翼流场的网格采用结构嵌套网格方法生成,无黏通量计算采用Roe格式,时间推进采用双时间法。旋翼结构分析中,考虑旋翼配平,基于Hamilton变分原理和20自由度Timoshenko梁模型求解弹性旋翼非线性运动方程。分别对CSD和CFD方法进行验证,在此基础上,计算了SA349/2旋翼桨叶在前飞状态下的非定常气动力、挥舞弯矩、摆振弯矩和扭转力矩,并与飞行测试数据进行了对比。计算表明:CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,精确捕捉桨叶表面压强峰值、激波位置等,表明本文发展的旋翼CFD/CSD耦合方法可以有效地运用到旋翼气动弹性载荷的预测分析中。     

基于CFD和混合配平算法的直升机旋翼地面效应模拟

建立了一套基于非结构网格技术和动量源模型的直升机旋翼计算流体力学(CFD)方法,用来模拟贴地飞行时直升机旋翼非定常气动特性。其中,控制方程采用三维Navier-Stokes方程,空间方向上采用Jameson格式,时间方向上采用五步Runge-Kutta迭代法,选用Spalart-Allmaras湍流模型。旋翼对流场的作用采用动量源项模拟,为更真实地模拟地面效应(IGE)的作用,采用了"移动地面"的物面边界来代替常规的"固定地面"边界,并对旋翼附近及旋翼与地面之间的网格进行加密处理,以提高地面涡的捕捉精度。考虑实际飞行环境下旋翼的运动和操纵,在流场计算时考虑旋翼配平特性。其中,配平方程的旋翼气动力通过CFD方法和动量-叶素组合理论模型的耦合计算给出,为了提高配平方法的鲁棒性和效率,提出并建立了基于遗传算法/拟牛顿法的高效混合迭代算法。运用所建立的方法,首先,选用有试验结果可供对比的算例计算了地面效应作用下的旋翼拉力增益、功率变化,验证了计算方法的有效性,解决了涡流理论方法较难模拟的"小速度前飞旋翼需用功率突增"问题。然后,着重研究了UH-60A直升机旋翼在不同离地高度、不同前进比状态,旋翼需用功率、诱导速度、地面涡及旋翼操纵的变化规律。计算结果表明:地面涡出现在较小的前进比范围内,随前进比的增大,地面涡在纵向平面将顺来流方向移动,在轴向方位靠近地面方向移动,直至最后不断减弱消失。     

基于非线性耦合本构关系的改进边界条件

非线性耦合本构关系（NCCR）模型是在Eu的广义流体动力学方程（GHE）基础上,通过绝热假设、Eu封闭和Myong简化推导出的关于非守恒量（黏性应力与热流）的非线性代数方程,有效拓展了线性的纳维-斯托克斯-傅里叶（NSF）本构模型在非平衡流动中的模拟能力,为快速准确模拟连续与稀薄耦合流动问题提供了强有力的理论工具。针对该模型开展滑移边界条件研究,结合努森层内物理量非线性分布的特点,提出一套在物面处与模型精度相一致的非线性修正滑移边界条件。在有限体积框架下,采用AUSMPW⁺格式和LU-SGS方法以及NCCR的完整耦合求解算法,对不同稀薄程度的高超声速单原子氩气圆柱绕流和平板绕流问题进行数值模拟。研究结果表明,基于NCCR模型的修正边界条件准确刻画出物面努森层内流动的非线性特点,有效提高了固壁滑移边界的精度。采用非线性修正边界的NCCR模型准确预测了连续流、滑移流和过渡流域的物面压力、摩阻与热流系数。     

超声速条件下内埋式武器分离特性的数值分析

采用非结构动网格技术,耦合求解欧拉方程和六自由度(6-DOF)运动方程,对内埋式导弹在亚声速和超声速条件下的分离轨迹和姿态变化进行了比较,分析了不同马赫数下气动力对导弹分离特性的影响,重点讨论了超声速条件下导弹分离过程中相对于载机的位置和导弹的姿态变化情况,分析了超声速飞行环境下导弹安全分离的条件.     

建立在视线方程基础上的空战动力学模型

首先定义并介绍了视线动力学方程,在此基础上讨论视线坐标系上的机动策略问题。介绍视线坐标系上的空战策略的一种决策方法,以及视线坐标系上的空战策略向机体坐标系上的转化问题。最后建立在机体系上的力和力矩方程,用来说明如何操控飞机达到所需要的空战策略。     

二维H2/O2混合和流动的直接数值模拟

用高阶精度的紧致差分方法,对低速的随时间发展的二维可压缩H-2/O-2混合层流动编制了直接数值模拟的计算机程序,对雷诺数Re为400、来流马赫数Ma分别为0.4(上部的H₂)和1.57(下部的O₂)的H₂/O₂亚跨超声速混合层流动进行了计算,得到了相应的计算结果,并对结果进行了分析和研究.     

Euler computation of the nearfield wake vortex of an aircraft in take-off configuration

Numerical simulations based on the three-dimensional Euler equations are used to investigate the predictive capability of an Euler code for calculations of the nearfield wake of a narrow-body airliner wind tunnel model in take-off configuration up to a half span behind the wing tip trailing edge. Simulation results on both structured and unstructured grids are presented. The results on the block-structured grid were obtained within the scope of the EU-project EUROWAKE. The simulation quality of the vortex formation and spatial development is analysed by comparison to wind tunnel measurements of the spanwise lift distribution available from the EU-project DUPRIN II and to experimental PIV data available from the EU-project EUROWAKE.     

三维Euler方程自适应八叉树结构各向异性直角网格算法

采用各向异性直角网格对相应的三维复杂流场进行了模拟 ,提出了物面三角化方法及网格与物面的相交判断过程。用以中心差分为基础的 Jameson的有限体积法对 Euler方程进行求解。计算结果表明 ,各向异性直角网格算法对相应流场进行模拟 ,计算时间较短而效果较好