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 577 records · Page 32

Integrated modeling of advanced tokamaks for tearing mode avoidance on DIII-D

Understanding the physics of the onset and evolution of tearing modes (TMs) in tokamak plasmas is crucial for high-performance steady-state operations. The onset condition of (m,n) = (2,1) tearing stability, $Δ′$ > $Δ_c^′$> 0, has been studied in DIII-D steady-state hybrid discharges with accurate equilibrium reconstruction and well-measured plasma profiles. Here, m is the poloidal mode number and n is the toroidal mode number. The onset of n = 1 TMs is carefully identified using the Mirnov array to determine in-phase signals for the onset time and using electron cyclotron emission measurements to trace back from the saturated island width in the experiments for both initial island width and its growth rate. The tearing stability index $Δ′$ at TM onset is calculated in multiple ways, using the modified Rutherford equation and using the PEST-III or resistive DCON magnetohydrodynamic codes, and the results are in reasonable agreement with each other. The calculated $Δ′$ is compared with an analytical formula of the tearing stability threshold $Δ_c^′$ to determine the condition for the mode onset, $Δ′$ > $Δ_c^′$. Finally, the local gradient effect of the plasma current and the pressure profiles on the tearing stability is investigated in terms of the difference between $Δ′$ and $Δ_c^′$ through a series of modeled equilibria with the parametric variations. A database of experimental mode instability onsets is compatible with the simulated regions of instability obtained from the parametric variations.

Linear stability analysis↗

PARFUME/BISON Fission Product Release Predictions versus AGR-3/4 Heating Test Measurements

The fuel performance modeling codes PARFUME (PARticle FUel ModEl) and BISON were used to predict the release of fission products silver, cesium, and strontium from as-irradiated fuel compacts containing tristructural isotropic (TRISO) coated particles during heating tests post irradiation. The AGR-3/4 fuel compacts were irradiated as part of the third and fourth series of planned experiments to support the Advanced Gas Reactor (AGR) program. The heating tests were conducted at temperatures between 1200°C and 1700°C to simulate reactor accident conditions. The measured fission product release fractions from the heating tests were compared to modeling predictions calculated by PARFUME and BISON to evaluate how the codes compare to experimental results. Comparisons between the experimental measured fission product release fractions from silver, cesium and strontium indicate that both modeling codes overpredict the fission product release fractions demonstrating that the diffusivities used in the codes are overestimated. This results in a conservative estimate predicted by the codes when evaluating the fission product release relative to experimental data as it pertains to silver, cesium, and strontium.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Program Verification for Extreme-Scale Applications (Final Scientific/Technical Report)

The project seeks to develop tools and techniques to help software developers verify the correctness and accuracy of their code. The target domain is scientific software of the kind widely used and developed in the Department of Energy research community, with a particular focus on "extreme scale" programs - those that are expected to involve possibly millions of parallel threads of execution. The report covers the University of Delaware contribution to the collaborative project. The project had a number of successful outcomes, especially regarding the development and extension of the CIVL software verification framework. CIVL is a verification tool for C or Fortran programs that use MPI, OpenMP, CUDA, and/or Pthreads for parallelization. CIVL went through vast improvements and extensions, and was successfully applied to a number of challenging codes. It found a subtle bug in the Devito PDE framework. It was able to verify the functional correctness of a conjugate gradient solver using a novel probabilistic technique with vanishingly small chance of error. CIVL was used very successfully in verification competitions, and the CIVL solutions to the competition challenges were presented and published. CIVL was also used successfully by other researchers on a computational chemistry kernel.

97 MATHEMATICS AND COMPUTING↗

Introduction to an MCGIDI Mini-App and Performance Comparisons with XSBench

