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 469 records · Page 26

Advanced Accelerators: Particle, Photon and Plasma Wave Interactions

The overall objective of this project was to study the acceleration of electrons to very high energies over very short distances based on trapping slowly moving electrons in the fast moving potential wells of large amplitude plasma waves, which have relativistic phase velocities. These relativistic plasma waves, or wakefields, are the basis of table-top accelerators that have been shown to accelerate electrons to the same high energies as kilometer-length linear particle colliders operating using traditional decades-old acceleration techniques. The accelerating electrostatic fields of the relativistic plasma wave accelerators can be as large as GigaVolts/meter, and our goal was to study techniques for remotely measuring these large fields by injecting low energy probe electron beams across the plasma wave and measuring the beam’s deflection. Our method of study was via computer simulations, and these results suggested that the deflection of the probe electron beam was directly proportional to the amplitude of the plasma wave. This is the basis of a proposed diagnostic technique, and numerous studies were performed to determine the effects of changing the electron beam, plasma wave and laser beam parameters. Further simulation studies included copropagating laser beams with the relativistic plasma waves. New interesting results came out of these studies including the prediction that very small scale electron beam bunching occurs, and an anomalous line focusing of the electron beam occurs under certain conditions. These studies were summarized in the dissertation of a graduate student who obtained the Ph.D. in physics. This past research program has motivated ideas for further research to corroborate these results using particle-in-cell simulation tools which will help design a test-of-concept experiment in our laboratory and a scaled up version for testing at a major wakefield accelerator facility.

43 PARTICLE ACCELERATORS↗

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗

Towards Predictive Plasma Science and Engineering through Revolutionary Multi-Scale Algorithms and Models (Final Report)

This report describes the high-level accomplishments from the Plasma Science and Engineering Grand Challenge LDRD at Sandia National Laboratories. The Laboratory has a need to demonstrate predictive capabilities to model plasma phenomena in order to rapidly accelerate engineering development in several mission areas. The purpose of this Grand Challenge LDRD was to advance the fundamental models, methods, and algorithms along with supporting electrode science foundation to enable a revolutionary shift towards predictive plasma engineering design principles. This project integrated the SNL knowledge base in computer science, plasma physics, materials science, applied mathematics, and relevant application engineering to establish new cross-laboratory collaborations on these topics. As an initial exemplar, this project focused efforts on improving multi-scale modeling capabilities that are utilized to predict the electrical power delivery on large-scale pulsed power accelerators. Specifically, this LDRD was structured into three primary research thrusts that, when integrated, enable complex simulations of these devices: (1) the exploration of multi-scale models describing the desorption of contaminants from pulsed power electrodes, (2) the development of improved algorithms and code technologies to treat the multi-physics phenomena required to predict device performance, and (3) the creation of a rigorous verification and validation infrastructure to evaluate the codes and models across a range of challenge problems. These components were integrated into initial demonstrations of the largest simulations of multi-level vacuum power flow completed to-date, executed on the leading HPC computing machines available in the NNSA complex today. These preliminary studies indicate relevant pulsed power engineering design simulations can now be completed in (of order) several days, a significant improvement over pre-LDRD levels of performance.

42 ENGINEERING↗

Numerical Assessment of the Diagnostic Capabilities of the Instrumented Calorimeter for SPIDER (STRIKE)

An important feature of the ITER project is represented by additional heating via injection of neutral beams from accelerated negative ions. To study and optimise their production, the SPIDER test facility is under construction in Padova, with the aim of testing beam characteristics and to verify the source proper operation.STRIKE (Short-Time Retractable Instrumented Kalorimeter Experiment) is a diagnostic to characterise the SPIDER negative ion beam during short operation (several seconds). During long pulse operations, STRIKE is parked off-beam in the vacuum vessel. The most important measurements are beam uniformity, beamlet divergence and stripping losses. STRIKE is directly exposed to the beam and is formed of 16 tiles, one for each beamlet groups. The measurements are provided by thermal cameras, current sensors, thermocouples and electrostatic sensors. This paper presents the investigation of the influence on the response of STRIKE of: thermal characteristics of the tile material, exposure angle, features of some dedicated diagnostics. The uniformity of the beam will be studied by measurements of the current flowing through each tile and by thermal cameras. Simulations show that it will be possible to verify experimentally whether the beam meets the ITER requirement about the maximum allowed beam non-uniformity (below {+-}10%). In the simulations also the influence of the beam halo has been included; the effect of off-perveance conditions has been studied. To estimate the beamlet divergence, STRIKE can be moved along the beam direction at two different distances from the accelerator. The optimal positions have been defined taking into account design constraints. The effect of stripping on the comparison between currents and heat loads has been assessed; this will allow to obtain an experimental estimate of stripping. Electrostatic simulations have provided the suitable tile biasing voltage in order to reabsorb secondary particles into the same tile as the one where they were emitted from.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

