Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma fusion”

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 289 records · Page 16

A novel photomultiplier tube neutron time-of-flight detector

A traditional neutron time-of-flight (nTOF) detector used in inertial confinement fusion consists of a scintillator coupled with a photomultiplier tube (PMT). The instrument response function (IRF) of such a detector is dominated by the scintillator-light decay. At DT neutron yields larger than 1013 a novel detector consisting of a microchannel-plate photomultiplier tube in a housing without a scintillator (PMT nTOF) can be used to measure DT yield, ion temperature, and neutron velocity. Most of the neutron signal in PMT nTOF is produced from neutron interaction with a PMT window. The direct interaction of neutrons with MCP has negligible contribution. The elimination of the scintillator removes the scintillator decay from the instrument response function and makes the IRF of the PMT nTOF faster, which makes the ion temperature and neutron velocity measurements more accurate. Three PMT nTOFs were deployed on the OMEGA Laser System for the first time to diagnose inertial confinement fusion plasma. Here, the design details, characteristics, and calibration results of these detectors in DT implosions on OMEGA are presented. Recommendations on the use of different PMTs for specific applications are provided.

47 OTHER INSTRUMENTATION↗

New, improved analysis of correlation ECE data to accurately determine turbulent electron temperature spectra and magnitudes (invited)

Turbulent electron temperature fluctuation measurement using a correlation electron cyclotron emission (CECE) radiometer has become an important diagnostic for studying energy transport in fusion plasmas, and its use is widespread in tokamaks (DIII-D, ASDEX Upgrade, Alcator C-Mod, Tore Supra, EAST, TCV, HL-2A, etc.). The CECE diagnostic typically performs correlation analysis between two closely spaced (within the turbulent correlation length) ECE channels that are dominated by uncorrelated thermal noise emission. This allows electron temperature fluctuations embedded in the thermal noise to be revealed and fluctuation level and spectra determined. We have demonstrated a new, improved CECE coherency-based analysis for calculating the temperature fluctuation frequency spectrum and level, which has been verified both numerically through the simulation of synthetic ECE radiometer data and through analysis of experimental data from the CECE system on DIII-D. The new formulation places coherency-based analysis on a firm foundational footing and corrects some currently published methodologies. This new method accurately accounts for bias error in the coherence function and correctly calculates noise levels for a fixed data record length. It provides excellent accuracy in determining temperature fluctuation level (e.g., <10% error) even for a small realization number in the ensemble average. The method also has a smaller uncertainty (i.e., error bar) in the power spectrum when compared to the more standard cross-power method when evaluated at low coherency. Direct calculation of system noise level using correlation between randomized intermediate frequency signals is recommended.

Wang, G. (ORCID:0000000225739827)↗

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗

Inferring performance metrics for laser direct drive experiments on OMEGA

Quantifying performance improvements on the OMEGA laser facility requires robust inference of established no-alpha performance metrics, which requires, at minimum, a model to infer the shocked fuel mass and pressure of the confined fusion plasma. In this work, we describe the methodology used to infer performance metrics on OMEGA and present the current state-of-the art model used to infer these metrics from OMEGA experiments. In particular, since neutron images of cryogenic implosions are not available on OMEGA at present, we present how x-ray sizes are determined on OMEGA using a Gaussian Process regression model and how the neutron production region's size is inferred from them. As a result, we end by benchmarking the model using synthetic data and 1-D LILAC simulations and test its experimental self-consistency across available x-ray diagnostic channels.

