膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

面向非结构混合网格高精度阻力预测的梯度求解方法

针对非结构混合网格的特点,通过改进传统的Green-Gauss梯度求解方法,提出了一种可提高非结构混合网格黏性计算精度的节点型Green-Gauss梯度求解方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析。改进后的梯度求解方法残差收敛更好,下降量级更多,阻力系数和试验吻合更好,激波区域压力分布和分离区域流场细节的模拟更精确,说明改进后的梯度求解方法有效提高了程序的鲁棒性和阻力预测精度,验证了方法的有效性。采用改进后的方法对第5届AIAA阻力预测研讨会的通用研究模型（CRM）进行了详细的模拟分析,结果表明：改进的梯度求解方法更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

T型尾翼颤振模型光学测量实验与弯扭特性解算

在FL-26跨声速风洞半模试验段进行了某高速飞机T型尾翼颤振模型的光学测量实验,并依据测量结果解算了尾翼颤振模型的弯扭特性。颤振模型表面用白色圆点进行标记,用于记录模型表面的位移变化,两台固定在风洞试验段上壁板观察孔旁肋板上的400万像素工业相机用来采集图像,采集到的图像通过自主开发的图像解算软件进行图像的识别与求解,计算出尾翼颤振模型表面标记点的三维坐标。模型表面标记点的三维坐标通过坐标变化转换到风洞气流坐标系中,利用不同时刻模型表面坐标的变化计算模型剖面扭角和弹性轴位移的分布。T型尾翼右平尾图像采集实验与弯扭特性计算结果表明,非接触光学测量技术可以用于高速颤振试验的定量分析中。     

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

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

振动环境中相机位置坐标与姿态角解算的实验研究

我国高速暂冲式风洞洞体振动较大,故振动环境中的相机位置坐标与姿态角解算技术对模型变形视频测量(VMD)至关重要.为此,建立振动模拟实验平台,研究相机位置坐标与姿态角求解方法,结果表明:基于蒙特-卡洛法的相机位置坐标与姿态角求解法最适于振动环境中VMD测量的相机位置坐标与姿态角求解.     

风洞模型静弹性变形对气动力影响研究

介绍了一种基于模型变形视频测量系统和计算空气动力学研究静弹性变形对气动力影响的方法.利用模型变形视频测量系统获取模型在气动载荷作用下的静弹性变形,驱动模型表面网格运动,得到模型变形后的表面CFD计算网格.CFD计算变形前后网格外形下的气动力,研究模型变形对模型气动特性的影响.对一大展弦比连接机翼的测量与计算结果进行了分析,分析结果表明:模型变形对升力系数影响最大发生在升力线性变化的最大迎角附近,模型变形对阻力系数影响最大发生在失速迎角附近,模型静弹性变形对气动力的最大影响量远远超出风洞测力实验的精度指标,因此开展风洞模型静弹性变形影响研究与修正是十分必要的.     

2．4m跨声速风洞的模型位移视频测量精度研究

风洞试验中模型的位置和变形测量对试验数据精准度至关重要。为此,创建2．4m跨声速风洞的模型位移视频测量系统,提出度量其测量误差的方法,并实验研究其测量精度。研究发现,试验中的振动对测量精度影响极大,采用振动环境中相机位、姿解算方法后,试验段底部的编码标记点的测量误差从22．80-48．48mm降至0．03～0．64mm。     

火箭子级垂直回收布局气动特性及发动机喷管影响

垂直回收可重复使用运载火箭是运载火箭发展的一个重要方向,大长细比的火箭子级垂直再入过程属于典型的非规则钝头体绕流,与传统低阻力流线体飞行器气动特性差异较大。采用风洞试验辅以数值仿真分析的方法,对基于栅格舵的火箭子级垂直回收构型基本气动特性和非规则钝头体绕流情况进行了研究,获得了发动机外露喷管和栅格舵对火箭子级垂直回收气动特性的影响规律,给出了火箭子级垂直回收布局设计建议。结果表明：火箭子级倒飞状态下肩部区域会在小迎角下产生大分离流动,外露发动机喷管左右两侧诱导出较强的分离涡结构,与火箭尾翼、肩部大分离流动相互作用;垂直回收构型在超声速下阻力会一直处于较高的水平,不同马赫数下压心移动量较大,倒飞时发动机外露喷管会产生较大的干扰静不稳定力矩,其量值与栅格舵提供的静稳定控制力矩基本相当,在火箭子级垂直回收方案设计时需要引起注意。     

基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用

为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进一步提高GMRES的计算效率,发展了方程组右端项残值雅可比精确计算方法,针对无黏通量Roe格式和黏性通量BR2(Bassi Rebay 2)黏性计算方法,分别解析给出其对守恒变量多项式自由度的雅可比矩阵。基于建立的方法首先采用NACA0012翼型研究了GMRES的重启次数及收敛参数对方法收敛性影响,然后采用无黏及黏性算例对比研究了基于雅可比矩阵不同计算方法的GMRES计算效率,同时对比研究了雅可比矩阵完全近似求解下GMRES和LU-SGS(Lower Upper-Symmetric Gauss-Seidel)的计算效率。结果表明,建立的基于右端项残值雅可比矩阵精确求解的GMRES方法能够大幅提高不同精度DG方法的CFL(CourantFriedrichs-Lewy)数,相比前面提到的其它方法具有更高的计算效率,其收敛速度实现量级以上的提高。     

激波/湍流边界层干扰物面剪切应力统计特性

为了揭示激波/湍流边界层干扰区内物面剪切应力统计特性的演化规律，采用直接数值模拟方法对来流马赫数2.9、12°激波角的入射激波与平板湍流边界层相互作用问题进行了研究。通过与风洞试验数据的比较分析，验证了计算结果的可靠性。系统地探究了干扰区内物面剪切应力的典型特征，如预乘谱、概率密度分布和相干结构等。研究结果表明，分离激波的低频振荡运动对流向及展向分量的预乘谱均没有实质影响，其脉动仍以高频特征为主，低频能量变化较小。干扰区内流向剪切应力概率密度函数分布变化剧烈，分离泡内对数正态分布规律不再适用，而展向剪切应力在干扰区内与正态分布较为接近。相较于上游湍流边界层，分离泡内物面剪切应力矢量夹角与幅值的联合概率密度变化显著，峰值概率降低，峰值范围增大。此外，流向剪切应力脉动场的本征正交分解分析指出，主能量模态与分离激波的低频振荡以及下游再附边界层内的Görtler-like流向涡结构密切相关。