Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “C 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.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cyclic Ramberg–Osgood Parameters of 316SS Weld Metal Under In-Air and PWR-Water Environments

We report at present, the available Ramberg–Osgood (R–O) parameters for different metals (e.g., in ASME code and other literatures) are static (generally based on a tensile curve). These static R–O parameters cannot accurately model the cyclic plasticity behavior. This work presents the cyclic R–O material hardening parameters for 316 stainless steel similar metal welds. The parameters were estimated under various conditions (in-air at room temperature, 300 °C in-air, and in-air at primary water conditions for a pressurized water reactor (PWR)). It is anticipated that the reported results would be useful for computational mechanics-based shakedown analysis and fatigue life estimation of PWR components.

42 ENGINEERING↗
DOE OSTI2024

Common Electric Power Transmission System Model JSON Schema Specification

The Common Electric Power Transmission System Model (CTM) is an intuitive, extensible, language-agnostic, and error-resistant specification of electric power network components parameter names and units, and relation between components, intended for use by the research community developing new computational methods for power systems operations and simulation. Power system datasets following the CTM specification can be read as dictionaries and manipulated in that form in most programming languages (e.g., Python, Julia, C++). This standard data structure in CTM makes it easy to work in multiple power systems domains (e.g., economic operation, reliability assessment, electricity markets, stability assessment, etc.) without requiring conversions between use-case-specific file formats with information loss in the process. This repository specifies CTM as a JSON Schema, provides documentation, derivate (code-generated) implementations of CTM, and example data and usage of the schema for important use cases.

Aravena Solis, Ignacio↗
DOE OSTI2022

Time-explicit Darwin PIC algorithm

A new approach to Darwin particle-in-cell plasma simulation is described. Using a finite-element approach, the vector potential is assured to be exactly solenoidal. This allows writing the system of multi-species particles as an action-at-a-distance Hamiltonian system. Applying recently-developed explicit symplectic methods to this gives a time-explicit algorithm with all desired properties. Elliptic systems for the electrostatic and magneto-static portions of the problem are inverted efficiently by an algebraic multi-grid algorithm. The algorithm is implemented in a two-dimensional Cartesian-geometry code and tested by application to Weibel instability and to Alfvén wave dynamics. Results show the effectiveness of this approach. Extensions to arbitrary meshes and to a partially time-implicit scheme are briefly discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The MCNP ® 6 code: A decade of progress

After several years of effort involved in merging the Los Alamos National Laboratory MCNP5 and MCNPX codes, in 2013 the first production release of version 6 of the Monte Carlo N-Particle ® , or MCNP ® , code MCNP6.1 was distributed publicly. Since then, three significant releases have been issued: MCNP6.1.1beta in 2014, MCNP6.2 in 2018, and MCNP6.3 in 2023. While each release always contains new features, code enhancements, and bug fixes, each version has had a different primary focus, ranging from improved calculational efficiency to new powerful utilities and tools, to software modernization of the code base. With all that has been learned over the first decade of the MCNP6 code, continuous progress is being made toward a modernized, general-purpose Monte Carlo radiation transport code that remains a trusted resource for the global community of practitioners. This paper describes these first 10+ years of the MCNP6 code and its continually improving data libraries, and gives some insight into how the next decade is expected to unfold.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

First coupled GENE–XGC microturbulence simulations