In high-performance Monte Carlo (MC) radiation transport codes, cross section lookups often account for the majority of computational expens. In response to this observation, a number of nuclear data mini-apps have been developed to profile and optimize the lookup process. A review of these mini-apps can be found in Ref. In contrast to a fully fledged MC code, nuclear data mini-apps do not feature particle tracking; instead, the cross section lookup process is considered in isolation, a simplification that allows for easier implementation and evaluation of various lookup schemes and techniques. GIDI+, as described on its public repository1, is ...a collection of C++ libraries for accessing evaluated and processed nuclear data stored in the Generalized Nuclear Database Structure (GNDS). In addition to reading GNDS files, GIDI+ has functions to sum and collapse multi-group data as needed by deterministic transport codes, and to sample GNDS data as needed by Monte Carlo transport codes. Certain modes of Mercury, a LLNL-developed MC radiation transport code, utilize GIDI+ for cross section lookups. Specifically, Mercury uses the MCGIDI library contained within GIDI+ to perform lookups on GPUs. A mini-app exercising MCGIDI’s cross section lookup capabilities would prove useful for optimization, which in turn could improve the performance of client codes such as Mercury. Such a mini-app could also be used to assess the performance of MCGIDI compared to other nuclear data mini-apps. In this document, we introduce a nuclear data mini-app built with the MCGIDI library. The MCGIDI mini-app (MCGIDI-MA) was developed to closely resemble the methodology present in XSBench, a nuclear data mini-app developed at Argonne National Laboratory. The remainder of this section provides an overview on XSBench. Section 2 describes MCGIDI-MA and its usage. Section 3 contains performance results of MCGIDI-MA on LLNL’s Lassen and Quartz compute platforms, as well as comparisons to XSBench performance where appropriate. In Section 4, we present a summary of our findings.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Verification of Reduced Order Modeling Based Uncertainty/Sensitivity Estimator (ROMUSE)

This article presents a number of verification case studies for a recently developed sensitivity/uncertainty code package. The code package, ROMUSE (Reduced Order Modeling based Uncertainty/Sensitivity Estimator) is an effort to provide an analysis tool to be used in conjunction with reactor core simulators, in particular the Virtual Environment for Reactor Applications (VERA) core simulator. ROMUSE has been written in C++ and is currently capable of performing various types of parameter perturbations and associated sensitivity analysis, uncertainty quantification, surrogate model construction and subspace analysis. The current version 2.0 has the capability to interface with the Design Analysis Kit for Optimization and Terascale Applications (DAKOTA) code, which gives ROMUSE access to the various algorithms implemented within DAKOTA, most importantly model calibration. The verification study is performed via two basic problems and two reactor physics models. The first problem is used to verify the ROMUSE single physics gradient-based range finding algorithm capability using an abstract quadratic model. The second problem is the Brusselator problem, which is a coupled problem representative of multi-physics problems. This problem is used to test the capability of constructing surrogates via ROMUSE-DAKOTA. Finally, light water reactor pin cell and sodium-cooled fast reactor fuel assembly problems are simulated via SCALE 6.1 to test ROMUSE capability for uncertainty quantification and sensitivity analysis purposes.

97 MATHEMATICS AND COMPUTING↗

AMReX: a framework for block-structured adaptive mesh refinement

AMReX is a C++ software framework that supports the development of block-structured adaptive mesh refinement (AMR) algorithms for solving systems of partial differential equations (PDEs) with complex boundary conditions on current and emerging architectures. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving the local descriptions of different physical processes in complex multiphysics algorithms. Current AMReX-based application codes span a number of areas, including atmospheric modeling, astrophysics, combustion, cosmology, fluctuating hydrodynamics, multiphase flows, and particle accelerators. In particular, the AMReX-Astro GitHub repository holds a number of astrophysical modeling tools based on AMReX. The origins of AMReX trace back to the BoxLib software framework. AMReX supports a number of different time-stepping strategies and spatial discretizations. Solution strategies supported by AMReX range from level-by-level approaches (with or without subcycling in time) with multilevel synchronization to full-hierarchy approaches, and any combination thereof. User-defined kernels that operate on patches of data can be written in C++ or Fortran; there is also a Fortran-interface functionality which wraps the core C++ data structures and operations in Fortran wrappers so that an application code based on AMReX can be written entirely in Fortran.

