Search NASA⌕ Search

SEARCH · Search NASA

Results for “C Code”

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 991 records · Page 55

A future of inertial confinement fusion without laser-plasma instabilities

From the beginning of inertial confinement fusion (ICF) research, laser-plasma instabilities excited by narrowband lasers have limited the hydrodynamic design space of all laser-based approaches to inertial fusion energy (IFE). With advances in broadband laser technologies, the next generation of ICF drivers will likely have large bandwidth that is engineered to mitigate laser-plasma instabilities, thereby expanding the hydrodynamic design space to include both robust high yields (>200−MJ) with large-energy laser systems (>4 MJ) and high gains (>10) with moderate-energy laser facilities (<2 MJ). State-of-the-art simulations indicate that laser bandwidths of a few percent are required to mitigate instabilities for IFE-relevant conditions. To test these models and demonstrate that high-bandwidth lasers mitigate laser-plasma instabilities, the Fourth-generation Laser for Ultrabroadband eXperiments will be used with the OMEGA Laser System. The goal is to provide the community with the confidence to invest in a multiple-beam high-bandwidth laser facility that will demonstrate the necessary ablation pressures for robust direct-drive ignition without detrimental levels of hot electrons that degrade the implosion performance. It is important to recognize that the highest performing ICF implosions will likely never be completely LPI free because of the significant advantages to maximizing the laser intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatiotemporal evolution of runaway electrons from synchrotron images in Alcator C-Mod

We report in the Alcator C-Mod tokamak, relativistic runaway electron (RE) generation can occur during the flattop current phase of low density, diverted plasma discharges. Due to the high toroidal magnetic field (B 0 = 5.4 T), RE synchrotron radiation is measured by a wide-view camera in the visible wavelength range (λ ≈ 400–900 nm). In this paper, a statistical analysis of over one thousand camera images is performed to investigate the plasma conditions under which synchrotron emission is observed in C-Mod. In addition, the spatiotemporal evolution of REs during one particular discharge is explored in detail via a thorough analysis of the distortion-corrected synchrotron images. To accurately predict RE energies, the kinetic solver CODE is used to evolve the electron momentum-space distribution at six locations throughout the plasma: the magnetic axis and flux surfaces q = 1, 4/3, 3/2, 2, and 3. These results, along with the experimentally-measured magnetic topology and camera geometry, are input into the synthetic diagnostic SOFT to simulate synchrotron emission and detection. Interesting spatial structure near the surface q = 2 is found to coincide with the onset of a locked mode and increased MHD activity. Furthermore, the RE density profile evolution is fit by comparing experimental to synthetic images, providing important insight into RE spatiotemporal dynamics.

synthetic diagnostic↗

JETSCAPE framework: p+p results

The JETSCAPE framework is a modular and versatile Monte Carlo software package for the simulation of high energy nuclear collisions. In this work we present a new tune of JETSCAPE, called PP19, and validate it by comparison to jet-based measurements in p + p collisions, including inclusive single jet cross sections, jet shape observables, fragmentation functions, charged hadron cross sections, and dijet mass cross sections. These observables in p+p collisions provide the baseline for their counterparts in nuclear collisions. Quantifying the level of agreement of JETSCAPE results with p+p data is thus necessary for meaningful applications of JETSCAPE to A+A collisions. The calculations use the JETSCAPE PP19 tune, defined in this paper, based on version 1.0 of the JETSCAPE framework. For the observables discussed in this work calculations using JETSCAPE PP19 agree with data over a wide range of collision energies at a level comparable to standard Monte Carlo codes. These results demonstrate the physics capabilities of the JETSCAPE framework and provide benchmarks for JETSCAPE users.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modifying PyUltraLight to model scalar dark matter with self-interactions

