Search NASA⌕ Search

SEARCH · Search NASA

Results for “O 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 469 records · Page 26

Phase Separation in Ultramassive White Dwarfs

We report ultramassive white dwarfs are extreme endpoints of stellar evolution. Recent findings, such as a missing multi-Gyr cooling delay for a number of ultramassive white dwarfs and a white dwarf with a quasi-Chandrasekhar mass, motivate a better understanding of their evolution. A key process still subject to important uncertainties is the crystallization of their dense cores, which are generally assumed to be constituted of 16 O, 20 Ne, and a mixture of several trace elements (most notably 23 Na and 24 Mg). In this work, we use our recently developed Clapeyron integration technique to compute accurate phase diagrams of three-component mixtures relevant to the modeling of O/Ne ultramassive white dwarfs. We show that, unlike the phase separation of 22 Ne impurities in C/O cores, the phase separation of 23 Na impurities in O/Ne white dwarfs cannot lead to the enrichment of their cores in 23 Na via a distillation process. This severely limits the prospect of transporting large quantities of 23 Na toward the center of the star, as needed in the white dwarf core-collapse mechanism recently proposed by Caiazzo et al. We also show that despite representing ≈10% of the ionic mixture, 23 Na and 24 Mg impurities only have a negligible impact on the O/Ne phase diagram, and the two-component O/Ne phase diagram can be safely used in white dwarf evolution codes. We provide analytic fits to our high-accuracy O/Ne phase diagram for implementation in white dwarf models.

79 ASTRONOMY AND ASTROPHYSICS↗

The SENSEI Generic In Situ Interface: Tool and Processing Portability at Scale [Book Chapter]

One key challenge when doing in situ processing is the investment required to add code to numerical simulations needed to take advantage of in situ processing. Such instrumentation code is often specialized, and tailored to a specific in situ method or infrastructure. Then, if a simulation wants to use other in situ tools, each of which has its own bespoke API [4], then the simulation code team will quickly become overwhelmed with having a different set of instrumentation APIs, one per in situ tool or method. In an ideal situation, such instrumentation need happen only once, and then the instrumentation API provides access to a large diversity of tools. In this way, a data producer’s instrumentation need not be modified if the user desires to take advantage of a different set of in situ tools. The SENSEI generic in situ interface addresses this challenge, which means that SENSEI-instrumented codes enjoy the benefit of being able to use a diversity of tools at scale, tools that include Libsim, Catalyst, Ascent, as well as user-defined methods written in C++ or Python. SENSEI has been shown to scale to greater than 1M-way concurrency on HPC platforms, and provides support for a rich and diverse collection of common scientific data models. Furthermore, this chapter presents the key design challenges that enable tool and processing portability at scale, some performance analysis, and example science applications of the methods.

Bethel, E. Wes↗

Conjugate Heat Transfer Modeling of Salt-Filled Fuel Pins for Stable Salt Reactor Safety Analysis

The Stable Salt Reactor (SSR) combines the proven structural design of light water reactor fuel assemblies with the inherent safety and fuel-cycle advantages of molten salt technology. In its fast reactor configuration, the SSR utilizes recycled nuclear waste as fuel, sealed within narrow salt-filled fuel pins and cooled by a surrounding liquid salt coolant. Reliable transfer of heat from the molten fuel salt through the cladding to the external coolant is essential for both reactor safety and performance. This work investigates conjugate heat transfer (CHT) in the SSR’s salt-filled fuel pins using NekRS, a high-fidelity spectral element computational fluid dynamics (CFD) solver. The analyses capture internal natural convection within the molten fuel salt and external forced convection in the coolant, under steady-state and transient operating conditions. Parametric studies evaluate how variations in reactor power and coolant flow rate influence heat transfer distributions and system response. The high-fidelity CFD results are time-averaged and post-processed for direct comparison with moderate-fidelity Reynolds-averaged Navier–Stokes (RANS) models, and for the development of reduced-order models within the SAM system code. These validated models support fast-running safety analyses of normal and off-normal transients, improving predictive capability for key safety margins. By integrating advanced CFD with system-level safety tools, this study strengthens the modeling framework for SSR design, reduces uncertainty in molten salt CHT simulations, and accelerates the engineering and licensing of next-generation nuclear reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Performance-portable Binary Neutron Star Mergers with AthenaK

