Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variable Grid Method”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗
DOE OSTI2020

Hierarchical, Grid-Aware, and Economically Optimal Coordination of Distributed Energy Resources in Realistic Distribution Systems

Renewable portfolio standards are targeting high levels of variable solar photovoltaics (PV) in electric distribution systems, which makes reliability more challenging to maintain for distribution system operators (DSOs). Distributed energy resources (DERs), including smart, connected appliances and PV inverters, represent responsive grid resources that can provide flexibility to support the DSO in actively managing their networks to facilitate reliability under extreme levels of solar PV. This flexibility can also be used to optimize system operations with respect to economic signals from wholesale energy and ancillary service markets. Here, we present a novel hierarchical scheme that actively controls behind-the-meter DERs to reliably manage each unbalanced distribution feeder and exploits the available flexibility to ensure reliable operation and economically optimizes the entire distribution network. Each layer of the scheme employs advanced optimization methods at different timescales to ensure that the system operates within both grid and device limits. The hierarchy is validated in a large-scale realistic simulation based on data from the industry. Simulation results show that coordination of flexibility improves both system reliability and economics, and enables greater penetration of solar PV. Discussion is also provided on the practical viability of the required communications and controls to implement the presented scheme within a large DSO.

distribution network↗
DOE OSTI2020

Embedding physics domain knowledge into a Bayesian network enables layer-by-layer process innovation for photovoltaics

Process optimization of photovoltaic devices is a time-intensive, trial-and-error endeavor, which lacks full transparency of the underlying physics and relies on user-imposed constraints that may or may not lead to a global optimum. Herein, we demonstrate that embedding physics domain knowledge into a Bayesian network enables an optimization approach for gallium arsenide (GaAs) solar cells that identifies the root cause(s) of underperformance with layer-by-layer resolution and reveals alternative optimal process windows beyond traditional black-box optimization. Our Bayesian network approach links a key GaAs process variable (growth temperature) to material descriptors (bulk and interface properties, e.g., bulk lifetime, doping, and surface recombination) and device performance parameters (e.g., cell efficiency). For this purpose, we combine a Bayesian inference framework with a neural network surrogate device-physics model that is 100× faster than numerical solvers. With the trained surrogate model and only a small number of experimental samples, our approach reduces significantly the time-consuming intervention and characterization required by the experimentalist. As a demonstration of our method, in only five metal organic chemical vapor depositions, we identify a superior growth temperature profile for the window, bulk, and back surface field layer of a GaAs solar cell, without any secondary measurements, and demonstrate a 6.5% relative AM1.5G efficiency improvement above traditional grid search methods.

14 SOLAR ENERGY↗
DOE OSTI2024

ARM shortwave spectrometers to study the clear-cloud transition zone and mixing processes

The proposed research is a collaborative effort between NASA/Goddard Space Flight Center and the Hebrew University of Jerusalem. While the NASA team focused on analyzing ground-based hyper-spectral radiance observations to understand cloud edge properties and their connection to mixing processes, the Hebrew University team tackled the problem through cloud modeling activities. By approximating the shortwave spectra in the cloud-clear transition zone as a linear combination of purely clear and purely cloudy spectra we can characterize the variations of cloud optical thickness and cloud droplet effective radius in the transition zone. When applying this method to the measurements of a ground-based shortwave spectroradiometer at the ARM’s SGP site, representing continental conditions, and MAGIC field campaign between Log Angeles, California and Honolulu, Hawaii, representing maritime scenarios, we found that cloud optical depth consistently decreases in both cases, but droplet size decreases much more substantially for the continental regime, suggesting different mixing processes for the continental and maritime conditions. The investigation and measurements of radiation clouds were coupled with a unique cloud modeling. A novel spectral bin microphysics was developed and implemented to the System of Atmospheric Modeling (SAM). In order to resolve cloud transition zones with high spatial gradients of microphysical variables a unique high resolution (10 m) was used in simulations. The model calculates droplet size distributions in each grid point. The model output was transferred to the NASA/GSFC team for utilization in radiative calculations and testing of both radiative algorithm and model representation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2022

