南京地区斜向加热电离层效应数值模拟

与垂直加热相比,斜向加热电离层具有灵活性高、影响范围广和便于实际操作等优势.在非偏区考虑欧姆吸收,基于电子能量方程和连续性方程构建电波斜向加热低电离层的物理模型,并利用IRI-2007和NRLMSISE-00经验模型提供的背景参数对南京地区斜向加热低电离层进行数值模拟,对比不同加热条件下电子温度和电子密度的扰动情况.研究结果表明:电波加热效果随入射仰角和有效辐射功率的增大而增大;电子温度和电子密度增幅随电波频率增大而减小;X波模比O波模造成的电子温度扰动幅度和电子密度扰动幅度更大,同时X波模比O波模能更快地使电子温度和电子密度达到稳定状态;一定范围内较小仰角、较低频率、较大有效辐射功率的电波能使电子密度更快达到稳定,后两者还能加快电子温度达到稳定的过程;电子温度达到稳定所需时间随入射电波仰角呈单峰变化,仰角为62°时达到最大.      

基于线性调频波的有源对消隐身仿真及分析

有源对消隐身技术由于其适用性广、灵活主动及不破坏目标外形等特点而受到重视。针对其技术要求中的快速反应、延时较短的特点,提出了一种基于线性调频（LFM）波的对消系统方案。该方案对入射波采样后进行相位和频率的调制,然后直接转发,使其匹配滤波输出与原入射波的回波输出进行对消;探讨了有源对消技术的电磁学原理并以理想导电无限长圆柱为散射源,对平面波与柱面波的对消效果进行了仿真验证,得出不同视角、平面波与柱面波的夹角都会对对消效果产生影响的结论,且在两列波夹角为0°时达到最佳;推导了进行相位与频率调制的对消后信号输出公式,并利用MATLAB对其进行了仿真验证,结果表明,对消后输出信号峰值由时延脉宽之比决定,带宽对输出信号主峰的宽窄程度有直接影响。      

不同入射角度下等离子鞘套中C波段电磁波传输特性

为解决临近空间通信黑障问题,研究了不同入射角对等离子体鞘套中C波段（4~8GHz）电磁波传输特性的影响,分别采用均匀等离子体模型与高斯分布非均匀等离子体模型,根据电磁波在分层介质面的反射和透射传播原理,对不同角度下C波段通信电磁波的传输性能进行计算和对比分析.结果表明,入射波频率的增加以及入射角的减小有利于降低C波段的衰减值,提高透射率,使通信电磁波更有效地穿透等离子体鞘套,为解决临近空间通信黑障问题提供了理论依据.      

Oblique sounding of the ionosphere by powerful wave beams

The article is devoted to modeling the impact on the ionosphere powerful obliquely incident wave beam. The basis of this analysis will be orbital variational principle for the intense wave beams-generalization of Fermat’s principle to the case of a nonlinear medium (, ,  and ). Under the influence of a powerful wave beam appears manageable the additional stratification of the ionospheric layer F2. Explicit expressions show how the properties of the test beam, with a shifted frequency, released in the same direction as the beam depend on the intensity of a powerful beam and the frequency shift.     

VLF波从大气层到低电离层的传输特性分析

基于传播矩阵法计算了均匀半空间电离层的反透射系数, 同时解Booker复系数四次方程得到电离层的复折射指数, 分别研究了电离层反射透射系数及折射指数随VLF频段入射电波频率、入射角和地磁倾角、电离层电子浓度及碰撞频率的变化规律. 计算结果表明, 在VLF频段, 垂直电偶极子辐射的横磁(Traverse Magnetic)波更易透射进入电离层, 而水平电偶极子辐射的横电(Traverse Electric) 波易被限制在地-电离层波导内来回反射. 电离层电子密度较低时 (如夜间), 在高纬度地区, 观测到地震电离层VLF异常的概率更大. 当考虑地磁场 的影响时, 电离层将允许地震辐射的超低频(Ultra Low Frequency ,ULF)/甚低频 (Very Low Frequency, VLF)部分的电磁波透射进入电离层, 这一点已有很多卫星观测事实为证, 但其进一步的物理机制尚需深入研究.      

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

