反演航天器在轨瞬态外热流的导热反问题

获得航天器在轨飞行过程中的外热流数据对于研究热控涂层在轨退化规律、各种空间因素对热控产品的影响以及航天器姿轨控发动机羽流热效应都有非常重要的意义,然而直接测量热流存在很多困难,因此可以通过求解导热反问题得到满足一定精度的结果.首先,通过研究利用航天器设备在轨遥测温度值反演出航天器在轨瞬态外热流的导热反问题方法,建立了反演航天器在轨瞬态外热流的数学模型,采用共轭梯度法求解导热反问题并从物理概念角度改进了共轭梯度法的迭代过程以增加其抗不适定性;然后构造了两组能够代表目前大多数地球轨道航天器以及深空探测航天器在轨吸收外热流变化的数值试验对共轭梯度法的反演效果进行了检验.除阶跃变化位置以外反演值与真实值的最大相对偏差为2.9%,反演效果非常好;对于阶跃变化位置的吸收外热流在对反演结果进行分析处理后也能够得到较好的反演结果.      

强磁暴、能量粒子暴与热层大气密度涨落之间的相关关系

利用1997-2007年由GOES8, GOES11和GOES12星载高能粒子探测器在地球同步轨道高度上所探测到的高能质子和高能电子通量探测数据以及高度560km左右星载大气密度探测器所得的热层大气密度探测数据, 统计分析了强地磁扰动、高能粒子通量跃变和热层大气密度涨落之间的相关关系, 初步获得强地磁扰动期间, 地球同步轨道(外辐射带外环)均出现了增幅大于三个数量级的高能质子通量(尤其是E>1MeV)强增强现象, 随后热 层大气密度强烈上涨, 表明三者之间是正相关关系. 在时间上地球同步轨道高能质子通量强增强现象先于日均Ap值(地磁活动程度)上涨约一天左右, 而热层大气密度强涨落现象又明显滞后于强地磁扰动事件.      

基于神经网络及深层充电的电子通量反演模型

为了实现对空间高能电子通量的估计及航天器深层充放电的风险评估,基于深层充电和空间电子环境的关联性,利用人工神经网络（ANN）建立了一种由深层充电反演空间高能电子环境的模型。以深层充电探测电流密度及电子能量作为模型输入,电子通量作为模型输出,使用AE9模型对神经网络进行训练,将神经网络的MSE降低到了0.04122,并使用Giove A卫星的深层充电探测数据及GOES卫星的电子通量探测数据验证了模型反演电子环境的准确性。同时对由探测电流计算航天器典型介质材料最大内电场的神经网络模型进行了研究,以实现对航天器内充电风险实时评估。     

FY-3A卫星与NOOA系列卫星高能带电粒子实测结果的比较

FY-3A卫星是运行于830 km高度的太阳同步轨道气象卫星,其搭载的空间环境监测器可观测3～300 MeV的高能质子和0.15～5.70 MeV的高能电子.FY-3A卫星在轨工作期间,太阳活动处于由谷年向峰年过渡期,空间环境非常平静,探测结果显示3～300 MeV的高能质子分布主要集中在南大西洋辐射带异常区,0.15～5.70 MeV的高能电子分布区域除南大西洋异常区外,还分布在南北两极高纬区域.FY-3A与NOAA卫星测量结果反映出带电粒子强度及分布区域随投掷角变化的空间各向异性特征.本文在充分考虑了带电粒子时间、空间分布差异以及比对探测器之间自身设计差异的前提下,经过归一化处理后,首次对两颗卫星同期探测结果进行相关性分析,验证了两颗卫星相同时空条件下高能带电粒子通量分布的一致性;说明FY-3A空间环境监测器不仅具备空间带电粒子辐射监测能力,且探测结果有效可靠,可用作辐射带环境数据源的组成部分,为发展新的模型,深入研究辐射带高能粒子的分布、起源和传输等提供新的观测依据.     

FY-3A卫星与NOAA系列卫星高能带电粒子实测结果的比较

