Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation framework”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

NEAMS-IPL MOOSE Framework Activities 2017

The MOOSE Framework is a modular pluggable framework for building complex simulations. The ability to add new objects with custom syntax is a core capability that makes MOOSE a powerful platform for coupling multiple applications together within a single environment. The creation of a new, more standardized JSON syntax output improves the external interfaces for generating graphical components or for validating input file syntax automatically without the need for explicit developer interaction. The design of this interface and the requirements it satisfies are covered in this short report.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A New Capability of E4D For 3D Parallel Joint Inversion of DC Resistivity And Traveltime Data on Unstructured Mesh

A major challenge in interpreting geophysical data is how to derive consistent three-dimensional (3D) earth models of different physical properties from spatially and temporally limited measurements. Joint inversion with cross-gradient constraints is an approach to find such models by imposing structural similarities between different physical parameters. We have developed a parallel distributed-memory joint inversion code for direct-current (DC) resistivity and traveltime data using the cross-gradient constraint on unstructured mesh. The code utilizes existing E4D framework for parallel forward simulation, distributed storage and computation of the Jacobian matrix of forward operator, and parallel execution of matrix-vector multiplication during inversion. Besides, the joint inversion is solved by nonlinear conjugate gradient algorithm parallelized for DC resistivity and traveltime data. The joint inversion capability of E4D was tested using synthetic data from cross-borehole DC resistivity and traveltime data. The results indicate that the shape and size of the anomalies from the joint inversion are more reliable than those from separate inversions.

58 GEOSCIENCES↗
DOE OSTI2020

A Study on Modeled Wind Speed Errors Using the U.S. Department of Energy Buoys

Pacific Northwest National Laboratory (PNNL) operates two AXYS WindSentinel lidar buoys for the U.S. Department of Energy’s Wind Energy Technologies Office. The purpose of these buoys is to collect hub-height winds and supporting meteorological and oceanographic information to facilitate the development of wind energy in the U.S. waters. The first deployment for one buoy was off the coast of Virginia from December 2014 to May 2016, and the first deployment for the other buoy was off the coast of New Jersey from November 2015 until February 2017. This report describes recent analysis of data collected during these first two deployments. Specifically, we compare hub-height wind speed estimates using Monin-Obukhov Similarity Theory (MOST) to the lidar measurements, and examine how those errors are affected by wind direction, atmospheric stability, wind-wave direction differences, and various measures of the wave-state. The comparisons are done using standard similarity functions based on MOST; including the Businger - Dyer, the Beljaars & Holtslag and the Vickers & Mahrt similarity functions. All models produce large errors over the range of atmospheric stabilities that were observed, with the largest errors occurring for stable flows. The Vickers & Mahrt function resulted in the largest overall bias and standard deviation, while Beljaars & Holtslag function gave the smallest bias and standard deviation due to its better performance under stable conditions. The models perform best under unstable conditions, but even in this regime there is a consistent overestimation of the wind speed of between roughly 0 to 1 ms -1 compared to the lidar measurements. We identify specific metocean conditions (i.e. stability and wind and wave directions) at each of the deployment locations that lead to large errors in MOST predictions. Finally, a coupled ocean-atmosphere model framework was investigated to simulate large errors in weather research forecasting (WRF).

17 WIND ENERGY↗
DOE OSTI2022

Dayflow: CONUS Daily Streamflow Reanalysis, Version 1 (V1)

Dayflow V1 is a historical streamflow reanalysis dataset reconstructed for a 36-year period (1980-2015). The dataset provides both daily and monthly scale streamflow information at about 2.7 million NHDPlusV2 stream reaches in the conterminous US (CONUS). Dayflow is the result of a nationally scalable modeling framework that integrates the simulated runoff from the Variable Infiltration Capacity (VIC) model with the Routing Application for Parallel computatIon of Discharge (RAPID) routing model. Two types of streamflow products, simulated streamflow with or without assimilation of historic US Geological Survey (USGS) streamflow observations, are provided in Dayflow V1. A comprehensive evaluation at 7,526 USGS National Water Information System (NWIS) gauges is performed for both types of streamflow products. The resulting key evaluation metrics are also included in the Dayflow V1 Dataset.

13 HYDRO ENERGY↗
DOE OSTI2020

Experimental protocol for sea level projections from ISMIP6 stand-alone ice sheet models

