Negative VTEC with spherical harmonic expansion model under different solar activity conditions

The Earth’s ionosphere can be described by a spherical harmonic (SH) expansion up to a specific degree. However, there exist negative vertical total electron content (VTEC) values in the global ionosphere map (GIM) with the SH expansion model. In this contribution, we specifically investigated the negative VTEC values that are induced by the SH expansion model and validated the performance of the inequality-constrained least squares (ICLS) method in eliminating the negative VTEC values. The GPS data from 2004 to 2017 was selected to cover one solar cycle and the experiments under different solar activity conditions were analyzed. The results in our work show that the occurrence of the negative VTEC values is attributed to the deficiency of the SH expansion model when the VTEC itself is small instead of the unevenly distribution of the GNSS stations. The negative VTEC values appear periodically in the temporal domain, showing apparently one year and half year periods. During one year, two peaks in June and December can be observed in the time series of the negative VTEC values. The number of negative VTEC values in June is obvious larger than that in December. During one solar cycle, the number of negative VTEC values under quiet solar activity condition is obvious larger than that under strong solar activity condition. In the spatial domain, the appearance of the negative VTEC values is strongly related with the movement of the subsolar point. In the latitude of the subsolar point has the largest magnitude, the negative values will appear on the opposite hemisphere and the further from the subsolar point the more negative values. The maximum number of the negative VTEC values in the southern hemisphere appears in June, while the peak value in the northern hemisphere appears in December. The maximum number of negative VTEC values in the southern hemisphere is generally larger than that in the northern hemisphere. In addition, the negative VTEC values are distributed both at middle latitude and high latitude in the southern hemisphere, while they are mainly distributed at high latitude in the northern hemisphere. When the ICLS method is used, the negative VTEC values can be eliminated efficiently and it has nearly no influence on the positive VTEC values. The ICLS method can also improve the receiver’s differential code bias (DCB) and significantly decrease the unreasonable negative slant TEC (STEC) values along the lines of sight. Using the final GIM product of the Jet Propulsion Laboratory (JPLG) as a reference, the root mean square (RMS) of the ICLS solution shows maximum 25%, 20% and 45% improvement relative to the least squares (LS) solution at northern high latitude, southern middle latitude and southern high latitude, respectively.     

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.     

化学非平衡流动对超燃冲压发动机尾喷管性能的影响

采用RNGk-ε湍流模型和有限体积方法数值求解有组分的守恒形式Navier-Stokes方程,针对设计的超声速燃烧冲压发动机单斜面膨胀喷管,采用氢氧七组元八反应模型,数值模拟不同高度、不同来流马赫数条件下的喷管流场和性能.计算结果表明,喷管内的进一步燃烧对燃烧室起到了一定的补燃作用,对喷管性能影响较大.化学非平衡状态下,喷管推力特性和升力特性显著提高.超燃冲压发动机尾喷管的性能研究,应考虑化学非平衡流动的影响.      

基于模型的转子系统不平衡量的估计

采用基于模型的方法对转子系统的不平衡量进行估计,首先记录初始状态下的振动信号,然后根据实时测得的不平衡响应信号,得出差值振动信号,联立转子系统的动力学方程,得出不平衡等效载荷,根据等效载荷,得出不平衡的位置、相位及大小.采用模态扩展技术实现由有限几个点的振动响应获取所需点上的响应.由于噪声的影响,提出了基于高斯小波的降噪方法.仿真与实验结果表明,这一方法能够比较准确地估计出转子系统的不平衡量.      

Propagation of non-linear ion-sound pulses in dusty plasma waveguides

There are a lot of objects in space associated with dusty plasma inclusions. Such inclusions may bear a prolonged shape and behave as waveguides for ion-sound waves. In the case of space plasmas, the dust particles can possess both negative charge, due to electron attachment, and positive one, due to photoionization. In this paper the propagation of linear and non-linear ion-sound wave pulses in the dusty plasma waveguides, possessing positive charge, is studied. It has been demonstrated that non-linear dynamics of baseband pulse propagation in plasma waveguide possesses essentially non-solitonic behavior. Namely, propagation of a long ion-sound pulse leads to an excitation of a shock-like wave but not a stable localized nonlinear pulse. Also, when a Korteveg–de Vries (KdV) soliton is incident onto the dusty plasma waveguide, some part of the soliton energy is captured by the waveguide and transformed into a multi-pulse structure. Additionally, an interaction of dusty plasma inclusions with KdV soliton can lead to the occurrence of transverse instabilities of the soliton and its eventual destruction.     

