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.     

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.     

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.     

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.     

Numerical study of the generation of linear energy transfer spectra for space radiation applications

In analyzing charged particle spectra in space due to galactic cosmic rays (GCR) and solar particle events (SPE), the conversion of particle energy spectra into linear energy transfer (LET) distributions is a convenient guide in assessing biologically significant components of these spectra. The mapping of LET to energy is triple valued and can be defined only on open energy subintervals where the derivative of LET with respect to energy is not zero. Presented here is a well-defined numerical procedure which allows for the generation of LET spectra on the open energy subintervals where, in spite of their singular nature, the spectra are integrable. The efficiency of the numerical procedures is demonstrated by providing examples of computed differential and integral LET spectra and their equilibrium components for historically large SPEs and 1977 solar minimum GCR environments. Due to the biological significance of tissue, all simulations are done with tissue as the target material.