脉冲爆震发动机管外复杂波系的数值计算

为了研究脉冲爆震发动机外流场复杂波系的变化特点,用时空守恒元与求解元(CE/SE)方法对单爆震管和三爆震管脉冲爆震发动机外流场进行了计算。三个爆震管并排排列,管内填充按化学当量比混合的乙烯氧气混和物。通过计算,获得了脉冲爆震发动机外流场的压力分布、温度分布、激波与膨胀波的变化特点。计算结果表明该时空守恒元与求解元方法是一种较好的爆震波模拟方法,能有效地捕捉激波等强间断;三爆震管脉冲爆震发动机相对于单爆震管脉冲爆震发动机外流场复杂,存在多道相互作用、相互影响的激波和膨胀波。     

SE(3)上航天器姿轨耦合固定时间容错控制

针对相对运动航天器姿轨一体化控制问题,考虑执行器故障和控制输入饱和的影响,提出了一种基于滑模的模糊自适应固定时间容错控制方法。首先,在李群SE（3）的框架下建立并推导相对运动航天器姿轨一体化误差动力学模型;其次,引入执行器故障和控制输入饱和的问题,采用双幂次快速终端滑模面,并结合模糊自适应方法设计了固定时间稳定的容错控制器,可以实现执行器故障情况下相对运动航天器的高精度快速跟踪控制;然后,运用Lyapunov方法证明了系统的稳定性,该控制器不仅能不依赖于系统的初始状态实现滑模趋近和到达阶段的固定时间稳定性,而且由于采用模糊逼近方法结合自适应更新策略可以实时高精度地估计系统的总扰动信息,因此可以达到快速高精度的容错控制目标;最后,对所提出的的控制方法进行数值仿真分析,结果验证了该方法的有效性和可行性。     

基于CE/SE法的翼型流场激波捕捉计算

在二维四边形结构网格下,针对时空守恒元和解元(CE/SE, space-time Conservation Element and Solution Element)方法捕捉激波需要双倍空间网格点的问题,重新对守恒元和解元进行了定义,计算点同时包含单元中心和网格节点,并推导得到了新的CE/SE方法的计算公式.以此为基础,结合非敏感克朗数(CNIS, Courant Number Insensitive Scheme)计算格式和当地时间步长方法,对激波翼型流场进行数值模拟,并与原方法及AGARD报告进行对比.结果表明,对原CE/SE方法的改进是有效的,可明显提高激波分辨能力,并且激波前后无明显的数值震荡发生,新方法适合应用于翼型流场中捕捉激波.     

水下爆轰燃气泡形态与激波传播过程研究*(爆震专刊)

针对脉冲爆轰发动机在水下工作过程中形成的燃气射流问题,搭建了水下爆轰燃气实验系统,研究了第一个爆轰循环在水下的燃气泡发展变化过程。建立了基于气液两相双流体模型的脉冲爆轰发动机水下喷射模型,采用时-空守恒元和求解元方法,模拟了爆轰波退化为激波在水中的传播及衰减过程。研究结果表明：燃气泡前期受外界水环境阻滞作用呈现“豌豆状”形态,充分发展阶段气液交界面逐渐失稳,在达到最大尺寸后开始收缩并在中心轴线位置出现凸出的射流；水下爆轰燃气射流发展过程中同时存在脱离燃气泡的水中前导激波和管口燃气泡内的高压区两部分,水中前导激波在传播过程中压力迅速衰减至常压量级,而管口燃气泡内则一直保持较高压力；中心轴线区域气液交界面处反射激波的回传使管口附近出现回击现象,并导致前导激波波阵面上压力峰值逐渐出现在30°方向上。     

地球边缘大气密度时空特性分析

分析了地球边缘大气密度的时变特性及空间分布特性。从物理原因出发并基于NRLMSISE-00和MET-99大气模型计算随不同时空变量变化的相对大气密度,以说明各时空因素对地球边缘大气密度的影响程度。分析结果表明地球边缘大气密度具有较强的时变特性,昼夜和季节均对大气密度产生一定影响。空间特性表现为纬度对大气密度的影响明显,而经度影响微弱。季节不同导致纬度的影响程度不同,纬度因素在夏季和冬季的影响强于春秋。随着纬度的增加,季节对地球边缘大气密度的影响逐渐增强,而当地时间的影响呈减弱趋势。相比较而言,低纬区的昼夜波动强于季节波动,高纬区的季节波动强于昼夜波动。     

