Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation package”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

End-Use Savings Shapes Upgrade Package Documentation: Wall and Roof Insulation and New Windows

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock and ComStock models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) upgrade measures, or upgrades. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility upgrade applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on an upgrade package of three end-use savings shapes upgrades - Window Replacement, Exterior Wall Insulation, and Roof Insulation, which we will refer to collectively as the "High Efficiency Envelope" package. More details on the individual upgrades can be found on the ComStock Measures Documentation page. An upgrade package applies two or more EUSS upgrades to a single building model simulation. Since ComStock is a bottom-up physics-based model, an upgrade package will go beyond aggregating or summing the individual upgrade results and produce novel results by simulating interactions between the upgrades. For example, pairing an envelope upgrade with an electrification upgrade would likely result in higher savings results than the sum of these upgrades individually, and the size of the heating, ventilating, and air conditioning (HVAC) equipment may be reduced if the envelope upgrade reduces the loads significantly.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

End-Use Savings Shapes Upgrade Package Documentation: LED Lighting, HP-RTU and ASHP-Boiler

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock and ComStock models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) upgrade measures, or upgrades. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility upgrade applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on an upgrade package of three end-use savings shapes upgrades - Light Emitting Diode (LED) Lighting, Heat Pump Rooftop Unit (RTU) (HP-RTU), and Air-Source Heat Pump (ASHP) Boiler, which we will refer to collectively as the "Interior Lighting and Heat Pump" package. More details on the individual upgrades can be found on the ComStock Measures Documentation page. An upgrade package applies two or more EUSS upgrades to a single building model simulation. Since ComStock is a bottom-up physics-based model, an upgrade package will go beyond aggregating or summing the individual upgrade results and produce novel results by simulating interactions between the upgrades. For example, pairing an envelope upgrade with an electrification upgrade would likely result in higher savings results than the sum of these upgrades individually, and the size of the heating, ventilating, and air conditioning (HVAC) equipment may be reduced if the envelope upgrade reduces the loads significantly.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

ZERNIPAX: A fast and accurate Zernike polynomial calculator in Python

Zernike polynomials serve as an orthogonal basis on the unit disc, and have proven to be effective in optics simulations, astrophysics, and more recently in plasma simulations. Unlike Bessel functions, Zernike polynomials are inherently finite and smooth at the disc center (r=0), ensuring continuous differentiability along the axis. This property makes them particularly suitable for simulations, requiring no additional handling at the origin. We developed ZERNIPAX, an open-source Python package capable of utilizing CPU/GPUs, leveraging Google's JAX package and available on GitHub as well as the Python software repository PyPI. Furthermore, our implementation of the recursion relation between Jacobi polynomials significantly improves computation time compared to alternative methods by use of parallel computing while still performing more accurately for high-mode numbers.

Astrophysics↗
DOE OSTI2022

MPLBM-UT: Multiphase LBM library for permeable media analysis

MPLBM-UT is a specialized lattice-Boltzmann library that makes running single- and two-phase flow simulations in porous media accessible to everyone. We provide a suite of tools to pre-process computational domains for simulation, to set up custom boundary conditions, to run simulations, to post-process simulation outputs, and to visualize simulation results and data. All of these tools are easily accessible to users through the mplbm_utils Python package included in and automatically installed with MPLBM-UT. The high-performance, highly parallel library Palabos is used as the solver backend. MPLBM-UT is easily deployed in a variety of systems, from laptops to supercomputer clusters. MPLBM-UT also features multiple examples and benchmark templates that allow for fast prototyping of different porous media problems. We also provide an interface for reading in different file types and downloading domains from the Digital Rocks Portal to perform simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Predictive Contaminant Transport Simulation with P2R Model for the Composite Analysis Limited Source Sensitivity Case.

The Plateau to River (P2R) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by the Central Plateau Cleanup Company at the Hanford Site in Washington State. Figure 1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the Composite Analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3 documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application. Simulations conducted to support the dose calculations required by the CA are documented in ECF-HANFORD-19-0119, Predictive Flow Simulation with the P2R Model for the Composite Analysis Base Case and ECF-HANFORD-19-0120, Contaminant Transport Simulation with the P2R Model for the Composite Analysis Base Case.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

