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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Holographic measurement and bulk teleportation

Holography has taught us that spacetime is emergent and its properties depend on the entanglement structure of the dual theory. In this paper, we describe how changes in the entanglement due to a local projective measurement (LPM) on a subregion A of the boundary theory modify the bulk dual spacetime. We find that LPMs destroy portions of the bulk geometry, yielding post-measurement bulk spacetimes dual to the complementary unmeasured region A c that are cut off by end-of-the-world branes. Using a bulk calculation in AdS 3 and tensor network models of holography (in particular, the HaPPY code and random tensor networks), we show that the portions of the bulk geometry that are preserved after the measurement depend on the size of A and the state we project onto. The post-measurement bulk dual to A c includes regions that were originally part of the entanglement wedge of A prior to measurement. This suggests that LPMs performed on a boundary subregion A teleport part of the bulk information originally encoded in A into the complementary region A c . In semiclassical holography an arbitrary amount of bulk information can be teleported in this way, while in tensor network models the teleported information is upper-bounded by the amount of entanglement shared between A and A c due to finite-N effects. When A is the union of two disjoint subregions, the measurement triggers an entangled/disentangled phase transition between the remaining two unmeasured subregions, corresponding to a connected/disconnected phase transition in the bulk description. Our results shed new light on the effects of measurement on the entanglement structure of holographic theories and give insight on how bulk information can be manipulated from the boundary theory. They could also represent a first step towards a holographic description of measurement-induced phase transitions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Temperature response and brine availability to heated boreholes in bedded salt

Abstract There is a growing need for disposal of high‐level nuclear waste. To reduce uncertainty associated with brine availability to repository excavations in salt formations, a collaboration between Sandia, Los Alamos, and Lawrence Berkeley National Laboratories is performing a series of borehole‐scale coupled process tests. Here, we report on the first round of the Brine Availability Test in Salt (BATS) project, a “shakedown” experiment called Phase 1s. Experimental testing included placing a resistive heater, a 260‐W radiative heater, and a 750‐W radiative heater within previously drilled horizontal boreholes at the Waste Isolation Pilot Plant (WIPP) while monitoring temperature and water inflow. The experiments successfully achieved the targeted temperature of 120 °C when using the 750‐W radiative heater. Simulations using FEHM (Finite Element Heat and Mass transfer code) and TOUGH‐FLAC (Transport Of Unsaturated Groundwater and Heat–Fast Lagrangian Analysis of Continua) were able to accurately predict the coupled thermo–hydro–mechanical–chemical response of salt, matching the observed temperature and brine production. Due to the extremely low permeability of salt, these systems take many years to reach steady state when perturbed by mining activities. Long‐term numerical simulations are used to develop the initial pressure and saturation conditions. The inclusion of a damaged rock zone with higher permeability around the borehole also affects the saturation and pressure distributions and plays an important role in dissipating the potential for thermal pressurization. Knowledge gained from this round of experimentation and simulation will be used to conduct the next BATS project experiment in new boreholes at WIPP.

58 GEOSCIENCES↗
DOE OSTI2022

Verifying raytracing/Fokker–Planck lower-hybrid current drive predictions with self-consistent full-wave/Fokker–Planck simulations

Raytracing/Fokker–Planck (FP) simulations used to model lower-hybrid current drive (LHCD) often fail to reproduce experimental results, particularly when LHCD is weakly damped. A proposed reason for this discrepancy is the lack of ‘full-wave’ effects, such as diffraction and interference, in raytracing simulations and the breakdown of the raytracing approximation. Previous studies of LHCD using non-Maxwellian full-wave/FP simulations have been performed, but these simulations were not self-consistent and enforced power conservation between the FP and full-wave code using a numerical rescaling factor. Here, we have created a fully self-consistent full-wave/FP model for LHCD that is automatically power conserving. This was accomplished by coupling an overhauled version of the non-Maxwellian TORLH full-wave solver and the CQL3D FP code using the Integrated Plasma Simulator. We performed converged full-wave/FP simulations of Alcator C-Mod discharges and compared them with raytracing. We found that excellent agreement in the power deposition profiles from raytracing and TORLH could be obtained, however, TORLH had somewhat lower current drive efficiency and broader power deposition profiles in some cases. This discrepancy appears to be a result of numerical limitations present in the TORLH model and a small amount of diffractional broadening of the TORLH wave spectrum. Our results suggest full-wave simulation of LHCD is likely not necessary as diffraction and interference represented only a small correction that could not account for the differences between simulations and experiment.

Physics↗
DOE OSTI2023

