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 109 records · Page 6

Progress on Optimizing Wind Farms and Rotor Designs Using Adjoints

Modern wind plants are increasingly tasked with multiple performance objectives. In addition to designing plants that maximize power output and minimize the levelized cost of energy (LCOE), the design and operation of wind plants is increasingly influenced by challenges regarding grid integration of variable generation renewables. This places a growing emphasis on making wind plants more controllable and predictable. WindSE is a Reynolds-averaged Navier-Stokes (RANS) model designed around analytical gradient and adjoint methods, with the ability to capture terrain-induced effects, as shown in Figure 1. The recent addition of an unsteady solver with an actuator line method (ALM) and ongoing work to enable massively parallel optimizations gives it a unique niche to explore coupled plant-level controls and design problems. This code is an open source python package built on the FEniCS framework that utilizes fast, parallel PETSc solvers to model fluid flow throughout wind-farm scale domains. Two recent studies performed using WindSE demonstrate the capability to optimize under a wide variety of flow conditions and objective functions. In the first, we present an optimization focused on modifying the layout of a wind farm with a fixed number of turbines for maximum total power output [1]. This study highlights the ability to quickly perform simulations using the steady Navier-Stokes solver combined with rotors represented as actuator disks while also stressing the importance of capturing terrain-induced effects. Gradient-based optimization using the RANS equations is viable due to the inclusion of efficiently computed adjoint derivatives. We interpret the physical results of the optimal layout and also discuss the computational cost of scaling to larger problems. In the second study, we present the capabilities of the unsteady Navier-Stokes solver, where rotor-blade profiles represented by actuator lines are optimized to enhance wake steering effects and overall power production [2]. We quantify the wind plant performance gains obtained from this type of simultaneous control co-design optimization as compared to optimizing the blade design and yaw independently. Figure 2 shows the differences between a baseline two-turbine system and an optimized system where we fine-tune the blade chord profile. Results and challenges from each study are quickly summarized and used to motivate the current development efforts within WindSE. Current and future work is focused on enabling higher-resolution studies with more degrees of freedom through parallelization of both the simulation and optimization algorithms. We present benchmarking results to show that WindSE performs well in both weak- and strong-scaling tests and further demonstrate that the optimizer obtains the same convergence rates in both shared- and distributed-memory environments. Using larger wind farms, we can study deep-array effects within an optimization context, allowing the use of objective functions that have been previously unstudied. As an example, we present ongoing work on a blockage metric which characterizes the loss of available kinetic energy due to wake effects from multiple upstream turbines.

adjoint optimization↗

Identifying snowfall elevation patterns by assimilating satellite-based snow depth retrievals

Precipitation in mountain regions is highly variable and poorly measured, posing important challenges to water resource management. Traditional methods to estimate precipitation include in-situ gauges, Doppler weather radars, satellite radars and radiometers, numerical modeling and reanalysis products. Each of these methods is unable to adequately capture complex orographic precipitation. Here, we propose a novel approach to characterize orographic snowfall over mountain regions. We use a particle batch smoother to leverage satellite information from Sentinel-1 derived snow depth retrievals and to correct various gridded precipitation products. This novel approach is tested using a simple snow model for an alpine basin located in Trentino Alto Adige, Italy. Here, we quantify the precipitation biases across the basin and found that the assimilation method (i) corrects for snowfall biases and uncertainties, (ii) leads to cumulative snowfall elevation patterns that are consistent across precipitation products, and (iii) results in overall improved basin-wide snow variables (snow depth and snow cover area) and basin streamflow estimates.

54 ENVIRONMENTAL SCIENCES↗

Characterizing uncertainties in ground “truth” of precipitation over complex terrain through high-resolution numerical modeling

Ground observation of precipitation over complex terrain is subject to large uncertainties due to inadequate sampling. This study explores a method that combines limited gauge data and a high-resolution numerical simulation to quantify the precipitation uncertainties in central Himalaya. Specifically, the Coefficient of spatial Variability (CV) of precipitation and the minimum Number of Required Stations (NRS) to obtain areal-mean precipitation ground truth values within a 0.25-degree area are investigated using fine-scale meteorological simulation at 1.5 km grid spacing. Evaluation over a densely-gauged area demonstrates comparable CV and NRS values between station observations and simulations. The simulated CV and NRS values in a larger area show a strong and positive dependence on each other and an expected positive (negative) correlation with topographic complexity (temporal scale). The proposed method sheds lights on evaluating precipitation products and holds promise for informing the layout of rain gauge networks in complex terrain.

