高超声速钝头体边界层转捩试验

为研究高超声钝头体边界层转捩以及头部钝度对转捩的影响,在FD-14和FD-14A两座激波风洞中开展了热流、压力扰动和高速纹影显示等综合测量。试验结果表明,转捩雷诺数关于钝度雷诺数的变化显示出转捩反转的趋势。压力扰动的功率谱密度（PSD）分析结果以流向离散分布云图形式显示,边界层高速纹影图像显示了第二模态波的发展、湍流的生成和熵层对边界层结构的显著影响。大头部钝度带来的强熵梯度熵层流动对边界层压力扰动频谱特性和流动结构影响显著,在转捩反转机理中起到重要作用。此外,马赫数对转捩的影响不容忽视。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

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

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

“蚌式”进气道附面层扫除特性风洞试验

以某低、高速风洞为平台,设计搭建了“蚌式”进气道附面层扫除特性测量试验系统,进行了不同流量系数和来流马赫数下进气道鼓包表面附面层扫除特性的风洞试验,通过对试验数据的整理、计算和对比分析同型号的飞行试验结果,研究了“蚌式”进气道鼓包表面附面层扫除特性。研究结果表明:在相同的来流马赫数下,随着流量系数的增大,鼓包表面附面层的扫除能力逐渐减弱;在亚声速工况的绝大多数流量范围内,鼓包表面压力系数沿鼓包中心线对称分布、压力梯度变化明显,且在不同截面沿主流方向具有增大的特征,鼓包构型对附面层扫除效果较强;超声速工况下具有明显附面层扫除能力的流量范围明显小于亚声速工况,进气道唇口形成的弓形激波是影响鼓包表面不同位置压力梯度变化的主要因素,进而决定着附面层扫除特性。在接近来流马赫数18及以上飞行工况下,附面层的扫除能力减弱,附面层分离加强,进而会造成较大的进气压力损失和畸变。     

带非线性支撑的转子有限元模型求解方法

用数值方法研究了非线性支撑的柔性转子系统的动学行为,提出了一种将有限元与非线性支撑结合的模型和求解方法。利用有限元法(FEM)构建转轴和转盘部分的模型,通过矩阵进行组合;利用离散元方法对包含滚动轴承和挤压油膜阻尼器(SFD)的支撑部分进行建模,此部分包含4个单元,分别为轴承内圈、外圈、SFD内圈和支撑鼠笼。有限元部分和离散元部分通过轴端节点相连,仿真过程中轴端位移传递给非线性支撑部分,支撑部分通过位移计算得到的非线性力反过来作用于有限元转子轴端部分。为了耦合求解有限元转子和非线性支撑组成的数学模型,提出了一种综合的迭代求解方法,克服传统的有限元求解方法对轴端隐性非线性支撑的求解局限性。由于转轴部分采用了Timoshenko梁单元建模,对比与简单转子模型,可以考虑陀螺力矩和轴的柔性特征,更能体现非线性支撑对振动真实影响。在建立的20个轴单元的有限元转子模型中,非线性响应更多体现在靠近非线性支撑的节点1和节点21处,响应频谱中靠近轴端的节点能体现出滚动轴承的2倍和3倍变柔振动频率。     

考虑可压缩与热传导的壁面函数边界条件及其应用

考虑到可压缩和热传导效应的壁面函数边界条件,被耦合到了采用k-ω两方程湍流模型、用有限体积法求解N-S方程的程序中。壁面函数基于耦合的速度和温度型,并且在边界层内的粘性子区和对数区内一致有效。引入壁面函数边界条件后,通过算例验证在y <100的范围内,得到的物面压力、摩阻、热流与实验结果比较,结果可靠。而无壁面函数边界条件时,要得到相同精度的结果,要求y ≈1。壁面函数的引入,为工程上准确预测飞行器在湍流流动中表面受力与气动热提供了保障。     

日照边缘区域电离层对耀斑的响应特点研究

利用MSIS模型和背景太阳辐射谱模型,在一定大耀斑辐射谱假设的前提下,计算了耀斑期间日照边缘区域的电子产生率,分析了这一区域电离层电子密度的变化特点.结果表明,大耀斑期间在日照边缘区域,甚至大于太阳天顶角90°的区域都有明显的电子产生率的增加.从不同太阳天顶角处的电子产生率剖面的形态来看,随着天顶角的增加最大电离率减少,但高度增加.计算还显示了在太阳天顶角小于90°的区域内电子产生率的垂直分布有明显的双峰结构,这种结构对应着电离层的E区和F区,但在天顶角大于90°区域,F区的电子产生率要大得多.考虑到离子和电子的复合过程,这一区域的总电子含量的增加主要产生在高F区.      

A swept fin-induced flow field with different height mounting gaps

In order to apply the air fin successfully and ensure the maneuverability of hypersonic vehicle, a key problem to be studied urgently is the heat flux brought by the fin mounting gap. The appearance of mounting gap and fin shaft can induce many complex flow structures which need more attentions to be investigated. Under Ma 6, Nano-tracer-based Planar Laser Scattering (NPLS) and Temperature Sensitive Paints (TSP) were applied to visualize and measure transient flow structures and heat flux distribution of a swept fin-induced flow field with different height mounting gaps. Complementarily, Reynolds-averaged N-S equations were solved with k-ω SST turbulent model. The heat flux distribution results of numerical simulation and TSP observed the change of high heat flux region with different mounting gap, both in position and magnitude. The streamlines based on Computational Fluid Dynamics (CFD) and flow visualization results obtained by NPLS revealed the cause of high heat flux region. The high heat flux region in this flow field is mainly related to the reattachment of vortex and flow stagnation. The increase of gap height can lead to stronger gap overflow and shaft-induced horseshoe vortex, which are source of the high heat flux around the fin. The case with the highest mounting gap (4 mm) en-counters the most severe aerodynamic heating, both on the surface of fin and plate. Thus, under the premise of ensuring the flexibility of the fin, the gap should be set as small as possible.     

潜入式喷管喉衬界面间隙优化设计

为了解决固体火箭发动机潜入式喷管喉衬在热载荷与内压联合作用下界面间隙优化问题,数值模拟方法求解了喷管热结构响应.采用流体仿真软件,计算了潜入式喷管流动过程的稳态流场,获取了高温燃气流动参数与固体壁面对流换热系数.并采用三维有限元结构计算平台,编写了非均布壁面压力载荷与非均布对流换热系数子程序,求解了喉衬前后搭接界面间隙...     

模态综合的直接变换法

基于固定界面模态综合的Craig-Bampton方法,提出了模态综合的直接变换方法,该法在形式上更简洁,适用于大规模动力学系统的模态分析和动力学模型减缩。对采用直接变换法模态综合过程的计算误差做了分析,给出了计算误差界,阐明了模态综合后总体结构的固有频率与子结构保留主模态个数之间的关系。针对子结构界面自由度太大的情况,提出了直接变换法的分部算法,不需组集系统总刚度、总质量矩阵,仅在子结构的基础上完成系统的减缩。计算实例表明该算法的有效性及计算的稳定性。