A symbolic framework to obtain mid-fidelity models of flexible multibody systems with application to horizontal-axis wind turbines

Abstract. The article presents a symbolic framework (also called computer algebra program) that is used to obtain, in symbolic mathematical form, the linear and nonlinear equations of motion of a mid-fidelity multibody system including rigid and flexible bodies. Our approach is based on Kane's method and a nonlinear shape function representation for flexible bodies. The shape function approach does not represent the state of the art for flexible multibody dynamics but is an effective trade-off to obtain mid-fidelity models with few degrees of freedom, taking advantage of the separation of space and time. The method yields compact symbolic equations of motion with implicit account of the constraints. The general and automatic framework facilitates the creation and manipulation of models with various levels of complexity by adding or removing degrees of freedom. The symbolic treatment allows for analytical gradients and linearized equations of motion. The linear and nonlinear equations can be exported to Python code or dedicated software. There are multiple applications, such as time domain simulation, stability analyses, frequency domain analyses, advanced controller design, state observers, and digital twins. In this article, we describe the method we used to systematically generate the equations of motion of multibody systems and present the implementation of the framework using the Python package SymPy. We apply the framework to generate illustrative land-based and offshore wind turbine models. We compare our results with OpenFAST simulations and discuss the advantages and limitations of the method. The Python implementation is provided as an open-source project.

Branlard, Emmanuel (ORCID:0000000277506128)↗
DOE OSTI2023

pyTDGL: Time-dependent Ginzburg-Landau in Python

Time-dependent Ginzburg-Landau (TDGL) theory is a phenomenological model for the dynamics of superconducting systems. Due to its simplicity in comparison to microscopic theories and its effectiveness in describing the observed properties of the superconducting state, TDGL is widely used to interpret or explain measurements of superconducting devices. Here, we introduce pyTDGL, a Python package that solves a generalized TDGL model for superconducting thin films of arbitrary geometry, enabling simulations of vortex and phase dynamics in mesoscopic superconducting devices. pyTDGL can model the nonlinear magnetic response and dynamics of multiply connected films, films with multiple current bias terminals, and films with a spatially inhomogeneous critical temperature. We demonstrate these capabilities by modeling quasi-equilibrium vortex distributions in irregularly shaped films, and the dynamics and current-voltage-field characteristics of nanoscale superconducting quantum interference devices (nanoSQUIDs).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Large-scale application of free energy perturbation calculations for antibody design

Abstract Alchemical free energy perturbation (FEP) is a rigorous and powerful technique to calculate the free energy difference between distinct chemical systems. Here we report our implementation of automated large-scale FEP calculations, using the Amber software package, to facilitate antibody design and evaluation. In combination with Hamiltonian replica exchange, our FEP simulations aim to predict the effect of mutations on both the binding affinity and the structural stability. Importantly, we incorporate multiple strategies to faithfully estimate the statistical uncertainties in the FEP results. As a case study, we apply our protocols to systematically evaluate variants of the m396 antibody for their conformational stability and their binding affinity to the spike proteins of SARS-CoV-1 and SARS-CoV-2. By properly adjusting relevant parameters, the particle collapse problems in the FEP simulations are avoided. Furthermore, large statistical errors in a small fraction of the FEP calculations are effectively reduced by extending the sampling, such that acceptable statistical uncertainties are achieved for the vast majority of the cases with a modest total computational cost. Finally, our predicted conformational stability for the m396 variants is qualitatively consistent with the experimentally measured melting temperatures. Our work thus demonstrates the applicability of FEP in computational antibody design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Correcting aberrations of a transverse-field neutron resonance spin echo instrument

Neutron resonance spin echo (NRSE) technique has the potential to increase the Fourier time and energy resolution in neutron scattering by using radio frequency (rf) neutron spin-flippers. However, aberrations arising from variations in the neutron path length between the rf flippers reduce the polarization. Here, we develop and test a transverse static-field magnet, a series of which are placed between the rf flippers, to correct for these aberrations. The prototype correction magnet was both simulated in an NRSE beamline using McStas, a Monte Carlo neutron ray-tracing software package, and measured using neutrons. The results from the prototype demonstrate that this static-field design corrects for transverse-field NRSE aberrations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Tungsten–carbon surface evolution and erosion modeling for a small angle slot divertor in DIII-D