The WRF-Solar Ensemble Prediction System: Development, Test, and Validation

Providing reliable probabilistic solar radiation information is needed to improve management of the uncertainty and variability of solar generation. Thus, guidance on how to develop skillful and accurate ensemble forecasts is essential and it will ultimately contribute to integration of high amounts of solar energy on the grid. A team from the National Renewable Energy Laboratory and the National Center for Atmospheric Research had been collaborating to develop the WRF-Solar ensemble prediction system (WRF-Solar EPS) in the past three years to produce probabilistic solar irradiance forecasts and better predict solar energy by quantifying forecast uncertainty. The WRF-Solar EPS basically generates ensemble members for solar irradiance based on stochastic perturbations to provide intraday and day-ahead probabilistic forecasts. This study will present main research steps in developing the WRF-Solar EPS including: (a) tangent linear analysis for identifying key input variables of six WRF-Solar modules significantly related to predicting of cloud and solar irradiance, (b) combining stochastic perturbation technique with the WRF-Solar model, and (c) ensemble calibration method to decrease error and uncertainty of ensemble-based solar forecasts. The capability of WRF-Solar EPS is now updated to the most recent version of standard WRF model. This presentation will summarize comprehensive results from the evaluation of forecasts against the National Solar Radiation Data Base as well as ground-measured observations. Moreover, we will introduce the user's guide for WRF-Solar EPS (e.g., parameters to configure stochastic perturbations) and future extension of this research.

day-ahead forecast↗
DOE OSTI2025

In situ Detection of Plasma Induced Surface Interaction based on Deep Learning based Visual Diagnostics (Technical Report)

It is characteristic for many plasma devices to undergo plasma-material interaction leading to surface erosion. These processes, often not easily detectable, lead to changes in device performance and lifespan. State-of-the-art lifetime tests and wear experiments require over 1000s hours. A self-consistent model for accurately predicting the erosion's effects is not available. In situ detection of these processes is not a trivial task since the surface variations at the early stages have a micron scale. Such limitations not only restrict testing and prediction capabilities but also slow the development of new thrusters and limit mission duration. To address these challenges, an in-situ diagnostic for real-time erosion assessment has been developed, aiming to expedite lifetime testing and broaden experimental campaigns. Several works were dedicated to real-time and in situ monitoring of material erosion during plasma exposure using laser holography, microscopy, and with telemicroscopes. However, the applicability of these approaches is limited due to complexity, cost and less flexibility as they often require placing diagnostic equipment inside the vacuum chamber. In collaboration with Princeton Collaborative Research Facility (PCRF), Princeton Plasma Physics Laboratory (PPPL), a new diagnostic approach is developed, where geometry modifications to the ceramic channel walls were introduced that would result in accelerated channel erosion. We employed Long-distance microscope (LDM) imagery, combined with Deep-Learning based Shape from focus or depth from focus (DFF or SFF) approach, that provides an accessible and cost-effective solution. LDM employs focus variation techniques to continuously capture multiple images of the target object at distinct focal planes. DFF, an optical focus variation method, generates a 3D topographical surface depth map from a sequence of variably focused images. Combined with the developed diagnostic, this approach offers a controllable means to study erosion under accelerated conditions. In this work, we develop Neural Network-based DFF algorithm applicable for LDM data to quantitatively evaluate plasma induced surface modification from LDM data. Next, we develop Deep Learning-based super-resolution depth map image reconstruction technique to increase the resolution of depth maps obtained from DFF algorithm to improve the accuracy of erosion measurements. Thirdly, we develop several image processing techniques to remove noise and improve the quality of depth map image. Here we report the results of initial tests for this approach. An experimental setup designed and built in PPPL was employed that consists of a 3-cm gridded ion source that produces a neutralized argon beam with energies up to 600 eV. A hexagonal boron nitride (h-BN) ceramic target, designed based on computational predictions, was used. Tests were conducted to reconstruct the complex geometry of the target under the lighting conditions of the operated ion source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Myriad World Baseline: Global Geodemographic Estimates

