SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Dipole antenna radiation in compressible anisotropic electron plasma overlying imperfectly conducting half space, solving for lunar environment
The removal of undesired beam nonlinear energy chirp (time-energy correlation) or linearization of the beam longitudinal phase space (LPS) is crucial for high-brightness linac-based scientific applications, such as x-ray free-electron lasers. In this paper, we propose that a low-density hollow channel plasma can be used as a near-ideal passive linearizer to significantly linearize the beam LPS and, at the same time, preserve the beam emittance. Physically, the passage of the beam through the hollow plasma channel excites a strong quasicosinoidal longitudinal wakefield that acts to mitigate the beam nonlinear energy chirp by superimposing a reverse chirp on the beam. The theoretical analyses and large-scale three-dimensional start-to-end simulations confirm that the beam longitudinal phase space can be almost completely linearized without noticeable beam emittance growth. In conclusion, application of such a near-ideal linearizer may significantly improve the performance of numerous accelerator-based applications.
The electrodynamic tether consists of a satellite deployed to a distance of 20 km by an electrically conducting tether. The space station hardware consists of a 12 meter deployment boom, satellite cradle, tether reel and motor, and other tether support systems. The electrodynamic tether will be used to perform a variety of wave experiments by exciting a wide spectrum of low frequency waves in the ionospheric plasma. The system can also be used to artificially generate and study field aligned currents and associated plasma effects. Hydromagnetic waves generated by the passage of the system through the space plasma are of particular interest in space plasma research.
The large-scale dynamics of most conventional space and astrophysical plasmas are predominantly governed by Alfvén modes, which are low-frequency magnetohydrodynamic modes existing in magnetized media. At scales smaller than the ion gyroscale or frequencies exceeding the ion cyclotron frequency, the Alfvén modes transform into kinetic-Alfvén or whistler modes that significantly contribute to plasma dynamics. However, this scenario reverses in nonneutral pair plasmas, such as those found in the magnetospheres of pulsars and magnetars, around rotating black holes, and in their relativistic jets, as well as in certain laboratory plasmas. In these systems, the large-scale dynamics are governed by hybrid whistler–Alfvén modes, which transform into pure Alfvén modes at smaller scales. We derive the nonlinear equations that describe the dynamics of whistler–Alfvén modes in ultrarelativistic nonneutral magnetically dominated pair plasma and discuss the spectrum of turbulence governed by these equations.
Hypervelocity impacts on solid surfaces generate and propagate shock waves and create extreme density, pressure, and temperature conditions. Understanding hypervelocity impacts requires a multi-physics approach at the intersection of fluid dynamics, solid mechanics, and quantum theory, which has the capability of spanning multiple different dominating physics regimes. These strong shock conditions found in hypervelocity impacts are studied in the context of high-speed collisions between spacecraft and meteoroids and debris. We study the problem of hypervelocity impacts on metals with the lens of understanding the processes that lead to the creation of an impact generated plasma. In order to characterize the impacting particle as well as its threat to spacecraft from the subsequent radiating plasma an analysis was performed that highlights how the kinetic energy of the impactor is deposited into a target, how the shock wave and shocked material in the target evolve, and how these effects can lead towards the creation of an impact plasma through heating. The analysis suggests that at low velocities the majority of the energy from a hypervelocity impact is stored in reversible elastic energy and does not contribute to strong heating. At higher impact velocities, nucleus and electron thermal components dominate, increasing the energy available for ionization. Overall, this work provides the complete set of state variable initial conditions for any study of hypervelocity impact plasma expansion and its effects.
Imbalanced Alfvénic turbulence is a universal process playing a crucial role in energy transfer in space, astrophysical, and laboratory plasmas. A fundamental and long-lasting question about the imbalanced Alfvénic turbulence is how and through which mechanism the energy transfers between scales. Here, we show that the energy transfer of imbalanced Alfvénic turbulence is completed by coherent interactions between Alfvén waves and co-propagating anomalous fluctuations. These anomalous fluctuations are generated by nonlinear couplings instead of linear reflection. We also reveal that the energy transfer of the waves and the anomalous fluctuations is carried out mainly through local-scale and large-scale nonlinear interactions, respectively, responsible for their bifurcated power-law spectra. This work unveils the energy transfer physics of imbalanced Alfvénic turbulence, and advances the understanding of imbalanced Alfvénic turbulence observed by Parker Solar Probe in the inner heliosphere.
Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.
The tritium aspects of the DT fuel cycle embody some of the most challenging feasibility and attractiveness issues in the development of fusion systems. The review and analyses in this paper provide important information to understand and quantify these challenges and to define the phase space of plasma physics and fusion technology parameters and features that must guide a serious R&D in the world fusion program. We focus in particular on components, issues and R&D necessary to satisfy three 'principal requirements': (1) achieving tritium self-sufficiency within the fusion system, (2) providing a tritium inventory for the initial start-up of a fusion facility, and (3) managing the safety and biological hazards of tritium. A primary conclusion is that the physics and technology state-of-the-art will not enable DEMO and future power plants to satisfy these principal requirements. We quantify goals and define specific areas and ideas for physics and technology R&D to meet these requirements. A powerful fuel cycle dynamics model was developed to calculate time-dependent tritium inventories and flow rates in all parts and components of the fuel cycle for different ranges of parameters and physics and technology conditions. Dynamics modeling analyses show that the key parameters affecting tritium inventories, tritium start-up inventory, and tritium self-sufficiency are the tritium burn fraction in the plasma (f b ), fueling efficiency (η f ), processing time of plasma exhaust in the inner fuel cycle (t p ), reactor availability factor (AF), reserve time (tr) which determines the reserve tritium inventory needed in the storage system in order to keep the plant operational for time t r in case of any malfunction of any part of the tritium processing system, and the doubling time (t d ). Results show that η f f b > 2% and processing time of 1–4 h are required to achieve tritium self-sufficiency with reasonable confidence. For η f f b = 2% and processing time of 4 h, the tritium start-up inventory required for a 3 GW fusion reactor is ~11 kg, while it is <5 kg if η f f b = 5% and the processing time is 1 h. To achieve these stringent requirements, a serious R&D program in physics and technology is necessary. The EU-DEMO direct internal recycling concept that carries fuel directly from the plasma exhaust gas to the fueling systems without going through the isotope separation system reduces the overall processing time and tritium inventories and has positive effects on the required tritium breeding ratio (TBR R ). A significant finding is the strong dependence of tritium self-sufficiency on the reactor availability factor. Simulations show that tritium self-sufficiency is: impossible if AF < 10% for any η f f b , possible if AF > 30% and 1% ≤ η f f b ≤ 2%, and achievable with reasonable confidence if AF > 50% and η f f b > 2%. These results are of particular concern in light of the low availability factor predicted for the near-term plasma-based experimental facilities (e.g. FNSF, VNS, CTF), and can have repercussions on tritium economy in DEMO reactors as well, unless significant advancements in RAMI are made. There is a linear dependency between the tritium start-up inventory and the fusion power. The required tritium start-up inventory for a fusion facility of 100 MW fusion power is as small as 1 kg. Since fusion power plants will have large powers for better economics, it is important to maintain a 'reserve' tritium inventory in the tritium storage system to continue to fuel the plasma and avoid plant shutdown in case of malfunctions of some parts of the tritium processing lines. But our results show that a reserve time as short as 24 h leads to unacceptable reserve and start-up inventory requirements. Therefore, high reliability and fast maintainability of all components in the fuel cycle are necessary in order to avoid the need for storing reserve tritium inventory sufficient for continued fusion facility operation for more than a few hours. The physics aspects of plasma fueling, tritium burn fraction, and particle and power exhaust are highly interrelated and complex, and predictions for DEMO and power reactors are highly uncertain because of lack of experiments with burning plasma. Fueling by pellet injection on the high field side of tokamak has evolved to be the preferred method to fuel a burning plasma. Extrapolation from the DIII-D penetration scaling shows fueling efficiency expected in DEMO to be <25%, but such extrapolations are highly uncertain. The fueling efficiency of gas in a reactor relevant regime is expected to be extremely poor and not very useful for getting tritium into the core plasma efficiently. Gas fueling will nonetheless be useful for feedback control of the divertor operating parameters. Extensive modeling has been carried out to predict burn fraction, fueling requirements, and fueling efficiency for ITER, DEMO, and beyond. The fueling rate required to operate Q = 10 ITER plasmas in order to provide the required core fueling, helium exhaust and radiative divertor plasma conditions for acceptable divertor power loads was calculated. If this fueling is performed with a 50–50 DT mix, the tritium burn fraction in ITER would be ~0.36%, which is too low to satisfy the self-sufficiency conditions derived from the dynamics modeling for fusion reactors. Extrapolation to DEMO using this approach would also yield similarly low burn fraction. Extensive analysis presented shows that specific features of edge neutral dynamics in ITER and fusion reactors, which are different from present experiments, open possibilities for optimization of tritium fueling and thus to improve the burn fraction. Using only tritium in pellet fueling of the plasma core, and only deuterium for edge density, divertor power load and ELM control results in significant increase of the burn fraction to 1.8–3.6%. These estimates are performed with physics models whose results cannot be fully validated for ITER and DEMO plasma conditions since these cannot be achieved in present tokamak experiments. Thus, several uncertainties remain regarding particle transport and scenario requirements in ITER and DEMO. The safety standard requirements for protection of the public and release guidelines for tritium have been reviewed. General safety approaches including minimizing tritium inventories, reducing tritium permeation through materials, and decontaminating material for waste disposal have been suggested.
Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady energy transfer rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. Finally, these results provide explicit evidence that reconnection dynamics fundamentally involves turbulence-like energy transfer.
Noncollisional current sheets that form during the nonlinear development of spontaneous magnetic reconnection are characterized by a small thickness, of the order of the electron skin depth. They can become unstable to the formation of plasmoids, which allows the magnetic reconnection process to reach high reconnection rates. In this work, we investigate the marginal stability conditions for the development of plasmoids when the forming current sheet is purely collisionless and in the presence of a strong guide field. We analyze the geometry that characterizes the reconnecting current sheet, and what promotes its elongation. Once the reconnecting current sheet is formed, we identify the regimes for which it is plasmoid unstable. Finally, our study shows that plasmoids can be obtained, in this context, from current sheets with an aspect ratio much smaller than in the collisional regime, and that the plasma flow channel of the marginally stable current layers maintains an inverse aspect ratio of 0.1.
The mechanisms that generate “seed” magnetic fields in our Universe and that amplify them throughout cosmic time remain poorly understood. By means of fully kinetic particle-in-cell simulations of turbulent, initially unmagnetized plasmas, here we study the genesis of magnetic fields via the Weibel instability and follow their dynamo growth up to near-equipartition levels. In the kinematic stage of the dynamo, we find that the rms magnetic field strength grows exponentially with rate γ B ≃0.4u rms /L, where L/2π is the driving scale and u rms is the rms turbulent velocity. In the saturated stage, the magnetic field energy reaches about half of the turbulent kinetic energy. Here, magnetic field growth is balanced by dissipation via reconnection, as revealed by the appearance of plasmoid chains. At saturation, the integral-scale wave number of the magnetic spectrum approaches k int ≃12π/L. Our results show that turbulence—induced by, e.g., the gravitational buildup of galaxies and galaxy clusters—can magnetize collisionless plasmas with large-scale near-equipartition fields.
Not Available
For over a decade, the SpacePy project has contributed open-source solutions for the production and analysis of heliophysics data and simulation results. Here we introduce SpacePy’s functionality for the scientific user and present relevant design principles. We examine recent advances and the future of SpacePy in the broader scientific Python ecosystem, concluding with some of the work that has used SpacePy.
The X8.2-class limb flare on 2017 September 10 is among the best studied solar flare events owing to its great similarity to the standard flare model and the broad coverage by multiple spacecraft and ground-based observations. These multiwavelength observations indicate that electron acceleration and transport are efficient in the reconnection and flare looptop regions. However, there lacks a comprehensive model for explaining and interpreting the multi-faceted observations. In this work, we model the electron acceleration and transport in the early impulsive phase of this flare. We solve the Parker transport equation that includes the primary acceleration mechanism during magnetic reconnection in the large-scale flare region modeled by MHD simulations. We find that electrons are accelerated up to several MeV and fill a large volume of the reconnection region, similar to the observations shown in microwaves. The electron spatial distribution and spectral shape in the looptop region agree well with those derived from the microwave and hard X-ray emissions before magnetic islands grow large and dominate the acceleration. Future emission modelings using the electron maps will enable direct comparison with microwave and hard X-ray observations. These results shed new light on the electron acceleration and transport in a broad region of solar flares within a data-constrained realistic flare geometry.
We present a spatial filtering (or coarse-graining) analysis on 3D magnetized magnetohydrodynamic (MHD) turbulence simulations. The filtered compressible MHD formulae show transfer of kinetic and magnetic energies from large to small scales, as well as energy conversion between kinetic, magnetic, and thermal energies. The anisotropic filtering enables separate analyses of the energy flows perpendicular and parallel to the global mean magnetic field. Anisotropy in energy cascade is demonstrated by the larger perpendicular energy cascade rate and also the larger perpendicular wavenumbers associated with the peak energy transfer rate. We also find that the “inertial range” along the parallel (perpendicular) direction in the anisotropic energy cascade formulation is no longer strictly dissipation-free, because it includes the dissipation in the perpendicular (parallel) direction. A change in the driving force (kinetic only versus kinetic and magnetic) affects the energy conversion between kinetic and magnetic energies. While the compressibility of the driving force changes the partition of different channels of energy transfer and conversion, and also increases the total energy transfer rate, the global energy flow remains unaffected by compressibility qualitatively. Our analysis can be applied to multispacecraft observations of turbulence in the solar wind or a planetary magnetosphere.
The structure of the magnetic reconnection boundary, particularly the presence of slow-mode shocks in the near-Earth magnetotail was studied by using magnetospheric multiscale (MMS) observations and 2.5D hybrid simulations. A total of 51 crossings of MMS from 2017 to 2021 were analyzed. We found that the detection percentage of slow-mode shocks in the near-Earth magnetotail is 41%–55%. Previous studies have only reported one slow-mode shock event in the near-Earth magnetotail and a slow-mode shock detection percentage of 10% or lower in the mid-to-distant magnetotail. It was observed that if the high-energy beam region data is removed from the slow-mode shock downstream observations then the detection of slow-mode shocks reduces, implying that the kinetic effects play an important role in the detection of slow-mode shocks. For the crossings where the interface was not identified as a slow-mode shock, it was found that the turbulence in those crossings can change the mass flux values and disrupt the detection of slow-mode shock. However, the macroscopic slow-mode shock-like structure stably exists around the magnetic reconnection interface, as most of the conditions for slow-mode shocks were satisfied. This result suggests that slow-mode shocks are a general feature of magnetic reconnection geometry. We find that the lack of detection of slow-mode shocks in previous observations and simulations can be explained by taking into account the kinetic structure of slow-mode shocks and the presence of turbulence.
The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.
Recent observations and simulations indicate that solar flares undergo extremely complex 3D evolution, making 3D particle transport models essential for understanding electron acceleration and interpreting flare emissions. In this study, we investigate this problem by solving Parker’s transport equation with 3D MHD simulations of solar flares. By examining energy conversion in the 3D system, we evaluate the roles of different acceleration mechanisms, including reconnection current sheet (CS), termination shock (TS), and supra-arcade downflows (SADs). We find that large-amplitude turbulent fluctuations are generated and sustained in the 3D system. The model results demonstrate that a significant number of electrons are accelerated to hundreds of keV and even a few MeV, forming power-law energy spectra. These energetic particles are widely distributed, with concentrations at the TS and in the flare looptop region, consistent with results derived from recent hard X-ray (HXR) and microwave (MW) observations. By selectively turning particle acceleration on or off in specific regions, we find that the CS and SADs effectively accelerate electrons to several hundred keV, while the TS enables further acceleration to MeV. However, no single mechanism can independently account for the significant number of energetic electrons observed. Instead, the mechanisms work synergistically to produce a large population of accelerated electrons. Our model provides spatially and temporally resolved electron distributions in the whole flare region and at the flare footpoints, enabling synthetic HXR and MW emission modeling for comparison with observations. These results offer important insights into electron acceleration and transport in 3D solar flare regions.