Here we introduce a modification of the pysiultralight code that models the dynamical evolution of ultralight axionlike scalar dark matter fields. Our modified code, pysiultralight, adds a quartic, self-interaction term to reflect the one which arises naturally in axionlike particle models. Using a particle mass of 10−22 eV/c 2 , we show that pysiultralight produces spatially oscillating solitons, exploding solitons, and collapsing solitons which prior analytic work shows will occur with attractive self-interactions. Using our code we calculate the oscillation frequency as a function of soliton mass and equilibrium radius in the presence of attractive self-interactions. We show that when the soliton mass is below the critical mass ($M_c$ = $\frac{\sqrt{3}}{2} M_{max}$) described by Chavanis and the initial radius is within a specific range, solitons are unstable and explode. We test the maximum mass criteria described by Chavanis and Chavanis and Delfini for a soliton to collapse when attractive self-interactions are included. We also analyze both binary soliton collisions and a soliton rotating around a central mass with attractive and repulsive self-interactions. We find that when attractive self-interactions are included, the density profiles get distorted after a binary collision. We also find that a soliton is less susceptible to tidal stripping when attractive self-interactions are included. We find that the opposite is true for repulsive self-interactions in that solitons would be more easily tidally stripped. Including self-interactions might therefore influence the survival timescales of infalling solitons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Phase diagram of the three-dimensional subsystem toric code

Subsystem quantum error-correcting codes typically involve measuring a sequence of noncommuting parity check operators. They can sometimes exhibit greater fault tolerance than conventional codes, which use commuting checks. However, unlike subspace codes, it is unclear if subsystem codes—in particular their advantages—can be understood in terms of ground-state properties of a physical Hamiltonian. In this paper, we address this question for the three-dimensional subsystem toric code (3D STC), as recently constructed by Kubica and Vasmer [], which exhibits single-shot error correction. Motivated by a conjectured relation between single-shot properties and thermal stability, we study the zero- and finite-temperature phases of an associated noncommuting Hamiltonian. By mapping the Hamiltonian model to a pair of 3D Z 2 gauge theories coupled by a kinetic constraint, we find various phases at zero temperature, all separated by first-order transitions: There are 3D toric code-like phases with deconfined point-like excitations in the bulk, and there are phases with a confined bulk supporting a 2D toric code on the surface when appropriate boundary conditions are chosen. The latter is similar to the surface topological order present in 3D STC. However, the similarities between the single-shot correction in 3D STC and the confined phases are only partial: they share the same sets of degrees of freedom, but they are governed by different dynamical rules. Instead, we argue that the process of single-shot error correction can more suitably be associated with a path (rather than a point) in the zero-temperature phase diagram, a perspective, which inspires alternative measurement sequences enabling single-shot error correction. Moreover, since none of the above-mentioned phases survives at nonzero temperature, the single-shot error-correction property of the code does not imply thermal stability of the associated Hamiltonian phase. Published by the American Physical Society 2024

Li, Yaodong (ORCID:0000000337421944)↗

An ab-initio Investigation: The Physical Properties of ScIr{sub 2} Superconductor

Using ab initio technique the physical properties of ScIr{sub 2} superconductor have been investigated with T{sub c} 1.03 K with a MgCu{sub 2}-type structure. We have carried out the plane-wave pseudopotential approach within the framework of the first-principles density functional theory (DFT) implemented within the CASTEP code. The calculated structural parameters confirm a good agreement with the experimental and other theoretical results. Using the Voigt-Reuss-Hill (VRH) averaging scheme the most important elastic properties including the bulk modulus B, shear modulus G, Young’s modulus E and Poisson’s ratio ν of ScIr{sub 2} are determined. At ambient condition, the values of Cauchy pressure and Pugh’s ratio exhibit ductile nature of ScIr{sub 2}. The electronic and optical properties of ScIr{sub 2} were investigated for the first time. The electronic band structure reveals metallic conductivity and the major contribution comes from Ir-5d states. In the ultraviolet region the reflectivity is high up to 50 eV as evident from the reflectivity spectrum.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

KinBot

KinBot 2.0 searches for chemical reaction pathways for gas-phase systems containing C, H, and O atoms. It determines the structure of key intermediates and products, characterizes their molecular structure and related properties, and creates an output that can be conveniently used to calculate rate coefficients by other codes.

Zador, Judit↗

Corrosion of Structural Materials for Advanced Supercritical Carbon- Dioxide Brayton Cycle

