Engineering topics

Dorelli, John C.

Publications and source records attributed to Dorelli, John C..

14 records

NASA NTRS2019

Does the Solar Wind Electric Field Control the Reconnection Rate at Earth's Subsolar Magnetopause?

We demonstrate that the reconnection rate at the subsolar magnetopause is stronglycontrolled by the solar wind electric field and depends weakly on the local properties of the dissipationregion. Our approach is to match the solar wind and magnetospheric states in an internal boundarylayer described by the Cassak and Shay (2007, https://doi.org/10.1063/1.2795630) expression fortwo-dimensional asymmetric reconnection. Faraday's law along the Sun-Earth line determinesthe variation of the solar wind electric field from the bow shock to the magnetopause. While themagnetospheric plasma exerts some control over the reconnection rate, magnetic flux pileup in the sheathpartially compensates for any local reduction in the reconnection rate. For a fixed magnetospheric state,the reconnection rate is shown to be directly proportional to the solar wind electric field, thus explainingwhy the solar wind electric field correlates well with geomagnetic indices.

Dorelli, John C.↗

NASA NTRS2017

Spacecraft and Instrument Photoelectrons Measured by the Dual Electron Spectrometers on MMS

Secondary electrons are continuously generated via photoemission from sunlit spacecraft and instrument surfaces. These particles can subsequently contaminate low-energy channels of electron sensors. Spacecraft photoelectrons are measured at energies below that of a positive spacecraft potential and can be removed at the expense of energy resolution. However, fluxes of photoelectrons generated inside electron instruments are independent of spacecraft potential and must be fully characterized in order to correct electron data. Here we present observations of spacecraft and instrument photoelectron populations measured with the Dual Electron Spectrometers (DES) on NASA's Magnetospheric Multiscale (MMS) mission. We leverage observations from Earth's nightside plasma sheet taken during MMS commissioning and develop an empirical model of instrument photoelectrons. This model is used with DES velocity distribution functions to correct plasma moments and has been made publicly available on the MMS science data center for use by the scientific community.

photoelectrons↗

NASA NTRS2016

Ion-Scale Structure in Mercury's Magnetopause Reconnection Diffusion Region

The strength and time dependence of the electric field in a magnetopause diffusion region relate to the rate of magnetic reconnection between the solar wind and a planetary magnetic field. Here we use approximately 150 milliseconds measurements of energetic electrons from the Mercury Surface, Space Environment, GEochemistry, and Ranging (MESSENGER) spacecraft observed over Mercury's dayside polar cap boundary (PCB) to infer such small-scale changes in magnetic topology and reconnection rates. We provide the first direct measurement of open magnetic topology in flux transfer events at Mercury, structures thought to account for a significant portion of the open magnetic flux transport throughout the magnetosphere. In addition, variations in PCB latitude likely correspond to intermittent bursts of approximately 0.3 to 3 millivolts per meter reconnection electric fields separated by approximately 5 to10 seconds, resulting in average and peak normalized dayside reconnection rates of approximately 0.02 and approximately 0.2, respectively. These data demonstrate that structure in the magnetopause diffusion region at Mercury occurs at the smallest ion scales relevant to reconnection physics.

Gershman, Daniel J.↗

NASA NTRS2016

Electron Dynamics in a Subproton-Gyroscale Magnetic Hole

Magnetic holes are ubiquitous in space plasmas, occurring in the solar wind, downstream of planetary bow shocks, and inside the magnetosphere. Recently, kinetic-scale magnetic holes have been observed near Earth's central plasma sheet. The Fast Plasma Investigation on NASA's Magnetospheric Multiscale (MMS) mission enables measurement of both ions and electrons with 2 orders of magnitude increased temporal resolution over previous magnetospheric instruments. Here we present data from MMS taken in Earth's nightside plasma sheet and use high-resolution particle and magnetometer data to characterize the structure of a subproton-scale magnetic hole. Electrons with gyroradii above the thermal gyroradius but below the current layer thickness carry a current sufficient to account for a 10-20 depression in magnetic field magnitude. These observations suggest that the size and magnetic depth of kinetic-scale magnetic holes is strongly dependent on the background plasma conditions.

Gershman, Daniel J.↗

NASA NTRS2016

Extended Magnetohydrodynamics with Embedded Particle-in-Cell Simulation of Ganymede's Magnetosphere

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

Toth, Gabor↗

NASA NTRS2015

In Flight Calibration of the Magnetospheric Multiscale Mission Fast Plasma Investigation

The Fast Plasma Investigation (FPI) on the Magnetospheric Multiscale mission (MMS) combines data from eight spectrometers, each with four deflection states, into a single map of the sky. Any systematic discontinuity, artifact, noise source, etc. present in this map may be incorrectly interpreted as legitimate data and incorrect conclusions reached. For this reason it is desirable to have all spectrometers return the same output for a given input, and for this output to be low in noise sources or other errors. While many missions use statistical analyses of data to calibrate instruments in flight, this process is insufficient with FPI for two reasons: 1. Only a small fraction of high resolution data is downloaded to the ground due to bandwidth limitations and 2: The data that is downloaded is, by definition, scientifically interesting and therefore not ideal for calibration. FPI uses a suite of new tools to calibrate in flight. A new method for detection system ground calibration has been developed involving sweeping the detection threshold to fully define the pulse height distribution. This method has now been extended for use in flight as a means to calibrate MCP voltage and threshold (together forming the operating point) of the Dual Electron Spectrometers (DES) and Dual Ion Spectrometers (DIS). A method of comparing higher energy data (which has low fractional voltage error) to lower energy data (which has a higher fractional voltage error) will be used to calibrate the high voltage outputs. Finally, a comparison of pitch angle distributions will be used to find remaining discrepancies among sensors.

Heliophysics↗

NASA NTRS2015

The Parameterization of Top-Hat Particle Sensors with Microchannel-Plate-Based Detection Systems and its Application to the Fast Plasma Investigation on NASA's Magnetospheric MultiScale Mission

The most common instrument for low energy plasmas consists of a top-hat electrostatic analyzer geometry coupled with a microchannel-plate (MCP)-based detection system. While the electrostatic optics for such sensors are readily simulated and parameterized during the laboratory calibration process, the detection system is often less well characterized. Furthermore, due to finite resources, for large sensor suites such as the Fast Plasma Investigation (FPI) on NASA's Magnetospheric Multiscale (MMS) mission, calibration data are increasingly sparse. Measurements must be interpolated and extrapolated to understand instrument behavior for untestable operating modes and yet sensor inter-calibration is critical to mission success. To characterize instruments from a minimal set of parameters we have developed the first comprehensive mathematical description of both sensor electrostatic optics and particle detection systems. We include effects of MCP efficiency, gain, scattering, capacitive crosstalk, and charge cloud spreading at the detector output. Our parameterization enables the interpolation and extrapolation of instrument response to all relevant particle energies, detector high voltage settings, and polar angles from a small set of calibration data. We apply this model to the 32 sensor heads in the Dual Electron Sensor (DES) and 32 sensor heads in the Dual Ion Sensor (DIS) instruments on the 4 MMS observatories and use least squares fitting of calibration data to extract all key instrument parameters. Parameters that will evolve in flight, namely MCP gain, will be determined daily through application of this model to specifically tailored in-flight calibration activities, providing a robust characterization of sensor suite performance throughout mission lifetime. Beyond FPI, our model provides a valuable framework for the simulation and evaluation of future detection system designs and can be used to maximize instrument understanding with minimal calibration resources.

Instrument-model↗

NASA NTRS2015

Performance of a Discrete Wavelet Transform for Compressing Plasma Count Data and its Application to the Fast Plasma Investigation on NASA's Magnetospheric Multiscale Mission

Plasma measurements in space are becoming increasingly faster, higher resolution, and distributed over multiple instruments. As raw data generation rates can exceed available data transfer bandwidth, data compression is becoming a critical design component. Data compression has been a staple of imaging instruments for years, but only recently have plasma measurement designers become interested in high performance data compression. Missions will often use a simple lossless compression technique yielding compression ratios of approximately 2:1, however future missions may require compression ratios upwards of 10:1. This study aims to explore how a Discrete Wavelet Transform combined with a Bit Plane Encoder (DWT/BPE), implemented via a CCSDS standard, can be used effectively to compress count information common to plasma measurements to high compression ratios while maintaining little or no compression error. The compression ASIC used for the Fast Plasma Investigation (FPI) on board the Magnetospheric Multiscale mission (MMS) is used for this study. Plasma count data from multiple sources is examined: resampled data from previous missions, randomly generated data from distribution functions, and simulations of expected regimes. These are run through the compression routines with various parameters to yield the greatest possible compression ratio while maintaining little or no error, the latter indicates that fully lossless compression is obtained. Finally, recommendations are made for future missions as to what can be achieved when compressing plasma count data and how best to do so.

Particles↗

NASA NTRS2015

In Flight Calibration of the Magnetospheric Multiscale Mission Fast Plasma Investigation

The Fast Plasma Investigation (FPI) on the Magnetospheric Multiscale mission (MMS) combines data from eight spectrometers, each with four deflection states, into a single map of the sky. Any systematic discontinuity, artifact, noise source, etc. present in this map may be incorrectly interpreted as legitimate data and incorrect conclusions reached. For this reason it is desirable to have all spectrometers return the same output for a given input, and for this output to be low in noise sources or other errors. While many missions use statistical analyses of data to calibrate instruments in flight, this process is difficult with FPI for two reasons: 1. Only a small fraction of high resolution data is downloaded to the ground due to bandwidth limitations and 2: The data that is downloaded is, by definition, scientifically interesting and therefore not ideal for calibration. FPI uses a suite of new tools to calibrate in flight. A new method for detection system ground calibration has been developed involving sweeping the detection threshold to fully define the pulse height distribution. This method has now been extended for use in flight as a means to calibrate MCP voltage and threshold (together forming the operating point) of the Dual Electron Spectrometers (DES) and Dual Ion Spectrometers (DIS). A method of comparing higher energy data (which has low fractional voltage error) to lower energy data (which has a higher fractional voltage error) will be used to calibrate the high voltage outputs. Finally, a comparison of pitch angle distributions will be used to find remaining discrepancies among sensors.

Detection-system↗

NASA NTRS2013

A Simple GPU-Accelerated Two-Dimensional MUSCL-Hancock Solver for Ideal Magnetohydrodynamics

We describe our experience using NVIDIA's CUDA (Compute Unified Device Architecture) C programming environment to implement a two-dimensional second-order MUSCL-Hancock ideal magnetohydrodynamics (MHD) solver on a GTX 480 Graphics Processing Unit (GPU). Taking a simple approach in which the MHD variables are stored exclusively in the global memory of the GTX 480 and accessed in a cache-friendly manner (without further optimizing memory access by, for example, staging data in the GPU's faster shared memory), we achieved a maximum speed-up of approx. = 126 for a sq 1024 grid relative to the sequential C code running on a single Intel Nehalem (2.8 GHz) core. This speedup is consistent with simple estimates based on the known floating point performance, memory throughput and parallel processing capacity of the GTX 480.