From Interstellar Cloud to Star to Laboratory: Frontier HEDP Studies of Magnetized Colliding Plasma Flows with Strong Radiative Cooling

The project investigated a number of phenomena related to the interaction of highly supersonic, strongly cooled flows where magnetic field and the development of plasma instabilities play a role in the system evolution. Colliding plasma jets and flow interactions are commonplace in astrophysical systems which are strongly collisional, and radiatively cooled despite their large size. Many observations require the testing of hypotheses which can lead to a greater understanding of the dominant plasma condition and processes. We have developed experiments to investigate counter-propagating, high Mach number flows where the radiative loss was rate varied. The development of extended shock structure at the intersection of these flows was observed and characterized over a range of parameter space. Of particular note, was the fragmentation of shocks in collisional, radiatively cooled flows which show strong indications of the development of a turbulent power density spectrum driven by cooling instabilities. This is significant in that similar processes are believed to develop in stellar nurseries where gravity play little role. We have demonstrated the efficacy of novel experiments where a high Mach number flow interacts with a magnetized target. The target B-field dominates the morphology of the shock evolution, in contract to the purely hydrodynamic case. Additionally, the interaction of 2 bow shocks was shown to result in regular reflection mid-plane, in line with detailed simulation work. This is directly relevant to the often-complicated structures formed at the head of HH jets in which multiple shocks can generate new structures in observations as they interact. Research work from all the above has been published in peer-reviewed journals and presented at international conferences, regularly by graduate and undergraduate students.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Anomalous Retrograde Drifts in Obstructed Magnetron Microdischarges: a Consequence of a Field Reversal in the Anode Sheath?

The goal of this project was to carry out experiments at the Plasma Collaborative Research Facility (Princeton Plasma Physics Laboratory) aimed at confirming or refuting the hypothesis that an electric field reversal occurs in the near-anode region of a strongly obstructed planar magnetron microdischarges. Magnetron discharges are magnetized plasmas with E x B field topologies, and variants of such discharges are widely used in materials processing. The topologies are also similar to that of Hall thrusters used in space propulsion. Understanding the structure of these discharges will aid in developing reliable simulations that are intended to facilitate the further development of such technologies, as the field reversals, presumed to be due to strong axial-driven electron transport, provides further insight into transport properties. The possible existence of such a field reversal is hinted at because of the retrograde-propagating (i.e., negative E x B) spokes or gradient-driven drift instabilities seen in the high speed videos and segmented anode probing of the plasma. The presence of this field reversal should be reflected in the ion dynamics, particularly, in an ion drift towards the anode. The proposed experiments to be conducted at the PCRF are to characterize this anode-directed ion drift using state-of-the-art laser-induced fluorescence velocimetry. In particular, we seek to see if the near-anode ion velocity reverses direction when the external circuit is varied to alter the electron current demand imposed at the anode. Such an ion velocity reversal would confirm the mechanism, which is that the highly obstructed (small gapped) discharge produces strong axial gradients in plasma density, driving more electrons towards the anode then that demanded by the external circuit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dust and Ions: Self Organization and Stability (Final Report)

This project explores the stability and structure of systems with non-reciprocal interactions, challenging the traditional understanding based on Newton's third law, which states that every action has an equal and opposite reaction. Reciprocal forces are fundamental to the stability of systems ranging in size from atomic nuclei to galactic clusters. Our research investigates what happens when the forces between two objects are not equal and opposite. We used dusty plasmas as a model system to study non-reciprocal interactions. In a plasma chamber, micron-sized dust particles acquire a negative charge and form 2D planar "dust crystals" when levitated by the electric field present in the plasma sheath at the interface between the plasma and the lower surface of the chamber. This electric field also drives a vertical ion flow, creating a positively charged "plasma wake" downstream of the dust grains. While horizontally aligned dust grains interact reciprocally, a slight vertical displacement causes non-reciprocal interactions due to the attraction of the lower dust grain to the upper dust grain’s ion wake. Our experiments investigated the range of plasma conditions (gas pressure and system power) where stable dusty plasma structures are able to self-organize, aided by the ion wake. We studied systems ranging from pairs of dust particles to large 2D crystals, providing insights into the conditions that lead to stable or unstable structures. We used numerical simulations to investigate how ion wakes changed in response to changes in the operating conditions as well as how the wakes of separate grains interact when dust grains are in close proximity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Fusion Nuclear Science Pathways Assessment