We introduce an extension to the AthenaK code for general-relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes using a 3+1 conservative Eulerian formulation. Like the fixed-spacetime GRMHD solver, we use standard finite-volume methods to evolve the fluid and a constrained-transport scheme to preserve the divergence-free constraint for the magnetic field. We also utilize a first-order flux correction (FOFC) scheme to reduce the need for an artificial atmosphere and optionally enforce a maximum principle to improve robustness. We demonstrate the accuracy of AthenaK using a set of standard tests in flat and curved spacetimes. Using a SANE accretion disk around a Kerr black hole, we compare the new solver to the existing solver for stationary spacetimes using the so-called "HARM-like" formulation. We find that both formulations converge to similar results. We also include the first published binary neutron star (BNS) mergers performed on graphical processing units (GPUs). Thanks to the FOFC scheme, our BNS mergers maintain a relative error of $\mathcal{O}$(10 –11 ) or better in baryon mass conservation up to collapse. Finally, we perform scaling tests of AthenaK on OLCF Frontier, where we show excellent weak scaling of ≥80% efficiency up to 32,768 GPUs and 74% up to 65,536 GPUs for a GRMHD problem in dynamical spacetimes with six levels of mesh refinement. AthenaK achieves an order-of-magnitude speedup using GPUs compared to CPUs, demonstrating that it is suitable for performing numerical relativity problems on modern exascale resources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Connect microstructure evolution to property degradation with validated simulation

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

Implications of a two-fluid heat diffusion model for RF condensation effects in magnetic islands

Theoretical calculations have predicted an “RF condensation” effect through which the temperature perturbation and current drive may be nonlinearly enhanced in electron cyclotron (EC) heated magnetic islands. Previously, two-fluid effects had been neglected in the OCCAMI code used to model RF condensation. For realistic experimental conditions, we show that two-fluid effects in an EC heated plasma may have a major impact on RF condensation. Less collisional plasmas, where electron and ion temperatures are decoupled, exhibit an enhanced linear temperature perturbation in the island, and an even stronger nonlinear enhancement of the perturbed temperature. This can occur even when nonlinear effects are small in the single-fluid limit. We also calculate corresponding hysteresis effects. Finally, we demonstrate that other physical effects, including the rotation of the magnetic island and the Gaussian spread of the beam from a gyrotron, can significantly decrease the predicted observable RF condensation effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Prospects for H-mode inhibition in negative triangularity tokamak reactor plasmas

Instability to high toroidal mode number ($\textit{n}$) ballooning modes has been proposed as the primary gradient-limiting mechanism for tokamak equilibria with negative triangularity (NT) (δ) shaping, preventing access to strong H-mode regimes when δ is sufficiently negative. To understand how this mechanism extrapolates to reactor conditions, we model the infinite-$\textit{n}$ ballooning stability as a function of internal profiles and equilibrium shape using a combination of the CHEASE and BALOO codes. While the critical δ required for avoiding 2nd stability to high-$\textit{n}$ modes is observed to depend in a complicated way on various shaping parameters, including the equilibrium inverse aspect ratio, elongation and squareness, equilibria with NT are robustly prohibited from accessing the 2nd stability region, offering the prediction that NT reactors should maintain L-mode-like operation. In order to access high-$\textit{n}$ 2nd stability, the local shear over the entire bad curvature region must be sufficiently negative to overcome curvature destabilization on the low field side. Finally, scalings of the ballooning-limited pedestal height are provided as a function of plasma parameters to aid future scenario design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurement and modeling of proton-induced reactions on arsenic from 35 to 200 MeV

