SEARCH · Search NASA

Results for “Open source”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Emulation of radiation transport in 3D stochastic media using 1D planar Monte Carlo stochastic media radiation transport algorithms

A subset of stochastic media radiation transport problems involves those in which radiation is incident on a thin slab of stochastic material. Particle tracking in 3D for such problems is expensive, and 1D planar models lack accuracy because they only allow the material to change in one dimension. Therefore, we propose dimensional emulation, which through a slight modification allows existing 1D planar geometry stochastic media radiation transport models to reproduce results from the equivalent 3D models by allowing the material to change in all three dimensions, reproducing the fidelity of the 3D model for the low computational cost of the 1D planar model. In this work, we apply dimensional emulation to three Monte Carlo stochastic media radiation transport models: Chord Length Sampling (CLS), the Local Realization Preserving method (LRP), and a variant of Conditional Point Sampling (CoPS). For a common Markovian benchmark set, the 3D emulation variants of these algorithms are numerically verified to reproduce the results of the 3D variants within statistics while running 1.3 to 2 times faster in the implementation within Sandia National Laboratories open-source research code PlaybookMC. The 3D emulation variants are also shown to yield a 72%–92% reduction in error for the thin slab problems in comparison to the 1D benchmark. As a result, the 3D emulation variant of CLS and CoPS-1 are shown to reproduce 3D CLS results that were used to approximate results for a 3D spherical inclusion geometry benchmark set.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Direct sensitivity analysis on the parameterization of crystal plasticity models

Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.

Elementary Effects Test↗

DOE OSTI2025

GMFOLD: Subgraph matching for high-throughput DNA-aptamer secondary structure classification and machine learning interpretability

Aptamers are oligonucleotide receptors that bind to their targets with high affinity. Here, we consider aptamers comprised of single-stranded DNA that undergo target-binding-induced conformational changes, giving rise to unique secondary and tertiary structures. Given a specific aptamer primary sequence, there are well-established computational tools (notably mfold) to predict the secondary structure via free energy minimization algorithms. While mfold generates secondary structures for individual sequences, there is a need for a high-throughput process whereby thousands of DNA structures can be predicted in real-time for use in an interactive setting, when combined with aptamer selections that generate candidate pools that are too large to be experimentally interrogated. We developed a new Python code for high-throughput aptamer secondary structure determination (GMfold). GMfold uses subgraph matching methods to group aptamer candidates by secondary structure similarities. We also improve an open-source code, SeqFold, to incorporate subgraph matching concepts. We represent each secondary structure as a lowest-energy bipartite subgraph matching of the DNA graph to itself. These new tools enable thousands of DNA sequences to be compared based on their secondary structures, using machine-learning algorithms. This process is advantageous when analyzing sequences that arise from aptamer selections via systematic evolution of ligands by exponential enrichment (SELEX). This work is a building block for future machine-learning-informed DNA-aptamer selection processes to identify aptamers with improved target affinity and selectivity and advance aptamer biosensors and therapeutics.

Aptamer↗

DOE OSTI2025

Thermodynamic properties of ZrCl 4 with LiCl, NaCl, KCl, CsCl, MgCl 2 , and UCl 3 for molten salt reactor applications

Here, a set of self-consistent Gibbs energy functions has been developed to describe the thermochemical behavior of the major reactive fission product ZrCl 4 with key components in chloride molten salt reactors (MSRs): LiCl–ZrCl 4 , NaCl–ZrCl 4 , KCl–ZrCl 4 , CsCl–ZrCl 4 , MgCl 2 –ZrCl 4 and UCl 3 –ZrCl 4 . Low ZrCl 4 concentration phase equilibria most relevant to MSR applications have been confirmed via differential scanning calorimetry for NaCl–ZrCl 4 , and X-ray diffraction analysis performed on equilibrated samples of NaCl–ZrCl 4 , KCl–ZrCl 4 , MgCl 2 –ZrCl 4 , and UCl 3 –ZrCl 4 . Within the framework of the modified quasi-chemical model in the quadruplet approximation, extrapolations of pseudo-binary models were generated to represent KCl–MgCl 2 –ZrCl 4 , KCl–NaCl–ZrCl 4 , and MgCl 2 –NaCl–ZrCl 4 , which show agreement with phase equilibria data where available without the use of ternary interaction parameters. The optimized thermodynamic descriptions for these systems and others are available in the open-source compendium Molten Salt Thermal Properties Database – Thermochemical (MSTDB–TC).

