镜模波识别方法研究及其在火星磁鞘中的应用

镜模波是温度各向异性等离子体中的一种波动结构,根据磁场和离子分布及波动特性可以进行识别.本文对比了只使用磁场数据与同时使用磁场及离子数据两种识别方法,分析了两类方法的特点.只使用磁场数据的方法基于磁场强度变化大、方向沿背景磁场的特征,通常使用磁场强度的波动幅度ΔB/|B|以及磁场变化方向与背景磁场的夹角θ_min,θ_max作为参数;同时使用磁场及粒子数据的方法利用的是磁场纵波特性、总压平衡和波动在等离子体坐标系下静止的特征.使用两种方法对MAVEN卫星在火星磁鞘内的数据进行识别,结果表明在某些情况下,只使用磁场数据会导致对镜模波的误判.通过研究改变上述参数阈值时识别结果的变化,发现当θ_min> 40°,θ_max < 40°,ΔB/|B|> 80%时,只用磁场数据可取得较好的识别效果.      

压电柔性结构振动的鲁棒控制

针对具有多面体型参数不确定性和外界扰动的柔性悬臂梁系统,研究其振动的鲁棒控制问题.基于压电柔性悬臂梁系统的状态方程,利用一种扩展LMI(Linear Matrix Inequality),对柔性悬臂梁的振动提出了改进的 H₂/H_∞ 综合状态反馈鲁棒控制方法.此方法针对多面体型参数不确定性系统引入不同的Lyapunov函数,降低了原来 H₂/H_∞ 综合鲁棒控制方法中由于使用一个公共Lyapunov函数而产生的保守性.通过Matlab分别采用改进和未改进的 H₂/H_∞ 鲁棒控制方法对柔性悬臂梁系统的振动控制进行了数值仿真.仿真结果表明,2种方法都使柔性悬臂梁的振动得到了良好控制,而改进方法的控制效果更好.     

基于准静态拉伸试验的临界CTOA测量

为了寻求更为简单、准确的临界裂纹尖端张开角（CTOA_c）测量方法以推广CTOA准则的应用,本文对材料CTOA_c的测量方法进行了研究。根据CTOA的性质及G?teborgs Kungliga Segel S?llskap（GKSS）提出的CTOA_c估算方法导出了两种不同裂纹尖端张开位移（CTOD）定义之间的关系,并进一步得到一种基于准静态拉伸试验的CTOA_c测量方法,该方法与GKSS提出的估算方法相等效但更加简单。为了验证准静态拉伸试验法的准确性,对7B04 M（T）试样进行了静态拉伸试验,并使用上述两种方法估算7B04的CTOA_c,结果表明两种方法估算出的CTOA_c仅相差1%。同时分别采用平面应变核模型、平面应力模型对试样进行了弹塑性有限元模拟,模拟结果与试验结果相吻合。      

喷流制冷的高超声速后台阶流场气动光学效应研究

近年来,临近空间的军事应用价值受到了越来越多的关注。为了研究高超声速飞行器典型后台阶结构的气动光学问题,采用研究稀薄气体动力学的直接模拟蒙特卡洛算法(DSMC)研究了临近空间不同喷流控制下的后台阶流动特性,并采用光线追踪方法研究了不同流动控制下流场的气动光学效应。计算结果表明,随着飞行高度增加,壁面热流将增加,同时气动光效应将变弱。在四种喷流冷却方式中,喷流方式A的冷却效果最好,同时没有增加气动光学效应,B、D方式具有一定的冷却效应,但使得气动光学效应更加严重。     

不完全加速退化的三参数幂模型线性化方法

针对加速度计零偏和标度因数不完全加速退化数据的非线性预测问题,提出一种加速度计加速退化的三参数幂函数模型线性化方法.该方法假设三参数幂函数退化轨迹模型的幂指数与加速应力无关,考虑到不同应力水平的试验时间相当时,较高应力下的参数退化速率较快,其退化轨迹与较低应力水平的相比,退化量更加接近性能参数的失效阈值,而低应力水平下退化量与失效阈值差距较大,因此最高应力水平下的退化轨迹更能反映加速度计性能参数的长期非线性退化规律.据此进一步提出了加速退化数据三参数幂函数模型幂指数确定的均方相关系数最大化原则,该原则通过高应力水平试验数据确定模型的幂指数,由此实现三参数幂函数模型的线性化,克服由于初值敏感而导致模型预测时稳健性较差的问题.为加速退化轨迹非线性模型的预测提供了一种工程化方法.     

基于线性回归近似替代模型的分布式卫星系统不确定性分析方法