The LandScan Myriad World Baseline (MWB) method produces global, residential (nighttime/home-location) gridded geodemographic estimates based on 5-year age/gender cohorts—at 30-arcsecond (≈1 km) resolution. MWB is designed to fill gaps where detailed, georeferenced survey data (e.g., Demographic and Health Surveys (DHS)) are missing or outdated, and to provide a baseline that can support human security analysis, including consequence assessment, “patterns of life” modeling, and scenario-based population futures. MWB’s workflow spatializes household-level age/gender characteristics from the GLOPOP-S dataset by conflating household and gridded expected relative wealth adapted from Global Gridded Relative Deprivation Index (GRDI), then adjusts them to a target year of interest. Age/gender estimates are then applied to harmonize lowest-administrative-level statistics with LandScan residential counts, yielding final geodemographic estimates. Two validation case studies are presented: Ghana (2021) and Tokyo/Kanagawa, Japan (2020), illustrating spatial variability in demographic cohorts and comparing MWB outputs to official gridded statistics. Results show close overall alignment relative to validation criteria including population pyramids and age-dependency ratios.

Tuccillo, Joe [ORNL] (ORCID:0000000259300943)↗
DOE OSTI2023

Sensitivity of the pseudo-global warming method under flood conditions: a case study from the northeastern US

Abstract. Intensified extreme precipitation and corresponding floods are the most relevant consequences of climate change over the northeastern US (NEUS). To evaluate the impacts of climate change or certain climate perturbations on future extreme weather events which are dynamically similar to historic analogs, the pseudo-global warming (PGW) method has been frequently employed; however, this method lacks precise definition and guidelines, thus limiting its application. More specifically, three key questions related to the application of the PGW method remain unanswered: at what spatial scale should climate perturbations be applied? Among the different meteorological variables available, which ones should be perturbed? And will PGW projections vary significantly when different perturbations are applied? To address these questions, we examine the sensitivity and robustness of conclusions drawn from the PGW method over the NEUS by conducting multiple PGW experiments with varied perturbation spatial scales and choice of perturbed meteorological variables. The results show that the projections of precipitation and other essential variables at the regional mean scale are consistent across the PGW simulations, with a relative difference of much less than 10 %; however, different perturbation modifications can cause significant displacements of the storm events being simulated. Several previously assumed advantages of modifying only the temperature at regional mean scale, such as the preservation of geostrophic balance, do not appear to hold. Also, for these experiments, we find the regional mean perturbation produces a positive precipitation bias because it ignores the land–ocean warming contrast, which is a robust regional response to global warming. Overall, PGW experiments with perturbations from temperature or the combination of temperature and wind at the grid point scale are both recommended, depending on the research questions. The first approach can isolate the spatially dependent thermodynamic impact, and the latter incorporates both the thermodynamic and dynamic impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Gradient Based Subgrid-Scale Parameterization for Ocean Mesoscale Eddies

Mesoscale eddies play an important role in transport of heat and biogeochemical tracers in the global ocean circulation. Resolving these energetic eddies, however, is challenging in ocean general circulation models (OGCM) because it requires a horizontal grid spacing of ≲1/8° that is computationally expensive. As a result, we are required to parameterize mesoscale eddy effects on large-scale ocean flows. In this work, we introduce a new subgrid-scale (SGS) model that is developed based on a Taylor series expansion of resolved variables to parameterize subgrid mesoscale eddy transports and momentum fluxes in OGCM. We have performed an a priori study to evaluate the performance of our new gradient model using high-resolution ocean simulations. Our results show that the gradient model well predicts the actual SGS thickness fluxes in the zonal and meridional directions in coarse-resolution simulations with the grid spacing ≳1/4°. The unresolved kinetic energy at the ocean surface is also skillfully estimated. More importantly, unlike current mesoscale eddy parameterizations, which are mainly developed based on an assumption of flat bottom ocean, our new SGS model can capture the structure of unresolved standing meanders at the ocean surface. We have also developed a dynamic procedure for setting in non-dimensional parameters in our new parameterization through a non-ad hoc and tuning-free method. Overall, this work suggests that implementing the gradient model in OGCM can improve the model accuracy with an affordable computational cost in eddy-permitting and non-eddying simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Improving Prediction of Surface Solar Irradiance Variability by Integrating Observed Cloud Characteristics and Machine Learning