Height and time variability of planetary wave activity

The height–season and year-to-year regularities of parameters of first and second spatial harmonics determine the structure of the stratosphere and mesosphere circulation and its variability. In the period 1992–2002 at heights 0–55 km, the amplitudes and phases of the first and second spatial harmonics in the field of temperature, geopotential height, zonal and meridional wind were calculated by the method of harmonic decomposition. Dispersion (standard or mean square deviation) of their day-to-day and year-to-year variations was calculated by their wavelength constants. Height and season patterns of variability have been estimated. The difference in height–longitude variability for wave numbers m = 1 and 2 has been discovered. At the same time, the intensity of wave disturbances for m = 1 is less than for m = 2 excluding the polar areas, where a significant variability appears at the heights 0–55 km. There is also a tendency for the intensity of year-to-year variations to decrease in comparison with day-to-day variations. In cold and warm periods the amplitude of perturbation waves with m = 2 both for day-to-day and year-to-year variations is greater than for waves with m = 1. Transient height areas in the interval of 20–30 km are more distinct for day-to-day variations of polar area.     

Spectroscopic Constraints on Models of Ion-cyclotron Resonance Heating in the Polar Solar Corona

Using empirical velocity distributions derived from UVCS and SUMER ultraviolet spectroscopy, we construct theoretical models of anisotropic ion temperatures in the polar solar corona. The primary energy deposition mechanism we investigate is the dissipation of high frequency (10-10000 Hz) ion-cyclotron resonant Alfvén waves which can heat and accelerate ions differently depending on their charge and mass. We find that it is possible to explain the observed high perpendicular temperatures and strong anisotropies with relatively small amplitudes for the resonant waves. There is suggestive evidence for steepening of the Alfvén wave spectrum between the coronal base and the largest heights observed spectroscopically. Because the ion-cyclotron wave dissipation is rapid, even for minor ions like O⁵⁺, the observed extended heating seems to demand a constantly replenished population of waves over several solar radii. This indicates that the waves are generated gradually throughout the wind rather than propagated up from the base of the corona. This revised version was published online in June 2006 with corrections to the Cover Date.     

Low-Frequency Waves and Instabilities in Collisionless Plasmas

We present a theoretical overview of low-frequency waves and instabilities in collisionless, multi-component plasmas with gyrotropic ( ) thermal pressure. We show that the complete dispersion relation can be obtained in the framework of a mixed magnetohydrodynamic (MHD)-kinetic formalism, which uses the MHD mass, momentum, and induction equations, together with the kinetically corrected version of the double-adiabatic equations of state. The complete dispersion relation contains not only the three standard modes (fast, slow, and Alfvén) from double-adiabatic MHD, but also the mirror mode from kinetic theory. We examine the stability properties of these four modes, firstly in the case of a uniform medium, and secondly in the case of a stratified and rotating medium. We also discuss the connections with the quasi-interchange modes (interchange and translation) often referred to in the context of magnetospheric physics.     

三维四向编织复合材料等效热特性数值分析和试验研究

 采用“米”字型枝状胞体有限元模型和试验方法对三维四向编织复合材料的整体等效热膨胀系数和等效热传导系数进行了分析并将计算结果和试验值进行了比较。研究表明枝状胞体模型能较为真实的反映三维四向编织复合材料的结构构形;有限元方法在分析热膨胀、热传导方面具有较好的精度;三维四向编织复合材料编织方向的热膨胀系数随着纤维体积比的增加而减小,随着编织角的增加而减小;热传导系数随着纤维体积比的增加而增大,随着编织角的增大而增大。     

爆炸冲击载荷下典型机身结构动响应及破坏

确定满足适航要求的最小风险炸弹位置,必须研究爆炸冲击载荷下机身结构动响应及破坏模式。以某型飞机典型机身结构为研究对象,采用LS-DYNA商用软件,建立了爆炸冲击载荷下机身典型结构动响应数值模型。采用控制变量法分析了爆炸物当量、爆炸冲击距离以及爆炸冲击位置对典型机身结构动响应及破坏模式的影响,同时研究了损伤后典型机身结构的剩余强度。研究结果表明,造成机身结构有效破坏的爆炸物临界当量与爆炸冲击距离密切相关;爆炸冲击距离对典型机身结构损伤及剩余强度影响不明显;典型机身结构筋条位置对剩余强度影响较大。在此基础上,提出了表征剩余强度的无量纲系数,并建立了剩余强度无量纲系数与爆炸物当量及爆炸冲击距离之间的函数关系。