分布式卫星系统在空间科学领域具有广泛的应用前景,同时也对概念设计阶段的分析工作提出了更高的要求.在分布式卫星系统概念设计阶段,评价不确定性参数对于最终探测效能的影响有着重要的工程应用价值.传统的不确定性分析方法存在解析困难、数值模拟计算效率低且耦合关系表征不明确的缺点.本文结合概念设计阶段不确定性分析的特点和需求,提出一种线性回归近似替代模型,将不确定性参数对于探测效能影响的结果计算降维为一个线性组合系数的求解问题.利用一个空间多点探测物理场分界面任务作为典型事例进行仿真验证,结果表明,相比传统数值模拟方法,本文方法在计算效率上具有明显优势且对关键不确定性参数有一定识别作用.      

低纬(海南)地区电离层偶发E层的潮汐类周期研究

基于低纬(海南)台站(19.5°N,109.1°E)电离层测高仪在2004年观测到的偶发E层(Es)临界频率(f₀E_s)和虚高(h'E_s)数据,利用谱分析方法,给出了中国低纬地区Es层在全年和不同季节内的短周期变化特性.对全年数据谱分析表明,Es层除具有非常强的24h周期外,还具有显著的12h和8h周期.对不同季节内两个参数的谱分析进一步表明,8h周期主要出现在f₀E_s参数的春秋分季和夏末及在h'E_s参数的几乎所有季节;6h周期主要出现在f₀E_s参数的11-12月.海南低纬Es层中具有显著的6h和8h周期特性,特别是这两个周期还能同时出现在冬季,这是以往研究中不曾发现的新结果.初步分析认为,Es层的6h和8h周期极可能是分别由1/4日和1/3日潮汐所造成的.      

发动机燃气喷流红外辐射场的数值模拟

在高温流动中热辐射与流场是耦合的,但在发动机燃气喷流流场的温度范围内可以在流场计算中忽略热辐射的影响,于是流场与辐射场的模型方程可解耦.为此提出了模拟飞行器绕流流场红外辐射的解耦模型以及相应的算法.该算法首先采用总变差减小TVD(Total Variation Diminishing)格式对流场进行模拟,再引入有限体积概念,应用所取得的流场参数,在同一组网格上计算流场中的红外辐射特征.以轴对称喷管内N₂,O₂,CO₂,H₂O,CO,HCl,H₂等7组分高度欠膨胀冻结流与均匀绕流干扰流场为算例进行了验证性的模拟计算,并与已有试验和计算数据进行了对比,表明该算法是可靠的,可在较大幅度地降低计算量的情况下给出满足工程需要的结果.     

一种基于白谱法的电离层天气扰动指数

基于一种电离层扰动提取方法——白谱法,利用IGS提供的电离层TEC网格数据,获得电离层J_s指数、J_r指数和J_p指数,分别反映单站、纬度圈（沿经度积分）及行星际尺度下的电离层天气扰动状态.在2015年3月的一次磁暴过程中,J_s指数、J_r指数及J_p指数均很好地反映出电离层响应地磁暴的过程,磁暴前后J_p指数与Dst指数相关系数达到-0.72;J_s图从二维角度很好地表征了电离层天气的扰动过程.在此基础上,统计分析了2011——2014年J_p指数与Dst指数的相关性,结果表明:限定J_p≥2,J_p指数与对应时间Dst指数的相关系数为-0.67;限定J_p≥3,二者相关系数更高,达到-0.87.通过分析不同J_p指数阈值下不同等级磁暴的次数,发现J_p指数可以很好地反映磁暴下的电离层整体扰动,为指示电离层天气状态提供了可能的参数.      

数字通信系统信噪比参数的测量方法比较

针对数字通信系统中常见的信噪比参数E_b/N₀,本文从功率测量、功率谱密度测量及矢量分析三个角度出发,分别给出了该参数的测量原理与方法,并通过具体测量实验比较了三种测量方法的各优缺点,对于该参数的实际工程测量具有一定指导意义。     

Inconsistency of magnetic field and plasma velocity variations in the distant plasma sheet: Violation of the “frozen-in” criterion?

Measurements of the magnetic field and low energy plasma by the GEOTAIL spacecraft have been used to study the relationship between variations of the plasma velocity and of the magnetic field in the distant (100–200 R_E) and middle (40–80 R_E) tail. The analysis was carried out separately for the tail lobes and the plasma sheet. It is shown that the absolute values of the magnetic field and plasma velocity, as well as their corresponding components (V_X and B_X, V_Y and B_Y, V_Z and B_Z), are linearly connected in the tail lobes. In the plasma sheet, however, the plasma velocity and the magnetic field do not seem to be related to one another. The distant plasma sheet seems to be in a regime of turbulence. The diffusion coefficients estimated from our data set of the velocity parameters in the plasma sheet are in good agreement with the theoretical predictions of Antonova and Ovchinnikov (1996, 1999).     

Steady and oscillatory convection in vertical cylinders heated from below. Numerical simulation of asymmetric flow regimes

Three-dimensional buoyancy driven flows in vertical cylinders are simulated with a finite difference technique. Asymmetric convective flows are studied for three values of the (length-to-radius) aspect ratio A. One result concerns the occurrence of asymmetric convection in a flat A=0.485 cylinder for Pr=6.7, near the threshold and up to 3Ra_c (where Ra_c is the critical Rayleigh number and Pr is the Prandtl number). Complex (steady and time-dependent) supercritical regimes have been simulated in cylinders of aspect ratios A=2 and A=4 for Pr=0.02 and in ranges of Ra up to 8Ra_c and 6Ra_c, respectively. The flow patterns are analysed graphically and discussed with respect to a weakly nonlinear analysis and to experiments. In a particular case, A=4 and Pr=0.02, an oscillatory motion has been obtained and some features of the flow pattern are shown during a period.     

Plasma convection in the low-latitude mantle of the near-earth magnetosphere

A detailed analysis of Geotail observations on June 1, 1997 provides evidence for the appearance of the plasma mantle at low-latitudes, Z_GSM 3 R_E for southward interplanetary magnetic field (IMF), and clarifies the characteristics of the plasma flow in the low-latitude mantle. The tailward flow speed in the mantle is 80 160 km/s, smaller than the flow speed (300 km/s) in the adjacent magnetosheath. The component of the velocity perpendicular to the magnetic field lines is predominantly poleward at speeds up to 60 km/s, and the convection tends to be stronger for |;B_Y|<|;B_Z| than for IMF |B_Y|>|B_Z|;. This polewardconvection can be interpreted as being driven by a J × B force operating on the field lines.     

Interplanetary shocks and geomagnetic activity during solar maximum (2000) and solar minimum (1995–1996)

Plasma and magnetic field parameter variations through fast forward interplanetary shocks were correlated with the peak geomagnetic activity index Dst in a period from 0 to 3 days after the shock, during solar maximum (2000) and solar minimum (1995–1996). Solar wind speed (V) and total magnetic field (B_t) were the parameters with higher correlations with peak Dst index. The correlation coefficients were higher during solar minimum (r² = 56% for V and 39% for B_t) than during solar maximum (r² = 15% for V and 12% for B_t). A statistical distribution of geomagnetic activity levels following interplanetary shocks was obtained. It was observed that during solar maximum, 36% and 28% of interplanetary shocks were followed by intense (Dst  −100 nT) and moderate (−50  Dst < −100 nT) geomagnetic activity, whereas during solar minimum 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity. It can be concluded that the upstream/downstream variations of V and B_t through the shocks were the parameters better correlated with geomagnetic activity level, and during solar maximum a higher relative number of interplanetary shocks can be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50% to be followed by intense/moderate geomagnetic activity.     

Collisionless reconnection in the magnetotail

Recent developments regarding collisionless reconnection in current sheets with a finite normal magnetic field component (B_z) are reviewed. In 2-D x, z configurations the ion tearing mode is stabilized by the electron compressibility. When the y dependence is included, cross-field current instabilities can be excited. Of these, the drift kink mode appears to be particularly important. 3-D electromagnetic particle simulations indicate that this mode can act as the precursor to the growth of tearing modes and subsequent reconnection.     

IMF effect on ionospheric trough occurrence at equinoxes

Previous observations have shown that there is a relationship between the F region trough and both B_z and B_y components of the interplanetary magnetic field (IMF). Since IMF governs the polar cap convection, we investigate here if this relationship can be explained by means of polar cap convection. The study is limited to equinox seasons. The poleward and equatorward edges of the trough are determined from satellite tomographic observations and their locations are plotted in magnetic coordinates together with the convection pattern given by Papitashvili and Rich [Papitashvili, V.O., Rich, F.J. High-latitude ionospheric convection models derived from DMSP ion drift observations and parameterized by the IMF strength and direction. J. Geophys. Res. 107, 2002, doi:10.1029/2001JA000264] using IMF measurements coincident with trough observations. The results indicate a close relationship between the troughs and convection. Most of the troughs are seen within the dusk cell and the pattern of trough observations rotates with the convection pattern, when B_y changes its sign. More dayside troughs are observed when B_z is negative than in the opposite case, i.e. fast convective flow favours the dayside trough occurrence. Nightside troughs are observed more frequently when B_y is negative. In both evening and morning sectors the trough is situated close to the edges of convection cells, which partly contradicts previous results showing that the troughs are associated with the convection reversal. It is concluded that plasma convection has an important role in trough generation, although the effect of a strong electric field and other mechanisms like precipitation certainly have a role of their own.     

The model of turbulent plasma sheet during IMF Bz > 0

Theory of the plasma sheet with medium-scale developed turbulence gives the possibility to explain the main processes of plasma sheet bifurcation and theta-aurora formation during IMF B_z > 0. The model suggests that during IMF B_z > 0 small bulge structure in the plasma sheet center is formed. The polarization of the bulge due to dawnward electron motion and duskward ion motion decreases the large-scale electric field in the bulge region. The decrease of the large-scale field in the conditions of constant coefficient of diffusion leads to the bulge growth. The results of plasma sheet bifurcation and theta-aurora formation modelling are presented.     

The collective gyration of a heavy ion cloud in a magnetized plasma

In both the ionospheric barium injection experiments CRIT I and CRIT II, a long-duration oscillation was seen with a frequency close to the gyro frequency of barium and a time duration of about one second. A model for the phenomenon which was proposed for the CRIT 1 experiment is here compared to the results from CRIT II which made a much more complete set of measurements. The model follows the motion of a low-β ion cloud through a larger ambient plasma. The internal field of the model is close to antiparallel to the injection direction v_i but slightly tilted towards the self-polarization direction E_P = −v_ixB. As the ions move across the magnetic field, the space charge is continuously neutralized by magnetic-field aligned electron currents from the ambient ionosphere, drawn by the divergence in the perpendicular electric field. These currents give a perturbation of the magnetic field related to the electric field perturbation by ΔE/ΔB ≈ V_a. The model predictions agree quite well with the observed vector directions, field strengths, and decay times of the electric and magnetic fields in CRIT II. The possibility to extend the model to the active region, where the ions are produced in this type of self-ionizing injection experiments, is discussed.     

INTERBALL statistical study of ion flow fluctuations in the plasma sheet

We use ion distribution measurements with CORALL instrument on-board the INTERBALL/Tail spacecraft to study plasma flows in the mid-tail (−9> X> −27 R_E) plasma sheet. Three velocity components computed every 2 minutes exhibit two types of velocity variations: Earthward bursty bulk flows (BBFs) and random flow fluctuations. Their properties are in a good agreement with the observations of the ISEE-2 spacecraft (Borovsky et al., 1997). The INTERBALL/Tail spacecraft configuration favors measurements of V_z component, in contrast to previous experiments in which only V_x and V_y were measured reliably. In the outer part of the plasma sheet V_y and V_z fluctuations were close to each other (variances σ(V_y) and σ(V_z) were about 160 and 110 km/s, respectively), but in the inner part at the dusk flank amplitude of V_y fluctuations increased and was 2 times higher than that of V_z component. This asymmetry of fluctuations should be taken into account during modern theoretical analysis and simulations.     

Large geomagnetic storms of extreme solar event periods in solar cycle 23

During extreme solar events such as big flares or/and energetic coronal mass ejections (CMEs) high energy particles are accelerated by the shocks formed in front of fast interplanetary coronal mass ejections (ICMEs). The ICMEs (and their sheaths) also give rise to large geomagnetic storms which have significant effects on the Earth’s environment and human life. Around 14 solar cosmic ray ground level enhancement (GLE) events in solar cycle 23 we examined the cosmic ray variation, solar wind speed, ions density, interplanetary magnetic field, and geomagnetic disturbance storm time index (D_st). We found that all but one of GLEs are always followed by a geomagnetic storm with D_st  −50 nT within 1–5 days later. Most(10/14) geomagnetic storms have D_st index  −100  nT therefore generally belong to strong geomagnetic storms. This suggests that GLE event prediction of geomagnetic storms is 93% for moderate storms and 71% for large storms when geomagnetic storms preceded by GLEs. All D_st depressions are associated with cosmic ray decreases which occur nearly simultaneously with geomagnetic storms. We also investigated the interplanetary plasma features. Most geomagnetic storm correspond significant periods of southward B_z and in close to 80% of the cases that the B_z was first northward then turning southward after storm sudden commencement (SSC). Plasma flow speed, ion number density and interplanetary plasma temperature near 1 AU also have a peak at interplanetary shock arrival. Solar cause and energetic particle signatures of large geomagnetic storms and a possible prediction scheme are discussed.