On the validity of the aluminum equivalent approximation in space radiation shielding applications

The origin of the aluminum equivalent shield approximation in space radiation analysis can be traced back to its roots in the early years of the NASA space programs (Mercury, Gemini and Apollo) wherein the primary radiobiological concern was the intense sources of ionizing radiation causing short term effects which was thought to jeopardize the safety of the crew and hence the mission. Herein, it is shown that the aluminum equivalent shield approximation, although reasonably well suited for that time period and to the application for which it was developed, is of questionable usefulness to the radiobiological concerns of routine space operations of the 21st century which will include long stays onboard the International Space Station (ISS) and perhaps the moon. This is especially true for a risk based protection system, as appears imminent for deep space exploration where the long-term effects of Galactic Cosmic Ray (GCR) exposure is of primary concern. The present analysis demonstrates that sufficiently large errors in the interior particle environment of a spacecraft result from the use of the aluminum equivalent approximation, and such approximations should be avoided in future astronaut risk estimates. In this study, the aluminum equivalent approximation is evaluated as a means for estimating the particle environment within a spacecraft structure induced by the GCR radiation field. For comparison, the two extremes of the GCR environment, the 1977 solar minimum and the 2001 solar maximum, are considered. These environments are coupled to the Langley Research Center (LaRC) deterministic ionized particle transport code High charge (Z) and Energy TRaNsport (HZETRN), which propagates the GCR spectra for elements with charges (Z) in the range 1 ? Z ? 28 (H–Ni) and secondary neutrons through selected target materials. The coupling of the GCR extremes to HZETRN allows for the examination of the induced environment within the interior of an idealized spacecraft as approximated by a spherical shell shield, and the effects of the aluminum equivalent approximation for a good polymeric shield material such as generic polyethylene (PE). The shield thickness is represented by a 25 g/cm² spherical shell. Although, one could imagine the progression to greater thickness, the current range will be sufficient to evaluate the qualitative usefulness of the aluminum equivalent approximation. Upon establishing the inaccuracies of the aluminum equivalent approximation through numerical simulations of the GCR radiation field attenuation for PE and aluminum equivalent PE spherical shells, we further present results for a limited set of commercially available, hydrogen rich, multifunctional polymeric constituents to assess the effect of the aluminum equivalent approximation on their radiation attenuation response as compared to the generic PE.     

A model for predicting the radiation exposure for mission planning aboard the international space station

The International Space Station Cosmic Radiation Exposure Model (ISSCREM) has been developed as a possible tool for use in radiation mission planning as based on operational data collected with a tissue equivalent proportional counter (TEPC) aboard the ISS since 2000. It is able to reproduce the observed trapped radiation and galactic cosmic radiation (GCR) contributions to the total dose equivalent to within ±20% and ±10%, respectively, as would be measured by the onboard TEPC at the Zvezda Service Module panel 327 (SM-327). Furthermore, when these contributions are combined, the total dose equivalent that would be measured at this location is estimated to within ±10%. The models incorporated into ISSCREM correlate the GCR dose equivalent rate to the cutoff rigidity magnetic shielding parameter and the trapped radiation dose equivalent rate to atmospheric density inside the South Atlantic Anomaly. The GCR dose equivalent rate is found to vary minimally with altitude and TEPC module location however, due to the statistics and data available, the trapped radiation model could only be developed for the TEPC located at SM-327. Evidence of the variation in trapped radiation dose with detector orientation and the East–West asymmetry were observed at this location.     

Method for the prediction of the effective dose equivalent to the crew of the International Space Station

This paper describes a methodology for assessing the pre-mission exposure of space crew aboard the International Space Station (ISS) in terms of an effective dose equivalent. In this approach, the PHITS Monte Carlo code was used to assess the particle transport of galactic cosmic radiation (GCR) and trapped radiation for solar maximum and minimum conditions through an aluminum shield thickness. From these predicted spectra, and using fluence-to-dose conversion factors, a scaling ratio of the effective dose equivalent rate to the ICRU ambient dose equivalent rate at a 10 mm depth was determined. Only contributions from secondary neutrons, protons, and alpha particles were considered in this analysis.     

An analytical model for the prediction of a micro-dosimeter response function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (δ ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCRs) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 s/ion. The ionizing radiation environment at LEO is represented by O’Neill’s GCR model (2004), covering charged particles in the 1 ? Z ? 28 range. O’Neill’s free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.     