CALPHAD↗

DOE OSTI2025

PHASE: Personalized Head-based Automatic Simulation for Electromagnetic properties in 7T MRI

Accurate and individualized human head models are becoming increasingly important for electromagnetic (EM) simulations. These simulations depend on precise anatomical representations to realistically model electric and magnetic field distributions, particularly when evaluating Specific Absorption Rate (SAR) within safety guidelines. State of the art simulations use the Virtual Population due to limited public resources and the impracticality of manually annotating patient data at scale. Here, this paper introduces Personalized Head-based Automatic Simulation for EM properties (PHASE), an automated open-source toolbox that generates high-resolution, patient-specific head models for EM simulations using paired T1-weighted (T1w) magnetic resonance imaging (MRI) and computed tomography (CT) scans with 14 tissue labels. To evaluate the performance of PHASE models, we conduct semi-automated segmentation and EM simulations on 15 real human patients, serving as the gold standard reference. The PHASE model achieved comparable global SAR and localized SAR averaged over 10 grams of tissue (SAR-10g), demonstrating its potential as a promising tool for generating large-scale human model datasets in the future. The code and models of PHASE toolbox have been made publicly available: https://github.com/hrlblab/PHASE.

Deep learning↗

DOE OSTI2026

Evaluation of simulated HPGe detector efficiencies in OpenMC compared to MCNP

In late 2023, the open-source radiation transport code OpenMC introduced a pulse-height tally (PHT) feature, enabling users to track the total energy deposited by individual photons in OpenMC cells. This function represents an important improvement of the OpenMC code because it allows users to simulate the response function of a gamma detector without relying on closed-source alternatives. Despite this, limited work has been published evaluating OpenMC as a radiation transport code for simulating gamma spectroscopy experiments. This study attempts to demonstrate the usefulness of OpenMC in this space by directly comparing its PHT output to MCNP (a trusted industry-standard Monte Carlo code). In the first half of this study, a radiation transport experiment is described in which a detector with a complex internal geometry is exposed to various gamma-emitting isotopes ( 133 Ba, 137 Cs, and 60 Co) over a range of distances. The setup of this experiment was modeled in MCNP with sufficient detail to capture the efficiency characteristics of a high purity germanium (HPGe) detector for the primary gammas of each isotope. After good agreement between the MCNP model and experiment was achieved, an identical model was produced in OpenMC to allow direct comparison between the full energy peak (FEP) values produced in OpenMC and MCNP. The results show strong agreement between OpenMC and MCNP across the full range of tested energies, with each model’s FEP values typically within 2% of each other and most FEP areas within 4% of experimental data. Worse agreement was seen between the Monte Carlo codes and experiment below 300 keV (an expected result). For the 662 keV line of 137 Cs, both codes were found in poor agreement with experiment and each other over the full range of distances tested (possibly indicating an error with the 137 Cs experimental data). Ignoring the anomalous results of the 662 keV line, all other data show good qualitative and quantitative agreement between MCNP and OpenMC. This result demonstrates the accuracy of OpenMC’s PHT feature for spectroscopic applications in which detector efficiency is a primary concern.

07 - ISOTOPES AND RADIATION SOURCES↗

DOE OSTI2024

Nuclear safety Enhanced: A Deep dive into current and future RAVEN applications