Ouyang, Lin↗

Optimal interpolation of global dissolved oxygen: 1965–2015

Oxygen inventory of the global ocean has declined in recent decades potentially due to the warming-induced reduction in solubility as well as the circulation and biogeochemical changes associated with ocean warming and increasing stratification. Earth System Models predict continued oxygen decline for this century with profound impacts on marine ecosystem and fisheries. Observational constraint on the rate of oxygen loss is crucial for assessing the ability of models to accurately simulate these changes. There are only a few observational assessments of the global oceanic oxygen inventory reporting a range of oxygen loss. This study develops a gridded data set of dissolved oxygen for the global oceans using optimal interpolation method. The resulting gridded product includes full-depth map of dissolved oxygen as 5-year moving average from 1965 to 2015 with uncertainty estimates. The uncertainty can come from unresolved small-scale and high-frequency variability and mapping errors. The multi-decadal trend of global dissolved oxygen is in the range of -281 to -373 Tmol/decade. This estimate is more conservative than previous works. In this study, the grid points far from the observations are essentially set equal to zero anomaly from the climatology. Calculating global inventory with this approach produces a relatively conservative estimate; thus, the results from this study likely provide a useful lower bound estimate of the global oxygen loss.

54 ENVIRONMENTAL SCIENCES↗

Stochastic multiscale flux basis for Stokes-Darcy flows

Three algorithms are developed for uncertainty quantification in modeling coupled Stokes and Darcy flows. The porous media may consist of multiple regions with different properties. The permeability is modeled as a non-stationary stochastic variable, with its log represented as a sum of local Karhunen-Loève (KL) expansions. The problem is approximated by stochastic collocation on either tensor-product or sparse grids, coupled with a multiscale mortar mixed finite element method for the spatial discretization. A non-overlapping domain decomposition algorithm reduces the global problem to a coarse scale mortar interface problem, which is solved by an iterative solver, for each stochastic realization. In the traditional domain decomposition implementation, each subdomain solves a local Dirichlet or Neumann problem in every interface iteration. To reduce this cost, two additional algorithms based on deterministic or stochastic multiscale flux basis are introduced. The basis consists of the local flux (or velocity trace) responses from each mortar degree of freedom. It is computed for each subdomain independently before the interface iteration begins. The use of the multiscale flux basis avoids the need for subdomain solves on each iteration. The deterministic basis is computed at each stochastic collocation and used only at this realization. The stochastic basis is formed by further looping over all local realizations of a subdomain's KL region before the stochastic collocation begins. It is used over multiple realizations. Numerical tests are presented to illustrate the performance of the three algorithms, with the stochastic multiscale flux basis showing significant savings in computational cost compared to the other two algorithms.

97 MATHEMATICS AND COMPUTING↗

Scattering of radiofrequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radiofrequency (RF) electromagnetic waves that propagate through this region. It is important to know the spectral properties of these scattered RF waves, whether used for diagnostics or for heating and current drive. The spectral changes influence the interpretation of the obtained diagnostic data, and the current and heating profiles. A full-wave, three-dimensional (3-D) electromagnetic code ScaRF has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method used for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma by the anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can also be used for the study of the scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U and DIII-D. For the case of medium-sized tokamaks, there is experimental evidence that drift waves and rippling modes are present in the edge region. Hence, we have studied the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes. The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion method in conjunction with the Smolyak sparse-grid integration, which is known as the PCE-SG method. The PCE-SG method is proven to be accurate and more efficient compared with alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhancing Active Distribution Systems Resilience by Fully Distributed Self-Healing Strategy

