液氢液氧发动机喷管化学反应流数值模拟

用时间相关的半隐格式求解液氢液氧发动机喷管中的化学反应流动，考虑了６种燃烧产物，８个基元反应，得到了流动参数和产物的质量分数在流场中的分布，为发动机化学动力学损失计算和喷管热结构设计提供了必要数据。计算结果表明，采用半隐格式解决组分连续方程刚性引起的数值不稳定是有效的。     

一种模拟喷管二维流场的新方法

采用ＳＩＭＰＬＥ方法对固体火箭发动机喷管进行了纯气相流场的数值模拟。推导了建立在同位网格之上利用协变物理速度分量为计算变量的离散方程，采用协变物理速度分量推导压力修正方程，对密度采用一阶迎风格式插值，将ＳＩＭＰＬＥ方法扩大到计算可压流动。计算结果表明：本方法计算精度高，收敛速度快，程序编写简单。     

Magnetosphere Imaging Instrument (MIMI) on the Cassini Mission to Saturn/Titan

The magnetospheric imaging instrument (MIMI) is a neutral and charged particle detection system on the Cassini orbiter spacecraft designed to perform both global imaging and in-situ measurements to study the overall configuration and dynamics of Saturn’s magnetosphere and its interactions with the solar wind, Saturn’s atmosphere, Titan, and the icy satellites. The processes responsible for Saturn’s aurora will be investigated; a search will be performed for substorms at Saturn; and the origins of magnetospheric hot plasmas will be determined. Further, the Jovian magnetosphere and Io torus will be imaged during Jupiter flyby. The investigative approach is twofold. (1) Perform remote sensing of the magnetospheric energetic (E > 7 keV) ion plasmas by detecting and imaging charge-exchange neutrals, created when magnetospheric ions capture electrons from ambient neutral gas. Such escaping neutrals were detected by the Voyager l spacecraft outside Saturn’s magnetosphere and can be used like photons to form images of the emitting regions, as has been demonstrated at Earth. (2) Determine through in-situ measurements the 3-D particle distribution functions including ion composition and charge states (E > 3 keV/e). The combination of in-situ measurements with global images, together with analysis and interpretation techniques that include direct “forward modeling’’ and deconvolution by tomography, is expected to yield a global assessment of magnetospheric structure and dynamics, including (a) magnetospheric ring currents and hot plasma populations, (b) magnetic field distortions, (c) electric field configuration, (d) particle injection boundaries associated with magnetic storms and substorms, and (e) the connection of the magnetosphere to ionospheric altitudes. Titan and its torus will stand out in energetic neutral images throughout the Cassini orbit, and thus serve as a continuous remote probe of ion flux variations near 20R _S (e.g., magnetopause crossings and substorm plasma injections). The Titan exosphere and its cometary interaction with magnetospheric plasmas will be imaged in detail on each flyby. The three principal sensors of MIMI consists of an ion and neutral camera (INCA), a charge–energy–mass-spectrometer (CHEMS) essentially identical to our instrument flown on the ISTP/Geotail spacecraft, and the low energy magnetospheric measurements system (LEMMS), an advanced design of one of our sensors flown on the Galileo spacecraft. The INCA head is a large geometry factor (G ∼ 2.4 cm² sr) foil time-of-flight (TOF) camera that separately registers the incident direction of either energetic neutral atoms (ENA) or ion species (≥5^∘ full width half maximum) over the range 7 keV/nuc < E < 3 MeV/nuc. CHEMS uses electrostatic deflection, TOF, and energy measurement to determine ion energy, charge state, mass, and 3-D anisotropy in the range 3 ≤ E ≤ 220 keV/e with good (∼0.05 cm² sr) sensitivity. LEMMS is a two-ended telescope that measures ions in the range 0.03 ≤ E ≤ 18 MeV and electrons 0.015 ≤ E≤ 0.884 MeV in the forward direction (G ∼ 0.02 cm² sr), while high energy electrons (0.1–5 MeV) and ions (1.6–160 MeV) are measured from the back direction (G ∼ 0.4 cm² sr). The latter are relevant to inner magnetosphere studies of diffusion processes and satellite microsignatures as well as cosmic ray albedo neutron decay (CRAND). Our analyses of Voyager energetic neutral particle and Lyman-α measurements show that INCA will provide statistically significant global magnetospheric images from a distance of ∼60 R _S every 2–3 h (every ∼10 min from ∼20 R _S). Moreover, during Titan flybys, INCA will provide images of the interaction of the Titan exosphere with the Saturn magnetosphere every 1.5 min. Time resolution for charged particle measurements can be < 0.1 s, which is more than adequate for microsignature studies. Data obtained during Venus-2 flyby and Earth swingby in June and August 1999, respectively, and Jupiter flyby in December 2000 to January 2001 show that the instrument is performing well, has made important and heretofore unobtainable measurements in interplanetary space at Jupiter, and will likely obtain high-quality data throughout each orbit of the Cassini mission at Saturn. Sample data from each of the three sensors during the August 18 Earth swingby are shown, including the first ENA image of part of the ring current obtained by an instrument specifically designed for this purpose. Similarily, measurements in cis-Jovian space include the first detailed charge state determination of Iogenic ions and several ENA images of that planet’s magnetosphere.This revised version was published online in July 2005 with a corrected cover date.     