Covering the core and the edge region of a tokamak, respectively, the two gyrokinetic turbulence codes Gyrokinetic Electromagnetic Numerical Experiment (GENE) and X-point Gyrokinetic Code (XGC) have been successfully coupled by exchanging three-dimensional charge density data needed to solve the gyrokinetic Poisson equation over the entire spatial domain. Certain challenges for the coupling procedure arise from the fact that the two codes employ completely different numerical methods. This includes, in particular, the necessity to introduce mapping procedures for the transfer of data between the unstructured triangular mesh of XGC and the logically rectangular grid (in a combination of real and Fourier space) used by GENE. Constraints on the coupling scheme are also imposed by the use of different time integrators. First, coupled simulations are presented. We have considered collisionless ion temperature gradient turbulence, in both circular and fully shaped plasmas. Coupled simulations successfully reproduce both GENE and XGC reference results, confirming the validity of the code coupling approach toward a whole device model. Here, many lessons learned in the present context, in particular, the need for a coupling procedure as flexible as possible, should be valuable to our and other efforts to couple different kinds of codes in pursuit of a more comprehensive description of complex real-world systems and will drive our further developments of a whole device model for fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Bulk viscosity and cavitation in heavy ion collisions

Relativistic heavy ion collisions generate nuclear-sized droplets of quark-gluon plasma (QGP) that exhibit nearly inviscid hydrodynamic expansion. Smaller collision systems such as p + Au, d + Au, and 3 He + Au at the Relativistic Heavy Ion Collider, as well as p + Pb and high-multiplicity p + p at the Large Hadron Collider may create even smaller droplets of QGP. If so, the standard time evolution paradigm of heavy ion collisions may be extended to these smaller systems. These small systems present a unique opportunity to examine pre-hydrodynamic physics and extract properties of the QGP, such as the bulk viscosity, where the short lifetimes of the small droplets makes them more sensitive to these contributions. Here we focus on the influence of bulk viscosity, its temperature dependence, and cavitation effects on the dynamics in small and large systems using the publicly available hydrodynamic codes SONIC and MUSIC. In conclusion, we also compare pre-hydrodynamic physics in different frameworks including AdS/CFT strong coupling, IP-GLASMA weak coupling, and free streaming or no coupling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

MELCOR Code Change History (Revision 14959 to 18019)

This document summarily provides brief descriptions of the MELCOR code enhancement made between code revision number 14959and 18019. Revision 14959 represents the previous official code release; therefore, the modeling features described within this document are provided to assist users that update to the newest official MELCOR code release, 18019. Along with the newly updated MELCOR Users Guide and Reference Manual, users are aware and able to assess the new capabilities for their modeling and analysis applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

The Fracture Toughness of Nuclear Graphite Grades, Revision 1

New measurements of graphite mode I critical stress intensity factor, K Ic (commonly referred to as the fracture toughness) and the mode II critical shear stress intensity, K IIc , are reported and compared with prior data for K Ic and K IIc . The new data are for graphite grades PCEA, IG-110 and 2114. Variations of K Ic and acoustic emission (AE) data with graphite texture are reported and discussed. The Codes and Standards applications of fracture toughness, K Ic , data are also discussed. A specified minimum value for nuclear graphite K Ic is recommended.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

ORNL Neutron Cross Section Measurements of 90 Zr NCSP ND-1 Task [Slides]

Luiz Leal's idea to integrate the experimental uncertainties, covariances, and sensitivities in this workflow and centralize the process at one location led to the creation of a cohesive tool for data testing. It has been practiced for more than two decades by the NCS and ND groups at ORNL. The development and use of the ORNL codes TSUNAMI and TSURFER within SCALE made it possible. Today, it has been extended to SAMINT which enables the coupling of differential and integral data evaluation in a continuous-energy framework.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Grid Optimization Competition Challenge 3 Problem Formulation

This report contains the problem formulation for the Grid Optimization (GO) Competition Challenge 3. The Grid Optimization Competition is run by a team of researchers from a number of organizations, including the sponsor Advanced Research Projects Agency - Energy (ARPA-E), lead organization Pacific Northwest National Laboratory (PNNL), and technical contributors from Los Alamos National Laboratory (LANL), National Renewable Energy Laboratory (NREL), Texas A&M University (TAMU), Georgia Institute of Technology (GT), University of Wisconsin (UW), and others. The GO Competition poses challenge problems in the field of power grid management, invites entrants to develop solvers for these problems, invokes the solvers on a set of problem instances using common hardware, ranks the solvers according to their performance, and awards prizes according to the rankings. The overall goal of the GO Competition is to spur innovative research on high impact and computationally challenging problems in power grid management from initial development through commercial deployment. Complete information about the GO Competition can be found online at [2]. The webpage covers previous Challenges, rules, timeline, registration information, data formats, scoring methods, computational platform information, information on supported solvers and languages, sponsor information, frequently asked questions, administrator contact information, publicly available problem instances, computer code for reading and evaluating problem and solution data, a sandbox for testing solvers, and a solver submission interface, results, and publications.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