As the horizon of nuclear energy expands with the advent of small modular reactors, IV generation reactors, and fusion reactors, there is a growing perspective that the licensing process could benefit from a more comprehensive approach. Moving beyond traditional deterministic and PRA analysis might pave the way for a novel safety analysis paradigm propelled by the increasing computational power at our disposal. This paper explores different methodologies that can improve the outcomes of nuclear safety analysis. These range from uncertainty quantification techniques, aimed at enhancing the precision of safety margins, to deploying dynamic event trees by driving system code simulations, capturing the potential evolutions of severe accidents. These methodologies introduce innovative dimensions to safety analysis, considering the consequences of postulated events and the dynamics of accident sequences. However, they also bring forth challenges, especially in managing the complexity and sheer volume of potential scenarios. The paper touches upon some strategies to counter these challenges, emphasizing the importance of adaptability and continuous evolution in the face of emerging nuclear safety concerns. Additionally, the paper sheds light on the need for advanced tools to apply these methodologies. Among these tools is RAVEN, an open-source software designed for parametric and probabilistic analyses. Its core components, including distribution, sampler, and reduced order model, enable various applications, from risk assessment and mitigation to dynamic learning and plant control logic simulations.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Alchemy: A Model-Based Approach for 2D to 3D Autonomous Nuclear System Design

Engineering design of nuclear power plant (NPP) piping and equipment systems frequently bypasses crucial 2D system planning, instead moving straight to 3D modeling. This often leads to designs that exceed building envelope constraints, forcing expensive and time-consuming redesigns. When 2D modeling is employed, it typically involves labor-intensive manual workflows that convert 2D drawings into 3D models, resulting in inefficiencies and errors across design iterations. These workflows further suffer from poor software interoperability and dependence on proprietary software ecosystems, thereby contributing to schedule delays and cost overruns. This paper presents Alchemy, an autonomous framework that transforms 2D system definitions into Industry Foundation Classes (IFC)-compliant 3D building information models (BIMs) for expediting nuclear facility design at the conceptual preliminary phase. Using a model-based approach, the framework treats the 2D system diagram as the central reference model employed to automatically generate all subsequent outputs, ensuring consistency between the system definition and the resulting physical design. A web-based interface enables engineers to define hierarchical system topologies including associated equipment, geometric properties, and connectivity requirements. A two-phase equipment layout optimization algorithm automatically computes collision-free spatial configurations within predefined building envelopes. An artificial intelligence (AI)-assisted pipe routing module then generates orthogonal, collision-free routing paths, allowing the user to select either an A* search-based method or an Ant Colony Optimization (ACO)-based method. All outputs are authored natively in IFC format, relying on open-source technologies and standardized formats in order to ensure extensibility and eliminate proprietary software dependencies. The proposed framework is validated on two representative pressurized-water reactor (PWR)-based case studies, for which it autonomously generates IFC-compliant 3D models in minutes, drastically reducing workflows that typically require hours of manual effort. The generated model demonstrates topologically correct equipment placement, physically plausible spatial relationships, and collision-free pipe routing consistent with known PWR loop configurations. This work represents a foundational step toward digital engineering for nuclear facility preliminary design, with future ongoing development targeting design code compliance and expanded system complexity.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Slender-body approach for computing second-order wave loads in the frequency domain

This work presents a slender-body approach to evaluate the second-order wave loads acting on a floating structure in the frequency domain. The approach is in the same spirit as the common use of Morison’s equation to approximate the wave loads without solving the radiation/diffraction problem. To do so, we employ Rainey’s equation, which can be seen as an extension of the inertial part of Morison’s equation to include nonlinear effects. We introduce modifications to Rainey’s formulation in order to evaluate wave kinematics at the mean body position instead of the original approach of considering instantaneous displacements. We also propose a simple approximation to partially account for wave scattering effects on the second-order loads based on the analytical solution of a surface-piercing bottom-mounted vertical circular cylinder. Though limited to structures composed of cylinders, this slender-body approach is orders of magnitude faster than computing second-order wave coefficients with a radiation/diffraction code. We implemented this approach for difference-frequency (slow drift) loads in an open-source frequency-domain floating wind turbine model. We present comparisons against results obtained with radiation/diffraction theory for three reference floating wind turbine designs: the OC3-Hywind spar, the OC4-DeepCwind semisubmersible, and the VolturnUS-S semisubmersible. In general, the results show that the proposed slender-body approach with the correction to approximate wave scattering effects provides useful estimations of the difference-frequency wave loads and the resulting motions of the floater.

