屏蔽材料对木星轨道卫星内部介质充电效应的影响研究

在木星轨道的空间辐射环境中,占主导地位的粒子是能量大于1MeV（甚至高于100MeV）的高能电子,这可能会产生卫星内部介质充电效应。在卫星的防辐射设计中,通常需要一定厚度的材料来屏蔽这些电子,使得进入卫星内部的电子通量达到安全的水平。利用所建立的GEANT4-RIC（radiation induced conductivity）方法,研究了运行于木星轨道的卫星对高能电子的最佳屏蔽材料设计。研究了铝、钛、铁、铜、钽和铅作为卫星屏蔽材料的可能性。研究结果表明,在木星探测任务中,为了减轻卫星内部介质充电效应,高原子序数材料比低原子序数材料在相同质量下提供的屏蔽效果更好。因此,用钛、铁、铜、钽或铅代替地球轨道卫星上常用的屏蔽材料铝,可以节省屏蔽质量。     

基于数据网格化方法的低轨辐射带建模技术

地球辐射带中的高能带电粒子是引起航天器材料和器件性能退化甚至失效的主要空间环境因素.因此,航天器设计中所采用的辐射带模型的准确程度对于航天器的生存能力和航天任务的完成质量至关重要.在利用我国自主辐射带高能粒子探测数据进行的辐射带建模中,探测数据的空间网格化是一项非常重要的工作.介绍了我国辐射带探测数据的情况,以及辐射带建模的方法和步骤;重点研究了不同插值方法在低地球轨道(LEO)空间辐射带建模数据网格化中的应用,并开展了误差分析.研究结果表明:在各种常用的插值方法中,反距离加权法、自然邻点法和最近邻点法适合工程化应用.其中,反距离加权法生成的数据网格对粒子通量的反演结果精度最高,该方法采用低阶距离时得到的反演结果更为合理.      

大气中子实时计算模式

为实时评估0~100km高度范围内的大气中子全球分布,对宇宙线在地磁场和大气中的传输过程进行了分析.利用蒙特卡罗方法工具包Geant4,预先计算不同能量的粒子在大气层中产生的次级粒子能谱分布,形成大气次级粒子数据库,并与相关模型进行对比,验证了该数据库的有效性和可靠性.以实测或预报的空间环境参数作为输入,计算同步轨道银河宇宙线和太阳质子事件能谱以及100km高度上的地磁垂直截止刚度,最终得到大气层顶上的粒子能谱.通过对大气次级粒子数据库的线性插值,实现1h分辨率的大气中子能谱和辐射剂量全球分布的实时计算.      

航天员受银河宇宙线辐射的剂量计算

在近地空间(LEO)和深空探测中,航天员遭受的辐射风险主要来自于银河宇宙线(GCR)照射.银河宇宙线的辐射剂量是航天员辐射风险评价的基础.国际放射防护委员会(ICRP)于2013年提出了新的航天员空间辐射剂量估算方法,以更准确给出空间重离子辐射的剂量.基于此方法,开发了宇宙线粒子在物质中输运的蒙特卡罗程序,并在程序中实现用中国成年男性人体数字模型来仿真航天员.采用该程序计算了粒子(Z=1~92)各向同性照射航天员时器官的通量-器官剂量转换因数,并估算出航天员在近地轨道空间受银河宇宙线辐射的剂量.     

脉冲激光模拟空间载荷单粒子效应研究进展

脉冲激光作为模拟测试空间探测载荷半导体器件的单粒子效应现象的一种较新型手段,具有可以定位器件对单粒子效应敏感的具体单元以及动态测试电路系统对单粒子效应的时间响应特性的特点,能够满足工程部门、器件研发部门的不同需求。通过实验与理论研究,建立单粒子锁定与翻转效应的激光阈值能量与重离子LET值的对应关系,解决了脉冲激光模拟测试的激光结果如何定量的关键问题,据此可以定量摸底评估器件的单粒子效应敏感度,使脉冲激光测试载荷的结果更具评价以及指导意义,这对建立统一的脉冲激光单粒子效应评估试验标准以及对脉冲激光试验的推广具有重要意义。空间探测载荷发生单粒子效应后器件功能特性及电路系统的影响、防范单粒子效应电路条件影响的手段下电路系统的抗单粒子效应设计措施是的有效性,以及为空间探测专门研制的抗辐射ASIC电路评价,都需要更加精细的单粒子效应测试方法。通过建立便捷、低成本的脉冲激光定量试验的手段,解决了空间探测载荷上述单粒子效应试验的问题。     

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

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

天基微小空间碎片探测研究

随着空间碎片数量的不断增多,天基微小空间碎片探测已经成为一个热点．首先介绍了空间碎片在低地球轨道上的分布情况以及它对于航天器的危害,然后介绍了国外微小碎片探测器的基本情况,并在这些探测结果的基础上提出了一个探测器方案．这种探测器的传感器采用了新型的压电材料聚偏二氟乙稀(PVDF),使用了飞行时间法(TOF)准确测定空间碎片的飞行速度,以及快脉冲分析系统分析碎片的质量．     