With the strong commitment of the US to the success of the ITER burning plasma mission, and the project overall, it is prudent to consider how to take the most advantage of this investment. The production of energy from fusion has been a long sought goal, and the subject of several programmatic investigations and time line proposals. The nuclear aspects of fusion research have largely been avoided experimentally for practical reasons, resulting in a strong emphasis on plasma science. Meanwhile, ITER has brought into focus how the interface between the plasma and engineering/technology, presents the most challenging problems for design. In fact, this situation is becoming the rule and no longer the exception. ITER will demonstrate the deposition of 0.5 GW of neutron heating to the blanket, deliver a heat load of 10-20 MW/m 2 or more on the divertor, inject 50-100 MW of heating power to the plasma, all at the expected size scale of a power plant. However, in spite of this, and a number of other technologies relevant power plant, ITER will provide a low neutron exposure compared to the levels expected to a fusion power plant, and will purchase its tritium entirely from world reserves accumulated from decades of CANDU reactor operations. Such a decision for ITER is technically well founded, allowing the use of conventional materials and water coolant, avoiding the thick tritium breeding blankets required for tritium self-sufficiency, and allowing the concentration on burning plasma and plasma-engineering interface issues. The neutron fluence experienced in ITER over its entire lifetime will be ~ 0.3 MW-yr/m 2 , while a fusion power plant is expected to experience 120-180 MW-yr/m 2 over its lifetime. ITER utilizes shielding blanket modules, with no tritium breeding, except in test blanket modules (TBM) located in 3 ports on the midplane, which will provide early tests of the fusion nuclear environment with very low tritium production (a few g per year).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Note on numerical study of the beam energy spread in NDCX-I

The kinetic energy spread (defined here as the standard deviation of the beam particle energies) sets the ultimate theoretical limit on the longitudinal compression that can be attained on NDCX-I and NDCX-II. Experimental measurements will inevitably include the real influences on the longitudinal phase space of the beam due to injector and accelerator field imperfections1. These induced energy variations may be the real limit to the longitudinal compression in an accelerator. We report on a numerical investigation of the energy spread evolution in NDCX-I; these studies do not include all the real imperfections, but rather are intended to confirm that there are no other intrinsic mechanisms (translaminar effects, transverse-longitudinal anisotropy instability, etc.) for significant broadening of the energy distribution. We have performed Warp simulations that use a realistic Marx voltage waveform which was derived from experimental measurements (averaged over several shots), a fully-featured model of the accelerating and focusing lattice, and new diagnostics for computing the local energy spread (and temperature) that properly account for linear correlations that arise from the discrete binning along each physical dimension (these capabilities reproduce and extend those of the earlier HIF code BPIC). The new diagnostics allow for the calculation of multi-dimensional maps of energy spread and temperature in 2-D axisymmetric or 3-D Cartesian space at selected times. The simulated beam-line was terminated at z = 3 m by a conducting plate, so as to approximately reproduce the experimental conditions at the entrance of the spectrometer that was used for mapping the longitudinal phase space. Snapshots of the beam projection and current, as well as the Marx waveform and history of beam kinetic energy collected at the end plate, are shown in Fig. 1. A two-dimensional axisymmetric map of energy spread from simulations of a typical NDCX-I configuration is shown in Fig. 2 (a). The energy spread starts at 0.1 eV at the source and rapidly rises to a few eV, then fluctuates between a fraction of an eV and tens of eV, ending near the exit in a range of a few eV at the outer edge of the beam to a few tens of eV near the axis. The higher value on-axis is associated with greater numerical noise there, due to the axisymmetric geometry of the calculation, resulting in poorer simulation-particle statistics at small radius. A scatter plot of the macroparticles kinetic energy (KE) versus radius (R) and longitudinal position (0.28 m < z < 3 m) colored by local energy spread is shown in Fig. 2 (b). As expected, there is a correlation of the kinetic energy with radius that is clearly visible at z = 2.8 m and vanishes at the metal plate at z = 3 m. More snapshots from simulations varying the time step, grid resolution and number of macroparticles are given in Appendix II. The macro-particles were collected at the exit plate and their kinetic energy history is plotted in Fig. 3 (left) and contrasted to an experimental measurement using a streak camera shown in Fig. 3. For some types of measurements, averaging over several pulses to improve signal-to-noise will contribute an additional spread that may not be present on any single beam pulse. The upper bound for the energy spread is in the range of a few 100 eV for the experiment while in the range of a few eV for the reported Warp simulations. The Marx voltage exhibits variations in the range of up to several hundreds of volts, playing a significant role in the experimentally measured energy spread, which may account for the difference between the experimental and the simulated bounds.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of a colinear Second-Harmonic Orthogonal Polarization (SHOP) interferometer for electron areal density measurements in Magnetically Insulated Transmission Lines (MITLs)

