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.     

Cosmic ray hits in the central nervous system at solar maximum.

It has been suggested that a manned mission to Mars be launched at solar maximum rather than at solar minimum to minimize the radiation exposure to galactic cosmic rays. It is true that the number of hits from highly ionizing particles to critical regions in the brain will be less at solar maximum, and it is of interest to estimate how much less. We present here calculations for several sites within the brain from iron ions (z = 26) and from particles with charge, z, greater than or equal to 15. The same shielding configurations and sites in the brain used in an earlier paper for solar minimum are employed so that direct comparison of results between the two solar activity conditions can be made. A simple pressure-vessel wall and an equipment room onboard a spacecraft are chosen as shielding examples. In the equipment room, typical results for the thalamus are that the probability of any particles with 7 greater than or equal to 15 and from 2.3 percent to 1.3 percent for iron ions. The extra shielding provided in the equipment room makes little difference in these numbers. We conclude that this decrease in hit frequency (less than a factor of two) does not provide a compelling reason to avoid solar minimum for a manned mission to Mars. This conclusion could be revised, however, if a very small number of hits is found to cause critical malfunction within the brain.     

Global ena Image Simulations

The energetic neutral atom (ENA) images obtained by the ISEE and POLAR satellites pointed the way toward global imaging of the magnetospheric plasmas. The Imager for Magnetopause to Aurora Global Exploration (IMAGE) is the first mission to dedicate multiple neutral atom imagers: HENA, MENA and LENA, to monitor the ion distributions in high-, medium- and low-energy ranges, respectively. Since the start of science operation, HENA, MENA and LENA have been continuously sending down images of the ring current, ionospheric outflow, and magnetosheath enhancements from high pressure solar wind. To unfold multiple-dimensional (equal or greater than 3) plasma distributions from 2-dimensional images is not a trivial task. Comparison with simulated ENA images from a modeled ion distribution provides an important basis for interpretation of features in the observed images. Another approach is to develop image inversion methods to extract ion information from ENA images. Simulation studies have successfully reproduced and explained energetic ion drift dynamics, the transition from open to closed drift paths, and the magnetosheath response to extreme solar wind conditions. On the other hand, HENA has observed storm-time ion enhancement on the nightside toward dawn that differs from simple concepts but can be explained using more sophisticated models. LENA images from perigee passes reveal unexpected characteristics that now can be interpreted as evidence for a transient superthermal exospheric component that is gravitationally-influenced if not bound. In this paper, we will report ENA simulations performed during several IMAGE observed events. These simulations provide insight and explanations to the ENA features that were not readily understandable previously.     

Electrical Charging of Volcanic Plumes

Many explosive terrestrial volcanic eruptions are accompanied by lightning and other atmospheric electrical phenomena. The plumes produced generate large perturbations in the surface atmospheric electric potential gradient and high charge densities have been measured on falling volcanic ash particles. The complex nature of volcanic plumes (which contain gases, solid particles, and liquid drops) provides several possible charging mechanisms. For plumes rich in solid silicate particles, fractoemission (the ejection of ions and atomic particles during fracture events) is probably the dominant source of charge generation. In other plumes, such as those created when lava enters the sea, different mechanisms, such as boiling, may be important. Further charging mechanisms may also subsequently operate, downwind of the vent. Other solar system bodies also show evidence for volcanism, with activity ongoing on Io. Consequently, volcanic electrification under different planetary scenarios (on Venus, Mars, Io, Moon, Enceladus, Tethys, Dione and Triton) is also discussed.     

Helium-3: Status and Prospects

We report on our continuing efforts to determine 3He abundances in H II regions and planetary nebulae. Our detections of 3He in some PNe show that some stars produce large amounts of 3He. However the H II region abundances show no evidence for this production. From our sample of > 40 H II regions, the subsample which should yield the most reliable abundances has 3He/H abundances which scatter between 1-2 × 10-5. There is no trend with either galactocentric distance or metallicity. Even if we do not understand the underlying mechanisms, we see empirically that stars neither produce nor destroy 3He in a major way. We thus suggest that the level of the "3He Plateau" (3He/H = 1.5 _-0.5 ^+1.0 × 10^-5) is a reasonable estimate for the primordial 3He.     