ELM sputter erosion modeling of a tungsten coated small angle slot divertor in DIII-D

We modeled plasma edge localized mode (ELM) sputter erosion for a Small Angle Slot divertor with a tungsten coated region (SAS-VW), designed for experiments in the DIII-D tokamak, and proposed for use in future advanced tokamaks. The simulations use a free-streaming, 1000 eV, C +6 and D +1 ELM impingement model, with SOLPS-ITER, ITMC-DYN, and REDEP/WBC code packages for background plasma, material response, and erosion/redeposition respectively. The results show ELM'ing plasma gross and net tungsten erosion fluxes of the mixed-material C/W surface peaking at the slot entrance region, and an order of magnitude higher than for non-ELMs. The per-pulse erosion, however, remains low, of order 0.5 nm, due to expected moderate ELM frequencies and duration in DIII-D. Here, the ELMs result in a ~25x higher peak sputtered W current leaving the divertor slot region, towards the core plasma, compared to the ELM-free plasma case. The time-integrated escape current, however, may not significantly affect core plasma high-Z contamination concerns, for a 1% ELM duty factor, but may be an issue for higher frequency ELMs. In general, the modeling results appear favorable for effective testing of the SAS-VW divertor in DIII-D, and extrapolation to innovative divertor designs in future ITER-like and DEMO fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Sphingobium lignivorans sp. nov., isolated from river sediment downstream of a paper mill

Here, a bacterial isolate, B1D3A T , was isolated from river sediment collected from the Hiwassee River near Calhoun, TN, by enrichment culturing with a model 5–5' lignin dimer, dehydrodivanillate, as its sole carbon source. B1D3A T was also shown to utilize several model lignin-derived monomers and dimers as sole carbon sources in a variety of minimal media. Cells were Gram-stain-negative, aerobic, motile, rod-shaped and formed yellow/cream-coloured colonies on rich agar. Optimal growth occurred at 30°C, pH 7–8, and in the absence of NaCl. The major fatty acids of B1D3A T were C 18:1 ω7c and C 17:1 ω6c. The predominant hydroxy fatty acids were C 14: 0 2-OH and C 15:0 2-OH. The polar lipid profile consisted of a mixture of phosphatidylethanolamine, phosphatidylglycerol, diphosphatidylglycerol, phosphatidyldimethylethanolamine and sphingoglycolipid. B1D3A T contained spermidine as the only major polyamine. The major isoprenoid quinone was Q-10 with minor amounts of Q-9 and Q-11. The genomic DNA G+C content of B1D3A T was 65.6mol%. Phylogenetic analyses based on 16S rRNA gene sequences and coding sequences of 49 core, universal genes defined by Clusters of Orthologous Groups gene families indicated that B1D3A T was a member of the genus Sphingobium. B1D3A T was most closely related to Sphingobium sp. SYK-6, with a 100% 16S rRNA gene sequence similarity. B1D3A T showed 78.1–89.9%average nucleotide identity and 19.5–22.2% digital DNA–DNA hybridization identity with other type strains from the genus Sphingobium. On the basis of phenotypic and genotypic properties and phylogenetic inference, strain B1D3A T should be classified as representing a novel species of the genus Sphingobium, for which the name Sphingobium lignivorans sp. nov. is proposed. The type strain is strain B1D3A T (ATCC TSD-279 T =DSM 111877 T ).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A simple model for jet quenching in Heavy Ion Physics

This is a simple python code to calculate jet quenching for nuclear modification factors in relativistic heavy-ion collisions. This code can be used to generate all results specified in RR0010127, to be submitted to Physical Review C.

Soltz, Ron↗
DOE OSTI2023

SALSA_python (SALSython) v1

SALSA_python (Semi-Analytical Leakage Solutions for Aquifers) is a software that computed semi-analytical solutions for hydraulic head and brine leakage in multilayered aquifer–aquitard systems with geologic pressure forcing. It can simulate brine leakage into aquifers in a multi-aquifer-aquitard system with multiple injection and leaky wells. This situation is encountered in underground CO2 storage wherein brine leakage from pressurized reservoirs into aquifers is of concern. SALSA_python calls the original SALSA[1] subroutines in python by using the salsa2.so library. This enables the incorporation and coupling of SALSA computations into existing python-based codes and tools. SALSA_python is available for Linux and Mac operating systems. [1] Cihan, A., Oldenburg, C. M., & Birkholzer, J. T. (2022). Leakage from coexisting geologic forcing and injection-induced pressurization: A semi-analytical solution for multilayered aquifers with multiple wells. Water Resources Research, 58, e2022WR032343.

