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

基于直接蒙特卡洛(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%则会发生熄火。     

航天器运输包装箱仿真验证平台技术研究

针对数量有限的物理跑车试验无法满足减振与保温性能测试需求的问题，提出一套航天器运输包装箱动力学与热学仿真验证方法，包括：建立适用于包装箱系统的刚柔耦合多体动力学系统，通过结合线路条件测试生成的动力学系统外部激励，实现减振性能虚拟跑车测试；建立基于计算流体力学的包装箱热学模型，通过模拟自然对流和空调控制，实现包装箱保温性能虚拟跑车测试；基于C/S架构和导航式流程设计思想，建立航天器运输包装箱仿真验证平台，通过实际案例证明该平台仿真结果与实际跑车测试数据具有较高的一致性。     

空天飞行器返回滑翔段在线制导方法

针对空天往返飞行器的返回滑翔段在线制导问题，设计了一种新的滑翔段飞行剖面，实现了滑翔段终端交班高度、位置和倾角约束的自动满足，减少了在线制导算法中需处理的约束数量。推导了滑翔段运动状态、过程约束和性能指标的解析表达式，获得了剩余航程和终端速度间的函数关系。在此基础上，提出了一种双层在线制导方法：内层解析重构飞行剖面，同时通过解析确定路径点来改变剩余航程的变化率，进而对终端交班速度进行控制；外层借助解析表达式，使用粒子群优化算法(PSO)和改进共轭梯度法(CGM)优化飞行剖面，从而满足过程约束和指标要求。最后通过数学仿真验证了方法的正确性。     

多臂空间机器人的视觉伺服与协调控制

针对多臂空间机器人自主目标抓捕任务，首先建立多臂空间机器人的运动模型和其与目标的相对运动模型，采用Kane方法建立多臂空间机器人的动力学模型；其次，研究基于视觉伺服的机械臂在线轨迹规划算法，并引入零反作用机动，消除机械臂运动对平台姿态的扰动；再次，在不使用零反作用机动功能时，分别使用基于角动量前馈补偿的协调控制算法和逆动力学方法设计了协调控制器，在机械臂运动时保持平台姿态和相对目标的位置。最后，开发了基于Matlab的仿真软件MASS（多臂空间机器人仿真），仿真结果校验了上述方法的有效性。     

Evaluation of SARAL/AltiKa performance using GNSS/IMU equipped buoy in Sajafi,Imam Hassan and Kangan Ports

Performance of SARAL/AltiKa mission has been evaluated within 2016 altimeter calibration/validation framework in Persian Gulf through three campaigns conducted in the offshore waters of Sajafi, Imam Hassan and Kangan Ports, while the altimeter overflew the passes 470, 111 and 25 on 13 Feb, 7 March and 17 June 2016, respectively. As the preparation, a lightweight buoy was equipped with a GNSS receiver/choke-ring antenna and a MEMS-based IMU to measure independent datasets in the field operations. To obtain accurate sea surface height (SSH) time series, the offset of the onboard antenna from the equilibrium sea level was predetermined through surveying operations as the buoy was deploying in the onshore waters of Kangan Port. Accordingly, the double-difference carrier phase observations have been processed via the Bernese GPS Software v. 5.0 so as to provide the GNSS-derived time series at the comparison points of the calibration campaigns, once the disturbing effects due to the platform tilt and heave have been eliminated. Owing to comparing of the SSH time series and the associating altimetry 1?Hz GDR-T datasets, the calibration/validation of the SARAL/AltiKa has been performed in the both cases of radiometer and ECMWF wet troposphere corrections so as to identify potential land contamination. An agreement of the present findings in comparison with those attained in other international calibrations sites confirms the promising feasibility of Persian Gulf as a new dedicated site for calibration/validation of ongoing and future altimetry missions.     

Assessment of between-floor structural and topological properties on cognitive map development in multilevel built environments

The present study investigated cognitive map development in multilevel built environments. Three experiments were conducted in complex virtual buildings to examine the effects of five between-floor structural factors that may impede the accuracy of humans’ ability to build multilevel cognitive maps. Results from Experiments 1 and 2 (of three experiments) revealed that difficulties in developing multilevel cognitive maps are not solely caused by the z-axis offset, as is suggested in the literature, but are due to the factorial combination of a between-floor overlap and a z-axis offset. Results from Experiment 2 showed that this process becomes substantially more difficult when the reference directions between different floors have an angular offset from each other. Finally, results from Experiment 3 demonstrated that confusing between-floor heading shifts in aligned buildings did not make it reliably harder to build multilevel cognitive maps. The implications of these findings are discussed in terms of theories of mental representations in multilayered three-dimensional spaces, as well as for architectural design.     

Simulating planetary wave propagation to the upper atmosphere during stratospheric warming events at different mountain wave scenarios

Parameterization schemes of atmospheric normal modes (NMs) and orographic gravity waves (OGWs) have been implemented into the mechanistic Middle and Upper Atmosphere Model (MUAM) simulating atmospheric general circulation. Based on the 12-members ensemble of runs with the MUAM, a composite of the stratospheric warming (SW) has been constructed using the UK Met Office data as the lower boundary conditions. The simulation results show that OGW amplitudes increase at altitudes above 30 km in the Northern Hemisphere after the SW event. At altitudes of about 50 km, OGWs have largest amplitudes over North American and European mountain systems before and during the composite SW, and over Himalayas after the SW. Simulations demonstrate substantial (up to 50–70%) variations of amplitudes of stationary planetary waves (PWs) during and after the SW in the mesosphere-lower thermosphere of the Northern Hemisphere. Westward travelling NMs have amplitude maxima not only in the Northern, but also in the Southern Hemisphere, where these modes have waveguides in the middle and upper atmosphere. Simulated variations of PW and NM amplitudes correspond to changes in the mean zonal wind, EP-fluxes and wave refractive index at different phases of the composite SW events. Inclusion of the parameterization of OGW effects leads to decreases in amplitudes (up to 15%) of almost all SPWs before and after the SW event and their increase (up to 40–60%) after the SW in the stratosphere and mesosphere at middle and high northern latitudes. It is suggested that observed changes in NM amplitudes in the Southern Hemisphere during SW could be caused by divergence of increased southward EP-flux. This EP-flux increases due to OGW drag before SW and extends into the Southern Hemisphere.     

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.