Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid impact analysis”

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 775 records · Page 43

Impact of Frequency-Watt Control on the Dynamics of a High DER Penetration Power System

Inverter-coupled generation is slowly displacing synchronous generators in the power system, resulting in reduced system inertia and primary frequency response, which makes the system susceptible to larger frequency deviations. Such deviations may result in load shedding or tripping of inverters, thus reducing the reliability. Frequency-watt control of inverters is an autonomous approach to support grid frequency at a fast timescale by adjusting the inverter output with changes in frequency. This study investigates the dynamics of a power system with a high penetration of distributed PV generation with frequency-watt droop control enabled. A simple aggregated reduced-order two-source model of the O'ahu Island (Hawaii) power system is developed and validated against a full dynamic model, and the small-signal stability of the reduced model is then examined the under various PV penetrations. The impact of frequency-watt parameters on the system frequency dynamics is investigated by observing trajectories of small- signal eigenvalues. The results of the eigenvalue analysis are validated through time-domain simulation. Considerations for frequency-watt parameter selection are discussed.

14 SOLAR ENERGY↗

The Technical and Economic Potential of the H2@Scale Hydrogen Concept within the United States

The U.S. energy system is evolving as society and technologies change. Renewable electricity generation - especially from wind and solar - is growing rapidly, and alternative energy sources are being developed and implemented across the residential, commercial, transportation, and industrial sectors to take advantage of their cost, security, and health benefits. Systemic changes present numerous challenges to grid resiliency and energy affordability, creating a need for synergistic solutions that satisfy multiple applications while yielding system-wide cost and emissions benefits. One such solution is an integrated hydrogen energy system. This is the focus of H2@Scale - a U.S. Department of Energy (DOE) initiative led by the Office of Energy Efficiency and Renewable Energy’s Hydrogen and Fuel Technologies Office. H2@Scale brings together stakeholders to advance affordable hydrogen production, transport, storage, and utilization in multiple energy sectors. The H2@Scale concept involves hydrogen as an energy intermediate. Hydrogen can be produced from various conventional and renewable energy sources including as a responsive load on the electric grid. Hydrogen has many current applications and many more potential applications, such as energy for transportation - used directly in fuel cell electric vehicles (FCEVs), as a feedstock for synthetic fuels, and to upgrade oil and biomass - feedstock for industry (e.g., for ammonia production, metals refining, and other end uses), heat for industry and buildings, and electricity storage. Owing to its flexibility and fungibility, a hydrogen intermediate could link energy sources that have surplus availability to markets that require energy or chemical feedstocks, benefiting both. This document builds upon a growing body of analyses of hydrogen as an energy intermediate by reporting the results from our initial analysis of the potential impacts of the H2@Scale vision by the mid-21 st century for the 48 contiguous U.S. states. Previous estimates have been based on expert elicitation and focused on hydrogen demands. We build upon them, first, by estimating hydrogen’s serviceable consumption potential for possible hydrogen applications and the technical potential for producing hydrogen from various resources. We define the serviceable consumption potential as the quantity of hydrogen that would be consumed to serve the portion of the market that could be captured without considering economics (i.e., if the price of hydrogen were $0/kg over an extended period); thus, it can be considered an upper bound for the size of the market. We define the technical potential as the resource potential constrained by real-world geography and system performance, but not by economics. We then compare the cumulative serviceable consumption potential with the technical potential of a number of possible sources. Second, we estimate economic potential: the quantity of hydrogen at an equilibrium price at which suppliers are willing to sell and consumers are willing to buy the same quantity of hydrogen. We believe this method provides a deeper understanding than was available in the previous analyses. We develop economic potentials for multiple scenarios across various market and technology-advancement assumptions.

08 HYDROGEN↗

The global gridded crop model intercomparison: Data and modeling protocols for Phase 1 (v1.0)

We present protocols and input data for Phase 1 of the Global Gridded Crop Model Intercomparison, a project of the Agricultural Model Intercomparison and Improvement Project (AgMIP). The project consist of global simulations of yields, phenologies, and many land-surface fluxes using 12–15 modeling groups for many crops, climate forcing data sets, and scenarios over the historical period from 1948 to 2012. The primary outcomes of the project include (1) a detailed comparison of the major differences and similarities among global models commonly used for large-scale climate impact assessment, (2) an evaluation of model and ensemble hindcasting skill, (3) quantification of key uncertainties from climate input data, model choice, and other sources, and (4) a multi-model analysis of the agricultural impacts of large-scale climate extremes from the historical record.