Computational methods for the HZETRN code.

Asymptotic expansion has been used to simplify the transport of high charge and energy ions for broad beam applications in the laboratory and space. The solution of the lowest order asymptotic term is then related to a Green's function for energy loss and straggling coupled to nuclear attenuation providing the lowest order term in a rapidly converging Neumann series for which higher order collisions terms are related to the fragmentation events including energy dispersion and downshift. The first and second Neumann corrections were evaluated numerically as a standard for further analytic approximation. The first Neumann correction is accurately evaluated over the saddle point whose width is determined by the energy dispersion and located at the downshifted ion collision energy. Introduction of the first Neumann correction leads to significant simplification of the second correction term allowing application of the mean value theorem and a second saddle point approximation. The regular dependence of the second correction spectral dependence lends hope to simple approximation to higher corrections. At sufficiently high energy nuclear cross-section variations are small allowing non-perturbative methods to all orders and renormalization of the second corrections allow accurate evaluation of the full Neumann series.     

New ways to rotation rates

Progress in determining rotation rates of Algol primary stars from light curves is assessed. Included are elementary ideas on why the new method works, how to carry out solutions, and how and when to apply the constraint on rotational lobe filling. Some comments are made on desirable improvements for measuring rotation from line broadening. A few remarks are made about the problem of identifying likely fast rotating Algols.     

Track Initiation of Occasionally Unresolved Radar Targets

In a dense detection environment track-while-scan algorithms will introduce many false tracks when processing is performed on a scan-by-scan basis. The maximum likelihood solution involving all detections on all scans is formulated and evaluated for the initiation problem consisting of false alarms, missed detections, and unresolved detections. A comparison is made between the maximum likelihood solutions including and not including the resolution likelihood and it is shown that the resolution likelihood must be included if the data contain unresolved detections. While the maximum likelihood method (with unresolved detections) does not appear to be implementable in real time without very high speed integrated circuit (VHSIC) technology, it can be used as a standard to which more implementable methods can be compared.     

Observations of the profiles of solar UV emission lines and their analysis in terms of the heating and production of the corona

Observations of the solar spectrum have been made between 1200–2200 with high spectral resolution. The results were obtained with an all-reflecting echelle spectrograph carried by a stabilized Skylark rocket launched in April 1970. Measurements of the profiles of a number of emission lines due to Siii, Cii, Siiii and Civ formed in the temperature range 10⁴-10⁵ K, indicate ion energies which are considerably in excess of the electron temperatures derived from the ionization balance. Since the ion/electron relaxation time is very short the observed ion energies cannot correspond to an ion temperature and hence a non-thermal mechanical energy component exists in the transition zone.It is postulated that the non-thermal energy component represents the actual mechanical energy responsible for the heating of the corona, and, that, it is propagated as an acoustic wave. On this basis and with a preliminary estimate of the reflection from the transition zone, a flux of 3 × 10⁵ erg cm ^-2 s ^-1 is established as entering the corona. This value is in agreement with estimates of the total energy loss from the corona due to conduction, radiation and the solar wind, thus establishing a gross energy balance.Theoretical calculations are currently underway to establish the physical nature of the atmosphere which would result from such a propagating flux. At the present time this has been carried out for an atmosphere in hydrostatic equilibrium and the energy balance equation solved. A preliminary temperature structure which results is shown in Figure 1, together with the derived distribution in electron density. This gives a corona of the right temperature and density but the observed structure deviates in detail from those derived from an analysis of the solar XUV spectrum.     

Analysis of Limit Cycles in a Two-Transistor Saturable-Core Parallel Inverter

A familiar two-transistor saturable-core parallel inverter is modeled as a nonlinear negative resistance in parallel with energy-storage elements. The techniques of singular-point analysis are combined with piecewise linear techniques to permit determination of solution trajectories on the phase plane. Clear insight is provided, not only into steady-state oscillation, but also into transient behavior of the circuit. Experimental results confirming the analytical model are included.