The supercritical carbon-dioxide (referred to as SC-CO 2 hereon) Brayton cycle is being considered for power conversion systems for a number of nuclear reactor concepts, including the sodium fast reactor (SFR), fluoride saltcooled high temperature reactor (FHR), and high temperature gas reactor (HTGR), and several types of small modular reactors (SMR). The SC-CO 2 direct cycle gas fast reactor has also been recently proposed. The SC-CO 2 Brayton cycle (discussed in Chapter 1) provides higher efficiencies compared to the Rankine steam cycle due to less compression work stemming from higher SC-CO 2 densities, and allows for smaller components size, fewer components, and simpler cycle layout. For example, in the case of a SFR using a SC-CO 2 Brayton cycle instead of a steam cycle would also eliminate the possibility of sodium-water interactions. The SC-CO 2 cycle has a higher efficiency than the helium Brayton cycle, with the additional advantage of being able to operate at lower temperatures and higher pressures. In general, the SC-CO 2 Brayton cycle is well-suited for any type of nuclear reactor (including SMR) with core outlet temperature above ~ 500°C in either direct or indirect versions. In all the above applications, materials corrosion in high temperature SC-CO 2 is an important consideration, given their expected lifetimes of 20 years or longer. Our discussions with National Laboratories and private industry early on in this project indicated materials corrosion to be one of the significant gaps in the implementation of SC-CO 2 Brayton cycle. Corrosion can lead to a loss of effective load-bearing wall thickness of a component and can potentially lead to the generation of oxide particulate debris which can lead to three-body wear in turbomachinery components. Another environmental degradation effect that is rather unique to CO 2 environment is the possibility for simultaneous occurrence of carburization during oxidation of the material. Carburization can potentially lead to embrittlement of structural alloys in SC-CO 2 Brayton cycle. An important consideration in regards to corrosion is that the temperatures can vary widely across the various sections of the SC-CO 2 Brayton cycle, from room temperature to 750°C, with even higher temperatures being desirable for higher efficiencies. Thus the extent of corrosion and corrosion mechanisms in various components and SC-CO 2 Brayton cycle will be different, requiring a judicious selection of materials for different sections of the cycle. The goal of this project was to address materials corrosion-related challenges, identify appropriate materials, and advance the body of scientific knowledge in the area of high temperature SC-CO 2 corrosion. The focus was on corrosion of materials in SC-CO 2 environment in the temperature range of 450°C to 750°C at a pressure of 2900 psi for exposure duration for up to 1000 hours. The Table below lists the materials tested in the project. The materials were selected based on their high temperature strength, their code certification status, commercial availabilities, and their prior or current usage in the nuclear reactor industry. Additionally, pure Fe, Fe-12%Cr, and Ni-22%Cr were investigated as simple model materials to more clearly understand corrosion mechanisms. This first phase of the project involved testing in research grade SC-CO 2 (99.999% purity). Specially designed autoclaves with high fidelity temperature, pressure, and flow control capabilities were built or modified for this project.

36 MATERIALS SCIENCE↗

MCNP Version 6.2 Release Notes

Monte Carlo N-Particle or MCNP ® is a general-purpose Monte Carlo radiation-transport code designed to track many particle types over broad ranges of energies. This MCNP Version 6.2 follows the MCNP6.1.1 beta version and has been released in order to provide the radiation transport community with the latest feature developments and bug fixes for MCNP. Since the last release of MCNP major work has been conducted to improve the code base, add features, and provide tools to facilitate ease of use of MCNP version 6.2 as well as the analysis of results. These release notes serve as a general guide for the new/improved physics, source, data, tallies, unstructured mesh, code enhancements and tools. For more detailed information on each of the topics, please refer to the appropriate references or the user manual which can be found at http://mcnp.lanl.gov. This release of MCNP version 6.2 contains 39 new features in addition to 172 bug fixes and code enhancements. There are still some 33 known issues the user should familiarize themselves with (see Appendix).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

User's Guide for RESRAD-ONSITE Code