54 ENVIRONMENTAL SCIENCES↗

Numerical study of a delta planform with multiple jets in ground effect

The flow past a 60-deg delta wing equipped with two thrust-reverser jets near the inboard trailing edge has been analyzed by numerical solution of the 3D thin-layer Navier-Stokes equations. An implicit, partially flux-split, approximately-factored Navier-Stokes solver coupled with a multiple grid embedding scheme has been adapted to this problem. Studies of the impact of numerical parameters (e.g., grid refinement and dissipation levels), and flow-field parameters such as the height of the delta wing above the ground plane and the jet size on the solution, were performed. Results of these numerical studies indicate some challenges in the accurate resolution of complex 3D free shear layers and jets. Nevertheless, flow features such as jet deformation and ground vortex formation observed in experimental flow visualizations are captured. Further, comparisons with experimental data confirm the ability to simulate the loss of wing-borne lift, commonly referred to 'suckdown, as the delta planform flies at slow speeds in close proximity to the ground. Detailed analysis of the numerical results has also given additional insight into the structure of the ground vortex and the mechanisms of lift loss.

Chawla, K.↗

Load Flexibility Analysis of Residential HVAC and Water Heating and Commercial Refrigeration

The emergence of connected devices and the Internet of Things has opened the door to the concept of widespread use of building equipment for load-side flexibility. This stems from the increased use of distributed and clean energy resources for electric generation and the growing demand for a more resilient and reliable grid. A significant challenge associated with the use of building equipment as a distributed energy resource is quantifying the magnitude of power or energy flexibility. Unlike batteries that have well-defined and relatively constant capacity and maximum charge/discharge rates, building equipment often has characteristics that vary based on external factors (e.g., weather or occupant behavior). This makes quantifying the availability and flexibility of these resources difficult. Controlling a fleet of equipment while minimizing any adverse impacts to the occupants (e.g., discomfort or food spoilage) is a nontrivial matter. Although there are many different control methodologies available, this study focuses on two: priority-based control (PBC) and model-predictive control (MPC). These were selected because they represent two ends of the control spectrum. PBC is a relatively simple control strategy that requires only the current state of the devices to be controlled and an estimate of how much power each device consumes when turned on. MPC is much more complex and requires a model of the system and forecasts of any disturbances that influence the system. Because MPC has a predictive element, it can prepare for a desired response before the event and therefore is expected to provide a more optimal control relative to PBC. Both MPC and PBC are compared with a baseline simulation that uses feedback control typical of the thermostatically controlled devices being evaluated. This report investigates the load flexibility of three building systems: residential air conditioners, residential water heaters, and commercial refrigeration. Air heating and cooling and water heating account for 45% of residential electricity use, making them good candidates for providing load flexibility. Refrigeration in commercial buildings accounts for 16% of commercial building electricity use, making it another good candidate. Simplified models of each system were generated and validated based on measured data from actual equipment. These models were used to simulate the thermal and power response of the systems using PBC and MPC for a fleet of devices trying to achieve a variety of different responses, including hourly load shed, hourly load up, and daily peak demand reduction. Significant findings from each piece of equipment follow.

42 ENGINEERING↗

Impact of Wind, Solar, and Other Factors on Wholesale Power Prices: An Historical Analysis—2008 through 2017

Wholesale power markets have evolved. Some of the most prominent changes over the last decade in the United States include growth in wind and solar, a reduction in the price of natural gas, weakened load growth, and an increase in the retirement of thermal power plants. Here we empirically assess the degree to which wind and solar—among other factors—have influenced wholesale electricity prices. We show that wind and solar have contributed to reductions in overall average annual wholesale electricity prices since 2008, but that natural gas prices have had the largest impact. More notable is that expansion of variable renewable energy has led to significant changes in locational, time of day, and seasonal pricing patterns in some regions. These altered pricing patterns reflect a fundamental shift, and hold important implications for the grid-system value of wind and solar, and for other electric-sector planning and operating decisions.

14 SOLAR ENERGY↗

