Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Space plasma physics”

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.

At least 883 records · Page 49

Efficient Nonthermal Ion and Electron Acceleration Enabled by the Flux-Rope Kink Instability in 3D Nonrelativistic Magnetic Reconnection

The relaxation of field-line tension during magnetic reconnection gives rise to a universal Fermi acceleration process involving the curvature drift of particles. However, the efficiency of this mechanism is limited by the trapping of energetic particles within flux ropes. Using 3D fully kinetic simulations, we demonstrate that the flux-rope kink instability leads to strong field-line chaos in weak-guide-field regimes where the Fermi mechanism is most efficient, thus allowing particles to transport out of flux ropes and undergo further acceleration. As a consequence, both ions and electrons develop clear power-law energy spectra that contain a significant fraction of the released energy. The low-energy bounds are determined by the injection physics, while the high-energy cutoffs are limited only by the system size. Furthermore, these results have strong relevance to observations of nonthermal particle acceleration in space and astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

SOHO Ultraviolet Coronagraph Spectrometer (UVCS) Mission Operations and Data Analysis

The scientific goal of UVCS is to obtain detailed empirical descriptions of the extended solar corona as it evolves over the solar cycle and to use these descriptions to identify and understand the physical processes responsible for coronal heating, solar wind acceleration, coronal mass ejections (CMEs), and the phenomena that establish the plasma properties of the solar wind as measured by "in situ" solar wind instruments. This report covers the period from 15 February 2003 to 14 April 2004. During that time, UVCS observations have consisted of three types: 1) standard synoptic observations comprising, primarily, the H I Lyalpha line profile and the 0 VI 103.2 and 103.7 nm intensity over a range of heights from 1.5 to about 3.0 solar radii and covering 360 degrees about the Sun, 2) sit and stare observations for major flare watches, and 3) special observations designed by the UVCS Lead Observer of the Week for a specific scientific purpose. The special observations are often coordinated with those of other space-based and ground-based instruments and they often are part of SOHO joint observation programs and campaigns. Lead observers have included UVCS Co-Investigators, scientists from the solar physics community and several graduate and undergraduate level students. UVCS has continued to achieve its purpose of using powerful spectroscopic diagnostic techniques to obtain a much more detailed description of coronal structures and dynamic phenomena than existed before the SOHO mission. The new descriptions of coronal mass ejections (CMEs) and coronal structures from UVCS have inspired a large number of theoretical studies aimed at identifying the physical processes responsible for CMEs and solar wind acceleration in coronal holes and streamers. UVCS has proven to be a very stable instrument. Stellar observations have demonstrated its radiometric stability. UVCS has not required any flight software modifications and all mechanisms are operational. The UVCS 0 VI Channel with its redundant optical path for wavelengths near H I Lyalpha is capable of observing the entire UVCS wavelength range. The regions of the detector currently being used require different grating angles for direct OVI observations and redundant path H I Lyalpha observations, and so those can no longer be observed simultaneously. Since December 1998, the 0 VI Channel has been used for all UVCS observations. Although the H I Lyalpha Channel and detector are still operational, increases in the dark count up to about 5x10(exp 4) counts/sec/pixel and an increase in high voltage current to within a factor of two of the maximum used in the laboratory before flight led to the decision to not use that detector after 1998. The visible light channel functioned nominally during the reporting period. UVCS data, data analysis software, calibration files and the mission log are available from the SOHO archive and SAO. All UVCS data are now available within three months of the observations to scientists and the general public via the SOHO Data Archive and SAO. UVCS has resulted in 33 scientific papers in 2003. There were numerous presentations at scientific meetings. UVCS Education and Public Outreach activities involved nine members of the UVCS team. During the reporting period, there were over a dozen events directed at students and teachers, museum audiences, and public audiences via the mass media, internet and educational literature.

Gurman, Joseph↗

Computing the shape gradient of stellarator coil complexity with respect to the plasma boundary