Distributed restoration can exploit smart grid technologies to enhance the resilience of active distribution networks toward a self-healing smart grid. However, the large number of decision variables, especially the binary ones for reconfiguration, bring challenges to developing scalable distributed distribution service restoration (DDSR) strategies. This paper proposes a fully distributed solution procedure based on the alternating direction method of multipliers (ADMM) for mixed-integer programming problems and applies to develop the DDSR framework. The method consists of relax-drive-polish phases, 1) relaxing binary variables, and applying the convex ADMM as a warm start; 2) driving the solutions toward Boolean values through a proximal operator; 3) fixing the obtained binding binary variables and solving the rest of the problem to polish results and achieve a high-quality suboptimal solution. Then, an autonomous clustering strategy and consensus ADMM are integrated with the proposed method to realize the fully distributed cluster-based framework of DDSR. This framework can first determine DER scheduling and switch status for reconfiguration to energize the out-of-service areas from local faults, and then provide the load restoration solution in a distributed manner for total blackouts in large-scale distribution networks. Furthermore, the effectiveness and scalability of the proposed DDSR framework are demonstrated through testing on the IEEE 123-node, IEEE 8500-node, and synthetic 100k-node test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SCOPE Climate: a 142-year daily high-resolution ensemble meteorological reconstruction dataset over France

SCOPE Climate (Spatially COherent Probabilistic Extended Climate dataset) is a 25-member ensemble of 142-year daily high-resolution reconstructions of precipitation, temperature, and Penman–Monteith reference evapotranspiration over France, from 1 January 1871 to 29 December 2012. SCOPE Climate provides an ensemble of 25 spatially coherent gridded multivariate time series. It is derived from the statistical downscaling of the Twentieth Century Reanalysis (20CR) by the SCOPE method, which is based on the analogue approach. SCOPE Climate performs well in comparison to both dependent and independent data for precipitation and temperature. The ensemble aspect corresponds to the uncertainty related to the SCOPE method. SCOPE Climate is the first century-long gridded high-resolution homogeneous dataset available over France and thus has paved the way for improving knowledge on specific past meteorological events or for improving knowledge on climate variability, since the end of the 19th century. This dataset has also been designed as a forcing dataset for long-term hydrological applications and studies of the hydrological consequences of climate variability over France. SCOPE Climate is freely available for any non-commercial use and can be downloaded as NetCDF files from https://doi.org/10.5281/zenodo.1299760 for precipitation, https://doi.org/10.5281/zenodo.1299712 for temperature, and https://doi.org/10.5281/zenodo.1251843 for reference evapotranspiration.

54 ENVIRONMENTAL SCIENCES↗

Modeling the Grid Integration of Variable Renewable Energy (VRE) Resources

The presentation was given at the November 21-23, 2018 'CEM Days - Integration of Renewables in the Electric Sector: Paths and Challenges to Energy Planning' in Rio de Janeiro, Brazil. It provides a summary of NREL's current modeling approaches, datasets, and considerations for integrating renewables into the power system. It also covers current model method improvement efforts, including our vision for a truly integrated modeling approach.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Advances in solar forecasting: Computer vision with deep learning

Renewable energy forecasting is crucial for integrating variable energy sources into the grid. It allows power systems to address the intermittency of the energy supply at different spatiotemporal scales. To anticipate the future impact of cloud displacements on the energy generated by solar facilities, conventional modeling methods rely on numerical weather prediction or physical models, which have difficulties in assimilating cloud information and learning systematic biases. Augmenting computer vision with machine learning overcomes some of these limitations by fusing real-time cloud cover observations with surface measurements acquired from multiple sources. This Review summarizes recent progress in solar forecasting from multisensor Earth observations with a focus on deep learning, which provides the necessary theoretical framework to develop architectures capable of extracting relevant information from data generated by ground-level sky cameras, satellites, weather stations, and sensor networks. Overall, machine learning has the potential to significantly improve the accuracy and robustness of solar energy meteorology; however, more research is necessary to realize this potential and address its limitations.

14 SOLAR ENERGY↗

Effect of the Grell-Freitas Deep Convection Scheme in Quasi-Uniform and Variable-Resolution Aquaplanet CAM Simulations

The Grell-Freitas deep convection scheme (GF) has been implemented in the Community Atmosphere Model (CAM) with the nonhydrostatic Model for Prediction Across Scales (MPAS) dynamical core. Aquaplanet simulations have been performed with a global quasi-uniform resolution mesh with ~120 km grid spacing and a variable resolution mesh employing a circular refined region with ~30 km grid spacing centering at the equator and ~120 km elsewhere. GF produces overall patterns of the general circulations comparable to those of other aquaplanet experiments, including those produced using the default Zhang-McFarlane convection scheme (ZM) in CAM Version 5.4. GF alleviates the overestimated occurrence of weak intensity precipitation characteristic of ZM due to GF’s different closure methods and deep-convection trigger conditions. Compared to ZM, the GF characteristics of the equatorial Kelvin waves are closer to observations. GF also better simulates precipitation at variable resolution, with no unrealistic wave distortion by the refinement. Overall, GF shows less sensitivity of convective precipitation and heating/moistening tendencies to changes in horizontal resolution compared to ZM. However, vertically-overlapped total cloud fraction simulated with GF in the tropics is more sensitive to the resolution than ZM due to the larger resolution-sensitivity of high-level cloud fraction in the tropics.

