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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

CICE on a C-grid: new momentum, stress, and transport schemes for CICEv6.5

Abstract. This article presents the C-grid implementation of the CICE sea ice model, including the C-grid discretization of the momentum equation, the boundary conditions (BCs), and the modifications to the code required to use the incremental remapping transport scheme. To validate the new C-grid implementation, many numerical experiments were conducted and compared to the B-grid solutions. In idealized experiments, the standard advection method (incremental remapping with C-grid velocities interpolated to the cell corners) leads to a checkerboard pattern. A modal analysis demonstrates that this computational noise originates from the spatial averaging of C-grid velocities at corners. The checkerboard pattern can be eliminated by adjusting the departure regions to match the divergence obtained from the solution of the momentum equation. We refer to this novel approach as the edge flux adjustment (EFA) method. The C-grid discretization with edge flux adjustment allows for transport in channels that are one grid cell wide – a capability that is not possible with the B-grid discretization nor with the C-grid and standard remapping advection. Simulation results match the predicted values of a novel analytical solution for one-grid-cell-wide channels.

Lemieux, Jean-François (ORCID:0000000320845759)↗
DOE OSTI2024

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram (j - α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j - α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

ASME Code Revisions to Incorporate 316H and Alloy 617 Viscoplastic Constitutive Models to Section III, Division 5 and Code Case N-898

This report provides a final status update on work to develop and implement two new constitutive models for 316H stainless steel and the Ni-based Alloy 617 in Nonmandatory Appendix Z to the Section III, Division 5, Subsection HB, Subpart B ASME Boiler & Pressure Vessel Code rules covering the design and construction of Class A high temperature nuclear reactor components. This report summarizes the objections of the overall project and provides the final versions of the constitutive models proposed for incorporation into the ASME Code. The report also provides an update on the balloting status at ASME of the two proposed constitutive models. As of the time of writing (July 2022) the models are on-track to be approved by ASME after the August 2022 Code Week. If so, this will mean the 316H model will be published in the 2023 edition of the Code and the A617 model available as part of a revised Code Case immediately

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simplified inelastic constitutive models for ASME Section III, Division 5 design by inelastic analysis

This report describes the development of simplified, universal constitutive model that captures the high temperature monotonic and cyclic behavior of a range of commonly-used high temperature materials. The goal of the work is to provide a simple, universal constitutive model to replace the current bespoke models for Grade 91, 316H, and Alloy 617 included in Nonmandatory Appendix HBB-Z of the ASME Boiler & Pressure Vessel Code, and to extend this model to cover Alloy 800H. We initiated this work in response to feedback from reactor vendors and other Code users requesting simplified models, compared to the current models, that are easier to implement and use in commercial finite element analysis software. This report describes the completion of this effort by developing a model to correct the defects in standard model forms presently used for high temperature material modeling, described in past work, developing and implementing new numerical methods to train this model against test data, and then actually training the model for the four materials. The report provides a complete mathematical description of the model along with the tabulated material coefficients for the four materials. The final step will be to formulate an ASME Code change to introduce the new models into the Code.

36 MATERIALS SCIENCE↗
DOE OSTI2024

VARI3D & PERSENT: Perturbation and Sensitivity Analysis (Revision 5)

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry-based transport code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Design Basis Document / Owner’s Technical Specification for Nitrate Salt Systems in CSP Projects (Final Technical Report)

The number of commercial coal, gas, and nuclear projects built over the past 100 years number in the thousands. As such, there is a large database, available to a wide range of commercial engineering contractors, on proven designs. In essence, unsuccessful designs have, through generations of iterations, been identified and then deleted from further consideration. In contrast, the number of commercial parabolic trough projects using nitrate salt for the thermal storage media is perhaps 60. Further, the number of commercial central receiver projects using nitrate salt as the working fluid is on the order of 20, including estimates for China. Given the relative immaturity of salt technology, and commercial pressures to successfully bid new solar projects into a mature electricity market, solar projects often promise more than has been delivered. The Design Basis Document / Owner’s Technical Specification is a first step in the iteration process. The report describes the successful features of commercial projects, outlines a range equipment and system failures in projects that didn’t operate as intended, and provides a draft set of design changes intended to correct the known problems. The product of the study is 3 volumes of technical material; one volume is on parabolic trough technologies; a second is on central receiver technologies, and the third is on potential design changes to parabolic trough and central receiver projects. The 3 volumes, which total some 590 pages, can be found at https://www.solardynllc.com/csp-plant-technologies. One of the principal topics in the report is the use of functional or prescriptive specifications. Functional specifications describe what the equipment needs to do, consistent with the minimum legal requirements of the local jurisdictions. The details of how this is to be accomplished is developed by the engineering contractor. Prescriptive specifications, which are developed by the Owner, prescribe to the engineering contractor how the functional requirements are to be met. This arrangement ensures that the favorable experience from a previous project is repeated. One example is the design code for the hot salt tank in central receiver projects. The closest design basis is API Standard 650 Welded Steel Tanks for Oil Storage. However, the maximum design temperature in API 650 is 260 °C. As such, solar projects have typically adopted a hybrid Code approach, in which allowable material stresses are taken from ASME Section II Materials. Further, since the tanks experience daily changes in temperature and in (static) pressure, and since portions of the tank can operate at stresses beyond the elastic range, the low cycle fatigue life of the tank is conducted using the rules of Section VIII Division 2. However, in a recent study by NREL, the principal damage mechanism was identified as creep rather than fatigue. Further, design stresses permitted under Section VIII Division 2, corresponding to a fatigue life of 30 years, result in projected creep lifetimes of only 2 to 5 years. An alternate design approach, prescribed by the Owner, would be based on Code sections intended for high temperature service in the creep regime. A candidate is Section III Division 5. Granted, this is a nuclear code section, and it’s use would not likely be mandated by local jurisdictions. However, the effects of creep have been deemed to be of sufficient importance that one nuclear project developer, and one central receiver project developer, have stipulated in the tank design specification that the equipment be designed to the requirements of Section III Division 5.

14 SOLAR ENERGY↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4

A translation of the PENELOPE physics subroutines to C++, designed as an extension of the GEANT4 toolkit, is presented. The Fortran code system PENELOPE performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. PENELOPE implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of GEANT4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original PENELOPE programs. The combined code, named PENG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of GEANT4.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Evaluation of Increased Peak Temperatures for Spent Fuel Cladding Performance during Dry Storage (Rev.1)

NUREG-2214, “Managing Aging Processes in Storage (MAPS) Report,” (2019) documents the technical basis for age-related degradation mechanisms on spent nuclear fuel assemblies during dry storage operations. The technical bases have been reviewed, and hoop cladding stress calculations were conducted using the Fuel Analysis under Steady-state and Transients (FAST) code at elevated normal condition temperatures of 425°C and 450°C, and at an off-normal and accident conditions temperature of 600°C. The NUREG-2214 conclusions are substantiated by the referenced literature. The calculations demonstrate that the aging mechanisms primarily driven by hoop stress such as hydride reorientation, delayed hydride cracking, thermal and athermal (low-temperature) creep, and localized mechanical overload are either not credible or do not compromise the fuel assembly’s performance for increased peak normal condition temperatures of 425°C or 450°C during the up to 60-year dry storage period. Other aging mechanisms considered are also not credible for compromising the fuel assembly’s performance during 60-year dry storage with a peak temperature of 400ºC as specified in ISG – 11 Rev. 3, and the same is indicated for the increased peak normal condition temperatures of 425°C or 450°C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Fuel Performance Evaluation of THOR-C Experiments

The Temperature Heatsink Overpower Response Commissioning (THOR-C) and THOR-Metal (THOR-M) experiments will be performed as part of an ongoing project for testing sodium fast reactor fuels with the Japan Atomic Energy Agency (JAEA). The THOR-C experiments consist of fresh metallic fuel pins and have been analyzed using the ABAQUS, Ansys codes and the BISON fuel performance code. THOR-M-Loss of Flow-1 (THOR-M-LOF-1) is designed to test an EBR-II irradiated fuel pin under LOF conditions. Simulation of the THOR-MLOF-1 experiment required first simulating the base irradiation of the fuel pin in EBR-II. MFUEL module of SAS4A/SASSYS-1 [1] is a physics-based metallic fuel performance model applicable to the normal operation, transient scenarios and fuel failure modeling including scenarios with bulk fuel melting. The model has been validated using EBR-II normal operation, separate effect transient tests as well as TREAT M-Series transient tests [2]. In this study, MFUEL models has been utilized together with a new capsule heat transfer model developed in this project. The new heat transfer model was necessary due to (1) significant amount of heat losses that required 2D heat transfer, (2) the presence of a titanium heat sink, rejecting a significant amount of heat, and (3) stagnant coolant conditions, which are inconsistent with SAS4A/SASSYS-1 (SAS) heat transfer model. Updates to SAS4A/SASSYS-1 and MFUEL has been described below, followed by a preliminary validation effort using the results from THOR-C-2 fresh fuel capsule experiment. A previous study for THOR-C-2 analysis using BISON code is also utilized in this study to model this test [3]. [1] D. O’Grady, A. J. Brunett, L. Ibarra, A. Karahan, T. Kim, T. S. Sumner, R. Thomas, T. H. Fanning, “The SAS4A/SASSYS-2 Version 5.7 Safety Analysis Code System,” Argonne National Laboratory,ANL/NSE-SAS/5.7, (2023). [2] A. Karahan, T. Kim, T. Fanning, D. O’Grady, “Validation of MFUEL Metal Fuel Performance Models of SAS4A/SASSYS-1,” Argonne National Laboratory, ANL/NSE-23/11, (2023). [3] M. Mihelish, A. Zabriskie, K. Paaren, P. Medvedev, C. Jensen, “Fuel Performance Predictions for the TREAT THOR-C Experiments,” Idaho National Laboratory, INL/RPT-23-73397, Revision 0, (2023)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Neotropical termite microbiomes as sources of novel plant cell wall degrading enzymes

In this study, we used shotgun metagenomic sequencing to characterise the microbial metabolic potential for lignocellulose transformation in the gut of two colonies of Argentine higher termite species with different feeding habits, Cortaritermes fulviceps and Nasutitermes aquilinus. Our goal was to assess the microbial community compositions and metabolic capacity, and to identify genes involved in lignocellulose degradation. Individuals from both termite species contained the same five dominant bacterial phyla (Spirochaetes, Firmicutes, Proteobacteria, Fibrobacteres and Bacteroidetes) although with different relative abundances. However, detected functional capacity varied, with C. fulviceps (a grass-wood-feeder) gut microbiome samples containing more genes related to amino acid metabolism, whereas N. aquilinus (a wood-feeder) gut microbiome samples were enriched in genes involved in carbohydrate metabolism and cellulose degradation. The C. fulviceps gut microbiome was enriched specifically in genes coding for debranching- and oligosaccharide-degrading enzymes. These findings suggest an association between the primary food source and the predicted categories of the enzymes present in the gut microbiomes of each species. To further investigate the termite microbiomes as sources of biotechnologically relevant glycosyl hydrolases, a putative GH10 endo-β-1,4-xylanase, Xyl10E, was cloned and expressed in Escherichia coli. Functional analysis of the recombinant metagenome-derived enzyme showed high specificity towards beechwood xylan (288.1 IU/mg), with the optimum activity at 50 °C and a pH-activity range from 5 to 10. These characteristics suggest that Xy110E may be a promising candidate for further development in lignocellulose deconstruction applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Gapped lineon and fracton models on graphs

We introduce a $\mathbb{Z}_N$ stabilizer code that can be defined on any spatial lattice of the form $\Gamma\times C_{L_z}$, where $\Gamma$ is a general graph. We also present the low-energy limit of this stabilizer code as a Euclidean lattice action, which we refer to as the anisotropic $\mathbb{Z}_N$ Laplacian model. It is gapped, robust (i.e., stable under small deformations), and has lineons. Its ground state degeneracy (GSD) is expressed in terms of a "mod $N$-reduction" of the Jacobian group of the graph $\Gamma$. In the special case when space is an $L\times L\times L_z$ cubic lattice, the logarithm of the GSD depends on $L$ in an erratic way and grows no faster than $O(L)$. Furthermore, we also discuss another gapped model, the $\mathbb{Z}_N$ Laplacian model, which can be defined on any graph. It has fractons and a similarly strange GSD.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