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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of MHD instabilities with fluid runaway electron model in M3D- C 1

Runaway electrons are generated in a tokamak during the start up, during normal operation and during a plasma disruption. During a disruption, runaway electrons can be accelerated to high energies, potentially damaging the first wall. To predict the consequences of runaway generation during a disruption, it is necessary to consider resonant interactions of runaways with the bulk plasma. Here we consider the interactions of runaways on low mode number tearing modes. We have developed a fluid runaway electron model for the 3D MHD code M3D-C1[Jardin,et al. J Comput. Sci Discovery 6 014002 (2012)]. To benchmark, we have reproduced the MHD linear tearing mode results (with runaway electrons) in a circular cylinder presented in previous analytic studies[[Helander, P., et al, Phys. Plasmas 14 144102 (2007)] and extended here with a numerical eigenvalue calculation. Furthermore, we find that the low mode number tearing mode has a rotation caused by the MHD - runaways interaction, and the toroidal current scale length is much smaller with runaways than that for without and decreases as the runaway speed increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Urban microclimate and its impact on building performance: A case study of San Francisco

Urban microclimate exerts an increasing influence on urban buildings, energy, and sustainability. This study uses 10-year measured hourly weather data at 27 sites in San Francisco, California, to (1) analyze and visualize the urban microclimate patterns and urban heat island effect; (2) simulate annual energy use and peak electricity demand of typical large office buildings and large hotels to investigate the influence of urban microclimate on building performance; (3) simulate indoor air temperature of a single-family house without air-conditioning during the record three-day heatwave of 2017, to quantify the divergence of climate resilience due to urban microclimate effect. Results show significant microclimate effects in San Francisco with up to 11 °C outdoor air temperature difference between the coastal and downtown areas on September 1, 2017, during the record three-day heatwave. The simulated energy results of the prototype large office and large hotel buildings using the 2017 weather data show over 100% difference in annual heating energy use and 65% difference in annual cooling energy use across different stations; as well as up to 30% difference in peak cooling electricity demand. The impacts on annual site or source energy use are minimal (less than 5%) as cooling and heating in a mild climate are a relatively small portion of overall building energy use in San Francisco. Results also show the microclimate effects influence indoor air temperature of unconditioned homes by up to 5 °C. Newer buildings and homes are much less affected by microclimate effects due to more stringent performance requirements of the building envelope and energy systems. Furthermore, these findings inform that San Francisco microclimate variations should be considered in urban energy planning, building energy codes and standards, as well as heat resilience policymaking.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Steady-State Thermal-Hydraulic Analysis of the LEU-Fueled Dalat Nuclear Research Reactor

This paper presents results of steady-state thermal-hydraulic analysis for the designed working core of the Dalat Nuclear Research Reactor (DNRR) using the PLTEMP/ANL code. The core was designed to be loaded with 92 low-enriched uranium (LEU) VVR-M2 fuel bundles (FBs) and 12 beryllium rods surrounding a neutron trap at the core center, for replacement of the previous core with 104 high-enriched uranium (HEU) VVR-M2 FBs. Before using this code for thermohydraulic analysis of the designed LEU working core, it was validated by comparing calculation results with experimental data collected from the HEU working core of the DNRR. The discrepancy between calculated results and measured data was at the maximum about 0.8°C and 1.5°C of fuel cladding and outlet coolant temperatures, respectively. In the design calculation, thermohydraulic safety was confirmed through evaluation of the fuel cladding and coolant temperatures, as well as of other safety parameters such as Departure from Nucleate Boiling Ratio (DNBR) and Onset of Nucleate Boiling Ratio (ONBR). The calculation results showed that, in normal operation conditions at full nominal thermal power of 500 kW without uncertainty parameters, the maximum fuel cladding temperature of the hottest FB was about 90.4°C, which is lower than its limit value of 103°C, the minimum DNBR was 32.0, which is much higher than the recommended value of 1.5, and the minimum ONBR was 1.43, which is higher than the recommended value of 1.4 for VVR-M2 LEU fuel type. When the global and local hot channel factors were taken into account, the maximum temperature of fuel cladding at the hottest FB was about 98.4 °C, for global only, and 114.3°C, for global together with local hot channel factors. The calculation results confirm the safety operation of the designed LEU core loaded with 92 fresh VVR-M2 FBs.

Nguyen, Kien-Cuong↗
DOE OSTI2025

Real Vector Framework

SAND2025-11463O The Real Vector Framework (RVF) is a modern and flexible C++ vector math library for developing scientific computing software that involves vector computations. RVF allows an opt-in approach to functionality that parallels the familiar base-class and override structures of object-oriented programming. Users can reuse and customize the code without inheritance entanglements and dynamic dispatch, while enabling seamless interoperability between diverse container types. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

von Winckel, Gregory [Sandia National Lab. (SNL-CA↗
DOE OSTI2025

An update to δ-ray production in MCNP6

As originally introduced in MCNP6.2, the δ-ray production subroutine had a number of flaws related to treatment of particle or nuclear spin contributions to the production cross section. A recent update by the author to the MCNP6 code corrects these flaws. This report documents those changes for the benefit of MCNP6 users who use the delta ray production capability. Additionally, this report aims to provide a more comprehensive documentation of the delta ray production, building on the prior documentation by C. Anderson and coworkers. On this latter point, the goal is to save future MCNP developers from a visit to the library to dig out ancient texts from 1952, or at least to provide useful points of reference should such a visit be necessary.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Quantitative comparisons of electron-scale turbulence measurements in NSTX via synthetic diagnostics for high-k scattering.

Two synthetic diagnostics are implemented for the high-k scattering system in NSTX [D. R. Smith, E. Mazzucato, W. Lee, H. K. Park, C. W. Domier, and N. C. Luhmann Jr., Rev. Sci. Instrum., 79, 123501 (2008)] allowing direct comparisons between synthetic and experimentally detected frequency and wavenumber spectra of electron-scale turbulence fluctuations. Synthetic diagnostics were formulated in real-space and in wavenumber space, and were deployed to realistic electron-scale simulations carried out with the GYRO code [J. Candy and R. E. Waltz, J. Comput. Phys. 186 (2003) 545]. A highly unstable ETG regime (electron temperature gradient mode) in a modest beta NSTX NBI-heated H-mode discharge was chosen for the analysis. Mapping the measured wavenumbers to field aligned coordinates showed that the high-k system is sensitive to fluctuations that are closer to the spectral peak in the density fluctuation wavenumber spectrum (streamers) than originally predicted. Analysis of synthetic spectra show that the frequency response of the detected fluctuations is dominated by Doppler shift and insensitive to the turbulence drive. The shape of the high-k density fluctuation wavenumber spectrum is sensitive to the ETG turbulence drive conditions, and can be reproduced in a sensitivity scan of the most pertinent turbulent drive terms in the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Management of Risks Associated with Application of Novel Materials in Novel Operating Environments in Novel Reactor Designs

There is currently no widely agreed, detailed general method for licensing a novel plant incorporating novel materials (or materials being deployed in novel environments); in many such situations, there are no directly applicable engineering code cases for decision-makers (including regulators) to rely on. This paper discusses a framework for solving this problem that is based on the Reliability and Integrity Management (RIM) approach delineated in ASME BPVC Section XI Division 2. NRC Regulatory Guide 1.246, Rev. 0, endorses, with conditions, the subject portion of the ASME Code. The proposed framework is meant to support development of a licensing case by addressing certain technical challenges. The framework discussed here is compatible with the Licensing Modernization Project, but applying it in a specific case will call for advances in the state of practice, if not the state of the art. The RIM approach calls for applicants to (a) allocate reliability targets to plant structures, systems, and components (SSCs), (b) show that they are able to relate the currently observed physical condition of each SSC in the program to its failure probability well enough to determine whether the target reliability allocations are being satisfied, allowing for uncertainty related to the novelty of the materials/designs/operating environments, and (c) be able to demonstrate that the proposed program of surveillances will reliably detect unacceptable degradation of an SSC before SSC failure occurs. These challenges are discussed in the paper, and a potentially applicable modeling approach based on cumulative damage rather than failure rates is briefly illustrated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Simulation of multi-shell fullerenes using Machine-Learning Gaussian Approximation Potential

Multi-shell fullerenes ”buckyonions ” were simulated, starting from initially random configurations, using a density-functional-theory (DFT)-trained machine-learning carbon potential within the Gaussian Approximation Potential (GAP) Framework [Volker L. Deringer and Gábor Csányi, Phys. Rev. B 95, 094203 (2017)]. Fullerenes formed from seven different system sizes, ranging from 60 ~ 3774 atoms, were considered. The buckyonions are formed by clustering and layering starting from the outermost shell and proceeding inward. Inter-shell cohesion is partly due to interaction between delocalized π electrons protruding into the gallery. The energies of the models were validated ex post facto using density functional codes, VASP and SIESTA , revealing an energy difference within the range of 0.02 - 0.08 eV/atom after conjugate gradient energy convergence of the models was achieved with both methods.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Putting Humpty Dumpty Back Together Again: What Does Protein Quantification Mean in Bottom-Up Proteomics?

Bottom-up proteomics provides peptide measurements and has been invaluable for moving proteomics into large-scale analyses. Commonly, a single quantitative value is reported for each protein-coding gene by aggregating peptide quantities into protein groups following protein inference or parsimony. However, given the complexity of both RNA splicing and post-translational protein modification, it is overly simplistic to assume that all peptides that map to a singular protein-coding gene will demonstrate the same quantitative response. Here, by assuming that all peptides from a protein-coding sequence are representative of the same protein, we may miss the discovery of important biological differences. To capture the contributions of existing proteoforms, we need to reconsider the practice of aggregating protein values to a single quantity per protein-coding gene.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Evaluation of the excitation spectra with diffusion Monte Carlo on an auxiliary bosonic ground state

We aim to improve upon the variational Monte Carlo (VMC) approach for excitations replacing the Jastrow factor by an auxiliary bosonic (AB) ground state and multiplying it by a fermionic component factor. The instantaneous change in imaginary time of an arbitrary excitation in the original interacting fermionic system is obtained by measuring observables via the ground-state distribution of walkers of an AB system that is subject to an auxiliary effective potential. The effective potential is used to (i) drive the AB system’s ground-state configuration space toward the configuration space of the excitations of the original fermionic system and (ii) subtract from a diffusion Monte Carlo (DMC) calculation contributions that can be included in conventional approximations, such as mean-field and configuration interaction (CI) methods. In this novel approach, the AB ground state is treated statistically in DMC, whereas the fermionic component of the original system is expanded in a basis. The excitation energies of the fermionic eigenstates are obtained by sampling a fermion–boson coupling term on the AB ground state. We show that this approach can take advantage of and correct for approximate eigenstates obtained via mean-field calculations or truncated interactions. We demonstrate that the AB ground-state factor incorporates the correlations missed by standard Jastrow factors, further reducing basis truncation errors. Relevant parts of the theory have been tested in soluble model systems and exhibit excellent agreement with exact analytical data and CI and VMC approaches. In particular, for limited basis set expansions and sufficient statistics, AB approaches outperform CI and VMC in terms of basis size for the same systems. The implementation of this method in current codes, despite being demanding, will be facilitated by reusing procedures already developed for calculating ground-state properties with DMC and excitations with VMC.

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

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗
DOE OSTI2020

Neural-network accelerated coupled core-pedestal simulations with self-consistent transport of impurities and compatible with ITER IMAS

pedestal structure, current profile, and plasma equilibrium physics has been developed and tested against a DIII-D discharge. Here, key features of the achieved core-pedestal coupled workflow are its ability to account for the transport of impurities in the plasma self-consistently, as well as its use of machine learning accelerated models for the pedestal structure and for the turbulent transport physics. Notably, the coupled workflow is implemented within the OMFIT framework, and makes use of the ITER integrated modeling and analysis suite (IMAS) data structure for exchanging data among the physics codes that are involved in the simulations. Such technical advance has been facilitated by the development of a new numerical library named OMAS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

CLMM : a LSST-DESC cluster weak lensing mass modeling library for cosmology

ABSTRACT We present the v1.0 release of CLMM, an open source python library for the estimation of the weak lensing masses of clusters of galaxies. CLMM is designed as a stand-alone toolkit of building blocks to enable end-to-end analysis pipeline validation for upcoming cluster cosmology analyses such as the ones that will be performed by the Vera C. Rubin Legacy Survey of Space and Time-Dark Energy Science Collaboration (LSST-DESC). Its purpose is to serve as a flexible, easy-to-install, and easy-to-use interface for both weak lensing simulators and observers and can be applied to real and mock data to study the systematics affecting weak lensing mass reconstruction. At the core of CLMM are routines to model the weak lensing shear signal given the underlying mass distribution of galaxy clusters and a set of data operations to prepare the corresponding data vectors. The theoretical predictions rely on existing software, used as backends in the code, that have been thoroughly tested and cross-checked. Combined theoretical predictions and data can be used to constrain the mass distribution of galaxy clusters as demonstrated in a suite of example Jupyter Notebooks shipped with the software and also available in the extensive online documentation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A chromosome-scale reference genome of trifoliate orange ( Poncirus trifoliata ) provides insights into disease resistance, cold tolerance and genome evolution in Citrus

Trifoliate orange (Poncirus trifoliata), a deciduous close relative of evergreen Citrus, has important traits for citrus production, including tolerance/resistance to citrus greening disease (Huanglongbing, HLB) and other major diseases, and cold tolerance. It has been one of the most important rootstocks, and one of the most valuable sources of resistance and tolerance genes for citrus. Here we present a high-quality, chromosomescale genome assembly of P. trifoliata. The 264.9-Mb assembly contains nine chromosomal pseudomolecules with 25 538 protein-coding genes, covering 97.2% of the estimated gene space. Comparative analyses of P. trifoliata and nine Citrus genomes revealed 605 species-specific genes and six rapidly evolving gene families in the P. trifoliata genome. Poncirus trifoliata has evolved specific adaptation in the C-repeat/ DREB binding factor (CBF)-dependent and CBF-independent cold signaling pathways to tolerate cold. We identified candidate genes within quantitative trait loci for HLB tolerance, and at the loci for resistance to citrus tristeza virus and citrus nematode. Genetic diversity analysis of Poncirus accessions and Poncirus/ Citrus hybrids shows a narrow genetic base in the US germplasm collection, and points to the importance of collecting and preserving more natural genetic variation. Two phenotypically divergent Poncirus accessions are found to be clonally related, supporting a previous conjecture that dwarf Flying Dragon originated as a mutant of a non-dwarfing type. The high-quality genome reveals features and evolutionary insights of Poncirus, and it will serve as a valuable resource for genetic, genomic and molecular research and manipulation in citrus.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Carbon–Oxygen Classical Novae Are Galactic 7 Li Producers as well as Potential Supernova Ia Progenitors

Here, we report on studies of classical nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on carbon–oxygen (CO) white dwarfs (WDs). We vary both the mass of the WD (from 0.6 M ⊙ to 1.35 M ⊙ ) and the composition of the accreted material. Our simulations are guided by the results of multidimensional studies of TNRs in WDs, which find that sufficient mixing with WD core material occurs after the TNR is well underway, and levels of enrichment are reached that agree with observations of CN ejecta abundances. We use NOVA (our one-dimensional hydrodynamic code) to accrete solar matter until the TNR is ongoing and then switch to a mixed composition (either 25% WD material and 75% solar or 50% WD material and 50% solar). Because the amount of accreted material is inversely proportional to the initial 12 C abundance, by first accreting solar matter the amount of material taking part in the outburst is larger than in those simulations where we assume a mixed composition from the beginning. Our results show large enrichments of 7 Be in the ejected gases, implying that CO CNe may be responsible for a significant fraction (~100 M ⊙ ) of the 7 Li in the galaxy (~1000 M ⊙ ). Although the ejected gases are enriched in WD material, the WDs in these simulations eject less material than they accrete. We predict that the WD is growing in mass as a consequence of the accretion–outburst–accretion cycle, and CO CNe may be an important channel for SN Ia progenitors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Simulation of Peach Bottom turbine trip transient with RETRAN-3D using BXFGEN

The Peach Bottom 2 Turbine Trip transient is used to evaluate the capabilities and accuracy of coupled 3-D neutronics/thermal hydraulics best estimate calculations for support of safety analysis. Selected results, obtained with RETRAN-3D, are compared to available measurements and benchmark results. The comparisons are analyzed and the application of RETRAN-3D for a best estimate analysis approach within the safety analysis methodology is discussed. Overall, both the thermal-hydraulic and neutronic results show a good comparison to the experimental data, and the model is considered to accurately predict test data. In particular, the initial acoustic response of the main steam line and the subsequent main pressurization stages are well-predicted. During the main pressurization stage, the timing and magnitude of the in-core local pressurization and associated void compression are properly predicted based upon the agreement between local power range monitors (LPRM) predictions upwards, but the shift is slightly underpredicted. This indicates the methodology is well-posed for simulating this class of transient response.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cryogenic tracer irradiation facility at the university of Texas at Austin

We report on a cryogenic gas irradiation facility in the 1.1 MW TRIGA reactor at The University of Texas at Austin. The system was designed to produce radioactive xenon and argon for environmental studies, and it can be applied to produce other gaseous radiotracers. The system design includes modeling accident scenarios to ensure operation does not risk damage to the reactor and minimizes risk of release of radioactive material. In conclusion, the SCALE code was used to model a 30-day irradiation of one liter of Xe-126 and the predicted activity is 0.9 Ci of Xe-127 at irradiation end.

Ar-37↗