Bhuvankar, Pramod↗
DOE OSTI2025

TInCup

SAND2025-11464O TInCuP (Tag Invoked Customization Points) is a modern, header-only C++20 library that addresses the boilerplate problem in tag invoke based customization points. It provides comprehensive code generation and verification tools. 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 OSTI2023

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment. Dataset.

To simulate the effects of mycorrhizal fungi on soil organic matter cycling, we incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model to develop a new soil model Myco-CORPSE. The new model was calibrated and evaluated against soil measurements taken at temperate forests in New Hampshire (NH) and Georgia (GA). A series of scenario analysis were also conducted to explore the conditions under which ectomycorrhizal (ECM) N acquisition processes can induce different soil C accumulation in ECM systems compared to arbuscular (AM) systems.In this data package, we included:-The Python codes of the standard Myco-CORPSE model we developed: "Standard Myco_CORPSE python codes.zip". The main program is the "gradient_sim.py" which calculates the bulk soil microbes and CN content along a user defined gradient of clay, soil temperature, soil moisture and mycorrhizal dominance, and relies on two subprograms "CORPSE_deriv.py" and "CORPSE_integrate.py". "CORPSE_deriv.py" calculated the changes in all simulated soil stock within every time step and "CORPSE_integrate.py" integrate the changes in all simulated soil stock within simulated time period. The program "Plot.py" is used to plot the major outputs produced by the main program "gradient_sim.py".-The modified Python codes of Myco-CORPSE models with site-level environmental inputs (in NH and GA) used to conduct simulations in NH and GA sites: "NH_GA model simulations.zip". -The Python codes used to evaluate the Myco-CORPSE simulation outputs in NH and GA sites against site-level measurements: "Plot NH_GA simulation against measurements.zip". It includes both the evaluation Python code, the model outputs on NH and GA sites, and the measured soil properties in both sites.-The modified Python codes of Myco-CORPSE models "Scenario analysis_model simulations.zip" that is used to conduct scenario analysis of how different litter properties, mycorrhizal fungal traits, climate, and seasonal variation in temperature and vegetation phenology impact the mycorrhizal effects on soil CN properties. The sub file folder "Scenario analysis_litter traits" contains the codes for scenario analysis of different litter properties; The sub file folder "Scenario analysis_ECM types" contains the codes for scenario analysis of different ECM fungal traits; The sub file folder "Scenario analysis_climate&seasonality" contains the codes for scenario analysis of different climate and seasonalities;-"Scenario analysis_model results and plotting codes.zip" contains all the output files from the the scenario analysis of Myco-CORPSE model as described above and the plotting codes used to the generate the heatmaps and scatterplots shown in the manuscript "Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment"The majority of the model outputs did not have specific geographic information or temporal coverage because the analysis we conducted are mainly hypothetical model simulations. We only provided geographic description, coordinates and temporal coverage for those soil measurements which we used for model evaluations (included in the "Plot NH_GA simulation against measurements.zip").

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Data Repository for Multi-Objective Urban Observational Strategies: A risk-based framework for expanding flood sensor networks.

These data support the manuscript "Multi-Objective Urban Observational Strategies: A risk-based framework for expanding flood sensor networks." These data are generated to allow water managers to reason about optimal locations to expand a flood observation system from multiple perspectives, specifically focusing on flood hazards, and population exposure to flooding. The data included are a) a shapefile of individual sensor locations b) a shapefile of river reach catchments, c) raster of FEMA flood likelihood layers d) shapefile of population locations and population socioeconomic characteristics. The code is written in R and includes all files necessary to generate the figures for the associated manuscript. Interactive maps of the final calculated maps of hazard, vulnerability, exposure, and risk are also included as html files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Preliminary description of a new creep-fatigue design method that reduces over conservatism and simplifies the high temperature design process

The report provides the initial description of a new creep-fatigue design method for structural components in high temperature nuclear service. The new method is based on an integrated elastic-perfectly plastic (EPP) analysis and Simplified Model Test (SMT) approach that reduces over conservatism, improves the treatment of elastic follow up, and simplifies the design procedure, when compared with the current creep-fatigue design methods in ASME Boiler and Pressure Vessel Code. Developing the design charts for the EPP-SMT design method requires extrapolating SMT test data as a function of hold time and follow up factor. The report develops the preliminary design charts for Alloy 617 at temperatures between 800°C and 950°C by combining two extrapolation approaches developed in a previous work. The report also presents a comparative analysis between the EPP-SMT design method and the current ASME creep-fatigue design methods by evaluating design life of two sample geometries under different loading conditions. Results from the comparative analysis verify the EPP-SMT design charts but suggest the requirement of additional test data in the low strain range regime for improving the extrapolation procedure that will further reduce the over conservatism in the creep-fatigue damage evaluation. The report also concludes that the EPP-SMT design procedure can account for effect of primary load on creep-fatigue life by using a fixed, bounding value of follow up in constructing the design charts. The conclusions to this report describe the future work required to complete this new design method so it can be codified through a nuclear Code Case.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

