碳纳米管阵列推力器放电特性的数值模拟研究

基于直接蒙特卡洛(DSMC)算法、浸入式有限元算法(IFE)、粒子云及蒙特卡洛碰撞(PIC-MCC)算法以及Ammosov-Delone-Krainov(ADK)隧穿电离模型建立了碳纳米管阵列推力器工质气体和CNT场电离过程的三维仿真模型。利用文献中的实验数据对本文仿真模型的有效性进行了验证,结果表明本文仿真模型与实验结果基本一致。针对碳纳米管阵列推力器的放电特性进行了研究,并重点分析了离子与原子碰撞过程对碳纳米管阵列推力器放电特性的影响。结果表明CNT可以极大地增强局部电场,稳态时场增强因子约为1308;与原子的碰撞使得离子在CNT附近更加集中,降低了局部电场,进而降低了其场增强能力。     

Solar sailing in realistic mode,based on a locally-optimal control laws

Ballistic design of solar sailing missions in the solar system is composed of defining the design parameters, the control programs, and the trajectories that provide performance goals of a flight. The use of a solar sail spacecraft imposes specific restrictions on mission parameters that include the degradation limit on the flight duration, the maximum temperature of solar sail's surface, the minimum distance from the Sun, the maximum angular velocity of the spacecraft's rotation and others.Many authors considered the impact of these restrictions on the design of the mission separately, but they used a sophisticated method of finding the exact optimal motion control or applied the most straightforward laws of motion control. This paper uses local-optimal control laws at the complete mathematical models of motion and functioning of solar sail spacecraft to describe a technique of designing interplanetary missions. The described method avoids the need to obtain an accurate optimal solution to the control problem and does not cause significant computational difficulties.     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

天基轻气炮发射清除低轨碎片方案及关键技术分析

文章分析了现有的空间碎片清除方式,并以800~1200 km低地球轨道高度上1~10 cm量级的空间碎片为清除目标,提出了天基轻气炮清除碎片的新方法。首先分析了轻气炮有效载荷在典型参数下的弹丸加速能力;之后根据将碎片降轨使其坠入大气层烧毁的设想,提出天基轻气炮共面清除碎片的方式,并选择轨道高度800 km的圆轨道作为碎片运行轨道进行可行性分析。计算表明,对半径10 cm、厚度1 cm的铝合金圆板碎片(质量211.95 g),使用初速1 km/s、重10 g的黏性弹丸可按任务方案达到清除效果。此外,计算出该参数弹丸对轨道高度800~1200 km的圆轨道上可清除的最大碎片质量为500~825 g,证明轻气炮弹丸对1~10 cm的碎片具有较强的清除能力。最后,分析了以轻气炮为有效载荷的航天器在完成清除碎片任务时的关键技术。     

飞机液压系统柱塞泵热分析

随着飞机液压系统向高压化、大功率化方向发展，液压系统的发热及温升问题受到了广泛关注。柱塞泵作为液压系统中重要的能量转换装置，对液压系统的功率损失和温度有重大的影响，良好的柱塞泵的热分析对液压系统的热设计和温度控制有重要的意义。然而传统的柱塞泵热分析方法，如平均油温法、软件仿真法和热节点网络法，因不够准确、高效和简洁而不能满足现代飞机液压系统热分析和热设计的相关需求。本文针对飞机液压系统柱塞泵发热的相关问题，采用基于集总参数的热网络分析方法，建立了柱塞泵热网络模型，通过对柱塞泵不同工况下的仿真，验证了模型的有效性和准确性，为柱塞泵及液压系统的热分析提供了参考。     

Comparison of vertical E × B drift velocities and ground-based magnetometer observations of DELTA H in the low latitude under geomagnetically disturbed conditions

In the present work, we analyzed the daytime vertical E × B drift velocities obtained from Jicamarca Unattended Long-term Ionosphere Atmosphere (JULIA) radar and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca, and Piura in Peru for 22 geomagnetically disturbed events in which either SC has occurred or Dst_max < ?50 nT during the period 2006–2011. The ΔH component of geomagnetic field is measured as the differences in the magnitudes of horizontal H component between magnetometer placed directly on the magnetic equator and one displaced 6–9° away. It will provide a direct measure of the daytime electrojet current, due to the eastward electric field. This will in turn gives the magnitude of vertical E × B drift velocity in the F region. A positive correlation exists between peak values of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of the events. It was observed that 45% of the events have daytime vertical E × B drift velocity peak in the magnitude range 10–20 m/s and 20–30 m/s and 20% have peak ΔH in the magnitude range 50–60 nT and 80–90 nT. It was observed that the time of occurrence of the peak value of both the vertical E × B drift velocity and the ΔH have a maximum (40%) probability in the same time range 11:00–13:00 LT. We also investigated the correlation between E × B drift velocity and D_st index and the correlation between delta H and D_st index. A strong positive correlation is found between E × B drift and D_st index as well as between delta H and D_st Index. Three different techniques of data analysis – linear, polynomial (order 2), and polynomial (order 3) regression analysis were considered. The regression parameters in all the three cases were calculated using the Least Square Method (LSM), using the daytime vertical E × B drift velocity and ΔH. A formula was developed which indicates the relationship between daytime vertical E × B drift velocity and ΔH, for the disturbed periods. The E × B drift velocity was then evaluated using the formulae thus found for the three regression analysis and validated for the ‘disturbed periods’ of 3 selected events. The E × B drift velocities estimated by the three regression analysis have a fairly good agreement with JULIA radar observed values under different seasons and solar activity conditions. Root Mean Square (RMS) errors calculated for each case suggest that polynomial (order 3) regression analysis provides a better agreement with the observations from among the three.     

