各向异性等离子体多重电流片的演化

应用二维三维量磁流体动力学方程组模拟，数值研究了压力各向异性等离子体周期多重电流片的演化，发现在β＜＜1时，微弱的各向异性仅仅使电流片中磁场重联的速度加快，而当β≈1（即热压与磁压相当时，微弱的垂直各向异性不仅大大加速了磁场重联的速度，还使重联位置发生显著变化。初始反对称的磁场重联位形逐步转化为对称的磁场位形，进而再演化为与初始反相的反对称磁场位形，场向各向异性和强的垂直各向异性都导致不稳定性增长率成倍增加。但在场各向异下磁场位形基本不变，只在电流片中心出现不规则扰动，同样参数下向向异性等离子体多重电流片中磁场重联形成的磁岛比各向同性多层电流要小。     

电流片的自发重联与无力场──磁岛的形成

采用二维三分量磁流体力学数值模拟方法,研究了线性无力场对电流片的自发重联的影响。结果表明:电流片的背景场(无力场)对重联的拓扑位形有重要的影响,当无力场因子达到一定值时,在模拟区域内出现大范围的磁岛,磁岛厚度随时间的发展而增加,而无力场因子较小时则没有磁岛形成。当无力场因子为零时,即回到普通的位场情形的单x线重联,位场情形磁场拓扑位形与Scholar的结果一致。      

无力场的触发重联和太阳耀斑

采用多步隐格式,对在瞬间形成的电流片的触发下的高剪切无力场的磁重联过程进行了数值模拟。磁重联首先在交界面处的非中性电流片区出现,然后向无力场区蔓延。在磁重联过程中,在无力场区形成一高温环状结构,物质向光球层流动。在高温环内侧的新喷发场区,物质向上流动。磁重联主要集中在初始电流片外侧的高剪切无力场区,高温环顶部的温度最高,位置基本固定。在磁重联的过程中,剪切磁场分量的空间梯度减小,无力场因子下降。     

初始电流片长宽比对非对称多重X线重联的影响

采用可压缩电阻性磁流体力学模型,研究了初始电流片的长宽比Lx对非对称多重X线磁场重联的影响.研究结果发现,当Lx超过一定的阈值时,非对称磁场重联的演化过程中将伴有多重X线重联发生.进一步地研究结果表明,当Lx越大,多重X线重联就越容易出现,相同时间间隔内所诱发的次级磁岛的尺寸也随之相应变大,不仅如此,次级磁岛的尺寸占整...     

2.5维自适应磁场重联MHD模式

磁雷诺数(R_m)是影响磁场重联的重要因素. 真实的物理环境中R_m往往很高, 例如, 在行星际空间和太阳日冕中R_m通常大于104量级. 高R_m条件下的磁重联表现出很多异常特性, 然而高R_m条件下的磁场重联数值模拟需要很高的时空分辨率, 否则很难分辨出重联过程中形成的薄电流片. 本文基于自适应软件包PARAMESH将并行自适应网格技术引入磁场重联数值模拟, 建立了一个2.5维自适应磁场重联MHD模式, 研究高磁雷诺数条件下重联的动态演化过程, 进而将不同磁雷诺数的参数进行对比研究. 结果表明, 该模式可以自动捕捉到磁场重联产生的奇性电流片, 高磁雷诺数条件下产生的慢激波结构可提供一种快速磁能释放机制.      

日冕冲浪形成的磁流体动力学模拟

