Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Benchmarking Microscale Ductility Measurements (Final Report of the Project DE-NE0008799)

Conventional macroscale experimentation is generally considered to be straightforward with few limitations. Conversely, micro/nanoscale experimentation presents numerous challenges in loading device design, sample preparation and handling, as well as accurate understanding of grain size and local texture effects on recorded measurements. Despite these challenges, nanopillar compression, MEMs based micro-tension, and nanoindentation approaches have been able to provide fundamental contributions to the understanding of material behavior at small lengthscales. However, the overarching shortcoming of these micro/nanoscale experimentation approaches, is the inability to directly translate measurements evaluated at the nm and µm length scales (e.g., hardness) to macroscale tensile material behavior (i.e., elastic modulus, yield strength, and ductility). The objectives of the proposed study are, 1) to establish best practices for obtaining tensile microscale ductility measurements, and 2) to validate methodologies to for comparing microscale ductility measurements to macroscale ductility measurements. In order to achieve these objectives, a multi-lengthscale, multi-temperature testing protocol and simulation framework are executed first on copper as a model material to validate the following approach, and second on reactor grade Zircaloy-2. Experiments are conducted on specimens extracted from the same test piece to ensure nominally identical grain size and texture from specimens to specimen. Motivated by the need to isolate the contribution of size-effects on obtained mechanical property measurements, specimens are manufactured with thicknesses at the micro- (1-10 µm), meso- (10-100's µm), and macroscales (sub-sized ASTM E8). In-situ full-field deformation techniques (scanning electron microscopy (SEM) grid methods and optical DIC) are incorporated into testing at each specimen length-scale to capture plasticity localization and evolution. Experimental testing for all specimens is conducted at both room temperature and elevated temperatures to probe the role of thermal activation on plastic deformation accommodation processes. Simulation efforts focus on examining the mechanical behavior of microscale specimens using a finite element approach with explicitly resolved grain morphologies, and an embedded crystal plasticity model. The cost-efficient implementation method allows for the modeling of a statistically significant number of both real (i.e., digital twin) and generated microstructures to obtain an understanding of the interrelationships between specimen microstructure and geometric variables (grain size, texture, specimen geometry, etc.) on microscale mechanical behavior.

36 MATERIALS SCIENCE↗

Optimization and Multimachine Learning Algorithms to Predict Nanometal Surface Area Transfer Parameters for Gold and Silver Nanoparticles

Interactions between gold metallic nanoparticles and molecular dyes have been well described by the nanometal surface energy transfer (NSET) mechanism. However, the expansion and testing of this model for nanoparticles of different metal composition is needed to develop a greater variety of nanosensors for medical and commercial applications. In this study, the NSET formula was slightly modified in the size-dependent dampening constant and skin depth terms to allow for modeling of different metals as well as testing the quenching effects created by variously sized gold, silver, copper, and platinum nanoparticles. Overall, the metal nanoparticles followed more closely the NSET prediction than for Förster resonance energy transfer, though scattering effects began to occur at 20 nm in the nanoparticle diameter. To further improve the NSET theoretical equation, an attempt was made to set a best-fit line of the NSET theoretical equation curve onto the Au and Ag data points. An exhaustive grid search optimizer was applied in the ranges for two variables, 0.1≤C≤2.0 and 0≤α≤4, representing the metal dampening constant and the orientation of donor to the metal surface, respectively. Three different grid searches, starting from coarse (entire range) to finer (narrower range), resulted in more than one million total calculations with values C=2.0 and α=0.0736. The results improved the calculation, but further analysis needed to be conducted in order to find any additional missing physics. With that motivation, two artificial intelligence/machine learning (AI/ML) algorithms, multilayer perception and least absolute shrinkage and selection operator regression, gave a correlation coefficient, R2, greater than 0.97, indicating that the small dataset was not overfitting and was method-independent. This analysis indicates that an investigation is warranted to focus on deeper physics informed machine learning for the NSET equations.

Demers, Steven M. E. (ORCID:0000000192213246)↗

Evaluating fine-resolution, regional outputs of a variable resolution global climate model

Climate models have been used to study water resources and regional hydrologic responses to climate change, but climate model outputs must be downscaled to provide relevant regional data. However, the accuracy of this regional data is limited by uncertainties across and within downscaling methods, uncertainty across global outputs, and discontinuities at downscaled boundaries. A new alternative to traditional downscaling is a variable resolution model that incorporates fine-resolution regions directly into a coarse-resolution, global climate simulation in order to capture contiguous dynamics across resolution boundaries. In this study, we used the Variable-Resolution Community Earth System Model (VR-CESM) to generate one-eighth degree (14 km) fine-resolution outputs for the western U.S. and eastern China from 1970-2006. We focus our evaluation on precipitaiton, temperature, snow pack, solar radiation, and wind. We compare the model outputs with remote-sensing-based precipitation data, and both reanalysis and gridded weather station data for precipitation and temperature. VR-CESM generally has a cold bias in winter and a warm bias in summer in the western U.S., which compensate each other to reduce the annual bias. In eastern China, however, the sign of temperature biases are more consistent throughout the year with cold biases in the higher mountains and warm biases throughout most of the rest of the region. Precipitation biases are dependent upon reference data, and show slight overestimation in high mountain regions in both the U.S. and China with respect to gridded weather station data. Simulated snow cover in the western U.S. is reasonable compared to remote sensing data, but snow cover and snow water equivalent have larger biases when compared to reanalysis data. In eastern China there are widespread snow cover biases compared to remote sensing data. VR-CESM underestimates downward shortwave radiation to a greater degree in summer than in winter, and underestimates surface layer windspeed over mountains to a greater degree than in other areas. Comparison between VR-CESM and a coarser simulation (1-degree Beijing Climate Center model) shows reduced precipitation biases in the mountainous regions with finer resolution, indicating the value of variable-resolution modeling for reigonal studies.

Di Vittorio, Alan↗

Voltron Compatible Whole Building Root-Fault Detection and Diagnosis

Buildings consume more than 40% of primary energy in the U.S. Malfunctioning sensors, components, and control systems, as well as degrading components in Heating, Ventilating and Air-conditioning (HVAC) and lighting systems are main reasons for energy waste and unsatisfactory indoor environment. Field studies have shown that an energy saving of 5-30% of total building energy consumption and improved indoor air quality can be achieved by simply applying automated fault detection and diagnosis (AFDD), followed by corrections even if these are not done in real time. Extensive research has been made on the development of component level AFDD tools in the past two decades. However, in recent years, the development of method for a whole building level AFDD has received much attention. Here, a whole building fault refers to a fault happening in one subsystem but impact on more than one subsystems, or has a significant impact on building performances, such as building’s energy consumption. Existing component level AFDD tools could fail to detect or give false alarm for faults that trigger abnormalities in multiple subsystems. Isolating such whole building faults, which have impacts on multiple subsystems, is also challenging by simply using component based AFDD solutions. On the other hand, it might be more cost-effective for medium- or small-sized commercial buildings that have subsystems (e.g. primary cooling subsystem, air loop subsystem, heating subsystem, lighting subsystem, etc.) to install one whole building level AFDD tool, rather than many component level AFDD tools, to focus on significant faults. The focus of this project is firstly to provide training opportunities in building AFDD research area for undergraduate students and to develop a workforce pipeline in building sciences. The objective of the proposed research activities is to develop cost-effective and scalable AFDD solutions for whole building level faults. The developed AFDD solution is VOLTTRONTM compatible to further increase its plug-n-play capacity and the market penetration. VOLTTRONTM is an open-source innovative distributed control and sensing software platform developed by the Pacific Northwest National Laboratory. A VOLTTRONTM compatible AFDD tool brings the benefits such as low-cost deployment, and better building to grid integration, etc. The performance of the developed whole building AFDD tools is evaluated using data collected from a mixed-use medium-sized university campus building. Another significant goal of this project is to engage undergraduate students in research activities and expose them to building control and AFDD fields. A demonstration building – Nesbitt Hall at Drexel University is identified for this study. This demonstration building has a typical HVAC system configuration, which includes a water-cooled chiller system, three variable air volume (VAV) air handling unit (AHU) systems, and a hydronic heating system. Faults, that are expected to have a whole building level impact, are artificially implemented in this demonstration building for three different seasons. Building data have been collected in this project, which include fault free building automation system (BAS) data and BAS data that contain system behaviors caused by artificially implemented and naturally occurred faults. These collected data are used to evaluate the developed whole building AFDD tools. A data-driven method, which includes a weather and schedule-based pattern matching (WPM) method and Feature based Principal Component Analysis (FPCA) method, is developed for whole building level fault detection. Parameter sensitivity tests, including snapshot window size, data sample searching pool size, etc., are implemented to evaluate the impacts of various WPM-FPCA parameters. A data-driven and expert knowledge/rule-based method using Bayesian Network (BN) and WPM is developed for whole building level fault diagnosis. The developed WPM-BN method, which includes a two-layer BN structure model and BN parameters, are either learned from baseline data or developed from expert knowledge. The developed whole building AFDD tools are evaluated using BAS data collected from the demonstration building. Using the evaluation data set, the WPM-FPCA fault detection method is evaluated under two different Principle Component (PC) retention rates, i.e., 0.65 retention rate and 0.95 retention rate. It is observed that when using 0.65 PC retention rate, the WPM-FPCA fault detection method achieves a 79% fault detection rate (this rate reaches to 85% when including the six make-up fault test cases discussed in the next bullet), while when using 0.95 PC retention rate, the WPM-FPCA fault detection method achieves 100% fault detection rate (this rate remains to be 100% when including the six make-up fault test cases). Under both PC retention rates, the WPM-FPCA fault detection method demonstrates 0% false alarm rates. Using the evaluation data set, the developed WPM-BN fault diagnosis method has successfully diagnosed thirteen (13) faults out of the fourteen (14) fault test cases. The one mis-diagnosed case exhibits fewer fault symptoms in different subsystems, when compared with other fault test cases. Overall speaking, the developed whole building AFDD solutions demonstrate satisfactory fault detection/diagnosis accuracy with very low false alarm rate. Lastly, a survey and market study has been performed to understand the existing AFDD market, including its gaps and needs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Unified 0.25-degree gridded infrastructure-critical extreme weather for the United States from 1979 to 2100

Extreme weather events can severely disrupt critical infrastructure, triggering cascading effects on power, transportation, and essential services. However, standard weather and climate datasets often lack specialized variables necessary for hazard assessments. We present a unified dataset of infrastructure-critical weather and climate variables across the United States at 0.25° resolution, covering daily or sub-daily intervals from 1979 to 2100. The dataset includes temperature, dew point, wind gusts, precipitation partitioned by rain, snow, and freezing rain or ice pellets, lightning, and wildfire metrics. Historical conditions (1979-2023) are synthesized from observations and reanalysis products, while future projections are derived from 14 CMIP6 global climate models (historical, SSP245, and SSP585 experiments). Physically based and data-driven methods are used to estimate variables not directly provided by existing models. By integrating these variables into a single unified dataset, we enable consistent, high-resolution assessments of weather-related infrastructure risks across past and future periods, supporting wide-ranging applications in energy, transportation, water resources, emergency management, and beyond.

Climate and Earth system modelling↗

Data-Driven Day-Ahead PV Estimation Using Hybrid Deep Learning

Ongoing smart grid activities and associated automation resulted in rich set of data. Here, these data can be utilized for monitoring and estimation of real time photovoltaic (PV) generation. Inherent variability in PV and related impact on power systems is a challenging problem. Improving the accuracy of PV generation estimation is beneficial for both the PV owners and the grid operators. Recently, deep learning algorithms possible by the availability of data have shown its advantages for time series estimation; however, its application on PV generation estimation is still in the early stage. In this paper, a hybrid estimation model with a combination of long-short-term-memory network (LSTM) and persistence model (PM) is developed to provide day-ahead PV estimation at 15-minute time interval with high accuracy and robustness. Simulation results show the superior performance of the proposed method over existing methods for most of the test cases.

24 POWER TRANSMISSION AND DISTRIBUTION↗

PATCHWORK: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe Fluid Simulations

Here, we present a "multipatch" infrastructure for the numerical simulation of fluid problems in which subregions require different grid scales, different grid geometries, different physical equations, or different reference frames. Its key element is a sophisticated client–router–server framework for efficiently linking processors supporting different regions ("patches") that must exchange boundary data. This infrastructure may be used with a wide variety of fluid dynamics codes; the only requirement is that their primary dependent variables be the same in all patches, e.g., fluid mass density, internal energy density, and velocity. Its structure can accommodate either Newtonian or relativistic dynamics. The overhead imposed by this system is both problem and computer cluster architecture dependent. Compared to a conventional simulation using the same number of cells and processors employed on a problem not requiring multipatch methods, the cell update per processor rate decreases by an amount that can range from negligible to a factor of a few; however, even in these problems, the infrastructure can permit substantial decreases in the total number of cell updates required.

97 MATHEMATICS AND COMPUTING↗

Development of observation-based global multilayer soil moisture products for 1970 to 2016

Abstract. Soil moisture (SM) datasets are critical to understanding the global water, energy, and biogeochemical cycles and benefit extensive societal applications. However, individual sources of SM data (e.g., in situ and satellite observations, reanalysis, offline land surface model simulations, Earth system model – ESM – simulations) have source-specific limitations and biases related to the spatiotemporal continuity, resolutions, and modeling and retrieval assumptions. Here, we developed seven global, gap-free, long-term (1970–2016), multilayer (0–10, 10–30, 30–50, and 50–100 cm) SM products at monthly 0.5∘ resolution (available at https://doi.org/10.6084/m9.figshare.13661312.v1; Wang and Mao, 2021) by synthesizing a wide range of SM datasets using three statistical methods (unweighted averaging, optimal linear combination, and emergent constraint). The merged products outperformed their source datasets when evaluated with in situ observations (mean bias from −0.044 to 0.033 m3 m−3, root mean square errors from 0.076 to 0.104 m3 m−3, Pearson correlations from 0.35 to 0.67) and multiple gridded datasets that did not enter merging because of insufficient spatial, temporal, or soil layer coverage. Three of the new SM products, which were produced by applying any of the three merging methods to the source datasets excluding the ESMs, had lower bias and root mean square errors and higher correlations than the ESM-dependent merged products. The ESM-independent products also showed a better ability to capture historical large-scale drought events than the ESM-dependent products. The merged products generally showed reasonable temporal homogeneity and physically plausible global sensitivities to observed meteorological factors, except that the ESM-dependent products underestimated the low-frequency temporal variability in SM and overestimated the high-frequency variability for the 50–100 cm depth. Based on these evaluation results, the three ESM-independent products were finally recommended for future applications because of their better performances than the ESM-dependent ones. Despite uncertainties in the raw SM datasets and fusion methods, these hybrid products create added value over existing SM datasets because of the performance improvement and harmonized spatial, temporal, and vertical coverages, and they provide a new foundation for scientific investigation and resource management.

54 ENVIRONMENTAL SCIENCES↗

An Uncertainty Management Framework for Integrated Gas-Electric Energy Systems

In many parts of the world, electric power systems have seen a significant shift toward generation from renewable energy and natural gas. Because of their ability to flexibly adjust power generation in real time, gas-fired power plants are frequently seen as the perfect partner for variable renewable generation. However, this reliance on gas generation increases interdependence and propagates uncertainty between power grids and gas pipelines and brings coordination and uncertainty management challenges. To address these issues, we propose an uncertainty management framework for uncertain, but bounded gas consumption by gas-fired power plants. The admissible ranges are computed based on a joint optimization problem for the combined gas and electricity networks, which involves chance-constrained scheduling for the electric grid and a novel robust optimization formulation for the natural-gas network. This formulation ensures feasibility of the integrated system with a high probability, while providing a tractable numerical formulation. A key advance with respect to existing methods is that our method is based on a physically accurate, validated model for transient gas pipeline flows. Our case study benchmarks our proposed formulation against methods that ignore how reserve activation impacts the fuel use of gas power plants and only consider predetermined gas consumption. Here, the results demonstrate the importance of considering uncertainty to avoid operating constraint violations and curtailment of gas to the generators.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Economic analysis of a nuclear hybrid energy system in a stochastic environment including wind turbines in an electricity grid

The electricity market in the U.S. is moving toward an energy mix including more variable renewable energy (VRE) sources like wind and solar. This increasing VRE penetration is altering the profile of the net electricity demand. Normally, flexible generators (e.g. gas turbines) are used to absorb fluctuations in the electricity demand. Therefore, an increase in VRE also leads to an increase in installed capacity for such flexible generators which, in turn, may increase the price of electricity. Nuclear-Renewable Hybrid Energy Systems (N-R HES) are seen as a solution to this problem. The N-R HES contains a nuclear power plant electrically and/or thermally linked to an industrial process (IP). This paper demonstrates a new methodology and presents the corresponding RAVEN/Modelica-based software framework to evaluate the financial performance of N-R HES. The novelty of the method is that it explicitly incorporates the stochastic nature of N-R HES inputs such as electricity demand and wind speeds. The method is applied to a generic (rather than region-specific) N-R HES. The analysis shows that under the right market conditions, economic benefits can be achieved by adding a suitable IP to the base-load generator. In particular, there exists an optimal relationship between the IP capacity and the demand. Furthermore, global trends indicate VRE penetration raises the cost of electricity due to the added volatility, an issue only mitigated when operating with high profit margins for the IP.

97 MATHEMATICS AND COMPUTING↗

Optimization of turbine design in wind farms with multiple hub heights, using exact analytic gradients and structural constraints

Wind farms are generally designed with turbines of all the same hub height. If wind farms were designed with turbines of different hub heights, wake interference between turbines could be reduced, lowering the cost of energy (COE). This paper demonstrates a method to optimize onshore wind farms with two different hub heights using exact, analytic gradients. Gradient-based optimization with exact gradients scales well with large problems and is preferable in this application over gradient-free methods. Our model consisted of the following: a version of the FLOw Redirection and Induction in Steady-State wake model that accommodated three-dimensional wakes and calculated annual energy production, a wind farm cost model, and a tower structural model, which provided constraints during optimization. Structural constraints were important to keep tower heights realistic and account for additional mass required from taller towers and higher wind speeds. We optimized several wind farms with tower height, diameter, and shell thickness as coupled design variables. Our results indicate that wind farms with small rotors, low wind shear, and closely spaced turbines can benefit from having two different hub heights. A nine-by-nine grid wind farm with 70-meter rotor diameters and a wind shear exponent of 0.08 realized a 4.9% reduction in COE by using two different tower sizes. If the turbine spacing was reduced to 3 diameters, the reduction in COE decreased further to 11.2%. Allowing for more than two different turbine heights is only slightly more beneficial than two heights and is likely not worth the added complexity.

17 WIND ENERGY↗

Optimizing the physical design and layout of a resilient wind, solar, and storage hybrid power plant

We report as wind and solar technologies improve and their costs decrease, the share of power produced by these sources will increase. As the market penetration increases, these power sources will need to provide grid services, such as dispatchability, in addition to providing energy. One way to reduce variability, provide higher quality power to the grid, and address local grid stability issues is through colocating wind and solar power plants. In addition to operating reliably during normal operating conditions, in scenarios with high penetrations of renewable generation, it is important that these hybrid plants can withstand production disruptions and continue to supply power despite prolonged resource reduction, extreme weather events, or other disruptions. In this paper, we present a methodology to optimize a wind-solar-battery hybrid power plant down to the component level that is resilient against production disruptions and that can continually produce some minimum required power. We introduce the models and assumptions we used to simulate a hybrid power plant as well as the design variable parameterization and specific methods we used to optimize the plant. We demonstrate the performance of our method by comparing a plant optimized for different objectives, generation outage durations, minimum power requirements, and power purchase agreements. Although the plant design is sensitive to model parameters and various other assumptions, our results demonstrate some of the optimal designs that occur in different scenarios and what one should expect when designing a hybrid wind-solar-storage power plant.

14 SOLAR ENERGY↗

Challenges and alternatives to empirical orthogonal functions for earth system data

Empirical orthogonal functions (EOFs) applied to gridded Earth system data enables users to diagnose modes of variability with relative ease. Yet, many challenges to interpretation exist such that they must be used with awareness and intention when applied to gridded climate data, especially with large ensembles. Utilizing data from two different Earth system modelling large ensemble frameworks, the Energy Exoscale Earth System Model and the Community Earth System Model, as well as reanalysis data, common EOF pitfalls are summarized and discussed. Challenges include erroneous mode swapping, sign flipping, and the temporal variability of the centers of action. For modes of variability with similar contribution to variance, mode swapping is not uncommon. Sign flipping can occur with almost any mode where the pattern is correct, but the sign is arbitrary. Although the variability of the center of action is not necessarily problematic, it potentially complicates interpretation over multi-century timescales. A wide variety of alternative methods to EOFs exist, but fitness-for-purpose must be evaluated. Additionally, illustrations of alternative methods and examples of proper use are provided. Alternative methods fit into three categories: EOF variants, linear methods, and multilinear methods.

54 ENVIRONMENTAL SCIENCES↗

Optimal Control of Differentially Private EV Charging: A Scalable Learning Approach Under Uncertainty

Internet of Things (IoT)-enabled electric vehicles (IoEVs) enable intelligent charging coordination that accounts for grid congestion. However, increased data exchange raises privacy concerns, as charging patterns can reveal sensitive driver behavior to grid operators. Here, we propose a differentially private (DP) EV charging framework that enables coordinated control while protecting driver data with theoretical privacy guarantees. Nevertheless, integrating DP inevitably introduces uncertainty into the control strategy for EVs, which can lead to infeasible solutions. To tackle this challenge, we develop a feasible and scalable control algorithm based on constrained reinforcement learning (CRL) and convex hulls. While our framework is designed to handle the uncertainty introduced by DP, it is general and also applicable to other sources of uncertainty in EV charging, such as the stochastic nature of driver behavior and renewable variability. This ensures feasible and privacy-preserving coordination of EV charging at scale. Our method constructs convex hulls within the action space to guarantee feasibility under stochastic constraints and incorporates constraint reduction techniques to improve scalability. Case studies based on IEEE benchmark systems demonstrate that the proposed approach effectively balances feasibility under uncertainty, scalability, and privacy in large-scale EV charging control.

Engineering - Power transmission and distribution↗

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 1—Theory

Summary Core-scale measurements are considered the ground truth that oil and gas or electric utility companies use to predict the migration of fluids such as oil, natural gas, carbon dioxide, or brine deep underground during their extraction or injection operations. To provide a greater understanding of petrophysical properties of low-permeability geologic formations such as shales and tight gas sandstones, this study introduces a novel core-analysis procedure. The technique follows conventional pressure-pulse-decay permeametry, where the pressure in an inlet chamber adjacent to a cylindrical core plug undergoes a rapid pressurization, the system is shut in, and the pressure reaches a new equilibrium. However, unlike a standard unidirectional pressure-pulse decay, the full-immersion pressure-pulse decay (Hannon 2019) applies a pressure disturbance to the entire outer surface area of the sample. This article covers the numerical simulator designed to model flow through an anisotropic porous sample in this scenario. The model assumes distinct but uniform permeabilities along the radial and axial directions of the cylindrical plug sample. When extracting a plug vertically (or perpendicular to bedding), the permeability along the radial direction associates with the horizontal permeability (i.e., parallel to bedding), whereas in the axial direction, flow occurs perpendicular to bedding (a vertical permeability). Investigations of these model outputs demonstrate an approximately 20-fold decrease in time to complete a full-immersion experiment compared with conventional pressure-pulse decay. Furthermore, the pressure-decay curves resulting from the full-immersion method have slightly different shapes than those resulting from other unidimensional transient methods. These differences begin to demonstrate that, under achievable experimental conditions, the analysis of pressure data from one full-immersion test could enable the simultaneous estimation of the apparent permeabilities parallel and perpendicular to bedding of a cylindrical sample in addition to its porosity. A follow-up article finalizes the proof of this capability with a parameter-estimation procedure and presents experimental verification through a proof-of-concept study. Described in greater detail in that article, the parameter estimator requires multiple forward simulations to analyze the data, which behooves minimizing the compute time for each simulation. By using an alternating direction implicit time-marching scheme and a structured but variably spaced grid, the numerical simulator built for this purpose provides a forward model output with suitable accuracy in approximately 0.5 seconds.

Energy & Fuels↗

Management of Risk and Uncertainty Through Optimized Co-Operation of Transmission Systems and Microgrids With Responsive Loads (Final Report)

The evolution of the power system to the reliable, efficient and sustainable system of the future will involve development of both demand- and supply-side technology and operations. Ambitious national and state-level goals around the decarbonization of electricity relies on the integration of very high levels of renewable resources, most of which are variable and intermittent. The use of demand response is an ideal approach to counterbalance the intermittency of renewable generation and brings the consumer into the spotlight. Though individual consumers are interconnected at the low-voltage distribution system, these resources are typically modeled as variables at the transmission network level. Demand-side participation cannot be leveraged effectively without explicitly including the distribution system dynamics in the optimization-based wholesale market operations. This project grew from a vision for co-optimized interaction of distribution systems, or microgrids, with the high-voltage transmission system. In this framework, microgrids encompass consumers, distributed renewables and storage. The energy management system of the lower voltage system (distribution or microgrid) can also sell (buy) excess (necessary) energy from the transmission system. Until recently, very little research had been conducted on the co-optimization of these two systems due to computational limitations. However, advances in computational capabilities, and the judicious use of decomposition methods and innovative approximation methods for high-dimension dynamic programming made this goal a viable objective for this project, leading to a fundamental shift in the ability to integrate and fully utilize demand-side resources. To this end, the modeling framework developed introduces a novel co-optimization framework, to include the operations of both the transmission and distribution systems (or microgrids) in operational decision making. This framework was used to analyze renewable and distributed generation along with responsive demand and to compare the capability of co-optimized systems to perform with higher levels of variable renewables. An ideal microgrid is defined as an electric entity capable of operating in both interconnected (with the high-voltage grid) and islanded mode. As such, the microgrid should incorporate generating units (traditional units and intermittent) and if needed, exchange power with the high-voltage grid. The interplay between the microgrid and high-voltage grid motivated the development of the co-optimization approach to ensure efficient performance of the interconnected network. Results show that the use of a bi-level optimization approach is an appropriate structure, capable of co-optimizing a transmission system with multiple distribution systems and microgrids. While increasing the number of connected systems provides increasing flexibility for renewables integration this can also the economic benefits to the low-voltage subsystems with each additional system connected. Comparison of a traditional single-level decision structure with the co-optimization approach illustrates a reduction in overall system cost under co-optimization, while specific cost allocations to transmission and distribution systems are changed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