Experimental measurements of low density plasmas forming in Magnetically Insulated Transmission Line (MITL) regions are desired to improve our understanding of current loss and power flow. Therefore, a new optical interferometer diagnostic was commissioned via this LDRD project. To measure the expected 10 13 - 10 17 cm -3 electron densities inside the 0.5 - 6 mm Anode-Cathode (A-K) gaps, a colinear SHOP interferometer diagnostic was constructed. The diagnostic was initially fielded on the University of New Mexico (UNM) Helicon-Cathode (HelCat) plasma device which provided a highly repeatable and well understood plasma source for which the colinear SHOP interferometer’s functionality could be verified and measured. Utilizing the highly repeatable plasma source and shot averaging, the interferometer was able to achieve an areal density sensitivity of 1×10 14 cm -2 . This work at UNM lead to a Review of Scientific Instruments (RSI) publication [20], DOI:10.1063/5.0101687. After the diagnostic’s capability was proven at UNM, the colinear SHOP interferometer was commissioned for use on the Sandia National Laboratories (SNL) Mykonos accelerator. Here, it provided the first temporal areal density measurements of plasma formation in a parallel plate MITL. The diagnostic was able to achieve a single shot (no multi-shot averaging like at UNM) areal density sensitivity of 1×10 15 cm -2 along a ~ 2mm probing path length, which provided adequate capability to conduct fundamental physics research of MITL plasma formation. CHICAGO and ALEGRA simulations support the diagnostics experimental findings. More experimental and computational work will continue, likely leading to another publication(s). The smaller scale Mykonos accelerator work has also provided justification that the colinear SHOP interferometer is a capable diagnostic for measuring plasma areal densities in the inner MITL and convolute regions of larger TW-class accelerators like SNL’s Z machine.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An improved model for the accurate calculation of parallel heat fluxes at the JET bulk tungsten outer divertor

Parallel heat flux calculations at the JET divertor have been based on the assumption that all incoming heat is due to the projection of the heat flux parallel to the magnetic line, q ∥ , plus a constant background. This simplification led to inconsistencies during the analysis of a series of dedicated tungsten melting experiments performed in 2013, for which infrared (IR) thermography surface measurements could not be recreated through simulations unless the parallel heat flux was reduced by 80% for L-mode and 60% for H-mode. We give an explanation for these differences using a new IR inverse analysis code, a set of geometrical corrections, and most importantly an additional term for the divertor heat flux accounting for non-parallel effects such as cross-field transport, recycled neutrals or charge exchange. This component has been evaluated comparing four different geometries with impinging angles varying from 2 to 90°. Its magnitude corresponds to 1.2%–1.9% of q ∥ , but because it is not affected by the magnetic projection, it accounts for up to 20%–30% of the tile surface heat flux. The geometrical corrections imply a further reduction of 24% of the measured heat flux. In addition, the application of the new inverse code increases the accuracy of the tile heat flux calculation, eliminating any previous discrepancy. The parallel heat flux computed with this new model is actually much lower than previously deduced by inverse analysis of IR temperatures—40% for L-mode and 50% for H-mode—while being independent of the geometry on which it is measured. This main result confirms the validity of the optical projection as long as a non-constant and non-parallel component is considered. For a given total heating power, the model predicts over 10% reduction of the maximum tile surface heat flux compared to strict optical modelling, as well as a 30% reduced sensitivity to manufacturing and assembling tolerances. Furthermore these conclusions, along with the improvement in the predictability of the divertor thermal behaviour, are critical for JET future DT operations, and are also directly applicable to the design of the ITER divertor monoblocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

