新建正激变换器运行理论

用两种与传统不同的电路模型表征正激变换器.其一用LC串联谐振实现磁复位功能(这一模型中没有磁复位绕组,并确认开关管存在与其并联的电容).其二主要用磁复位绕组实现磁复位功能(这一模型中有磁复位绕组,但也确认开关管存在与其并联的电容).用相平面分析法分析这些模型,得到了正激变换器运行时电路和磁路特征参量的解析表达式.本文的分析和解析结果以及后续对这两种正激变换器的全面分析被称为"新建正激变换器运行理论".     

磁壳参数L与磁暴Dst指数和行星际条件的关系

用磁坐标L-/A来描述地球近地空间粒子特性和卫星位置是近年来空间物理探测研究和数据分析中的一个新趋势.利用T96磁场模型计算了L值,并比较了在地球磁层剧烈活动期间和不同行星际条件下,用偶极子模型,国际参照磁场(IGRF)模型,和T96磁场模型这三种方式计算得到的地球表面L-A磁坐标之间的区别.在地磁纬度大于30°时偶极子近似和IGRF磁场模型计算得到的L值差别开始增大.在地磁纬度大于50°时,用IGRF磁场模型和T96磁场模型计算得到的L值差别开始增大.由于T96磁场模型引入了行星际磁场南北和东西分量,计算的L值包含了行星际条件的影响,并具有了随地方时变化的特性.本工作对于辐射带粒子动态模式的建立,以及正确理解卫星磁坐标位置等具有重要意义.     

脉宽调制器中一类振荡器的全面解析结果

用分析法全面讨论脉宽调制器（IC：1842/3/4/5,1842A/3A/4A/5A）中振荡器的振荡特征参数（占空比D和频率f）与外接电阻R和电容C的关系.我们得到以R、C为自变量表达D和f的函数,以D、f为自变量求R、C的函数.并根据这些函数的计算结果绘制描述D、f、R、C（1F）4个参数之间的关系曲线.这些结果,不仅提升了对这类振荡器的认知,而且可供设计时直接引用和参考.     

中国地区MAGSAT卫星标量和矢量磁异常图

使用MAGSAT卫星资料,编绘中国及邻近地区的卫星磁异常图(10°N—60°N,70°E—140°E)。为了提取地壳异常场,必须从观测资料中消除主磁场、磁层场,感应场和电离层场。本文选用GSFC(12/83)地磁场模型,消除主磁场,磁层场和感应场。为了进一步消除电离层场的影响,改善相邻轨道磁异常的一致性,用纬度的二次多项式拟合每一条轨道的初始地壳异常,并从初始地壳磁异常值减去这一拟合值,从而获得最后的地壳异常值(△X,△Y,△Z,△F)。把研究地区分成1°×1°的网格,将每个网格内的磁异常进行平均,以获得网格中心点的异常值。根据这些网点值绘制中国及邻近地区的卫星标量和矢量磁异常图。等值线间隔为2nT。为了检验卫星磁异常图的可靠性,将卫星资料分成黎明组、黄昏组和联合组,分别进行处理和绘制卫星磁异常图。结果表明:根据上述资料绘制的卫星磁异常图有很好的一致性。中国地区的卫星磁异常值位于—10nT—10nT。卫星磁异常与区域构造特征有较好的一致性;塔里木地台、扬子地台和中朝地台是正磁异常区,西藏高原是负磁异常区。卫星磁异常反映出下地壳磁化强度的横向不均匀性。     

磁云边界层与间断的观测研究

利用WIND飞船的磁场和等离子体观测数据,分析了1995年2月至2003年8月之间82个磁云中的边界层事件.我们认为(1)磁云边界层中方向间断DD(Directional Discontinuity)类型中RD(Rotational Discontinuity),TD(Tangential Discontinuity),ED(Either Discontinuity),ND(NeitherDiscontinuity)的分布为37:18:44:1(%),与背景太阳风中的分布51:12:35:2(%)不同;主要区别在于RD与TD的比例变化. (2)磁云边界层的内外边界切向间断的比例很高,占总数的20%,而且两侧的密度和温度平均相对变化分别为|△N/N|=0.24和|△T/T|=0.19,大于边界层内部的平均值,显示了更多的切向间断特征;此外,磁云边界层中方向间断的出现频次约为太阳风中的1.87倍. (3)磁云边界层中方向间断的法向分布在θ-φ平面中不是随机分布,而是以θ=-24.90°,φ=217.49°为几何中心,主要是指向远离太阳的方向,而不是简单的各向同性分布.初步结果表明,间断是磁云边界层中的重要结构,它有着不同于背景太阳风间断类型比,为诊断磁云边界层的形成机理提供了依据.     