72 As is a promising positron emitter for diagnostic imaging that can be employed locally using a 72 Se generator. However, current reaction pathways to 72 Se have insufficient nuclear data for efficient production using regional 100–200 MeV high-intensity proton accelerators. In order to address this deficiency, stacked-target irradiations were performed at LBNL, LANL, and BNL to measure the production of the 72 Se/ 72 As positron emission tomography (PET) generator system via 75 As (p, x) between 35 and 200 MeV. This work provides the most well-characterized excitation function for 75 As ( p , 4 n ) 72 Se starting from threshold. Additional focus was given to report the first measurements of 75 As (p, x) 68 Ge and bolster an already robust production capability for the highly valuable 68 Ge/ 68 Ga PET generator. Thick target yield comparisons with prior established formation routes to both generators are made. In total, high-energy proton-induced cross sections are reported for 55 measured residual products from 75 As, nat Cu , and nat Ti targets, where the latter two materials were present as monitor foils. These results were compared with literature data as well as the default theoretical calculations of the nuclear model codes TALYS, COH, EMPIRE, and ALICE. Reaction modeling at these energies is typically unsatisfactory due to few prior published data and many interacting physics models. Therefore, a detailed assessment of the talys code was performed with simultaneous parameter adjustments applied according to a standardized procedure. Particular attention was paid to the formulation of the two-component exciton model in the transition between the compound and preequilibrium regions, with a linked investigation of level density models for nuclei off of stability and their impact on modeling predictive power. This paper merges experimental work and evaluation techniques for high-energy charged-particle isotope production in an extension to an earlier study of this kind.

43 PARTICLE ACCELERATORS↗

Enhancements supporting IC usage of PEM libraries on next-gen platforms

This milestone reports on the culmination of several years of effort by multiple PEM support software development teams to provide capabilities for use in LLNL-developed integrated codes on next-gen ASC platforms, including GPU support. We will provide a survey of relevant Application Program Interfaces (API) that are required to support LLNL IC code capability on relevant architectures, with a focus on Sierra and El Capitan. We will identify and summarize all dependencies between PEM supported libraries and IC supported physics codes. We will provide an assessment of algorithmic improvements that have been deployed, as well as future developments that are required to complete the GPU porting efforts. This assessment will include a description of programming models adopted by each of the PEM projects, distinct algorithmic challenges for each of the capabilities, and information about sharing GPU memory between the APIs and host codes. We will develop targeted test problems to assess computational performance. Finally, this milestone will result in identification of gaps in our effort to assist the LLNL ASC program in prioritization of effort for porting software to El Capitan.

97 MATHEMATICS AND COMPUTING↗

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Simulation of full wave, without paraxial approximation, high-resolution solution of wave equations in frequency domain in the electron cyclotron resonance (ECR) frequency range for realistic Tokamak plasma parameters became possible by using recently formulated hybrid iterative algorithm [Svidzinski et al., Phys. Plasmas 25, 082509 (2018)] for numerically solving discretized wave equations. This approach combines time evolution and iterative relaxation techniques into iteration cycles. This algorithm is implemented in 2D code FullWave, solving wave equations in Tokamaks in cold and hot plasma models, and it has been tested in 3D full wave iterative RF beams simulation tool, which is presently being developed to model 3D ECRH RF beams in fusion devices using dynamic grid adaptation. The results of 2D full wave modeling, assuming specified toroidal mode number, of ECRH RF beams in DIII-D plasma, performed in the cold and hot plasma models for outboard and top launch scenarios using FullWave are presented. Nonlocal hot plasma response model, based on accurate numerical solution of linearized Vlasov equation, is used to model beam propagation and absorption in the 2nd electron cyclotron harmonic region. Demonstration of capability of the hybrid iterative algorithm to model ECRH RF beams in 3D is made by simulating a substantial part of realistic beam in DIII-D, launched from outboard side of the machine. All relevant physics of RF beam propagation, most of which is not captured in paraxial approximation, such as beam's divergence, interference between the X and O modes in the beam, X-O mode conversion, beam splitting into the X and O mode beams, transformation of beam's cross section, and absorption at the 2nd electron cyclotron harmonic, is captured in the simulations. A numerical technique to find an optimal beam polarization at the launcher to launch a nearly pure X or O mode beam in plasma is developed and tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Monarch Initiative in 2024: an analytic platform integrating phenotypes, genes and diseases across species

