Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wind direction variability”

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.

552 records · Page 31

Data and scripts from: “Denoising autoencoder for reconstructing sensor observation data and predicting evapotranspiration: noisy and missing values repair and uncertainty quantification”

This data package includes data and scripts from the manuscript “Denoising autoencoder for reconstructing sensor observation data and predicting evapotranspiration: noisy and missing values repair and uncertainty quantification”.The study addressed common challenges faced in environmental sensing and modeling, including uncertain input data, missing sensor observations, and high-dimensional datasets with interrelated but redundant variables. Point-scaled meteorological and soil sensor observations were perturbed with noises and missing values, and denoising autoencoder (DAE) neural networks were developed to reconstruct the perturbed data and further predict evapotranspiration. This study concluded that (1) the reconstruction quality of each variable depends on its cross-correlation and alignment to the underlying data structure, (2) uncertainties from the models were overall stronger than those from the data corruption, and (3) there was a tradeoff between reducing bias and reducing variance when evaluating the uncertainty of the machine learning models.This package includes:(1) Four ipython scripts (.ipynb): “DAE_train.ipynb” trains and evaluates DAE neural networks, “DAE_predict.ipynb” makes predictions from the trained DAE models, “ET_train.ipynb” trains and evaluates ET prediction neural networks, and “ET_predict.ipynb” makes predictions from trained ET models.(2) One python file (.py): “methods.py” includes all user-defined functions and python codes used in the ipython scripts.(3) A “sub_models” folder that includes five trained DAE neural networks (in pytorch format, .pt), which could be used to ingest input data before being fed to the downstream ET models in ‘ET_train.ipynb” or ‘ET_predict.ipynb’.(4) Two data files (.csv). Daily meteorological, vegetation, and soil data is in “df_data.csv”, where “df_meta.csv” contains the location and time information of “df_data.csv”. Each row (index) in “df_meta.csv” corresponds to each row in “df_data.csv”. These data files are formatted to follow the data structure requirements and be directly used in the ipython scripts, and they have been shuffled chronologically to train machine learning models. The meteorological and soil data was collected using point sensors between 2019-2023 at(4.a) Three shrub-dominated field sites in East River, Colorado (named “ph1”, “ph2” and “sg5” in “df_meta.csv”, where “ph1” and “ph2” were located at PumpHouse Hillslopes, and “sg5” was at Snodgrass Mountain meadow) and(4.b) One outdoor, mesoscale, and herbaceous-dominated experiment in Berkeley, California (named “tb” in “df_meta.csv”, short for Smartsoils Testbed at Lawrence Berkeley National Lab).- See "df_data_dd.csv" and "df_meta_dd.csv" for variable descriptions and the Methods section for additional data processing steps. See "flmd.csv" and "README.txt" for brief file descriptions.- All ipython scripts and python files are written in and require PYTHON language software.

54 ENVIRONMENTAL SCIENCES↗

VALENTInE: A Concept for a New Frontiers Class Long Duration In-Situ Balloon Mission to Venus