The advanced evolution of massive stars

The nuclear rates for reactions involving 12 C and 16 O are key to computing the energy release and nucleosynthesis evolution of massive stars during their advanced burning phases. Ultimately, these burning rates shape the stellar structure and evolution and influence the nature of the compact objects produced at the end of the stellar life. We explore the implications of new nuclear reaction rates from both experimental and theoretical studies for 12 C(α, γ) 16 ​O, 12 C+ 12​ C, 12 C+ 16 ​O, and 16 O+ 16 ​O reactions for massive stars. Our goal is to investigate how the chemical structure and nucleosynthesis evolve from the He-exhaustion stage to the O-burning phase and how these processes influence the ultimate stellar fate. We computed rotating and non-rotating models for stars of different masses at solar metallicity. We used the stellar evolution code GENEC, which includes a large network of nuclear reactions and isotopes involved in advanced phases, as well as updated rates for 12 C(α,γ) 16 O. For the three fusion reactions involving 12 C and 16 O, we considered new rates following a data-driven fusion suppression scenario (hereafter HIN(RES)) and new theoretical rates obtained with time-dependent Hartree-Fock (TDHF) calculations. The updated 12 C(α, γ) 16 ​O rates mainly impact the chemical structure evolution changing the 12 C/ 16 O ratio at He-exhaustion and have little effect on the CO core mass. This variation in the 12 C/ 16 O ratio is in some cases critical for predicting the final fate of the model, which is very sensitive to 12 C abundance, and in particular the 20 M ⊙ remnant may change from a black hole to a neutron star. The He-burning (C-burning) lifetime is also decreased (increased) by about −2% (+15%). The combined new rates for 12 C+ 12 ​C and 16 O+ 16 ​O fusion reactions according to the HIN(RES) model lead to shorter C- and O-burning lifetimes by ≈ − 10%, and −50%, respectively, and shift the ignition conditions to higher temperatures and densities. In contrast, the theoretical TDHF rates primarily affect C-burning, increasing its duration by about 30% and lowering the ignition temperature. These changes modify the chemical structure of the core, the size and duration of C-burning shells, and hence their compactness. They also impact the central and shell nucleosynthesis (by ±1 dex and by factors of ±2–10, respectively), while 12 C+ 16 ​O reaction rates variations remain the least important.

abundances↗
DOE OSTI2026

AstraAI v1

AstraAI is an open-source, structure-aware AI coding agent designed for large scientific and DOE-HPC codebases such as AMReX-based applications. Unlike general-purpose coding assistants, AstraAI combines retrieval-augmented generation (RAG) with compiler-level Abstract Syntax Tree (AST) analysis to perform precise, scope-constrained code modifications. It identifies exact function spans, enforces locality of edits, and maintains cross-file invariants, enabling deterministic and build-safe transformations in complex C++/GPU environments. AstraAI is intended for developers working on large, evolving HPC frameworks where correctness, reproducibility, and structural integrity are critical. Typical use cases include modifying physics kernels, updating GPU device lambdas, and performing multi-file refactors without breaking compilation or runtime semantics. Compared to conventional LLM-based coding agents - even those with repository access - AstraAI provides structural guarantees rather than free-form text patches. It minimizes unintended diffs, prevents scope drift, preserves formatting and build stability, and reduces structural hallucinations. By integrating compiler tooling directly into the generation loop, AstraAI transforms AI-assisted coding from probabilistic text editing into deterministic, structure-preserving program transformation suitable for mission-critical scientific software.