Motion of ejecta around a binary system with an ellipsoid and a sphere

A kinetic impact occurs when an asteroid is moving. The impact may be caused by other small celestial body or an artificial object. Ejecta is produced from an impact. Dynamics of ejecta near a binary system which contains an ellipsoid and a sphere is analyzed. A phase diagram which comes from numerical simulation is shown in this work. The phase diagram shows a clear structure. Some special trajectories are also shown, which indicates a potential source of asteroid orbiting objects.     

Multi-objective optimal formation reconfiguration with scalarization and stochastic boundaries

A new method is proposed for multi-objective optimal control of satellite formation reconfiguration. First, necessary optimality conditions of single-objective fuel-optimal reconfiguration are studied. Problem of initial guess of the costate and control discontinuity are addressed by using a new developed stochastic based method. Terminal conditions of Two Point Boundary Value Problem (TPBVP) of the optimal control equations are relaxed using augmented Gaussian. Variances or bandwidths of the terminal conditions are set to decision variables to be optimized. Responses of varied initial costate are modeled and propagated using Riccati equation. Using quadratic convolution of relaxed terminal conditions, objective function of the fuel-optimal problem is reformulated into a quadratic equation, and a gradient based optimizer is used to search the optimal initial guess. In the second part of the paper, optimal control of format flying reconfiguration with two competitive objectives:fuel and time of flight (ToF) is studied. With Pascoletti-Serafini scalarization, the Multi-objective Optimal Control Problem (MOCP) is scalarized into a set of weighted Single-objective Optimal Control Problems (SOCP). Hamiltonian and switch function for the weighted SOCP are developed and constructed. New developed stochastic based approach is then used to search the optimal Pareto solutions. Verification and performances of proposed algorithms are demonstrated through numerical simulations of single-objective and multi-objective optimal control of low-thrust elliptical orbit formation reconfiguration. Potential applications of the algorithms to some other preliminary space mission designs are also discussed.     

铝粉浓度对汽油-纳米铝粉悬浮液滴多相爆轰影响的一维数值计算

采用带化学反应的一维多相流模型,运用CE/SE数值计算方法和处理源项的四阶龙格-库塔方法对爆轰管内汽油/纳米铝粉悬浮液滴、空气混合物的多相爆轰过程进行了数值模拟,讨论了悬浮液滴中纳米铝粉浓度对爆轰参数及燃烧转爆轰过程的影响。数值结果表明,添加纳米铝粉后,稳定传播的爆轰波速度变大;随着纳米铝粉浓度的增加,爆轰波压力峰值与...     

Protel99在模拟电子线路中的仿真方法

文章介绍了当今比较先进的电子线路EDA辅助设计软件Protel99－电路仿真部分，阐明了电路仿真原理，特点和仿真应用技术，并给出一个实例。     

Sq(H) field variations at Sonmiani geomagnetic observatory,Pakistan for solar cycle 24

In this study, Sq(H) field results for Sonmiani geomagnetic observatory (SON), Pakistan are presented first time for solar cycle i.e., Solar Cycle 24. Sonmiani observatory was established in 2008 and was included in the list of INTERMAGNET Magnetic Observatory (IMO) in 2012. The analysis of Sq(H) field is performed to examine diurnal, annual and seasonal behaviour of Sq(H) field at Sonmiani along with Honolulu which is used as reference observatory. Both the observatories lie at the Equatorial Ionisation Anomaly (EIA) crest. In general strong dependency of Sq(H) field on solar cycle has been observed. First peak of SC24 was visible in Sq(H) field at both observatories. However, at HON maximum Sq(H) was observed in 2015 instead of 2014, a year later than the year of solar maximum. Prominent longitudinal difference has been noted between both observatories. At HON, maximum Sq(H) was observed during equinox while at SON it was noted during equinox and summer as a consequence of shift in latitude of Sq focus. Phase shift of Sq(H) at SON followed a general trend, that is maximum Sq(H) shifted to later hours in solar maximum as compared to solar minimum. Whereas, an opposite trend was noted at HON which might be due to coastal effect. In case of seasonal phase shift, at both observatories maximum Sq(H) occurred at later hours in summer than in winter which is opposite to the results obtained by many workers. These points need further investigation and geomagnetic as well as wind, and electric field data of other observatories situated at the EIA crest region is required to interpret the Sq(H) phase shift extensively.