FY-3A卫星是运行于830 km高度的太阳同步轨道气象卫星, 其搭载的空间环境监测器可观测3~300 MeV的高能质子和0.15~5.70 MeV的高能电子. FY-3A卫星在轨工作期间, 太阳活动处于由谷年向峰年过渡期, 空间环境非常平静, 探测结果显示3~300 MeV的高能质子分布主要集中在南大西洋辐射带异常区, 0.15~5.70 MeV的高能电子分布区域除南大西洋异常区外, 还分布在南北两极高纬区域. FY-3A与NOAA卫星测量结果反映出带电粒子强度及分布区域随投掷角变化的空间各向异性特征. 本文在充分考虑了带电粒子时间、空间分布差异以及比对探测器之间自身设计差异的前提下, 经过归一化处理后, 首次对两颗卫星同期探测结果进行相关性分析, 验证了两颗卫星相同时空条件下高能带电粒子通量分布的一致性; 说明FY-3A空间环境监测器不仅具备空间带电粒子辐射监测能力, 且探测结果有效可靠, 可用作辐射带环境数据源的组成部分, 为发展新的模型, 深入研究辐射带高能粒子的分布、起源和传输等提供新的观测依据.      

基于BP神经网络的GEO等离子体环境参数反演分析

空间等离子体环境诱发的表面充电效应会对航天器运行产生干扰,严重时将导致太阳电池等部件失效。通过神经网络反演方法,以GEO环境中介质表面充电电位曲线作为输入,在双峰麦克斯韦分布假设下,可以逆向得到高能峰的等离子体参数。分析了GEO等离子体环境参数对表面充电电位曲线的影响,表明高能峰在充电过程中起决定性作用;其次通过MATLAB搭建BP神经网络,采用COMSOL计算得到多组充电曲线进行网络训练和反演计算,得到等离子体密度反演的平均相对误差为0.42%,温度反演的平均相对误差为0.03%,整体误差在0.1%～5.6%。结果表明,采用神经网络对等离子体环境进行反演具有可行性,该方法可以作为空间等离子体环境探测结果的对比参考和航天器非探测点表面电位计算的输入条件。     

由地面宇宙线中子强度反演空间质子和重离子的分布

采用线性回归的方法研究了太阳活动平静时期空间高能质子和重离子微分通量与地面OULU宇宙线台站中子强度之间的相关性,利用地面中子强度数据来反演空间高能质子和重离子通量.从质子和重离子的能谱出发,结合OULU台站中子强度的数据,提出了由地面中子强度数据反演空间高能质子和重离子微分通量的新方法.文中以GOES卫星上350～420 MeV,420～510 MeV,510～700 MeV,＞700 MeV(P8～P11)四能道的高能质子和ACE卫星上的元素C为例,并将反演的数据与测量结果比较,二者符合较好.为了探索这种反演方法在空间辐射环境描述中的应用,同时利用2006年11月和12月质子通量的反演结果对反演方法进行了验证,证明利用这一反演方法可以从地面中子强度的数据很好地反演出空间质子和重离子的微分通量.      

探测-1发现陌生空间现象

□□中国探测-1卫星于1月3日开始向地面传回数据。科学家对数据分析后发现,探测-1首次探测到了一种陌生的空间现象,即在其轨道运行空间内有一种新的高能电子分布规律,这在以往的空间探测中从未发现过。探测数据显示,在距离地球比较远的时候,高能电子探测器探测到的电子分布状况     

VAP发现内辐射带相对论电子数量低于预期

<正>NASA网站3月15日报道,利用范艾伦探测器(VAP)数据,科学家发现内辐射带的相对论电子数量低于预期,可能有助于设计更廉价的航天器。相关论文发表在Journal of Geophysical Research:Space Physics上。以往的空间任务无法区分内辐射带中的高能质子和电子。利用VAP搭载的磁电子和离子光谱仪(MagEIS),科学家首次排除了高能质子的影响,并发现内辐射带大多数时候不包含相对论电子,这与     

基于嫦娥一号高能粒子数据的地球磁层屏蔽效应研究

月球绕地球运行轨道约有1/4位于地球磁层内,因此,地球磁层是否会为月球轨道附近高能粒子提供足够的磁场屏蔽对于探索月球活动具有重要影响.嫦娥一号是中国首颗绕月人造卫星,其绕月飞行的工作轨道距离月球表面200 km.通过对嫦娥一号高能粒子探测器(HPD)的探测数据进行分析,比较了当月球位于地球磁层内外6个不同能道(能量范围4～400 MeV)时质子通量的变化,发现当月球位于地球磁层内时,这些能道的质子通量并没有发生显著减少,结果表明地球磁层不能为月球轨道附近高能粒子提供显著的磁屏蔽.     