Venus and Earth are similar in bulk composition, size, density, and approximate distance from the Sun, yet Venus’s modern-day climate and surface geology is distinctly different [1]. Previous missions to Venus revealed unusual volcanic features, possible continental crust, widespread volcanic plains, a weak magnetic field [2], and insights into the Venusian atmosphere [3,4]. Unfortunately, data from these missions were limited in spatial and temporal resolution and global extent. A future mission to Venus is critical to address fundamental questions surrounding the chemical composition and dynamics of the Venusian atmosphere [5], its geologic history [6,7], its internal structure [8,9], and its habitability throughout time [10,11]. We present the Venus Air and Land Expedition: a Novel Trailblazer for In situ Exploration (VALENTInE) mission to meet this need. VALENTInE is a variable altitude balloon that will passively float in Venus’s atmosphere between 45 and 55 km altitude. VALENTInE will acquire atmospheric data at varying latitudes and longitudes in addition to mapping the surface geomorphology and mineralogy across multiple terrains. Mission Objectives: The VALENTInE mission concept is driven by four main science objectives: 1. Determine whether the driving force of the superrotation of Venus’s atmosphere is caused by horizontal or vertical momentum transport. 2. Determine whether the atmospheric composition and noble gas inventory of the Venusian atmosphere is a product of outgassing from the initial protoplanetary source or if there are significant contributions from exogenic sources. 3. Determine whether the tesserae regions (particularly Aphrodite Terra) are felsic and relatively older than surrounding regions. 4. Determine if there is any evidence of a recent dynamo preserved in the rock record of Venus. Mission Overview: A balloon architecture provides a robust system that can survive long-term in the Venusian environment while taking accurate measurements of the lower cloud deck and surface. Two prior balloon missions to Venus, VEGA 1 and VEGA 2 in 1984, have demonstrated the potential for such planetary exploration; however, these missions were short-lived (46 hr), with limited range (54 km). Our spacecraft concept consists of a battery powered balloon with a gondola and flyby carrier stage (Fig. 1). Our novel design allows for at least 15 days of atmospheric exploration, including multiple ascents and descents in the Venus atmosphere. VALENTInE is baselined to launch in 2032, cruise on solar power for 128 days, enter, descend, and inflate (EDI) into Venus’s atmosphere above Aphrodite Terra, and float in Venus’s atmosphere for a total nominal mission duration of 15 days. The balloon, able to control its altitude by changing its buoyancy, will make one full circumnavigation every 4–8 days and move poleward ~1 o latitude per day; it will be passively directed by the horizontal air currents on Venus. The balloon can only be controlled in the z direction, and the expected latitude range is ±10° from EDI. The balloon itself will be a tracer for atmospheric structures such as zonal winds [5, 12]. During the 15-day operational period, there will be five dips to 45 km for a compositional study of the lower atmosphere and geological and magnetic mapping of the surface. Dips are used to take images and measurements closer to the surface and to obtain in situ vertical atmospheric profiles between 45 and 55 km. Instruments: The VALENTInE instrument payload will allow for extensive study of the geology, atmosphere, and interior of Venus. The instrument suite consists of six instruments. The mission profile for each instrument is shown in Fig. 2. Lower Atmosphere Meteorology Analyzer (LAMA) is an atmospheric structure instrument consisting of a thermometer, barometer, and accelerometer to continuously measure temperature and pressure as a function of altitude, longitude, and latitude. TracE and Noble Gas Investigator (TENGI) is aneutral mass spectrometer used to sample the dense atmosphere. TENGI will operate continuously at 45 km and 55 km and will measure D/H ratios, as well as Ne, Ar, and O isotopic ratios. Kilometer Scale Spectral Imager (KSSI) is amultispectral imager (850-1150 nm) used to resolve surface features on the order of km. Near InfraRed Multispectral Photometer (NIRMP) is an Infrared (IR) photometer used to image the surface through the clouds to characterize mineral assemblages at five different areas of the surface. NIRMP will be operated during dwell. ELevation REConnaissance (ELREC) is a light detection and ranging (LIDAR) instrument used for measuring topography at five different areas of the surface. The ELREC data will be combined with those of NIRMP to determine mineral assemblages and how they correlate to topography. Magnetic Exploration and Interior Detective (MEID) is a magnetometer used to detect any near-surface magnetic anomalies. The magnetometer will continuously operate at all altitudes. Mission Design: The spacecraft’s propulsion system will have a wet mass of 930 kg and will be launched on an Atlas V rocket. The spacecraft will jettison the payload upon arrival in Venus’s atmosphere before decelerating to orbital velocity. Figure 2. Notional power-dip profile while the balloon is on the dayside. Risks and Challenges: Raising and lowering the spacecraft to dip beneath the haze layer require a large amount of energy. This issue was partly compensated by limiting the number of dips over the 15 day prime mission. Power needs for pumping helium are reduced by slowing the descent speed, which consequently increases spatial resolution of in situ measurements between 45 and 55 km. The main limitation for the duration of the mission is the need to carry 15 days’ worth of batteries. The bus experiences external Venusian temperatures ranging from ~27°C at 55 km altitude to ~110°C at 45 km altitude. These high atmospheric temperatures require thermal protection for the bus (instruments, electronics, and flight systems) during its descent, dwell at 45 km, and ascent. The bus is maintained at a mechanically safe temperature range of -10°C to 50°C using white external paint, multi-layer insulation (MLI), mechanical/thermal isolation (e.g. Ti, composites), and ~57 kg of Phase Change Materials (PCM). Images and spectra taken below the haze dominate the available data transmission regardless of the time spent at 45 km. Therefore, less time at the lower altitude mainly reduces the coverage of in situ measurements there. Further, we assumed an unrealistically low-density material for the helium storage tanks on the gondola. However, this issue is partially resolved if we were to jettison 75% of the spent helium storage tanks. After the entry process, the balloon volume remains mostly inflated and storage tanks are required only for reducing balloon volume in dipping to 45 km. This mission was planned against an expected New Frontiers 5 (NF5) cost cap of $1B, as the NF5 call had not yet been released. The mission we describe fits within the predicted cost cap if the flyby carrier stage, responsible for powering the spacecraft during cruise and for relaying in situ measurements back to Earth, can be contributed by another space agency. Acknowledgments: We thank the JPL Planetary Science Summer School, especially A. Nash, J. Scully, K. Mitchell, L. Lowes, and J. Armijo, and our mentors from JPL Team X. Thank you to our review panel for their insightful review. References: [1] Kane, S. et al. (2019) JGR:Ps, 124 , 2015–2028. [2] Phillips, J. L., & Russell, C. T. (1987). JGR: Space Physics , 92 (A3), 2253-2263. [3] Nakamura, M. et al. (2018) Earth Planets Space, 70 (1), 144. [4] Svedhem, H. et al. (2009) JGR:P 114 (E5). [5] Limaye et al. (2009) Decadal Survey White Paper . [6] Ivanov, M. and Head, J. (2011) PSS, 59 (13) , 1559-1600. [7] Smrekar, S. et al. (2018) Space Sci Rev, 214 (5), 88. [8] O’Rourke, J. et al. (2018) EPSL, 502 , 46-56. [9] O’Rourke J. et al. (2019) GRL , (46), 5768–5777 [10] Way, M. et al. (2016) GRL, 43 (16), 8376-8383. [11] Way, M. et al. (2020) JGR:P, 125 (5), e2019JE006276. [12] Preston, R. A. et al. (1986) Science, 231 (4744), 1414-1416.

Mission Concept↗

Application of Remote Sensing to Assess the Impact of Short Term Climate Variability on Coastal Sedimentation

The purpose of this joint University of Wisconsin (UW) and Louisiana State University (LSU) project has been to relate short term climate variation to response in the coastal zone of Louisiana in an attempt to better understand how the coastal zone is shaped by climate variation. Climate variation in this case largely refers to variation in surface wind conditions that affect wave action and water currents in the coastal zone. The primary region of focus was the Atchafalaya Bay and surrounding bays in the central coastal region of Louisiana. Suspended solids in the water column show response to wind systems both in quantity (through resuspension) and in the pattern of dispersement or transport. Wind systems associated with cold fronts are influenced by short term climate variation. Wind energy was used as the primary signature of climate variation in this study because winds are a significant influence on sediment transport in the micro-tidal Gilf of Mexico coastal zone. Using case studies, the project has been able to investigate the influence of short term climate variation on sediment transport. Wind energy data, collected daily for National Weather Service (NWS) stations at Lake Charles and New Orleans, LA, were used as an indicator of short term climate variation influence on seasonal time scales. A goal was to relate wind energy to coastal impact through sediment transport. This goal was partially accomplished by combining remote sensing and wind energy data. Daily high resolution remote sensing observations are needed to monitor the complex coastal zone environment, where winds, tides, and water level all interact to influence sediment transport. The NASA Earth Observing System (EOS) era brings hope for documenting and revealing response of the complex coastal transport mosaic through regular high spatial resolution observations from the Moderate resolution Imaging Spectrometer (MODIS) instrument. MODIS observations were sampled in this project for information content and should continue to be viewed as a resource for coastal zone monitoring. The project initialized the effort to transfer a suspended sediment concentration (SSC) algorithm to the MODIS platform for case 2 waters. MODIS enables monitoring of turbid coastal zones around the globe. The MODIS SSC algorithm requires refinements in the atmospheric aerosol contribution, sun glint influence, and designation of the sediment inherent optical properties (IOPs); the framework for continued development is in place with a plan to release the algorithm to the MODIS direct broadcast community.