Alfvén gravity waves in viscous solar plasma: Effect of the background flow

We study the Propagation and damping properties of Alfvén gravity waves in the presence of the vertical magnetic field in the viscous solar plasma under influence of the background flow by deriving a fourth order dispersion relation in terms of the Doppler shifted frequency. We derive the dispersion relation under WKB and Boussinesq approximation. We study the damping of Alfvén gravity waves for the wave frequencies less than and greater than the Brunt-Väisälä frequency. We find that the Brunt-Väisälä frequency divides the frequency ranges where the weakly or strongly damped oscillations occur. The background flow exhibits a strong effect on weakly damped oscillations and a weak effect on the strongly damped oscillations. We also notice that the damping of both the strong and weakly damped oscillations depend on the Brunt-Väisälä frequency and wave number. The effect of the background flow is also being governed by the Brunt-Väisälä frequency and wave number. We also study the properties of gravity wave mode after filtering the Alfvén wave mode by minimizing the magnetic field and noticed that the background flow shows a very strong effect on the gravity wave mode.     

左旋寻常的极光千米波与辐射带高能电子相互作用的参数化研究

参数化研究了在L=4.5外辐射带区域左旋寻常（Left-hand Ordinary mode, L-O模）的极光千米波与高能电子的相互作用，定量计算了左旋寻常的极光千米波在不同峰值频率、传播角分布和纬度分布条件下的弹跳平均投掷角扩散系数、动量扩散系数和投掷角–动量交叉扩散系数。计算结果表明，动量扩散系数一般比投掷角扩散系数大100倍左右，说明动量扩散在L-O模与电子间的相互作用中起主导作用。随着传播角范围的改变，扩散系数均未发生明显变化，这说明L-O模与电子相互作用的扩散系数对传播角范围的依赖性很小。此外，扩散系数会因为L-O模纬度分布的改变而发生剧烈变化，该结果与之前得到的纬度分布对右旋奇异的极光千米波的影响是一致的。通过参数化研究结果表明，适当条件下，L-O模可能会显著影响外辐射带高能电子的动力学过程。      

Study of waves in the magnetotail region with cluster and DSP

The study of the neutral sheet is of fundamental importance in understanding the dynamics of the Earth’s magnetosphere. From the earliest observation of the magnetotail, it has been found that the neutral sheet frequently appears to be in motion due to changing solar wind conditions and geomagnetic activity. Multiple crossings of the neutral sheet by spacecraft have been attributed to a flapping motion of the neutral sheet in the north–south direction, a wavy profile either along the magnetotail or the dawn–dusk direction. Cluster observations have revealed that the flapping motions of the Earth’s magnetotail are of internal origin and that kink-like waves are emitted from the central part of the tail and propagate toward the tail flanks. This flapping motion is shown here to propagate at an angle of ∼45° with x_GSM. A possible assumption that the flapping could be created by a wake travelling away from a fast flow in the current sheet is rejected. Other waves in the magnetotail are found in the ULF range. One conjunction event between Cluster and DoubleStar TC1 is presented where all spacecraft show ULF wave activity at a period of approximately 5 min during fast Earthward flow. These waves are shown to be Kelvin–Helmholtz waves on the boundaries of the flow channel. Calculations show that the conversion of flow energy into magnetic energy through the Kelvin–Helmholtz instability can contribute to a significant part of flow breaking between Cluster and DoubleStar TC1.     

多离子空间等离子体中快磁声波与粒子的回旋共振特性