临近空间中子辐射环境分析及其引起的单粒子效应预计研究

使用Space Radiation 7.0工具分析临近空间中子辐射环境，研究其与海拔高度、太阳活动和经纬度的关系及内在原因.在此基础上，提出了一种基于蒙特卡罗方法的大气中子实时错误率预计方法，并以航空电子系统关键集成电路FPGA为例，预计其单粒子翻转敏感模块包括配置存储器、块存储器和用户触发器，单粒子功能中断敏感模块包括上电复位电路、SelectMAP接口等的实时飞行错误率.结果表明，配置存储器中发生的单粒子翻转达到总单粒子翻转率的77%，而上电复位电路和SelectMAP接口中发生的单粒子功能中断各占总单粒子功能中断率的36%.根据RTCA DO-254对飞行系统失效等级的划分，该FPGA器件不可用于航空电子系统关键位置.      

电离总剂量复合屏蔽模拟仿真及验证试验

空间环境中辐射粒子的电离总剂量效应对卫星电子器件危害严重,需要采用合适的材料进行屏蔽防护.本文采用蒙特卡罗方法模拟材料对电子的屏蔽,将双层复合屏蔽方法与单质屏蔽方法进行对比,结果表明对电子而言,复合屏蔽在屏蔽厚度足够大时比单质屏蔽效果更好.利用90Sr-90Y电子放射源进行了复合屏蔽效果的验证试验,试验结果与模拟结果规律相符,研究结果可为辐射防护的优化设计提供参考.      

2 AU以内的“渐进型”太阳高能粒子事件模拟

太阳高能粒子（Solar Energetic Particle,SEP）事件是影响地球空间以及深空辐射环境的主要因素之一。“渐进型”太阳高能粒子事件中的高能粒子主要来自于日冕物质抛射（Coronal Mass Ejection,CME）所驱动的激波扩散加速（Diffusive Shock Acceleration,DSA）过程。CME驱动的激波在行星际的传播过程中,其结构不断演化,进而影响到高能粒子的加速过程。本文利用二维太阳高能粒子加速和传播模型,对发生于2014年4月18日的太阳高能粒子事件实例进行了数值模拟。模型考察了黄道面上2 AU的距离以内包含地球所在位置的4个不同点,分别计算了每个点上高能粒子的通量。数值模拟的结果表明：黄道面内不同位置的观察点,与激波波前的磁力线连接不同,从而导致观察点处高能粒子的通量有着显著的差异。该模型的计算结果可以为深空探测计划开展辐射环境研究提供必要的输入。     

Monte Carlo calibration of the response of the University of Chicago’s Cosmic Ray Nuclei Experiment (CRNE) on IMP-8 to electrons above 0.5 MeV

Modern instrument-simulation techniques offer the possibility of increasing the scientific yield from archival space datasets. In this paper, we report on a simulation of the electron response of the University of Chicago’s Cosmic Ray Nuclei Experiment (CRNE) instrument on the IMP-8 satellite. IMP-8/CRNE returned data from 1973 to 2006. The CRNE particle telescope was designed to measure the isotopic composition of Galactic cosmic-ray (GCR) nuclei and has also been used in many studies of protons and ions above 10 MeV/nucleon from solar energetic particle (SEP) events. But CRNE also functions as a highly-capable detector for solar electrons above 0.5 MeV, an energy range that has not been extensively studied. Utilization of the CRNE electron data has heretofore been limited by the fact that CRNE was never calibrated for electrons. We have therefore used the GEANT4 Monte Carlo simulation package to model the CRNE response to electrons and (separately) protons for multiple energies and incident angles. The results were used to compute the energy- and angle-dependence of the effective area and the energy-dependence of the geometric factor. The response to protons, which was already well understood, was used to verify the mass model, the simulation settings, and the post-processing software. Our simulation of the IMP-8/CRNE electron response now allows analysis of hundreds of relativistic solar electron events observed by CRNE over the years, including studies of evolution of electron energy spectra with high time resolution. We show examples of these results and briefly discuss potential applications to future scientific investigations.     

Simulation of ASTROD I test mass charging due to solar energetic particles and interplanetary electrons

As ASTROD I travels through space, its test mass will accrue charge due to exposure of the spacecraft to high-energy particles. This test mass charge will result in Coulomb forces between the test mass and the surrounding electrodes. In earlier work, we have used the GEANT 4 toolkit to simulate charging of the ASTROD test mass due to cosmic-ray protons of energies between 0.1 and 1000 GeV at solar maximum and at solar minimum. Here we use GEANT 4 to simulate the charging process due to solar energetic particle events and interplanetary electrons. We then estimate the test mass acceleration noise due to these fluxes. The predicted charging rates range from 2247 e⁺/s to 47,055 e⁺/s, at peak intensity, for the four largest SEP events in September and October 1989. Although the noise due to charging exceeds the ASTROD I budget for the two larger events, it can be suppressed through continuous discharging. The acceleration noise during the two small events is well below the design target. The charging rate of the ASTROD I test mass due to interplanetary electrons in this simulation is about −11% of the cosmic-ray protons at solar minimum, and over −37% at solar maximum. In addition to the Monte Carlo uncertainty, an error of ±30% in the net charging rates should be added to account for uncertainties in the spectra, physics models and geometry implementations.     

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

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