Abstract We modeled tungsten–carbon mixed surface evolution, sputtering erosion, and transport for the tungsten coated region of a small angle slot (SAS) divertor design for the DIII-D tokamak. This divertor concept aims to achieve a closed slot dissipative plasma to minimize heat load and surface erosion, and to study high- Z material performance. Our advanced simulations use coupled ITMC-DYN material mixing/response and 3D full kinetic REDEP/WBC erosion/redeposition code packages, with divertor plasma solution from the SOLPS-ITER package with 4 MW power input. The SAS design geometry and resulting in-slot plasma parameters cause significant differences in predicted sputter/transport from a conventional divertor. For 2% C/D incident plasma ratio, an equilibrium mixed C/W surface is attained at ∼30 s of discharge, from wall sputtered carbon transported to the 10 cm long tungsten divertor region. Tungsten remains exposed to the plasma, but the evolved surface composition varies with different C/D ratios. Tungsten is primarily sputtered from the mixed surface by impinging carbon ions in the +1 to +4 charge states, with some self-sputtering. Redeposition of sputtered tungsten to the divertor is significant, ∼80% along the higher plasma temperature attached plasma SAS entrance region, but this is less than the typically near-unity values for a conventional divertor. Plasma-incident carbon is highly backscattered (∼50%) from the mixed surface, with little redeposition (<10%); this helps maintain tungsten in the surface sputter zone. Carbon is mainly sputtered from the mixed surface by D + ions, also with low redeposition (∼10%–30%). Finally, the modeling shows non-zero but low sputtered tungsten current from the divertor to the core plasma direction. These results appear favorable for effective testing of a tungsten-containing SAS divertor in DIII-D, and extrapolation of mixed-material evolution/response findings to the analogous low- Z /high- Z , Be/W, ITER plasma facing system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Calculation of Groundwater Pathway Radiological Dose for the Hanford Site Composite Analysis Null-Space Model Carlo Flow Model Set

The purpose of this environmental calculation file (ECF) is to present the results of the exposure route-specific and total radiological dose assessments for the groundwater pathway based on the null space Monte Carlo (NSMC) groundwater concentrations as a part an uncertainty analysis for the updated Hanford Site Composite Analysis (CA). The Plateau-to River (P2R) Groundwater Model (CP-57037, Model Package Report: Plateau to River Groundwater Model Version 8.3) is the CA base case that simulates the fate and transport of radiological contaminants within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. An NSMC analysis was performed to identify and quantify the potential uncertainties associated with the P2R Model. The result of the NSMC analysis is a set of flow and transport simulations that provide an estimated range of possible outcomes that are used to quantify the uncertainty associated with the simulated base case concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Mechanistic Modeling of Creep and Aging for Time Extrapolation in LPBF 316 SS - LANL

A polycrystal constitutive model was extended to quantify the effects heterogeneity in the microstructure, porosity, and internal stresses on the mechanical response (creep, tensile) of 316H steel processed by additive manufacturing. The starting point is a recently developed and fully calibrated advanced microstructuresensitive mechanistic constitutive model fit to literature data pertaining to conventionally processed steels. The model is embedded within the full field fast Fourier transform framework (LApx) and several hundred simulations were performed. Synthetic microstructures with different grain aspect ratios were generated using the Neper package to study the role of grain morphology on the mechanical responses. The effect of initial cell interior dislocation density, cell-wall density, precipitate content, porosity, and internal stresses on creep rates are quantified for different polycrystalline microstructures. The study allowed to prospectively quantify the role of distinct features of the microstructure on the mechanical response under tension and creep loading scenarios. Critically it is found that minor anisotropies in the tensile response of samples, emanating from the aspect ratio of grains correlated into significant differences in creep rate as a function of loading directions. The study also allows to identify initial microstructure porosities below which the materials response only moderately degraded by the initial microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hands-on with OWL: the Oak–Ridge Wang–Landau Monte Carlo software suite