97 MATHEMATICS AND COMPUTING↗

Data for “Fast-decaying plant litter enhances soil carbon in temperate forests, but not through microbial physiological traits”

This data package contains data and code used in the paper “Fast-decaying plant litter enhances soil carbon in temperate forests, but not through microbial physiological traits”. This paper details results from two studies: 1) a laboratory leaf litter incubation experiment (lab experiment) and 2) a multi-site observational field study (field study). Both studies were designed to test the relationships among litter quality, microbial physiological traits, and mineral-associated soil carbon (C). In the lab experiment, we incubated 16 temperate tree litters of differing chemical quality with isotopically distinct soil, measuring microbial physiological traits and the flow of litter-derived C into the mineral-associated soil C pool. In the field study, we sampled soils (0-5 cm) across six eastern US temperate forests, measuring microbial physiological traits, soil abiotic properties, and leaf litter chemistry. The analysis data are provided in two .csv files corresponding to either the lab experiment or field study. Both files contain data on litter chemistry, microbial physiological traits (growth and turnover rates, and carbon use efficiency [CUE]), and the mineral-associated soil C pool. The lab experiment dataset additionally contains litter-derived versus soil-derived soil C, respiration, and litter decomposition parameters. The field study dataset additionally contains site-level climatic information, ectomycorrhizal dominance of plots, and plot-level soil properties. Also provided is the R code and output reproducing the results in the related publication. Analyses were originally performed using R version 3.6.1 using the package "lavaan" and the packages listed on lines 5-7 in the file "data.R".

54 ENVIRONMENTAL SCIENCES↗

Standardized Protocol for Real-Time APIs as Required by Title 23 CFR 680.116(c)

Improving the ability of drivers to easily locate working and available chargers is key to improving the public charging experience. Electric vehicle charging providers who are recipients of federal funds through the National Electric Vehicle Infrastructure (NEVI) Formula Program, Charging and Fueling Infrastructure (CFI) Discretionary Grant Program, and other funding programs as identified under Title 23 of the U.S. Code must deploy and maintain an application programming interface (API) to access information about charging stations they operate.1 This includes information about individual charging ports, pricing, and availability in accordance with the Federal Highway Administration’s National Electric Vehicle Infrastructure Standards and Requirements, 23 CFR 680.116(c), herein referred to as the minimum standards (Federal Highway Administration 2023). Specifically outlined in the minimum standards, states and other designated recipients are required to ensure that charging station information including location, connector type, power level, real-time status, and real-time price to charge are available free of charge to third-party software developers through an API. These requirements are intended to enable effective communication with consumers about available charging stations and help consumers make informed decisions about trip planning, including when and where to charge. This document provides a standardized protocol for how to structure data, data update frequency, and practices for making the data required to be shared via API usable for improving public transparency and the customer experience. These are recommendations only and do not modify the Federal Highway Administration’s minimum standards in any way.

33 ADVANCED PROPULSION SYSTEMS↗

Optimization of Condensed Matter Physics Application with OpenMP Tasking Model

The Density Matrix Renormalization Group (DMRG++) is a condensed matter physics application used to study superconductivity properties of materials. It’s main computations consist of calculating hamiltonian matrix which requires sparse matrix-vector multiplications. This paper presents task-based parallelization and optimization strategies of the Hamiltonian algorithm. The algorithm is implemented as a mini-application in C++ and parallelized with OpenMP. The optimization leverages tasking features, such as dependencies or priorities included in the OpenMP standard 4.5. The code refactoring targets performance as much as programmability. The optimized version achieves a speedup of 8.0× with 8 threads and 20.5× with 40 threads on a Power9 computing node while reducing the memory consumption to 90 MB with respect to the original code, by adding less than ten OpenMP directives.