应用二维时变可压缩磁流体动力学模拟,数值研究了双极-单极磁场中电阻撕裂模不稳定性引起的磁场重联过程,用于模拟日冕冲浪的形成.结果表明,在包含有三区——双极场、电流片和单极场的磁静力平衡初态下,双极场和单极场中的磁力线将会直接重联,磁场演变成鞭状(whip)结构.由弯曲磁力线支撑的等离子体团向上运动到最高位置后,逐渐下落和弥散.等离子体团上升速度可达到0.10v_A(v_A为双极场中的Alfv'én速度).模拟结果证实日冕冲浪的形成可能与双极-单极场中的磁场重联密切相关.      

电流片的形成

采用二维三分量的耗散磁流体力学模型,对由于光球层的剪切运动引起的四极磁场内电流片的形成过程进行数值模拟。结果表明,磁分隔线在磁场的演化过程中起关键作用。电流主要在分隔面和根部剪切梯度较大的磁环中形成,其中在分隔线上最强。磁分隔线可以被拉长,形成强电流片,随着速度剪切位置的不同,电流片的形成模式也不同,既可以形成与光球层垂直的电流片,也可以形成平行的电流片。      

等离子体彗尾动力学的数值研究

本文基于可压缩磁流体动力学模型,数值研究了尾瓣巾具有超Alfven速流动的等离子体彗尾的动力学特征。结果表明,等离子体片和尾瓣之间的剪切等离子体流动将会激发流动撕裂模不稳定性,引起彗尾等离子体片中发生磁场重联,形成磁岛和高密度的等离子体团。进而模拟了太阳风引起的局部驱动力对等离子体彗尾中磁场重联的影响,其特征时间远大于流动撕裂模。我们认为一些观测到的等离子体彗尾中的四块和彗尾截断事件可能主要与彗尾中剪切等离子体流动所引起的流动撕裂模不稳定性有关。      

2001年9月15日Cluster卫星对电流片的观测研究

2001年9月15日0430-0515 UT期间,Cluster卫星多次穿越磁尾电流片.由FGM、CIS等仪器获得了电流片磁场,粒子速度等数据变化情况,并得到了磁尾高速流的两次明显反转.本文采用GSM坐标系.利用求磁场空间梯度张量的方法获得了越尾电流,其电流密度的峰值为28nA/m2.并对电流片在垂直方向上的摆动和在晨昏方向的波动现象进行了分析.数据显示此时电流片为薄电流片并有一个变薄的趋势,其厚度大约为0.2-0.3 Re.磁场重联时地向流与尾向流均超过了1000 km/s,并测得了电流片的法线方向和运动速度,从而得到了磁尾电流片的结构和运动情况.      

空间等离子体压力各向异性对磁场重联的影响

基于二维时变可压缩磁流体动力学模拟,数值研究了等离子体压力各向异性对磁场重联的影响,发现一个小的压力各向异性(P_⊥=1.02P_//)即可大大加速磁场重联的发展,这可能是由于磁镜不稳定性与撕裂模不稳定性共同起作用.在P_⊥P_//的情况下,撕裂模受到抑制,电流片中不能形成大型磁岛.      

基于Object-FMA的软件代码审查方法

针对传统代码审查方法中使用的检查单具有的两个不足:缺少从代码到检查单缺陷记录的映射方式,检查单记录的代码缺陷类型不全而导致审查时代码缺陷的遗漏,提出了一种将代码进行层次划分,并使用对不同层次代码对象进行失效模式分析来代替使用传统检查单的代码审查方法.此外,提出了一种对象失效模式分析方法(Object-FMA,Object Fail-ure Modes Analysis),用于系统地分析代码对象的失效模式,为不同层次的代码对象建立失效模式库.在某型号软件的代码审查中对该方法进行了实例应用,结果表明此方法提供了一种从审查代码对象到失效模式库的映射方式,提高了代码审查的效率,审查时不仅能发现传统检查单中记录的代码词法、语法和常见语义的缺陷,还能有效地发现代码中人为疏忽或算法设计错误而引入的代码缺陷.     

Three-dimensional global simulation of multiple ICMEs’ interaction and propagation from the Sun to the heliosphere following the 25–28 October 2003 solar events

This study performs simulations of interplanetary coronal mass ejection (ICME) propagation in a realistic three-dimensional (3D) solar wind structure from the Sun to the Earth by using the newly developed hybrid code, HAFv.2+3DMHD. This model combines two simulation codes, Hakamada–Akasofu–Fry code version 2 (HAFv.2) and a fully 3D, time-dependent MHD simulation code. The solar wind structure is simulated out to 0.08 AU (18 Rs) from source surface maps using the HAFv.2 code. The outputs at 0.08 AU are then used to provide inputs for the lower boundary, at that location, of the 3D MHD code to calculate solar wind and its evolution to 1 AU and beyond. A dynamic disturbance, mimicking a particular flare’s energy output, is delivered to this non-uniform structure to model the evolution and interplanetary propagation of ICMEs (including their shocks). We then show the interaction between two ICMEs and the dynamic process during the overtaking of one shock by the other. The results show that both CMEs and heliosphere current sheet/plasma sheet were deformed by interacting with each other.     

First results of operational ionospheric dynamics prediction for the Brazilian Space Weather program

It is shown the development and preliminary results of operational ionosphere dynamics prediction system for the Brazilian Space Weather program. The system is based on the Sheffield University Plasmasphere–Ionosphere Model (SUPIM), a physics-based model computer code describing the distribution of ionization within the Earth mid to equatorial latitude ionosphere and plasmasphere, during geomagnetically quiet periods. The model outputs are given in a 2-dimensional plane aligned with Earth magnetic field lines, with fixed magnetic longitude coordinate. The code was adapted to provide the output in geographical coordinates. It was made referring to the Earth’s magnetic field as an eccentric dipole, using the approximation based on International Geomagnetic Reference Field (IGRF-11). During the system operation, several simulation runs are performed at different longitudes. The original code would not be able to run all simulations serially in reasonable time. So, a parallel version for the code was developed for enhancing the performance. After preliminary tests, it was frequently observed code instability, when negative ion temperatures or concentrations prevented the code from continuing its processing. After a detailed analysis, it was verified that most of these problems occurred due to concentration estimation of simulation points located at high altitudes, typically over 4000 km of altitude. In order to force convergence, an artificial exponential decay for ion–neutral collisional frequency was used above mentioned altitudes. This approach shown no significant difference from original code output, but improved substantially the code stability. In order to make operational system even more stable, the initial altitude and initial ion concentration values used on exponential decay equation are changed when convergence is not achieved, within pre-defined values. When all code runs end, the longitude of every point is then compared with its original reference station longitude, and differences are compensated by changing the simulation point time slot, in a temporal adjustment optimization. Then, an approximate neighbor searching technique was developed to obtain the ion concentration values in a regularly spaced grid, using inverse distance weighting (IDW) interpolation. A 3D grid containing ion and electron concentrations is generated for every hour of simulated day. Its spatial resolution is 1° of latitude per 1° of longitude per 10 km of altitude. The vertical total electron content (VTEC) is calculated from the grid, and plotted in a geographic map. An important feature that was implemented in the system is the capacity of combining observational data and simulation outputs to obtain more appropriate initial conditions to the ionosphere prediction. Newtonian relaxation method was used for this data assimilation process, where ionosonde data from four different locations in South America was used to improve the system accuracy. The whole process runs every day and predicts the VTEC values for South America region with almost 24 h ahead.     

卫星扩频应答机的伪码快捕方法

在低地球轨道通信系统中,由于系统的高动态,使扩频应答机接收的信号附加了几十千赫兹的多普勒频移,传统的伪码捕获方法很难完成捕获.对此提出了一种部分相关和FFT(Fast Fourier Transform)相结合的捕获方法,利用FFT补偿多普勒造成的相关峰损失,在搜索码相位的同时得到多普勒频移的估计,减少了捕获时间.给出了捕获系统的结构,分析了FFT对相关峰的补偿性能,给出了部分相关和FFT的参数设计方法,使用了自适应门限技术,以满足输入信号载噪比大幅变化的情况下实现伪码的捕获判决.理论分析和仿真结果表明:该方法可以在高多普勒频移和低载噪比环境下实现伪码的快速捕获,比其他快捕方法占用更少的FPGA(Field Programmable Gate Arrays)资源.     

利用GPS实现高轨卫星定位的抗远近效应算法

针对GPS应用于高轨卫星定位面临的远近效应问题, 提出利用正交相关的联合极大似然估计算法来对GPS信号进行二维搜索, 可以得到正确的码延时及多普勒估计. 极大似然估计算法首先利用滑行相关法对仿真信号中的强信号进行搜索, 得到强信号模型, 然后利用正交相关去掉强信号, 最终实现对弱信号的正确捕获. 本文对高轨卫星接收到的GPS信号进行了功率分析, 并建立了信号模型, 进行了算法仿真. 结果表明, 这种估计算法能够提高二维搜索性能, 有效解决远近效应问题.      

空间辐射场蒙特卡罗模拟计算方法研究

对利用蒙特卡罗方法对由银河宇宙射线引起的空间辐射场各成分进行计算的方法进行了调研,对计算模型的建立以及计算过程中通常使用的方差减小技术进行分析,给出了美国的Roesler等人利用FLUKA程序以及加拿大Anid等人利用MCNPX程序计算得到的由银河宇宙射线引起的空间辐射场各量值及其与实验结果的比较,验证了计算方法与计算模型的可靠性。对任意航线空间辐射场剂量分布预评估方法进行分析,给出了由银河宇宙射线引起的空间辐射场的基本特征。     

Analysis of the code and phase between-receiver inter-system biases of the overlapping frequencies for GPS/Galileo/BDS

The overlapping-frequency signals from different GNSS constellations are interoperable and can be integrated as one constellation in multi-GNSS positioning when inter-system bias (ISB) is properly disposed. The look-up table method for ISB calibration can enhance the model strength, maximize the number of integer-estimable ambiguities, and thus is preferred. However, the characteristics and magnitudes of the receiver code ISB and phase fractional ISB (F-ISB) are not well known and the wrong values of the biases can seriously degrade the positioning results. In this contribution, we first estimate the between-receiver code ISB and phase F-ISB of hundreds of the baselines up to around 25km in the European Permanent GNSS Network (EPN) and the Multi-GNSS Experiment (MGEX) for the overlapping frequencies L1-E1 (L1), L5-E5a (L5) and E5b-B2b (L7). The data collected from 1st January 2016 to 1st January 2019. Second, the receiver-type and firmware-version combinations for the receivers of Trimble, Leica, Javad, Septentrio and NovAtel are carefully classified. Results show that the Septentrio receivers have consistent code and phase ISB values for the three overlapping frequencies i.e. only one value for each frequency and no receivers are different. The Leica, Trimble and Javad receivers have two or more ISB values for at least one of the three frequencies. A few receivers with biases to the groups are also found and listed. Third, the code ISB and phase F-ISB of the groups are adjusted by the least-squares method. The root mean square errors (RMSE) of the least square adjustment are 0.240 m, 0.250 m and 0.200 m for code of L1, L5 and L7 frequencies, respectively, and are 0.0009 m, 0.0015 m and 0.0031 m for phase of L1, L5 and L7 frequencies, respectively. Finally, the effects of code ISB errors on code positing are investigated with the zero-baseline MAT1_MATZ. The distance root mean square error (DRMS) of L1-E1 code positioning can be reduced by 48.2% with 5 GPS and Galileo satellites and the DRMS degrades quickly when the code ISB error is larger.     

空间光学遥感器真空热试验工装模块化设计

针对空间光学遥感器在空间环境地面模拟试验中舱板温度分布不均的问题,运用模块化思想对某型号空间光学遥感器工装的舱板结构开展了优化设计。将舱板按照温度分布进行分块,提出模块化拼装和模块化热控2种设计方案,模块化拼装是对舱板及其表面加热片进行独立分块划分,模块化热控则是将舱板视为整体,仅对其表面的加热片进行分块设计。结果表明:模块化热控的舱板平均温度偏差为0.205 K,低于未模块化设计的0.87 K和模块化拼装的0.30 K,提高了舱板的温度均匀性。同时,模块化拼装改善了舱板温度分布,使得符合热控要求的测点比例由34.8%提高到96.7%,但独立划分的模块之间仍存在一定温差;模块化热控则消除了模块间的温差,将符合热控要求的测点比例进一步提高到100%,完全满足热控要求。      

Extending multipath hemispherical model to account for time-varying receiver code biases

Multipath effects on code observables account for one of the major error sources in high-accuracy Global Positioning System (GPS)-based positioning, atmosphere sounding and timing applications. The multipath hemispherical model (MHM) represents one of the most widely used methods of mitigating code multipath effects by taking advantage of their spatial repeatability. The use of MHM usually assumes that the receiver code biases (RCBs) are time-invariant; however, this assumption is not always valid, as RCBs and linear combinations thereof (differential code biases, for instance) have long been found to be time-varying over a period of one day. In this contribution, we propose an extended multipath hemispherical model (EMHM) that is capable of mitigating the code multipath effects in the presence of time-varying RCBs. Consequently, the proposed EMHM has two advantages. First, the EMHM gives rise to code multipath corrections with improved reliability because it addresses the intraday variability of RCBs. Second, more interestingly, the EMHM allows easy and effective calibration of short-term temporal variations, if any, in the RCB on each frequency. These advantages are hopeful to benefit GPS code-related applications.     

能量粒子进入磁层的数值模拟研究

来自行星际的能量粒子进入磁层的过程一直是磁层研究的热点,也是空间天气研究关注的焦点之一,而以往的工作大都集中在讨论不同行星际条件下粒子进入磁层的统计规律.本文在T89磁层模式的基础上,建立了太阳能量粒子进入磁层运动的程序,计算了几种不同能量的质子从不同方位角入射到磁层的运动轨迹.模拟结果表明,能量粒子沿着磁力线方向才可能穿越磁层到达地球表面,越偏离这个方向,则越早被反弹.能量越高的粒子进入到地球磁层上空的角度范围越大,但仍然只有沿着磁力线入射的粒子才能到达地球表面.这些结果与理论预言是一致的.