Projection of the contribution of ice sheets to sea level change as part of the Coupled Model Intercomparison Project Phase 6 (CMIP6) takes the form of simulations from coupled ice sheet–climate models and stand-alone ice sheet models, overseen by the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). This paper describes the experimental setup for process-based sea level change projections to be performed with stand-alone Greenland and Antarctic ice sheet models in the context of ISMIP6. The ISMIP6 protocol relies on a suite of polar atmospheric and oceanic CMIP-based forcing for ice sheet models, in order to explore the uncertainty in projected sea level change due to future emissions scenarios, CMIP models, ice sheet models, and parameterizations for ice–ocean interactions. We describe here the approach taken for defining the suite of ISMIP6 stand-alone ice sheet simulations, document the experimental framework and implementation, and present an overview of the ISMIP6 forcing to be used by participating ice sheet modeling groups.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Analysis of Power Electronic Solutions for Wide-Range Wave Energy Input

Wave energy is a largely untapped energy source with the potential to provide 290 TWh per year to the continental United States. As a new technology, wave energy converters are less efficient and reliable compared to established energy sources, leaving the vast energy of the sea largely unutilized. Wave energy is uniquely inconsistent, with large fluctuations. Paired with harsh operating environments and difficulty in repair, the power electronics designed to harvest wave energy need to be reliable, cost-effective, and able to work in a wide operating band. The literature on existing power electronics in wave energy converters and lessons learned across other industries provide a framework of topologies to simulate. WEC-Sim, an open-source hydrodynamic platform in MATLAB Simulink developed by the National Laboratory of the Rockies, was used to simulate a wide range of power electronics under a variety of wave conditions. This paper provides a comprehensive investigation into existing wave energy power electronics converters, state-of-the-art topologies that can be applied to wave energy, and novel solutions to wave conversion.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

NMPC for Mode-Switching Operation of Reversible Solid Oxide Cell Systems

Solid oxide cells (SOCs) are a promising dual-mode technology that generates hydrogen through high-temperature water electrolysis and generates power through a fuel cell reaction that consumes hydrogen. Reversible operation of SOCs requires a transition between these two modes for hydrogen production setpoints as the demand and price of electricity fluctuate. Moreover, a well-functioning control system is important to avoid cell degradation during mode-switching operation. In this work, we apply nonlinear model predictive control (NMPC) to an SOC module and supporting equipment and compare NMPC performance to classical proportional integral (PI) control strategies, while ramping between the modes of hydrogen and power production. While both control methods provide similar performance in many metrics, NMPC significantly reduces cell thermal gradients and curvatures (mixed spatial-temporal partial derivatives) during mode switching. A dynamic process flowsheet of the reversible SOC system was developed in the open-source, equation-based IDAES modeling framework. Our IDAES dynamic simulation results show that NMPC can ramp the SOC system between hydrogen and power production targets within short mode-switching times. Moreover, NMPC can comply with operating limits in the SOC system more effectively than PI, and only NMPC can directly enforce user-specified limits for mixed spatial-temporal partial derivatives of temperature. This allows for management of the trade-off be-tween operating efficiency and cell degradation, which is dependent on these temperature curvatures.

Li, Mingrui↗
DOE OSTI2019

Electromagnetics Simulations with Vector-valued Finite Elements in MOOSE

Vector-valued finite elements have become an essential part of modern electromagnetic (EM) numerical simulation. This talk will describe the development and initial validation of a MOOSE-based EM solver application, ELK (Electromagnetics Library for Kinetics and fluids), that utilizes the recently added vector-valued finite element system. The primary motivation for the development of ELK was to study EM wave propagation in radio-frequency (rf) driven systems, such as those encountered in the study of antenna structures, waveguides, and plasma physics. In simulating these scenarios with traditional scalar-valued finite elements, boundary conditions on tangential field components can sometimes be difficult to define, and field discontinuities encountered with certain geometric or structural features (e.g., abrupt corners and inhomogeneous materials) can often lead to non-physical results. Vector-valued elements can largely mitigate and even eliminate these difficulties. ELK verification and validation efforts thus far have focused on standard wave propagation benchmarks for antennas and waveguides, and the impact of vector-valued elements on the implementation and results of these simulations will be described. Progress using ELK for plasma physics simulation, where inhomogeneous simulation domains are commonplace, will also be discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Time-Disciplined Non-PLL Active Synchronization for Grid-Forming Inverters: Preprint