基于ASIC芯片的多探测单元设计与验证

以空间辐射环境中的粒子为研究目标, 研制了用硅微条传感器作为探头、用集成芯片IDE3160进行信号处理的空间粒子探测系统. 该系统采用两片硅微条传感器组成的硅微条探测阵列作为前端探头, 并应用数字化的信号处理方法, 获取空间粒子入射的位置及在硅微条内单位长度沉积的能量(线性能量传递, LET). 从诱发单粒子效应的物理机制角度, 对比重离子和脉冲激光在硅半导体中所产生物理效应的不同. 采用1.064 μm脉冲激光开展系统测试, 获得良好的LET线性结果: 该系统数据采集所需时间为2.47 ms, 可探测到的LET阈值约为0.1 MeV·cm²·mg^–1, 皮尔逊相关系数(PCC)达0.998, 表明系统测量结果与理论设计符合性良好. 该系统动态范围宽、线性度好, 具有较高的集成度、可拓展性及可移植性, 可以搭载在各种空间探测卫星上.     

利用质子强度变化对太阳质子事件进行短期预报方法初步研究

采用GOSE-10卫星4～9 MeV(P2),9～15 MeV(P3),15～40 MeV(P4),40～80 MeV(P5)能段上的质子通量数据,结合质子能谱,对太阳质子事件发生前各能谱参数的变化特征进行分析,详细介绍利用能谱参数的变化特征及能量E>10 MeV的质子通量数据对太阳质子事件进行预报的新方法,并运用这种方法对2002-2006年期间太阳质子事件进行了预报.预报结果显示,预报提前量最多达到100 h以上,对质子事件的报准率达97.5%,预报方法具备一定的有效性和实用性.     

An improved neutron transport algorithm for HZETRN

Long-term human presence in space requires the inclusion of radiation constraints in mission planning and the design of shielding materials, structures and vehicles. It is necessary to expose the numerical tools commonly used in radiation analyses to extensive verification, validation and uncertainty quantification. In this paper, the numerical error associated with energy discretization in HZETRN is addressed. An inadequate numerical integration scheme in the transport algorithm is shown to produce large errors in the low energy portion of the neutron and light ion fluence spectra. It is further shown that the errors result from the narrow energy domain of the neutron elastic cross section spectral distributions and that an extremely fine energy grid is required to resolve the problem under the current formulation. Since adding a sufficient number of energy points will render the code computationally inefficient, we revisit the light ion and neutron transport theory developed for HZETRN and focus on neutron elastic interactions. Two numerical methods (average value and collocation) are developed to provide adequate resolution in the energy domain and more accurately resolve the neutron elastic interactions. An energy grid convergence study is conducted to demonstrate the improved stability of the new methods. Based on the results of the convergence study and the ease of implementation, the average value method with a 100 point energy grid is found to be suitable for future use in HZETRN.     

Pion and electromagnetic contribution to dose: Comparisons of HZETRN to Monte Carlo results and ISS data

Recent work has indicated that pion production and the associated electromagnetic (EM) cascade may be an important contribution to the total astronaut exposure in space. Recent extensions to the deterministic space radiation transport code, HZETRN, allow the production and transport of pions, muons, electrons, positrons, and photons. In this paper, the extended code is compared to the Monte Carlo codes, Geant4, PHITS, and FLUKA, in slab geometries exposed to galactic cosmic ray (GCR) boundary conditions. While improvements in the HZETRN transport formalism for the new particles are needed, it is shown that reasonable agreement on dose is found at larger shielding thicknesses commonly found on the International Space Station (ISS). Finally, the extended code is compared to ISS data on a minute-by-minute basis over a seven day period in 2001. The impact of pion/EM production on exposure estimates and validation results is clearly shown. The Badhwar–O’Neill (BO) 2004 and 2010 models are used to generate the GCR boundary condition at each time-step allowing the impact of environmental model improvements on validation results to be quantified as well. It is found that the updated BO2010 model noticeably reduces overall exposure estimates from the BO2004 model, and the additional production mechanisms in HZETRN provide some compensation. It is shown that the overestimates provided by the BO2004 GCR model in previous validation studies led to deflated uncertainty estimates for environmental, physics, and transport models, and allowed an important physical interaction (π/EM) to be overlooked in model development. Despite the additional π/EM production mechanisms in HZETRN, a systematic under-prediction of total dose is observed in comparison to Monte Carlo results and measured data.     

