微波等离子推力器喷管流动数值模拟

为了准确把握微波等离子推力器喷管流动的机理与特点,采用冻结流和非平衡流两种模型对其进行了对比数值模拟分析。非平衡流模型考虑了流动过程中的分解反应、电离反应和复合反应,化学动力学模型为4组分、4反应的有限速率化学反应模型,采用二阶精度NND格式数值求解耦合化学反应源项的N S方程组。数值模拟的结果揭示了喷管内的流场结构,反映了喷管内的离解电离状况,得到了推力器的推力和比冲。分析表明,数值方法有效,计算结果合理,具有工程应用价值,能为相关研究提供参考。     

月球探测软着陆与采样返回段弹道确定

针对月球探测中软着陆与采样返回段弹道计算问题,提出用数值逼近弹道确定方法。通过B样条对探测器状态进行建模,进而综合全弧段数据进行统计定轨的方法。由于样条法良好的数值逼近性能,使得该方法对探测器弹道异常复杂情况下的状态确定较为有效。对嫦娥三号探测器动力软着陆弧段进行了仿真与实测数据处理。分析了采样返回段的基本动力学与控制特征,为后续的嫦娥五号探测器的软着陆及其采样返回提供初步的可行弹道计算方法。在嫦娥三号探测器动力落月段实测数据处理中,通过评估,该段弹道确定精度优于100 m,其弹道末点与NASA的月球勘测轨道器（LRO）给出的结果差异优于50 m,证实了文章提出的软着陆弹道确定方法的有效性。     

直升机行星架疲劳裂纹扩展寿命预测

直升机行星齿轮保持架（简称行星架）是传动系统的重要部件,其可靠性对于直升机的飞行安全至关重要。行星架疲劳裂纹故障的发生、发展受多种因素影响,其故障诊断和寿命预测都有难度。为准确预测行星架疲劳裂纹寿命,研究了其裂纹故障发生、发展规律,提出了基于对数线性的方法,将裂纹扩展过程离散化处理,采用Paris公式,定量描述裂纹扩展速率,结合Miner准则,累积其疲劳损伤过程,最终得到随裂纹长度变化的行星架剩余使用寿命数值。利用有限元软件ANSYS Workbench的疲劳寿命模块对计算得到的数值结果进行了仿真对比,验证了方法的精确性和有效性。      

一种基于数字聚束的极坐标格式算法波前弯曲补偿

提出了一种基于数字聚束技术的PFA波前弯曲误差补偿新方法,该方法首先对回波数据进行两维空域滤波和降采样,将原始的宽波束大场景数据分解成若干窄的子波束子场景数据,再结合PFA的非共面运动补偿能力,能够有效解决在雷达任意航迹条件下的PFA波前弯曲误差补偿问题,显著提高PFA有效成像场景范围。仿真数据处理表明了算法的有效性。     

非线性广义时变时滞系统鲁棒模糊控制

针对一类参数不确定模糊广义时变时滞系统,研究该系统的鲁棒稳定性控制问题。选取特殊的Lyapunov函数,给出了该系统稳定的充分条件,采用线性矩阵不等式技术,将时变时滞系统稳定性条件转化为求解一组线性矩阵不等式问题,并设计出状态反馈控制器,仿真结果说明该方法的可行性。     

随机事元及其在解决矛盾问题中的应用

通过建立随机事件的可拓模型,给出了随机事元、随机事元的概率、随机事元集和可拓事件集的概念。在对随机事元可拓性研究的基础上给出了随机事元的拓展概率的概念。利用随机事元的多特征性,讨论了当随机事元的某一个或若干个特征及其量值改变时,由可拓变换的传导性,将在随机事元的各要素之间产生传导效应,从而导致随机事元概率的改变,进而为涉及随机事件发生的概率的矛盾问题的解决提供了有效途径。最后举例说明了利用概率的改变处理矛盾问题的方法。     

微波等离子体推力器流动模拟(英文)

采用N-S方程求解了100 W微波等离子体推力器(MPT)选用不同推进工质时的性能参数;并采用直接蒙特卡洛模拟方法(DsMC)对MPT羽流进行了数值模拟.结果表明,几种工质的推力变化不大,氮气为23.6 mN,氮气为24.8mN,氩气为24.8 nuN;但比冲区别较大,氮气为565.2 s,氮气为243.7 8,氢气为180.2 s.羽流场中,密度、压强及温度沿轴向和径向均逐渐减小;轴向速度在轴线附近变化不大,采用氩气工质时,约1 700 m/s,在远离轴线区域,沿流动方向逐渐增大,沿径向逐渐减小;径向速度沿轴向变化不大,沿径向逐渐增大,并在接近流动区域边界时迅速减小.     

膏体冲压发动机构型对燃烧效率的影响研究

膏体冲压发动机是在膏体发动机和冲压发动机基础上提出的一种新型组合动力装置.探讨了膏体冲压发动机的多种构型.并对其反应流场进行了数值模拟和燃烧效率分析.结果表明,采用将膏体富燃料直接喷射到补燃室的方式燃烧效率较高;随着进气道轴向距离的增加,补燃室头部的燃烧效率有所降低;随着进气道夹角的增加,补燃室燃烧效率略有增加.     

狭小空间内气膜孔流量系数的数值模拟

针对多层层板冷却、双层壁冷却等冷却技术,利用数值模拟分析了这些冷却结构形成的狭小冷却通道中,气膜孔附近的流动特性,重点研究了气膜孔流量系数C_{d随吹风比M(0.5～2.0)、气膜孔雷诺数Re(5000～10000)以及冷却通道高度H和气膜孔直径D之比H/D(0.33～1.0)等参数的变化规律.数值计算中湍流模型为Realizable k-ε模型,近壁处采用非平衡壁面函数,利用SIMPLE算法和二阶迎风格式进行离散求解.计算结果表明:气膜孔流量系数随吹风比的增大而增大,在吹风比M小于1时,影响尤为明显.研究中同时发现在相同吹风比的条件下,Cd随着气膜孔雷诺数的增大而减小,但变化的幅度不大;在相同的气膜孔雷诺数下,Cd随着H/D的减小而降低,特别是当H/D小于0.75时,随着流体在进入气膜孔前逐步受到限制,Cd随着H/D的减小而快速降低.}      

A Simulation of the Mid-and Low-latitude Ionospheric Electric Fields

A theoretical model of ionospheric electric fields at mid- and low-latitudes is developed. In the geomagnetic dipolar coordinate system, the ionospheric dynamo equations were solved, and the ionospheric electric potential and electric field were derived respectively. Major parameters for the model inputs, such as the neutral winds, the densities and temperatures of electron, ions and neutrals, are obtained from empirical models. The global ionospheric electrical potential and field at mid- and low-latitudes derived from our model are largely in agreement with the results presented by other authors and the empirical model. Using our model, it is found that the diurnal component of the HWM93 wind mainly contributed to the formation of the vertical electric field, while the semidiurnal component mainly contributed to the zonal electric field. Finally, by adjustment of the input F region winds and conductivities, most discrepancies between our model and the empirical one can be eliminated, and it is proved that the F region dynamo is the most significant contribution to the electric fields.