Menzel, W. Paul↗

Automated point dendrometer, soil moisture and temperature, and meteorological variables datasets, Oct 2024 – Nov 2025, G.A. Pearson Natural Area, Flagstaff, AZ, USA

This data package includes parsed, cleaned, and calibrated data from 48 TOMST automated point dendrometers, 48 TOMST 15 cm soil moisture sensors, and 12 TOMST 30 cm soil moisture sensors. The point dendrometers were cleaned with the “dendRoAnalyst” package in RStudio. The soil sensors were cleaned and calibrated for volumetric water content (VWC) with the “myClim” package in RStudio using the soil texture of the site (sandy clay loam). Additionally, this data package also includes raw data from 2 METER weather stations. Dendrometers and soil sensors have both their sensor ID, as well as the ID for the specific tree they were instrumented on at the G.A. Pearson Natural Area (GPNA) site and their experimental group. The purpose of these data is to understand how ponderosa pine trees in restored (thinned and burned) vs. unrestored (no treatment) areas are responding to drought and seasonal precipitation. These data use radial growth and soil moisture data to answer the following question: how are active season length, growth on different time scales (weekly, monthly, seasonally, and annually), growth during dry periods and after precipitation events, and environmental and biological drivers of radial growth different between restored versus unrestored areas?

Air temperature↗

Chapter 10 - Remote Sensing Measurements of Aerosol Properties

Satellite instruments have proven especially capable at monitoring the quantity of airborne particles in columns of atmosphere, globally. This chapter describes the principles of satellite measurements and retrieval algorithms, and surveys current instruments and their capabilities. We outline the issues associated with retrieval algorithms, such as surface characterization and aerosol proximity to clouds, and the challenges with interpretation of the results. The relationship between measured aerosol properties and climate-relevant aerosol properties simulated in models is outlined, as well as how measurements are used to evaluate models. Most space-based aerosol instruments are passive sensors that measure reflected sunlight at multiple wavelengths, some at multiple viewing angles. A few are active sensors that send out their own laser light and measure the returned signal. Except when clouds are present, the excess amount of light scattered back to space, beyond that expected from the surface and atmospheric gas, is attributed to aerosol. Satellite measurements are used in many ways in aerosol research. They often provide the only method for monitoring hazardous phenomena such as major wildfire and volcanic eruption plumes, especially in remote areas. Stable, long-term, near-global-scale satellite data records make it possible to identify regional and global aerosol trends. Aerosol radiative effects on climate can be quantified on a near-global scale and used to estimate the strength of aerosol–radiation and aerosol–cloud interactions as well as to evaluate climate model simulations of these interactions. Aerosol-type mapping from satellite imagery is helpful for source attribution, model validation, and to constrain particle light-absorption properties that are essential for radiative forcing calculations. The range of aerosol properties retrieved from satellite observations has grown considerably since the first global estimates of aerosol optical depth (τ a) over ocean were made in the late 1970s. Methods for retrieving particle size and light-absorption properties were explored in the 1990s using multispectral, multi-angle observations, and polarization in visible and near-infrared wavelengths. Sensitivity to particle light absorption, primarily from black or brown carbon content, improved with the inclusion of UV channels, and sensitivity to very thin aerosol layers in the upper troposphere and lower stratosphere was advanced with the use of limb-sounding instruments and active sensors. There are limitations to every measurement technique, including satellite aerosol remote sensing. For wide-swath, passive instruments, aerosol retrievals near clouds can present substantial challenges as far as 15 km away due to cloud-scattered light contaminating the signal. In nearly all cases, retrievals over bright snow and ice surfaces are precluded because surface reflectance uncertainties can overwhelm the aerosol signal. Similarly, meteorological cloud is identified and masked out where possible. Data from passive sensors also lack vertical resolution except those that view toward the limb or where multi-angle imagery is acquired over plumes from wildfires, erupting volcanoes, and wind-blown dust. Yet, passive sensors provide vastly more coverage than the active instruments that mitigate these issues. Particle microphysical information is qualitative from all remote sensing techniques, relying on proxies to infer particle composition, hygroscopicity, and the amount of light-absorbing material. Further, particles smaller than about 200 nm diameter cannot be distinguished from atmospheric gas molecules with remote sensing, which hinders studies of cloud condensation nuclei and their effects on clouds. Most satellite instruments dedicated to aerosol observations are in low-Earth, near-polar, sun-synchronous orbits, which means they cross the equator at the same local time each day. Most are set on cycles that repeat approximately every 16 days, which makes it difficult to monitor aerosol evolution locally. Geostationary satellites make it possible to observe changes occurring from minutes to hours over regions up to 8000 km in size, but lack coverage of high latitudes, and often provide more limited constraints on aerosol properties. Ground-truth data are vital for satellite aerosol-retrieval validation. The AErosol RObotic NETwork (AERONET) of sun photometers was created in 1993 and has become an established global network of over 350 instruments for validating satellite measurements. The network, as well as global networks of ground-based lidars, solar flux radiometers and other sun photometers, are widely used for evaluating global satellite retrievals and model simulations. NASA's Earth Observing System (EOS) program beginning in 1999 led to improvements in reliability, spatial resolution, and spectral resolution (and hence, to improved particle size discrimination and light absorption properties). Satellite payloads include advanced broad-swath and multi-angle imagers, along with the first space-based active sensor focused largely on long-term aerosol monitoring. Since about 2002, Europe's SENTINEL and operational meteorological satellite fleets are also providing sustained aerosol observations, with planned continuation until at least 2030. Satellite remote sensing instruments offer valuable data for evaluating aerosol representations in global climate models. They have been used to assess aerosol optical and physical properties, trends and distributions, and are applied increasingly as direct model constraints in data assimilation to create global aerosol reanalysis products. Aerosol optical depth is the most common quantity adopted for routine model evaluation, including multiwavelength data to loosely constrain particle-size distributions. These evaluations of multiple models have revealed general biases in their regional aerosol amounts and seasonal patterns of transport and removal. Although satellite measurements have near-global coverage, substantial errors can be introduced into the model observation comparison unless attention is paid to spatial and temporal collocation, cloud screening, subgrid-scale variability, and measurement uncertainties that vary with retrieval conditions.