Geometric error analysis of an over-constrained parallel tracking mechanism using the screw theory

This paper deals with geometric error modeling and sensitivity analysis of an overconstrained parallel tracking mechanism. The main contribution is the consideration of overconstrained features that are usually ignored in previous research. The reciprocal property between a motion and a force is applied to tackle this problem in the framework of the screw theory. First of all, a nominal kinematic model of the parallel tracking mechanism is formulated. On this basis, the actual twist of the moving platform is computed through the superposition of the joint twist and geometric errors. The actuation and constrained wrenches of each limb are applied to exclude the joint displacement. After eliminating repeated errors brought by the multiplication of wrenches, a geometric error model of the parallel tracking mechanism is built. Furthermore,two sensitivity indices are defined to select essential geometric errors for future kinematic calibration. Finally, the geometric error model with minimum geometric errors is verified by simulation with SolidWorks software. Two typical poses of the parallel tracking mechanism are selected, and the differences between simulation and calculation results are very small. The results confirm the correctness and accuracy of the geometric error modeling method for over-constrained parallel mechanisms.     

飞秒激光系统图像调焦装置以及方法

飞秒激光直写技术在复杂三维微结构加工领域具有显著优势，而调焦是否精准直接影响了所加工结构的完整度.提出了在光路中临时置入调焦光源和物的图像调焦技术，通过调节物的位置使其成像面与激光聚焦面一致，从而通过清晰可分辨的成像状态间接反映激光聚焦状态.利用Zemax软件模拟分析了原飞秒激光光路与加入调焦光源和物的调焦光路，二者可实现相同加工物镜后工作距离与良好成像质量，证明了该方法的可行性.通过分析得到该过程的成像误差主要由成像镜头焦深(3.9 μm)引起，我们获得的理想调焦精度可达到1/2焦深以内.设计了单层高度为5 μm的二层圆柱结构，通过多次实验验证了所加工元件高度误差在1.5 μm范围以内，与理论分析一致，满足飞秒激光系统的调焦要求.     

超临界环境下煤油替代物液滴燃烧特性数值研究

为研究液氧/煤油火箭发动机燃烧室内经喷注形成的煤油液滴的燃烧过程，基于实际气体状态方程、高压热物性修正、高压气液平衡和详细化学反应动力学，建立一维的全瞬态液滴燃烧模型，对超临界环境下两组分煤油替代物液滴的燃烧特性及液滴初始直径的影响进行仿真研究。结果表明，在超临界环境下，相比于煤油液滴纯蒸发过程，煤油液滴燃烧过程的迁移时刻大大提前；煤油液滴着火之后很快进入超临界燃烧阶段，此时液滴燃烧过程可以看成中心附近的燃料高浓度区与外侧氧气高浓度区之间的扩散燃烧过程；煤油液滴的火焰半径先增大，达到最大值之后开始减小，并减小为零，火焰温度在着火之后快速上升至最大值，并基本保持不变，在火焰半径减小为零之后开始降低；随着液滴初始直径的增大，火焰特性以及液滴中心参数变化曲线趋势不变、整体延迟，着火时间、迁移时间和液滴寿命增大。     

火箭发动机超声速过膨胀射流气动噪声特性研究

采用大涡模拟LES方法计算了火箭发动机超声速过膨胀射流形态及近场声压分布,研究了入口温度与环境温度的比值(温度比)对声场的影响;将声源分解,基于Ffowcs Williams-Hawkings (FW-H)方程获取了不同位置噪声源的远场噪声,并根据声压级频谱和湍流形态分析了超声速射流噪声的产生机理。研究表明,超声速过膨胀射流气动噪声由湍流混合噪声和宽频激波噪声组成,近场噪声源以马赫波形式向大方位角辐射中高频噪声,下游大尺度湍流向低方位角范围辐射低频噪声,声压级峰值频率随观测角度增大而升高;随温度比升高,马赫波辐射角度增大,噪声指向性发生改变。该研究可为运载火箭发动机地面试车或火箭发射段声学环境设计提供参考。