从属函数在地磁扰动预报研究中的初步应用

根据1966-1982年期间有关太阳耀斑、行星际激波和地磁扰动的观测资料而建立的从属函数,对1984-1985年间的行星际闪烁观测中能证认出的耀斑-激波所引起的地磁扰动作了预报试验。结果表明:(1)磁扰开始时间预报的相对误差,δT/T≤10%的事件数为20个,占总事件数的50%,δT/T≤20%的事件占总事件数的70%以上;(2)磁扰幅度(ΣKp)大小的预报,其相对误差δΣKp/ΣKp≤30%的事件数为32个,占总事件数的80%,而δΣKp/ΣKp≥60%仅占15%.本文方法显示了一定潜力,有待从聚类分析方面进一步深入。      

基于NOAA卫星数据的辐射带质子投掷角分布研究

利用NOAA-12卫星数据对空间环境平静时期太阳同步轨道处辐射带质子投掷角分布进行了研究. 根据投掷角分布的经验公式,计算出90°投掷角的质子方向强度和各向异性指数n. 质子投掷角分布按n的取值范围可分为三类,即90°峰值分布、平顶分布和蝴蝶形分布. 观测证实,对于辐射带质子,三种分布类型均存在并且具有明显的空间区域特征. 在内辐射带边缘地区90°峰值分布占主要优势;在外辐射带高L值区域,90°峰值分布明显减少,平顶分布和蝴蝶形分布逐渐增多. 针对90°峰值分布,研究了质子强度各向异性的区域分布特征,对于内辐射带区域,n值随L值的增大而增大,对于外辐射带,n值表现为逐渐下降的趋势. 为了研究质子投掷角分布对磁地方时的依赖关系,分析了能量为250～800keV的质子在两个不同磁地方时范围的投掷角分布规律. 结果显示,在内辐射带,质子强度的投掷角分布相对稳定,随磁地方时的变化并不显著;而在外辐射带的高L值区域,质子强度的投掷角分布随磁地方时变化明显,与磁地方时之间有明显的依赖关系.     

磁纬10°以下哨声路径纬度和所需抬升因子积分中值的估算方法

本文采用国际地磁参考场(1980IGRF,n=8)磁力线经验公式和三层模式电离层,结合哨声波沿磁力线传播所需电子浓度的横向梯度,给出了磁纬10°以下哨声路径纬度和所需抬升因子积分中值的估算方法;从而得出:(1)哨声色散值D不仅与电离层N_mF2正相关,而且与h_mF2也有明显相关性.当路径纬度ф₉₀≤10.5°(IGRF,n=8,下同)时,D与h_mF2负相关;>12时,正相关;在10.5°—12°之间时,正或负相关取决于路径顶点高度与h_m之差值.估算表明,三亚(磁纬7.04°N)收到的哨声中路径纬度ф₉₀≤10.5。的占94.5%,D与h_mF2负相关,这与观测结果一致.(2)所需抬升因子积分中值一般在7—31%之间.      

用于GMA的新型永磁偏置闭合磁路

利用ANSYS有限元软件进行仿真模拟,对超磁致伸缩执行器(GMA,Giant Magnetostrictive Actuator)的磁路结构进行了分析和优化.根据磁致伸缩棒内磁场均匀性好,场外漏磁小的设计目标,利用极性相反的双永磁补偿原理,提出了一种新型的内置式永磁组合偏置结构,使得沿磁致伸缩棒长方向磁场分布均匀的同时,偏置磁场和激励磁场都形成良好的闭合磁路.ANSYS分析结果表明,棒中偏置磁场不均匀度为3.51%,激励磁场不均匀度为8.73%,均满足不均匀度小于10%的设计目标;在执行器指定位置(距离GMA 7 cm)的总漏磁场强度为30.4 A/m,符合实际应用的标准要求(<80 A/m).     

NASA成功发射Lucy小行星探测任务