The Oak–Ridge Wang–Landau (OWL) package is an open-source scientific software specialized for large-scale, Monte Carlo simulations for the study of materials properties at finite temperature. In this paper, we discuss the main features and capabilities of OWL, followed by detailed descriptions of building and running the code. The readers will be guided through the usage and functionality of the code with a few hands-on examples. This paper is based on a tutorial on OWL given at the 32nd Center for Simulational Physics Workshop on Recent Developments in Computer Simulation Studies in Condensed Matter Physics.

Li, Ying Wai↗
DOE OSTI2020

A scalable matrix-free spectral element approach for unsteady PDE constrained optimization using PETSc/TAO

In this work, we provide a new approach for the efficient matrix-free application of the transpose of the Jacobian for the spectral element method for the adjoint-based solution of partial differential equation (PDE) constrained optimization. This results in optimizations of nonlinear PDEs using explicit integrators where the integration of the adjoint problem is not more expensive than the forward simulation. Solving PDE constrained optimization problems entails combining expertise from multiple areas, including simulation, computation of derivatives, and optimization. The Portable, Extensible Toolkit for Scientific computation (PETSc) together with its companion package, the Toolkit for Advanced Optimization (TAO), is an integrated numerical software library that contains an algorithmic/software stack for solving linear systems, nonlinear systems, ordinary differential equations, differential algebraic equations, and large-scale optimization problems and, as such, is an ideal tool for performing PDE-constrained optimization. This paper describes an efficient approach in which the software stack provided by PETSc/TAO can be used for large-scale nonlinear time-dependent problems. Time integration can involve a range of high-order methods, both implicit and explicit. The PDE-constrained optimization algorithm used is gradient-based and seamlessly integrated with the simulation of the physical problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

ELM model simulations of Plum Island Ecosystems LTER low marsh site 2018-2020

Model simulations using the Department of Energy's Energy Exascale Earth System Model (E3SM) land model (ELM) with improved capabilities to represent vegetation response to salinity and inundation. The simulations were conducted for a tidal salt marsh at Plum Island Ecosystems Long Term Ecological Research (LTER) site near Rowley, Massachusetts, USA; the site is a low marsh dominated by Spartina alterniflora. The model was forced with site-specific meteorology, salinity and tidal cycles from 2018-2020. Four sets of model simulations are included and described below:1. Parameterization of the salinity response function. These simulations tested different combinations of values for optimal salinity and salinity tolerance.2. Model evaluation. This comparison conducted simulations using the default model, the salinity function only, the submergence function only, and both the salinity and submergence functions. 3. Salinity scenarios. These simulations used the 2018 salinity input data varied by -5 to +10 ppt salinity.4. Water level scenarios. These simulations used the tide height varied by -10 to +50 cm. These simulations were used to demonstrate how the salinity and submergence functions better represent carbon uptake by tidal salt marshes.The data package includes netCDF files used as forcing files for tide height and salinity, one for each year 2018-2020 at observed salinity concentrations, and an additional three forcing files in which salinity concentrations were varied 5 ppt lower, 5 ppt higher, and 10 ppt higher than the measured 2018 time series. Also included are python scripts for creating forcing files, plain text parameter and command files for running simulations, model outputs in netCDF format, and python scripts for visualizing outputs. Code for the modified E3SM model is archived in Sulman et al 2023 at doi:10.15485/1991625. More detail about files is provided in the README.md file.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

MPC4CLR (Model-Predictive-Control-for-Critical-Load-Restoration-in-Power-Distribution-Systems) [SWR-22-24]