A bi-level advanced control framework for large-scale control of buildings with system-level impact

Increased electricity consumption combined with new forms of generation is testing the reliability of our grid infrastructure. This work describes a method to improve the reliability of the grid through large-scale advanced building control. This paper develops a bi-level distributed control framework to shift the load of 153 buildings to achieve a system-level objective of tracking a power reference signal. This bi-level control is based on the previously-developed ANPV-MPC, a predictive controller that uses a Bayesian neural network to generate an accurate control model and adapt to changing conditions over time. By shifting the building electricity demand to better match the available power, the grid system supplying the buildings is more reliable as evidenced by the analysis of node voltages across an IEEE 13-bus distribution system. The proposed bi-level control framework tracks the system-level power reference with enough accuracy to regulate node voltages across the IEEE 13-bus distribution system within ANSI limits of ±5%. Additionally, the adaptive nature of ANPV-MPC allows each building across the system to adapt to changing conditions, further amplifying the system-level reliability.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The climatological distribution of extreme Arctic winds and implications for ocean and sea ice processes

Abstract Some of the strongest near‐surface winds on Earth form in the Arctic and sub‐Arctic due to intense midlatitude cyclones and mesoscale processes, and these strong surface winds have important impacts on ocean and sea ice processes. We examine the climatological distribution of over‐ocean, near‐surface wind speeds within a Pan‐Arctic domain for 18 years (1990–2007) in four gridded data sets: the European Centre for Medium‐Range Weather Forecasts Interim reanalysis (ERA‐I), the Climate Forecast System Reanalysis, version 2 of the Common Ocean‐Ice Reference Experiment data set, and a regional climate simulation generated using the Weather Research and Forecasting (WRF) model run at 50 km (WRF50) horizontal resolution with ERA Interim as lateral boundary conditions. We estimate probability density functions, the annual cycle, and map the 50th and 99th percentile winds. We then perform the same statistical analysis of winds for 2 years when 10 km WRF data are available (June 2005 to May 2007); despite the much shorter time period, the Pan‐Arctic statistics are very similar to those from the 18 year analysis. We repeat the wind speed statistical analysis within a subdomain surrounding Greenland and find that WRF10 has consistently larger maximum wind speeds, but this difference only appears at wind speed percentiles higher than 99%. Differences in the 99th percentile wind speeds are spatially heterogeneous. An investigation of surface fluxes within WRF50 and WRF10 reveals unrealistically large sensible heat fluxes along the sea ice edge, and the geographic distribution and magnitude of these fluxes is shown to be sensitive to sea ice representation in WRF.

Hughes, Mimi↗

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↗

Global land use for 2015–2100 at 0.05° resolution under diverse socioeconomic and climate scenarios

Abstract Global future land use (LU) is an important input for Earth system models for projecting Earth system dynamics and is critical for many modeling studies on future global change. Here we generated a new global gridded LU dataset using the Global Change Analysis Model (GCAM) and a land use spatial downscaling model, named Demeter, under the five Shared Socioeconomic Pathways (SSPs) and four Representative Concentration Pathways (RCPs) scenarios. Compared to existing similar datasets, the presented dataset has a higher spatial resolution (0.05° × 0.05°) and spreads under a more comprehensive set of SSP-RCP scenarios (in total 15 scenarios), and considers uncertainties from the forcing climates. We compared our dataset with the Land Use Harmonization version 2 (LUH2) dataset and found our results are in general spatially consistent with LUH2. The presented dataset will be useful for global Earth system modeling studies, especially for the analysis of the impacts of land use and land cover change and socioeconomics, as well as the characterizing the uncertainties associated with these impacts.

54 ENVIRONMENTAL SCIENCES↗

Hy-DAT: A Tool to Address Hydropower Modeling Gaps Using Interdependency, Efficiency Curves, and Unit Dispatch Models

As the power system continues to be flooded with intermittent resources, it becomes more important to accurately assess the role of hydro and its impact on the power grid. While hydropower generation has been studied for decades, dependency of power generation on water availability and constraints in hydro operation are not well represented in power system models used in the planning and operation of large-scale interconnection studies. There are still multiple modeling gaps that need to be addressed; if not, they can lead to inaccurate operation and planning reliability studies, and consequently to unintentional load shedding or even blackouts. As a result, it is very important thathydropower is represented correctly in both steady-state and dynamic power system studies. In this paper, we discuss the development and use of the Hydrological Dispatch and Analysis Tool (Hy-DAT) as an interactive graphical user interface (GUI), that uses a novel methodology to address the hydropower modeling gaps like water availability and inter-dependency using a database, and algorithms, to generate accurate representative models for power system simulation.

