高超声速进气道快速破膜开启的流动特性

进气道处于起动状态是保证超燃冲压发动机正常工作的前提,进气道帽罩快速开启时的非定常效应可以有效提高进气道的起动能力。采用非定常数值计算深入研究了唇口帽罩不同安装位置开启时的非定常效应对进气道起动过程的影响,分析了不同帽罩安装位置开启时进气道流场的演化过程,并揭示了喉道分离泡的形成机理。研究结果表明,当帽罩上游不存在分离泡时,破膜非定常激波在压缩面反射,与上游复杂波系作用形成沿壁面的低速流,在唇口激波作用下在喉道形成分离泡。帽罩安装靠近唇口可通过缩短激波/边界层作用距离减小低速流动区范围,进气道临界起动内压比随之增大;而当帽罩上游出现大分离泡时,分离泡会先演变为低速流,之后在唇口激波作用下重新聚集形成大尺度分离,进气道临界起动内压比显著降低。     

连续旋转爆轰发动机参数特性的三维数值模拟

采用一步化学反应模型,基于任意坐标系下Euler控制方程,对连续旋转爆轰发动机(CRDE)进行了三维数值模拟研究, 详细分析了来流总压对CRDE参数特性的影响.研究发现:随着来流总压的增大,爆轰波峰值随之增大,但是燃烧室头部Laval型喷注段的壅塞比基本不变.不同来流总压下,可燃气体在燃烧室头部均以亚声速入射.随着来流总压的增大,燃烧室内流场的平均压强增大,但是其平均轴向流速基本不变.CRDE的流量、推力和比冲均随着来流总压的增加而变大;但是,不同来流总压下,CRDE的流量与Laval型喷注段最大流量之比不变,并且此比值约等Laval型喷注段的壅塞比.Laval型喷注段最大流量与来流总压成正比,因此以上分析从理论上进一步解释了CRDE的流量与来流总压成正比的原因.      

喷注壁面楔状流层流边界层速度与温度相似解

建立了壁面有喷注的楔状流层流边界层冷却数学模型,求解经相似变换得到的描述无量纲流函数和无量纲温度的常微分方程,获得了楔状流层流边界层无量纲速度和温度的相似解,给出了考虑壁面喷注边界条件的形式简洁、物理意义明确的无量纲流函数拟合式;用Runge-Kutta法求解无量纲速度与温度的常微分方程,获得了壁面有喷注的楔状流层流边界层无量纲速度和温度随楔形角、吹风参数、冷却介质温度的变化规律.通过计算0°,18°和36°楔形角的楔状流层流边界冷却速度与温度分布得到以下结论:楔状流速度边界层和温度边界层厚度都随楔形角的增大而变小,随着吹风参数的增大而增大;当冷却介质温度越低、吹风参数越小、楔形角越大时,靠近壁面处楔状流层流边界层内温度梯度越大.      

一种稀薄两相流动的数值模拟方法

从描述粒子运动的微观层次出发,采用双向耦合技术,建立了一种适用于稀薄条件下两相流动的DSMC数值模拟方法。对相间相互作用进行解耦处理,实现了气固两相间动量和能量相互作用的模拟。相关算例验证了方法的可靠性。结果表明:在大样本统计平均条件下,该模拟方法能够保证相间动量和能量交换时守恒,固体颗粒尺寸、气体稀薄程度等条件对相间的相互作用有很大的影响。     

喷注壁面外掠平板层流边界层速度和温度的HPM-Padé解析解

运用HPM-Padé(同伦摄动-帕德逼近)法导出了当喷注流体速度与距平板前缘距离的平方根成反比时,不可压缩常物性流体外掠平板层流边界层内无量纲速度和无量纲温度解析表达式,无量纲流函数1阶导数的HPM-Padé解析解与4阶龙格库塔法的数值解结果一致.利用HPM-Padé解析解研究了喷注系数和流体Prandtl数对速度和温度分布的影响.结果表明:Prandtl数越大,温度边界层越薄,壁面温度梯度越大;喷注系数越大,速度边界层厚度和温度边界层厚度越大,壁面速度梯度和壁面温度梯度越小,壁面喷注对平板有隔热作用;当喷注系数为0.619时,壁面的速度梯度和温度梯度为零,高温来流向壁面的传热被喷注流体完全阻隔.     