The RESRAD-ONSITE code is equipped with a user-friendly interface that has many features to facilitate its use and comprehension. The design of the interface provides various navigation options, including entering data, performing calculations, and displaying the results of the calculations. General and context-specific help are available, offering information on editing and viewing the radionuclide database, the definitions of input parameters and their use in the calculations, and the selection and output of the results of the calculations into other applications. Two types of sensitivity analysis are supported by the code, deterministic and probabilistic, which can be used to study the influence of input parameters on the calculated results.

54 ENVIRONMENTAL SCIENCES↗

High Brightness Superconducting RF Gun Cavity

The Phase I goals were to numerically optimize a TM020/TM010 SRF gun with a demountable back wall, using SLAC’s state-of-the-art massive parallel ACE3P suite of 3D Finite Element Analysis codes available at the National Energy Research Scientific Computing Center (NERSC) , which includes: a) TM020 cell shape, including impact of RF focusing on beam dynamics; b) optimum number of TM010-cells to achieve low emittances; c) multipacting simulations back wall opening; d) dark current investigations; e) studies on LOM and HOM damping with waveguides/port on TM020-cells using wakefield solver can be conducted as well but are not required stringently; f) conceptual design of a power input coupler; g) multiphysics simulations to optimize the gun thermal management and cooling requirements for a superconducting joint; h) multivariate beam dynamical optimization in conjunction with the gun optimization; and i) finalize an engineering design for vertical high power tests at JLab.

97 MATHEMATICS AND COMPUTING↗

Spatiotemporal evolution of runaway electrons from synchrotron images in Alcator C-Mod

In the Alcator C-Mod tokamak, relativistic runaway electron (RE) generation can occur during the flattop current phase of low density, diverted plasma discharges. Due to the high toroidal magnetic field (B0 = 5.4 T), RE synchrotron radiation is measured by a wide-view camera in the visible wavelength range (~400-900 nm). In this paper, a statistical analysis of over one thousand camera images is performed to investigate the plasma conditions under which synchrotron emission is observed in C-Mod. In addition, the spatiotemporal evolution of REs during one particular discharge is explored in detail via a thorough analysis of the distortion-corrected synchrotron images. To accurately predict RE energies, the kinetic solver CODE [Landreman et al 2014 Comput. Phys. Commun. 185 847-855] is used to evolve the electron momentum-space distribution at six locations throughout the plasma: the magnetic axis and flux surfaces q = 1, 4/3, 3/2, 2, and 3. These results, along with the experimentally-measured magnetic topology and camera geometry, are input into the synthetic diagnostic SOFT [Hoppe et al 2018 Nucl. Fusion 58 026032] to simulate synchrotron emission and detection. Interesting spatial structure near the surface q = 2 is found to coincide with the onset of a locked mode and increased MHD activity. Furthermore, the RE density profile evolution is fit by comparing experimental to synthetic images, providing important insight into RE spatiotemporal dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards Generic Parallel Programming in Computer Science Education with Kokkos

Parallel patterns, views, and spaces are promising abstractions to capture the programmer's intent as well as the contextual information that can be used by an underlying runtime to efficiently map software to parallel hardware. These abstractions can be valuable in cases where an algorithm must accommodate requirements of code and performance portability across hardware architectures and vendor programming models. Kokkos is a parallel programming model for host- and accelerator architectures that relies on these abstractions and targets these requirements. It consists of a pure C++ interface, a specification, and a programming library. The programming library exposes patterns and types and maps them to an underlying abstract machine model. The abstract machine model offers a generic view of parallel hardware. While Kokkos is gaining popularity in large-scale HPC applications at some DOE laboratories, we believe that the implemented concepts are of interest to a broader audience including academia as they may contribute to a generic, vendor, and architecture-independent education of parallel programming. In this work, we give an insight into the design considerations of this programming model and list important abstractions. Further, we document best practices obtained from giving virtual classes on Kokkos and give pointers to resources that the reader may consider valuable for a lecture on generic parallel programming for students with preexisting knowledge on this matter.

Ciesko, Jan↗

ECP Milestone Memo WBS 2.3.4.13 ECP/VTK-m FY18Q4 [MS-18/09-10] Dynamic Types / Rendering Topologies STDA05-16