Abstract Bridging the gap between genetic variations, environmental determinants, and phenotypic outcomes is critical for supporting clinical diagnosis and understanding mechanisms of diseases. It requires integrating open data at a global scale. The Monarch Initiative advances these goals by developing open ontologies, semantic data models, and knowledge graphs for translational research. The Monarch App is an integrated platform combining data about genes, phenotypes, and diseases across species. Monarch's APIs enable access to carefully curated datasets and advanced analysis tools that support the understanding and diagnosis of disease for diverse applications such as variant prioritization, deep phenotyping, and patient profile-matching. We have migrated our system into a scalable, cloud-based infrastructure; simplified Monarch's data ingestion and knowledge graph integration systems; enhanced data mapping and integration standards; and developed a new user interface with novel search and graph navigation features. Furthermore, we advanced Monarch's analytic tools by developing a customized plugin for OpenAI’s ChatGPT to increase the reliability of its responses about phenotypic data, allowing us to interrogate the knowledge in the Monarch graph using state-of-the-art Large Language Models. The resources of the Monarch Initiative can be found at monarchinitiative.org and its corresponding code repository at github.com/monarch-initiative/monarch-app.

60 APPLIED LIFE SCIENCES↗

Highly Extensible X-ray Diffraction Toolkit

This code has been previous released as open source under LLNL-CODE-529294. The HEXRD program comprises a library of X-ray diffraction analysis tools and a stand-alone GUI application. While there are many software packages that provide similar capabilities, HEXRD was designed to provide an extremely flexible and extensible basis for describing a broad array of X-ray diffraction instruments; this includes both mono- and poly-chromatic modalities as well as poly- and single-crystal samples. The high-level objects include Material: an abstraction for defining crystalline materials, including unit cell parameters and anisotropic elastic moduli Imageseries: an interface for series of diffraction images containing one or more frames, including provisions for metadata. Instrument: an abstraction of an X-ray diffraction instrument including support for the definition of multiple detector elements as well as powder, Laue, and mono-chromatic single/multi-grain diffraction analysis modalities (including attendant calibration routines). The hexrdgui package is a standalone PySide2/Qt5 GUI built on top of hexrd. The 0.8.1 release includes various bug fixes as well as some new functionality in the GUI, including: o Enhancements to the import tools for LLNL diffraction instruemtns (PXRDIP and TARDIS) o The ability to save 1-d LeBail fits of integrated diffraction spectra.

Barton, NathanR.↗

Research Plan and Preliminary Results in Developing the Fabrication Parameters for Alloy 709 in Different Product Forms―Grain Coarsening Temperature Evaluation

The Advanced Reactor Technologies (ART) Program has established a multi-year plan to develop Alloy 709 advanced stainless steel (A709), generate the data package and develop material-specific design parameters in qualifying it as a new structural material for Class A component design in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code, Section III, Division 5, High Temperature Reactors. In collaboration with material vendors, the Advanced Materials Development activities under ART have successfully scaled the A709 plate form production from a laboratory heat of 500 pounds to commercial heats totaling 133,000 pounds of A709 plate fabricated from three heats. The goal of the overall A709 development program is to establish the necessary microstructural and mechanical properties relationship for A709 to ultimately develop fabrication parameters for other product forms such as bars, piping and forging using the available ART A709 materials. The objective of this A709 development work at ORNL in FY2023 is to experimentally determine grain coarsening behavior for the A709 heats and to experimentally generate the continuous cooling precipitation (CCP) diagram for A709 using the as-rolled commercial heat plate materials. Integral to this work is the characterization of the as-rolled materials and the determination of an effective solution annealing process. This report summarizes the work performed to identify the grain coarsening temperature for commercial heat 58776-3RB fabricated by G. O. Carlson and heat 529900-02 fabricated by Allegheny Technologies Incorporated (ATI) Flat Rolled Products.

36 MATERIALS SCIENCE↗

Gamma-Ray Imaging to Monitor Accumulating Holdup

Over time, nuclear materials can accumulate in process equipment and keeping track of such material is of interest to many stakeholder communities including safeguards, process monitoring, criticality safety, and radiological safety. Passive gamma-ray imaging provides a visual and quantitative tool that can be applied to this problem, allowing one to determine the amount of material present from the images. In this work we emulate the accumulation of highly enriched uranium (HEU) in a large rectangular stainless-steel duct by stacking 5, 27 × 50 cm2 cards, each of which has approximately 11 g of 235U in oxide form distributed across its surface, behind each other. A series of images were collected using a coded-aperture gamma-ray imager, starting with a single card, adding a second card, etc. Long integrations were collected to allow evaluating the performance of the approach as a function of integration time. The images are analyzed using first principles to determine the amount of material in each image. The results are compared to the known amount of material in each configuration, and as a function of time. The analysis includes a first-principles error propagation and to further evaluate the uncertainty, a boot-strap sampling approach is also applied to the data. The work offers a proof of concept for unattended monitoring of uranium holdup with a quantitative gamma-ray imaging technique.

