俯仰振荡引起的二元高超声速进气道不起动/再起动特性

针对高超声速飞行器受到扰动后俯仰姿态可能会瞬时大幅度改变或振荡的问题,为了分析其对进气道起动/不起动特性的影响,在Ma=4.03条件下对二元高超声速进气道俯仰振荡的流场进行了非定常数值模拟研究。研究结果表明:大幅度俯仰振荡的二元高超声速进气道会出现不起动和再起动现象,进入不起动/再起动状态时的流场特征和性能参数均发生剧烈变化;不同折合频率俯仰动态变化时,进气道的性能参数及流场特征存在明显差异,气动性能存在迟滞现象;折合频率越大,进气道发生不起动时的攻角值越大,再起动时的攻角值越小;折合频率增加到一定程度时,进气道可能出现全程起动或不起动现象,初始攻角将是决定进气道处于全程起动或不起动状态的关键因素。     

流向磁场作用下二维磁流体槽道湍流直接模拟

对低磁雷诺数近似下流向磁场作用的二维磁流体槽道湍流进行了直接数值模拟(DNS,Direct Numerical Simulation),给出了 Re=10 000时不同磁相互作用数下二维磁流体槽道湍流的近壁速度分布、湍流脉动速度均方根、雷诺应力等统计量的变化,并与中性流体二维槽道湍流进行了比较.结果表明:流向磁场作用会导致对数区上移;雷诺应力最大值随磁相互作用数呈线性变化;随着磁相互作用数的增大,下壁面平均涡量的时间演化由拟周期性向周期性转变,且脉动周期逐渐增大并当流动层流化后下壁面平均涡量成为常值.      

对流层雷暴源所产生的柱形重力波

已有充分的证据表明, 大气对流层的雷暴迹象是大气中间层重力波活动的显著代表源.在雷暴迹象的上方, 通过火箭已观测到大气中间层出现的热效应, 也已通过雷 达探测到大气平流层出现的上行重力波, 从地面和卫星平台上观察到夜间气辉有序而 成环状的重力波波形. 所有这些实验结果都与位于观测点下方的雷暴活动有紧密联系. 此类雷暴通常主要集中在中国东部沿海以及地球其他沿海海湾地域. 关于此类雷暴对大 气中间层的影响尚未被充分研究和了解. 为能有效地探究其成因, 利用所开发的一个二维计算机数值模型模拟和研究大气对流层的雷暴源所引发的上行重力波, 进而揭示 此类重力波产生的基本物理机理, 及其在大气中间层的能量耗散. 通过模拟研究发现, 雷暴源可以大面积高强度地聚集和释放积雨云的能量, 当这种周期性的对流变化引发大气对流层的不稳定性后, 就会有圆柱体重力波的产生和传播.      

Numerical modeling and optimization of bundled fiber optic FTIR probes for spectroscopy of small targets

Fiber optic probes coupled to a Fourier Transform Infrared (FTIR) spectrometer allow sample analysis without the need to place the sample directly into the instrument. This is particularly attractive for spacecraft lander/rover applications where the additional mechanics, mass and cost associated with excavating and placing a small sample in a sample bay is not desirable. This paper uses numerical modeling to analyze fiber bundle array configurations and fiber numerical apertures to determine their light gathering efficiency and sample size as a function of probe-to-sample distance. The model was validated against experimental measurements for a selection of off the shelf fiber optic probes. The collection efficiency of a probe can be optimized by using a hexagonal packing method, and by using fibers with high numerical apertures. By increasing the numerical aperture of the fibers used in our model probes from 0.2 to 0.5, we were able to increase the collection efficiency of a seven fiber probe from 4.5% to 9.5%, with the peak collection efficiency corresponding to a spot size of diameter 3 mm.     

磁流体管流准一维模型的数值验证与性能分析

为了快速分析磁流体管流,采用涡轮喷气发动机模型对磁流体管流进行类比分析,引入变截面管流理论和修改过的多变效率表达式,推导了磁流体发电器和加速器的准一维模型,采用二维数值方法验证了准一维模型的可行性,利用该准一维模型对磁流体管流进行了性能分析,结果表明,多变效率的引入使准一维模型更加完整,当负载系数不变时,磁流体装置中的平均马赫数越低,其等熵效率越高。     

间隙高度对某涡轮叶尖间隙泄漏流影响的研究

