Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic load model”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multiphysics Simulation of Supercritical CO2 Gasification for Hydrogen Production

A bench-scale supercritical CO2 gasifier is being considered by the U.S. Department of Energy’s National Energy Technology Laboratory (NETL). A reacting computational fluid dynamics (CFD) model was developed using ANSYS/Fluent in order to investigate the impact of various operating parameters on syngas composition. The model was validated using available data in the literature. Simulations were conducted over a range of coal-slurry loadings, reaction temperatures, system pressures, and oxygen mass flowrates with the goal of optimizing hydrogen production for the considered reactor. The results from the simulations identified a positive correlation between slurry loading and hydrogen production, while reactor temperature and pressure had a limited impact on H2 production at the conditions of interest. Additionally, the simulations indicated that oxygen mass flowrate has less of an impact on hydrogen production at higher slurry loadings.

Prlina, Alexander↗
DOE OSTI2022

Numerical modeling and parametric study of a dual purpose underground thermal battery

We report integrating thermal energy storage with building energy systems can enable flexible building electric demands at buildings to help mitigate the mismatch between electricity supply and demand. A novel building heating and cooling system that integrates a dual-source heat pump with hybrid thermal storage named dual-purpose underground thermal battery (DPUTB) has been developed for reshaping building electric demands. The proposed DPUTB integrated geothermal heat pump system is an original innovation that enables Grid-interactive Energy Efficient buildings. This paper focuses on the study of the novel DPUTB. The DPUTB works as both a thermal storage tank (an inner tank) and a ground heat exchanger (an outer tank separated from the inner tank by the insulation material). High fidelity and computationally effective models are needed to predict the performance of the novel DPUTB. This study has developed a simplified dynamic model for the DPUTB according to heat transfer and energy conservation principles and validated it by using experimental data obtained from testing a small-scale DPUTB apparatus. A parametric study was conducted to identify a design that can achieve the target thermal storage performance of load shift and energy efficiency. The parametric study results show that the inner tank shell thermal conductivity and the phase change material's melting point are the two most influencing factors on the performance of the DPUTB. One single full-size DPUTB with the identified design could provide 1-ton cooling (3.51 kW) with the supply temperature lower than 11 °C for 4 hours in summer after being fully charged in 8 hours. The inner tank filled with phase change material is for cooling thermal storage as a latent tank in the design. However, its capacity can be as high as 60 MJ as a sensible water tank for heating storage in winter. In the future study, the DPUTB model will be incorporated into the dual-source heat pump system for evaluating the overall system performance of demand side management in the long term.

25 ENERGY STORAGE↗
DOE OSTI2021

An asymptotically compatible treatment of traction loading in linearly elastic peridynamic fracture

Meshfree discretizations of state-based peridynamic models are attractive due to their ability to naturally describe fracture of general materials. However, two factors conspire to prevent meshfree discretizations of state-based peridynamics from converging to corresponding local solutions as resolution is increased: quadrature error prevents an accurate prediction of bulk mechanics, and the lack of an explicit boundary representation presents challenges when applying traction loads. Herein, we develop a reformulation of the linear peridynamic solid (LPS) model to address these shortcomings, using improved meshfree quadrature, a reformulation of the nonlocal dilatation, and a consistent handling of the nonlocal traction condition to construct a model with rigorous accuracy guarantees. In particular, these improvements are designed to enforce discrete consistency in the presence of evolving fractures, whose a priori unknown location render consistent treatment difficult. In the absence of fracture, when a corresponding classical continuum mechanics model exists, our improvements provide asymptotically compatible convergence to corresponding local solutions, eliminating surface effects and issues with traction loading which have historically plagued peridynamic discretizations. When fracture occurs, our formulation automatically provides a sharp representation of the fracture surface by breaking bonds, avoiding the loss of mass. We provide rigorous error analysis and demonstrate convergence for a number of benchmarks, including manufactured solutions, free-surface, nonhomogeneous traction loading, and composite material problems. Finally, we validate simulations of brittle fracture against a recent experiment of dynamic crack branching in soda-lime glass, providing evidence that the scheme yields accurate predictions for practical engineering problems.

42 ENGINEERING↗
DOE OSTI2020

Thermal Analysis Evaluations Using the Dry Cask Simulator

Different thermal analysis models were developed to simulate the dry cask simulator (DCS). The DCS is an experiment designed to simulate dry storage of a single boiling water reactor fuel assembly under a variety of heat loads and internal pressures. The DCS was set up and tested in both a vertical and horizontal configuration to determine cladding temperatures in vertical and horizontal dry cask storage systems. The models included a detailed STAR-CCM+ model with the fuel assembly geometry explicitly modeled, a porous STAR-CCM+ model with the fuel assembly geometry modeled as porous media region with calculated effective properties, and a COBRA-SFS model. COBRA-SFS is a thermal-hydraulic code developed for steady-state and transient analysis of multi assembly spent-fuel storage and transportation systems. STAR-CCM+ is a commercial computation fluid dynamics (CFD) code. Both a detailed and porous STARCCM+ model were developed to look at the effective thermal conductivity (keff) approach to modeling a fuel assembly. A keff fuel model is typically model in CFD thermal analyses due to its significantly lower computational costs. The models were run for a combination of low and high canister pressures (100 kPa and 800 kPa) and low and high internal heat loads (0.5 kW and 5 kW). Results from all three models were compared against experimental data taken from the DCS for the peak cladding temperature (PCT) and inlet air mass flow.