Results from the ESA SREM monitors and comparison with existing radiation belt models

The Standard Radiation Environment Monitor (SREM) is a simple particle detector developed for wide application on ESA satellites. It measures high-energy protons and electrons of the space environment with a 20° angular resolution and limited spectral information. Of the ten SREMs that have been manufactured, four have so far flown. The first model on STRV-1c functioned well until an early spacecraft failure. The other three are on-board, the ESA spacecraft INTEGRAL, ROSETTA and PROBA-1. Another model is flying on GIOVE-B, launched in April 2008 with three L-2 science missions to follow: both Herschel and Planck in 2008, and GAIA in 2011). The diverse orbits of these spacecraft and the common calibration of the monitors provides a unique dataset covering a wide range of B-L* space, providing a direct comparison of the radiation levels in the belts at different locations, and the effects of geomagnetic shielding. Data from the PROBA/SREM and INTEGRAL/IREM are compared with existing radiation belt models.     

SHAKING: Adjusted spherical harmonics adding KrigING method for near real-time ionospheric modeling with multi-GNSS observations

Precise positioning based on Global Navigation Satellite System (GNSS) technique requires high accuracy ionospheric total electron content (TEC) correction models to account for the ionospheric path delay errors. We present an adjusted Spherical Harmonics Adding KrigING method (SHAKING) approach for regional ionospheric vertical TEC (VTEC) modeling in real time. In the proposed SHAKING method, the VTEC information over the sparse observation data area is extrapolated by the Adjusted Spherical Harmonic (ASH) function, and the boundary distortion in regional VTEC modeling is corrected by the stochastic VTEC estimated using Kriging interpolation. Using real-time GPS, GLONASS and BDS-2/3 data streams of the Crust Movement Observation Network of China (CMONOC), the SHAKING-based regional ionospheric VTEC maps are re-constructed over China and its boundary regions. Compared to GNSS VTECs derived from the independent stations, the quality of SHAKING solution improves by 13–31% and 6–33% with respect to the ASH-only solution during high and low geomagnetic periods, respectively. Compared to the inverse distance weighting (IDW) generated result, significant quality improved of SHAKING-based VTEC maps is also observed, especially over the edge areas with an improvement of 60–80%. Overall, the proposed SHAKING method exhibits notable advantage over the existing regional VTEC modeling techniques, which can be used for regional TEC modeling and associated high-precision positioning applications.     

基于SPH方法的球形贮箱液体晃动动力学分析

针对发生液体破碎现象的航天器球形贮箱液体晃动问题,基于光滑粒子流体动力学方法（SPH方法）对N S方程进行离散近似,分别对核函数、应力张量、边界问题、邻域粒子搜索等方面进行选择确定或计算处理,最终得到基于SPH方法的流体动力学模型。通过建立包括不同燃料液体体积及不同挡板安装位置的仿真模型,分析了航天器球形贮箱的受力及局部压强。通过分析贮箱受力情况及设置测量点压强情况,得到运动参数对晃动的影响规律,为航天器球形贮箱的设计提供了理论依据。     

航天器表面污染物质沉积变化和控制因子评估

航天器进入空间环境以后,空间环境分子污染和颗粒污染形成了航天器表面污染层,从而对航天器的各技术分系统产生不同的负面影响.介绍了中外中轨道航天器表面污染物质沉积变化在轨探测结果,并对污染物质沉积量变化和控制因子做初步评估.结果表明,污染物质沉积量在航天器入轨初期的1~2年内受航天器自身出气物质量、放气速率、表面温度及所处的气流方向等因子所控制.初期沉积量大,正是受到航天器入轨后自身出气量大、放气速率较高的控制,同时迎风面比背风面沉积量大.入轨后期表面沉积量长期变化呈现出明显的降变或缓慢涨落,而且具有全向性特征,因此探讨了具有全向性影响能力的控制因子相关特性,其中高能粒子通量和太阳紫外辐射通量变化可能是主要控制因子.      

A model of radiation conditions during spacecraft flights in the interplanetary space and in the Earth's magnetosphere.

Based on the available measurement data, simulations of radiation conditions during spacecraft flights in the interplanetary space and in the Earth's and Jupiter's radiation belts has been carried out. The > or = 10 MeV and > or = 30 MeV solar flare proton fluence forecast has been proposed for Cycle 22. Radiation conditions due to both magnetospheric electrons and protons and to solar flare protons, magnetic rigidity cutoff being taken into account, have been evaluated on spacecraft trajectories in the Earth's and Jupiter's magnetospheres.     

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

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