基于SIMION软件的空间等离子体探测器的数值仿真

空间等离子体探测器采用带偏转板的柱形静电分析器作为传感器,具有大于320°探测视场,可探测约10eV～30keV的热等离子体。利用SIMION软件采用均匀随机抽样建立粒子源的方法仿真探测器的基本特性参数:静电分析器因子、能量分辨率、偏转板因子和偏转角。仿真结果与实测结果吻合很好。静电分析器因子的仿真与实测结果偏差为3.2%,能量分辨率仿真与实测结果偏差为1.5%,偏转角-偏转板因子的关系式与实测结果的相关系数为0.999 4。仿真结果可以有效地应用于该类仪器的设计和定标过程中。     

离子推力器羽流热效应仿真分析

离子推力器工作时向外喷出的羽流与航天器表面碰撞,会引起敏感材料热变形等热效应,严重时会导致航天任务失败。针对兰州空间技术物理研究所研制的LIPS-200型离子推力器羽流热效应进行了仿真分析。仿真中,使用粒子网格（PIC）方法处理等离子体运动,使用直接模拟蒙特卡罗（DSMC）方法处理粒子间碰撞,使用Maxwell模型处理粒子与壁面的能量交换,对电推进羽流热效应测量中的部分测点进行了数值模拟。结果表明,仿真结果与实验数据符合较好,离子推力器出口轴线上滞止热流仿真值与实验测量值误差小于17.0%。此外,热流计对流场的影响主要集中在热流计附近0.1 m范围内,对整体流场影响较小。      

Comparisons of several transport models in their predictions in typical space radiation environments

We have used several transport codes to calculate dose and dose equivalent values as well as the particle spectra behind a slab or inside a spherical shell shielding in typical space radiation environments. Two deterministic codes, HZETRN and UPROP, and two Monte Carlo codes, FLUKA and Geant4, are included. A soft solar particle event, a hard solar particle event, and a solar minimum galactic cosmic rays environment are considered; and the shielding material is either aluminum or polyethylene. We find that the dose values and particle spectra from HZETRN are in general rather consistent with Geant4 except for neutrons. The dose equivalent values from HZETRN and Geant4 are not far from each other, but the HZETRN values behind shielding are often lower than the Geant4 values. Results from FLUKA and Geant4 are mostly consistent for considered cases. However, results from the legacy code UPROP are often quite different from the other transport codes, partly due to its non-consideration of neutrons. Comparisons for the spherical shell geometry exhibit the same qualitative features as for the slab geometry. In addition, results from both deterministic and Monte Carlo transport codes show that the dose equivalent inside the spherical shell decreases from the center to the inner surface and this decrease is large for solar particle events; consistent with an earlier study based on deterministic radiation transport results. This study demonstrates both the consistency and inconsistency among these transport models in their typical space radiation predictions; further studies will be required to pinpoint the exact physics modules in these models that cause the differences and thus may be improved.     

A detailed FLUKA-2005 Monte-Carlo simulation for the ATIC detector

We have performed a detailed Monte-Carlo (MC) simulation for the Advanced Thin Ionization Calorimeter (ATIC) detector using the MC code FLUKA-2005 which is capable of simulating particles up to 10 PeV. The ATIC detector has completed two successful balloon flights from McMurdo, Antarctica lasting a total of more than 35 days. ATIC is designed as a multiple, long duration balloon flight, investigation of the cosmic ray spectra from below 50 GeV to near 100 TeV total energy; using a fully active Bismuth Germanate (BGO) calorimeter. It is equipped with a large mosaic of silicon detector pixels capable of charge identification, and, for particle tracking, three projective layers of x–y scintillator hodoscopes, located above, in the middle and below a 0.75 nuclear interaction length graphite target. Our simulations are part of an analysis package of both nuclear (A) and energy dependences for different nuclei interacting in the ATIC detector. The MC simulates the response of different components of the detector such as the Si-matrix, the scintillator hodoscopes and the BGO calorimeter to various nuclei. We present comparisons of the FLUKA-2005 MC calculations with GEANT calculations and with the ATIC CERN data.     

旋转调制器分析方法的研究

本文研究了为探测150keV左右的硬X射线源而设计的旋转调制器的数据分析方法。改进Durouchoux的方法使之适用于无位置灵敏的镶嵌式旋转调制器系统,并在计算中采用了我们自己导出的响应函数,为验证其正确性,我们做了Monte Carlo模拟加以校验。