aerosol properties↗

Machine learning model inputs, outputs, and scripts associated with “Artificial intelligence-guided iterations between observations and modeling significantly improve environmental predictions”

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “Artificial intelligence-guided iterations between observations and modeling significantly improve environmental predictions” (Malhotra et al., in prep). This effort was designed following ICON (integrated, coordinated, open, and networked) principles to facilitate a model-experiment (ModEx) iteration approach, leveraging crowdsourced sampling across the contiguous United States (CONUS). New machine learning models were created every month to guide sampling locations. Data from the resulting samples were used to test and rebuild the machine learning models for the next round of sampling guidance. Associated sediment and water geochemistry and in situ sensor data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1923689, https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1729719, and https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1603775. This data package is associated with two GitHub repositories found at https://github.com/parallelworks/dynamic-learning-rivers and https://github.com/WHONDRS-Hub/ICON-ModEx_Open_Manuscript. In addition to this readme, this data package also includes two file-level metadata (FLMD) files that describes each file and two data dictionaries (DD) that describe all column/row headers and variable definitions. This data package consists of two main folders (1) dynamic-learning-rivers and (2) ICON-ModEx_Open_Manuscript which contain snapshots of the associated GitHub repositories. The input data, output data, and machine learning models used to guide sampling locations are within dynamic-learning-rivers. The folder is organized into five top-level directories: (1) “input_data” holds the training data for the ML models; (2) “ml_models” holds machine learning (ML) models trained on the data in “input_data”; (3) “examples” contains files for direct experimentation with the machine learning model, including scripts for setting up “hindcast” run; (4) “scripts” contains data preprocessing and postprocessing scripts and intermediate results specific to this data set that bookend the ML workflow; and (5) “output_data” holds the overall results of the ML model on that branch. Each trained ML model resides on its own branch in the repository; this means that inputs and outputs can be different branch-to-branch. There is also one hidden directory “.github/workflows”. This hidden directory contains information for how to run the ML workflow as an end-to-end automated GitHub Action but it is not needed for reusing the ML models archived here. Please see the top-level README.md in the GitHub repository for more details on the automation. The scripts and data used to create figures in the manuscript are within ICON-ModEx_Open_Manuscript. The folder is organized into four folders which contain the scripts, data, and pdf for each figure. Within the “fig-model-score-evolution” folder, there is a folder called “intermediate_branch_data” which contains some intermediate files pulled from dynamic-learning-rivers and reorganized to easily integrate into the workflows. NOTE: THIS FOLDER INCLUDES THE FILES AT THE POINT OF PAPER SUBMISSION. IT WILL BE UPDATED ONCE THE PAPER IS ACCEPTED WITH ANY REVISIONS AND WILL INCLUDE A DD/FLMD AT THAT POINT. We thank the United States Forest Service, Washington Department of Fish and Wildlife, Washington Department of Natural Resources, Cowiche Canyon Conservatory, Washington State Parks and Recreation Commission (Scientific Research Permit #210901), and the Confederated Tribes and Bands of the Yakama Nation for access to field locations where the samples labeled “SSS” were collected. We also thank the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview. WHONDRS consortium members were asked to provide any acknowledgments for the collection of samples labeled “CM” and the following is a list of acknowledgments that were submitted with their corresponding Site IDs: (MART) Research activities were conducted in part on the Wind River Experimental Forest within the Gifford Pinchot National Forest; (MP- 100379) Philadelphia is part of Lenapehoking, the ancestral homelands of the Lenape peoples; (MP-102398) Land surveyed is the ancestral homelands of the Nookhose'iinenno (Arapaho), Tsis tsis'tas (Cheyenne), and Nuuchu (Ute); (MP-100749 and MP- 100747) Georgia Coastal Ecosystem LTER, OCE-1832178; (SP-70 and SP-72) Eastern Shoshone, Shoshone-Bannock; (MP- 102944) Funded by Oregon Watershed Enhancement Board. On the traditional lands of the Confederated Tribes of the Siletz, Confederated Tribes of the Grand Rhonde, and the Clatsop-Nehalem Confederated Tribe; (MP- 100607) Holiday Creek is located on the traditional territory of the Monacan Indian Nation; (SP-45) Lafayette Blue Springs State Park; (MP-102420) NSF DEB-2016749; (MP-100019) New Hampshire Agriculture Experiment Station; (SP-35) Rayonier (land owner; https://www.rayonier.com/); (MP- 101276) US Department of Energy, Office of Science, Biological and Environmental Research, Subsurface Biogeochemical Research, Watershed Dynamics and Evolution SFA at ORNL; (MP- 103224) Watershed Dynamics and Evolution SFA at ORNL; (MP- 101584) Traditional lands of the Oceti Sakowin (Dakota, Lakota, Nakoda) and Anishinaabe Peoples.