17 WIND ENERGY↗

DOE OSTI2025

Stabilizing dynamic subsea power cables using Bi-stable nonlinear energy sinks

This study investigates vibration mitigation of dynamic subsea cables through passive bi-stable nonlinear energy sinks (B-NESs). These devices suppress vibration energy in a broadband way, and can be regarded as extensions of classical linear tuned mass dampers (TMDs) which are narrowband devices. Through the open-source MoorDyn library, we simulated the vibrations of a vertical subsea cable equipped with a set of B-NESs. Multi-objective optimization was performed to detect the B-NES parameters for optimal mitigation of the cable vibrations. Advanced signal processing verified the efficacy of the optimized B-NESs not only to rapidly absorb and locally dissipate vibration energy, but also to nonlinearly scatter vibration energy from low to high frequencies within the cable itself. This last feature is especially beneficial for vibration mitigation of the undersea cable, as at higher frequencies the cable vibrations exhibit drastically reduced amplitudes and are more effectively dissipated by inherent structural damping and hydrodynamic radiation damping. This contrasts with traditional TMDs which can mitigate vibration only at a single frequency. Furthermore, our robustness study confirms the B-NES's effectiveness under even varying environmental conditions. Overall, the B-NES's capacity for broadband vibration mitigation renders it a promising retrofit solution for improving the performance and operational safety of dynamic power cables in offshore wind farms and other marine applications.

17 WIND ENERGY↗

DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

3D Poisson equation↗

DOE OSTI2024

Asynchronous GPU-based DEM solver embedded in commercial CFD software with polyhedral mesh support

A novel graphical processing unit-based discrete element method solver is introduced to improve stability, performance, and provide seamless integration into commercial or open-source computational fluid dynamics software. A key innovation is eliminating a need for network communication between solvers, which was previously required for cross-platform coupling. This is accomplished by a direct coupling method that employs dynamic-linked libraries. Furthermore, the solver optimizes memory usage by streamlining the particle-cell search algorithm by eliminating the cells' searching grid. This ensures the solver is compatible with a wide range of mesh types, providing high geometric flexibility. The approach simplifies the simulation process by directly incorporating computational fluid dynamics mesh information into the discrete element method solver. The performance analysis indicates about sixteen times boost in computational speed compared to benchmark central processing unit-based solvers. Finally, the solver's compatibility with polyhedral meshes, a vital advantage for complex geometries, is tested against a referenced study regarding the simulation of an immersed-tube fluidized bed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

On the optimal sizing of power take-off systems for wave energy converters

The power take-off (PTO) system in a wave energy converter (WEC) is the means by which the energy in an ocean wave is converted into useful energy. There has been a general lack of technology convergence in the wave energy industry and PTO design has been largely device-specific with the design process determined by the WEC developer. An optimized PTO design can significantly affect the efficiency, reliability, performance, and overall cost and viability of a device. There remains a need in the marine energy industry for best practices and efficient design processes to be identified as it could greatly benefit developers and researchers in designing an optimal WEC PTO. This research utilizes the open-source Matlab-Simulink based software, WEC-Sim, which was developed by the National Renewable Energy Laboratory, and the Eagle, a high-performance computing system, to model and optimize two WECs of different archetypes and then investigates if any trends exist that can be exploited for greater efficiency in optimal PTO design. The results from this research indicate that there may be an optimal sizing for WEC PTOs for both power and PTO force rating and that substantial reductions can be made in the cost of a device without significant loss in the energy produced. Furthermore this optimal sizing may be independent of WEC type and deployment location. With an optimal, or near optimal, approach to WEC PTO sizing demonstrated, a methodology is proposed to address the challenge of nameplate ratings in the wave energy industry.

13 HYDRO ENERGY↗

DOE OSTI2025

DFAT: A web-based toolkit for estimating demand flexibility in building-to-grid integration

