Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fusion Simulation Project”

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 343 records · Page 19

Characterization of an electrothermal plasma source for fusion transient simulations

The realization of fusion energy requires materials that can withstand high heat and particle fluxes at the plasma material interface. Here in this work, an electrothermal (ET) plasma source has been designed as a transient heat flux source for a linear plasma material interaction device. An ET plasma source operates in the ablative arc regime driven by a DC capacitive discharge. The current channel width is defined by the 4 mm bore of a boron nitride liner. At large plasma currents, the arc impacts the liner wall, leading to high particle and heat fluxes to the liner material, which subsequently ablates and ionizes. This results in a high density plasma with a large unidirectional bulk flow out of the source exit. The pulse length for the ET source has been optimized using a pulse forming network to have durations of 1 and 2 ms. The peak currents and maximum source energies seen in this system are 1.9 kA and 1.2 kJ for the 2 ms pulse and 3.2 kA and 2.1 kJ for the 1 ms pulse, respectively. This work is a proof of the principal project to show that an ET source produces electron densities and heat fluxes comparable to those anticipated in transient events in large future magnetic confinement fusion devices. Heat flux, plasma temperature, and plasma density were determined for each shot using infrared imaging and optical spectroscopy techniques. This paper will discuss the assumptions, methods, and results of the experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Prostate implant reconstruction from C-arm images with motion-compensated tomosynthesis

Purpose: Accurate localization of prostate implants from several C-arm images is necessary for ultrasound-fluoroscopy fusion and intraoperative dosimetry. The authors propose a computational motion compensation method for tomosynthesis-based reconstruction that enables 3D localization of prostate implants from C-arm images despite C-arm oscillation and sagging. Methods: Five C-arm images are captured by rotating the C-arm around its primary axis, while measuring its rotation angle using a protractor or the C-arm joint encoder. The C-arm images are processed to obtain binary seed-only images from which a volume of interest is reconstructed. The motion compensation algorithm, iteratively, compensates for 2D translational motion of the C-arm by maximizing the number of voxels that project on a seed projection in all of the images. This obviates the need for C-arm full pose tracking traditionally implemented using radio-opaque fiducials or external trackers. The proposed reconstruction method is tested in simulations, in a phantom study and on ten patient data sets. Results: In a phantom implanted with 136 dummy seeds, the seed detection rate was 100% with a localization error of 0.86 {+-} 0.44 mm (Mean {+-} STD) compared to CT. For patient data sets, a detection rate of 99.5% was achieved in approximately 1 min per patient. The reconstruction results for patient data sets were compared against an available matching-based reconstruction method and showed relative localization difference of 0.5 {+-} 0.4 mm. Conclusions: The motion compensation method can successfully compensate for large C-arm motion without using radio-opaque fiducial or external trackers. Considering the efficacy of the algorithm, its successful reconstruction rate and low computational burden, the algorithm is feasible for clinical use.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

On the specificity between mapping of initial and final states of the magneto Rayleigh-Taylor instability

In this LDRD we investigated the application of machine learning methods to understand dimensionality reduction and evolution of the Rayleigh-Taylor instability (RTI). As part of the project, we undertook a significant literature review to understand current analytical theory and machine learning based methods to treat evolution of this instability. We note that we chose to refocus on assessing the hydrodynamic RTI as opposed to the magneto-Rayleigh-Taylor instability originally proposed. This choice enabled utilizing a wealth of analytic test cases and working with relatively fast running open-source simulations of single-mode RTI. This greatly facilitated external collaboration with URA summer fellowship student, Theodore Broeren. In this project we studied the application of methods from dynamical systems learning and traditional regression methods to recover behavior of RTI ranging from the fully nonlinear to weakly nonlinear (wNL) regimes. Here we report on two of the tested methods SINDy and a more traditional regression-based approach inspired by analytic wNL theory with which we had the most success. We conclude with a discussion of potential future extensions to this work that may improve our understanding from both theoretical and phenomenological perspectives.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Edge-localized-modes in tokamaks