快磁声波是空间等离子体中一种接近垂直传播的右旋极化电磁波,能够在等离子体层内外传播.快磁声波与带电粒子的回旋共振相互作用能够导致高能电子随机加速和投掷角扩散、能量质子投掷角扩散等,从而影响辐射带高能带电粒子的动态过程.分别基于完整的色散关系和高密度近似的色散关系,在不同空间等离子体条件下研究多离子空间等离子体中不同传播角的快磁声波色散曲线,并计算了快磁声波与H+,He+和O+离子的最小共振能量.结果表明,当传播角较小时,采用高密度近似与采用完整色散关系计算的离子最小共振能量没有太大差别.在中低密度中强磁场空间等离子体中,传播角≥ 88°时高密度近似色散关系会带来很大的误差,因此应利用完整色散关系计算最小共振能量.      

非均匀双流等离子体中波传播耦合波动方程

我们利用麦克斯韦方程组和矩方程组,研究了非均匀双流等离子体在电流源作用下,激发波的耦合过程。给出了三模[P_e,P_i,E_o]耦合的普遍方程▽2[P]-[C]·▽[P]+[K2][P]=[D][▽·(ω_p2E_o)]。作为一特例,本文首先研究了均匀、有损失(v≠0)双流等离子体中,等离子体波的耦合过程。指出v≠0,即使在均匀等离子体中,也将产生电子等离子体波和离子等离子体波间耦合。但不出现光学模和等离子体模间耦合。      

Electron acceleration and parallel electric fields due to kinetic Alfvén waves in plasma with similar thermal and Alfvén speeds

We investigate electron acceleration due to shear Alfvén waves in a collissionless plasma for plasma parameters typical of 4–5R_E radial distance from the Earth along auroral field lines. Recent observational work has motivated this study, which explores the plasma regime where the thermal velocity of the electrons is similar to the Alfvén speed of the plasma, encouraging Landau resonance for electrons in the wave fields. We use a self-consistent kinetic simulation model to follow the evolution of the electrons as they interact with a short-duration wave pulse, which allows us to determine the parallel electric field of the shear Alfvén wave due to both electron inertia and electron pressure effects. The simulation demonstrates that electrons can be accelerated to keV energies in a modest amplitude sub-second period wave. We compare the parallel electric field obtained from the simulation with those provided by fluid approximations.     

哨声波的有质动力

本文导出哨声波有质动力的一个公式。该公式表明,当ω<Ω_e/2时,有质动力倾向于把电子推向强波场区,从而在波的传播路径上形成密度脊;当ω>Ω_e/2时情况正相反,有质动力倾向于在波的传播路径上形成密度槽。这就表明,weibel的结果[1]只是当ω《Ω_e时本文所得到的结果的特殊情形。      

Case study of inclined sporadic E layers in the Earth’s ionosphere observed by CHAMP/GPS radio occultations: Coupling between the tilted plasma layers and internal waves

We have used the radio occultation (RO) satellite data CHAMP/GPS (Challenging Minisatellite Payload/Global Positioning System) for studying the ionosphere of the Earth. A method for deriving the parameters of ionospheric structures is based upon an analysis of the RO signal variations in the phase path and intensity. This method allows one to estimate the spatial displacement of a plasma layer with respect to the ray perigee, and to determine the layer inclination and height correction values. In this paper, we focus on the case study of inclined sporadic E (E_s) layers in the high-latitude ionosphere based on available CHAMP RO data. Assuming that the internal gravity waves (IGWs) with the phase-fronts parallel to the ionization layer surfaces are responsible for the tilt angles of sporadic plasma layers, we have developed a new technique for determining the parameters of IGWs linked with the inclined E_s structures. A small-scale internal wave may be modulating initially horizontal E_s layer in height and causing a direction of the plasma density gradient to be rotated and aligned with that of the wave propagation vector k. The results of determination of the intrinsic wave frequency and period, vertical and horizontal wavelengths, intrinsic vertical and horizontal phase speeds, and other characteristics of IGWs under study are presented and discussed.     

自适应非奇异快速终端滑模固定时间收敛制导律