ORNL Neutron Cross Section Measurements of 90 Zr NCSP ND-1 Task [Slides]

Luiz Leal's idea to integrate the experimental uncertainties, covariances, and sensitivities in the ND workflow and to centralize the process at one location created a cohesive tool for data testing. It has been practiced for more than two decades by the NCS and ND groups at ORNL. The development and use of the ORNL codes TSUNAMI and TSURFER within SCALE made it possible. Today, it has been extended to SAMINT, which enables the coupling of differential and integral data evaluation in a continuous-energy framework.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Characterization of modeling and experimental data inconsistencies from burst testing for high-burnup commercial fuel rod applications

Fuel fragmentation, relocation, and dispersal are some of the largest issues remaining in the nuclear industry before rod-average burnup can be increased beyond 62 GWd/tU. The issue is primarily related to the potential for fuel to be dispersed into the reactor primary system, which may increase public risk. One way to prevent dispersal is to avoid cladding burst. The objective of this work is to support the high burnup safety case by evaluating cladding burst under high-burnup, full-length fuel rods and to identify uncertainties that could improve model predictions. The results of this analysis will evaluate realistic, prototypic loss-of-coolant accident (LOCA) conditions; support future cladding burst test designs; and inform the development of mechanistic material models. Realistic high-burnup operating conditions were implemented in the BISON fuel performance code to simulate steady-state and LOCA transient fuel rod evolution to the point at which cladding burst occurred. Parametric studies are performed to assess code response to changes in rod internal pressures and heating rates. Results were compared with simulated LOCA experiments to identify inconsistencies between commercial fuel rod analysis and experimental validation. The representative full-length fuel rod LOCA simulation results did not agree with cladding burst tests. Cladding burst tests indicated burst occurring well below (100–150 °C) those calculated in the full-length fuel rod LOCA analysis. Further investigation indicated that the cladding burst tests do not appear to be representative for full-length fuel rods. The inconsistency investigated in this work showed that the differences between the BISON simulation and the experiment's cladding burst conditions arise from an incomplete characterization of the cladding surface temperature, detailed rodlet characterization, lack of cladding strain measurements, and uncertainty in the cladding creep and failure models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Uncovering supramolecular chirality codes for the design of tunable biomaterials

In neurodegenerative diseases, polymorphism and supramolecular assembly of β-sheet amyloids are implicated in many different etiologies and may adopt either a left- or right-handed supramolecular chirality. Yet, the underlying principles of how sequence regulates supramolecular chirality remains unknown. Here, we characterize the sequence specificity of the central core of amyloid-β 42 and design derivatives which enable chirality inversion at biologically relevant temperatures. We further find that C-terminal modifications can tune the energy barrier of a left-to-right chiral inversion. Leveraging this design principle, we demonstrate how temperature-triggered chiral inversion of peptides hosting therapeutic payloads modulates the dosed release of an anticancer drug. These results suggest a generalizable approach for fine-tuning supramolecular chirality that can be applied in developing treatments to regulate amyloid morphology in neurodegeneration as well as in other disease states.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Reconstruction of beam parameters and betatron radiation spectra measured with a Compton spectrometer

