Search NASA⌕ Search

SEARCH · Search NASA

Results for “modelling studies”

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

Impacts of Spatial Resolution in a High-Fidelity Capacity Expansion Model: An ERCOT Case Study

Capacity expansion models are important tools in examining the evolution of the electric power sector. Embedded in these tools are many modeling choices with consequential impacts on computational burden and associated analysis. In this study, we adjust the spatial resolution of the Regional Energy Deployment System (ReEDS) to understand the implications of higher-fidelity modeling on energy system projections and model solve times. The native ReEDS regions capture the contiguous United States in 134 balancing areas whereas the regions in the higher-resolution version are defined by over 3,000 U.S. counties. Using both resolutions, we conduct a case study of the Texas Interconnection (The Electric Reliability Council of Texas [ERCOT]) to explore differences in model projections and to inform appropriate applications of high spatial resolution in a large-scale, applied capacity expansion model.

county↗
DOE OSTI2023

Moiré Kanamori-Hubbard model for transition metal dichalcogenide homobilayers

Ab initio and continuum model studies predicted that the Γ valley transition metal dichalcogenide (TMD) homobilayers could simulate the conventional multiorbital Hubbard model on the moiré honeycomb lattice. Here, in this work, we perform the Wannierization starting from the continuum model and show that a more general moiré Kanamori-Hubbard model emerges, beyond the extensively studied standard multiorbital Hubbard model, which can be used to investigate the many-body physics in the Γ valley TMD homobilayers. Using the unrestricted Hartree-Fock and Lanczos techniques, we study these half-filled multiorbital moiré bands. By constructing the phase diagrams we predict the presence of an antiferromagnetic state and in addition we found unexpected and dominant states, such as a S = 1 ferromagnetic insulator and a charge density wave state. Our theoretical predictions made using this model can be tested in future experiments on the Γ valley TMD homobilayers.

2-dimensional systems↗
DOE OSTI2025

Quantifying Uncertainties in Modeling Wind Resource Data from Different PBL Schemes in the WRF Model: A Case Study Over the Puerto Rico Region

This study examines the modeling uncertainty of wind resource data stemming from the use of various planetary boundary layer (PBL) parameterizations available in the Weather Research and Forecasting (WRF) model. WRF-based wind simulations spanning 20 years at 3-km resolution using 11 different PBL schemes are used to objectively investigate the uncertainty in modeling wind speed for land-based wind (LBW) and offshore wind (OSW) locations in Puerto Rico. The uncertainty in the wind modeling for the 20-year dataset is quantified using the spread index (SI) and standard deviation (SD). For virtual LBW and OSW sites, the SI and SD values are analyzed as calculated across various spatial and temporal scales. Because the PBL's atmospheric stability conditions can be characterized into two dominant categories, the study focuses on analyzing the SI and SD for daytime (mainly unstable PBL conditions) and nighttime (mainly stable PBL conditions). For wind shear (10 m-200 m) at the OSW and LBW sites, WRF-based numerical experiments indicate the following SI (or SD) ranges: 39%-94% (0.74 m/s-1.44 m/s) during the daytime for OSW, 50%-75% (0.68 m/s-1.19 m/s) during the daytime for LBW, 37%-60% (0.73 m/s-1.12 m/s) during the nighttime for OSW, and 57%-143 % (0.65 m/s-1.43 m/s) during the nighttime for LBW. While a high SI is observed when modeling LBW during the nighttime, there are notable modeling uncertainties during the daytime on the leeward side of the orographic barriers for Puerto Rico.

17 WIND ENERGY↗
DOE OSTI2022

Surrogate Modeling of Nonlinear Dynamic Systems: A Comparative Study