Grid Forming Inverters (GFMIs) are a critical component of the future Grid. One major challenge with GFMIs is synchronization. Currently, most research in GFMI synchronization is focused around specific transition types but this is only a subcategory of the possible events facing GFMIs. In this paper, we propose categorizing all possible synchronization events into four categories: (1) Initial Grid Synchronization, (2) Ride-Through, (3) Blackstart Coordination, (4) Grid Reconnection Coordination. By categorizing events, we can better compare and integrate different control strategies under specific test conditions as demonstrated in simulation. This paper’s framework helps guide GFMI primary-level control design to further research towards the 100% inverter-dominant Grid.

blackstart↗
DOE OSTI2021

Time Disciplined Non-PLL Active Synchronization for Grid Forming Inverters

Grid forming inverters (GFMIs) are a critical component of the future grid. One major challenge with GFMIs is synchronization. Currently, most research in GFMI synchronization is focused around specific transition types but this is only a subcategory of the possible events facing GFMIs. In this paper, we propose categorizing all possible synchronization events into four categories: (1) initial grid synchronization, (2) ride-through, (3) blackstart coordination, (4) grid reconnection coordination. By categorizing events, we can better compare and integrate different control strategies under specific test conditions as demonstrated in simulation. This paper's framework helps guide GFMI primary-level control design to further research towards the 100% inverter-dominant grid.

blackstart↗
DOE OSTI2025

Towards An Updated Simulation of the Booster Neutrino Beam

Precise, accurate neutrino flux predictions for neutrino beam experiments are crucial for physics results. Flux predictions for the Booster Neutrino Beam at Fermilab were first published in 2009 by the MiniBooNE collaboration. It is no longer possible to run the simulation used to create these predictions due to outdated software versions. In this exciting period for the Short-Baseline Neutrino Program at Fermilab, with both the far detector (ICARUS) and the near detector (SBND) operating, an updated flux model for the BNB is necessary. For the purpose of using a dynamic, up-to-date beam model, the SBN program has created a new simulation named G4BNB. This framework contains new features, such as a full neutrino ancestry to handle hadron production systematics with more precision, and the inclusion of neutral mesons from the proton-beryllium scatter to allow the study of exotic BSM scenarios.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Petri Nets for Adversarial Models using Monte Carlo Simulation

We describe an object-oriented framework and user interface for developing models of adversarial scenarios, implemented as stochastic Petri nets. The framework is implemented using the statistical computing language R, and provides facilities for eliciting models from subject matter experts, displaying models graphically, simulating scenario execution, and presenting scenario outcomes and statistics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Optimizing district energy systems by integrating Borehole Thermal Energy Storage Using a Mixed-Integer Linear Programming g-function framework with a Multi-Timescale Rolling Horizon method

Shallow geothermal has gained increasing attention in recent years; however, a reliable framework for its accurate incorporation into large-scale energy system optimization remains lacking. This study proposes a Mixed-Integer Linear Programming (MILP) framework combined with the g-function approach to integrate Borehole Thermal Energy Storage (BTES) technology into energy system optimization. Validation against a Modelica-based reservoir network simulation demonstrates that the proposed framework effectively captures the ground thermal response under varying energy loads and accurately estimates the borefield energy supply. To enhance scalability, a Rolling Horizon with Multi-Timescale (RH-MTS) method is further introduced, reducing computational time by 73 % for the 1-year optimization model with only minor loss of optimality. The framework is demonstrated through the case study of the UC Berkeley campus. Results indicate that BTES is a cost-effective and low-carbon solution: two borefields comprising 382 boreholes can meet 8.0 % and 6.6 % of the total campus heating and cooling demand, respectively, at an average energy rate of 0.70–0.77 USD/kWh and carbon intensity of 0.54 kg-CO2/kWh. Short-term analysis reveals a 35%–65% decline in BTES energy flow after 3–6 months of continuous heating/cooling operation, while long-term simulation shows that annual energy production of BTES can vary by up to 12.0 % after four years before stabilizing. Overall, this study develops a novel optimization framework that couples physics-based g-function method with MILP optimization framework, thereby advancing methodological development for shallow-geothermal integration and providing actionable guidance for BTES deployment in district-energy systems.

Yang, Jiahui↗
DOE OSTI2021

Search for efficient formulations for Hamiltonian simulation of non-Abelian lattice gauge theories