54 ENVIRONMENTAL SCIENCES↗

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Variability and Anomalous Trends in the Global Sea Ice Cover

The advent of satellite data came fortuitously at a time when the global sea ice cover has been changing rapidly and new techniques are needed to accurately assess the true state and characteristics of the global sea ice cover. The extent of the sea ice in the Northern Hemisphere has been declining by about -4% per decade for the period 1979 to 2011 but for the period from 1996 to 2010, the rate of decline became even more negative at -8% per decade, indicating an acceleration in the decline. More intriguing is the drastically declining perennial sea ice area, which is the ice that survives the summer melt and observed to be retreating at the rate of -14% per decade during the 1979 to 2012 period. Although a slight recovery occurred in the last three years from an abrupt decline in 2007, the perennial ice extent was almost as low as in 2007 in 2011. The multiyear ice, which is the thick component of the perennial ice and regarded as the mainstay of the Arctic sea ice cover is declining at an even higher rate of -19% per decade. The more rapid decline of the extent of this thicker ice type means that the volume of the ice is also declining making the survival of the Arctic ice in summer highly questionable. The slight recovery in 2008, 2009 and 2010 for the perennial ice in summer was likely associated with an apparent cycle in the time series with a period of about 8 years. Results of analysis of concurrent MODIS and AMSR-E data in summer also provide some evidence of more extensive summer melt and meltponding in 2007 and 2011 than in other years. Meanwhile, the Antarctic sea ice cover, as observed by the same set of satellite data, is showing an unexpected and counter intuitive increase of about 1 % per decade over the same period. Although a strong decline in ice extent is apparent in the Bellingshausen/ Amundsen Seas region, such decline is more than compensated by increases in the extent of the sea ice cover in the Ross Sea region. The results of analysis of MODIS, AMSR-E and SSM/I data reveal that the sea ice production rate at the coastal polynyas along the Ross Ice Shelf has been increasing since 1992. This also means that the salinization rate and the formation of bottom water in the region are going up as well. Simulation studies indicate that the stronger production rate is likely associated with the ozone hole that has caused a deepening of the lows in the West Antarctic region and therefore stronger winds off the Ross Ice Shelf. Stronger winds causes larger coastal polynyas near the shelf and hence an enhanced ice production in the region during the autumn and winter period. Results of analysis of temperature data from MODIS and AMSR-E shows that the area and concentration of the sea ice cover are highly correlated with surface temperature for both the Arctic and Antarctic, especially in the seasonal regions where the correlation coefficients are about 0.9. Abnormally high sea surface temperatures (SSTs) and surface ice temperatures (SITs) were also observed in 2007 and 2011when drastic reductions in the summer ice cover occurred, This phenomenon is consistent with the expected warming of the upper layer of the Arctic Ocean on account of ice-albedo feedback. Changes in atmospheric circulation are also expected to have a strong influence on the sea ice cover but the results of direct correlation analyses of the sea ice cover with the Northern and the Southern Annular Mode indices show relatively weak correlations, This might be due in part to the complexity of the dynamics of the system that can be further altered by some phenomena like the Antarctic Circumpolar Wave and extra polar processes like the El Nino Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (POD),

Comiso, Josefino C.↗

Data-driven Community-centered Resilient Assessment and Planning Toolkit for Nexus of Energy and Water (DCRAPT-NEW)

Urban areas, including Detroit and Pittsburgh, have suffered significant dual outages of the electrical and water infrastructure in the past decade due, in part, to the increasing number of extreme weather events. With increasing temperatures and rainfall intensity, these regions need to prepare for increasing extreme events through community-based energy and water resilience analysis, planning, and enhancement. This project developed a suite of open-source, open-access, community-centered, data-driven assessment and distributed energy resource (DER) and planning tools for energy and water resilience enhancement in urban areas. Through establishing a multi-level community awareness and engagement mechanism and a comprehensive collection of power outage and flooding data, an innovative group of community energy and water resilience assessment and planning tools have been developed for a wide range of users with differing and variable sets of data available to them. The developed tools include (1) DOE EAGLE-I data-driven, deep-learning assisted resilience assessment and DER planning tools at the county level with socioeconomic factors incorporated; (2) Utility annual power outage data-driven tools for long term resilience assessment and DER planning and 15-min power outage data-driven tools for short term resilience assessment and planning; (3) Detailed engineering tools for energy and water systems resilience assessment and planning when the system topology and component fragility curves are available; (4) Alternative Resiliency Metric Calculation that extracts and separates outage and restoration processes; and (5) Co-optimization tools that evaluate the resilience of the power and sewage system and allow users to conduct joint planning with energy and wastewater systems. The developed tools provide planners, decision-makers, and stakeholders with powerful capabilities to systematically evaluate system/community resilience and optimal and actionable guidance for enhancing resilience while prioritizing DER investments. The tools have been used and validated in Detroit and Pittsburgh and can be used in other areas of the nation. In addition, this project will (1) advance the knowledge and applications of machine-learning methods in analyzing and fusing different layers of information and generating meaningful data points such as generating rare weather events; (2) significantly improve the energy and water resilience of the identified communities in Detroit and Pittsburgh and prepare for more frequent and severe weather conditions; (3) help communities assess extreme weather event impacts and address short-term and long-term resilience-related issues The developed tools have been made public via GitHub and demonstrated to community stakeholders and utility companies via the two annual workshops and numerous community engagement meetings. The project outcomes are also disseminated through publications in various journals and conference proceedings, and presentations at top conferences.

13 HYDRO ENERGY↗

Flight Mechanics Modeling and Simulation of the Earth Entry System