What next: Further implosion space exploration on the path to NIF extended yield capability

We present quantitative motivations and assessments of various proposed and ongoing directions to further improving yields and target gain of igniting indirect-drive implosions at the National Ignition Facility (NIF). These include increasing compression and confinement time, improving hohlraum and ablator efficiency, and further increasing peak power and laser energy. 1D hydroscaled simulations, augmented by analytic 1D theory, have been used to project yield improvements for each of these implosion optimization tracks, normalized to the best current performing 4 MJ shot. At current NIF capabilities of 2.2 MJ, 450 TW, we project several paths could reach 15 MJ yield levels. We also expect several key implosion physics questions will be addressed in attempting to reach this yield level. These include demonstrating to what extent lower adiabat designs leading to higher compression will increase gain and efficiency, and whether we can reduce residual kinetic energy and ablator-fuel mix that is probably limiting the current burn-up fraction. For an envisaged NIF upgrade to EL = 3 MJ at fixed 450 TW peak power, scaling capsule size and fuel thicknesses faster than pure hydroscaling should allow for yields that could reach up to 60–80 MJ, depending on the efficiency gains realized in increasing deuterium-tritium fuel thickness, reducing hohlraum losses, and switching to lower Z ablators. The laser-plasma instability and beam transmission scaling in these larger hohlraums is shown to be favorable if the spot size is increased with hohlraum scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Examples of Use of SINBAD Database for Nuclear Data and Code Validation

Benchmark experiments can serve several purposes and the availability of suitable experimental data is fundamental to ensure that modelling tools meet the requirements of nuclear industry and licensing needs. Essentially, they are used for: - validation of nuclear cross-section data, - validation of nuclear codes and models used, - qualification and training of users through benchmarking and code comparison exercises. The Shielding Integral Benchmark Archive and Database (SINBAD) project was started in the early 1990's at the OECD Nuclear Energy Agency Data Bank (NEA DB) and the Radiation Safety Information Computational Center (RSICC) as a collection of quality experiments for validating and benchmarking computer codes and data for radiation transport and shielding applications. Currently, the SINBAD database comprises over 100 benchmark compilations and evaluations of relevance to: - reactor shielding, pressure vessel dosimetry (48); - fusion blanket neutronics (31); - accelerator shielding (23). In order to facilitate their use, the data are organized in the standardized format structure. The index and the master document (called abstract) use the Hyper Text Markup Language (HTML) format with links to text files, figures, Portable Document Format (pdf) documents and computer code inputs. The format allows easy access and forward/backward navigation among the data and has proved to be suitable, general enough and above all easy (and low-cost) to maintain. In addition to the characterization of the radiation source, description of the shielding set-up, the instrumentation and the relevant detector measurements, most sets in SINBAD contain also the deterministic or Monte Carlo radiation transport computer model used for the interpretation of the experiment and, where available, results from uncertainty analysis. The set of primary documents used for the benchmark compilation and evaluation are provided in computer readable form. The database is intended for different types of users, including nuclear data evaluators, computer code developers, experiment designers and university students. SINBAD is available from the NEA Data Bank and RSICC (see https://www.oecdnea.org/science/wprs/shielding/sinbad/). Up to now several hundreds of SINBAD packages were distributed from the NEA DB and RSICC. Since 2011 the maintenance and further development of the database is coordinated in the scope of the OECD NEA Nuclear Science Committee (NSC) Working Party on Scientific Issues of Reactor Systems (WPRS) Expert Group on Radiation Transport and Shielding (EGRTS). A key objective of the group is to identify, evaluate and preserve experimental data on shielding benchmarks. The Expert Group maintains close links with the International Reactor Physics Experiment Evaluation (IRPhE) and International Criticality Safety Benchmark Evaluation (ICSBEP) projects, as well as in co-ordination with the JEFF and EFF NEA Data Bank projects and with the SATIF community of accelerator shielding specialists. In the last years, much effort was devoted to the quality review of the available benchmark information.

42 ENGINEERING↗

Overview of Nuclear Integration Activities for Mitigation of Radiation Loads to ITER Superconducting Magnets

