Search NASA⌕ Search

Engineering topics

Cunningham, Gregory Scott

Publications and source records attributed to Cunningham, Gregory Scott.

Resolution of a few problems in the application of quasilinear theory to calculating diffusion coefficients in heliophysics [Slides]

Three problems with the historical literature were discussed. Kennel and Engelmann (1966): ‘notational error’ that implies that diffusion coefficients depend on the distribution of wave power with respect to ψ. This dependence is unphysical and is not present in Lyons+ 1971. Kennel and Engelmann (1966) assume that only a single frequency can be associated with each wave vector: not necessarily true for EMIC, e.g. Lyons (1974a, 1974b) incorrectly used the Jacobian to relate the power spectral density in (ω,χ) space to that into (κ ⊥ , κ ∥ space. The correct use of the Jacobian: (1) Eliminates undesirable properties associated with N(ω) : lack of integrability, lack of superposition; and (2) Can produce orders of magnitude changes to the diffusion coefficients, but depends on the wave distribution: bigger effect for more oblique waves, and depends on the plasma parameter: bigger effect for smaller $ω_{pe}/ω_{ce}$.

79 ASTRONOMY AND ASTROPHYSICS↗

Scattering by whistler-mode waves during a quiet period perturbed by substorm activity

We study the dynamics of radiation belt electrons during a 10-day quiet period perturbed by substorm activity and preceding a high-speed stream (HSS), aiming at a global description of the radiation belts in L-shell, L in [2, 6], and energy [0.1, 10] MeV. We combine Van Allen Probes observations and Fokker-Planck numerical simulations of pitch-angle diffusion. The Fokker-Planck model uses event-driven pitch angle diffusion coefficients from whistler-mode waves, built from the wave properties and the ambient plasma density measurements from the Van Allen Probes. We first find this event has some similar characteristics to regular quiet times previously studied; a widely extended plasmasphere within which we observe strong and varying whistler-mode waves. These ambient conditions lead to strong pitch-angle scattering, which contributes to the creation of a wide slot region as well as a significant decay of the outer radiation belts, which are observed and qualitatively well simulated. In addition, we find the substorm activity causes short duration (within ± 4h) decay of the plasma density and a lowering amplitude of the whistler-mode waves within the plasmasphere, both causing opposite effects in terms of pitch angle diffusion. This leads to a diminution of pitch-angle diffusion at the time of the main substorm activity. Conversely, whistler-mode waves become enhanced in the time periods between the substorm injections. All effects cumulated, we find an enhancement of pitch angle diffusion by whistler-mode waves above L~4.7 during the 10-day period. This directly relates to the combination of quietness and substorm activity which allows pitch angle diffusing of up to 1 MeV electrons in the outer belt. Relativistic electrons of 1–2 MeV remain trapped in the outer belt, from L~4.7 to L~5.2, forming, in both the observations and the simulations, a distinct pocket of remnant electrons.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Radiation Effects of Strontium Ions on Satellite Solar Cells in Low Earth Orbits

This study focuses on the radiation effects 19 of Sr + ions—generated from high-altitude nuclear explosions (HANE)—on satellite solar cells in low-Earth orbits (LEO). Along four selected satellite orbits, ion fluences are sampled inside the evolving Sr + distributions for days, determined from our newly developed HANE environment model. These fluences, along with the help of radiation transport codes including the MULASSIS and SRIM models, enable us to quantify the radiation damages by determining the values of total ionizing doses and the equivalent 1 MeV electron fluences for displacement damages. Comparing the dose values to existing experimental data, we conclude that HANE-generated Sr + ions have limited darkening effects to quartz solar cell coverglasses in LEO with apogees of 100s to 1000 km. In addition, with the extremely high equivalent fluences, we also conclude that these Sr+ ions may cause severe or even fatal displacement damage to exposed solar photovoltaic (PV) cells on satellites in LEO. The radiation effects of Sr + ions are much less significant for the orbits with high apogees beyond ten thousand km. We also conducted model parameter sensitivity studies on the charge exchange cross-sections, neutral atmosphere density profiles and explosion local time positions, and the above conclusions stay unchanged. The methodology developed in this study can be extended to other HANE-generated heavy ion species in the future.

14 SOLAR ENERGY↗

Analysis of Electric and Magnetic Lightning-Generated Wave Amplitudes Measured by the Van Allen Probes

Herein we provide a statistical analysis of both electric and magnetic field wave amplitudes of very low frequency lightning-generated waves (LGWs) based on the equivalent of 11.5 years of observations made by the Van Allen Probes encompassing ~24.6 × 10 6 survey mode measurements. We complement this analysis with data from the ground-based World Wide Lightning Location Network to explore differences between satellite and ground-based measurements. LGW mean amplitudes are found to be low compared with other whistler mode waves (1 ± 1.6 pT and 19 ± 59 μV/m). Extreme events (1/5,000) can reach 100 pT and contributes strongly to the mean power below $L$ = 2. We find excellent correlations between World Wide Lightning Location Network-based power and wave amplitudes in space at various longitudes. We reveal strong dayside ionospheric damping of the LGW electric field. LGW amplitudes drop for $L$ < 2, contrary to the Earth's intense equatorial lightning activity. We conclude that it is difficult for equatorial LGW to propagate and remain at $L$ < 2.

Electric wave power↗

Response to Comment on “Radiation-Belt Remediation Using Space-Based Antennas and Electron Beams” by G. Ganguli and C. Crabtree

Ganguli and Crabtree have written a comment about a recent article by the authors listed above on radiation-belt remediation. They have objected to our evaluation of the Naval Research Laboratory’s chemical release concept which states that this concept may be impractical due to an apparently low overall efficiency. In their comment, they provide a scientific argument and refer to the published literature to counter our statement. Here, in this response, we provide more details on our numerical calculations and experimental results which led to this evaluation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Local and statistical maps of lightning-generated wave power density estimated at the Van Allen Probes footprints from the World-Wide Lightning Location Network database

We introduce a new method that uses the World-Wide Lightning Location Network (WWLLN) to estimate both the local and the drift lightning power density at the Van Allen Probes footprints during 4.3 years (~2x10 8 strokes.). The ratio of the drift power density to the local power density defines a time-resolved WWLLN-based model of lightning-generated wave (LGW) power density ratio, $\mathcal{R_{WWLLN}}$. $\mathcal{R_{WWLLN}}$ is computed every ~34s. This ratio multiplied by the time-resolved LGW intensity measured by the Probes allows direct computation of pitch angle diffusion coefficients used in radiation belts codes. According to Statistical analysis, the median power density ratio is $\mathcal{\hat{R}_{WWLLN}}$ = 0.3 - 4 over the Americas. Elsewhere, $\mathcal{\hat{R}_{WWLLN}}$ > 1 in general. Over oceans, $\mathcal{\hat{R}_{WWLLN}}$ is larger than ~10 and varies with season (x2.5 from winter to summer). The yearly median $\mathcal{\hat{R}_{WWLLN}}$ decays as $\mathcal{\hat{R}_{WWLLN}}$~9.9/L 0.91 . The strong geographical and temporal variation should be kept in assessing effects in space. $\mathcal{R_{WWLLN}}$> 1 suggests significant LGW effects in the inner belt.

58 GEOSCIENCES↗

Incorporation of drift loss cone effects on the trapping of an artificial radiation belt into LANL's modeling capability

Los Alamos National Laboratory has developed an end-to-end model that quantifies the total ionizing dose that a shielded component on a space asset will accrue due to interactions with trapped beta-decay electrons produced after a high-altitude nuclear explosion. The model consists of three major components: initial trapping, spatio-temporal evolution of the trapped population, and interaction of the evolving trapped population with a space asset to produce total ionizing dose. The initial trapping module calculates the fraction of betas produced by radioactive decay of nuclear fission products that will be trapped in the magnetic field surrounding Earth for at least one drift around Earth. Previously, the initial trapping calculation assumed a centered dipole magnetic field model whereas we know that Earth’s magnetic field is offset from its geographic center and has higher-order non-dipole terms. Both the offset and the higher-order terms are captured in the International Geomagnetic Reference Field (IGRF) internal field model. In a centered dipole, the electrons that are trapped on a magnetic field line at one longitude will be trapped at all other longitudes as they drift around Earth due to azimuthal symmetry of the field, but the offset and higher-order terms of the IGRF field affects the drift of the trapped electrons and causes some that are trapped at one longitude to be lost at other longitudes. These electrons are called "quasi-trapped" and are said to be in the drift loss cone since they will be lost during their first drift around Earth. In this paper we report on the effect of incorporating the drift of electrons around Earth using the IGRF field on the trapping fraction as compared to using a centered dipole. Here, we find that the new, more accurate, methodology produces higher trapping fractions for bursts above the South Atlantic Anomaly, a region of low magnetic field strength, and lower trapping fractions for bursts in regions with higher magnetic field strength than expected from a centered dipole.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On the time needed to reach an equilibrium structure of the radiation belts

In this paper, we complement the notion of equilibrium states of the radiation belts with a discussion on the dynamics and time needed to reach equilibrium. We solve for the equilibrium states obtained using 1D radial diffusion with recently developed hiss and chorus lifetimes at constant values of Kp = 1, 3 and 6. We find that the equilibrium states at moderately low Kp, when plotted vs L-shell (L) and energy (E), display the same interesting S-shape for the inner edge of the outer belt as recently observed by the Van Allen Probes. The S-shape is also produced as the radiation belts dynamically evolve toward the equilibrium state when initialized to simulate the buildup after a massive dropout or to simulate loss due to outward diffusion from a saturated state. Physically, this shape, intimately linked with the slot structure, is due to the dependence of electron loss rate (originating from wave-particle interactions) on both energy and L-shell. Equilibrium electron flux profiles are governed by the Biot number (τ Diffusion /τ loss ), with large Biot number corresponding to low fluxes and low Biot number to large fluxes. The time it takes for the flux at a specific (L, E) to reach the value associated with the equilibrium state, starting from these different initial states, is governed by the initial state of the belts, the property of the dynamics (diffusion coefficients), and the size of the domain of computation. Its structure shows a rather complex scissor form in the (L, E) plane. The equilibrium value (phase space density or flux) is practically reachable only for selected regions in (L, E) and geomagnetic activity. Convergence to equilibrium requires hundreds of days in the inner belt for E > 300 keV and moderate Kp (≤3). It takes less time to reach equilibrium during disturbed geomagnetic conditions (Kp ≥ 3), when the system evolves faster. Restricting our interest to the slot region, below L = 4, we find that only small regions in (L, E) space can reach the equilibrium value: E ~ [200, 300] keV for L = [3.7, 4] at Kp = 1, E ~ [0.6, 1] MeV for L = [3, 4] at Kp = 3, and E ~ 300 keV for L = [3.5, 4] at Kp = 6 assuming no new incoming electrons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