Natarajan, Mahesh [Lawrence Berkeley National Labo↗
DOE OSTI2021

RMACXX: An Efficient High-Level C++ Interface over MPI-3 RMA

Parallel scientific applications can benefit from de- coupling communication and synchronization. One-sided pro- gramming abstractions, which separate communication from syn- chronization, have in fact served as a motivation for partitioned global address space (PGAS) models. However, the use of PGAS models in application codes in a manner that fully exploits the benefit of these programming models requires significant development effort. Meanwhile, a vast majority of scientific codes already use the Message Passing Interface (MPI) and need convenient features to support application-specific one-sided communication scenarios. MPI Remote Memory Access (RMA) can be employed for this purpose. MPI is a low-level API, however, and developing applications with MPI RMA requires programmers to be well versed in its nuances. We present RMACXX, a compact set of C++ bindings to MPI-3 RMA, to ease the use of MPI RMA. Unlike other PGAS models, which may have interoperability issues with MPI, RMACXX is written on top of MPI and uses the same runtime as MPI. The basic functionality of RMACXX adds only a relatively small number of extra instructions (about 20) to the critical communication path. Moreover, RMACXX provides an intuitive API for building a wide variety of scientific applications while enjoying performance matching handwritten MPI-3 RMA codes.

MPI-3 RMA, one-sided communication, PGAS, C++, exp↗
DOE OSTI2023

Community College Energy Code Training Program

The U.S. Department of Energy (DOE) Office of Energy Efficiency and Renewable Energy (EERE) funded a series of research and training programs under the Advanced Building Construction with Energy Efficient Technologies and Practices (ABC) initiative. The aim is to develop new training programs to support the development of a skilled workforce in the building energy efficiency industry. In response, the University of Illinois Smart Energy Design Assistance Center (SEDAC), in collaboration with State Energy Offices in Illinois, Hawaii, and Nevada, developed the Community College Energy Code Training Program. The program is designed to prepare a new generation of energy-literate professionals that understand energy efficiency and its relation to energy codes for new construction and existing buildings through: a) the development of a curriculum and innovative pedagogical materials, b) dissemination through a national network of state energy offices: c) outreach targeting Community College students and instructors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

FLASSH 1.0: Thermal scattering law evaluation and cross section generation for reactor physics applications

The Full Law Analysis Scattering System Hub (FLASSH) is a modern, advanced code which evaluates the thermal scattering law (TSL) along with accompanying cross sections. FLASSH features generalized methods which accommodate any material structure. Historical approximations including the incoherent and cubic approximations have been removed. Instead, the latest release of FLASSH features advanced physics options including distinct corrections (1-phonon contributions) and non-cubic formulations. The non-cubic elastic and inelastic contributions are necessary to accurately evaluate 1-phonon contributions. Both non-cubic and 1-phonon calculations require high-density sampling of the various scattering directions. Optimization and parallelization of these routines were therefore necessary to produce results in a reasonable timeframe. With these notable improvements to the generalized TSL, FLASSH 1.0 meets benchmark requirements, demonstrating noticeable agreement with experiment for both TSLs and the resulting integrated cross sections. Additional features including a graphical user interface (GUI), plotting diagnostics, and formatted output options including ACE files allow users to complete a TSL evaluation with minimal input and maximum flexibility. The user GUI creates input files for FLASSH, reducing user error and also providing built-in error checks. Autofill options and suggested input values help make TSL evaluation accessible to novice users. The FLASSH code is compiled to run on both Windows and Linux platforms with automatic parallelization. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Validation of heat transport modeling using directly driven beryllium spheres

Recent experiments involving directly driven beryllium spheres are reported. Plasma conditions are measured using Thomson scattering with the probe beam pointed 200, 300, and 400 μ m from the surface of the sphere. Laser coupling is assessed using calorimeters that collect scattered light placed at various locations within the target chamber. Laser intensities of 1 0 14 W / c m 2 and 2.5 × 1 0 14 W / c m 2 are chosen to minimize unmodeled laser-plasma interactions (LPIs) that lead to laser-target decoupling. Two-dimensional simulations are compared to the interpreted data using the radiation-hydrodynamics code Lasnex. Heat transport is simulated using flux-limited Spitzer–Harm with both high ( f = 0.15) and low ( f = 0.03) flux limiters and the Schurtz–Nicolai-Busquet (SNB) model. At 1 0 14 W / c m 2 , all three heat transport models agree well with the measurement, demonstrating that the heat flux is local at low intensities near the measurement locations. At 2.5 × 1 0 14 W / c m 2 , the SNB and high flux model roughly match the plasma conditions but predict 2% uncoupled light compared to 10% measured. Additionally, the use of drive multipliers to match the measured coupled light does not alter the agreement between measured and simulated plasma conditions, suggesting that decoupling due to LPI is unlikely to alter this agreement. The low flux model cannot match the plasma conditions and results in 19% scattered light. The use of a resonant absorption model can be used to bring the simulated scattered light into agreement, but the simulated plasma conditions are still in disagreement with the measurement. For this reason, the low flux model is rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