Criado, Joel↗

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Continuous-Energy ENDF/B-VIII.0 Cross Section and SCALE 6.2.4 Performance for Nuclear Criticality Safety Applications: 1 H, C, 58,60 Ni, 182,183,184,186 W, 235,238 U, 239 Pu

As part of the nuclear data evaluation and validation cycle, the ENDF/B-VIII.0 cross-section library released in 2018 requires testing to determine areas of improvement and deterioration. Previous work by the authors investigated the performance of 16 O, 56 Fe, and 63,65 Cu cross sections, with this study acting as an extension of the prior work. In addition to the isotopes and nuclear criticality safety benchmarks of interest to the prior work, benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook were selected for their k eff sensitivity to 1 H, C, 58,60 Ni, 182,183,184,186 W, 235,238U, or 239 Pu cross sections and were modeled in the SCALE code system maintained by Oak Ridge National Laboratory. In total, 253 benchmark configurations were selected for their sensitivities and modeled using SCALE 6.2.4 Criticality Safety Analysis Sequences (CSAS) continuous-energy Monte Carlo k eff calculations. This collection includes and expands upon the 99 benchmarks in the prior work. The AMPX-processed ENDF/B-VIII.0 library was decomposed into individual ENDF/B-VIII.0 datum libraries for each isotope of interest. Doing so allowed for the individual substitution of an ENDF/B-VIII.0 cross section in the place of ENDF/B-VII.1, determining isotope-specific effects of ENDF/B-VIII.0 relative to ENDF/B-VII.1. Full library calculations with entirely ENDF/B-VII.1 data or entirely ENDF/B-VIII.0 data were also executed. As a measure of performance, the average relative deviation was determined as the ratio of the deviation between calculated and experimental keff to the propagated calculational and experimental uncertainty. With calculated full library and isotope-specific ENDF/B-VIII.0 k eff ’s, an optimized combination of data libraries was estimated and confirmed with SCALE calculations. This showed that reverting 239 Pu, 58 Ni, 16 O, and 65 Cu cross sections to ENDF/B-VII.1 resulted in improved performance relative to the full ENDF/B-VIII.0 library. Across all 253 benchmarks, the average relative deviation was 1.29σ for the full ENDF/B-VII.1 library, 1.17σ for the full ENDF/B-VIII.0 library, and 0.97σ for the optimized combination. The reversion of 239 Pu, 58 Ni, 16 O, and 65 Cu cross sections to ENDF/B-VII.1 in the 99 benchmarks of the prior work resulted in further improved experimental agreement compared to the previously reported improvement from 16 O and 65 Cu alone. Therefore, it is suggested that applications with significant sensitivities to 239 Pu, 58 Ni, 16 O, and 65 Cu consider their choice of nuclear data library.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Structural basis for preferential binding of human TCF4 to DNA containing 5-carboxylcytosine

The psychiatric risk-associated transcription factor 4 (TCF4) is linked to schizophrenia. Rare TCF4 coding variants are found in individuals with Pitt-Hopkins syndrome—an intellectual disability and autism spectrum disorder. TCF4 contains a C-terminal basic-helix-loop-helix (bHLH) DNA binding domain which recognizes the enhancer-box (E-box) element 5'-CANNTG-3' (where N = any nucleotide). A subset of the TCF4-occupancy sites have the expanded consensus binding specificity 5'-C(A/G)-CANNTG-3', with an added outer Cp(A/G) dinucleotide; for example in the promoter for CNIH3, a gene involved in opioid dependence. In mammalian genomes, particularly brain, the CpG and CpA dinucleotides can be methylated at the 5-position of cytosine (5mC), and then may undergo successive oxidations to the 5-hydroxymethyl (5hmC), 5-formyl (5fC), and 5-carboxyl (5caC) forms. We find that, in the context of 5'- 0 CG- 1 CA- 2 CG- 3 TG-3'(where the numbers indicate successive dinucleotides), modification of the central E-box 2 CG has very little effect on TCF4 binding, E-box 1 CA modification has a negative influence on binding, while modification of the flanking 0 CG, particularly carboxylation, has a strong positive impact on TCF4 binding to DNA. Crystallization of TCF4 in complex with unmodified or 5caC-modified oligonucleotides revealed that the basic region of bHLH domain adopts multiple conformations, including an extended loop going through the DNA minor groove, or the N-terminal portion of a long helix binding in the DNA major groove. The different protein conformations enable arginine 576 (R576) to interact, respectively, with a thymine in the minor groove, a phosphate group of DNA backbone, or 5caC in the major groove. The Pitt-Hopkins syndrome mutations affect five arginine residues in the basic region, two of them (R569 and R576) involved in 5caC recognition. Our analyses indicate, and suggest a structural basis for, the preferential recognition of 5caC by a transcription factor centrally important in brain development.