The STDA05-16 milestone comprises the following 3 distinct deliverables. OpenMP VTK-m currently supports three types of devices: serial CPU, TBB, and CUDA. To run algorithms on multicore CPU-type devices (such as Xeon and Xeon Phi), TBB is required. However, there are known issues with integrating a software product using TBB with another one using OpenMP. Therefore, we will add an OpenMP device to the VTK-m software. When engaged, this device will run parallel algorithms using OpenMP directives. This will mesh more nicely with other code also using OpenMP. Rendering Topological Entities VTK-m currently supports surface rendering by tessellation of data structures,and rendering the resulting triangles. We will extend current functionality to include face, edge, and point rendering. Better Dynamic Types Impl For the best efficiency across all platforms, VTK-m algorithms use static typing with C++ templates. However, many libraries like VTK, ParaView, and Visit use dynamic types with virtual functions because data types often cannot be determined at compile time. We have an interface in VTK-m to merge these two typing mechanisms by generating all possible combinations of static types when faced with a dynamic type. Although this mechanism works, it generates very large executables and takes a very long time to compile. As we move forward, it is clear that these problems will get worse and become infeasible at exascale. We will rectify the problem by introducing some level of virtual methods, which require only a single code path, within VTK-m algorithms. This first milestone produces a design document to propose an approach to the new system.

97 MATHEMATICS AND COMPUTING↗

Evaluation of the Effect of Burnup Acceleration on UO 2 Microstructure Evolution

Accelerated fuel qualification (AFQ) is a methodology by which new nuclear fuels are developed in an accelerated time frame compared with historical fuel qualification approaches. AFQ generally relies on high-fidelity physics-based modeling and simulation tools to adequately describe fuel performance as well as on revolutionary methods to accelerate burnup accumulation and collect relevant data more quickly. This report summarizes the use of advanced fuel modeling and simulation tools to evaluate microstructures from commercially irradiated fuel and microstructures from proposed MiniFuel irradiations, in which burnup accumulation is accelerated while prototypic temperature conditions are maintained. In this milestone, we used the mesoscale fuel performance code MARMOT to model the evolution of irradiated UO 2 microstructures and their potential restructuring at high burnup. The simulation conditions were informed by BISON models of both commercially irradiated fuel and MiniFuel. A first set of simulations investigated the recrystallization behavior of fully dense microstructures and showed full recrystallization at burnups as low as 52 MWd/kgU at 950°C. However, these simulations did not account for the presence of fission gas bubbles (predicted by BISON). Therefore, a second set of simulations including fission gas bubbles was performed and indicated that at the lowest temperature considered (650°C), the porous UO 2 microstructures have the highest total Gibbs free energies and are likely to recrystallize earlier than higher temperature cases (800 and 950°C), which agrees with high-burnup fuel characterization data. The results also showed that at lower temperature (650°C), the total free energies of the PWR fuel and MiniFuel microstructures are not significantly different. However, at the highest temperature (950°C), MiniFuel microstructures have a lower free energy than that of the PWR fuel microstructure. The competing effects between the temperature-dependent grain nucleation rate and the reduction of the free energy of the microstructure at higher temperature as a result of diffusion indicated that restructuring may occur at even higher temperatures than those considered in this study. In addition to the microstructure evolution modeling efforts, the burnup gradient across a single fuel specimen was also considered. This evaluation was for the VXF-15 position of the High Flux Isotope Reactor (HFIR) using the code suite HFIRCON, which was developed to automate the workflow for evaluating targets and fuel as they are irradiated in HFIR. The burnup gradient evaluation showed a dependence on both the axial and radial locations within the specimen, with a maximum difference of 1.7 between the inner and outermost radial layers. This relationship was further supported by considering the fission product speciation with respect to location within the specimen, which showed a higher concentration of 239 Pu, 240 Pu, and 241 Pu on the outer radial locations of the specimen than the center. The findings of the burnup and speciation evaluation show that some amount of self-shielding is occurring in the specimen when irradiated in the high-flux environment of HFIR; however, this impact is more pronounced for natural uranium when compared to 6% enrichment due to the higher ratio of 238 U in the specimen. Further analyses are required to understand the sensitivity of this gradient to spatial mesh and enrichment of the specimen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Numerical study of transition between even and odd toroidal Alfvén eigenmodes on EAST