Introduction: The Mars Sample Return (MSR) Campaign being planned by NASA and ESA has the ambitious goal to return Mars samples back to Earth. This international collaboration had developed a concept of operations that included a ESA-designed Earth Return Orbiter (ERO) and NASA-designed Capture, Containment, and Return System (CCRS). The Earth Entry System (EES), consisting of a protective aeroshell that houses the samples as well as sample containment vessels, would conduct entry, descent, and landing (EDL) on a direct Earth trajectory. The EES would enter on a spin-stabilized ballistic trajectory with the goal to passively achieve aerodynamic stability throughout all regions of flight. The EDL sequence would end with the EES impacting the soft playa soil of the Utah Test and Training Range (UTTR). As of the submission of this abstract, the MSR campaign is undergoing a re-architecture leading to a pause in EES development. However, the novel approaches developed in flight mechanics modeling and simulation can significantly benefit the greater IPPW community in the development of Earth return vehicles. This paper will present the latest state of EES flight mechanics modeling and simulation. The paper will highlight the simulation architecture developed and key lessons learned from understanding of EDL trajectory sensitivities. Modeling and Simulation: Figure 1 provides a high-level concept of operations for the approach, entry, descent, and landing (AEDL) phase of the CCRS-portion of MSR. The objective of EES flight mechanics is to model and simulate the EES trajectory from ERO separation to ground impact at UTTR. A variety of flight mechanics simulation models were utilized to model both exo-atmopsheric and atmospheric portions of flight. 42, a 6-DOF simulation developed at Goddard Space Flight Center, is utilized for propagating the attitude of EES during exo-atmospheric flight. 42 allows for a variety of spin eject mechanism scenarios to be simulated for analysis. 10 minutes prior to entry, the 42 states are handed off to the EDL sims. The prime EDL sim utilized by EES is the Program to Optimize Simulated Trajectories II (POST2), a 6-DOF sim developed at Langley Research Center, and the independent verification and validation EDL sim utilized is DSENDS, a 6-DOF sim developed at Jet Propulsion Laboratory. Figure 2 provides a visualization of the flight mechanics simulation model flow through various points in the AEDL phase. Due to the existence of a variety of sim models, the EES flight mechanics team developed processes for data hand-off. These processes included the development of a centralized coordinate frame document, utilization of a single, centralized simulation input document for all sims to reference, and hand-off files containing both the technical data to be ingested by other flight mechanics sims as well as annotations of modeling assumptions utilized to generate the data. Figure~\ref{fig:post2simarchitecture} provides an overview of the POST2 sim architecture wherein POST2 ingests numerous subsystem models and input files. The dispersed state file generated by MONTE provides the position/velocity state of the trajectory while the 42 Handoff file provides the attitude. The aerodynamics database, delivered by the EES aeroscience team, is utilized to simulate the aerodynamic forces and moments experienced during EDL. A custom atmosphere model, developed by EES atmosphere team, is utilized to simulate the anticipated atmosphere environment around the region of Earth through which the EES trajectory flys. These inputs and subsystem models can be varied depending on the AEDL flight mechanics scenario being simulated. Monte Carlo simulations are utilized to generate statistical AEDL performance metrics in the form of scorecards and violin plots. Furthermore, outputs from the POST2 simulation are utilized for follow-on analyses including aerothermal and landing performance. \section{Flight Mechanics Lessons Learned} Though the EES flight mechanics team uncovered a variety of lessons learned through the analysis conducted to support CCRS through preliminary design review, this paper will highlight the most important lessons. A key AEDL performance goal is to ensure the landing footprint of EES remains on the UTTR south range. A common modeling strategy used in EDL analysis is One-Variable-At-a-Time (OVAT). OVAT analysis provides insight into the key drivers that affect AEDL performance metrics. Figure 3 shows the landing ellipses for single dispersion sources as compared to the baseline aggregate of all dispersions. The figure shows that atmosphere winds alone dominate the size of the footprint ellipse (note: EES does not use a parachute unlike previous Earth-return missions and is in wind-driven free fall for ~5min). The significance of the wind led the EES flight mechanics team to pursue the development of a Custom Atmosphere Model [4], in lieu of EarthGRAM [1], built on actual radiosonde wind measurements around the UTTR-region. This decision was driven by the realism in the generated footprint ellipses and lessons-learned from Stardust [5]. These findings will be invaluable for future Earth-return missions in providing an early understanding of the key drivers affecting footprint size and modeling considerations for which to account. Another lesson learned is tied to the AEDL performance goal of achieving passive stability throughout all regions of flight. It is well understood that blunt-body aeroshells are less stable as they transition from supersonic to subsonic. Eliminating a backshell does help improvestability; however, other phenomena such as roll-induced instability during terminal descent can still arise. The EES flight mechanics team developed stability metrics as tools to better understand the causes of and better predict the onset of dynamic instability. These tools were built upon analytical models developed by Jaffe [3] and Murphy [2]. The tools were shown to both be very accurate in correlation with actual unstable cases and useful in developing stability margin policies based on the vehicle design and simulation considerations (e.g. sphere-cone angle change, mass change, wind turbulence). These tools allowed for the current EES design to demonstrate the ability to achieve passive stability and can be an invaluable tool for consideration in the design of parachute-less Earth-return vehicles.

Rohan Deshmukh↗

Liquid Salt Combined-Cycle Pilot Plant Design