A 5-year, 1-minute resolution observational dataset of clouds and solar radiation was produced that includes two metrics of the variability in surface solar irradiance due to cloud type and fractional sky cover. Multiple regression models were trained to fit observations of surface solar irradiance variability from those two cloud property predictors. We found that ensemble tree-based methods, Random Forest and Gradient Boosting Machine, have the least overfitting issues and showed the best performance with an R2 of 0.42. While the observational data trained in this study was only from one site, the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site in Oklahoma, initial comparisons of the seasonality of the statistics suggest that these results are relatively weather regime independent; the generality of such a finding across sites will be tested in future work. The observational data and developed machine learning model are being used to create a numerical weather prediction model parameterization to enable day-ahead solar variability prediction in a computationally efficient way. This is a first step towards creating a new paradigm of predicting day-ahead variability with the potential to provide a new tool to improve grid operation, planning, and resilience.

Riihimaki, Laura↗
DOE OSTI2024

Comparing Gravity Waves in a Kilometer‐Scale Run of the IFS to AIRS Satellite Observations and ERA5

Abstract Atmospheric gravity waves (GWs) impact the circulation and variability of the atmosphere. Sub‐grid scale GWs, which are too small to be resolved, are parameterized in weather and climate models. However, some models are now available at resolutions at which these waves become resolved and it is important to test whether these models do this correctly. In this study, a GW resolving run of the European Center for Medium‐Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS), run with a 1.4 km average grid spacing (TCo7999 resolution), is compared to observations from the Atmospheric Infrared Sounder (AIRS) instrument, on NASA's Aqua satellite, to test how well the model resolves GWs that AIRS can observe. In this analysis, nighttime data are used from the first 10 days of November 2018 over part of Asia and surrounding regions. The IFS run is resampled with AIRS's observational filter using two different methods for comparison. The ECMWF ERA5 reanalysis is also resampled as AIRS, to allow for comparison of how the high resolution IFS run resolves GWs compared to a lower resolution model that uses GW drag parametrizations. Wave properties are found in AIRS and the resampled models using a multi‐dimensional S‐Transform method. Orographic GWs can be seen in similar locations at similar times in all three data sets. However, wave amplitudes and momentum fluxes in the resampled IFS run are found to be significantly lower than in the observations. This could be a result of horizontal and vertical wavelengths in the IFS run being underestimated.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Impact of FERC Order 2222 on DER Participation Rules in US Electricity Markets

Electricity markets in the bulk grid are beginning to implement market mechanisms that support the procurement of flexible capabilities from wide range of technologies, including distributed energy resources (DERs). The flexibility of these resources will help counterbalance supply uncertainties from large-scale integration of variable renewable generation. To encourage development of distributed and aggregated market participants, FERC Order 2222 was issued in September 2020 to require each Independent System Operator (ISO) in the US to implement rules that enable broader participation from aggregations of DERs in the bulk market. The following paper first describes the generic design of ISO markets before introducing the new market participation rules that ISOs have proposed for compliance with Order 2222. The paper then describes how software performance issues may continue to affect the eligibility requirements and offer structures for DER aggregations participating in ISOs, noting that continued research on computational methods may help reduce burdens for DER integration. The prospects for transmission and distribution system coordination is second major issue discussed, which will require minor changes to existing processes in the short term. In the longer term, there is more opportunity for more wide-ranging reforms, such as the development of a Distribution System Operator (DSO) framework. Newly proposed market rules may affect how Transactive Energy Systems (TES) will help facilitate efficient formation of DER aggregations and operation of the individual DERs within an aggregation. Within the TES context, the challenge is to fully understand how resource eligibility and operational and planning coordination methods will affect the design and implementation of TES.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Uncertainty in Atlantic Multidecadal Oscillation derived from different observed datasets and their possible causes

