OA212翼型主动流动控制的数值模拟研究

采用数值模拟的方法,探讨了基于零质量射流的主动流动控制技术对OA212旋翼翼型动态失速的控制效果和控制特性.以积分形式雷诺平均Navier-Stokes(N-S)方程为控制方程,采用格心有限体积法进行求解.空间离散采用AUSM~+-up格式,时间推进采用含牛顿型LU-SGS子迭代的全隐式双时间法,且引入了预处理方法和多重网格方法加速收敛.通过在喷口上施加非定常边界条件来模拟射流对翼型绕流的影响.研究了不同类型射流、不同位置射流以及不同控制参数(频率、相位、偏角、动量系数等)对动态失速控制效果的影响.研究表明:零质量射流和传统的定常射流均可减小动态失速迟滞环的回线面积,但在提高最大升力方面零质量射流明显优于定常射流;在12%c和62%c处施加组合零质量射流的控制效果最为明显.     

基于控制理论的旋翼翼型优化设计

基于粘性Navier-Stokes方程和最优控制理论原理,研究了气动/几何约束条件下多设计变量的旋翼翼型气动优化设计问题。根据给定目标函数的表达形式,在计算坐标系下推导出了相应的伴随方程和边界条件,以及梯度方程的数学表达式,并对流动控制方程和伴随方程进行了有效数值求解。综合流动方程、伴随方程、目标函数敏感性导数和优化算法,发展了一种高效旋翼翼型气动优化方法。通过典型旋翼翼型的算例验证,表明方法的有效性。     

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

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

基于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模型准确预测了连续流、滑移流和过渡流域的物面压力、摩阻与热流系数。     

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

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

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

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

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

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

基于CFD的静气动弹性优化设计方法

以三维Navier-Stokes方程为控制方程,数值计算弹性机翼跨音速气动力,耦合结构静平衡方程,研究弹性机翼的静气动弹性变形和真实载荷分布,并在此基础上,对机翼的型架外形进行基于静气动弹性的多学科的优化设计(MDO),以及对机翼进行基于总升力不变的飞行姿态确定,以满足弹性机翼在飞行时的品质设计要求.以某后掠机翼为例,设计结果达到预期目标.