59 BASIC BIOLOGICAL SCIENCES↗

Real-Time Interactive 4D-STEM Phase-Contrast Imaging From Electron Event Representation Data: Less computation with the right representation

The arrival of direct electron detectors (DED) with high frame-rates in the field of scanning transmission electron microscopy has enabled many experimental techniques that require collection of a full diffraction pattern at each scan position, a field which is subsumed under the name four dimensional-scanning transmission electron microscopy (4D-STEM). DED frame rates approaching 100 kHz require data transmission rates and data storage capabilities that exceed commonly available computing infrastructure. Current commercial DEDs allow the user to make compromises in pixel bit depth, detector binning or windowing to reduce the per-frame file size and allow higher frame rates. This change in detector specifications requires decisions to be made before data acquisition that may reduce or lose information that could have been advantageous during data analysis. The 4D Camera, a DED with 87 kHz frame-rate developed at Lawrence Berkeley National Laboratory, reduces the raw data to a linear-index encoded electron event representation (EER). Here we show with experimental data from the 4D Camera that linear-index encoded EER and its direct use in 4D-STEM phase contrast imaging methods enables real-time, interactive phase-contrast from large-area 4D-STEM datasets. Furthermore, we detail the computational complexity advantages of the EER and the necessary computational steps to achieve real-time interactive ptychography and center-of-mass differential phase contrast using commonly available hardware accelerators.

4D-STEM↗

Computer Science Research Needs for Parallel Discrete Event Simulation (PDES)

Historically, scientific computing efforts have demonstrated the clear need for, and effective use of, supercomputing with traditional time-stepped simulations. Nevertheless, there are several areas in the mission spaces of the U.S. Department of Energy and other agencies waiting to tap advanced computing research using a different, discrete event style of modeling, simulation, and analysis. These span a wide spectrum of applications including energy grid resilience, urban planning and policy, transportation science, building technologies, emergency response and planning, environmental impact analysis, computational epidemiology, Internet communications, cyber security, and cyber-physical systems, to name only a few. Even within traditional scientific applications, the role of discrete event modes of execution is increasing in the form of new event-based mathematical solvers such as quantized state integration methods and discrete-continuous hybrid system solvers. Co-design of advanced supercomputing hardware systems is another area that exploits discrete event simulation at its core for effective analyses. Complex systems, entity behaviors and interconnections play a significant role in all these applications, which are mapped to large-scale models with discrete event formulations. To make advancements in all the aforementioned scientific areas, many technical aspects need to be more thoroughly studied and deeply understood in parallel discrete event simulation (PDES). The unique dynamics inherent in a discrete event modeling approach, by their very nature, intersect and influence the entire stack of the computing system, including (a) the unique nature of the instruction sets exercised in PDES workloads without a predominance of high-precision floating point operations, (b) virtual time-constrained multi-threaded execution of many logical processes per processor, (c) extremely variable and difficult to predict network traffic characteristics, (d) interfaces and inter-dependencies with machine learning and artificial intelligence codes at higher software layers, and (e) highly challenging load balancing needs, especially in effectively accounting for accelerated/extremely heterogeneous computing in current and future high-performance computing systems. Efficient and accurate parallel execution of PDES workloads is also dominated by challenges in dealing with their asynchronous concurrency fundamentally present at the model level. Conservative synchronization, optimistic/speculative synchronization, and their hybrid schemes open new questions in fundamental computer science with respect to reversibility of computation and prediction (lookahead) of behaviors inherent within model codes. On the implementation front, there are relatively few scalable, general-purpose parallel discrete event simulators in the world, and even fewer have been studied on emerging hardware platforms. To enable scientific advances using PDES, the research needs in computer science must also be pursued and met in the intersection of the algorithmic and hardware-aware aspects of scalable PDES engines. This report is aimed at capturing a computer science-oriented view of this important area of research in PDES, presenting a sample of important applications with their inherent discrete event technology elements. Needs are outlined in core areas of parallel discrete event research as well as cross-cutting directions in computer science research that positively impact scientific advancements across several important application areas. A selection of priority research opportunities in advanced computing for PDES is identified to serve as reference for key research topics and their order of importance for scientific advancements.