Ziock, Klaus-Peter↗

PV System Owner's Guide to Identifying, Assessing, and Addressing Weather Vulnerabilities, Risks, and Impacts

The report is designed to help federal facility managers, engineers, contractors, and consultants identify the most common PV system vulnerabilities that could lead from mild to catastrophic loss of PV systems during severe weather events. The guide helps users understand the type and severity of severe weather that can occur at any given location, and provides pre- and post-storm operations and maintenance (O&M) measures that can reduce the potential for damage to a PV system during a severe weather event and help speed the recovery of a PV system damaged during a storm. The guide emphasizes safety in responding to damage from severe weather events.

14 SOLAR ENERGY↗

Nationwide HVAC Energy-Saving Potential Quantification for Office Buildings with Occupant-Centric Controls in Various Climates

The occupant-centric control (OCC) is receiving an increasing attention due to its ability to reduce building heating ventilation and air-conditioning (HVAC) system energy consumptions while not affecting the occupant thermal comfort. This paper aims to investigate and quantify the nationwide energy-saving potential of implementing the occupant-centric HVAC controls in typical office buildings using a whole building simulation software EnergyPlus. First, the medium office and large office from the Department of Energy (DOE) Commercial Prototype Building Models (CBPM) were enhanced to have detailed layouts and dynamic occupancy schedules. Then, a comprehensive simulation plan was created by incorporating the multiple zone-level and system-level occupant-centric building HVAC controls recommended by the updated ASHRAE Standard 90.1 – 2019 and ASHRAE Guideline 36 – 2018. Three control scenarios with different occupancy sensing methods were identified in this simulation plan. A nation-wide parametric analysis which includes two building types, three occupancy sensing scenarios, two building code versions, and 16 U.S. climate zones was carried out. The simulation results of the key control variables and HVAC energy consumption suggest that generally, both the occupancy presence sensor and occupant counting sensor could achieve energy savings for the office buildings in majority of the scenarios. However, compared with the occupancy presence sensor, which could support both the temperature setpoint reset and operational breathing zone airflow rate reset for the unoccupied zones, the occupant counting sensor only brings a marginal benefit. Besides, a higher HVAC energy-saving ratio could be achieved in the heating-dominated zone, since the energy reduction brought with the minimum outdoor airflow rate reset is stronger in the heating mode.

Pang, Zhihong↗

Data for Stetten et al. (2025), "Biogeochemical controls on iron speciation and cycling across upland to shoreline gradients in freshwater and estuarine coastal soils (Lake Erie and Chesapeake Bay, United States)"

Coastal environments are dynamic interfaces that mediate carbon and nutrient exchanges between terrestrial landscapes and open waters, but it is unclear how biogeochemical reactions, in particular iron (Fe) redox transformations, affect the understanding and prediction of coastal ecosystem functions. This dataset includes measurements from two freshwater sites in the Western and Central basins of Lake Erie (Ohio, United States) and two estuarine sites in the Chesapeake Bay (Maryland, United States); the analytical results were reported by Stetten et al. (2025) in Science of the Total Environment. It was produced as part of the COMPASS-FME project, which seeks to advance a scalable, predictive understanding of the fundamental biogeochemical processes, ecological structure, and ecosystem dynamics that distinguish coastal terrestrial-aquatic interfaces from the purely terrestrial or aquatic systems to which they are coupled. The sites were sampled in November 2022 (CRC), December 2022 (MSM), February 2023 (GCW), and March 2023 (OWC); site codes follow those used by Pennington et al. (2025).The dataset consists of the following soil data:- Solid data (Fe concentration, etc.)- Porewater data (sulfate, sulfide, etc.)- Linear combination fitting results of X-ray absorption near edge structure (XANES) spectra; i.e., quantitative results of the oxidation state of Fe, indicated as a proportion of pure Fe(III) and Fe(II) model compounds- Linear combination fitting results of EXAFS (extended X-ray absorption fine structure) spectra, indicated as proportion of of Fe-model compounds (illite, smectite, etc.)Each data type has a single file in comma-separated value (CSV) format. No special software is required to read it.

54 ENVIRONMENTAL SCIENCES↗