54 ENVIRONMENTAL SCIENCES↗

How do modeling decisions affect the spread among hydrologic climate change projections? Exploring a large ensemble of simulations across a diversity of hydroclimates

Methodological choices can have strong effects on projections of climate change impacts on hydrology. In this study, we investigate the ways in which four different steps in the modeling chain influence the spread in projected changes of different aspects of hydrology. To form the basis of these analyses, we constructed an ensemble of 160 simulations from permutations of two Representative Concentration Pathways, 10 global climate models, two downscaling methods, and four hydrologic model implementations. The study is situated in the Pacific Northwest of North America, which has relevance to a diverse, multinational cast of stakeholders. We analyze the effects of each modeling decision on changes in gridded hydrologic variables of snow water equivalent and runoff, as well as streamflow at point locations. Results show that the choice of representative concentration pathway or global climate model is the driving contributor to the spread in annual streamflow volume and timing. On the other hand, hydrologic model implementation explains most of the spread in changes in low flows. Finally, by grouping the results by climate region the results have the potential to be generalized beyond the Pacific Northwest. Future hydrologic impact assessments can use these results to better tailor their modeling efforts.

58 GEOSCIENCES↗

Deep Reinforcement Learning for Microgrid Cost Optimization Considering Load Flexibility

This paper proposes a novel Soft-Actor-Critic (SAC) based Deep Reinforcement Learning (DRL) method for optimizing the cost of microgrid operation by leveraging load flexibility. The proposed SAC-DRL method is designed to coordinate the control of distributed energy resources (DERs) and flexible load, addressing practical energy billing formation by power distribution utilities. Key contributions include an innovative reward function to mitigate sparse reward challenges and a mixed control strategy for discrete and continuous variables, ensuring radial network topology and minimizing power loss. We evaluate the proposed method on the model of a real microgrid located in Southern California, U.S.. The SAC-DRL model is tested to demonstrate its efficacy in reducing grid dependence, optimizing resource use, and minimizing costs. The results highlight the potential of DRL in modern energy systems, offering a sustainable and economically efficient solution for energy management in microgrids.

deep reinforcement learning↗

Feedback control of stored energy and rotation with variable beam energy and perveance on DIII-D

Design and testing of a new method for real-time control on DIII-D using the new voltage and perveance variation capability of DIII-D's neutral beam injection system is presented. The approach enables control of power, torque, stored energy, and rotation without requiring pulse-width modulation of the beams. Such modulations are perturbative to experiments, can cause beam reliability problems, and could potentially damage or fatigue accelerator grids, especially looking forward to active control of long-pulse devices and reactors. Here, the control algorithm accounts for beam line loss mechanisms and includes approaches that address the complexity of having three manipulated variables for each of the eight beams and compensate for the rate and magnitude limits affecting the voltage and perveance actuators. Results are presented that demonstrate the performance of the approach in the first feedback control experiments to use the expanded controllability of the DIII-D beams.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sulfur isotope analysis of microcrystalline iron sulfides using secondary ion mass spectrometry imaging: Extracting local paleo-environmental information from modern and ancient sediments