97 MATHEMATICS AND COMPUTING↗

Results of the Phenix Dissymmetric Test Benchmark Exercise

In the context of the European H2020 project SESAME, a work package was dedicated to the validation of code coupling for the thermal-hydraulics of pool type reactors. To this end, CEA proposed a benchmark exercise of the so-called dissymmetric test performed at the end of life of the PHENIX reactor. The PHENIX reactor was an experimental reactor operated by CEA in France from 1974 until 2009. The knowledge and data acquired over the lifetime of this reactor are of paramount importance for the development of future sodium-cooled generation IV reactors. At the beginning of this test, the reactor is operated at nominal power (350 MWth) with 3 primary pumps and 4 intermediate heat exchangers. During this test, a pump trip occurs at the secondary side of one of the two loops (one loop provides coolant to two intermediate heat exchangers). The chain reaction is then stopped while the pumps remain operating at nominal speed. Dissymmetric flow patterns are then observed in the cold pool as hot sodium leaves the intermediate heat exchangers that are no longer cooled. The failed secondary pump is then restarted while the rotational speed of the other is reduced. Eventually, a stable homogeneous state is reached. This transient is challenging for both modelers and codes since this exercise tackles an industrial scale reactor featuring complex flows and interactions between different scales. The physical phenomena observed include thermal stratification, flow mixing and buoyancy effects. In order to accurately model these phenomena for a reasonable computational costs, the partners used System Thermal Hydraulics codes coupled with Computational Fluid Dynamics codes and one even included a subchannel code. In this paper, the PHENIX reactor and its behavior during the dissymmetric test are described. Then, a comparison of the different modeling options selected by the partners is offered. Finally, the results are compared to the experimental data and discussed and modelling recommendations are proposed.

code coupling↗

Data for "Which plant traits increase soil carbon sequestration? Empirical evidence from a long-term poplar genetic diversity trial"

This archive contains all data and code used by the following publication: Field, J. L., Sloan, B. P., Craig, M. E., Calloway, P., Ottinger, S. L., Mead, T., Abramoff, R. Z., Venegas, M. P., Chhetri, H. B., Haiby, K., Kalluri, U. C., Muchero, W., Schadt, C. W., & Mayes, M. A. (2025). Which plant traits increase soil carbon sequestration? Empirical evidence from a long-term poplar genetic diversity trial (p. 2025.02.17.638464). bioRxiv. https://doi.org/10.1101/2025.02.17.638464 Our analysis combined several soil and root data sets collected by Oak Ridge National Laboratory (ORNL) researchers/collaborators from the Clatskanie Poplar Common Garden in Clatskanie, OR by from 2009-2024. The raw data data files are located */02-data/01-raw/* which we harmonized using the codes in */01-codes/01-harmonize-clatskanie-data-pub.qmd*. The final processed data set used in the paper is found at */02-data/02-processed/clatskanie-c-fit-data.csv* and its columns are described in the table below.

Sloan, Brandon [ORNL] (ORCID:0000000316304271)↗