Wang, Dewei↗

Integrated Security-Safety Consequence Study of a Heat Pipe Reactor

This report presents a preliminary demonstration of integrated security/safety consequence modeling to support potential physical security exemption justifications under the new Title 10 Code of Federal Regulations Part 53 licensing pathway, which allows exemptions if sabotage-induced radiological consequences do not exceed 25 rem at the site boundary two hours after release without mitigative actions. A hypothetical heat pipe microreactor subjected to direct sabotage by a 150 lbs. TNT-equivalent high-explosive device was analyzed using a simplified radionuclide release estimate based on DOE-HDBK-3010-94, implemented in the MELCOR severe accident analysis code to model release transport and attenuation through reactor/building compartments, with MELCOR outputs coupled to MACCS to compute offsite doses and the distance-to-threshold metric D 25rem . Six scenarios were evaluated to examine sensitivities to release pathway/building configuration and radionuclide form (vapor versus aerosols with different size distributions). Results show a positive correlation between leak path factor and D 25rem and indicate that reactor building compartmentalization can substantially reduce consequences (a single-story configuration reduced D 25rem by roughly an order of magnitude relative to a cavity-only configuration), while vapor versus aerosol assumptions had limited impact for the modeled two-story case on the spatial grid used.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Global land use for 2015-2100 at 0.05° resolution under diverse socioeconomic and climate scenarios

Global future land use (LU) is an important input for Earth system models for projecting Earth system dynamics and is critical for many modeling studies on future global change. Here we generated a new global gridded LU dataset using the Global Change Analysis Model (GCAM) and a land use spatial downscaling model, named Demeter, under the five Shared Socioeconomic Pathways (SSPs) and four Representative Concentration Pathways (RCPs) scenarios. Compared to existing similar datasets, the presented dataset has a higher spatial resolution (0.05°×0.05°) and spreads under a more comprehensive set of SSP-RCP scenarios (in total 15 scenarios), and considers uncertainties from the forcing climates. We compared our dataset with the Land Use Harmonization version 2 (LUH2) dataset and found our results are in general spatially consistent with LUH2. The presented dataset will be useful for global Earth system modeling studies, especially for the analysis of the impacts of land use and land cover change and socioeconomics, as well as the characterizing the uncertainties associated with these impacts.

54 ENVIRONMENTAL SCIENCES↗

New Approaches to Quantifying Transport Model Error in Atmospheric CO2 Simulations

In recent years, much progress has been made in observing CO2 distributions from space. However, the use of these observations to infer source/sink distributions in inversion studies continues to be complicated by difficulty in quantifying atmospheric transport model errors. We will present results from several different experiments designed to quantify different aspects of transport error using the Goddard Earth Observing System, Version 5 (GEOS-5) Atmospheric General Circulation Model (AGCM). In the first set of experiments, an ensemble of simulations is constructed using perturbations to parameters in the model s moist physics and turbulence parameterizations that control sub-grid scale transport of trace gases. Analysis of the ensemble spread and scales of temporal and spatial variability among the simulations allows insight into how parameterized, small-scale transport processes influence simulated CO2 distributions. In the second set of experiments, atmospheric tracers representing model error are constructed using observation minus analysis statistics from NASA's Modern-Era Retrospective Analysis for Research and Applications (MERRA). The goal of these simulations is to understand how errors in large scale dynamics are distributed, and how they propagate in space and time, affecting trace gas distributions. These simulations will also be compared to results from NASA's Carbon Monitoring System Flux Pilot Project that quantified the impact of uncertainty in satellite constrained CO2 flux estimates on atmospheric mixing ratios to assess the major factors governing uncertainty in global and regional trace gas distributions.

Ott, L.↗

Risk-Based Dynamic Contingency Analysis Applied to Puerto Rico Electric Infrastructure