2021年10月16日,NASA的Lucy小行星探测任务成功发射升空,将在未来12年探测一颗主带小行星和7颗特洛伊族小行星. Lucy是NASA发现计划(Discovery)的第13项独立任务,于2017年1月获批,主要目标是探测多个木星特洛伊族小行星.此类小天体是早期太阳系的遗迹,与木星共用轨道,分为位于木星轨道前方(L4)和后方(L5)的两群,围绕太阳运行.     

双轴模拟式太阳敏感器测试链误差分析及修正

基于四象限硅光电池,双轴模拟式太阳敏感器（简称双轴模太）可实现太阳矢量两轴方向角的同时测量。硅光电池各象限光生电流的采集精度直接决定了敏感器的工作性能。然而,测试链环节中存在的不一致性误差会引起双轴模太角度测量偏差。为此,基于数值仿真,分析并建立了硅光电池各象限响应率、各象限光生电流的电压转换、信号放大、A/D转换系数及暗电流等环节偏差的影响模型,进而对各环节的综合偏差进行标定和修正,从而在保证敏感器测量精度的同时,降低对测试链环节的要求。测试结果表明,在不改变测试链性能的情况下,双轴模太的测量精度由修正前的2.05° (α轴,3σ)和1.94° (β轴,3σ)提高到0.28° (α轴,3σ)和0.26° (β轴,3σ)。     

Estimation of the bias of the Minimum Variance technique in the determination of magnetic clouds global quantities and orientation

Magnetic clouds (MCs) are highly magnetized plasma structures that have a low proton temperature and a magnetic field vector that rotates when seen by a heliospheric observer. More than 25 years of observations of magnetic and plasma properties of MCs at 1 AU have provided significant knowledge of their magnetic structure. However, because in situ observations only give information along the trajectory of the spacecraft, their real 3D magnetic configuration remains still partially unknown. We generate a set of synthetic clouds, exploring the space of parameters that represents the possible orientations and minimum distances of the satellite trajectory to the cloud axis, p. The synthetic clouds have a local cylindrical symmetry and a linear force-free magnetic configuration. From the analysis of synthetic clouds, we quantify the errors introduced in the determination of the orientation/size (and, consequently, of the global magnetohydrodynamic quantities) by the Minimum Variance method when p is not zero.     

How to address the accuracy of empirical magnetic field models?

Empirical magnetic field models are compared with high-altitude magnetic field measurements and results from an MHD simulation. Comparison of the T96 model and observations from GOES-8 and GOES-9 shows that if the observed solar wind and IMF parameters are used to compute the model field, the model field is more stretched than the observed field. On the other hand, if measurements made by one spacecraft are used to find the model parameters that give a best-fit field at that location, the RMS error can be reduced also at the other spacecraft four hours away in local time. Comparison of T96 model and MHD simulation results shows that the empirical models have a thinner current sheet than the MHD simulation, but that the lobe field values are quite similar to each other. Furthermore, a comparison of an event-oriented, modified T89 model and MHD simulation by Pulkkinen et al. [2000] reveals that if the empirical model is constructed by fitting to in-situ measurements, the resulting model is very similar to the MHD simulation magnetic field. These results indicate that an efficient method of utilizing the present-day empirical models is to select model parameters based on measurements from a few individual points.     

Estimation of long-term cosmic ray intensity variation in near future and prediction of their contribution in expected global climate change

On the basis of results obtained in our paper [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], we determine: the dimension of the Heliosphere (modulation region), radial diffusion coefficient and other parameters of convection–diffusion; drift mechanisms of long-term variations of cosmic ray (CR) dependence on particle energy; level of solar activity (SA); and generally, the solar magnetic field. We obtain this important information on the basis of CR and SA data in the past, taking into account the theory of convection–diffusion and global drift modulation of galactic CR in the Heliosphere. By using these results and other regularly published predictions of expected SA variation in the near future, as well as predictions of the next SA cycle, we may make predictions of long-term cosmic ray intensity variation expected in the near future (up to 10–12 years). In [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], properties of connections between long-term variation in CR intensity and some part of a global climate change were estimated, controlled by solar activity through CR. We show that in this way we may make predictions of some part of a global climate change expected in the near future (up to 10–12 years and maybe more, depending upon the period during which definite predictions of SA can be made), controlled by solar activity through CR. In this case, estimations of expected long-term changes in the planetary distribution of cutoff rigidities, which also influence CR intensity, as well as CR-influenced effects on global climate variation, become important.     