Coil complexity is a critical consideration in stellarator design. The traditional two-step optimization approach, in which the plasma boundary is optimized for physics properties and the coils are subsequently optimized to be consistent with this boundary, can result in plasma shapes which cannot be produced with sufficiently simple coils. To address this challenge, we propose a method to incorporate considerations of coil complexity in the optimization of the plasma boundary. Coil complexity metrics are computed from the current potential solution obtained with the REGCOIL code (Landreman, Nucl. Fusion , vol. 57, 2017, 046003). While such metrics have previously been included in derivative-free fixed-boundary optimization (Drevlak et al. , Nucl. Fusion , vol. 59, 2018, 016010), we compute the local sensitivity of these metrics with respect to perturbations of the plasma boundary using the shape gradient (Landreman & Paul, Nucl. Fusion , vol. 58, 2018, 076023). We extend REGCOIL to compute derivatives of these metrics with respect to parameters describing the plasma boundary. In keeping with previous research on winding surface optimization (Paul et al. , Nucl. Fusion , vol. 58, 2018, 076015), the shape derivatives are computed with a discrete adjoint method. In contrast with the previous work, derivatives are computed with respect to the plasma surface parameters rather than the winding surface parameters. To further reduce the resolution required to compute the shape gradient, we present a more efficient representation of the plasma surface which uses a single Fourier series to describe the radial distance from a coordinate axis and a spectrally condensed poloidal angle. This representation is advantageous over the standard cylindrical representation used in the VMEC code (Hirshman & Whitson, Phys. Fluids , vol. 26, 1983, pp. 3553–3568), as it provides a uniquely defined poloidal angle, eliminating a null space in the optimization of the plasma surface. In comparison with previous spectral condensation methods (Hirshman & Breslau, Phys. Plasmas , vol. 5, 1998, p. 2664), the modified poloidal angle is obtained algebraically rather than through the solution of a nonlinear optimization problem. The resulting shape gradient highlights features of the plasma boundary that are consistent with simple coils and can be used to couple coil and fixed-boundary optimization.

Physics↗

A Comprehensive Review of Latent Space Dynamics Identification Algorithms for Intrusive and Non-Intrusive Reduced-Order-Modeling

Numerical solvers of partial differential equations (PDEs) have been widely employed for simulating physical systems. However, the computational cost remains a major bottleneck in various scientific and engineering applications, which has motivated the development of reduced-order models (ROMs). Recently, machine-learning-based ROMs have gained significant popularity and are promising for addressing some limitations of traditional ROM methods, especially for advection dominated systems. In this chapter, we focus on a particular framework known as Latent Space Dynamics Identification (LaSDI), which transforms the high-fidelity data, governed by a PDE, to simpler and low-dimensional latent-space data, governed by ordinary differential equations (ODEs). These ODEs can be learned and subsequently interpolated to make ROM predictions. Each building block of LaSDI can be easily modulated depending on the application, which makes the LaSDI framework highly flexible. In particular, we present strategies to enforce the laws of thermodynamics into LaSDI models (tLaSDI), enhance robustness in the presence of noise through the weak form (WLaSDI), select high-fidelity training data efficiently through active learning (gLaSDI, GPLaSDI), and quantify the ROM prediction uncertainty through Gaussian processes (GPLaSDI). We demonstrate the performance of different LaSDI approaches on Burgers equation, a non-linear heat conduction problem, and a plasma physics problem, showing that LaSDI algorithms can achieve relative errors of less than a few percent and up to thousands of times speed-ups.

