Engineering Papers⌕ Search

Engineering topics

Cowee, Misa M.

Publications and source records attributed to Cowee, Misa M..

8 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulating Long–Term Dynamics of Radiation Belt Electrons Using DREAM3D Model

Here, we compared the performance of DREAM3D simulations in reproducing the long-term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low-energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (DLL), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics inside L* = 6, with a better model performance at lower μ than higher μ, where μ is the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data-driven outer boundary condition is critical to the model performance, while adding in the data-driven seed population doesn't further improve the performance. (b) The model shows comparable performance with DLL from Brautigam and Albert (2000, https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014, https://doi.org/10.1002/2013ja019204), and Liu et al. (2016, https://doi.org/10.1002/2015gl067398), while with DLL from Ali et al. (2016, https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data-based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higher μ using the Wang et al. (2019, https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under-reproduced in the current model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

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↗
DOE OSTI2016

Scalings for the Alfvén-cyclotron instability: Linear dispersion theory and hybrid particle-in-cell simulations: Scalings of Alfven-Cyclotron Instability

The Alfvén-cyclotron instability is driven by an ion temperature anisotropy such that T ⊥ /T || >1 where ⊥ and || denote directions perpendicular and parallel to a uniform background magnetic field B o , respectively. The computations presented here consider a model of a magnetized, homogeneous, collisionless plasma. Two representations of the proton velocity distribution are considered: a single bi-Maxwellian and a magnetospheric-like configuration of two components, a more dense, relatively cool, isotropic component and a less dense, relatively hot, bi-Maxwellian component which drives the instability. Only wave propagation parallel to B o is considered. Here, using numerical solutions of the full kinetic linear dispersion equation, concise analytic expressions for the scaling of the dimensionless maximum instability growth rate and the corresponding dimensionless real frequency are derived as functions of three dimensionless variables: the hot proton temperature anisotropy, the relative hot proton density, and the hot proton β || . Furthermore, using one-dimensional hybrid particle-in-cell simulations of this same instability, a third relation for the scaling of the maximum amplitude of the dimensionless fluctuating magnetic field energy density is derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Predicting electromagnetic ion cyclotron wave amplitude from unstable ring current plasma conditions

Electromagnetic ion cyclotron (EMIC) waves in the Earth's inner magnetosphere are enhanced fluctuations driven unstable by ring current ion temperature anisotropy. EMIC waves can resonate with relativistic electrons and play an important role in precipitation of MeV radiation belt electrons. In this study, we investigate the excitation and saturation of EMIC instability in a homogeneous plasma using both linear theory and nonlinear hybrid simulations. We have explored a four-dimensional parameter space, carried out a large number of simulations, and derived a scaling formula that relates the saturation EMIC wave amplitude to initial plasma conditions. Finally, such scaling can be used in conjunction with ring current models like ring current-atmosphere interactions model with self-consistent magnetic field to provide global dynamic EMIC wave maps that will be more accurate inputs for radiation belt modeling than statistical models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Whistler anisotropy instabilities as the source of banded chorus: Van Allen Probes observations and particle-in-cell simulations

Magnetospheric banded chorus is enhanced whistler waves with frequencies ω r < Ω e , where Ω e is the electron cyclotron frequency, and a characteristic spectral gap at ω r ≃ Ω e /2. This paper uses spacecraft observations and two-dimensional particle-in-cell simulations in a magnetized, homogeneous, collisionless plasma to test the hypothesis that banded chorus is due to local linear growth of two branches of the whistler anisotropy instability excited by two distinct, anisotropic electron components of significantly different temperatures. The electron densities and temperatures are derived from Helium, Oxygen, Proton, and Electron instrument measurements on the Van Allen Probes A satellite during a banded chorus event on 1 November 2012. The observations are consistent with a three-component electron model consisting of a cold (a few tens of eV) population, a warm (a few hundred eV) anisotropic population, and a hot (a few keV) anisotropic population. The simulations use plasma and field parameters as measured from the satellite during this event except for two numbers: the anisotropies of the warm and the hot electron components are enhanced over the measured values in order to obtain relatively rapid instability growth. The simulations show that the warm component drives the quasi-electrostatic upper band chorus and that the hot component drives the electromagnetic lower band chorus; the gap at ~Ω e /2 is a natural consequence of the growth of two whistler modes with different properties.

58 GEOSCIENCES↗
DOE OSTI2012

Update on the Electron Source Model [Slides]

We summarize work done in FY12 on the Los Alamos Electron Source Model (ESM), which predicts the distribution of beta-decay electrons after a high altitude nuclear explosion (HANE) as a function of L, energy, and pitch angle. In the last year we have compared model results with data taken after the Russian 2 HANE test of 1962 and presented results at the HEART conference. We discuss our future plans to continue comparison with HANE data and to develop the code to allow a more complex set of initial conditions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2012

Conditions for Debris-Background Ion Interactions and Collisionless Shock Wave Generation

We use hybrid simulations and simple theoretical arguments to determine when debris ions streaming relative to background ions in a collisionless, magnetized plasma couple strongly enough to generate a magnetosonic shock wave. We consider three types of configurations: one-dimensional, the two-dimensional extension of the 1-D case, and a more complex 2-D geometry that contains some effects that would be found in a laser-produced, laboratory plasma. We show that the simulation results as well as previous Russian and LLNL results reduce to a simple condition (R m /ρ d = equal mass radius/debris ion gyroradius ≥ 0.7) for the generation of a shock wave. Strong debris interaction with the background is characterized by the formation of a magnetic pulse that steepens and speeds up as it encounters the debris ions deflected by the magnetic field. The pulse further evolves into a shock. As the earlier work has indicated, the process also involves the generation of a transverse electric field perpendicular to the flow and the magnetic field that accelerates the background ions radially outward, which in turn causes the speedup of the pulse. With electric and magnetic field probes, the UCLA laser experiments should be able to detect these signatures of coupling as well as the generation of the shock wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