The ORGANIC experiment on EXPOSE-R on the ISS: Flight sample preparation and ground control spectroscopy

In March of 2009, the ORGANIC experiment integrated into the European multi-user facility EXPOSE-R, containing experiments dedicated to Astrobiology, was mounted through Extra Vehicular Activity (EVA) externally on the International Space Station (ISS). The experiment exposed organic samples of astronomical interest for a duration of 97 weeks (∼22 months) to the space environment. The samples that were returned to Earth in spring 2011, received a total UV radiation dose during their exposure including direct solar irradiation of >2500 h, exceeding the limits of laboratory simulations. We report flight sample preparation and pre-flight ultraviolet–visible (UV–Vis) characterization of the ORGANIC samples, which include 11 polycyclic aromatic hydrocarbons (PAHs) and three fullerenes. The corresponding time-dependent ground control monitoring experiments for ORGANIC measured over ∼19 months are presented and the results anticipated upon return of the samples are discussed. We present the first UV–Vis spectrum of solid circobiphenyl (C₃₈H₁₆). Further, we present the first published UV–Vis spectra of diphenanthro[9,10-b′,10′-d]thiophene (C₂₈H₁₆S), dinaphtho[8,1,2-abc,2′,1′,8′-klm]coronene (C₃₆H₁₆), tetrabenzo[de,no,st,c′d′]heptacene (C₄₂H₂₂), and dibenzo[jk,a′b′]octacene (C₄₀H₂₂) in solid phase and in solution. The results of the ORGANIC experiment are expected to enhance our knowledge of the evolution and degradation of large carbon-containing molecules in space environments.     

基于多冲量能耗估算的小推力任务窗口搜索

针对小推力轨道设计中任务窗口参数具有范围跨度大、非连续性强、对设计结果影响显著等特点,提出了一种基于多冲量能耗估算的任务窗口搜索方法.该方法基于多冲量轨道和小推力轨道的能耗一致性,先由朗伯特两冲量轨道经过数次迭代扩展快速获得多冲量轨道,再通过多冲量轨道能耗估算对任务区间进行全局栅格搜索,最终得到任务窗口的可行区域.仿真结果表明:1)该方法同时适用于交会和飞越任务;2)同朗伯特两冲量相比,它与小推力轨道能耗的一致性更强;3)同Sims和Flanagan最优多冲量相比,计算时间缩短了1～2个量级,计算效率显著提高.地球—火星交会任务的仿真应用表明,该方法只需8.59s即可获得任务窗口的可行区域,进而快速给出小推力最优转移轨道,验证了该方法的有效性.     

Preliminary results of space experiment: Implications for the effects of space radiation and microgravity on survival and mutation induction in human cells

In view of the concern for the health of astronauts that may one day journey to Mars or the Moon, we investigated the effect that space radiation and microgravity might have on DNA damage and repair. We sent frozen human lymphoblastoid TK6 cells to the International Space Station where they were maintained under frozen conditions during a 134-day mission (14 November 2008 to 28 March 2009) except for an incubation period of 8 days under 1G or μG conditions in a CO₂ incubator. The incubation period started after 100 days during which the cells had been exposed to 54 mSv of space radiation. The incubated cells were then refrozen, returned to Earth, and compared to ground control samples for the determination of the influence of microgravity on cell survival and mutation induction. The results for both varied from experiment to experiment, yielding a large SD, but the μG sample results differed significantly from the 1G sample results for each of 2 experiments, with the mean ratio of μG to 1G being 0.55 for the concentration of viable cells and 0.59 for the fraction of thymidine kinase deficient (TK⁻) mutants. Among the mutants, non-loss of zygosity events (point mutations) were less frequent (31%) after μG incubation than after 1G incubation, which might be explained by the influence of μG on cellular metabolic or physiological function. Additional experiments are needed to clarify the effect of μG interferes on DNA repair.     

动磁式线性电机理论研究与测试技术

