Comparison of a collisional-radiative fluid model of H2 in UEDGE to the kinetic neutral code EIRENE
Not Available
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
SEARCH · Engineering Papers
Results for “D CODES”
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.
Benchmark between antenna code TOPICA, RAPLICASOL and Petra-M for the ICRH ITER antenna
ITER will be equipped with three plasma heating systems: neutral beam (NB), electron cyclotron (EC), and ion cy-clotron resonance heating (ICRH). The latter consists of two identical ICRH antennas to deliver 20 MW to the plasma (baseline, upgradable to 40 MW). ICRH will play a crucial role in the ignition and sustainment of burning plasmas in ITER. A high fidelity and robust modeling effort to understand the interaction of the IC waves with the scrape-off-layer (SOL) plasma is a very important aspect. Among the main important research topics, we have the assessment of the antenna loading for different plasma scenarios, the role of the lower hybrid resonance in front of the antenna and how to include it in our models, and the RF sheath boundary conditions to evaluate the antenna impurity generation. In this work, we tackle the first of these by reporting on ICRF simulations employing the Petra-M code, which is an electromagnetic simulation tool for modeling RF wave propagation based on MFEM [http://mfem.org] for the ITER ICRH antenna. Moreover, a benchmark between the well tested antenna codes TOPICA, RAPLI-CASOL, which is based on COMSOL [www.comsol.com], and the Petra-M code is also presented. In conclusion, S- and Z-matrices and wave electric field are compared showing an excellent agreement among these codes.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Neoclassical transport of impurities in tokamaks with non-axisymmetric perturbations
The effect of resonant magnetic perturbations (RMPs) on neoclassical transport of impurities is calculated in DIII-D and NSTX. Here, the neoclassical fluxes are evaluated using the NEO code with nonlinear one-fluid nonaxisymmetric equilibrium calculated using M3D-C1. Neoclassical fluxes of impurities show significant changes with RMPs if the impurity resides in low-collisionality regime, but are weakly affected by RMPs in the Pfirsch–Schlüter (P–S) regime. Charge number (Z) of the impurity affects the collisionality of impurity species, which then determines the effect of 3D fields on neoclassical impurity transport. This suggests that RMPs can be possibly used for exhaust of low-Z impurities in these facilities, but have little effect for impurities with charge number greater than 10 or so. Additionally, it is shown that it is the change in the convective fluxes which is a main contributor in modifying the neoclassical impurity fluxes in the presence of RMPs.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
LibERI—A portable and performant multi-GPU accelerated library for electron repulsion integrals via OpenMP offloading and standard language parallelism
A portable and performant graphics processing unit (GPU)-accelerated library for electron repulsion integral (ERI) evaluation, named LibERI, has been developed and implemented via directive-based (e.g., OpenMP and OpenACC) and standard language parallelism (e.g., Fortran DO CONCURRENT). Offloaded ERIs consist of integrals over low and high contraction s, p, and d functions using the rotated-axis and Rys quadrature methods. GPU codes are factorized based on previous developments with two layers of integral screening and quartet presorting. In this work, the density screening is moved to the GPU to enhance the computational efficacy for large molecular systems. Here, the L-shells in the Pople basis set are also separated into pure S and P shells to increase the ERI homogeneity and reduce atomic operations and the memory footprint. LibERI is compatible with any quantum chemistry drivers supporting the MolSSI Driver Interface. Benchmark calculations of LibERI interfaced with the GAMESS software package were carried out on various GPU architectures and molecular systems. The results show that the LibERI performance is comparable to other state-of-the-art GPU-accelerated codes (e.g., TeraChem and GMSHPC) and, in some cases, outperforms conventionally developed ERI CUDA kernels (e.g., QUICK) while fully maintaining portability.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Interpretive MHD modeling of dispersive shell pellet injection for rapid shutdown in tokamaks
Dispersive shell pellet (DSP) injection is modeled with the extended-MHD code NIMROD for interpretive insight into the results of recent DIII-D DSP experiments and to explore the dynamics of an inside-out thermal quench for disruption mitigation in tokamaks. Simulations of the pre-thermal quench (TQ) phase indicate that the upper bound for the quantity of ablated carbon shell material that will not perturb the flux surfaces is in the ballpark of, but somewhat below the experimental quantity. Even below this quantity, sufficient electrons are added to the plasma by the shell material to produce significant dilution cooling before the TQ is triggered. Simulations carried through the end of the TQ have very large amplitude MHD fluctuations (δB/B>10 -2 ) at the time of the plasma current spike associated with current profile redistribution. Finally, after the plasma current spike, which is of comparable amplitude to that measured in DIII-D experiments, none of the runaway electron test-particles whose orbits are tracked throughout the simulation remain confined.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Modeling the effect of MHD activity on runaway electron generation during SPARC disruptions
Magnetohydrodynamic (MHD) instabilities and runaway electrons (REs) interact in several ways, making it important to self-consistently model these interactions for accurate predictions of RE generation and the design of mitigation strategies, such as massive gas injection (MGI). Using M3D-C1 – an extended MHD code with a RE fluid model – we investigate the effects of 3-D nonlinear MHD activity, material injection, and 2-D axisymmetric vertical displacement events (VDEs) on RE evolution during disruptions on SPARC – a high-field, high-current tokamak designed to achieve a fusion gain Q > 1. Several cases, comprising different combinations of neon (Ne) and deuterium (D 2 ) injection, are considered. Here, our results demonstrate key effects that arise from the self-consistent RE + MHD coupling, such as an initial increase in RE generation due to MHD instability growth, decreased saturation energies of the m/n = 1/1 mode driving sawteeth-like activity, RE losses in stochastic magnetic fields, and subsequent RE confinement and plateau formation due to re-healing of flux surfaces. Large RE plateaus (>5 MA) are obtained with Ne-only injection (2-5 × 10 21 atoms), while combined D 2 + Ne injection (2 × 10 21 Ne atoms; 1.8 × 10 22 D 2 molecules) produces a lower RE current (<2 MA). With D 2 + Ne injection, a post thermal quench "cold" VDE terminates the RE beam, preventing a steady plateau. These simulations couple REs, 3-D MHD instabilities, MGI, and axisymmetric VDEs for the first time in SPARC disruption simulations and represent a crucial step in understanding RE generation and mitigation in high-current devices like SPARC.
Datta, Rishabh [Massachusetts Inst. of Technology ↗
TECHNICAL BACKGROUND FOR THE UPDATE AND EXTENSION OF THE TIME-DEPENDENT ALLOWABLE STRESSES OF TYPE 304H AND 316H STAINLESS STEELS
Type 304H and 316H stainless steels are codified in the ASME Boiler and Pressure Vessel Code for Section III, Division 5, Class A construction for up to 1500F (816C) and 300,000 hours. Extension of the design lifetimes to 500,000 hours has been undertaken by various research projects and by the ASME Section III, Division 5 Code Committees. Due to a long-standing issue related to non-classical creep behavior on the use of the time to the onset of tertiary creep as one of the time-dependent allowable stress criteria, little progress was made in the extension of the time-dependent allowable stresses for Type 304H and 316H stainless steels. A recent effort by Dabrow and Nestell (Impact of Tertiary Creep on Time Dependent Allowable Stresses for Type 304H and 316H Stainless Steels, Report 0300-0003-RPT-001, MPR Associates, Alexandria, VA, 2020) has provided a rational method for the treatment of the highly variable data for the time to onset of tertiary creep. In this paper, the background formulas for determining the time-dependent allowable stress values of Type 304H and 316H stainless steels, as assembled from relevant project work, are presented. These updates and extensions are being recommended to the ASME Code committees for approval, and are subject to change upon feedback from ASME Section III Code committees.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Evaluating Automated Face Identity-Masking Methods with Human Perception and a Deep Convolutional Neural Network
Face de-identification (or “masking”) algorithms have been developed in response to the prevalent use of video recordings in public places. Here, we evaluated the success of face identity masking for human perceivers and a deep convolutional neural network (DCNN). Eight de-identification algorithms were applied to videos of drivers’ faces, while they actively operated a motor vehicle. These masks were pre-selected to be applicable to low-quality video and to maintain coarse information about facial actions. Humans studied high-resolution images to learn driver identities and were tested on their recognition of active drivers in low-resolution videos. Faces in the videos were either unmasked or were masked by one of the eight algorithms. When participants were tested immediately after learning (Experiment 1), all masks reduced identification, with six of eight masks reducing identification to extremely poor performance. In a second experiment, two of the most effective masks were tested after a delay of 7 or 28 days. The delay did not further reduce identification of the masked faces. In all masked conditions, participants maintained stringent decision criteria, with low confidence in recognition, further indicating the effectiveness of the masks. Next, the DCNN performed an identity-matching task between high-resolution images and masked videos—a task analogous to that done by humans. The pattern of accuracy for the DCNN mirrored some, but not all, aspects of human performance, highlighting the need to test the effectiveness of identity masking for both humans and machines. The DCNN was also tested on its ability to match identity between masked and unmasked versions of the same video, based only on the face. DCNN performance for the eight masks offers insight into the nature of the information in faces that is coded in these networks.
97 MATHEMATICS AND COMPUTING↗
Lattice Optimization for Porous Electrode
This code performs an optimization over the structure of a porous electrode formed from a lattice of unit cells. This is a steepest descent version of the code used in: V. A. Beck, J. J. Wong, C. F. Jekel, D. A. Tortorelli, S. E. Baker, E. B. Duoss and M. A. Worsley, Computational design of microarchitected porous electrodes for redox flow batteries, J. Power Sources, 2021, 512, 230453.
Wong, JonathanJ↗
Computer Science Research Needs for Parallel Discrete Event Simulation (PDES)
Historically, scientific computing efforts have demonstrated the clear need for, and effective use of, supercomputing with traditional time-stepped simulations. Nevertheless, there are several areas in the mission spaces of the U.S. Department of Energy and other agencies waiting to tap advanced computing research using a different, discrete event style of modeling, simulation, and analysis. These span a wide spectrum of applications including energy grid resilience, urban planning and policy, transportation science, building technologies, emergency response and planning, environmental impact analysis, computational epidemiology, Internet communications, cyber security, and cyber-physical systems, to name only a few. Even within traditional scientific applications, the role of discrete event modes of execution is increasing in the form of new event-based mathematical solvers such as quantized state integration methods and discrete-continuous hybrid system solvers. Co-design of advanced supercomputing hardware systems is another area that exploits discrete event simulation at its core for effective analyses. Complex systems, entity behaviors and interconnections play a significant role in all these applications, which are mapped to large-scale models with discrete event formulations. To make advancements in all the aforementioned scientific areas, many technical aspects need to be more thoroughly studied and deeply understood in parallel discrete event simulation (PDES). The unique dynamics inherent in a discrete event modeling approach, by their very nature, intersect and influence the entire stack of the computing system, including (a) the unique nature of the instruction sets exercised in PDES workloads without a predominance of high-precision floating point operations, (b) virtual time-constrained multi-threaded execution of many logical processes per processor, (c) extremely variable and difficult to predict network traffic characteristics, (d) interfaces and inter-dependencies with machine learning and artificial intelligence codes at higher software layers, and (e) highly challenging load balancing needs, especially in effectively accounting for accelerated/extremely heterogeneous computing in current and future high-performance computing systems. Efficient and accurate parallel execution of PDES workloads is also dominated by challenges in dealing with their asynchronous concurrency fundamentally present at the model level. Conservative synchronization, optimistic/speculative synchronization, and their hybrid schemes open new questions in fundamental computer science with respect to reversibility of computation and prediction (lookahead) of behaviors inherent within model codes. On the implementation front, there are relatively few scalable, general-purpose parallel discrete event simulators in the world, and even fewer have been studied on emerging hardware platforms. To enable scientific advances using PDES, the research needs in computer science must also be pursued and met in the intersection of the algorithmic and hardware-aware aspects of scalable PDES engines. This report is aimed at capturing a computer science-oriented view of this important area of research in PDES, presenting a sample of important applications with their inherent discrete event technology elements. Needs are outlined in core areas of parallel discrete event research as well as cross-cutting directions in computer science research that positively impact scientific advancements across several important application areas. A selection of priority research opportunities in advanced computing for PDES is identified to serve as reference for key research topics and their order of importance for scientific advancements.
97 MATHEMATICS AND COMPUTING↗
Thermal Hydraulic Experimental Test Article - Fiscal Year 2022 Final Report
The Thermal Hydraulic Experimental Test Article (THETA) is a METL vessel experiment designed for testing and validating sodium fast reactor components and phenomena. THETA has been scaled using a non-dimensional Richardson number approach to represent temperature distributions during nominal and loss of flow conditions in a sodium fast reactor (SFR), this analysis was detailed in the THETA FY19 report. The facility is being constructed with versatility in mind, allowing for the installation of various immersion heaters, heat pipes, and heat exchangers without significant facility modification. THETA was designed in collaboration with systems code experts to inform the geometry and sensor placement to acquire the highest value code validation data.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Investigating Fission Reaction Rate Ratio Sensitivities [Abstract]
Reaction rate ratios are a measurable parameter for reactor and criticality applications. A number of foil irradiations and fission chamber measurements have been performed for critical assemblies at Los Alamos National Laboratory starting in the 1950’s including (i) Godiva, a bare HEU spherical assembly; (ii) Flattop-25, a spherical assembly consisting of an HEU core and a natural uranium reflector; (iii) Jezebel, a bare 239 Pu assembly; and (iv) Flattop-Pu, a spherical assembly consisting of a 239Pu core and a natural uranium reflector. Fission ratio data for 238 U(n,f)/ 235 U(n,f), 237 Np(n,f)/ 235 U(n,f), 233 U(n,f)/ 235 U(n,f) and 239 Pu(n,f)/ 235 U(n,f) were obtained and reported. The EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) project at Los Alamos National Laboratory (LANL) aims to constrain nuclear data by using a suite of measurement types beyond k-effective. Recent investigations include the use of pulsed spheres for nuclear data validation and other measurement methods of interest. One focus of the work is to determine if other methods are complimentary to the critical experiments utilized for nuclear data validation. It is anticipated the investigations will help inform methods that may be utilized in machine learning algorithms for nuclear validation. In order to use a measurement type for nuclear validation, it is necessary to obtain cross-section sensitivities for parameters. This work looks at reaction rate ratio sensitivities with SENSMG and Monte Carlo N-Particle R Code Version 6.21.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Integration of MELCOR for Fusion and TMAP4 for fusion reactor systems analysis and tritium inventory tracking
This is the poster presentation to be presented at a virtual poster presentation session at the 32nd SOFT conference. Discussed are the background and timeline of the development of the MELCOR-TMAP code developed at Idaho National Laboratory.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Out-of-Pile Testing and Instrumentation Transient Water Irradiation System
Current initiatives to increase the burnup of conventional nuclear fuels past the approximate 62 GWd/t limit have been spurred on by direct savings to refueling and waste storage. The technical justification for a new license limit requires extensive qualification through experimental testing. Unlike beginning-of-life fuels, high-burnup fuels are more susceptible to fuel fragmentation, relocation, and dispersal (FFRD), therefore more data is needed to characterize fuels under key accident scenarios. The Transient Reactor Test Facility (TREAT) located at the Idaho National Laboratory has developed a testing apparatus architecture to test fuels and claddings at prototypic conditions. The Transient Water Irradiation System (TWIST) is the latest iteration of a testing device capable of conducting loss of coolant accidents (LOCAs) in TREAT. The Out-of-Pile Testing and Instrumentation TWIST (OPTI-TWIST) is an electrically heated device that is analogous to TWIST. OPTI-TWIST allows for detailed instrumentation and thermal-hydraulic characterization. TWIST ultimately aims to conduct the most advanced in-situ diagnostics to evaluate FFRD in a prototypic LOCA. Moreover, it will explore the phenomenological bifurcation of a decay-energy heat up driven LOCA and a stored-energy heat up driven LOCA. The instrumentation suite includes conventional thermocouples and pressure transducers in addition to an electro impedance sensor, an acoustic emission sensor, an optical pressure sensor, and an optical pyrometer. Characterizing these instruments in OPTI-TWIST eliminates complications of irradiation effects while preserving extreme thermal-hydraulic conditions. Finally, benchmarking both devices to a thermal-hydraulic code like the Reactor Excursion and Leak Analysis Program (RELAP)5-3D provides a unique opportunity for iteration. Pre-test predictions and post-test interpretations inform the physical designs, operational procedures, test conditions, and instrumentation types and positions.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
The Impact of the Halden Reactor Project on the Development of the Bison Fuel Performance Code
The Bison fuel performance code has been under development at Idaho National Laboratory by various programs within the U.S. Department of Energy since 2009. As part of any fuel performance code development, validation is necessary. To this end, several Halden experiments have been used. Early on, the focus was on validating the thermal mechanical behaviour, including pellet cladding mechanical interaction and fission gas release models implemented in Bison. Support for design-basis accident capabilities began in approximately 2014. Also in 2014, the development of models for accident tolerant fuel (ATF) concepts also began. This includes dopedUO2 fuel for which Halden has some experiments with temperature and fission gas release measurements. The Bison team has been fortunate to have two senior developers be secondees to Halden for a year: the late Dr. Giovanni Pastore from May 2015 to April 2016 and the now retired Dr. Richard Williamson in 2018. The paper will provide examples of Bison validation using Halden data for standard UO2 fuel under normal operation and loss-of-coolant accident condition, and ATF concepts. A tribute to the legacy that both the late Giovanni and retired Richard left on the Bison code is provided.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Study of Codes & Standards for Energy Storage Systems: A Report to Congress
The Infrastructure Investment and Jobs Act (H.R. 3684, 2021) directed the Secretary of Energy to prepare a report identifying the existing codes and standards for energy storage technologies. The stated goals for the report are to enhance the safe development of energy storage systems by identifying codes that require updating and facilitation of greater conformity in codes across different types and usages of energy storage technologies. This paper will focus on the specific codes and standards for stationary energy storage systems (ESS).
25 ENERGY STORAGE↗
AMPX and SCALE Nuclear Data Libraries for Depletion [Slides]
This presentation touches on AMPX and SCALE nuclear data libraries for depletion. AMPX is a nuclear data processing code suite that generates continuous-energy and multigroup neutron and gamma libraries used in SCALE and VERA, covariance data used in S/U calculations, and depletion, activation and decay origin data for Oak Ridge Isotope Generation code (ORIGEN).
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Linear simulation of magnetohydrodynamic plasma response to three-dimensional magnetic perturbations in high-β P plasmas
In this work, we report the numerical analyses of linear magnetohydrodynamics (MHD) plasma response to applied three-dimensional magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β P (poloidal beta) plasmas, utilizing the extended-MHD code M3D-C1, with the purpose of realizing a better understanding of the existing experiment in which the n=3 MPs were applied to such high-β P plasmas attempting to control large amplitude type-I ELMs. Such high-β P plasmas obtained at the DIII-D tokamak feature an upper-biased double null configuration, a high edge safety factor q 95 ~6.4, and a stable internal transport barrier (ITB) leading to relatively high core pressures. Single-fluid simulations show that the plasma response to n=3 MPs, including both non-resonant/kinking and resonant components, is significantly weaker than that to n=1 or 2 MPs. To survey the impact of q 95 on plasma response to applied MPs, the SEGWAY (Self-consistent Equilibrium Generating Workflow for AnalYsis) module, developed in the OMFIT integrated modelling framework, is employed to generate a series of equilibria with a wide range of q 95 while other key parameters including the normalized beta, electron density at pedestal top, and plasma shape are kept fixed. Compared to the vacuum response, single-fluid M3D-C1 simulations predict a much more significant decrease of resonant plasma response to the applied n=3 MPs at the maximum penetration radii as q 95 increases. In contrast to single-fluid simulation results showing resonant penetration occurs only near the pedestal top where the E×B toroidal rotation frequency is zero, two-fluid simulations show two comparable resonant penetrations locating near the pedestal top and the ITB foot, where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot may be responsible for the observed formation of a staircase structure in both electron density and temperature profiles and thereby a considerable deterioration of global plasma performance when MPs are applied in high-β P plasmas. Motivated by this numerical work, we provide some ideas for the future research, with the purpose of realizing effective ELM control in such high-β P plasmas on the DIII-D and EAST devices.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