针对机动目标的末制导拦截问题,设计了一种带攻击角度约束的非奇异快速终端滑模固定时间收敛制导律。与有限时间收敛终端滑模制导律相比,所提制导律能够确保弹目视线（LOS）角和弹目视线角速率在固定时间内是收敛的,并且收敛时间是独立于制导系统初始条件的,可以根据制导律参数预先给定。构造了一种新型的非奇异快速终端滑模面,有效解决了奇异性问题,同时通过合理地改变滑模面与弹目视线角跟踪误差的趋近律指数,使得制导系统比现有的固定时间收敛控制具有更快的收敛速率。此外,设计了一种自适应律,针对目标机动引起的未知扰动进行估计,使得制导律的设计无需预先知道任何关于目标机动的信息。通过仿真实验验证了所提制导律能够使导弹成功拦截机动目标,并且与现有制导律相比,具有更快的系统收敛速率、更高的拦截精度及更短的拦截时间。      

波状扰动激发的重力波波包在中层大气中能量传播和转换特性的数值研究

采用谱配置方法对具有高斯分布的波状扰动激发重力波过程中的能量传播和转换特性进行数值研究．模拟结果表明：初始时刻单一的水平风、垂直风、密度或温度的波状扰动都可以激发出两支重力波波包，其中一文向上传播，而另一支则向下传播；在相同的波参数下，4种波状扰动所激发的重力波波包的能量传播路径几乎完全一致，与线性理论给出的射线路径却有些差异．进一步的数值分析则表明：4种波状扰动转化为波动的特征时间基本相等，约为2个波动周期，但是在这个转换过程中，不同扰动的能量转换效率却差别很大．     

基于传输矩阵法的等离子体鞘套传输特性分析

引入分层模型传输矩阵法,分析高速临近空间飞行器等离子体鞘套传输特性.通过对典型飞行器模型进行纳维-斯托克斯方程组数值求解,仿真了等离子体鞘套的流场电子密度分布,分析了电磁波在等离子体鞘套内传输时的衰减常数,给出信号传输的反射系数和透射系数.基于传输矩阵法的仿真结果与理论分析结果和时域有限差分法仿真分析结果有良好的一致性,同时,该方法有效降低了计算复杂度,适用于任意电子密度分布和厚度的等离子体鞘套传输特性的数值求解.      

带攻击角度约束的自适应终端滑模导引律

针对某些导弹在对目标进行打击时需要满足零脱靶量和攻击角度约束的要求,首先基于终端滑模控制和有限时间控制理论,改进了一种快速收敛的非奇异终端滑模函数,用于设计滑模面,结合自适应指数趋近律,提出了一种自适应非奇异终端滑模控制方法,解决了传统终端滑模控制中存在的奇异问题,并使状态变量在有限时间内快速收敛到平衡点。然后将所提方法用于导引律的设计,提出了一种带攻击角度约束的自适应非奇异和有限时间收敛导引律,实现了导弹对脱靶量和攻击角度约束的要求;采用有限时间控制理论对该导引律的收敛特性进行了分析,证明了制导系统状态的全局有限时间快速收敛特性。与传统的非奇异终端滑模导引律相比,本文所提导引律能够在更短的时间内以更小的脱靶量和更高精度的攻击角度对目标实施打击。最后进行了大量的对比仿真实验,仿真结果验证了所提导引律的有效性。      

Ion dynamics associated with Alfven wave in the near-Earth magnetotail: Two-dimensional global hybrid simulation

Ion dynamics in the near-Earth magnetotail region is examined during periods of fast Earthward flow with a two-dimensional (2-D) global-scale hybrid simulation. The simulation shows that shear Alfven waves are generated at x ∼ −10R_E, where the strong earthward flow is arrested by the dipole field, and propagate along field lines from the equator to both southern and northern polar ionosphere. Non-gyrotropic ion velocity distributions occur where the large-amplitude Alfven waves are dominant. The simulation indicates that the Alfven waves are generated by interaction of the fast earthward flow with the stationary near-Earth plasma. Beam ions are found to be pitch-angle scattered and trapped in the wave field, leading to the non-gyrotropic ion distributions in the high-latitude plasma sheet boundary. In addition, significant particle heating and acceleration are found to occur behind the dipolarization front due to the effect of wave turbulence.