Francis F. Badavi

A semiempirical nuclear fragmentation model

An abrasion/ablation model of heavy ion fragmentation is derived which includes a second order correction for the surface energy term and provides a reasonable representation of the present elemental fragmentation cross sections. The full development of the model must await the resolution of disagreement among different experiments and an expansion of the experimental data base to a broader set of projectile-target combinations.

Solar proton exposure of an ICRU sphere within a complex structure Part I: Combinatorial geometry.

The 3DHZETRN code, with improved neutron and light ion (Z≤2) transport procedures, was recently developed and compared to Monte Carlo (MC) simulations using simplified spherical geometries. It was shown that 3DHZETRN agrees with the MC codes to the extent they agree with each other. In the present report, the 3DHZETRN code is extended to enable analysis in general combinatorial geometry. A more complex shielding structure with internal parts surrounding a tissue sphere is considered and compared against MC simulations. It is shown that even in the more complex geometry, 3DHZETRN agrees well with the MC codes and maintains a high degree of computational efficiency.

Galactic cosmic ray simulation at the NASA Space Radiation Laboratory

Most accelerator-based space radiation experiments have been performed with single ion beams at fixed energies. However, the space radiation environment consists of a wide variety of ion species with a continuous range of energies. Due to recent developments in beam switching technology implemented at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory (BNL), it is now possible to rapidly switch ion species and energies, allowing for the possibility to more realistically simulate the actual radiation environment found in space. The present paper discusses a variety of issues related to implementation of galactic cosmic ray (GCR) simulation at NSRL, especially for experiments in radiobiology. Advantages and disadvantages of different approaches to developing a GCR simulator are presented. In addition, issues common to both GCR simulation and single beam experiments are compared to issues unique to GCR simulation studies. A set of conclusions is presented as well as a discussion of the technical implementation of GCR simulation.

3DHZETRN: Shielded ICRU spherical phantom

A computationally efficient 3DHZETRN code capable of simulating High (H) Charge (Z) and Energy (HZE) and light ions (including neutrons) under space-like boundary conditions with enhanced neutron and light ion propagation was recently developed for a simple homogeneous shield object. Monte Carlo benchmarks were used to verify the methodology in slab and spherical geometry, and the 3D corrections were shown to provide significant improvement over the straight-ahead approximation in some cases. In the present report, the new algorithms with well-defined convergence criteria are extended to inhomogeneous media within a shielded tissue slab and a shielded tissue sphere and tested against Monte Carlo simulation to verify the solution methods. The 3D corrections are again found to more accurately describe the neutron and light ion fluence spectra as compared to the straight-ahead approximation. These computationally efficient methods provide a basis for software capable of space shield analysis and optimization.

Dynamic/Anisotropic Low Earth Orbit Environmental Models

The first level of space infrastructure is represented by the International Space Station (ISS), which is now permanently occupied by space career astronauts. It is imperative that we understand the ISS exposures dynamically for career planning and insure that the regulatory requirement of keeping exposures as low as reasonably achievable (ALARA) is adequately implemented. This is especially true as ISS matures with increasing complexity resulting in a larger drag coefficient requiring operation at higher altitudes with increased exposure rates. In applying ALARA, it is desirable to augment the shielding in areas where the crew spends significant amounts of time. One such area is the crew quarter that has the added problem of being generally located against the pressure vessel wall where shielding is minimal and exposure rates tend to be high. A significant portion of the exposure results from trapped protons and is highly directional. In addition, some directions of approach of galactic cosmic rays are shielded below the Earth’s horizon and greater numbers of low energy particles enter from the west than the east. Clearly, the augmentation requires an understanding of this directionality to block major radiation components. We describe herein a new model of trapped and galactic cosmic radiation as seen in ISS operations in preparation for the design of an augmentation of the ISS crew quarter.

Solar proton exposure of an ICRU sphere within a complex structure part II: Ray-trace geometry.

A computationally efficient 3DHZETRN code with enhanced neutron and light ion (Z ≤ 2) propagation was recently developed for complex, inhomogeneous shield geometry described by combinatorial objects. Comparisons were made between 3DHZETRN results and Monte Carlo (MC) simulations at locations within the combinatorial geometry, and it was shown that 3DHZETRN agrees with the MC codes to the extent they agree with each other. In the present report, the 3DHZETRN code is extended to enable analysis in ray-trace geometry. This latest extension enables the code to be used within current engineering design practices utilizing fully detailed vehicle and habitat geometries. Through convergence testing, it is shown that fidelity in an actual shield geometry can be maintained in the discrete ray-trace description by systematically increasing the number of discrete rays used. It is also shown that this fidelity is carried into transport procedures and resulting exposure quantities without sacrificing computational efficiency.

Spectral Analyses and Radiation Exposures from Several Ground-Level Enhancement (GLE) Solar Proton Events: A Comparison of Methodologies

Several methods for analyzing the particle spectra from extremely large solar proton events, called Ground-Level Enhancements (GLEs), have been developed and utilized by the scientific community to describe the solar proton energy spectra and have been further applied to ascertain the radiation exposures to humans and radio-sensitive systems, namely electronics. In this paper 12 GLEs dating back to 1956 are discussed, and the three methods for describing the solar proton energy spectra are reviewed. The three spectral fitting methodologies are EXP [an exponential in proton rigidity (R)], WEIB [Weibull fit: an exponential in proton energy], and the Band function (BAND) [a double power law in proton rigidity]. The EXP and WEIB methods use low energy (MeV) GLE solar proton data and make extrapolations out to approx.1 GeV. On the other hand, the BAND method utilizes low- and medium-energy satellite solar proton data combined with high-energy solar proton data deduced from high-latitude neutron monitoring stations. Thus, the BAND method completely describes the entire proton energy spectrum based on actual solar proton observations out to ~10 GeV. Using the differential spectra produced from each of the 12 selected GLEs for each of the three methods, radiation exposures are presented and discussed in detail. These radiation exposures are then compared with the current 30-day and annual crew exposure limits and the radiation effects to electronics.

