带压缩性修正的离散涡方法在结冰数值模拟中的应用

离散涡方法（DVM）是一种无网格的涡运动算法,适用于解决易产生分离流的非定常问题,将其应用于结冰过程中的流场求解,在有效模拟分离流动的同时能避免冰形尖角对网格质量的影响。但离散涡方法基于不可压 N-S 方程,无法应用于预测可压缩流动下的结冰过程。本文在离散涡方法的基础上添加普朗特—格劳尔特压缩性修正,进行基于离散涡方法的可压缩流动下的数值模拟,并将其应用于翼型结冰预测;对流场分布、结冰冰形和结冰模型计算过程等仿真结果进行对比和分析。结果表明：引入压缩性修正后的离散涡方法能较好地模拟可压缩流动,与实验值相比,基于该方法得到的结冰数值模拟结果符合良好,对结冰数值模拟在工程上的应用提供了一定的参考。     

波瓣数对波瓣S型混合二元喷管气动热力性能影响

依据某型涡扇发动机波瓣S型混合二元喷管,保持波瓣混合器长度、内扩张角、外扩张角以及宽高比不变,依次取波瓣混合器波瓣数为12,14,16,18,20,建立了一组具有不同波瓣数的波瓣S型混合二元喷管模型.采用经过验证的CFD方法,研究了波瓣数对波瓣S型混合二元喷管气动热力性能的影响规律.结果表明:在波瓣尾缘截面至第1个S弯截面区域,波瓣数对流体混合程度产生很大影响,并且热混合效率近乎随波瓣数增加而增加.在第1个S弯截面至波瓣S型混合二元喷管出口截面区域,波瓣数为16的波瓣S型混合二元喷管模型的总压恢复系数始终最低,其余模型的总压恢复系数以及热混合效率没有明显差别.在波瓣S型混合二元喷管出口截面上,波瓣数为16的波瓣S型混合二元喷管模型的热混合效率最高,达到0.850,然而其总压恢复系数相对于该截面上最高值下降了0.289%.此外,波瓣S型混合二元喷管的渐缩型流道能够提高流向涡强迫混合效果,但同时也加速流向涡的耗散速率.      

基于POD曲线的预警机雷达检飞概率模型

针对预警机雷达检飞试验中距离取样间隔及检飞航线边界条件的确定问题,以保证发现概率(POD)曲线的覆盖率与精确性为目标,提出基于POD曲线的预警机雷达检飞概率模型.首先,根据POD曲线的覆盖率和精确性要求确定检飞试验拟合曲线所需的采样点数.然后,根据POD曲线的经验分布函数建立检飞试验中最小探测距离和最大探测距离的计算方法,并由此确定检飞航线上距离取样间隔的长度.最后,根据发现概率的统计特性得到距离取样间隔内采样点数以及航次数等试验参数,为雷达检飞试验设计与预警性能评定提供定量依据.      

基于涡系相交不稳定性的飞机尾流控制方法

在飞机飞行的过程中尾涡会伴随着升力产生,威胁后机的飞行安全.在简化机翼模型上添加扰流片,通过一个矩形翼以引入一个与主翼尾涡大小不同、方向相反的小涡,构建尾流自消散四涡系统,以期诱发尾涡的Rayleigh-Ludwig相交不稳定性.通过改变扰流片的大小形状,调整模型的攻角和拖曳速度,采用粒子图像速度场仪测量系统定量研究在低雷诺数下单主翼尾涡发展特性以及双涡相互作用特性.研究表明:在未添加扰流片时,尾涡环量在45个翼展内相对于初始环量基本保持不变;在添加扰流片的情况下尾涡的环量衰减可以达到35%~55%,而未添加的基本翼型的尾涡的环量则几乎保持不变,这说明添加适当的扰流片能诱发尾涡的Rayleigh-Ludwig相交不稳定性,加速尾涡的消散,当小涡和主涡的初始环量比为-0.489、初始距离比为0.5时,45个翼展范围内,尾涡环量衰减55.9%.本文系统性的实验结果可以为低尾流机翼的设计提供参考依据.      

高马赫数层流摩阻数值计算精度