本文以昆明物理研究所的SC100H线性斯特林制冷机为研究对象,对该制冷机采用的动磁式线性振荡电机磁路特性进行了理论与实验研究。建立了电机理论模型,对模型进行了结构简化以便于分析计算,利用等效磁路法和机电能量转换原理,分析电机推力特性与动子相对位置等的关系,利用拉力测试系统对电机在不同输入电流的推力特性进行了测试,测试结果与理论值很好的吻合,从而为电机结构参数和运行参数的设计优化提供了设计参考与依据。     

Statistical validation of HZETRN as a function of vertical cutoff rigidity using ISS measurements

Measurements taken in Low Earth Orbit (LEO) onboard the International Space Station (ISS) and transit vehicles have been extensively used to validate radiation transport models. Primarily, such comparisons were done by integrating measured data over mission or trajectory segments so that individual comparisons to model results could be made. This approach has yielded considerable information but is limited in its ability to rigorously quantify and differentiate specific model errors or uncertainties. Further, as exploration moves beyond LEO and measured data become sparse, the uncertainty estimates derived from these validation cases will no longer be applicable. Recent improvements in the underlying numerical methods used in HZETRN have resulted in significant decreases in code run time. Therefore, the large number of comparisons required to express error as a function of a physical quantity, like cutoff rigidity, are now possible. Validation can be looked at in detail over any portion of a flight trajectory (e.g. minute by minute) such that a statistically significant number of comparisons can be made. This more rigorous approach to code validation will allow the errors caused by uncertainties in the geometry models, environmental models, and nuclear physics models to be differentiated and quantified. It will also give much better guidance for future model development. More importantly, it will allow a quantitative means of extrapolating uncertainties in LEO to free space. In this work, measured data taken onboard the ISS during solar maximum are compared to results obtained with the particle transport code HZETRN. Comparisons are made at a large number (∼77,000) of discrete time intervals, allowing error estimates to be given as a function of cutoff rigidity. It is shown that HZETRN systematically underestimates exposure quantities at high cutoff rigidity. The errors are likely associated with increased angular variation in the geomagnetic field near the equator, the lack of pion production in HZETRN, and errors in high energy nuclear physics models, and will be the focus of future work.     

Interplanetary proton cumulated fluence model update

Solar particle events leading to important increase of particle fluxes at energies of order of magnitude ranging from MeV to GeV constitute an important hazard for space missions. They may lead to effects seen in microelectronics or damage to solar cells and constitute a potential hazard for manned missions. Cumulative damage is commonly expressed as a function of fluence which is defined as the integral of the flux over time. A priori deterministic estimates of the expected fluence cannot be made because over the time scale of a space mission, the fluence can be dominated by the contribution of a few rare and unpredictable high intensity events. Therefore, statistical approaches are required in order to estimate fluences likely to be encountered by a space mission in advance. This paper extends work done by Rosenqvist et al. [Rosenqvist, L., Hilgers, A., Evans, H., Daly, E., Hapgood, M., Stamper, R., Zwickl, R., Bourdarie, S., Boscher, D. Toolkit for updating interplanetary proton-cumulated fluence models. J. Spacecraft Rockets, 42(6), 1077–1090, 2005] to describe an updated predictive engineering model for the proton interplanetary fluence with energies >30 MeV. This model is derived from a complete list of solar proton fluences based on data from a number of calibrated sources covering almost three solar cycles.     

Analytical-HZETRN model for rapid assessment of active magnetic radiation shielding

The use of active radiation shielding designs has the potential to reduce the radiation exposure received by astronauts on deep-space missions at a significantly lower mass penalty than designs utilizing only passive shielding. Unfortunately, the determination of the radiation exposure inside these shielded environments often involves lengthy and computationally intensive Monte Carlo analysis. In order to evaluate the large trade space of design parameters associated with a magnetic radiation shield design, an analytical model was developed for the determination of flux inside a solenoid magnetic field due to the Galactic Cosmic Radiation (GCR) radiation environment. This analytical model was then coupled with NASA’s radiation transport code, HZETRN, to account for the effects of passive/structural shielding mass. The resulting model can rapidly obtain results for a given configuration and can therefore be used to analyze an entire trade space of potential variables in less time than is required for even a single Monte Carlo run. Analyzing this trade space for a solenoid magnetic shield design indicates that active shield bending powers greater than ∼15 Tm and passive/structural shielding thicknesses greater than 40 g/cm² have a limited impact on reducing dose equivalent values. Also, it is shown that higher magnetic field strengths are more effective than thicker magnetic fields at reducing dose equivalent.