低雷诺数下压气机的二次流旋涡结构

为了探究高空低雷诺数条件下跨声速压气机的流动规律,对NASA Rotor37进行单通道数值模拟,探索其在低雷诺数进气条件下二次流的旋涡结构.研究发现:马蹄涡压力面分支诱发压力面角区诱导涡,壁角涡形成了顺流和逆流的两段式结构,脱落涡由叶根角区发展起来后不断从尾缘脱落,泄漏涡近失速点仅局部破裂不是失稳触发的主要原因.通道中的激波系诱发了吸力面和压力面的两个径向涡,压力面径向涡构成闭合的气泡式分离,吸力面径向涡在叶顶的破碎诱导产生分离涡,触发了低雷诺数下压气机的失稳.流场旋涡结构由马蹄涡、壁角涡、径向涡、泄漏涡、分离涡、脱落涡6个大尺度旋涡以及其他小尺度旋涡组成.      

非线性负阻尼极限环型角位移振荡运动定量分析

从某飞行器两次飞行记录的俯仰角θ、偏航角Ψ观测值出发,分别采用对Van der Pol方程描述的θ和Ψ的角位移振荡运动微分方程进行气动参数辨识,和对观测值θ-狋、Ψ-t曲线的外包络线进行参数拟合两种分析方法,对该飞行器的角位移振荡运动特性进行了定量分析。两种分析方法取得了较为一致的结果,并证明该飞行器飞行中出现的锥形振荡运动是典型的非线性负阻尼极限环型的振荡运动,获得了非线性负阻尼极限环型气动阻尼力矩的典型表达式。分析结果表明非线性负阻尼极限环类型的气动阻尼能够导致飞行器出现动不稳定。     

涡轮叶片叶顶间隙变化减敏研究

采用数值方法联合标准k-ω两方程湍流模型求解雷诺平均Navier-Stokes方程组,研究了平顶和分别采用被动、主动以及复合间隙控制方法的凹槽顶、平顶喷气、凹槽顶喷气4种不同叶顶结构下涡轮性能对叶顶间隙变化的敏感性,分析了间隙变化对间隙流场结构及损失的影响,研究了涡轮性能随叶顶间隙变化的规律,并对采用不同间隙控制方法的叶顶结构所带来的效率收益对攻角变化敏感性进行了评估.结果表明:在不同间隙控制方法下叶顶泄漏涡强度随间隙的增加速度变缓,尤其在复合间隙控制方法下更为明显,从而使得在复合间隙控制方法下涡轮性能对间隙变化的敏感度最低;仅仅涡轮叶片前缘区域泄漏流动对攻角变化比较敏感,使得凹槽顶结构具有最佳的攻角变化适应性.      

球面收敛二元喷管电磁散射特性

采用迭代物理光学法和等效边缘电磁流法,对球面收敛二元喷管电磁散射特性进行研究.计算分析了该喷管电磁散射的构成及散射机理、散射场的空间分布特性、散射的频率特性,分别计算和对比分析了电磁波在俯仰平面和偏航平面入射时喷管在中间和最大两种状态下的电磁散射情况.研究表明该喷管在垂直极化下边缘绕射场对散射总场的贡献在大角度下体现.电磁波入射频率越高,该喷管后向散射越强烈,雷达散射截面分布对姿态角度越敏感.      

基于边界涡量流的二维叶型优化建模与求解

针对二维叶片整体造型优化问题,提出一种基于边界涡量流的二维叶型优化模型,并对该模型进行求解。建立以含参数的Bezier曲线定义的二维叶型型线为边界、以边界涡量流正峰值的最小化为目标函数的优化模型,通过分析边界涡量流的生成机制,把叶型型线边界信息代入边界涡量流中,将边界涡量流用叶型型线参数和型线上的气动参数表示,进而获得新的目标函数表达式。该目标函数表达式可直接对叶型型线参数求导,仅使用简单的定步长梯度算法即可对优化模型求解。结果表明:优化后的叶型型线上的边界涡量流的正峰值显著降低。