Dependence of the Martian radiation environment on atmospheric depth: Modeling and measurement

The energetic particle environment on the Martian surface is influenced by solar and heliospheric modulation and changes in the local atmospheric pressure (or column depth). The Radiation Assessment Detector (RAD) on board the Mars Science Laboratory rover Curiosity on the surface of Mars has been measuring this effect for over four Earth years (about two Martian years). The anticorrelation between the recorded surface Galactic Cosmic Ray-induced dose rates and pressure changes has been investigated by Rafkin et al. (2014) and the long-term solar modulation has also been empirically analyzed and modeled by Guo et al. (2015). This paper employs the newly updated HZETRN2015 code to model the Martian atmospheric shielding effect on the accumulated dose rates and the change of this effect under different solar modulation and atmospheric conditions. The modeled results are compared with the most up-to-date (from 14 August 2012 to 29 June 2016) observations of the RAD instrument on the surface of Mars. Both model and measurements agree reasonably well and show the atmospheric shielding effect under weak solar modulation conditions and the decline of this effect as solar modulation becomes stronger. This result is important for better risk estimations of future human explorations to Mars under different heliospheric and Martian atmospheric conditions.

3DHZETRN: Neutron leakage in finite objects

The 3DHZETRN formalism was recently developed as an extension to HZETRN with an emphasis on 3D corrections for neutrons and light ions. Comparisons to Monte Carlo (MC) simulations were used to verify the 3DHZETRN methodology in slab and spherical geometry, and it was shown that 3DHZETRN agrees with MC codes to the degree that various MC codes agree among themselves. One limitation of such comparisons is that all of the codes (3DHZETRN and three MC codes) utilize different nuclear models/databases; additionally, using a common nuclear model is impractical due to the complexity of the software. It is therefore difficult to ascertain if observed discrepancies are caused by transport code approximations or nuclear model differences. In particular, an important simplification in the 3DHZETRN formalism assumes that neutron production cross sections can be represented as the sum of forward and isotropic components, where the forward component is subsequently solved within the straight-ahead approximation. In the present report, previous transport model results in specific geometries are combined with additional results in related geometries to study neutron leakage using the Webber 1956 solar particle event as a source boundary condition. A ratio is defined to quantify the fractional neutron leakage at a point in a finite object relative to a semi-infinite slab geometry. Using the leakage ratio removes some of the dependence on the magnitude of the neutron production and clarifies the effects of angular scattering and absorption with regard to differences between the models. Discussion is given regarding observed differences between the MC codes and conclusions drawn about the need for further code development. Although the current version of 3DHZETRN is reasonably accurate compared to MC simulations, this study shows that improved leakage estimates can be obtained by replacing the isotropic/straight-ahead approximation with more detailed descriptions.

Low Earth orbit assessment of proton anisotropy using AP8 and AP9 trapped proton models

The completion of the International Space Station (ISS) in 2011 has provided the space research community with an ideal evaluation and testing facility for future long duration human activities in space. Ionized and secondary neutral particles radiation measurements inside ISS form the ideal tool for validation of radiation environmental models, nuclear reaction cross sections and transport codes. Studies using thermo-luminescent detectors (TLD), tissue equivalent proportional counter (TPEC), and computer aided design (CAD) models of early ISS configurations confirmed that, as input, computational dosimetry at low Earth orbit (LEO) requires an environmental model with directional (anisotropic) capability to properly describe the exposure of trapped protons within ISS. At LEO, ISS encounters exposure from trapped electrons, protons and geomagnetically attenuated galactic cosmic rays (GCR). For short duration studies at LEO, one can ignore trapped electrons and ever present GCR exposure contributions during quiet times. However, within the trapped proton field, a challenge arises from properly estimating the amount of proton exposure acquired. There exist a number of models to define the intensity of trapped particles. Among the established trapped models are the historic AE8/AP8, dating back to the 1980s and the recently released AE9/AP9/SPM. Since at LEO electrons have minimal exposure contribution to ISS, this work ignores the AE8 and AE9 components of the models and couples a measurement derived anisotropic trapped proton formalism to omnidirectional output from the AP8 and AP9 models, allowing the assessment of the differences between the two proton models. The assessment is done at a target point within the ISS-11A configuration (circa 2003) crew quarter (CQ) of Russian Zvezda service module (SM), during its ascending and descending nodes passes through the south Atlantic anomaly (SAA). The anisotropic formalism incorporates the contributions of proton narrow pitch angle (PA) and east-west (EW) effects. Within SAA, the EW anisotropy results in different level of exposure to each side of the ISS Zvezda SM, allowing angular evaluation of the anisotropic proton spectrum. While the combined magnitude of PA and EW effects at LEO depends on a multitude of factors such as trapped proton energy, orientation and altitude of the spacecraft along the velocity vector, this paper draws quantitative conclusions on the combined anisotropic magnitude differences within ISS SM target point between AP8 and AP9 models.