Global Response of Martian Plasma Environment to an Interplanetary Structure: From Ena and Plasma Observations at Mars

As a part of the global plasma environment study of Mars and its response to the solar wind, we have analyzed a peculiar case of the subsolar energetic neutral atom (ENA) jet observed on June 7, 2004 by the Neutral Particle Detector (NPD) on board the Mars Express satellite. The “subsolar ENA jet” is generated by the interaction between the solar wind and the Martian exosphere, and is one of the most intense sources of ENA flux observed in the vicinity of Mars. On June 7, 2004 (orbit 485 of Mars Express), the NPD observed a very intense subsolar ENA jet, which then abruptly decreased within ∼10 sec followed by quasi-periodic (∼1 min) flux variations. Simultaneously, the plasma sensors detected a solar wind structure, which was most likely an interplanetary shock surface. The abrupt decrease of the ENA flux and the quasi-periodic flux variations can be understood in the framework of the global response of the Martian plasma obstacle to the interplanetary shock. The generation region of the subsolar ENA jet was pushed towards the planet by the interplanetary shock; and therefore, Mars Express went out of the ENA jet region. Associated global vibrations of the Martian plasma obstacle may have been the cause of the quasi-periodic flux variations of the ENA flux at the spacecraft location.     

激光冲击残余应力场的有限元分析

建立模拟激光冲击残余应力场的非线性弹塑性有限元模型,预测AISI 304奥氏体不锈钢厚、薄靶材经激光冲击后表面以及内部的残余应力场分布.模拟过程中,考虑材料高应变率(106s-1)下的力学性能和激光诱导冲击波的精确加载,分析激光功率密度、激光光斑尺寸、激光脉宽、激光冲击次数、激光单面与双面冲击等激光冲击参数以及初始残余应力等因素对不锈钢靶材残余应力分布的影响.数值模拟结果与X射线衍射法测量值吻合较好.在有初始残余拉应力(250MPa)存在的情况下,激光冲击仍能使304奥氏体不锈钢靶材形成残余压应力层,说明激光冲击工艺可提高奥氏体不锈钢焊接构件的抗应力腐蚀开裂能力.     

高速列车进入带缓冲结构隧道的压力变化研究(Ⅰ)

给出了列车穿越带有缓冲结构的隧道压力变化的三维粘性流场数值模拟过程,控制方程为三维粘性、可压缩、非定常流的N-S方程,空间离散采用了中心有限体积法格式,时间采用预处理二阶精度多步后差分格式进行离散,对列车与隧道之间的相对运动采用移动网格技术处理。对不同的缓冲结构缓解隧道内瞬变压力及压力梯度的作用进行了研究。研究结果表明,缓冲结构的设置能够有效地降低隧道内的压力和压力梯度的最大值,其原因在于缓冲结构延长了压缩波压力上升的时间,降低列车突入隧道时所形成的最大压力梯度;另一方面由于压缩波在缓冲结构和列车、隧道之间多次的反射,也降低了压力峰值。     

冲击响应谱绝对校准研究中冲击响应谱的计算方法

介绍了应用于冲击响应谱绝对校准中的冲击响应谱的计算方法,该方法的实现基于MATLAB编程语言,即通过MATLAB中的仿真直接求解单自由度二阶微分方程的数值解以获得冲击响应谱的绝对复现.     

CE/SE法模拟爆震波点火及传播过程

运用时空守恒元和解元(CE/SE)法模拟一维爆震波的形成及传播过程,分析产生爆震波的不同机理和DDT(爆燃-爆震转捩)过程的影响因素.采取激波点火和热点点火两种点火方式,分别对应SDT(激波-爆震转捩)和DDT过程,热点点火初始场采用线性温度分布和阶梯性温度分布,整个管道均为常压.用隐式梯形方法实现刚性源项积分.结果表明,SDT和DDT对应于不同的爆震波机理,并且DDT受点火区不同参数的影响,为研究爆震波点火和DDT过程提供了理论依据.     

等离子体涡电磁散射特性及隐身性能

 针对非均匀等离子体在飞行器隐身中的应用,采用分段线性电流密度递归卷积时域有限差分（PLJERC-FDTD）方法计算等离子体涡及涡串电磁散射特性,分析等离子体涡对飞行器隐身性能影响。计算表明,等离子体涡在很大频率区间对电磁波吸收效果显著,RCS降低很大,具有明显的隐身效果。等离子体涡表现出一定规律性的极化特性,对L,S和C波段电磁波具有不同的吸收、反射特性。     

基于SIMULINK的超声速进气道仿真方法

将进气道一维流动的偏微分方程离散为一系列可积分的线性常微分方程,基于此创建了超声速进气道仿真模型,提出了一种基于SIMULINK的进气道仿真方法,模拟了来流参数和燃烧室反压扰动时流场内各个参数的变化规律。对一个一维变几何超声速进气道进行仿真,获得了结尾正激波的稳态特性及动态特性。该模型定义并输出了激波位置,有效地将一维超声速进气道的分布参数问题转换为零维问题,解决了系统分析与控制系统设计中无法实现集总参数的困难。     

旋转爆轰流场的数值模拟

基于带化学反应的二维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对充有当量比的氢气-空气预混气和空气的环形旋转爆轰流场,从点火燃烧到发展成旋转爆轰的过程进行了数值模拟。根据数值结果分析了波后流场中爆轰产物受激波、高温和离心力等作用而挤向外壁,形成有利于充入燃料,实现持续稳定旋转爆轰的流场特征。还讨论流场中爆轰波、激波与间断面和内外壁面反射或折射,从而形成多个激波相交的波系特征。为认识和理解旋转爆轰流场,开展旋转爆轰的实验研究等具有现实的指导意义。     

气流出口型面对凹面腔内激波聚焦起爆的影响

为提高两级脉冲爆震发动机工作的可靠性,探寻强化凹面腔内激波聚焦起爆爆震波的方法,通过数值模拟,采用氢气作为燃料、空气作为氧化剂,探讨了3种凹面腔气流出口型面及出口面积对激波聚焦起爆的影响。结果表明:采用垂直出口壁面有助于提高凹面腔内的激波聚焦起爆压力和温度,使起爆时刻提前,随着出口面积的减小,激波聚焦峰值压力有较大的提高;但效果不如垂直壁面,随着倾斜壁面出口面积减小,起爆点的压力和温度只有小幅度增加,起爆时刻基本无变化;采用弧形出口壁面,不能对向凹面腔底部传播的激波产生反射强化作用,不能有效起爆爆震波。对比3种出口壁面型面,垂直壁面能更有利于凹面腔内的能量聚集,强化激波聚焦起爆爆震波。      

异形腔体压电声衬声学性能

为提高压电声衬对低频噪声的抑制范围,对声衬腔体进行结构优化。利用平面波理论构建了两种曲线管道的声学物理模型,并分别建立了两种模型的传递矩阵,以此作为异形腔体亥姆霍兹共振器传递损失计算的理论依据,并通过仿真验证其正确性。结合压电振子的形变对声衬进行有限元仿真分析,结果表明:在压电振子施加500V驱动电压时,两种声衬频率偏移量分别为115Hz和120Hz。与圆柱形腔体声衬进行对比结果表明:在相同腔体厚度范围内,由曲率越大的曲线所生成的腔体,在相同驱动电压条件下,频率变化率越高,这为今后对声衬腔体结构优化提供一种有效的依据。