The photon flux resulting from high-energy electron beam interactions with high-field systems, such as those found in the upcoming FACET-II experiments at the SLAC National Accelerator Laboratory, yields deep insight into the electron beam’s underlying dynamics during the interaction. However, extracting this information is an intricate process. To demonstrate how to approach this challenge using modern methods, this paper utilizes simulated data that models plasma wakefield acceleration-derived betatron radiation in experiments to determine reliable methods of reconstructing key beam and beam-plasma interaction properties. For betatron radiation measurements, translating the observed 200⁢ keV to 30⁢ MeV photon double-differential energy-angle spectra obtained from an advanced Compton spectrometer requires testing multiple methods to optimize the pipeline from its response to incident electron beam information. The paper compares maximum likelihood estimation and machine learning to refine the translation of photon spectra into precise electron beam metrics, such as spot size, energy, and emittance, enhancing the understanding of beam behavior within these dense, high-field environments. We also introduce machine learning and the expected maximization algorithm to reconstruct the primary photon spectrum, employing a multilayer neural network for regression analysis of the energy and angle spectra. With appropriate modifications, the advanced methods reproduce relevant incident beam parameters with high accuracy, even for beam sizes in the <10 μ⁢m range. This capacity is critical to understanding intense beam propagation and its optimization in plasma.

Beam code development & simulation techniques↗
DOE OSTI2024

PV Generation and Load Forecasting for Adjuntas PR Community Microgrids

Existing frameworks to forecast time-series photovoltaic (PV) output power and consumer load for microgrid operations and controls assume a near-continuous availability of real-time input features from the field assets such as PV inverters, energy meters, and weather station. These incoming data points are used to periodically retrain models and update forecast snapshots over a moving horizon window, be it one hour-ahead, one-day ahead, or one-week ahead. However, such frameworks are not resilient to disruptions in data availability caused by losses in communications between the field sensors and data loggers. Hence, there is a need for programs that assume no availability of real-time microgrid asset data and still make reliable forecasts that can be used for decision-making. Such programs would be apt to function in extreme weather events such as hurricanes and would use lightweight recursive time-series models to independently forecast solar irradiance and ambient temperature, then compute PV power from those forecasts, as well as independently forecast consumer load. The codebase performs forecasting for the scenario of when the microgrid does not have a reliable access to forecasts or real-time observations of solar irradiance (I) and ambient temperature (AT) and load (Load) to be able to adequately forecast, in real-time, the PV power production or a business' load. In this case, using historical values of PV power and load, a univariate forecasting of generation and consumption are respectively made. The use-case in particular has two sub-scenarios: one, a normal 7-day ahead forecast where the unavailability of real-time data is assumed due to infrastructure issues such as loss of communication or sensor maintenance or service downtimes. Whereas a hurricane-caused unavailability of real-time data requires a second model trained specifically on historical hurricane days to be able to capture the extreme day behavior of generation in particular, and load if applicable. A gradient boosted regression tree comprises an ensemble of additive models that map between the input of historical values (be it irradiance, temperature, or load) and their corresponding output forecasts of a given horizon such that the individual learner predictions are summed up over the total number of such learners in the ensemble to produce an aggregate forecast. A weighting mechanism is applied to the training data in each iteration, where actual and forecast values are compared to penalize incorrect forecasts by increasing the weight and reducing it to reward correct forecasts. The code's benefits are that it: (a) accounts for a contingency where communication loss renders newly measured real-time data unavailable for model tuning and snapshot updates; (b) presents blind forecasting that recursively determines the next time-step value in a horizon using the forecast of the same attribute from a prior step; and (c) employs lightweight models that, once trained, can reliably generalize for different horizons, which make them suitable for enhancing the resilience of field microgrids prone to extreme events that encounter disruptions to data availability.