The Nuclear Integration Unit (NIU) is a trans-departmental team from the ITER Organization and Domestic Agencies, responsible for all nuclear integration aspects of the facility and for ensuring that systems and components are suitable for their radiation environment and satisfy all related performance requirements. This is a very challenging task namely due to the large number of systems, their still evolving designs, the intricate and far-reaching impacts of these in nuclear responses everywhere, and the complex world-wide distribution of procurement responsibilities. A sound systems engineering approach that involves all relevant experts and stakeholders, and trades-off between different performance requirements in an optimum way, needs to be taken and sustained during design evolution and construction, always with the objective of final overall project success. One of the most important current occupations of the NIU is the evaluation of, and protection from, radiation conditions to superconducting magnets during plasma operations, and particularly nuclear heating. This requires neutronic analyses in massive computer resources, with state-of-the-art methods and complex, extensive computer models. The main ITER superconducting magnet systems subject to significant nuclear loads are the toroidal field coils (TFCs) and the poloidal field coils (PFCs). These are illustrated in Fig. 1. Nuclear responses in these depend mostly on the design of all in-vessel (blankets, manifolds, divertor, in-vessel coils, port plugs), vessel and thermal shield components. Sources to consider are the 500 MW DT burning plasma (14.1 MeV neutrons) and the activated water in the in-vessel cooling system (photons up to 7 MeV from N-16 and neutrons up to 3.5 MeV from N-17). This paper discusses a snapshot of the present understanding of nuclear heat conditions in both the TFC and PFC magnets. This includes integral values and distributions, as well as best estimates from computations and conservative upper estimates by consideration of systematic effects and random uncertainties. The computer models, methods and assumptions used for evaluation of these responses are described. Nuclear integration aspects and ongoing efforts to study and implement additional shielding and other mitigating measures are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Studies of RF sheaths and diagnostics on IShTAR

IShTAR (Ion cyclotron Sheath Test ARrangement) is a linear magnetised plasma test facility for RF sheaths studies at the Max-Planck-Institut für Plasmaphysik in Garching. In contrast to a tokamak, a test stand provides more liberty to impose the parameters and gives better access for the instrumentation and antennas. The project will support the development of diagnostic methods for characterising RF sheaths and validate and improve theoretical predictions. The cylindrical vacuum vessel has a diameter of 1 m and is 1.1 m long. The plasma is created by an external cylindrical plasma source equipped with a helical antenna that has been designed to excite the m=1 helicon mode. In inductive mode, plasma densities and electron temperatures have been characterised with a planar Langmuir probe as a function of gas pressure and input RF power. A 2D array of RF compensated Langmuir probes and a spectrometer are planned. A single strap RF antenna has been designed; the plasma-facing surface is aligned to the cylindrical plasma to ease the modelling. The probes will allow direct measurements of plasma density profiles in front of the RF antenna, and thus a detailed study of the density modifications induced by RF sheaths, which influences the coupling. The RF antenna frequency has been chosen to study different plasma wave interactions: the accessible plasma density range includes an evanescent and propagative behaviour of slow or fast waves, and allows the study of the effect of the lower hybrid resonance layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High energy density Z-pinch plasmas using flow stabilization

The ZaP Flow Z-Pinch research project[1] at the University of Washington investigates the effect of sheared flows on MHD instabilities. Axially flowing Z-pinch plasmas are produced that are 100 cm long with a 1 cm radius. The plasma remains quiescent for many radial Alfvén times and axial flow times. The quiescent periods are characterized by low magnetic mode activity measured at several locations along the plasma column and by stationary visible plasma emission. Plasma evolution is modeled with high-resolution simulation codes – Mach2, WARPX, NIMROD, and HiFi. Plasma flow profiles are experimentally measured with a multi-chord ion Doppler spectrometer. A sheared flow profile is observed to be coincident with the quiescent period, and is consistent with classical plasma viscosity. Equilibrium is determined by diagnostic measurements: interferometry for density; spectroscopy for ion temperature, plasma flow, and density[2]; Thomson scattering for electron temperature; Zeeman splitting for internal magnetic field measurements[3]; and fast framing photography for global structure. Wall stabilization has been investigated computationally and experimentally by removing 70% of the surrounding conducting wall to demonstrate no change in stability behavior.[4] Experimental evidence suggests that the plasma lifetime is only limited by plasma supply and current waveform. The flow Z-pinch concept provides an approach to achieve high energy density plasmas,[5] which are large, easy to diagnose, and persist for extended durations. A new experiment, ZaP-HD, has been built to investigate this approach by separating the flow Z-pinch formation from the radial compression using a triaxial-electrode configuration. This innovation allows more detailed investigations of the sheared flow stabilizing effect, and it allows compression to much higher densities than previously achieved on ZaP by reducing the linear density and increasing the pinch current. Experimental results and scaling analyses will be presented. In addition to studying fundamental plasma science and high energy density physics, the ZaP and ZaP-HD experiments can be applied to laboratory astrophysics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Laser-plasmas in the relativistic-transparency regime: science and applications