graphics processing units↗

NASA NTRS2011

Porting a Hall MHD Code to a Graphic Processing Unit

We present our experience porting a Hall MHD code to a Graphics Processing Unit (GPU). The code is a 2nd order accurate MUSCL-Hancock scheme which makes use of an HLL Riemann solver to compute numerical fluxes and second-order finite differences to compute the Hall contribution to the electric field. The divergence of the magnetic field is controlled with Dedner?s hyperbolic divergence cleaning method. Preliminary benchmark tests indicate a speedup (relative to a single Nehalem core) of 58x for a double precision calculation. We discuss scaling issues which arise when distributing work across multiple GPUs in a CPU-GPU cluster.

Dorelli, John C.↗

NASA NTRS2010

Multiscale Simulations of Magnetic Island Coalescence

We describe a new interactive parallel Adaptive Mesh Refinement (AMR) framework written in the Python programming language. This new framework, PyAMR, hides the details of parallel AMR data structures and algorithms (e.g., domain decomposition, grid partition, and inter-process communication), allowing the user to focus on the development of algorithms for advancing the solution of a systems of partial differential equations on a single uniform mesh. We demonstrate the use of PyAMR by simulating the pairwise coalescence of magnetic islands using the resistive Hall MHD equations. Techniques for coupling different physics models on different levels of the AMR grid hierarchy are discussed.