Radiation analysis for manned missions to the Jupiter system.

An analysis for manned missions targeted to the Jovian system has been performed in the framework of the NASA RASC (Revolutionary Aerospace Systems Concepts) program on Human Exploration beyond Mars. The missions were targeted to the Jupiter satellite Callisto. The mission analysis has been divided into three main phases, namely the interplanetary cruise, the Jupiter orbital insertion, and the surface landing and exploration phases. The interplanetary phase is based on departure from the Earth-Moon L1 point. Interplanetary trajectories based on the use of different propulsion systems have been considered, with resulting overall cruise phase duration varying between two and five years. The Jupiter-approach and the orbital insertion trajectories are considered in detail, with the spacecraft crossing the Jupiter radiation belts and staying around the landing target. In the surface exploration phase the stay on the Callisto surface is considered. The satellite surface composition has been modeled based on the most recent results from the GALILEO spacecraft. In the transport computations the surface backscattering has been duly taken into account. Particle transport has been performed with the HZETRN heavy ion code for hadrons and with an in-house developed transport code for electrons and bremsstrahlung photons. The obtained doses have been compared to dose exposure limits.     

基于RQPSO-DMPC的多无人机编队自主重构控制方法

针对敌方防御区域内各种威胁,为了实现隐蔽突防并实施对敌有效打击,在突防过程中多无人机(UAV)编队需要进行重构控制,并且编队内的相互避碰问题与通信约束问题也需考虑。通过建立无人机虚拟领航编队模型并引入邻居集,采用分布式模型预测控制(DMPC)同时构建多无人机编队的重构代价函数,提出采用改进量子粒子群优化(RQPSO)算法进行求解,并将求解结果与采用粒子群优化算法的结果进行对比。仿真结果表明,本文算法能够有效控制多无人机编队完成自主重构,实现安全隐蔽突防任务。     

Survey of current situation in radiation belt modeling.

The study of Earth's radiation belts is one of the oldest subjects in space physics. Despite the tremendous progress made in the last four decades, we still lack a complete understanding of the radiation belts in terms of their configurations, dynamics, and detailed physical accounts of their sources and sinks. The static nature of early empirical trapped radiation models, for examples, the NASA AP-8 and AE-8 models, renders those models inappropriate for predicting short-term radiation belt behaviors associated with geomagnetic storms and substorms. Due to incomplete data coverage, these models are also inaccurate at low altitudes (e.g., <1000 km) where many robotic and human space flights occur. The availability of radiation data from modern space missions and advancement in physical modeling and data management techniques have now allowed the development of new empirical and physical radiation belt models. In this paper, we will review the status of modern radiation belt modeling.     

The bimodal magnetosphere and ring current, radiation belt, and tail transducers

It has generally been assumed that a geomagnetic storm is entirely driven by external forces—e.g., solar wind E_y = V_x × B_z, V_x, V²_x (where the components of the electric field, E, the magnetic field, B, and velocity, V, are given in GSE coordinates)—which would imply that particle injections in the ring current (RC) or outer radiation belts should be highly correlated. However the data from ISTP are showing that the magnetosphere can have at least two very different responses to the same solar wind (SW) conditions: a classic, enhanced RC with Dst response, or a 1000-fold increase in the outer radiation belt MeV electrons (ORBE). August 29, October 14 and 23, 1996 are examples of Dst storms, whereas April 15, 1996 and January 10, 1997 are examples of MeV storms. It is this second response that is so deadly to some geosynchronous spacecraft, whereas geomagnetic storms are categorized by the first response. Neither of these appear to be correlated to the SW conditions driving substorms. Why should the SW energy appear in the radiation belts or the ring current independently? We hypothesize that the RC couples to the electric power available (E_y), the ORBE couple to the mechanical power available (V_x), and the Tail couples to the magnetic energy (B_z) available in the SW. The transducer for RC may be subauroral parallel potentials, the transducer for ORBE may be the cusp, while the Tail substorm transducer is yet a third independent mechanism for extracting SW energy. Evidence for this theory comes from the novel POLAR satellite that traverses the cusp, the plasmasheet and the radiation belts.