Laser-plasma interactions in the novel regime of relativistically induced transparency (RIT) have been harnessed to generate intense ion beams efficiently with average energies exceeding 10 MeV/nucleon (>100 MeV for protons) at “table-top” scales in experiments at the LANL Trident Laser. By further optimization of the laser and target, the RIT regime has been extended into a self-organized plasma mode. This mode yields an ion beam with much narrower energy spread while maintaining high ion energy and conversion efficiency. This mode involves self-generation of persistent high magnetic fields (~10 4 T, according to particle-in-cell simulations of the experiments) at the rear-side of the plasma. These magnetic fields trap the laser-heated multi-MeV electrons, which generate a high localized electrostatic field (~0.1 TV/m). After the laser exits the plasma, this electric field acts on a highly structured ion-beam distribution in phase space to reduce the energy spread, thus separating acceleration and energy-spread reduction. Thus, ion beams with narrow energy peaks at up to 18 MeV/nucleon are generated reproducibly with high efficiency (≈5%). The experimental demonstration has been done with 0.12 PW, high-contrast, 0.6 ps Gaussian 1.053 μm laser pulses irradiating planar foils up to 250 nm thick at 2–8 × 10 20 W/cm 2 . These ion beams with co-propagating electrons have been used on Trident for uniform volumetric isochoric heating to generate and study warm-dense matter at high densities. These beam plasmas have been directed also at a thick Ta disk to generate a directed, intense point-like Bremsstrahlung source of photons peaked at ~2 MeV and used it for point projection radiography of thick high density objects. In addition, prior work on the intense neutron beam driven by an intense deuterium beam generated in the RIT regime has been extended. Neutron spectral control by means of a flexible converter-disk design has been demonstrated, and the neutron beam has been used for point-projection imaging of thick objects. Finally, we discuss the plans and prospects for further improvements and applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transport Properties of Magnetized High-Energy Density Plasma (Final Project Report)

This report summarizes results of project DE-SC0016159 "Transport Properties of Magnetized High-Energy-Density Plasma,'' which ran from 7/15/2016 - 7/14/2021. The primary objectives of the work, as stated in the original proposal, were to develop a theory that describes transport coefficients in magnetized high energy density plasmas and to test the theory with molecular dynamics (MD) simulations. The science challenge is that strong ion coupling, strong magnetization of electrons, and partial degeneracy of electrons, cause the system to be in a regime that is not well described by current theory. The particular processes that were to be investigated include: ion stopping power, diffusion, electron-ion temperature relaxation, thermal conduction and viscosity. All of the primary research objectives were accomplished during the course of this work. In addition, some unexpected results were discovered along the way that led to new and productive research directions. This work resulted in 17 publications in well-respected peer-reviewed journals (primarily Physics of Plasmas and Physical Review E), and 2 more are being prepared for publication. A few of these were selected as Editor's Picks and one was published as a Rapid Communication. Highlights of the research results include: Transport phase space: The first identification of the parameter space that defines fundamental transport regimes in terms of the Coulomb coupling and magnetization parameters. Mean force kinetic theory: Systematic derivation of a kinetic theory for strongly coupled plasmas based on a new expansion parameter of the BBGKY hierarchy. This provided the derivation of a theory that we had previously posed phenomenologically, and also revealed a new term that captures the equation of state properties at all coupling strengths. Transverse friction force: Discovery that the friction force on a test charge in a strongly magnetized plasma includes a component that is perpendicular to the motion of the test charge in the plane defined by the velocity and magnetic field vectors. Gyrofriction force: Discovery that the friction force on a test charge in a plasma that is both strongly magnetized and strongly coupled includes a component of the force in the direction of the Lorentz force. Kinetic theory for warm dense matter: Extension of the mean force kinetic theory to include partial electron degeneracy, so that it applies to warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