The U.S. Department of Energy’s (DOE) Office of Electricity and Office of Energy Efficiency & Renewable Energy have funded DOE National Laboratories to perform modeling, analysis, and high-level design of resilience-enhancement options for the power grid of the Commonwealth of Puerto Rico. The Pacific Northwest National Laboratory (PNNL) is one of the national laboratories contributing to the DOE effort. Under this funding, PNNL completed Phase I of their analysis in 2018, which identified high-priority transmission enhancements derived from detailed dynamic cascading analysis of severe contingencies, including a hurricane scenario example. This report describes additional analysis completed under Phase II, which was performed over 2019. The PNNL team applied decades of experience making complex power systems more resilient, reliable, secure, flexible, affordable and sustainable, with partners in government and industry. In this report, PNNL presents the Phase II analysis, in which a risk-based dynamic contingency analysis approach to evaluate impact of several hurricane scenarios was developed. This approach was used to identify high-priority enhancements and test resilience mitigation actions, including the evaluation of high-solar scenarios developed as a result of the 2019 Puerto Rico Integrated Resource Plan (IRP).

24 POWER TRANSMISSION AND DISTRIBUTION↗

Drought and Deluge: the Recurrence of Hydroclimate Extremes During the Past 600 Years in Eastern Australia’s Natural Resource Management (NRM) Clusters

Recent extremes of flood and drought across Australia have raised questions about the recurrence of such rare events and highlighted the importance of understanding multi-decadal climate variability. However, instrumental records over the past century are too short to adequately characterise climate variability on multi-decadal and longer timescales or robustly determine extreme event frequencies and their duration. Palaeoclimate reconstructions can provide much-needed information to help address this problem. Here, we use the 600-year hydroclimate record captured in the eastern Australian and New Zealand Drought Atlas (ANZDA) to analyse drought and pluvial frequency trends for East Australian Natural Resource Management (NRM) clusters. This partitioning of the drought atlas grid points into recognised biophysical areas (i.e. NRM clusters) revealed their differences and similarities in drought intensity and pluvial events over time. We find sustained multi-decadal periods of a wet–dry geographic 'seesaw' between eastern to central and southern NRMs (e.g. 1550–1600 CE and 1700–1750 CE). In contrast, other periods reveal spatially consistent wetting (e.g. 1500–1550 CE) or drying (e.g. 1750–1800 CE). Emerging hot spot analysis further shows that some areas that appear naturally buffered from severe drought during the instrumental period have a greater exposure risk when the longer 600-year record is considered. These findings are particularly relevant to management plans when dealing with the impacts of climate extremes developed at regional scales. Our results demonstrate that integrating and extending instrumental records with palaeoclimate datasets will become increasingly important for developing robust and locally specific extreme event frequency information for managing water resources.

Australian hydroclimate↗

Final Report for Project: Impacts of stratification and non-equilibrium winds and waves on hub-height winds