In this work, linear simulations of toroidal Alfvén eigenmodes (TAEs) driven by energetic particles (EPs) on EAST (Experimental Advanced Superconducting Tokamak) are performed using the hybrid-kinetic MHD model implemented in NIMROD code. Furthermore, the EAST equilibrium is reconstructed using the EFIT code based on the experimental measurement. The “slowing down” distribution is adopted for modeling the equilibrium distribution of the energetic ions from the deuterium neutral beam injection on EAST. The frequency, the dominant poloidal mode number, the radial location and the detailed 2D mode structure of the TAE/RSAE (Reverse Shear Alfvén Eigenmode)/EPM (energetic particle mode) modes are consistent between the eigenvalue analysis and the NIMROD simulation. As the β fraction of EP increases, a transition between even and odd TAEs occurs, along with that between the ballooning and antiballooning mode structures. It is demonstrated that when the β fraction of EP is close to the transition threshold, both types of TAEs coexist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlocal effects on thermal transport in hydrodynamic simulations of unmagnetized MagLIF-relevant gaspipes on NIF

We present simulations of heat flow relevant to gaspipe experiments on the National Ignition Facility to investigate kinetic effects on transport phenomena. D 2 and neopentane (C 5 H 12 ) filled targets are used to study the laser preheat stage of a MagLIF scheme where an axial magnetic field is sometimes applied to the target. Simulations were done with the radiation-MHD code HYDRA with a collision-dominated fluid model and the SNB nonlocal electron thermal conduction model. Using the SNB model to evolve the electron temperature increased the heat front propagation of neopentane gas targets compared to a local model by limiting radial heat flow. This increases electron temperature near the axis, which decreases laser absorption. We find that the effect of heat flow models on temperature profiles and laser propagation is modest. Beyond the SNB model, we utilize HYDRA to initialize plasma conditions for the Vlasov–Fokker–Planck K2 code. We run K2 until a quasi-steady state is reached and examine the impact of kinetic effects on heat transport. Although axial heat flow is well predicted by fluid models, the fluid model consistently overpredicts radial heat flow up to 150% in regions with the largest temperature gradient of D 2 filled gaspipes. On the other hand, the SNB nonlocal electron conduction model is found to be adequate for capturing kinetic heat flow in gaspipes.

Lau, Ryan Y. [Univ. of Colorado, Boulder, CO (Unit↗

Multiple Analytical Approach to Isotopic Transport Analysis in Magnetic Fusion Devices

High-Z impurities released from plasma-material interactions have been shown to limit the performance of fusion plasmas, and understanding these impurity transport mechanisms throughout the plasma scrape-off layer is a major challenge. Presented herein is a study of tungsten (W) erosion and transport by uniquely measuring absolute quantities of isotopic W in order to determine the source of natural and enriched 182W isotopes that have traveled throughout the tokamak discharges on the DIII-D National Fusion Facility at General Atomics. Two primary analysis methods have been implemented to characterize this W on graphite collector probes that were inserted into DIII-D’s outboard midplane. Results from experiments using Rutherford backscattering spectrometry (RBS) have measured W particle areal densities down the centerline of the probes as high as 6E14 atoms/cm2 with a detection limit of 1E12 atoms/cm2. Laser ablation inductively coupled plasma mass spectrometry (LAMS) has confirmed the elemental trends found with RBS and has provided additional insight into collector probe surface profiles. Two-dimensional elemental and isotopic maps from LAMS are used to reveal new collector probe features and further refine the source of collected W. Variations in isotopic profiles and total W content are coupled to (a) the face of the probe being analyzed, (b) the dimensions of the probe, and (c) the plasma pulse parameters that were used during probe exposure. These results provide one-of-a-kind empirical evidence that is now being utilized for validation of tokamak impurity transport through theoretical models and in codes such as 3D-LIM and OEDGE.

Duran, Jonah↗