Hamiltonian formulation of lattice gauge theories (LGTs) is the most natural framework for the purpose of quantum simulation, an area of research that is growing with advances in quantum-computing algorithms and hardware. It, therefore, remains an important task to identify the most accurate, while computationally economic, Hamiltonian formulation(s) in such theories, considering the necessary truncation imposed on the Hilbert space of gauge bosons with any finite computing resources. This paper is a first step toward addressing this question in the case of non-Abelian LGTs, which further require the imposition of non-Abelian Gauss’s laws on the Hilbert space, introducing additional computational complexity. Focusing on the case of SU(2) LGT in 1+1 dimensions coupled to one flavor of fermionic matter, a number of different formulations of the original Kogut-Susskind framework are analyzed with regard to the dependence of the dimension of the physical Hilbert space on boundary conditions, system’s size, and the cutoff on the excitations of gauge bosons. The impact of such dependencies on the accuracy of the spectrum and dynamics obtained from a Hamiltonian computation is examined, and the (classical) computational-resource requirements given these considerations are studied. Besides the well-known angular-momentum formulation of the theory, the cases of purely fermionic and purely bosonic formulations (with open boundary conditions), and the Loop-String-Hadron formulation are analyzed, along with a brief discussion of a Quantum-Link-Model formulation of the same theory. Clear advantages are found in working with the Loop-String-Hadron framework which implements non-Abelian Gauss’s laws a priori using a complete set of gauge-invariant operators. Although small lattices are studied in the numerical analysis of this work, and only the simplest algorithms are considered, a range of conclusions will be applicable to larger systems and potentially to higher dimensions. Future studies will extend this investigation to the analysis of resource requirements for quantum simulations of non-Abelian LGT, with the goal of shedding light on the most efficient Hamiltonian formulation of gauge theories of relevance in nature

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Structure and dynamics of aqueous solutions from PBE-based first-principles molecular dynamics simulations

Establishing an accurate and predictive computational framework for the description of complex aqueous solutions is an ongoing challenge for density functional theory based first-principles molecular dynamics (FPMD) simulations. In this context, important advances have been made in recent years, including the development of sophisticated exchange-correlation functionals. On the other hand, simulations based on simple generalized gradient approximation (GGA) functionals remain an active field, particularly in the study of complex aqueous solutions due to a good balance between the accuracy, computational expense, and the applicability to a wide range of systems. In such simulations we often perform them at elevated temperatures to artificially “correct” for GGA inaccuracies in the description of liquid water; however, a detailed understanding of how the choice of temperature affects the structure and dynamics of other components, such as solvated ions, is largely unknown. In order to address this question, we carried out a series of FPMD simulations at temperatures ranging from 300 to 460 K for liquid water and three representative aqueous solutions containing solvated Na + , K + , and Cl - ions. We show that simulations at 390–400 K with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional yield water structure and dynamics in good agreement with experiments at ambient conditions. Simultaneously, this computational setup provides ion solvation structures and ion effects on water dynamics consistent with experiments. These results suggest that an elevated temperature around 390–400 K with the PBE functional can be used for the description of structural and dynamical properties of liquid water and complex solutions with solvated ions at ambient conditions.

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

Multiscale Universal Interface: A concurrent framework for coupling heterogeneous solvers

Graphical abstract: - Abstract: Concurrently coupled numerical simulations using heterogeneous solvers are powerful tools for modeling multiscale phenomena. However, major modifications to existing codes are often required to enable such simulations, posing significant difficulties in practice. In this paper we present a C++ library, i.e. the Multiscale Universal Interface (MUI), which is capable of facilitating the coupling effort for a wide range of multiscale simulations. The library adopts a header-only form with minimal external dependency and hence can be easily dropped into existing codes. A data sampler concept is introduced, combined with a hybrid dynamic/static typing mechanism, to create an easily customizable framework for solver-independent data interpretation. The library integrates MPI MPMD support and an asynchronous communication protocol to handle inter-solver information exchange irrespective of the solvers' own MPI awareness. Template metaprogramming is heavily employed to simultaneously improve runtime performance and code flexibility. We validated the library by solving three different multiscale problems, which also serve to demonstrate the flexibility of the framework in handling heterogeneous models and solvers. In the first example, a Couette flow was simulated using two concurrently coupled Smoothed Particle Hydrodynamics (SPH) simulations of different spatial resolutions. In the second example, we coupled the deterministic SPH method with the stochastic Dissipative Particle Dynamics (DPD) method to study the effect of surface grafting on the hydrodynamics properties on the surface. In the third example, we consider conjugate heat transfer between a solid domain and a fluid domain by coupling the particle-based energy-conserving DPD (eDPD) method with the Finite Element Method (FEM)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