Edge-localized-modes (ELMs) are a ubiquitous feature of H-mode in tokamaks. When gradients in the H-mode transport barrier grow to exceed the MHD stability limit the ELM instability grows explosively, rapidly transporting energy and particles onto open field lines and material surfaces. Though ELMs provide additional particle and impurity transport through the H-mode transport barrier, enabling steady operation, the resulting heat flux transients to plasma facing surfaces project to large amplitude in future low collisionality burning plasma tokamaks. Measurements of the ELM heat flux deposition onto material surfaces in the divertor and main chamber indicate significant broadening compared to inter-ELM heat flux, with a timescale for energy deposition that is consistent with sonic ion flow and numerical simulation. Comprehensive ELM simulation is highlighting the important physics processes of ELM transport including parallel transport due to magnetic reconnection and turbulence resulting from collapse of the H-mode transport barrier. Encouraging prospects for ELM control and/or suppression in future tokamaks include intrinsic modes of ELM free operation, ELM triggering with frequent small pellet injection and the application of 3D magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Framework for International Collaboration on ITER Using Large-Scale Data Transfer to Enable Near-Real-Time Analysis

The global nature of the ITER project along with its projected ~ petabyte per day data generation presents a unique challenge, but also an opportunity for the fusion community to rethink, optimize, and enhance our scientific discovery process. Recognizing this, collaborative research with computational scientists was undertaken over the past several years to create a framework for large-scale data movement across wide-area networks (WANs), to enable global near-real time analysis of fusion data. This would broaden the available computational resources for analysis/simulation, and increase the number of researchers actively participating in experiments. An official demonstration of this framework for fast, large data transfer and real-time analysis was carried out between the KSTAR tokamak in Daejeon, Korea and PPPL in Princeton, USA. Streaming large data transfer, with near real-time movie creation and analysis of the KSTAR Electron Cyclotron Emission imaging (ECEI) data, was performed using the I/O framework ADIOS, and comparisons made at PPPL with simulation results from the XGC1 code. These demonstrations were made possible utilizing an optimized network configuration at PPPL, which achieved over 8.8 Gbps (88% utilization) in throughput tests from NFRI to PPPL. This demonstration showed the feasibility for large-scale data analysis of KSTAR data, and provides a nascent framework to enable use of globally distributed computational and personnel resources in pursuit of scientific knowledge from the ITER experiment.

43 PARTICLE ACCELERATORS↗

Relativistically Induced Transparency in Plasma (Final Research Report)

The major goal of this project was to study the process of relativistic transparency in plasmas. This process occurs when a very high intensity laser pulse interacts with a plasma. The electric field of the laser oscillates the electron speeds approaching the speed of light, increasing the effective mass of the electrons by the time-averaged Lorentz factor <γ>. The increase in the effective mass alters the critical density, the density at which the plasma becomes opaque to the electromagnetic wave of the laser. This project used experiments and particle-in-cell simulations to study this process in ultra-thin foil interactions. Experiments were performed through the LaserNetUS program on the Ohio State University laser system Scarlet. The transmitted and reflected light was studied as a function of the target thickness, the near-field profiles, the total energy and pulse duration were measured. Particle-in-cell simulations, using the OSIRIS code were used to model the interaction and better understand the dynamics. The results of both show current relativistic transparency theory is insufficient to predict the results, this is likely because the theories do not account for the laser energy absorption into the target and global plasma fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) - NIMROD DIII-D Shattered Pellet Injection Disruption Mitigation and RE Modeling. Final report

Final technical report for UCSD and subcontractor contributions to DE-SC0016420, Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM). Project description. MHD modeling includeing runaway electron drift orbit calculations will be carried out for DII-D disruption mitigation scenarios. The effects of pitch angle scattering on the energy distribution of runaway electron population will be evaluated using one or more collision models integrated / coupled with the runaway orbit following model in the nimrod code. An extension of the test particle model will be developed to predict runaway electron generation considering both primary and secondary sources of runaway electron production. Simulations simultaneously modelling the effects of generation, confinement, and moderation-space losses during a DII-D rapid shutdown (mitigated disruption) will be completed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collaborative Research: The Evolution of Magnetic Skeletons During 3D Reconnection. Final report

We investigated how the plasmoid instability affects the evolution of topological features of a magnetic field and the process of magnetic reconnection with three sub-projects. (i) We extensively studied the criterion for plasmoid-mediated current sheet disruption and the onset of fast reconnection. Scaling relations with respect to plasma parameters have been established, using both direct numerical simulations and analytic theories, for the disruption time, the current sheet width, the instability linear growth rate, and the dominant mode wavenumber. (ii) We applied the theories of plasmoid-mediated current sheet disruption to magnetohydrodynamic turbulence and predicted a steepening of the turbulence energy spectrum. This steepening has been confirmed by direct numerical simulations. (iii) We performed numerical simulations of ultraviolet burst events in the solar atmosphere with a spine-dome magnetic topology. The plasmoid instability leads to unsteady reconnection, which is consistent with solar observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗

Gyrokinetic Particle Simulation of Compressible Electromagnetic Turbulence in High-β Plasmas

Supported by this award, the PI and his research group at the University of California, Irvine (UCI) have carried out computational and theoretical studies of instability, turbulence, and transport in laboratory and space plasmas. Several massively parallel, gyrokinetic particle simulation codes have been developed to study electromagnetic turbulence in space and laboratory plasmas. In space plasma projects, the simulation codes have been successfully applied to study the spectral cascade and plasma heating in kinetic Alfven wave turbulence, the linear and nonlinear properties of compressible modes including mirror instability and drift compressional mode, and the stability of the current sheet instabilities with finite guide field in the context of collisionless magnetic reconnection. The research results have been published in 25 journal papers and presented at many national and international conferences. Reprints of publications, source codes, and other research-related information are also available to general public on the PI’s webpage (http://phoenix.ps.uci.edu/zlin/). Two PhD theses in space plasma physics are highlighted in this report.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Potential design problems for ITER fusion device

The international thermonuclear experimental reactor (ITER) is a worldwide project currently being built in France for the demonstration of the feasibility of thermonuclear technologies for future realization of successful commercial fusion energy. ITER is of the tokamak based design using strong magnetic fields to confine the very hot plasma needed to induce the fusion reaction. Tokamak devices are currently the front leading designs. Building a successful magnetic fusion device for energy production is of great challenge. A key obstacle to such design is the performance during abnormal events including plasma disruptions and so-called edge-localized modes (ELMs). In these events, a massive and sudden release of energy occurs quickly, due to loss of full or partial plasma confinement, leading to very high transient power loads on the reactor surface boundaries. A successful reactor design should tolerate several of these transient events without serious damages such as melting and vaporization of the structure. This paper highlights, through comprehensive state-of-the-art computer simulation of the entire ITER interior design during such transient events, e.g., ELMs occurring at normal operation and disruptions during abnormal operation, potential serious problems with current plasma facing components (PFCs) design. The HEIGHTS computer package is used in these simulations. The ITER reactor design was simulated in full and exact 3D geometry including all known relevant physical processes involved during these transient events. The current ITER divertor design may not work properly and may requires significant modifications or new innovative design to prevent serious damage and to ensure successful operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gyrokinetic Particle Simulation of Turbulent Transport in Burning Plasmas

The SciDAC project at the IFS advanced the state of high performance computing for turbulent structures and turbulent transport. The team project with Prof Zhihong Lin [PI] at Univ California Irvine produced new understanding of the turbulent electron transport. The simulations were performed at the Texas Advanced Computer Center TACC and the NERSC facility by Wendell Horton, Lee Leonard and the IFS Graduate Students working in that group. The research included a Validation of the electron turbulent transport code using the data from a steady state university experiment at the University of Columbia in which detailed probe measurements of the turbulence in steady state were used for wide range of temperature gradients to compare with the simulation data. These results were published in a joint paper with Texas graduate student Dr. Xiangrong Fu using the work in his PhD dissertation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma Interactions with Mixed Materials and Impurity Transport

The project brings together three discipline areas at LLNL to develop advanced capability to predict the impact of plasma/material interactions (PMI) on metallic surfaces in magnetic fusion energy (MFE) devices. These areas are (1) modeling transport of wall impurity ions through the edge plasma to the core plasma, (2) construction of a laser blow-off (LBO) system for injecting precise amounts of metallic atoms into a tokamak plasma, and (3) material science analysis of fundamental processes that modify metallic surfaces during plasma bombardment. The focus is on tungsten (W), which is being used for the ITER divertor and in designs of future MFE devices. In area (1), we have worked with the University of California, San Diego (UCSD) on applications of the UEDGE/DUSTT coupled codes to predict the influx of impurity ions from W dust through the edge plasma, including periodic edge-plasma oscillations, and revived a parallel version of UEDGE to speed up these simulations. In addition, the impurity transport model in the 2D UEDGE code has been implemented into the 3D BOUT++ turbulence/transport code to allow fundamental analysis of the impact of strong plasma turbulence on the impurity transport. In area (2), construction and testing of the LBO injection system has been completed. The original plan to install the LBO on the National Spherical Torus Experiment Upgrade (NSTX-U) at Princeton and its use to validate the impurity transport simulations is delayed owing to NSTX-U being offline for substantial magnetic coil repair period. In area (3), an analytic model has been developed to explain the growth of W tendrils (or fuzz) observed for helium-containing plasmas. Molecular dynamics calculations of W sputtering by W and deuterium (D) ions shows that a spatial blending of interatomic potentials is needed to describe the near-surface and deeper regions of the material.

36 MATERIALS SCIENCE↗

Development of an accelerator-based diagnostic for plasma-facing surfaces in magnetic confinement devices

The focus of this work was the development of the Accelerator-based In-situ Material Surveillance (AIMS) diagnostic. The understanding of plasma-material interaction processes in fusion devices has been greatly hampered by the lack of ability to measure the surface ~micron layers of solid surfaces that are in contact with and greatly modified by the boundary fusion plasma. AIMS was the first diagnostic developed and deployed that allowed non-destructive surface measurements with adequate spatial and time resolution. The deployment was primarily carried out on the Alcator C-Mod tokamak. AIMS diagnostic goal is to remotely generate isotopic concentration maps on a plasma shot-to-shot timescale that cover a large fraction of the plasma-facing surface inside of a magnetic fusion device without the need for vacuum breaks or physical access to the material surfaces. The diagnostic used a compact (~ 1 m), high-current (~ 1 milliamp) radio-frequency quadrupole accelerator to inject 0.9 MeV deuterons into the Alcator C-Mod tokamak. The tokamak magnetic fields were controlled – in between plasma shots – to steer the deuterons to material surfaces where the deuterons cause high-Q nuclear reactions with low-Z isotopes up to 5 microns into the material. The induced neutrons and gamma rays are measured with scintillation detectors; energy spectra analysis provides quantitative reconstruction of surface compositions. Experimental validation showed that low-Z isotopes such as deuterium and boron can be quantified on the material surfaces, and that magnetic steering provides access to different measurement locations. AIMS measurements of D retention on inner wall PFCs were acquired during diverted and limited plasma operations and during wall conditioning experiments om Alcator C-Mod. Intershot measurements demonstrate the local erosion and co-deposition of boron films on PFC surfaces with a constant D-B ratio. This is consistent with previous results suggesting that D co-=deposition with boron is insufficient to account for the net retention observed in Alcator C-Mod. Changes in deuterium concentration during boronization, electron cyclotron and glow cleanings were also measured. These initial results were limited to studying low-Z surface properties, because the Coulomb barrier precludes nuclear reactions between high-Z elements and the ~1 MeV AIMS deuteron beam. In order to measure the high-Z erosion of interest in Alcator C-Mod and future fusion devices, a new technique using deuteron-induced gamma emission and a low-Z depth marker was developed using implanted depth markers as a method of measuring bulk erosion and redeposition. This method is dubbed DEA: Depth marker Erosion with AIMS. DEA simulations and experiments were used to validate this technique, both in the laboratory and via experiments in the EAST tokamak. Based on these results, the DEA technique, using ion beam implanted depth markers as a reference to the surface in order to monitor surfaces for erosion and redeposition, is a viable option as a diagnostic, both in and ex situ. These markers show stability to temperature, the tokamak environment, and multiple probing beams. The ex-situ (eDEA) method, which uses resonances in the cross sections of gamma producing reactions to track the surface, was used to successfully measure erosion in samples that were exposed in the EAST tokamak. The in-situ (iDEA) method, which uses the gamma yield ratio to monitor changes to the surface, was used to differentiate between markers implanted to different depths. Both techniques were successfully simulated, and both presented unique challenges experimentally and in data analysis. In particular, the eDEA probing beam was found to have energy variations leading to initial experimental uncertainty, and the iDEA gamma spectra had conflicting gamma peaks which make data analysis difficult. However, the end result is a promising framework for a diagnostic technique that supplements the abilities of AIMS and expands the possibilities of surface diagnostics in fusion devices. While the results of this work are overwhelmingly encouraging for further development of AIMS in general, and the DEA technique specifically, several areas require significant development before this diagnostic is as robust and widely applicable as a surface diagnostic ought to be. Recommended future work includes the development of a cross section/yield measurement database, erosion tests with secondary verification of results using a fiducial mark from a focused ion beam, and acquiring a new accelerator with a more reliable high energy proton and deuterium beam The project resulted in 14 publications, including the PhD theses of three MIT graduate students, cited in the included bibliography.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SciDAC Center for Simulation of Wave-Plasma Interactions - Iterated Finite-Orbit Monte Carlo Simulations with Full-Wave Fields for Modeling Tokamak ICRF Wave Heating Experiments - Final Report

This final report describes the work performed under U.S. Department of Energy Cooperative Agreement DE-FC02-08ER54954 for the period April 1, 2011 through March 31, 2013. The goal of this project was to perform iterated finite-orbit Monte Carlo simulations with full-wall fields for modeling tokamak ICRF wave heating experiments. In year 1, the finite-orbit Monte-Carlo code ORBIT-RF and its iteration algorithms with the full-wave code AORSA were improved to enable systematical study of the factors responsible for the discrepancy in the simulated and the measured fast-ion FIDA signals in the DIII-D and NSTX ICRF fast-wave (FW) experiments. In year 2, ORBIT-RF was coupled to the TORIC full-wave code for a comparative study of ORBIT-RF/TORIC and ORBIT-RF/AORSA results in FW experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detection of Coherent Structures in Extreme-Scale Simulations

The analysis of coherent structures is a common problem in many scientific domains ranging from astrophysics to combustion, fusion, and materials science. The data from three-dimensional simulations are analyzed to detect the structures, extract statistics on them, and track them over time to gain insights into the phenomenon being modeled. This analysis is typically done off-line, using data that have been written out by the simulations. However, the move towards extreme scale architectures, with multi-core processors and graphical processing units, will affect how such analysis is done as it is unlikely that the systems will support the I/O bandwidth required for off-line analysis. Moving the analysis in-situ is a solution only if we know a priori what analysis will be done, as well as the algorithms used and their parameter settings. Even then, we need to ensure that this will not substantially increase the memory requirements or the data movement as the former will be limited and the latter will be expensive. In the Exa-DM project, a collaboration between Lawrence Livermore National Laboratory and University of Minnesota, we are exploring ways in which we can address the conflicting demands of coherent structure analysis of simulation data and the architecture of modern parallel systems, while enabling scientific discovery at the exascale. In this paper, we describe our work in two areas: the in situ implementation of an existing algorithm for coherent structure analysis and the use of graph-based techniques to efficiently compress the data.

79 ASTRONOMY AND ASTROPHYSICS↗

Cold Plasma Measurements

We have continued the simulation campaign in support of our ongoing magnetospheric cold plasma research project. This project aims to develop the next-generation particle instruments to measure the properties of the cold particle populations in the Earth’s magnetosphere. For this purpose, simulations have been performed with a Particle-In-Cell (PIC) code called the Curvilinear PIC (CPIC). The code is formulated in curvilinear geometry and couples the standard PIC algorithm with algorithms for the generation and adaptation of the underlaying computational mesh. It conforms to complex objects like spacecraft and it can place more grid points in regions where higher resolution is needed. The code also features a scalable solver based on the multigrid algorithm and it is fully parallelized via domain decomposition and MPI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimization of high heat flux components for DIII-D neutral beam upgrades

Upgrade of the DIII-D neutral beams leads to enhanced heat loads on many components, such as calorimeter, collimator, and pole shields which protect neutral beam magnets. Power increase from 2.6 MW to 3.2 MW per source leads to a normal heat flux loads of up to 55 MW/m2 for the calorimeter. The Princeton Plasma Physics Laboratory is responsible for the design and manufacturing of the upgrades of these components. Heat flux distribution on neutral beam components is very uneven and leads to significant thermal stresses. High heat flux density impact requires surface optimization to reduce surface heat flux projection, and avoid localized melting. Several new design features were introduced to accommodate increased heat loads, such as molybdenum inserts for the pole shields, two-dimensional shaping for the calorimeter, and three-dimensional shape optimization and replaceable copper inserts for the collimator. Also, all three components include an optimized cooling system design featuring peripheral cooling of copper components. The optimization process included applying analytical relations for the transient temperature distributions on the high heat flux components. These relations were confirmed by previous DIII-D experimental results. To confirm the designs, numerical simulations were performed. Results of the design optimization and numerical simulations will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