Gopalaswamy, V. [Laboratory for Laser Energetics, ↗

Doppler Backscattering Data Analysis and Integrated Modeling with OMFIT

One Modeling Framework for Integrated Tasks (OMFIT) is a widely used software tool in the magnetic fusion research community. OMFIT provides magnetic fusion energy researchers with a framework for the development of special-purpose physics modules. This paper describes an OMFIT physics module pertaining to the Doppler Backscattering (DBS) fusion plasma diagnostic. DBS measures density fluctuations and flow velocity through plasma scattering of electromagnetic waves. The OMFIT DBS module was developed to analyze experimental DBS data and facilitate modeling of DBS systems installed on multiple tokamak devices. The OMFIT DBS module is designed to support several analysis workflows: detailed analysis of experimental data, experimental planning, and theory-based synthetic diagnostic modeling. The DBS module uses integrated modeling by leveraging other OMFIT physics modules to perform tasks related to DBS, e.g. ray/beam–tracing simulations, edge-localized mode–synchronized data analysis, magnetic equilibrium reconstruction, and fitting kinetic profile data. Furthermore, this paper describes several supported workflows and serves a reference for the OMFIT DBS module.

Doppler backscattering↗

Stability analysis of Alfvén eigenmodes excited by ion cyclotron resonance heating on EAST

Alfvén wave instabilities driven by energetic particles are common in fusion devices. Understanding their behavior is essential for the good confinement of fusion plasma. For this purpose, a series of experiments are performed on EAST to investigate the excitation of toroidal Alfvén eigenmodes (TAEs). Experimentally, it was found that AEs with frequencies around 80, 134, 157 kHz are excited by ion cyclotron resonance heating (ICRH) hydrogen minority heating scenario. Moreover, the excitation of AEs is independent of the wall-coating materials, but determined with the generation and confinement of the fast ions. Statistical analysis suggests that the TAEs can only be excited when the ICRH power is larger than 2 MW and the H 98 factor is larger than 1.15. In line with the experiments, a set of simulations using TORIC, ASCOT, and NOVA-C is performed. The simulation results show that with the experimental density profile and safety factor, the measured TAE with frequency of 134 kHz is well reproduced by taking into account the plasma toroidal rotation frequency. However, the appearance of the 80 kHz mode cannot be captured in the simulations. The ICRH-generated fast hydrogen ions mainly have banana orbits and that their perpendicular velocity is much larger than that of the parallel ions. The stability analysis with NOVA-C suggests that the radial gradient of the fast H ion distribution is the only driving source of TAEs. For TAEs locate in the center of Alfvén continuum, the D ion Landau damping is the dominant damping term; For most other TAEs which have intersections with the continuum boundary, the continuum damping becomes the dominant damping mechanism.

Alfvén eigenmodes↗

Synergistic effects of W–Ta alloys exposed to various ELM-like fusion events

Abstract Recently tungsten alloys have shown capabilities in suppressing morphological changes when exposed to harsh fusion conditions. Morphological changes due to He + and D + ions at high temperatures (1000–2000 K) can lead to easier erosion and sputtering of tungten (W) particles towards quenching the fusion plasma. This study investigated W–3Ta alloy as potential plasma facing materials (PFMs) exposed to different edge-localized mode (ELM)-like energies and ion loadings. Performance of W–3Ta has been tested as a function of ELM-like energies (0.6–1.5 GW m −2 ) and varying He + and D + ion concentrations (100% He + , 50% He + :50% D + ,10% He + :90% D + (similar ELM-like reactor condition), and 100% D + ions). Early-stage fuzz formation and surface pores were observed during 100% He + ion irradiation and mixed ion loadings but not observed during 100% D + ion loadings. Pores decreased in quantity and increased in size with increasing transient heat loading. Grain boundaries acted as trapping sites of helium and deuterium atoms during dual beam ion irradiations, as observed through increased pore size along the grain boundary. Qualitative erosion estimates were measured using witness plates. Increased erosion was observed with increasing D + ion concentration, which is attributed to D supersaturation in the near surface region affecting thermal and mechanical properties. More erosion was observed for W–3Ta samples compared to pure W, evident through large (>10 µ m) cracks on the surface, and an increase in W metal x-ray photoelectron spectroscopy peaks under similar loading conditions. The slight enhancement in erosion is attributed to the slightly decreased thermal properties, ∼10% thermal diffusivity drop, of W–3Ta at high temperatures. The damage threshold was impacted by ion loading, where distinct cracking was observed for 100% He + at a lower heat loading. This matrix study using reactor-like loading supports the need for extensive studies on future W alloy PFMs before selecting one as an alternative.

Physics↗

Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s -1 , agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

43 PARTICLE ACCELERATORS↗

Asynchronous and Load-Balanced Union-Find for Distributed and Parallel Scientific Data Visualization and Analysis

We present a novel distributed union-find algorithm that features asynchronous parallelism and k-d tree based load balancing for scalable visualization and analysis of scientific data. Applications of union-find include level set extraction and critical point tracking, but distributed union-find can suffer from high synchronization costs and imbalanced workloads across parallel processes. In this study, we prove that global synchronizations in existing distributed union-find can be eliminated without changing final results, allowing overlapped communications and computations for scalable processing. We also use a k-d tree decomposition to redistribute inputs, in order to improve workload balancing. We benchmark the scalability of our algorithm with up to 1,024 processes using both synthetic and application data. Here, we demonstrate the use of our algorithm in critical point tracking and super-level set extraction with high-speed imaging experiments and fusion plasma simulations, respectively.

97 MATHEMATICS AND COMPUTING↗

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗

ASGarD: Adaptive Sparse Grid Discretization

Many areas of science exhibit physical processes that are described by high dimensional partial differential equations (PDEs), e.g., the 4D, 5D and 6D models describing magnetized fusion plasmas, models describing quantum chemistry, or derivatives pricing. Such problems are affected by the so-called “curse of dimensionality” where the number of degrees of freedom (or unknowns) required to be solved for scales as N D where N is the number of grid points in any given dimension D. A simple, albeit naive, 6D example is demonstrated in the left panel of Figure 1. With N = 1000 grid points in each dimension, the memory required just to store the solution vector, not to mention forming the matrix required to advance such a system in time, would exceed an exabyte - and also the available memory on the largest of supercomputers available today. The right panel of Figure 1 demonstrates potential savings for a range of problem dimensionalities and grid resolution. While there are methods to simulate such high-dimensional systems, they are mostly based on Monte-Carlo methods, which rely on a statistical sampling such that the resulting solutions include noise. Since the noise in such methods can only be reduced at a rate proportional to $\sqrt{N_p}$ where N p is the number of Monte-Carlo samples, there is a need for continuum, or grid/mesh-based methods for high-dimensional problems, which both do not suffer from noise and bypass the curse of dimensionality. We present a simulation framework that provides such a method using adaptive sparse grids.

97 MATHEMATICS AND COMPUTING↗

Leading magnetic fusion energy science into the big-and-fast data lane

To crack the code of unlimited clean energy from magnetic fusion, researchers collaborate among numerous facilities that are scattered around the globe. Probing hot fusion plasmas with sophisticated diagnostics, they routinely generate fast streams of high-dimensional time-series data. We present the adaptation of a production code-base for a plasma diagnostic to the big-and-fast data lane. The new code-base is used to stream measurements, made in Korea, into a Top-15 computing facility in the US for remote analysis. We discuss software and design choices made for this new HPC-enabled scientific workflow.

Kube, Ralph↗

Laboratory Study of Nonlinear QED in Intense Laser-Matter Interactions

The goal of the project was to investigate the interaction of powerful laser light with matter, which is of fundamental interest to electrodynamics, thermonuclear fusion plasmas, and high-energy-density physics for DOE. Experiments on highly nonlinear relativistic interactions were performed using the state-of-the-art DIOCLES laser at the Extreme Light Laboratory of the University of Nebraska–Lincoln (UNL). We have demonstrated a novel testbed for laser-plasma interactions in which laser pulses collided in plasma with femtosecond and micron-scale accuracy under extreme experiment conditions. This capability is crucial for future radiation reaction experiments and nonlinear QED studies. Moreover, this is also a milestone towards the efficient generation of GeV electron beams via travelling wave laser wakefield accelerator. We proposed a novel idea for measuring the radiation reaction process using travelling wave Thomson scattering, which was studied in numerical simulations. This method allows unambiguous measurement of radiation reaction in a single shot, therefore, surpassing the conventional methods, which suffer severely from the shot-to-shot fluctuations, due to the instability of laser wavefield acceleration.

74 ATOMIC AND MOLECULAR PHYSICS↗

Low-Jitter Clock Receivers for Fast Timing Applications

Precision timing is a key requirement for emerging 4D particle tracking, Positron Emission Tomography (PET), beam and fusion plasma diagnostics, and other systems. Time-to-Digital Converters (TDCs) are commonly used to provide digital estimates of the relative timing between events, but the jitter performance of a TDC can be no better than the performance of the circuits that acquire the pulses and deliver them to the TDC. Several clock receiver and distribution circuits were evaluated, and a differential amplifier with resistive loads driving a pseudo-differential clock distribution network, developed using design guidelines for radiation tolerance and cryogenic compatibility, was fabricated as part of three prototypes: an analog front-end testbed chip for high-precision timing pixel readout, a dedicated TDC evaluation chip, and a Low-Gain Avalanche Detector (LGAD) readout circuit. Based on TDC measurements of the prototypes, we infer that the jitter added by the clock receiver and distribution circuits is less than 2.25 ps-rms. This performance meets the requirements of many future precision timing systems. The clock receiver and on-chip pseudo-differential driver were fabricated in commercial 28-nm CMOS technology and occupy 2288 µm 2 .

47 OTHER INSTRUMENTATION↗

High resolution X-ray optics for astronomical and laboratory sources

In recent years X-ray telescopes have been used successfully for high resolution observations of astronomical objects, in particular the solar corona. At the present time X-ray microscopes are being developed for diagnostics of laser and electron beam induced fusion plasmas. This paper describes both types of imaging systems using as examples the AS&E X-ray telescope on Skylab and a 10X microscope being built for Los Alamos Scientific Laboratories. The designs, manufacturing tolerances, fabrication, testing, and performance are discussed. Steps for further development in both imaging areas are identified.

Silk, J. K.↗

Radiometry using synchrotron radiation

It is possible to use synchrotron radiation from electron synchrotrons and electron storage rings as an absolute source, especially in the VUV and soft X-ray regions where other standards are difficult to find. At the NBS, an electron storage ring (SURF-II) has been used to calibrate spectrometers and photometers utilized in solar and aeronomy research and in fusion plasma diagnostics. The radiation incident on these spectrometers can be calculated to uncertainties of 3%, and a technique to exactly determine the number of electrons orbiting in the ring is currently being developed to reduce this uncertainty. Detector calibrations between 5 and 55 nm are routinely performed at SURF-II and transfer standard detectors with 6-10% uncertainties over the range 5-254 nm are supplied.

Saloman, E. B.↗

A feasibility study and mission analysis for the Hybrid Plume Plasma Rocket

The Hybrid Plume Plasma Rocket (HPPR) is a high power electric propulsion concept which is being developed at the MIT Plasma Fusion Center. This paper presents a theoretical overview of the concept as well as the results and conclusions of an independent study which has been conducted to identify and categorize those technologies which require significant development before the HPPR can be considered a viable electric propulsion device. It has been determined that the technologies which require the most development are high power radio-frequency and microwave generation for space applications and the associated power processing units, low mass superconducting magnets, a reliable, long duration, multi-megawatt space nuclear power source, and long term storage of liquid hydrogen propellant. In addition to this, a mission analysis of a one-way transfer from low earth orbit (LEO) to Mars indicates that a constant acceleration thrust profile, which can be obtained using the HPPR, results in faster trip times and greater payload capacities than those afforded by more conventional constant thrust profiles.

Sullivan, Daniel J.↗