Sundararajan, Aditya [Oak Ridge National Laborator↗
DOE OSTI2024

Final CRADA Report – NFE-21-08693

TAE Technologies is developing a magnetic fusion energy concept known as the beam-driven field-reversed configuration (FRC) with the ultimate goal of developing a reactor for commercial electricity production capable of burning aneutronic pB11 fuel. To achieve the high plasma temperatures this requires, auxiliary radiofrequency (RF) heating will likely be needed. High Harmonic Fast Wave (HHFW) heating has been identified as a candidate RF heating scheme to overcome the unique challenges posed to RF heating by the FRC, including the large distance from the plasma edge to the last closed flux surface and a magnetic field profile with strength decreasing from edge to core and reversing sign at a null point inside the plasma. The purpose of this project was to develop the experimental capabilities to test HHFW on TAE’s C-2W device through the design of a phased array antenna and accompanying matching network. The design was performed by ORNL and informed by experiments with a prototype four-strap phased antenna-array that was manufactured and installed on the LArge Plasma Device (LAPD) at UCLA and simulations conducted with the Petra-M code under the purview of a previous INFUSE grant. The ORNL team completed the conceptual design of the antenna and matching network which was then handed off to the TAE Mechanical Design team. The design was then iterated on to ensure changes to the mechanical design did not interfere with the RF performance. This process is now complete, and, with mechanical design in hand, TAE is proceeding with plans for final integration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Beam impedance and heating analysis of the diagnostic stripline

In this work the heat transfer analysis has been performed for the stripline kicker, design of which is used at many accelerator facilities as a part of the transverse bunch-by-bunch feed-back system, including NSLS-II. The experimental data of the electrode temperature have been collected from a special diagnostic stripline designed with six infrared view ports. The temperature analysis is based on the power loss calculations and the loss factor computed by the 3-D electromagnetic simulation code GdfidL. To estimate the power loss, we analyze the real part of the longitudinal impedance and the loss factor for two different regimes. We use the ANSYS finite element code to correlate the temperature with the power loss. Lambertson and Shafer formalisms of the longitudinal and transverse beam impedances are discussed for the simplified stripline geometries. The required characteristic impedances and the geometric factors are determined analytically for the simplified geometries and compared with the results for the designed stripline geometry using the 2-D POISSON code by solving the Laplace’s equation in two dimensions. Lambertson’s equations are compared with numerical results. We introduce the frequency range, within which we can apply the analytical approach.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Thermo-Hydro-Chemical Simulation of Mid-Ocean Ridge Hydrothermal Systems: Static 2D Models and Effects of Paleo-Seawater Chemistry

Decades of research have resulted in characterization of the ocean floor manifestations of mid-ocean ridge (MOR) hydrothermal systems, yet numerical models accounting for the connections between heat transfer, hydrology and geochemistry have been slow to develop. The Thermo-hydro-chemical code ToughReact can be used to describe the coupled effects of fluid flow, heat transfer, and fluid-rock chemical interactions that occur in MOR systems. We describe the results of 2-dimensional simulations of steady state flow in fractured diabase with mineral-fluid chemical reactions. Basal heating and specified permeability yield maximum temperature of 400°C. Total fluid flux and high fracture flow velocities are in accord with observations. Fluid chemistry, mineralogical changes and 87 Sr/ 86 Sr ratios can be compared to observations to assess and calibrate models. Simulated high temperature fracture fluids have Mg and SO 4 near zero, elevated Ca and 87 Sr/ 86 Sr of about 0.7040. Total alteration is 10%–50% for simple models of spreading. Anhydrite forms mainly near the base of the upwelling zone and results in substantial local fracture porosity reduction. A calibrated model is used to predict how Sr isotopes and other features of altered oceanic crust would be different in the Cretaceous (95 Ma) early Proterozoic (1,800 Ma) and Archean (3,800 Ma), when seawater may have had high Ca and Sr concentrations, lower pH, higher temperature, and lower Na, Mg, and SO 4 . The simulations are offered as a start on what ultimately may require a longer-term community effort to better understand the role of MOR thermo-hydro-chemical systems in Earth evolution.

58 GEOSCIENCES↗