飞行器在高空高速飞行时黏性效应显著,摩阻预示精度对飞行器的关键气动性能意义重大。目前,摩阻预示主要依赖数值计算,但高马赫数层流摩擦阻力计算与试验测量结果仍存在差距。以具有高速飞行器典型部件特征的球锥、三角翼为对象,结合风洞试验摩阻测量结果,使用国家数值风洞（NNW）数值计算软件和自研CFD程序研究了数值计算中影响摩阻计算精度的格式数值耗散及壁面温度边界条件等重要因素。研究结果表明：数值耗散越小,表面摩阻的计算精度越高;在速度较低的边界层近壁区内关闭熵修正,将有助于提高表面摩阻的预示精度。此外,在高马赫数流动问题的数值模拟中,壁面温度条件对表面摩阻计算同样具有重要影响。最后,基于分析结果和工程需要提出了对高精度摩阻数值预示的研究需求。     

火星探测器降落伞拉直过程中的“绳帆”现象研究

针对火星探测器降落伞在拉直过程中出现的“绳帆”现象,以及火星探测器降落伞开伞前初始参数和大气密度与地球环境下的差异,建立了火星探测器降落伞拉直过程的数学模型,研究了火星环境下,伞包弹射速度、开伞前进入器的攻角、开伞马赫数以及大气密度对“绳帆”现象的影响。研究结果表明,选择较大的伞包弹射速度,并将开伞前进入器的攻角严格限定在较小范围内,将有利于避免或降低“绳帆”现象的发生。这一研究结果可为我国实施火星探测时减速着陆系统的设计分析提供一定参考。     

Analytical-HZETRN model for rapid assessment of active magnetic radiation shielding

The use of active radiation shielding designs has the potential to reduce the radiation exposure received by astronauts on deep-space missions at a significantly lower mass penalty than designs utilizing only passive shielding. Unfortunately, the determination of the radiation exposure inside these shielded environments often involves lengthy and computationally intensive Monte Carlo analysis. In order to evaluate the large trade space of design parameters associated with a magnetic radiation shield design, an analytical model was developed for the determination of flux inside a solenoid magnetic field due to the Galactic Cosmic Radiation (GCR) radiation environment. This analytical model was then coupled with NASA’s radiation transport code, HZETRN, to account for the effects of passive/structural shielding mass. The resulting model can rapidly obtain results for a given configuration and can therefore be used to analyze an entire trade space of potential variables in less time than is required for even a single Monte Carlo run. Analyzing this trade space for a solenoid magnetic shield design indicates that active shield bending powers greater than ∼15 Tm and passive/structural shielding thicknesses greater than 40 g/cm² have a limited impact on reducing dose equivalent values. Also, it is shown that higher magnetic field strengths are more effective than thicker magnetic fields at reducing dose equivalent.     

An empirical approach to predicting the key parameters for a sunspot number cycle

The common methodologies used to predict the smooth sunspot number (SSN) at peak (Rmax) and the rise time (Tr) for a cycle are noted. The estimates based on geomagnetic precursors give the best prediction of Rmax for five SSN cycles (20–24). In particular, an empirical technique invoking three-cycle quasi-periodicity (TCQP) in Ap index has made accurate predictions of Rmax and Tr for two consecutive SSN cycles (23 and 24). The dynamo theories are unable to account for TCQP. If it endures in the 21st century the Sun shall enter a Dalton-like grand minimum. It was a period of global cooling. The current status of the ascending phase of cycle 24 is described and the delayed reversal of the solar polar field reversal in the southern hemisphere in September 2013 is noted.     

铝合金薄板零件橡皮囊成形屈曲分析与实验

为了抑制铝合金薄板类零件在橡皮囊成形后出现屈曲的缺陷,采用薄板屈曲微分方程对零件成形后的屈曲进行理论分析,并通过橡皮囊成形实验研究自然时效时间和成形压力对屈曲缺陷的影响。理论分析结果表明：零件法兰屈曲波数与其内外径比有关,通过改变法兰宽度进行实验验证,实验结果与理论分析结果相符;而成形实验结果表明,屈曲缺陷高度随板料自然时效时间的增加而增加,随成形压力的增加而减少。因此,理论分析结果有效,并且可通过减少淬火后成形时间和增加成形压力来抑制屈曲缺陷。     

微分雷诺应力模型在激波分离流中的应用

由于Boussinesq涡粘性假设的湍流模型在一些流动中如激波分离流中预测失真,本文介绍一种基于完全雷诺应力方程的微分雷诺应力模型(DRSM),应用压缩性修正和改进耗散率方程构造出应用于可压流体的可实现DRSM,SSG压应力模型测试了两组参数。将其应用于ONERA M6机翼在跨声速大攻角的计算,结果表明DRSM能更好地模拟激波位置和流动分离。