The work described in this report is responsive to the Office of Fossil Energy program ‘Energy Storage for Fossil Power Generation.’ This Phase I report has been prepared by Pintail Power LLC, with support from Nexant ECA, Electric Power Research Institute (EPRI) and Southern Company Services as a deliverable for the U.S. Department of Energy for NETL Award DE-FE-00320016. The Liquid Salt Combined Cycle™ (LSCC™) technology provides large-scale energy storage integrated with Fossil Electric Generating Units (FEGUs) to meet critical needs in the energy transition by providing: • the lowest cost large-scale storage for time-shifting of renewable energy, • superior fuel efficiency to reduce GHGs from dispatchable resources, • flexible capacity and ramping to balance variability of wind and solar resources, • essential grid stability services to assure reliability of a low-carbon grid. The LSCC approach: • employs equipment that has already been proven in utility service, • uses safe, non-toxic, non-degrading, perpetual-life storage medium, • leverages and repurposes existing FEGU assets, • expands the value stack of energy storage to reduce market, financing, and commodity risks. Pintail Power has developed the LSCC technology to meet the need for reliable, efficient, and cost-effective integration of Variable Renewable Energy (VRE) into a low-carbon electric grid by coupling proven thermal energy storage with proven gas turbines, steam turbines, and heat transfer equipment. This novel approach is intended to address the key issues facing the grid and operators of renewable and fossil generating units including: • Overgeneration and curtailment of renewables, • Need for fast ramping dispatchable resources, • Improved efficiency and flexibility of fossil units, • Additional peaking capacity to support electrification of transportation and heating, • Provision of reliability services to support high penetration of VRE, especially synchronous inertia and fast frequency response. A Technology Readiness assessment by EPRI confirmed that LSCC technology consists of commercially proven hardware used in industrial and utility applications. Although the novel LSCC approach has not yet been demonstrated as a complete system, interfaces between major components have been conservatively specified. A Phase III pilot is planned to demonstrate equipment integration and operation. The patented innovation is removal of the evaporator section from the exhaust heat recovery system, with the evaporation performed by stored energy in a separate steam generator. This arrangement couples renewable and fossil power generation via long-duration energy storage to deliver cost, performance, and operational synergies, including superior charging and discharging flexibility, reduced fuel consumption and lower CO 2 emissions compared to conventional Combined Cycle Power Plants, and low-cost, large-scale energy storage. The LSCC technology is composed of proven equipment integrated with gas turbine exhaust heat in a novel system. During charging, electric heaters raise the salt temperature as it flows from the Cold Salt Tank to the Hot Salt Tank. During discharging, hot salt produces steam from feedwater that is heated with gas turbine exhaust, which also superheats steam to drive a steam turbine. LSCC technology can be added to any combustion-turbine to integrate renewable energy, provide needed grid services, and increase the value of fossil electric generating units based on the technology’s following attributes: • Long-duration storage enables time-shifting of VRE to avoid curtailment and impairment of renewable assets. • Long storage duration combined with fast-charging capability increases arbitrage opportunities by storing more energy when the price is low and discharging more hours when the price is high. • Long storage duration allows resource adequacy to be supplied across multiple days to increase reliability and reduce risk. • The stored energy reduces fuel heat rate and GHG emissions, and increases merit, so the LSCC dispatches earlier and longer to increase the plant’s capacity factor and asset value. • The stored energy enables pre-heating and startup of the steam cycle, without operating the gas turbine, to enable fast startup and ramping when dispatched for discharge. • The steam turbine can operate without the gas turbine so it can provide valuable synchronous inertia during charging without consuming fuel. • Fast frequency response and regulation services can be provided during charging using solid-state heater and pump controls to vary the charge power input in response to grid signals. • The LSCC system can be configured for resilience including black start, islanded/micro-grid operation, and even self-recharging of storage using either gas turbine power or gas turbine exhaust heat. The commercialization plan is to add LSCC technology to existing simple cycle gas turbine power plants with the 50MW GE LM6000 aero-derivative gas turbine as the reference design basis. A Techno-economic assessment of the reference design evaluated the benefits (Levelized Avoided Cost of Energy) and costs (Levelized Cost of Energy). The plant definition included all major systems and budgetary vendor quotes. Pintail Power and NexantECA developed the overall cost estimate for the LSCC plant up to the total plant cost level, following the DOE-NETL cost estimate guidelines at AACE Class 3 (-20%/+30%). This includes the equipment cost, bulk material, direct and indirect labor costs to arrive at the bare erected cost. Engineering costs are factored from the BEC and added to it to arrive at the EPC cost. Process and project contingencies were then factored from the EPC cost and rolled-up to yield the total plant cost of $\$$184 million for 1746 MWh of discharge electricity. • At $\$$105/kWh, the reference plant costs less than any of the Energy Storage Systems evaluated by PNNL in 2020 for the Energy Storage Grand Challenge. Operations and Maintenance cost estimates were scaled from combined cycle practice, assuming that the LSCC unit was co-located with and sharing some labor expense with other units, to arrive at $\$$2.2 million per year. Plant economics were evaluated using prices from the ERCOT Day-Ahead Market for calendar year 2019 (excluding the market disruptions from the COVID pandemic and the February 2020 deep freeze event). Assuming economic dispatch in the ERCOT Day-Ahead market, the reference plant capacity factor would have discharged for 2777 hours at 91.9 MW, a 31.66% capacity factor, with a marginal cost of $\$$25.59/MWh, and a LACE of $\$$82.41/MWh. Fixed charges were calculated according to EIA guidelines to arrive at an LCOE of $\$$83.48. The benefit-to-cost ratio of 0.99 suggests that the reference plant would have been cost-effective and competitive in the market. EPRI interviewed selected utilities to gauge the need for, applicability of and interest in the LSCC system. Several utilities are currently managing increased load growth along with the inclusion of increasing levels of renewable generation, putting pressure on conventional generation by requiring increased turndown requirements and ultimately lower capacity factors. All of the utilities interviewed have CO 2 reduction targets in the 2030-2050 timeframe that will severely limit the participation of fossil generation and require better utilization of carbon free generation. While there is limited opportunity for storage in the current markets, the utilities interviewed stated that there will be a substantial need for long duration energy storage in the future given the expected trends. Utilizing an energy storage system will generally be preferred over new gas capacity in some cases, with the capabilities of the LSCC system being a potential option for retrofit to existing simple cycle gas turbine units, allowing them to deliver greater participation in the market with lower carbon intensity. A technology gap assessment and technology maturation plan identified a pilot-scale demonstration as the final step before commercialization. Key gaps to be addressed during the Phase II FEED (Front-End Engineering Design) are component selection and design, commissioning procedures, and operational procedures and the control system for LSCC charging and discharging. The project team has been expanded to include Wood Group PLC as EPC. The proposed Phase II work leads to a pilot-scale engineering demonstration (TRL 6) to be conducted at Southern Company’s Plant Rowan, where the prototype system will perform “all the functions that will be required of the operational system.” The proposed pilot will facilitate commercialization (TRL-9) by scale-up to utility-scale systems integrated with peaking GTs or directly to facility scale systems using industrial GTs. The conceptual design for the pilot plant focuses on the novel integration aspects of LSCC technology. A slipstream of gas turbine exhaust will feed a waste heat recovery unit coupled to a molten salt steam generator heated by stored energy. The pilot is intended to demonstrate all key operating modes of the LSCC technology during charging, discharging and standby. The pilot equipment will be approximately one-seventh scale of the LM6000 commercial target and is expected to have commercial off-ramp potential for facility-scale applications.