Suffield, Sarah R.↗
DOE OSTI2022

Enhanced Inertial Support: Modeling Fast Frequency Response Controls for Energy Storage System Inverters

Reduction of power system inertia is a growing concern for grid operators given the transition under way from conventional fossil generators to inverter-based resources. One key concern is the rapid and severe dip in power system frequency following a large disturbance which could trigger remedial actions including operation of under frequency load shedding (UFLS) relays. A potential solution is to utilize battery energy storage systems (BESS) to provide fast frequency response(FFR) to mitigate the system frequency impact of low system inertia. In this work, we investigate different control parameters to achieve modeled FFR response from BESS in the simulated US Western Interconnection. First, a very low inertia model is developed from an existing positive sequence model of the Western Interconnection. Afterwards, we tune the BESS dynamic response parameters to eliminate any risks of reaching UFLS relay thresholds following the loss of large generators. Finally, we investigate the relative impact of having large BESS units in comparison of multiple smaller units in the same area. The results of the analyses indicate that 1 GW of transmission-connected BESS can help mitigate the frequency drop caused by large generators outages in the Western Interconnection, and significantly reduce the risk of under frequency load shedding.

power system dynamic performance, Low Inertia, Fas↗
DOE OSTI2024

HP-FLEX MPC v0.1.0

HP-FLEX MPC is control software developed by Lawrence Berkeley National Laboratory with support from the California Energy Commission (CEC) through EPIC-19-301. HP-FLEX aims to provide load flexibility for heat pumps (HPs) in response to dynamic grid signals (including Time-of-Use, Dynamic Pricing, and Critical Peak Pricing) while maintaining thermostat temperatures within user-specified bounds. The software includes a system-identification module, which models the dynamics of the building envelope with thermostat data, and a control module based on a model predictive controller (MPC) to make optimal decisions. HP-FLEX receives forecasts of outdoor air temperature, solar irradiation, and internal gain (if available), as well as trajectories of energy price, temperature lower and upper bounds over a prediction horizon. It then optimizes heating and cooling capacities to minimize energy cost and peak power (with a user-defined weight on peak power) over the prediction horizon, while maintaining room air temperature within the temperature constraints, and outputs the optimal thermostat setpoints.

Kim, Donghun↗
DOE OSTI2023

Scalable self attraction and loading calculations for unstructured ocean tide models

Self attraction and earth-loading effects are important for accurately modeling global tides. A common approach of handling this forcing is to expand mass anomalies into spherical harmonics, which are scaled by load Love numbers to account for elastic earth deformation. We investigate two different approaches to perform these calculations for ocean models that employ unstructured meshes and distributed memory parallelization. The first approach leverages a highly efficient spherical harmonics library, but requires all-to-one and one-to-all communications and interpolation operations between the unstructured and a structured mesh. This approach is compared to a parallel algorithm that computes the spherical harmonic transformations directly on the unstructured mesh with an all-reduce communication. Here, our results show that although the unstructured mesh calculations are more expensive, the scalability of the unstructured mesh approach allows for more efficient spherical harmonics transforms for high-resolution meshes and large processor counts. This methodology enables the efficient inclusion of tidal dynamics large-scale Earth system model simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dislocation density-based plasticity model from massive discrete dislocation dynamics database

In this work, we present a dislocation density-based strain hardening model for single crystal copper through a systematic coarse-graining analysis of more than 200 discrete dislocation dynamics (DDD) simulations of plastic deformation under uniaxial tension. The proposed constitutive model has two components: a generalized Taylor relation connecting resolved shear stresses to dislocation densities on individual slip systems, and a generalized Kocks-Mecking model for dislocation multiplication. The DDD data strongly suggests a logarithmic dependence of flow stress on the plastic shear strain rate on each slip system, and, equivalently, an exponential dependence of the plastic shear strain rate on the resolved shear stress. Hence the proposed generalized Taylor relation subsumes the Orowan relation for plastic flow. The DDD data also calls for a correction to the Kocks-Mecking model of dislocation multiplication to account for the increase of dislocation density on slip systems with negligible plastic shear strain rate. This is accomplished by allowing the multiplication rate on each slip system to include contributions from the plastic strain rates of the two coplanar slip systems. The resulting constitutive model successfully captures the strain hardening rate dependence on the loading orientation as predicted by the DDD simulations, which is also consistent with existing experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deformation and Delamination of Dynamically Bulging Bilayered Films