Estimating inter-system biases for tightly combined Galileo/BDS/GPS RTK

The overlapping carrier frequencies L1/E1, L5/E5a and B2/E5b from GPS/Galileo/BDS allow inter-system double-differencing of observations, which shows a clear advantage over differencing of the observations of each constellation independently. However, the inter-system biases destroy the integer nature of the inter-system double-differencing ambiguities. Two methods of direct rounding and parameter estimation are used to determine the ISB value. By analyzing data collected from Curtin University from 2015 to 2018, the phase and code inter-system bias (ISB) are related to the receiver type, firmware version and the selected overlapping frequency. Upgrade of receiver firmware version results in changes of ISB values. For example, the upgrade of Javad firmware in Dec, 15, 2017 causes the difference of 0.5 cycles ISB between BDS GEO and non-GEO satellites. By comparing the three dynamic models which include white noise process, random walk process, and random constant in the parameter estimation method, the ISB determined by the random constant model is consistent with the value obtained by the direct rounding method. After the calibration of ISBs, the performances of tightly combined positioning are assessed. The success rate of ambiguity resolution and accuracy of positioning for the tight combination (TC) are significantly improved in comparison with that for the loose combination (LC) over short baselines. For L5/E5a, on which only few satellites can be observed, the maximum increase in success rates of ambiguity resolution can reach 31.7%, i.e., from 54.9% of LC to larger than 86.6% of TC, and the positioning accuracies can even be increased by 0.13 m, i.e., from 0.208 of LC to 0.074 m of TC in East direction for the mix-receiver TRIMBLE NETR9-SEPT POLARX4 in 2018.     

Towards multi-dimensional space weather modeling for energetic oxygen ions in the earth''s inner magnetosphere: Equilibrium configuration

We report the first 3+1 dimensional model development for energetic atomic oxygen ions in the Earth's radiation belts. Energetic Oxygen ions cans be supplied to the Earth's Inner magnetosphere from the sun (as a component of solar wind and solar energetic particles), from anomalous cosmic rays, and from acceleration processes acting on ionospheric atomic oxygen ions. We have built a multi-dimensional oxygen ion model in the following free parameters: geomagnetic L-shell, the magnetic moment, the second adiabatic invariant, and the discrete charge state number. Quiet time, steady state oxygen ion distributions have been obtained numerically from an assumed outer radiation zone boundary condition at L=7, average values of the radial diffusion coefficients, and standard values for the exospheric neutral densities due to the MSIS-86 upper atmosphere and exosphere neutral thermal particle density model. Average distributions of free electrons in the plasmasphere were also assumed with a mean plasmapause location just beyond L=4. We included the six lowest ionic charge states of atomic oxygen (¹⁶O) based on an existing charge exchange cross section compilation by Spjeldvik and Fritz (1978). Computed oxygen ion distributions include the resulting equilibrium structure of energy oxygen ions between 10 KeV and 100 MeV.     

Modelling studies of ionospheric convection in northern and southern polar regions

The realistic model of Quegan et al. has been used to investigate the convection paths of ionospheric plasma at 300 km altitude, for different polar cap radii and in both hemispheres. Taking the Northern magnetic dip pole to be at a co-latitude of 11° and the Southern magnetic dip pole at a co-latitude of 23°, these paths are presented in a Sun-Earth frame, with the position of the Earth's axis fixed as it is on 21 March, as polar plots centred on the magnetic pole. There are marked hemispheric differences between 13 and 23 L.T., particularly near the stagnation region at 18 to 21 L.T., but only minor differences between 00 and 12 L.T., when the radius of the polar cap exceeds 12°. For a smaller polar cap, the differences between the hemispheres are small at all local times. The time taken to perform a complete circuit is most dependent on the polar cap radius, and most variable - between 15 and 36 h - for convection paths starting near 60° latitude. The time that plasma convecting from noon to near midnight across the Northern polar cap spends within the 10° co-latitude circle increases from 6 h, for a polar cap radius of 10°, to 11.5 h at 17°. These results are compared and contrasted with other model calculation results and with some ground-based and satellite observations of plasma densities at high latitudes.