As a leading mode of sea surface temperature (SST) variability over the North Atlantic in both observations and model simulations, the Atlantic Multidecadal Oscillation (AMO) can have a substantial influence on regional and global climate. By using Low-Frequency Component Analysis, we explore the uncertainties of the resulting AMO indices and the corresponding spatial patterns derived from three observational SST datasets. We found that the known coherent spatial pattern of the AMO at the basin scale over the North Atlantic appears in two out of the three datasets. Further analysis indicates that both the warming trend and the different techniques used to construct these observed gridded SSTs contribute to the AMO’s spatial coherence over the North Atlantic, especially during periods of sparse data sampling. The SST in the Extended Reconstructed SST dataset version 5 (ERSSTv5), changes from being systematically below the other datasets during the dense sampling periods on either side of the Second World War (WWII), to systematically above the other datasets during WWII, thereby introducing an artificial 10–20-year variability that affects the AMO’s spatial coherence. This coherence in the AMO’s spatial pattern is also affected by bias adjustment in ERSSTv5 at relative cool (i.e., non-summer) seasons, and by the heterogeneous North Atlantic warming pattern. The different AMO patterns can induce the different effects of wind, surface heat fluxes, and then drive ocean circulation and its heat transport convergence, especially for some seasons. For AMO indices, both the different detrending methods and different observational data result in uncertainty for the period 1935–1950. Such SST uncertainty is important to detect the relative role of the atmosphere and ocean in shaping the AMO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A non-hierarchical correlation discrete variable representation

The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. Here, the non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH 3 demonstrate the accuracy of the non-hierarchical CDVR approach.

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

Multi-Timescale Integrated Dynamics and Scheduling for Solar (MIDAS-Solar) (Final Technical Report)

Solar photovoltaic (PV) installations have experienced unprecedented growth in the United States. PV will become not only an energy producer but also a necessary provider of ancillary services at multiple timescales. Conventional methods to simulate power systems operations - such as long-term production simulation (which typically considers schedules from hours to minutes by using an optimization framework) and short-term transient studies (which simulate dynamics from seconds to sub-seconds using state variables and differential equations) - are not sufficient for studying the multiple-timescale variation of solar generation and its impact on system reliability. Long-term system economics and short-term system dynamics are highly coupled, particularly when the penetration level of renewable generation is extremely high, because the uncertainty and variability of solar generation will impact both power system steady-state and dynamic performance. This project helps meet and exceed the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Solar Energy Technologies Office goal of systems integration by directly addressing this stability and reliability challenge for power grid planning and operation. We have developed a temporally comprehensive, closed-loop simulation model, named Multi-timescale Integrated Dynamics and Scheduling (MIDAS), that seamlessly simulates power system operations from economic scheduling (day-ahead to hours) to dynamic response analysis (seconds to sub-seconds). For schedules with very high levels of inverter-based resources (IBRs), up to and including 100%, the stability of grid controls has been evaluated through electromagnetic transient (EMT) simulations and power-hardware-in-the-loop (PHIL) simulations of key transient events at key schedule points. Specifically, MIDAS provides: 1) a closed-loop simulation framework for simulating timescales from economic scheduling to dynamic stability analysis; 2) machine learning-based stability assessment; 3) EMT modeling and analysis for large-scale power systems; 4) MIDAS PHIL test bed. We worked with Hawaii Electric Companies to apply the MIDAS study framework to a Maui grid study. The entire island's transmission system was modeled in detail - from a yearly scheduling model, to a second-level frequency dynamic model, down to a sub-second-scale EMT model to address critical stability issues. The project demonstrated how MIDAS can help system planners and operators assess system reliability and stability while the power grid is marching toward a high-renewable, high-IBR future. In this Maui grid study, we found that 100% instantaneous IBR operation is achievable in EMT simulation and PHIL testing, and grid planners and operators might need new analysis/simulation tools to assess grid reliability and stability in the scheduling stage. MIDAS will bring Maui and other systems closer to 100% clean and stable energy futures. (In this study, we examined transient stability. Other topics necessary for 100% IBR operation, such as protection and resource adequacy, were not examined.)

100% Renewables↗
DOE OSTI2020

A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity

In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.

42 ENGINEERING↗