UKF方法及其在方位跟踪问题中的应用

采用UKF(Unscented Kalman Filter)方法处理了平面内地面站对目标的方位跟踪的估计问题。目标的位置和速度由选定的高斯分布采样点来近似，在每个更新过程中，采样点随着状态方程传播并随着非线性测量方程变换，由此不但得到目标位置和速度的均值及较高的计算精度，而且避免了对非线性方程的线性化过程。仿真结果表明，UKF方法比传统的扩展卡尔曼滤波(EKF)算法有更高的估计精度，并能有效地克服非线性严重时，方位跟踪问题中很容易出现的滤波发散问题。     

求解可压缩N－S方程的无粘－粘性分数步法

用按物理过程作时间分裂的无粘－粘性分数步法求解可压缩Ｎ－Ｓ方程，无粘步用高效欧拉算法解欧拉方程，粘性步用全隐格式解抛物方程；并采取了有效措施消除由分步法产生的时间和空间不相容性，文中给出了轴对称喷管内外流场的数值模拟结果。     

大瑞利数下二维Benard对流演化过程的DSMC仿真研究

利用DSMC（Direct Simulation Monte Carlo）方法研究了大Rayleigh（瑞利）数条件下稀薄气体的二维Benard对流的特性。考察了流动的非定常演化过程及若干流动特征量随时间的变化特征，得到了与实验较为一致的结论。计算结果表明，当流动控制参数不断增大时，二维Berard对流从热传导状态过渡到热对流状态表现出明显的非定常演化特征。最后当瑞利数增大到一定程度后，注动显示出混沌的特性。     

真实气体流动的相似规律

本文从具有化学反应的NS方程出发导出了真实气体流动的相似律,分别给出了高温空气非平衡流、平衡流和冻结流的相似参数。文中还讨论了航天飞机轨道器各种Mach数范围内的真实气体流动的相似参数和实验模拟问题。     

轴系稳定性裕度判据的理论研究

本文提出以系统最大能量比率作为衡量轴系稳定性裕度的判据,和传统的对数衰减率与系统阻尼判据相比,该判据综合了系统所有模态的信息,考虑了特征值与特征向量的共同影响。最大能量比率较小的轴系具有较大的稳定性裕度;最大能量比率较大的轴系具有较小的稳定性裕度。探讨了特征向量在评估稳定性裕度的重要作用,指出了传统的对数衰减率与系统阻尼判据在评估轴承稳定性裕度时的局限性。作为应用实例,用最大能量比率判据考察了一模型转了轴承系统的稳定性裕度,并与用对数衰减率及系统阻尼判据所得结论进行了对比。结果表明,用本文提出的判据所得的结论与实际情况比较一致。此外,应用模型转子的计算结果,进一步阐明了对数衰减率与系统阻尼判据的共同不足是采集的信息太少。     

航天飞行中的稀薄气体动力学

飞船、宇航探测器、航天飞机等复杂外形航天器给气体动力学，包括稀薄气体动力学提出了新的要求。本文简要介绍了为计算过渡领域中气动力与热而发展的基于位置元概念的DSMC方法的通用算法。该方法解决了计算物面通量量的技术难点并已用于模拟圆球、飞船、类航天飞机的绕流。正在进行的航天实践，如麦哲伦飞船对金星的探测、行星大气中的气动制动、伽利略飞船的木星之行、尾屏蔽在太空中获得高真空的实验等等提出了新的气动力问题，稀薄气体动力学和DSMC方法是有力的工具。     

燃烧室头部进气规律对煤粉燃烧效率的影响

用实验方法研究了火焰筒式煤粉燃烧室头部进气规律对燃烧效率的影响。通过对6种结构的头部配气方案的分析比较,得出了煤粉燃烧室合理的头部配气方案,为今后煤粉燃烧室的设计和工作状态的选择提供了重要的实验依据。