This project used a combination of turbulence-resolving large-eddy simulations, single-column modeling (where turbulence is parameterized), and currently available observations to improve, assess, and develop a parameterization of the impact of non-equilibrium wave states and stratification on the buoy-observed winds to establish reliable wind data at the turbine hub-height level. Analysis of turbulence-resolving simulations and observations illuminates the non-linear coupling between the atmosphere and the undulating sea surface. This analysis guides modification of existing boundary layer parameterizations to include wave influences for upward extrapolation of surface-based observations through the turbine layer. Our surface roughness modifications account for the interaction between stratification and the effects of swell’s amplitude and wavelength as well as swell’s relative motion with respect to the mean wind direction. The single-column version of the open source Weather and Research Forecasting (WRF) model (Skamarock et al., 2008) serves as our platform to test our proposed planetary boundary layer parameterization modifications that account for wave effects on marine atmospheric boundary layer flows. WRF has been widely adopted for wind resource analysis and forecasting. The single column version is particularly suitable to development, analysis, and testing of new boundary layer parameterizations. We utilize WRF’s single-column version to verify and validate our proposed modifications to the Mellor-Yamada-Nakanishi-Niino (MYNN) boundary layer parameterization (Nakanishi and Niino, 2004). We explore the implications of our modifications for two-way coupling between WRF and wave models (e.g.,Wavewatch III). The newly implemented parameterization accounting for marine atmospheric boundary layer-wave coupling is then tested in three-dimensional WRF simulations at grid sizes near 1 km. These simulations identify the behavior of simulated winds at the wind plant scale. Overall project conclusions include; In the presence of fast-moving swell (significant wave height Hs = 6.4 m, and phase speed cp = 18 ms -1 ), the atmospheric boundary layer grows more rapidly when waves propagate opposite to the winds compared to when winds and waves are aligned. Pressure drag increases by nearly a factor of 2 relative to the turbulent stress for the extreme case where waves propagate at 180° compared to the pressure gradient forcing. Net wind speed reduces by nearly 15% at hub-height for the 180°-case compared to the 0°-case, and turbulence intensities increase by nearly a factor of 2. These impacts diminish with decreasing wave age; Stratification increases hub height wind speeds and increases the vertical shear of the mean wind across the rotor plane. Fortuitously, this stability-induced enhanced shear does not influence turbulence intensity at hub height, but does increase (decrease) turbulence intensity below (above) hub height. Increased stability also increases the wave-induced pressure stress by ~ 10%; Off the East Coast of the United States during Coupled Boundary Layers Air-Sea Transfer - Low Wind (CBLAST-Low), cases with short fetch include thin stable boundary layers with depths of only a few tens of meters. In the coastal zone, the relationship between the mean wind and the surface fiction velocity (u*(V )) is significantly related to wind direction for weak winds but is not systematically related to the air sea difference of virtual potential temperature, δθv; since waves generally propagate from the south at the Air-Sea Interaction Tower (ASIT) tower, these results suggest that under weak wind conditions waves likely influence surface stress more than stratification does; and Winds and waves are frequently misaligned in the coastal zone. Stability conditions persist for long duration. Over a four year period, the Forschungsplattformen in Nord- und Ostsee Nr. 1 (FINO1) tower (a site with long fetch) primarily experienced weakly-unstable conditions, while stability at the ASIT tower (with a larger influence of offshore winds) experiences a mix of both unstable and stable conditions, where the summer months are predominantly stable. Wind-wave misalignment likely explains the large scatter in observed non-dimensional surface roughness under swell-dominated conditions. Andreas et al.’s (2012) relationship between u* and the 10-m wind speed under predicts the increased u* produced by wave-induced pressure drag produced by misaligned winds and waves. Incorporating wave-state (speed and direction) influences in parameterizations improves predictive skill. In a broad sense, these results suggest that one needs information on winds, temperature, and wave state to upscale buoy measurements to hub-height and across the rotor plane. Our parameterization of wave-state influences on surface drag has been submitted for inclusion in the next publicly available release. In combination, our project elucidates the impacts of two important physical processes (non-equilibrium wind/waves and stratification) on the atmosphere within which offshore turbines operate. This knowledge should help guide and inform manufacturers making critical decisions surrounding design criteria of future turbines to be deployed in the coastal zone. Reductions in annually averaged hub height wind speed error using our new wave-state-aware surface layer parameterization are relatively modest. However since wind turbine power production depends on the wind speed cubed, the error in estimated power production is close to 5%; which is significant and can substantially impact wind resource assessment and decision making with regards to the viability of particular location for a wind plant location. For a single 30-hour forecast, significant reductions in wind speed prediction errors can yield substantially improved wind power forecast skill, thereby mitigating costs and/or increasing revenue through improved; forecasting for maintenance operations and planning; day-ahead forecasting for power trading and resource allocation; and short-term forecasting for dispatch and grid balancing.

17 WIND ENERGY↗

Aerothermal Analysis of the Rocket Lab Venus Probe Heatshield