Surrogate models play a vital role in overcoming the computational challenge in designing and analyzing nonlinear dynamic systems, especially in the presence of uncertainty. This paper presents a comparative study of different surrogate modeling techniques for nonlinear dynamic systems. Four surrogate modeling methods, namely, Gaussian process (GP) regression, a long short-term memory (LSTM) network, a convolutional neural network (CNN) with LSTM (CNN-LSTM), and a CNN with bidirectional LSTM (CNN-BLSTM), are studied and compared. All these model types can predict the future behavior of dynamic systems over long periods based on training data from relatively short periods. The multi-dimensional inputs of surrogate models are organized in a nonlinear autoregressive exogenous model (NARX) scheme to enable recursive prediction over long periods, where current predictions replace inputs from the previous time window. Three numerical examples, including one mathematical example and two nonlinear engineering analysis models, are used to compare the performance of the four surrogate modeling techniques. The results show that the GP-NARX surrogate model tends to have more stable performance than the other three deep learning (DL)-based methods for the three particular examples studied. The tuning effort of GP-NARX is also much lower than its deep learning-based counterparts.

42 ENGINEERING↗
DOE OSTI2024

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗
DOE OSTI2020

Experimental validation of solid oxide fuel cell polarization modeling: An LSM-YSZ/YSZ/Ni-YSZ case study

Polarization modeling based on the measured I-V curves is a useful tool for analyzing SOFC performance and diagnosing performance limitations. Using the polarization model, the overpotential of a cell is usually separated into different polarization contributions, including area specific ohmic polarization, activation polarization, and anodic and cathodic concentration polarizations. In this study, we aim to experimentally validate the accuracy of the polarization model. As a case study, anode-supported cells with LSM-YSZ/YSZ/Ni-YSZ configuration are used. By curve-fitting the I-V curves of a cell tested under different cathode oxygen partial pressures into the polarization model, the contributions of different polarization losses are quantified. To validate the polarization model, independent experiments/studies including electrochemical impedance spectroscopy, cathode kinetics study, ex-situ diffusivity measurement, and tortuosity measurement using 3D reconstructed anode are conducted. The results from polarization modeling and from validation studies are in good agreement, thereby validating the polarization model for SOFC performance analysis.

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

CFD modeling and simulation for corrosive wear of refractory in molten slag

The development of high-wear resistant refractories having minimal production costs is facilitated by characterizing the wear mechanisms associated with their corrosive wear. Static cup testing is a commonly used method for comparing the corrosion resistance performance of two or more refractory materials. Although the static cup test conditions are not as severe as dynamic tests, this study shows that the thermal gradient present within the system during heating and cooling stages serves to generate movement of the slag leading to mechanical wear. The thermal gradient within the refractory, and between the slag and the refractory, occurs during the ramp stage of the test and lasts until the soaking stage is reached bringing the system to a thermal equilibrium. Using computational fluid dynamics (CFD) capabilities embedded within ANSYS software, this study modelled and quantified the convection currents within the slag and associated shear stresses generated on the refractory walls due to the thermal gradient. A traditional ladle furnace was employed as a case study to verify the results of the studied CFD model. The corrosion rate of the refractory lining was found to depend on the mass transfer coefficient of the refractory dissolution into the slag, and a velocity term which governs the extent of corrosion at any given location. This velocity term is a function of slag viscosity, as well as the concentration gradient and/or temperature gradient at the triple points. In this study, wall shear stress was used as a reliable proxy for identifying high-velocity regions prone to excessive corrosive wear. Elevated wall shear stress near the slag/air and slag/molten steel interfaces align with observed corrosion grooves, which reflects the intensified corrosive wear at these locations.