01 COAL, LIGNITE, AND PEAT↗

GC/MS Method Development for Separating Lunar Volatile Ice Simulant Headspace Gases

Various investigators propose the lunar surface contains widely distributed volatiles, especially water- like species, i.e. OH and H2O. Surface volatiles are theorized to exist as a hydrated regolith layer, concentrated in extremely cold polar permanently shadowed regions (PSR), and/or solar wind implantation reservoirs in lunar glasses. The proposed sources of lunar surface volatiles range from cometary impacts, solar wind, or a supply present during moon formation. Future Artemis missions aim to collect and return the samples containing volatiles collected near lunar polar craters or PSRs. We, as advanced curation scientists, are responsible for developing techniques and methodologies for preserving returned sample integrity as much as possible. Pristine volatile-bearing samples are invaluable to the scientific community seeking to unravel the history of the solar system. Realistically, a sample will experience alteration during collection, transportation back to earth, and storage. The Planetary Exploration and Astromaterials Research Lab (PEARL) seeks to understand temperature and pressure effects on high-fidelity volatile-containing regolith simulants, the foundation for the future of cold curation. This abstract outlines the separation, identification, and quantification of headspace gases over volatile ice feed stock material using gas chromatography/mass spectrometry (GC/MS). Preliminary objectives concentrated on sample handling, reproducibility, and understanding the elution characteristics for each analyte. Initial GC/MS method development experiments utilized diluted static headspace sample preparation. Diluted samples were used because sampling headspace gases directly from a vial containing liquid analyte resulted in overloading of the column and detector. Overloading is evident based on chromatogram peak shapes and instrument contamination, or carry over, between experiments. A mixture of three alcohols were used for a majority of the sample handling and reproducibility studies. Reproducibility was tested via multiple users, calibration curves, and check standards. Stock solutions of condensed lunar volatile analytes included methanol, ammonia in methanol, hydrogen sulfide in water, and an equal volume mixture of methanol, ethanol, and isopropanol. Current samples use room air as the headspace sample matrix, however future experiments will incorporate an inert purge gas, such as argon or nitrogen. Three mL of each analyte solution were capped in separate 20 mL crimp top GC vials. Dilutions were carried out by removing an aliquot of headspace gases with a calibrated 1 mL gastight syringe and immediately transferring to a 20 mL capped crimp top vial. The GC/MS is a Thermo Fisher Trace 1310/ISQ 7000 with a TriPlus RSH autosampler and split/splitless injector module. The experiments outlined in this abstract use the following hardware: a 2.5 mL gastight headspace syringe tool, 1 mm ID x 78.5 mm length ultra-inert straight injection liner, and a TG-BondQ 30 m × 0.32 mm × 10 μm column. Various parameters, such as hardware selection and the temperature, pressure, and split ratio set points, continue to evolve as the overall experiment is refined. Diluted headspace chromatograms were collected for the individual stock solutions. Retention times, peak shapes, and mass spectra were evaluated and added to the data processing method for each molecule of interest. Figure 1 shows the total ion chromatograms for the three major lunar volatile simulant stock solutions: methanol, 7 N ammonia in methanol, and 0.4% hydrogen sulfide in water. Tailing peak shapes for ammonia (2.98 min rt) and water (4.06 min rt) indicate the molecules are not properly eluting from the selected column with the current separation method. Additionally, hydrogen sulfide and ammonia have overlapping peak windows, which could impact quantification. Ongoing experiments aim to address the peak shape and overlapping via the separation method and hardware selection. Sample preparation reproducibility experiments used stock solution containing equal volumes of a non- interactive mixture of methanol, ethanol, and isopropanol. Mass spectrum ion traces were used to identify and quantify all three alcohols. Peaks were automatically detected, identified, and integrated through the mass spectra detection and processing parameters. Calibration response curves and check standards were used to evaluate the validity of the sample preparation procedure. Figure 2 shows the methanol chromatogram peak area versus total headspace dilution volume transferred from the alcohol mixture vial. The calibration response curves and check standards validate sample preparation procedure. Continuing data analysis efforts are working towards correlating the peak area and instrument response factor to the headspace analyte concentration and condensed phase composition. Static headspace gas chromatography theory relies on Dalton’s law, Raoult’s law, Henry’s Law, and the Kolb and Ettre equation to associate peak area to the analyte composition in a non-ideal solution. Equation 1 is a simplified expression derived from the aforementioned theories. Future experiments involve liquid injections of the individual stock solutions, liquid and headspace analysis of various stock solution combinations, and the addition of regolith simulants to the mixtures. Temperature is another variable expected to affect reaction rates and will be explored.

Cecilia L. Amick↗