Rationale: Sulfur isotope ratio measurements of bulk sulfide from marine sediments have often been used to reconstruct environmental conditions associated with their formation. In situ microscale spot analyses by secondary ion mass spectrometry (SIMS) and laser ablation multiple-collector inductively coupled plasma mass spectrometry (LA-MC-ICP-MS) have been utilized for the same purpose. However, these techniques are often not suitable for studying small (≤10 μm) grains or for detecting intra-grain variability. Methods: In this work, we present a method for the physical extraction (using lithium polytungstate heavy liquid) and subsequent sulfur isotope analysis (using SIMS; CAMECA IMS 7f-GEO) of microcrystalline iron sulfides. SIMS sulfur isotope ratio measurements were made via Cs + bombardment of raster squares with sides of 20–130 μm, using an electron multiplier (EM) detector to collect counts of 32 S - and 34 S - for each pixel (128 × 128 pixel grids) for between 20 and 960 cycles. Results: The extraction procedure did not discernibly alter pyrite grain-size distributions. The apparent inter-grain variability in 34 S/ 32 S in 1–4 μm-sized pyrite and marcasite fragments from isotopically homogeneous hydrothermal crystals was ~ ±2‰ (1σ), comparable with the standard error of the mean for individual measurements (≤ ±2‰, 1σ). In contrast, grain-specific 34 S/ 32 S ratios in modern and ancient sedimentary pyrites and marcasites can have inter- and intra-grain variability >60‰. The distributions of intra-sample isotopic variability are consistent with bulk 34 S/ 32 S values. Conclusions: SIMS analyses of isolated iron sulfide grains yielded distributions that are isotopically representative of bulk 34 S/ 32 S values. Populations of iron sulfide grains from sedimentary samples record the evolution of the S-isotopic composition of pore water sulfide in their S-isotopic compositions. These data allow past local environmental conditions to be inferred.

59 BASIC BIOLOGICAL SCIENCES↗

Scattering of radio frequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radio frequency (RF) electromagnetic waves that propagate through this region. It is important to know, whether used for diagnostics or for heating and current drive, the spectral properties of these scattered RF waves. The spectral changes influences the interpretation of the diagnostic-data obtained and the current and heating profiles. A full-wave, 3D electromagnetic code ScaRF (see Papadopoulos et al. 2019) has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma, anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can be for the study of scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U, DIII-D. For the case of medium-sized tokamaks, there's experimental evidence that drift waves and rippling modes are present in the edge region (see Ritz et al. 1984). Hence, we study the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes (see Papadopoulos et al. 2019). The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion (see Xiu & Karniadakis 2002) method in conjunction with the Smolyak sparse grid integration (see Papadopoulos et al. 2018) (PCE-SG). The PCE-SG method is proven accurate and much more efficient (roughly 2-orders of magnitude shorter execution time) compared to alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Review of data-driven models for quantifying load shed by non-residential buildings in the United States

Shifting and shedding power demand in buildings can be cost-effective techniques for grids to function reliably and for end users to earn compensation. Grid operators reimburse customers in proportion to the quantity of load shed. Simple data-driven methods are used to quantify this shed, which is the difference between a measured load during the event and modeled "baseline" that would have occurred in absence of the event. These methods have evolved over the years and in many cases have been integrated with building physics, to make them a hybrid between physics based and empirical models. However, there is no comprehensive analysis that provides guidance to building operators, grid operators and researchers in selecting appropriate models based on their specific needs and available data. Here, this work aims to fill this gap by critically assessing the performance of baseline models put forward from the year 2000 through 2023. The literature reviewed includes reports generated by grid operators, reports from national laboratories and academic journal articles. The work outlines modeling features like the inputs, training period, estimation method, adjustments to fine tune the predictions and metrics to evaluate the performance. A comprehensive list of 50 models has been provided. For each model, the study explores the applicability of the model to weather sensitive buildings, variability in the building profile, timing of the event, and whether the building reduces energy consumption before an event. The work identifies the situations in which a particular model works and draws lessons based on evidence of performance. Finally, recommendations to aid in model selection are given.

97 MATHEMATICS AND COMPUTING↗

Transmission Expansion Planning Test System for AC/DC Hybrid Grid With High Variable Renewable Energy Penetration

Test systems are important tools and benchmark for power system research. Currently, there is a lack of standard test systems for modern transmission expansion planning (TEP) research, especially under high variable renewable energy (VRE) penetrations. This article describes a 38-bus test system named the HRP-38 system dedicated to TEP. “HRP” stands for high renewable penetration. The objective of establishing such a system is to provide a consistent platform for different TEP methods to be tested and compared. The test system considers high VRE penetration of more than 30% energy share, which will be an important feature in many power systems in the future. Both AC and DC candidate lines are given to provide sufficient transmission planning alternatives. The complexity of the test system is well balanced considering calculation tractability and the ability to test the performance of the TEP planning model. The network topology, generation mix and load characteristics are described in detail. Year-round time series of VRE (wind power and PV) output and load curves for each bus are provided to facilitate a detailed evaluation of the transmission grid configuration under high VRE penetration. Finally, four benchmark results with different optimization settings are provided for the HRP-38 system.

AC/DC hybrid↗