Ramteke, Rajat Rajat Durgesh Ramteke [University o↗
DOE OSTI2024

Battery Degradation Modeling in Hybrid Power Plants: An Island System Unit Commitment Study: Preprint

As hybrid power plants (HPPs), such as photovoltaic (PV) and battery combinations, become increasingly important in power systems with high renewable energy penetration to address PV variability and ensure grid stability. This paper focuses on the urgent need to model the coordination between PV and battery systems in HPPs while accounting for battery degradation. We present a generation scheduling model that explicitly incorporates PV-battery hybridization in the unit commitment problem. Moreover, the cost function of the HPP scheduling problem endogenously considers battery degradation with adjustable weights to strike a balance between minimizing production costs and prolonging battery life, particularly when providing energy arbitrage and ancillary services. Using a realistic island system simulation, we demonstrate that accounting for battery degradation in the scheduling problem can significantly extend battery life with only minor additional production costs.

battery degradation↗
DOE OSTI2024

Fine-scale landscape characteristics, vegetation composition, and snowmelt timing control phenological heterogeneity across low-Arctic tundra landscapes in Western Alaska

The Arctic is warming at over twice the rate of the rest of the Earth, resulting in significant changes in vegetation seasonality that regulates annual carbon, water, and energy fluxes. However, a crucial knowledge gap exists regarding the intricate interplay among climate, permafrost, and vegetation that generates high phenology variability across extensive tundra landscapes. This oversight has led to significant discrepancies in phenological patterns observed across warming experiments, long-term ecological observations, and satellite and modeling studies, undermining our ability to understand and forecast plant responses to climate change in the Arctic. To address this problem, we assessed plant phenology across three low-Arctic tundra landscapes on the Seward Peninsula, Alaska, using a combination of in-situ phenocam observations and high-resolution PlanetScope CubeSat data. We examined the patterns and drivers of phenological diversity across the landscape by (1) quantifying phenological diversity among dominant plant function types (PFTs) and (2) modeling the interrelations between plant phenology and fine-scale landscape features, such as topography, snowmelt, and vegetation. Our findings reveal that both spring and fall phenology varied significantly across Arctic PFTs, accounting for about 25%–44% and 34%–59% of the landscape-scale variation in the start of spring [SOS] and start of fall [SOF], respectively. Deciduous tall shrubs (e.g. alder and willow) had a later SOS (~7 d behind the mean of other PFTs), but completed leaf expansion (within 2 weeks) considerably faster compared to other PFTs. We modeled the landscape-scale variation in SOS and SOF using Random Forest, which showed that plant phenology can be accurately captured by a suite of variables related to vegetation composition, topographic characteristics, and snowmelt timing (variance explained: 53%–68% for SOS and 59%–82% for SOF). Notably, snowmelt timing was a crucial determinant of SOS, a factor often neglected in most spring phenology models. Our study highlights the impact of fine-scale vegetation composition, snow seasonality, and landscape features on tundra phenological heterogeneity. Improved understanding of such considerable intra-site phenological variability and associated proximate controls across extensive Arctic landscapes offers critical insights for representation of tundra phenology in process models and associated impact assessments with climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Changes in dislocation punching behavior due to hydrogen-seeded helium bubble growth in tungsten

The accumulation of gas atoms in tungsten is a topic of long-standing interest to the plasma-facing materials community due the metal's use as a divertor material in some tokamak fusion reactors. The nucleation and growth of He/H gas bubbles (along with their isotopes) can result from impinging fluxes of these gases which give rise to damage at the W divertor surface. The inclusion of He or H in W has been studied extensively by the community, finding that He bubbles modify the surface through periodic dislocation punching and bursting mechanisms while H bubbles impact the metal through plastic-strain induced material failure. However, the mechanisms which are present during the combined flux of both He and H is not well-studied atomistically. Motivated by this, an atomistic modeling study is conducted using molecular dynamics to assess the behavior of mixed concentration He:H bubbles in W. Here we find that the introduction of H into a growing He bubble results in a dramatic change in the nature and presence of dislocation loops which are typically generated via dislocation punching in over-pressurized He bubbles. Most notably, at high H concentrations, there is a switchover in energetic favorability from glissile 1/2<111> dislocations to sessile <100> dislocations. This thermodynamic crossover could imply a significant reduction in W surface morphology changes than with pure He bubbles and, additionally, we show this to have implications on the trapping of H in the bubbles and their associated dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Heterogenous electromediated depolymerization of highly crystalline polyoxymethylene

Abstract Post-consumer plastic waste in the environment has driven the scientific community to develop deconstruction methods that yield valued substances from these synthetic macromolecules. Electrocatalysis is a well-established method for achieving challenging transformations in small molecule synthesis. Here we present the first electro-chemical depolymerization of polyoxymethylene—a highly crystalline engineering thermoplastic (Delrin®)—into its repolymerizable monomer, formaldehyde/1,3,5-trioxane, under ambient conditions. We investigate this electrochemical deconstruction by employing solvent screening, cyclic voltammetry, divided cell studies, electrolysis with redox mediators, small molecule model studies, and control experiments. Our findings determine that the reaction proceeds via a heterogeneous electro-mediated acid depolymerization mechanism. The bifunctional role of the co-solvent 1,1,1,3,3,3-hexafluoro-2-propanol (HFIP) is also revealed. This study demonstrates the potential of electromediated depolymerization serving as an important role in sustainable chemistry by merging the concepts of renewable energy and circular plastic economy.

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

Deploying a Model Predictive Traffic Signal Control Algorithm - A Field Deployment Experiment Case Study

This paper presents a field deployment experiment of a real-time traffic signal control algorithm. We implemented the model predictive control (MPC) algorithm based on the virtual phase-link (VPL) model. We selected the deployment locations and times based on an energy saving potential concept. We developed a set of experiment systems, which included sensing, processing, and actuating components, to enable field deployment. We tested the systems rigorously before the experiment days. We reported the key procedures on the experiment days, including the steps taken, the real-time control procedure, and the monitoring of the experiment. We evaluated the impact of the deployment by looking at the changes in delay and energy consumption.

deployment↗
DOE OSTI2021

Pairing correlations in the cuprates: A numerical study of the three-band Hubbard model

We study the three-band Hubbard model for the copper oxide plane of the high-temperature superconducting cuprates using determinant quantum Monte Carlo and the dynamical cluster approximation (DCA) and provide a comprehensive view of the pairing correlations in this model using these methods. Specifically, we compute the pair-field susceptibility and study its dependence on temperature, doping, interaction strength, and charge-transfer energy. Using the DCA, we also solve the Bethe-Salpeter equation for the two-particle Green's function in the particle-particle channel to determine the transition temperature to the superconducting phase on smaller clusters. Our calculations reproduce many aspects of the cuprate phase diagram and indicate that there is an “optimal” value of the charge-transfer energy for the model where T c is largest. Finally, these results have implications for our understanding of superconductivity in both the cuprates and other doped charge-transfer insulators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Opportunities and challenges in process modeling and simulation of electrochemical systems

Electrochemical technologies have garnered intense interest in both academic research and industrial applications due to their potential to increase energy efficiency and reduce carbon footprint. However, electrochemical process fundamentals have been absent from the chemical process simulators available to millions of chemical engineers worldwide. To expedite process research and development of electrochemical technologies in the chemical industry, it is imperative that essential electrochemical process fundamentals be incorporated into process simulators to support modeling and simulation of electrochemical processes. Here, this study examines three process fundamentals key to the research and development of electrochemical processes: electrochemical reaction kinetics, electrolyte thermodynamics, and heat and mass transfer. It further illustrates application of these process fundamentals with a case study modeling an electrochemical process for the conversion of acrylonitrile to adiponitrile. The modeling example highlights the roles of applied voltage on reaction rates, and electron flow rates on the performance of electrochemical processes. It suggests that the applied voltage and electron flow rates are two unique concepts that should be included in simulators for electrochemical processes.

09 BIOMASS FUELS↗
DOE OSTI2022

Spatial Microbial Respiration Variations in the Hyporheic Zones Within the Columbia River Basin

Abstract While the hyporheic zone (HZ) accounts for a significant portion of whole stream CO 2 concentrations, HZ respiration modeling studies are lacking in quantifying their contributions to the total CO 2 at large watershed/basin scales. Quantifying the contribution of anaerobic respiration is also underappreciated. This study used a carbon‐nitrogen‐coupled river corridor model to quantify HZ aerobic and anaerobic respiration and determined the key factors controlling their spatial variability within the Columbia River Basin (CRB). The modeled respiration patterns showed high spatial variability. Among the nine sub‐basins composing the CRB, the Lower Columbia and the Willamette, which receive higher precipitation, had higher respiration. Medium‐sized rivers (fourth to sixth orders) produced the highest aerobic and anaerobic respiration among reaches of different sizes. At the basin scale, aerobic respiration is dominant, representing approximately 98.7% of the total respiration across the CRB. While most of the reaches were dominant with aerobic respiration, reaches in agricultural land showed a relatively higher anaerobic respiration (18%) ratio. A variable importance analysis showed that hyporheic exchange flux controlled most of the spatial variability of HZ respiration, dominating over other physical variables such as residence time, stream dissolved organic carbon (DOC), nitrate, and dissolved oxygen (DO). The influence of substrate concentration (DOC and DO) is larger in modeling anaerobic respiration than aerobic respiration. Future efforts will focus on improving the estimation of the HZ exchange flux and the implementation of spatially explicit parameterizations for the reactions of interest to reduce model uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Moving-Wave Implementation in WRF to Study the Impact of Surface Water Waves on the Atmospheric Boundary Layer

Abstract While numerous modeling studies have focused on the interaction of ocean surface waves with the atmospheric boundary layer, most employ idealized waves that are either monochromatic or synthetically generated from a theoretical wave spectrum, and the atmospheric solvers are typically incompressible. To study wind–wave coupling in real-world scenarios, a model that can simulate both realistic meteorological and wave conditions is necessary. In this paper we describe the implementation of a moving bottom boundary condition into the Weather Research and Forecasting Model for large-eddy simulation applications. We first describe the moving bottom boundary conditions within WRF’s pressure-based vertical coordinate system. We then validate our code with idealized test cases that have analytical solutions, including flow over a monochromatic wave with and without viscosity. Finally, we present results from turbulent flows over a moving monochromatic wave with different wave ages, and demonstrate satisfactory agreement of the wave growth rate with results from the literature. We also compare atmospheric stress and wind parameters from two physically equivalent cases. The first specifies a wind moving in the same direction as a propagating wave, while the second involves a stationary wave with the wind adjusted such that the wind relative to the wave is the same as in the first case. Results indicate that the velocity and Reynolds stress profiles for the two cases match, further validating the moving bottom implementation.

17 WIND ENERGY↗
DOE OSTI2022

Rain-Induced Stratification of the Equatorial Indian Ocean and Its Potential Feedback to the Atmosphere

Surface freshening through precipitation can act to stably stratify the upper ocean, forming a rain layer (RL). RLs inhibit subsurface vertical mixing, isolating deeper ocean layers from the atmosphere. This process has been studied using observations and idealized simulations. The present ocean modeling study builds upon this body of work by incorporating spatially resolved and realistic atmospheric forcing. Fine-scale observations of the upper ocean collected during the Dynamics of the Madden-Julian Oscillation field campaign are used to verify the General Ocean Turbulence Model (GOTM). Spatiotemporal characteristics of equatorial Indian Ocean RLs are then investigated by forcing a 2D array of GOTM columns with realistic and well-resolved output from an existing regional atmospheric simulation. RL influence on the ocean-atmosphere system is evaluated through analysis of RL-induced modification to surface fluxes and sea surface temperature (SST). This analysis demonstrates that RLs cool the ocean surface on time scales longer than the associated precipitation event. A second simulation with identical atmospheric forcing to that in the first, but with rainfall set to zero, is performed to investigate the role of rain temperature and salinity stratification in maintaining cold SST anomalies within RLs. Approximately one third, or 0.1°C, of the SST reduction within RLs can be attributed to rain effects, while the remainder is attributed to changes in atmospheric temperature and humidity. The prolonged RL-induced SST anomalies enhance SST gradients that have been shown to favor the initiation of atmospheric convection. These findings encourage continued research of RL feedbacks to the atmosphere.

54 ENVIRONMENTAL SCIENCES↗