Computational Engineering, Finance, and Science (c↗

Exploring the parameter space of MagLIF implosions using similarity scaling. I. Theoretical framework

Magneto-inertial fusion concepts, such as the magnetized liner inertial fusion (MagLIF) platform, constitute an alternative path for achieving ignition and significant fusion yields in the laboratory. The space of experimental input parameters defining a MagLIF load is highly multi-dimensional, and the implosion itself is a complex event involving many physical processes. In the first paper of this series, we develop a simplified analytical model that identifies the main physical processes at play during a MagLIF implosion. Using non-dimensional analysis, we determine the most important dimensionless parameters characterizing MagLIF implosions and provide estimates of such parameters using typical fielded or experimentally observed quantities for MagLIF. Here, we then show that MagLIF loads can be “incompletely” similarity scaled, meaning that the experimental input parameters of MagLIF can be varied such that many (but not all) of the dimensionless quantities are conserved. Based on similarity-scaling arguments, we can explore the parameter space of MagLIF loads and estimate the performance of the scaled loads. Then, in the follow-up papers of this series, we test the similarity-scaling theory for MagLIF loads against simulations for two different scaling “vectors,” which include current scaling and rise-time scaling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct imaging of plasma waves using ultrafast electron microscopy

A femtosecond plasma imaging modality based on a new development of ultrafast electron microscope is introduced. We investigated the laser-induced formation of high-temperature electron microplasmas and their subsequent non-equilibrium evolution. Based on a straightforward field imaging principle, we directly retrieve detailed information about the plasma dynamics, including plasma wave structures, particle densities, and temperatures. We discover that directly subjected to a strong magnetic field, the photo-generated microplasmas manifest in novel transient cyclotron echoes and form new wave states across a broad range of field strengths and different laser fluences. Intriguingly, the transient cyclotron waves morph into a higher frequency upper-hybrid wave mode with the dephasing of local cyclotron dynamics. The quantitative real-space characterizations of the non-equilibrium plasma systems demonstrate the feasibilities of a new microscope system in studying the plasma dynamics or transient electric fields with high spatiotemporal resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neural network surrogate models for equations of state

Equation of state (EOS) data provide necessary information for accurate multiphysics modeling, which is necessary for fields such as inertial confinement fusion. Here, we suggest a neural network surrogate model of energy and entropy and use thermodynamic relationships to derive other necessary thermodynamic EOS quantities. We incorporate phase information into the model by training a phase classifier and using phase-specific regression models, which improves the modal prediction accuracy. Our model predicts energy values to 1% relative error and entropy to 3.5% relative error in a log-transformed space. Although sound speed predictions require further improvement, the derived pressure values are accurate within 10% relative error. Our results suggest that neural network models can effectively model EOS for inertial confinement fusion simulation applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conceptual Design of a Scintillator-Based Fast-Ion Loss Detector for the Wendelstein 7-X Stellarator

Here, a conceptual design of a scintillator-based fast-ion loss detector (FILD) has been developed for the Wendelstein 7-X stellarator (W7-X). Simulations using the Monte Carlo codes ASCOT5 and BEAMS3D have been performed to calculate the expected flux of neutral beam injection (NBI)-generated fast hydrogen ions onto the conceptual detector probe head. These fast-ion loss fluxes have been calculated for several magnetic field configurations as well as probe insertion positions. At the maximum insertion position, both co-and counter-going losses with high incident pitch angles are observed; however, at retracted positions, only co-going fast ions reach the probe head. The FILDSIM code has been used to optimize the geometry of the detector entrance and collimating elements to achieve a wide velocity space coverage as well as a high velocity–space resolution. A synthetic FILD signal is calculated for the expected loss distribution via forward modeling using the instrument response function. The synthetic signal is found to largely retain the velocity space features of the loss distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CUSP Energetic Particles: Confinement, Acceleration and Implications

The cusp energetic particle (CEP) event is a new magnetospheric phenomenon. The events were detected in the dayside cusp for hours, in which the measured helium ions had energies up to 8 MeV. All of these events were associated with a dramatic decrease and large fluctuations in the local magnetic field strength. During January 1999 - December 1999 covered by this report, I have studied the CEP events by analyzing the POLAR, GEOTAIL, and WIND particle and magnetic field data measured during the geomagnetic quiet periods in 1996 and one geomagnetic storm period in 1998. The simultaneous observations indicated that the ion fluxes in the CEP events were higher than that in both the upstream and the downstream from the bow shock. The pitch angle distribution of the helium ions in the CEP events was found to peak around 90 deg. It was found that the mirror parameter, defined as the ratio of the square root of the integration of the parallel turbulent power spectral component over the ultra-low frequency (ULF) ranges to the mean field in the cusp, is correlated with the intensity of the cusp MeV helium flux, which is a measure of the influence of mirroring interactions and an indication of local effect. It was also found that the turbulent power of the local magnetic field in the ultra-low frequency (ULF) ranges is correlated with the intensity of the cusp energetic helium ions. Such ULF ranges correspond to periods of about 0.33-500 seconds that cover the gyroperiods, the bounce periods, and the drift periods of the tens keV to MeV charged particles when they are temporarily confined in the high-altitude dayside cusp. These observations represent a discovery that the high-altitude dayside cusp is a new acceleration and dynamic trapping region of the magnetosphere. The cusp geometry is connected via gradient and curvature drift of these energized ions to the equatorial plasma sheet as close as the geostationary orbit at local midnight. It implies that the dayside cusp is potentially an important source of magnetospheric particles. The discovery of the CEP events has been recognized as one of the most significant results from POLAR. I was invited to give a talk at 1999 IUGG meeting to interpret the CEP events. This discovery has also been written into the web-based Space Physics Text Book (http://www.oulu.fi/- spaceweb/textbook/cusp.html).

Chen, Jiasheng↗

Atmospheric, Magnetospheric, and Plasmas in Space (AMPS) spacelab payload definition study, technical summary document

Some 60 instrument candidates and 80 possible science investigations were evaluated. The early analysis emphasized the science aspect in terms of the functional requirements for each of the potential experiments identified by the AMPS science working group. These requirements were then used for the grouping of instruments into practical payloads which would fit the capabilities of the Shuttle/Spacelab. This analysis resulted in the definition of eleven different AMPS configurations. The data were then used to define a typical set of requirements for a flexible AMPS laboratory. The data gathered to this point showed that a planned sequential buildup of the laboratory would be necessary to meet both physical and funding limitations. This led to the definition of five strawman payloads by the science working group, which were used to establish a conceptual laboratory and to define preliminary design of a configuration which could satisfy AMPS needs during the early program period.

Keeley, J. T.↗

Modeling laser-driven ion acceleration with deep learning

Developments in machine learning promise to ameliorate some of the challenges of modeling complex physical systems through neural-network-based surrogate models. High-intensity, short-pulse lasers can be used to accelerate ions to mega-electronvolt energies, but to model such interactions requires computationally expensive techniques such as particle-in-cell simulations. Multilayer neural networks allow one to take a relatively sparse ensemble of simulations and generate a surrogate model that can be used to rapidly search the parameter space of interest. In this study, we created an ensemble of over 1,000 simulations modeling laser-driven ion acceleration and developed a surrogate to study the resulting parameter space. A neural-network-based approach allows for rapid feature discovery not possible for traditional parameter scans given the computational cost. A notable observation made during this study was the dependence of ion energy on the pre-plasma gradient length scale. While this methodology harbors great promise for ion acceleration, it has ready application to all topics in which large-scale parameter scans are restricted by significant computational cost or relatively large, but sparse, domains.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗

A three dimensional dynamic study of electrostatic charging in materials

A description is given of the physical models employed in the NASCAP (NASA Charging Analyzer Program) code, and several test cases are presented. NASCAP dynamically simulates the charging of an object made of conducting segments which may be entirely or partially covered with thin dielectric films. The object may be subject to either ground test or space user-specified environments. The simulation alternately treats (1) the tendency of materials to accumulate and emit charge when subject to plasma environment, and (2) the consequent response of the charged particle environment to an object's electrostatic field. Parameterized formulations of the emission properties of materials subject to bombardment by electrons, protons, and sunlight are presented. Values of the parameters are suggested for clean aluminum, Al2O3, clean magnesium, MgO, SiO2 kapton, and teflon. A discussion of conductivity in thin dielectrics subject to radiation and high fields is given, together with a sample calculation.

Katz, I.↗

The Princeton Field-Reversed Configuration for Compact Nuclear Fusion Power Plants

The Princeton Field-Reversed Configuration (PFRC) nuclear fusion reactor concept is an innovative approach to fusion power generation prioritizing low neutron production and small size. Here, a combination of analytical modeling and numerical simulation shows that the novel heating approach generates an FRC with closed field lines. Simulation data from a single-particle Hamiltonian code predicts ms-scale plasma heating in reactor-scale conditions while PIC codes predict formation of warm FRC plasmas from initial mirror fields. The PFRC-1 and PFRC-2 experiments have heated electrons to energies well in excess of 100 eV and plasma durations to 300 ms, more than 10 4 times longer than the predicted tilt instability growth time. From these data, we have created a development plan and anticipated performance metrics for a fusion reactor based on the PFRC concept. The resulting 1–10 MW PFRC reactors would be suitable for diverse applications, from submarines to urban environments to space propulsion. PFRC is a steady-state, driven magnetic confinement device. Plasma, inside a cylindrical array of coils, is confined and heated by external RF antennae. PFRC would be ultra-low radiation due to both its fuel and small size. The choice of advanced fuels, deuterium and helium-3 (D– 3 He), may be enabled by the high-β FRC configuration. The small size of the reactor would enable rapid exhaust of the dangerous tritium ash. Low radiation would make the reactor safer to operate and, in combination with simple geometry and small size, dramatically lowers development and maintenance costs. This review paper gives an introduction to the physics of the PFRC and a summary of the PFRC-2 experiment results to date. It then discusses the future program plan and how PFRC reactors would be commercialized.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experimental study of hot electron generation in shock ignition relevant high-intensity regime with large scale hot plasmas

In the shock ignition (SI) laser fusion scheme, hot electrons generated by the laser spike pulse can either preheat the fuel or strengthen the ignition shock, depending on the hot electron characteristics. We conducted a planar target experiment on the OMEGA-EP laser facility and characterized the temperature and total energy of hot electrons generated from a kilojoule-class 100-ps infrared (IR) or a 1-ns ultraviolet (UV) laser interacting with a large (L n ~330-450 μm) and hot (T e ~1-2 keV) coronal plasma at the SI-relevant intensities (~10 16 W/cm 2 ). The IR laser converts ~2.5% energy into hot electrons with T hot ~ 60–90 keV, while the UV laser couples 0.8% ± 0.7% energy into T hot =27 ± 9 keV hot electrons. The IR-produced hot electrons yield five times higher Cu Kα emission than the UV case, confirming the higher electron conversion efficiency with the IR laser. The low energy conversion from the UV laser to hot electrons may be due to the refraction of the off-normal incident laser in the large coronal plasma. These findings are the first comparisons of hot electron generation between the IR and UV pulses at kilojoule scales in SI-relevant large-scale plasmas. The findings may expand the SI design space to include IR lasers as the possible spike lasers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Excitation of zonal flow by nonlinear geodesic acoustic mode

Here, the zero frequency zonal flow (ZFZF) excitation due to the nonlinear geodesic acoustic mode (GAM) is investigated in the framework of gyro-kinetic equations, which were first investigated. We show that after integrating over the velocity space, the first order finite orbit width (FOW) term of nonlinear GAM disappears, and hence, the second order FOW term should be taken into account to generate the ZFZF. The anisotropy of the equilibrium is also considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine Learning for Mapping Multipactor Susceptibility in RF Systems: Capabilities and Generalization Constraints

Multipactor is a surface-driven electron avalanche phenomenon that degrades the performance and reliability of radio-frequency (RF) systems in particle accelerator and vacuum electronics applications. Multipactor behavior in a given device structure is conventionally assessed through susceptibility charts, which provide a parameter-space characterization of the instability. In this work, we assess the capabilities of machine-learning (ML) models to learn and predict such susceptibility charts and analyze the constraints governing their generalization across materials. Using a simulation-derived dataset spanning six distinct secondary-electron-yield material profiles in a canonical two-surface planar geometry, we train supervised regression models and artificial neural networks to predict the time-averaged electron growth rate, δavg, across the relevant parameter space. Model performance is evaluated using metrics that explicitly probe the structure of susceptibility charts, including Intersection over Union, Structural Similarity Index, and correlation analysis. Tree-based ensemble models outperform neural-network models in reconstructing susceptibility regions and in generalizing across material domains. Principal-component analysis reveals disjoint material feature distributions, indicating that the piecewise mode structure of multipactor susceptibility is difficult to represent with a single global model and that generalization is constrained by data coverage rather than by model complexity. An exhaustive reduced-coverage study further shows that sparse material-space coverage can yield mean performance in the same general range but producing large variability in the susceptibility-region overlap. These results clarify the capabilities of ML-based surrogate models for parameter-space characterization of multipactor discharge. They also provide guidance for their appropriate use in RF system design.

43 PARTICLE ACCELERATORS↗

Localizing resonant magnetic perturbations for edge localized mode control in KSTAR

An external 3D magnetic perturbation typically drives a resonant response at the rational surfaces from the core to the edge of tokamak plasmas, due to strong mode coupling and amplification. In this paper, we present a method to isolate the edge from core resonant fields using the ideal perturbed equilibrium code and to design an edge-localized resonant magnetic perturbation (RMP) for effective edge localized mode (ELM) control. A robust feature of the edge-localized RMP is the curtailed response to the field at the low-field-side (LFS) midplane, as opposed to typical RMPs which strongly resonate with the LFS fields. This emphasizes the importance of off-midplane coils to improve ELM control without provoking a large core response that could lead to devastating instabilities. The conceptual design of new ELM control coils based on the edge-localized RMP in KSTAR shows how this new insight can be utilized to enhance the efficiency of our ELM suppression capabilities. Simple window-pane coils matching the edge-localized resonant mode structure substantially expand in the ELM suppression window beyond the existing coil. Further optimization using the flexible optimized coils using space-curves code leads to additional enhancement in the edge-localized control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