We employ finite element simulations to study the dynamic bulging of films under shock tube conditions and compare these responses with quasi-static bulging. Both scenarios subject the film to time-varying pressure. Two film types are considered: uniform material films and bilayers with a relatively soft interior layer and a hard exterior layer. The interface between the materials in the bilayer cases is modeled as a cohesive zone that follows a bilinear normal traction-separation cohesive law. Additionally, bilayer films with preexisting cracks in the interior layer are analyzed. Parametric studies are also conducted to explore the effects of loading rate and specimen size. The mechanical response of all materials is assumed to be strain-rate independent, focusing solely on the inertial effects in the response of the films under dynamic loading, and the absence of such effects under quasi-static loading. The results indicate that quasi-static bulging exhibits a predictable plate bending like deformation. In contrast, depending on the loading rate and specimen size, dynamic bulging can give rise to elastic waves, and a sequence of deformation processes including initial uniform acceleration, followed by bending-like deformation, and finally through-thickness biaxial stretching. Dynamic loading also leads to larger bulge formation and greater deformation than quasi-static loading. A key finding is the role of preexisting cracks in driving delamination under dynamic loading, which is absent in quasi-static loading. This work enhances our understanding of dynamic bulging and motivates further research on characterizing delamination in layered films using dynamic bulge tests.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A concurrent irradiation-mechanics multiscale coupling model

The vast majority of our current knowledge regarding the basic mechanisms controlling irradiation effect on mechanical properties is based almost entirely on the results of post-irradiation experiments or theoretical models. However, the concurrent effects of irradiation, mechanical stress, and thermal damage on the failure phenomena of materials and components remain largely unexplored due to its internal multiscale-multiphysics coupling nature. Here, we present here a concurrent irradiation-mechanics multiscale coupling model. The concurrent evolutions of nanoscale irradiation defect clusters, microscale dislocation configurations, and mechanical responses are well captured based on coupling cluster dynamics, discrete dislocation dynamics, and the finite element methods using an effective time marching scheme. Model predictions of defect densities and size are in general agreement with experimental observations. Irradiation hardening is shown to take place also in samples undergoing concurrent irradiation-mechanical loading, similar to samples tested post-irradiation. However, the occurrence of plastic flow localization and dislocation channel formation is not accompanied with apparent yield drop (softening) under concurrent irradiation-mechanical loading conditions, which is different from the post-irradiation case.

42 ENGINEERING↗
DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Helical Fins for Concentrated Solar Receivers: Design Optimization and Entropy Analysis

Concentrated solar power (CSP) with thermal energy storage (TES) has the potential to achieve grid parity. This can be realized by operating CSP systems at temperatures above 700 °C with high-efficiency sCO 2 power cycles. However, operating CSP systems at such temperatures poses several challenges, among which the design of solar receivers to accommodate increased thermal loads is critical. To this end, this work explores and optimizes various swirl-inducing internal fin designs for solar receiver tubes. These fin designs not only improve the thermal performance of receiver tubes but also levelize temperature unevenness caused by non-uniform thermal loading. In this work, the geometric parameters of the fin designs are optimized to maximize the Nusselt number with a constraint on the friction factor. This optimization, however, is computationally intensive, requiring hundreds of simulation calls to computational fluid dynamics (CFD) models. Here, to circumvent this problem, surrogate models are used to approximate the simulation outputs needed during the optimization. In addition, this study also examines the fin designs from an entropy generation perspective. To this end, the entropy contributions from thermal and viscous effects are quantitatively compared while varying the operational Reynolds number.

14 SOLAR ENERGY↗
DOE OSTI2024

A Dynamic Model of Refrigerator with Thermal Energy Storage for Demand Flexibility

A refrigerator is an essential appliance for both commercial and residential applications. Thermal (Cold) energy storage can effectively reduce the carbon emission from refrigerators. This preliminary study investigated the demand flexibility of a domestic refrigerator using thermal (cold) energy storage through dynamic modeling and controls. The models include both fresh and frozen compartments, a vapor compression cycle loop, and thermal storage panels. Results show that thermal energy storage can effectively achieve temperature control of the compartments within the 4-hour time window. The load reduction and cost savings (46.3%) were achieved through the TES using TOU rates.

Gao, Zhiming↗
DOE OSTI2024

Measurement of Wind Loading on Heliostats at the Crescent Dunes Power Plant: An Overview

The cost of solar collectors constitutes almost one third of the total cost of a CSP plant. One of the ongoing challenges in the design of these collectors is wind loading on mirrors, support structures, and drives. A particular challenge is dynamic wind loading, caused by the turbulent wind flow. To date, the design of solar collector structures has relied on wind tunnel experiments and numerical simulations that do not entirely capture the dynamic effects observed at scale. The CSP industry has shown increased interest in validating the idealized assumptions with measurements obtained in operational settings to improve wind load assumptions and increase reliability and cost-efficiency of the collector design. Further, performance models need realistic assumptions about wind loading and its impact on optical performance. In a parabolic trough field campaign, NREL successfully collected a wealth of long-term, high-resolution wind and loads data [2], that are publicly available and that can be used for the above-mentioned purposes. Heliostats are impacted differently by wind than parabolic troughs, due to their different shape, size, and field layout. To study the impact of wind and turbulence on heliostats, we initiated another field campaign in an operational power- tower plant, Crescent Dunes, in Nevada, USA. In this work, we present an overview of the measurements, first results, and potential implications of wind driven loads on Heliostats.

concentrated solar power↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