利用三维湍流数值模拟方法模拟某跨声速高压涡轮流场,研究了转子叶尖间隙高度对涡轮转子间隙泄漏流动的影响,详细描述了间隙流动结构随间隙高度的变化情况。研究结果表明:随着间隙高度的增加,间隙泄漏流流量线性增加,并导致转子叶尖附近流体偏转角度减小。间隙高度每增加1%相对叶高,间隙泄漏流量占总流量的比例约增加2.1%。间隙高度变化对间隙内部流动的影响明显:当间隙比较小时,随着间隙高度的增加,分离泡的大小迅速增加,回流区域减小。随着间隙高度的增加,间隙泄漏涡尺度迅速增加,同时由于泄漏流中高速流体流量增加,通道涡远离叶片吸力面并向靠近机匣的方向移动,泄漏涡也向远离吸力面方向运动。     

“2010.8.24”飞机积冰过程分析

通过对2010年8月24日积冰过程的分析得出:这次积冰发生地点附近为垂直运动上升区,配合云中密集的冰水混合物,是造成这次飞机严重积冰的主要原因。积冰指数对于积冰强度和区域的预报具有指示意义。WRF中尺度数值模式选用Eta Ferrier微物理方案的模式结果中的云冰、云水以及总凝结混合比,这对于积冰强度和范围有一定的预报能力。     

薄膜区传热传质蒸发对毛细蒸发弯月面稳定性的影响

蒸发毛细弯液面的稳定性对两相传热装置的运行有着重要意义。由于弯液面的蒸发50%以上发生在一个微小区域,即延展薄膜区,可以预期不稳定性将最先产生于这一区域。为对蒸发弯月面的稳定特性进行机理性研究,首次提出一种将线性稳定性分析与薄膜区传热传质数值分析巧妙结合的研究方法。该方法通过传热传质数值分析获得能够反映弯月面稳定特性的特征膜厚和长度,并将其应用于线性稳定性判据以获得弯月面的临界过热度。分析表明,两种方法的结合可以有效揭示蒸发弯月面的稳定机制,特别是毛细管半径、弯液面温度和工质的热物性对弯月面稳定性的影响。与相关文献报道的实验结果比较,证实了该方法的有效性。     

质子地磁垂直截止刚度在2004年11月7-8日两个时刻处的数值模拟

利用单粒子轨道理论, 在T04和IGRF2000磁场模式建立的磁层模型基础上, 应用四阶龙格库塔方法, 模拟计算宇宙线带电质子在地球磁层中的运动以及沿天顶方向入射到达地球磁场内的某一特定位置, 得到了距地球表面450km高度处全球质子垂直截止刚度在2004年11月7-8日中两个时刻的计算值. 根据计算得到的684个不同位置处的截止刚度值, 分析了同一时刻地磁垂直截止刚度随磁纬和磁经的变化. 与此同时, 模拟计算了相同时刻下磁纬为30°, 磁经为0°, 45°, 90°, 135°, 180°, 225°, 270°, 315°处, 地磁有效垂直截止刚度随高度的变化情况. 结果表明, 在磁纬30°处, 其质子有效垂直截止刚度随距离地心高度的变化沿向阳处向磁尾处方向, 变化越来越缓慢.      

Dynamic simulation of aerial towed decoy system based on tension recurrence algorithm

Dynamic model of aerial towed decoy system is established and simulations are performed to research the dynamic characteristics of the system. Firstly, Kinetic equations based on spinor are built, where the cable is discretized into a number of rigid segments while the decoy is modeled as a rigid body hinged on the cable. Then tension recurrence algorithm is developed to improve computational efficiency, which makes it possible to predict the dynamic response of aerial towed decoy system rapidly and accurately. Subsequently, the efficiency and validity of this algorithm are verified by comparison with Kane's function and further validated by wind tunnel tests.Simulation results indicate that the distance between the towing point and the decoy's center of gravity is suggested to be 5%–20% of the length of decoy body to ensure the stability of system.In up-risen maneuver process, the value of angular velocity is recommended to be less than0.10 rad/s to protect the cable from the aircraft exhaust jet. During the turning movement of aircraft, the cable's extent of stretching outwards is proportional to the aircraft's angular velocity.Meanwhile, the decoy, aircraft and missile form a triangle, which promotes the decoy's performance.