Dorelli, John C.↗

NASA NTRS2010

Does the Hall Effect Solve the Flux Pileup Saturation Problem?

It is well known that magnetic flux pileup can significantly speed up the rate of magnetic reconnection in high Lundquist number resistive MHD,allowing reconnection to proceed at a rate which is insensitive to the plasma resistivity over a wide range of Lundquist number. Hence, pileup is a possible solution to the Sweet-Parker time scale problem. Unfortunately, pileup tends to saturate above a critical value of the Lundquist number, S_c, where the value ofS_c depends on initial and boundary conditions, with Sweet-Parker scaling returning above S_c. It has been argued (see Dorelli and Bim [2003] and Dorelli [2003]) that the Hall effect can allow flux pileup to saturate (when the scale of the current sheet approaches ion inertial scale, di) before the reconnection rate begins to stall. However, the resulting saturated reconnection rate, while insensitive to the plasma resistivity, was found to depend strongly on the di. In this presentation, we revisit the problem of magnetic island coalescence (which is a well known example of flux pileup reconnection), addressing the dependence of the maximum coalescence rate on the ratio of di in the "large island" limit in which the following inequality is always satisfied: l_eta di lambda, where I_eta is the resistive diffusion length and lambda is the island wavelength.

Dorelli, John C.↗

NASA NTRS2009

The Role of Flux Transfer Events in Solar Wind-Magnetosphere Coupling

Recent global MHD simulations of Flux Transfer Events (FTE's)[Dorelli and Bhattachar ee, JGR, 114, 2009] demonstrate that they are topologically complex flux ropes which extend large distances away from the subsolar magnetopause. Thus, FTE's represent a significant perturbation to the magnetopause magnetic field topology. Specifically, prior to the formation of the first FTE of the simulation, the dayside magnetopause has a relatively simple topology consisting of a single magnetic separator draping northward (even for southward IMF) over the dayside magnetopause (in a manner similar to the "overdraping" which occurs under northward IMF conditions). After the first FTE forms, however, the magnetopause becomes topologically more complex, consisting of multiple separators and braided flux domains. One consequence of this topological complexity is that open flux tubes have direct access to the magnetosphere near the subsolar region (such direct access is not possible during steady separator reconnection). Since the magnetopause topology affects the mapping of the solar wind electric field to the magnetosphere, an interesting question arises: What is the effect of FTE generation on the global magnetospheric convection pattern? In this talk, we address this question by exploring the differences between magnetospheric and ionospheric convection before and after the formation of the first FTE of the simulation.

Dorelli, John C.↗