Demand Flexibility Assessment Tool (DFAT) is an open source web-based tool that estimates the demand flexibility potential of common control strategies in commercial buildings. The toolkit features a demand flexibility estimation tool that contains two calculators, basic and advanced, based on the level of input of customer data. The basic version calculates demand shed metrics for the control strategy “global temperature adjustment” and “cycle on/off compressors” using customer building information, local weather data, and electrical meter data. The advanced version, which uses detailed HVAC equipment data, calculates demand flexibility metrics for control strategies such as static pressure reset, global temperature adjustment, and cycle on/off compressors. In addition to the demand flexibility estimation tool, this toolkit offers a benchmarking tool that helps facility operators, aggregators, and utility resource managers assess demand flexibility opportunities, quantify/verify performance, and compare their performance against that of their peers.

Leong, Michael↗

DOE OSTI2025

NEML2: An efficient and modular multiphysics constitutive modeling library for hybrid computing environments

This paper presents NEML2, an open-source, high-performance library developed for constitutive material modeling, designed to support the flexible and modular development of models for complex material behavior. Building on the foundational structure of its predecessor, NEML, the NEML2 library introduces significant improvements, including enhanced vectorization, automatic differentiation, and seamless integration with PyTorch, facilitating the application of machine learning techniques in material simulations. NEML2 provides a C++ backend with Python bindings, enabling users to create custom material models that can be executed efficiently on both CPU and GPU platforms. The library also supports coupling with Multiphysics simulation frameworks like MOOSE, making it suitable for realistic simulations involving coupled physical processes. Rigorous quality assurance through unit and regression testing ensures the reliability of results, while the extensible, user-friendly design encourages collaboration and reproducibility across the scientific community. This paper provides an overview of NEML2’s architecture, core features, and applications, highlighting its impact on accelerating material qualification and advancing computational methods in materials science.

GPU↗

DOE OSTI2025

STREAM: A technology planning and capacity expansion model for the industrial sector

The Strategic Technology Roadmapping and Energy, Environmental, and Economic Analysis Model—STREAM—is an optimization-based modeling tool and analysis framework to assist with strategic planning and technology investments of the industrial sector. This open-source framework is written in Julia using the JuMP package, which enables users to model future “pathways” for incumbent and future production technologies, costs, fuels and energy carriers, and energy and non-energy environmental impacts from industries as they transform in pursuit of a robust and competitive manufacturing sector. The model starts with an initial stock of industrial production technologies and assets at a facility level and then determines pathways that minimize cost, subject to an array of possible constraints on demand, market shares, environmental flows, and other exogenously specified operational considerations such as capacity utilization rates or regional energy costs. Key features of the framework include flexibility to model a wide range of industries and industrial technologies/processes at varying levels of granularity, ability to perform parametric sensitivity analyses, and ability to visualize model results using visualization objects.

capacity expansion↗

DOE OSTI2025

EnergyPlus-MCP: A model-context-protocol server for ai-driven building energy modeling

Traditional building energy modeling with the EnergyPlus building performance simulation engine requires domain expertise, programming skills, and intensive manual efforts limiting its effective adoption. This paper introduces EnergyPlus-MCP, the first open-source Model Context Protocol (MCP) server specifically designed for EnergyPlus simulation workflows, establishing a new foundational infrastructure for AI-driven building energy modeling. The MCP server implements a layered architecture with 35 specialized tools spanning model management, editing and analysis, HVAC and other systems configuration inspection, and simulation execution, enabling Large Language Models to interact with EnergyPlus through conversational interfaces. The server addresses critical workflow barriers by automating model validation, streamlining energy efficiency measures modification, and providing intelligent output management with interactive visualization. Through practical demonstrations using a multi-zone building retrofit analysis, we show how the EnergyPlus-MCP server significantly reduces manual efforts while maintaining full simulation rigor. By providing accessible natural language interfaces to sophisticated building energy analysis, this approach enables scalable deployment of simulation expertise across public and private organizations, educational institutions, and research teams, fundamentally transforming traditional building energy modeling practices.

AI↗