基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

A progressive approach to predict shot peening process parameters for forming integral panel of Al7050-T7451

In this paper, a progressive approach to predict the multiple shot peening process parameters for complex integral panel is proposed. Firstly, the invariable parameters in the forming process including shot size, mass flow, peening distance and peening angle are determined according to the empirical and machine type. Then, the optimal value of air pressure for the whole shot peening is selected by the experimental data. Finally, the feeding speed for every shot peening path is predicted by regression equation. The integral panel part with thickness from 2 mm to 5 mm and curvature radius from 3200 mm to 16000 mm is taken as a research object, and four experiments are conducted. In order to design specimens for acquiring the forming data, one experiment is conducted to compare the curvature radius of the plate and stringer-structural specimens, which were peened along the middle of the two stringers. The most striking finding of this experiment is that the outer shape error range is below 3.9%, so the plate specimens can be used in predicting feeding speed of the integral panel. The second experiment is performed and results show that when the coverage reaches the limit of 80%, the minimum feeding speed is 50 mm/s. By this feeding speed, the forming curvature radius of the specimens with different thickness from the third experiment is measured and compared with the research object, and the optimal air pressure is 0.15 MPa. Then, the plate specimens with thickness from 2 mm to 5 mm are peened in the fourth experiment, and the measured curvature radius data are used to calculate the feeding speed of different shot peening path by regressive analysis method. The algorithm is validated by forming a test part and the average deviation is 0.496 mm. It is shown that the approach can realize the forming of the integral panel precisely.     

星载激光通信终端在轨机动对温度影响

对激光通信终端在轨瞬态温度变化开展了仿真，以期研究在轨机动的影响、热分析和热试验时机动模式的简化模拟。通过合理地分析与简化，建立终端的轨道热分析模型，准确模拟在轨机动，根据典型机动模式、外热流和涂层退化等因素设计了计算工况，得出了终端在轨运行过程中的温度场随时间和姿态变化规律。结果表明：不同机动模式下的温度变化存在差异，最大可达23.0℃，采用固定姿态或一维转动的简化热分析或热试验不能准确模拟实际飞行温度，甚至不能部分替代二维转动热分析或热试验；光学天线和反射镜的温度控制是制约终端工作的瓶颈，为终端设计合适的避光机动策略，可大幅度提高温度稳定度和均匀度。研究结果可以为光机电设备在轨机动策略的设计和改进提供参考。     

Topology optimization in lightweight design of a 3D-printed flapping-wing micro aerial vehicle

Topology optimization is an effective method to obtain a lightweight structure that meets the requirements of structural strength. Whether the optimization results meet the actual needs mainly depends on the accuracy of the material properties and the boundary conditions, especially for a tiny Flapping-wing Micro Aerial Vehicle (FMAV) transmission system manufactured by 3D printing. In this paper, experimental and numerical computation efforts were undertaken to gain a reliable topology optimization method for the bottom of the transmission system. First, the constitutive behavior of the ultraviolet (UV) curable resin used in fabrication was evaluated. Second, a numerical computation model describing further verified via experiments. Topology optimization modeling considering nonlinear factors, e.g. contact, friction and collision, was presented, and the optimization results were verified by both dynamic simulation and experiments. Finally, detailed discussions on different load cases and constraints were presented to clarify their effect on the optimization. Our methods and results presented in this paper may shed light on the lightweight design of a FMAV.     

Accurate approximation of in-ecliptic trajectories for E-sail with constant pitch angle

Propellantless continuous-thrust propulsion systems, such as electric solar wind sails, may be successfully used for new space missions, especially those requiring high-energy orbit transfers. When the mass-to-thrust ratio is sufficiently large, the spacecraft trajectory is characterized by long flight times with a number of revolutions around the Sun. The corresponding mission analysis, especially when addressed within an optimal context, requires a significant amount of simulation effort. Analytical trajectories are therefore useful aids in a preliminary phase of mission design, even though exact solution are very difficult to obtain. The aim of this paper is to present an accurate, analytical, approximation of the spacecraft trajectory generated by an electric solar wind sail with a constant pitch angle, using the latest mathematical model of the thrust vector. Assuming a heliocentric circular parking orbit and a two-dimensional scenario, the simulation results show that the proposed equations are able to accurately describe the actual spacecraft trajectory for a long time interval when the propulsive acceleration magnitude is sufficiently small.     

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

某型装备能源电池异常爆破故障分析

对一起因电池设计缺陷引发装备战训失利的故障进行分析,重点对电池工作机理和设计缺陷展开解析,明确故障定位,并对国内外装备弹上不同体系电池的优缺点进行比较,旨在为产品维护保障、国产化科研等工作提供参考。     

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

文章分析了现有的空间碎片清除方式,并以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的碎片具有较强的清除能力。最后,分析了以轻气炮为有效载荷的航天器在完成清除碎片任务时的关键技术。     

多种工质下螺旋波等离子体源波场结构数值模拟

针对螺旋波等离子体放电机理，开展了多种工质条件下的螺旋波放电等离子体内波场结构数值模拟研究。计算发现:氦气等离子体的Er分量在径向边界处的峰值更为突出，有利于等离子体在径向的输运，波电场径向分量决定了电流密度径向分量在内部的表现。在0.266 Pa和1.064 Pa两种气体压强条件下，通过波场结构验证了气压对于波阻尼影响的结论。波场结构是螺旋波在等离子体内传播以及能量沉积的微观体现，研究螺旋波波场结构是揭示其高电离效率的重要途径。初步探索了功率耦合机制，为实验系统优化及实验方案设定奠定理论基础。     

Unsteady wakes-secondary flow interactions in a high-lift low-pressure turbine cascade

Detailed experimental measurements were conducted to study the interactions between incoming wakes and endwall secondary flow in a high-lift Low-Pressure Turbine (LPT) cascade. All of the measurements were conducted in both the presence and absence of incoming wakes, and numerical analysis was performed to elucidate the flow mechanism. With increasing Reynolds number, the influence of the incoming wakes on suppressing the secondary flow gradually increased owing to the greater influence of incoming wakes on reducing the negative incidence angle at higher Reynolds numbers, leading to a lower blade loading near the leading edge and suppression of the Pressure Side (PS) leg of the horseshoe vortex. However, the effect of unsteady wakes on suppressing the profile losses gradually became weaker owing to the reduced size of the Suction Side (SS) separation bubble and increased mixing loss in the free-flow region at high Reynolds numbers. Incoming wakes clearly improved the aerodynamic performance of the low-pressure turbine cascade at low Reynolds numbers of 25,000 and 50,000. In contrast, at the high Reynolds number of 100,000, the profile loss at the midspan and mass-averaged total losses downstream of the cascade were higher in the presence of wakes than in the absence of wakes, and the unsteady wakes exerted a negative influence on the aerodynamic performance of the LPT cascade.