This document provides an overview of the aerothermodynamic analyses performed by the Aerothermodynamics Branch at NASA Langley Research Center for the Rocket Lab Venus Probe (RLVP). In addition to defining the baseline heating environment to the heatshield, this document pursues the experimental validation of key physical models at RLVP-relevant conditions. This experimental validation analysis, which captures the model form uncertainty, is used as one of two primary components of the margin assessment, where the other component is the parametric uncertainty. These model form (experimental) and parametric uncertainty components are used to construct a spatial and time varying margin for the heating to the RLVP heatshield. The margin is evaluated as the sum of the parametric and model form uncertainty components. The model form uncertainty is defined as the difference between the RLVP-relevant measurements and their simulations, using the upper limit uncertainty bounds for both the measurements and simulations in the comparisons. The differences in the dominant physics in the stagnation region and flank lead to the separate RLVP-relevant measurements for assessing the model form uncertainty in these two regions. These regions are addressed as follows: Stagnation Region Heating Environment: For the high-temperature stagnation-region, both the radiative heating and impact of blowing on convective heating are significant, while the impacts of turbulence and roughness are negligible. Coupled radiation and ablation LAURA/HARA solutions with ray-tracing provide the radiative heating over the entire vehicle, including the contributions from the Venus atmosphere and ablation species. Non-ablating LAURA simulations provide the convective heating. During the material-response computation typically used for TPS sizing, this non-ablating convective heating is corrected for the impact of ablation using the blowing correction. Coupled ablation LAURA simulations that capture finite-rate sur-face processes show that this blowing correction may be non-conservative over most of the heatshield. This non-conservatism is due to hydrogen recombination in the finite-rate surface model, which tends to increase the coupled ablation convective heating to near the non-ablating values, therefore making any reduction in the non-ablating value through the blowing correction non-conservative. This non-conservatism due to H catalysis is captured in the parametric component of the margin. The best available ground-test measurements that capture the impact of blowing on stagnation region convective heating, at RLVP-relevant conditions, indicate that the current blowing reduction model is non-conservative by up to 20% at RLVP-relevant blowing rates (the coupled blowing simulations were also non-conservative). Because of the relatively low velocity of the ground tests and the non-Venus atmospheric chemistry, these measurements do not capture the chemistry and therefore do not inform the uncertainty due to H catalysis. However, they do capture the fluid mechanics of blowing. The non-conservatism of the blowing correction implied by these measurements is covered by the model form component of the margin, which leads to total margin values over 50%. For the radiative heating, the shock-tube informed bias approach suggests a model form uncertainty of roughly 20%, while the parametric uncertainty analysis suggests values over 100%. The combined stagnation-point radiation margin of over 100% leads to peak margined radiative heating values of over300 W/cm2, which remains small relative to the peak margined convective heating of nearly 2000 W/cm2. Based on this analysis, at the stagnation point, the peak margined heat rate is 2203 W/cm2 and the margined total heat load is 31.5 kJ/cm2 for the current nominal trajectory. Flank Heating Environment: The forebody flank (and near-shoulder) heating environment is dominated by the impact of turbulence, roughness augmentation, and ablation on the convective heating. An extensive collection of ground test measurements with RLVP-relevant turbulence and roughness is studied to show that the maximum difference between the simulated and measured convective heating is 5%. However, with the exception of the Holden measurements from the 1980s, these measurements do not include roughness elements extending into the supersonic region of the boundary layer, which is likely to occur for RLVP (due to the 45 degree sphere-cone geometry). The interaction between the supersonic flow and roughness could cause convective heating augmentation fundamentally different than for locally subsonic flow. Although these Holden measurements are consistent with the other measurements considered, another path was pursued to assure that the rough-ness height extending into supersonic flow does not fundamentally change the roughness augmentation. This additional path was a computational effort to resolve the roughness elements in the CFD grid, so that the interaction be-tween the roughness elements and locally supersonic flow may be simulated in detail. This roughness-resolved CFD simulation is feasible because of the RLVP forebody TPS’s patterned roughness, which may be approximated analytically, and because of the axisymmetric nominal flow field, which allows a narrow surface region to be simulated and therefore make the computational expense feasible. These grid-resolved roughness simulations, which are performed at actual RLVP flight conditions, result in heating augmentation values that are below the design approach for roughness augmentation. This provides evidence that the design approach for RLVP roughness augmentation is sufficient. Based on this analysis, at this flank or near-shoulder location, the peak margined total heat rate is 2088 W/cm2and the margined total heat load is 26.0 kJ/cm2for the current nominal trajectory. Heat flux, shear, pressure and heat transfer coefficient at the RLVP stagnation point and near shoulder location are evaluated for the entire trajectory, and curved fit to a functional form of F=AρB∞UC∞. These simplified relationships for the nominal and margined aerothermal environments are referred to as aerothermal indicators, and presented at the end of this document.

Christopher O Johnston↗