LLNL Theory and Modeling highlights [Slides]
Abstract not provided.
SEARCH · Search NASA
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.
Abstract not provided.
Predictive simulations of heat and particle fluxes in the plasma boundary are a critical missing capability for magnetic fusion energy power plant design.
Understanding particulate formation in nuclear debris is critical for predicting fallout transport after a nuclear event. Improved characterization of fallout formation and transport could lead to better guidance for emergency response in a post-detonation scenario. By analyzing how U-oxides nucleate onto different forms of SiO 2 (crystalline and amorphous), we can gain insight into how entrained environmental materials may incorporate into fallout. In this experiment a Plasma Flow Reactor (PFR), was used to replicate the high temperature and extreme flow conditions in a nuclear fireball. Uranyl nitrate was injected into the plasma, and PFR-generated Uranium oxides were allowed to condense onto amorphous (nanoparticles of varying sizes) or crystalline (quartz) SiO 2 substrates. SiO 2 substrates were characterized before and after U-oxide deposition using Scanning Electron Microscopy (SEM) based techniques (i.e. EDS) in order to characterize how U-oxides may nucleate onto these substrates. After a collection time of 4 minutes at an RF coil distance of 25cm, both the amorphous SiO 2 nanoparticles and the crystalline structures demonstrated dendritic nucleation of U-oxide species as identified by SEM/EDS.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The WEST tokamak has recently been prepared for long pulse operation with a water-cooled full first wall. Heating is provided by radiofrequency systems, including Lower Hybrid Current Drive (LHCD). The Princeton Plasma Physics Laboratory has developed a multi-energy hard x-ray camera for profile measurements at WEST. The camera is based on a new generation of 2D pixel array detectors that allow the energy threshold to be independently set for each pixel. The diagnostic will provide spatial, temporal and energy resolved measurements of the hard x-ray emission from the full plasma cross-section, investigating several physical quantities such as the electron temperature from continuum emission, the fast electron tail density produced by radiofrequency current drive and runaway electrons, as well as characteristic tungsten x-rays due to beam-target emission at the edge. Overall, this work describes the engineering challenges that the WEST long pulse scenario poses for this diagnostic and how the design addresses and solves them. Vacuum, thermal stress and heat transfer calculations are presented and discussed.
The Large Helical Device (LHD), the largest superconducting helical system in the world, is equipped with advanced heating and diagnostic tools, facilitating plasma control and physics research. Data assimilation was employed for electron temperature control using a real-time Thomson scattering system and real time prediction code. A virtual LHD environment enabled visualization of escaping high-energy tritium ions and demonstrated that these ions impact the rear side of the divertor plate. Pioneering results crucial to plasma control have also been achieved. Real-time wall conditioning using Lithium granule dropping improved bulk ion energy and particle transport while simultaneously enhancing the heavy impurity transport. Progress has also been made in the investigation of turbulence-driven transport. At the confinement bifurcation, ion-scale turbulence decreased, while electron-scale turbulence increased. A change in the anisotropy of turbulent eddies was also observed at the confinement bifurcation. Coexistence of local and non-local turbulence was identified in electron-scale turbulence. Non-local turbulence exhibited the rapid spatial propagation of perturbations throughout the plasma, while local turbulence followed the temperature gradient. A transition between drift-wave turbulence and magnetohydrodynamics (MHD) turbulence was observed with the turbulence minimized at the transition condition. Machine learning analysis was employed to evaluate the temperate and density conditions of this turbulence transition. Then, real-time control of fueling and heating was applied to maintain the turbulence transition condition, improving the energy confinement enhancement factor by 20%. In addition, evidence was obtained for collisionless ion heating by energetic-ion-driven geodesic acoustic modes and MHD bursts. These achievements represent unique contributions to the development of fusion reactors.
Understanding complex physical systems often requires integrating data from multiple diagnostics, each with limited resolution or coverage. We present a machine learning framework that reconstructs synthetic high-temporal-resolution data for a target diagnostic using information from other diagnostics, without direct target measurements during the inference. This multimodal super-resolution technique improves diagnostic robustness and enables monitoring even in case of measurement failures or degradation. Applied to fusion plasmas, our method targets edge-localized modes (ELMs), which can damage plasma-facing materials. By reconstructing super-resolution Thomson Scattering data from complementary diagnostics, we uncover fine-scale plasma dynamics and validate the role of resonant magnetic perturbations (RMPs) in ELM suppression through magnetic island formation. The approach provides new observation supporting the plasma profile flattening due to these islands. Our results demonstrate the framework’s ability to generate high-fidelity synthetic diagnostics, offering a powerful tool for ELM control development in future reactors like ITER. The approach is broadly transferable to other domains facing sparse, incomplete, or degraded diagnostic data, opening new avenues for discovery.
Abstract Diagramming pathways of dimensionless power is a potent method for extrapolating between operating points on present-day tokamaks and future burning plasma devices. The heat transport power, current drive power, H-mode threshold power and other ‘plasma physics’ powers can be expressed in dimensionally correct (or normalized) form as Pa 3/4 , where P is the power and a is the plasma minor radius, with the relative gyroradius ( ρ * ) dependence ranging from gyro-Bohm-like for transport ( P a 3 / 4 ∝ ρ * − 3 / 2 ), Bohm-like for current drive ( P a 3 / 4 ∝ ρ * − 5 / 2 ) and worse than Bohm-like for H-mode threshold ( P a 3 / 4 ∝ ρ * − 3 ). The D–T fusion power cannot be normalized in the same fashion since it is governed by nuclear physics, but at fixed B T it scales like P a 3 / 4 ∝ ρ * − 9 / 2 . Other ‘mixed physics’ powers can be incorporated into the dimensionless power framework by holding B T fixed in the same manner. Diagramming these dimensionless powers vs ρ * shows how the pathway to a steady-state reactor can be optimized relative to various operational boundaries. Using a steady-state hybrid discharge with β N = 3.2 from DIII-D as the starting point, a multi-parameter optimization finds an attractive pathway to steady-state operation on ITER using 76 MW of current drive power (fusion gain of Q fus = 8), along with a pathway to Q fus = 20 in a JET-sized steady-state reactor with B T = 10 T.
A comprehensive suite of innovative diagnostic instrumentation and techniques has been developed by the Johns Hopkins collaboration for measurements of transport and MHD activity in the core and edge of the National Spherical Torus Experiment (NSTX) at the Princeton Plasma Physics Laboratory (PPPL). These include ultrafast multi-energy poloidal soft X-ray arrays, high spatial resolution multi-energy tangential soft X-ray arrays, high spectral resolution transmission grating imaging VUV spectrometers and ultrafast bolometric arrays, which measure over a broad range of spatial, spectral, and temporal scales the emission from ohmic, RF, and beam driven NSTX plasmas. The instrumentation was designed to fulfill the specific measurement needs of the Spherical Torus (ST) magnetic fusion concept and is accompanied by a comprehensive set of numerical tools for data analysis and modeling, aimed at determining as many as possible plasma parameters and at extracting as much as possible physics information from NSTX experiments. The Johns Hopkins diagnostics have been and continue to be instrumental for essentially all the NSTX experiments. In addition, our own research and experiments using these diagnostics has led over the years to fundamental conclusions for the ST fusion physics, such as the establishment of near neoclassical impurity transport in the NSTX core and edge plasma, and the critical role of the Global Alfven Eigenmode (GAE) MHD activity for the ST core electron thermal transport. The Johns Hopkins Diagnostics continue to be operational on the upgraded NSTX-U and will provide key plasma information in its upcoming high field, current and power experimental campaigns. The present report summarizes the main technical characteristics and measurement capabilities of the diagnostic instrumentation and methods developed within our research and describes a key physics result achieved during this research. Success in our research was made possible also by the excellent support and cooperation of the NSTX research team and of PPPL researchers, who are warmly acknowledged here.
A major goal of this project is to significantly increase our understanding of isochoric heating of matter using laser produced proton beams, and the associated high energy density (HED) and warm dense matter (WDM) regimes generated. This will benefit research fields such as planetary science, fusion energy, plasma physics, and material science. For example, it will enhance our understanding of WDM properties of iron and silica under conditions encountered in planetary interiors and diagnostic components in fusion devices exposed to high fluxes of energetic plasma ions. The project is motivated by recent experiments that irradiated Si targets with proton beams generated by the 20 TW-laser at the SLAC MEC end-station. The HED/WDM states are probed using the 50 fs hard X-rays available in the 3rd harmonic of the LCLS. As part of this project, results from the phase contrast X-ray imaging, which shows the generation of compression waves that produces rear surface spallation, are compared with results from the 3D multi-physics multi- material code, PISALE, that combines Arbitrary Lagrangian-Eulerian (ALE) hydrodynamics with Adaptive Mesh Refinement (AMR). This comparison required modifications to several physics models in the PISALE (Pacific Island Structured-AMR with ALE) code. An important aspect of this project is the continued training of graduate students in HED physics and in conducting complex multiphysics simulations.
The current manuscript is a final report on the activities carried out under the Project LDRD-CIS #226834. In scientific terms, the work reported in this manuscript is a continuation of the efforts started with Project LDRD-express #223796 with final report of activities SAND2021-11481, see [83]. In this section we briefly explain what pre-existing developments motivated the current body of work and provide an overview of the activities developed with the funds provided. The overarching goal of the current project LDRD-CIS #226834 and the previous project LDRD-express #223796 is the development of numerical methods with mathematically guaranteed properties in order to solve the Euler-Maxwell system of plasma physics and generalizations thereof. Even though Project #223796 laid out general foundations of space and time discretization of Euler-Maxwell system, overall, it was focused on the development of numerical schemes for purely electrostatic fluid-plasma models. In particular, the project developed a family of schemes with mathematically guaranteed robustness in order to solve the Euler-Poisson model. This model is an asymptotic limit where only electrostatic response of the plasma is considered. Its primary feature is the presence of a non-local force, the electrostatic force, which introduces effects with infinite speed propagation into the problem. Even though instantaneous propagation of perturbations may be considered nonphysical, there are plenty of physical regimes of technical interest where such an approximation is perfectly valid.
Fully relativistic particle-in-cell (PIC) simulations are crucial for advancing our knowledge of plasma physics. Modern supercomputers based on graphics processing units (GPUs) offer the potential to perform PIC simulations of unprecedented scale, but require robust and feature-rich codes that can fully leverage their computational resources. In this work, this demand is addressed by adding GPU acceleration to the PIC code Osiris. An overview of the algorithm, which features a CUDA extension to the underlying Fortran architecture, is given. Detailed performance benchmarks for thermal plasmas are presented, which demonstrate excellent weak scaling on NERSC's Perlmutter supercomputer and high levels of absolute performance. The robustness of the code to model a variety of physical systems is demonstrated via simulations of Weibel filamentation and laser-wakefield acceleration run with dynamic load balancing. Finally, measurements and analysis of energy consumption are provided that indicate that the GPU algorithm is up to ~14 times faster and ~7 times more energy efficient than the optimized CPU algorithm on a node-to-node basis. The described development addresses the PIC simulation community's computational demands both by contributing a robust and performant GPU-accelerated PIC code and by providing insight into efficient use of GPU hardware.
Significance The problem of solar wind heating is one of the long-standing fundamental problems of space plasma physics. We propose that the electron component of the solar wind plasma is heated by electrons streaming from the hot solar corona and slowly losing their energy due to weak Coulomb collisions. Our first-principle kinetic derivation predicts the scaling of the solar wind electron temperature with the heliocentric distance, which is in good agreement with observations. Our theory reveals the role of the non-Maxwellian electron distribution function in the heating mechanism, and it may be valuable for interpreting the results of in situ satellite measurements.
Both scientific and technical innovation is needed for the realization of an attractive engineering solution for a timely and cost-effective Pilot Plant, the design and construction of which is the overarching recommendation of the FESAC Long Range Plan, and the 2021 NASEM Pilot Plant reports, which underpin the Bold Decadal Vision. The two most significant plasma physics gaps to close for a Compact Pilot Plant (CPP) are core confinement improvement and heat flux mitigation, neither of which have been closed in an integrated fashion for any planned fusion power production device. High core confinement and stability are essential for producing majority self-driven plasmas in CPPs with reduced size and auxiliary heating power requirements, with an improvement in confinement being the major driver for cost reduction of a CPP. The National Spherical Tokamak Experiment - Upgrade (NSTX-U) is a unique low aspect ratio research facility that will address the fundamental challenge of developing the science and technology basis for a CPP design that integrates high core and edge confinement with the ability to mitigate very high incident heat fluxes. NSTX-U capabilities will enable the high performance, already achieved on NSTX, to extend into physics regimes much closer to those anticipated in Spherical Tokamak (ST)-based CPPs. These confinement and stability properties will be assessed by a full complement of diagnostics and analysis tools, which will also aid in the development of the underlying theory and predictive models needed for further optimization. Both conventional and transformative heat flux mitigation methods, such as liquid lithium plasma-facing components, will be developed and tested in-situ in NSTX-U at incident heat fluxes of ~100 MW/m 2 , and will inform plans and reduce risk for a subsequent major upgrade to the device to fully heated, high-Z wall and full liquid lithium divertor capability, a technology that potentially could then be implemented on any magnetic confinement device at any aspect ratio. NSTX-U research is fully complementary to programs performed on other STs, nationally and internationally. Furthermore, NSTX-U research has a direct connection to the private sector by informing design choices for future power production facilities being developed by these companies. The NSTX-U program will operate as a national User Facility, with collaborating researchers, engineers, and graduate students from 19 outside institutions, and open to participation and experiments led by researchers from both public and private entities. The research program will advance workforce development through training of young scientists, engineers, and technicians, and it will also serve for further diagnostic innovation, especially for high heat flux and high-Z wall environments, and implementation of advanced artificial intelligence (AI) for plasma and heat flux control.
Gaia has revealed a variety of substructures in the phase space of stars in the Solar neighbourhood, including the vertical ‘Snail’ $(z, v_z)$ in space. Such substructures are often interpreted as the incompletely phase-mixed response of the disc stars to a single perturbation, such as an impulsive encounter with a satellite galaxy. In this paper, we consider the possibility that such structures contain manifestations of phase-space echoes. First established in plasma physics in the 1960s, echoes arise when a collisionless system is perturbed twice: the macroscopic responses to both perturbations mix to small scales in phase space, whereupon they couple non-linearly, producing a third macroscopic ‘echo’ response without the need for a third perturbation. We derive the galactic analogue of the plasma echo theory using angle-action variables and apply it to a one-dimensional model of vertical motion in the Milky Way. We verify the predicted echo behaviour using idealized test particle simulations, both with and without the inclusion of diffusion through orbital scattering off molecular clouds. While we conclude that the Gaia Snail itself is unlikely a (pure) echo effect, the basic physics we uncover is sufficiently generic that we expect phase-space echoes to be common in disc galaxies.
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.
Unstructured mesh particle-in-cell, PIC, simulations executing on the current and next generation of massively parallel systems require new methods for both the mesh and particles to achieve performance and scalability on GPUs. The traditional approach to implementing PIC simulations defines data structures and algorithms in terms of particles with a full copy of the unstructured mesh on every process. To effectively scale the unstructured mesh and particles, mesh-based PIC uses the unstructured mesh as the predominant data structure with the particles stored in terms of the mesh entities. Here, this paper details the PUMIPic library, a framework for developing efficient and performance-portable mesh-based PIC simulations on GPU systems. A pseudo physics simulation based on a five-dimensional gyro-kinetic code for modeling plasma physics is used to examine the performance of PUMIPic. Scaling studies of the unstructured mesh partition and number of particles are performed up to 4096 nodes of the Summit system at Oak Ridge National Laboratory. The studies show that mesh-based PIC can utilize a partitioned mesh and maintain scaling up to system limitations.