Model predictive control (MPC) is a system or process control technique for making decisions under uncertainty via rolling look-ahead optimizations at each control step where only the current step decisions are applied, and the rest are discarded. In this work, we developed an MPC for a critical load restoration (CLR) in power distribution systems to recover system service (electricity delivery) following an extreme event-triggered substation outage. The method considers the problem of controlling distributed energy resources (DERs) of the distribution system with the objective of achieving maximum load pick up while satisfying distribution network flow and voltage constraints. A linearized optimal power flow (OPF) model is employed to represent the physics of the network. The problem formulation is augmented with a ramping (up) reserve product for the DERs to ensure improved and upward monotonic load restoration as time evolves. Simulation analysis and performance tests are performed using a modified IEEE 13-bus test feeder integrated with wind, solar, microturbine, and energy storage battery. The software is developed using various software packages in Julia and Python. The MPC model is implemented using the JuMP optimization language in Julia while the data analytics including renewable generation and load demand forecasts, running the MPC simulation and visualizations is performed in Python.

Eseye, Abinet Tesfaye↗
DOE OSTI2022

Achieving 50% Energy Savings in Chicago Homes: A Case Study for Advancing Equity and Climate Goals

Since 2019, the National Renewable Energy Laboratory (NREL) and Elevate have collaborated to identify pathways to deep energy retrofits in Chicago's housing stock, document equity implications and co-benefits of this transition, and validate the findings by implementing retrofits in real Chicago homes. This document summarizes our analysis process to model advanced retrofit packages that lead to greater than 50% energy savings in Chicago homes. Based on these findings, we have also developed a roadmap with the City to guide implementation, and are deploying the recommended retrofit packages in real Chicago homes to realize these energy savings. This work was developed in collaboration with two key stakeholders - the City of Chicago and Commonwealth Edison (ComEd) - and funded by the U.S. Department of Energy (DOE). NREL's Residential Buildings team maintains the best-in-class ResStock™ energy model of the U.S. residential building stock. For this work, we calibrated ResStock to Chicago's unique local housing stock to accurately simulate energy use in Chicago homes both for current conditions and with various retrofit scenarios. We simulated a wide range of potential building retrofits covering all aspects of residential energy use and then grouped these into packages based on energy and utility bill savings. ResStock can model diverse building types and housing characteristics, so we're able to observe the range of outcomes that might occur when these upgrades are deployed across the entire housing stock. We can then estimate potential energy savings from an advanced retrofit program on Chicago's housing stock by comparing the modeled energy use before versus after a retrofit. This novel version of ResStock, calibrated to Chicago with data from Elevate, can help City officials, ComEd, and other partners plan for community-scale decarbonization via residential retrofits. Specifically, this work contributes the following project goals: Develop a building retrofit prioritization strategy for Chicago single-family and 2- to 4-unit buildings; Identify neighborhoods and home types that have the highest potential for savings from electrification; and Assess the impact of advanced building retrofits on energy use, utility bills, and CO2 emissions at the city and building level. Although this study is specific to Chicago, its methods and learnings are applicable across the United States. These findings are especially notable for heat pumps and electrification retrofits in cold climates.

building energy modeling↗
DOE OSTI2021

Experimental inference of neutral and impurity transport in Alcator C-Mod using high-resolution x-ray and ultra-violet spectra

We present experimental inferences of cross-field impurity transport coefficients for AlcatorC-Mod plasmas without edge-localized modes, leveraging a novel forward model for the entire Ca Kα spectrum, including satellite lines within the spectral range, to compare to high-resolution x-ray imaging crystal spectroscopy (XICS). These measurements are complemented by extreme ultra-violet (EUV) spectroscopy that constrains transport closer to the edge. Using new atomic data sets for both XICS and EUV analysis has enabled consideration of line ratios across both spectral ranges and has increased the accuracy of inferred transport coefficients. Inclusion of charge exchange between edge thermal neutrals and impurities is shown to be extremely important in C-Mod pedestals. Atomic D neutral densities from experimental D Lyα measurements at the midplane are compared to SOLPS-ITER simulations, finding good agreement. Bayesian inferences of impurity transport coefficients are presented for L-, EDA H-, and I-mode discharges, making use of the Aurora package for forward modeling and combining our spectroscopic constraints. Experimentally inferred diffusion profiles are found to match turbulent transport models at midradius within uncertainties, using both quasilinear gyro-fluid TGLF SAT-1 and nonlinear ion-scale gyrokinetic CGYRO simulations. Furthermore, significant discrepancies in convection are observed in some cases, suggesting difficulties in predictions of flat or hollow impurity profiles.

neutrals↗