Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Prior Distribution”

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 1,513 records · Page 84

High Performance EVA Glove Collaboration: Glove Injury Data Mining Effort

Human hands play a significant role during Extravehicular Activity (EVA) missions and Neutral Buoyancy Lab (NBL) training events, as they are needed for translating and performing tasks in the weightless environment. Because of this high frequency usage, hand and arm related injuries are known to occur during EVA and EVA training in the NBL. The primary objectives of this investigation were to: 1) document all known EVA glove related injuries and circumstances of these incidents, 2) determine likely risk factors, and 3) recommend interventions where possible that could be implemented in the current and future glove designs. METHODS: The investigation focused on the discomforts and injuries of U.S. crewmembers who had worn the pressurized Extravehicular Mobility Unit (EMU) spacesuit and experienced 4000 Series or Phase VI glove related incidents during 1981 to 2010 for either EVA ground training or in-orbit flight. We conducted an observational retrospective case-control investigation using 1) a literature review of known injuries, 2) data mining of crew injury, glove sizing, and hand anthropometry databases, 3) descriptive statistical analyses, and finally 4) statistical risk correlation and predictor analyses to better understand injury prevalence and potential causation. Specific predictor statistical analyses included use of principal component analyses (PCA), multiple logistic regression, and survival analyses (Cox proportional hazards regression). Results of these analyses were computed risk variables in the forms of odds ratios (likelihood of an injury occurring given the magnitude of a risk variable) and hazard ratios (likelihood of time to injury occurrence). Due to the exploratory nature of this investigation, we selected predictor variables significant at p≤0.15. RESULTS: Through 2010, there have been a total of 330 NASA crewmembers, from which 96 crewmembers performed 322 EVAs during 1981-2010, resulting in 50 crewmembers being injured inflight and 44 injured during 11,704 ground EVA training events. Of the 196 glove related injury incidents, 106 related to EVA and 90 to EVA training. Over these 196 incidents, 277 total injuries (126 flight; 151 training) were reported and were then grouped into 23 types of injuries. Of EVA flight injuries, 65% were commonly reported to the hand (in general), metacarpophalangeal (MCP) joint, and finger (not including thumb) with fatigue, abrasion, and paresthesia being the most common injury types (44% of total flight injuries). Training injuries totaled to more than 70% being distributed to the fingernail, MCP joint, and finger crotch with 88% of the specific injuries listed as pain, erythema, and onycholysis. Of these training injuries, when reporting pain or erythema, the most common location was the index finger, but when reporting onycholysis, it was the middle finger. Predictor variables specific to increased risk of onycholysis included: female sex (OR=2.622), older age (OR=1.065), increased duration in hours of the flight or training event (OR=1.570), middle finger length differences in inches between the finger and the EVA glove (OR=7.709), and use of the Phase VI glove (OR=8.535). Differentiation between training and flight and injury reporting during 2002-2004 were significant control variables. For likelihood of time to first onycholysis injury, there was a 24% reduction in rate of reporting for each year increase in age. Also, more experienced crewmembers, based on number of EVA flight or training events completed, were less likely to report an onycholysis injury (3% less for every event). Longer duration events also found reporting rates to occur 2.37 times faster for every hour of length. Crewmembers with larger hand size reported onycholysis 23% faster than those with smaller hand size. Finally, for every 1/10th of an inch increase in difference between the middle finger length and the glove, the rate of reporting increased by 60%. DISCUSSION: One key finding was that the Series 4000 glove had a lower injury risk than the Phase VI, which provides a platform for further evaluation. General interventions that reduce hand overexertion and repetitive use exposure through tool development, procedural changes and shorter exposures may be one mitigation path, but due to the way the training event times were reported, we cannot provide a guideline for a specific event duration change. When the finger length was different from the glove length, the risk of injury increased indicating that the use of larger finger take-ups could be contributing to injury and therefore may not be recommended. Prior to this investigation, there was one previous investigation indicating hand anthropometry may be related to onycholysis. We found different hand anthropometry variables indicated by this investigation as compared to the prior, specifically differences in middle finger length compared to glove finger length, which point more towards a sizing issue than a specific anthropometry issue. Additionally, although this investigation has identified sizing as an issue, the force and environmental-related variables of the EVA glove that could also cause injury were not accounted for.

Reid, C. R.↗

Experiments on Nucleation in Different Flow Regimes

The vast majority of metallic engineering materials are solidified from the liquid phase. Understanding the solidification process is essential to control microstructure, which in turn, determines the properties of materials. The genesis of solidification is nucleation, where the first stable solid forms from the liquid phase. Nucleation kinetics determine the degree of undercooling and phase selection. As such, it is important to understand nucleation phenomena in order to control solidification or glass formation in metals and alloys. Early experiments in nucleation kinetics were accomplished by droplet dispersion methods. Dilitometry was used by Turnbull and others, and more recently differential thermal analysis and differential scanning calorimetry have been used for kinetic studies. These techniques have enjoyed success; however, there are difficulties with these experiments. Since materials are dispersed in a medium, the character of the emulsion/metal interface affects the nucleation behavior. Statistics are derived from the large number of particles observed in a single experiment, but dispersions have a finite size distribution which adds to the uncertainty of the kinetic determinations. Even though temperature can be controlled quite well before the onset of nucleation, the release of the latent heat of fusion during nucleation of particles complicates the assumption of isothermality during these experiments. Containerless processing has enabled another approach to the study of nucleation kinetics. With levitation techniques it is possible to undercool one sample to nucleation repeatedly in a controlled manner, such that the statistics of the nucleation process can be derived from multiple experiments on a single sample. The authors have fully developed the analysis of nucleation experiments on single samples following the suggestions of Skripov. The advantage of these experiments is that the samples are directly observable. The nucleation temperature can be measured by noncontact optical pyrometry, the mass of the sample is known, and post-processing analysis can be conducted on the sample. The disadvantages are that temperature measurement must have exceptionally high precision, and it is not possible to isolate specific heterogeneous sites as in droplet dispersions. Levitation processing of refractory materials in ultra high vacuum provides an avenue to conduct these kinetic studies on single samples. Two experimental methods have been identified where ultra high vacuum experiments are possible; electrostatic levitation in ground-based experiments and electromagnetic processing in low earth orbit on TEMPUS. Such experiments, reported here, were conducted on zirconium. Liquid zirconium is an excellent solvent and has a high solubility for contaminants contained in the bulk material as well as those contaminants found in the vacuum environment. Oxides, nitrides, and carbides do not exist in the melt, and do not form on the surface of molten zirconium, for the materials and vacuum levels used in this study. Ground-based experiments with electrostatic levitation have shown that the statistical nucleation kinetic experiments are viable and yield results which are consistent with classical nucleation theory. The advantage of low earth orbit experiments is the ability to vary the flow conditions in the liquid prior to nucleation. The put-pose of nucleation experiments in TEMPUS was to examine.

Bayuzick, R. J.↗

Validation of Long Bone Mechanical Properties from Densitometry

The objective of this study was to assess whether cross-sectional areal properties, calculated from densitometry, correlate to the true flexural properties. Right and left male embalmed tibiae were used in the study. Prior to scanning, the proximal end of each tibia was potted in a fixture with registration pins, flushed thoroughly with water under pressure to remove trapped air, and then placed in a constant thickness water bath attached to a precision indexer. Two sets of three scans of the entire tibia were taken with an Hologic QDR 1000/W densitometer at rotations of 0, 45, and 90 degrees about the tibia long axis. An aluminum step phantom and a bone step phantom, machined from bovine cortical bone, were also in the bath and scanned separately. Pixel attenuation data from the two sets of scans were averaged to reduce noise. Pixel data from the high energy beam were then converted to equivalent thicknesses using calibration equations. Cross-sectional areal properties (centroid, principal area moments and principal angle) along the length were computed from the three registered scans using methods developed in our laboratory. Flexural rigidities. Four strain gages were bonded around the circumference of each of 5 cross-sections encompassing the entire diaphysis. A known transverse load was then applied to the distal end and the bone was rotated 360 degrees in eight increments of 45 degrees each. Strains from the eight orientations were analyzed along with the known applied bending moments at each section to compute section centroids, curvatures, principal flexural rigidities and principal angle. Reference axes between the two methods were maintained within +/- 0.5 degrees using an electronic inclinometer. Principal angles (flexural - areal) differed by -2.0 +/- 4.0 degrees, and 1.0 +/- 2.5 degrees for the right and left tibia, respectively. Section principal flexural rigidities were highly correlated to principal areal moments (right: r(sup 2)= 0.997; left: r(sup 2)= 0.978) indicating a nearly constant effective flexural modulus. Right and left tibia exhibited a very high degree of symmetry when comparing either flexural or areal properties. To our knowledge this is the first study to validate the use of densitometry (DXA) to predict three dimensional structural properties of long bones. Our initial results support the conclusion that bone mineral and its distribution are the primary determinants of flexural modulus and rigidity.

Whalen, R.↗

Irrigation Characterization Improved by the Direct Use of SMAP Soil Moisture Anomalies Within a Data Assimilation System

Prior soil moisture data assimilation (DA) efforts to incorporate human management features such as agricultural irrigation has only shown limited success. This is partly due to the fact that observational rescaling approaches for bias correction used in soil moisture DA systems are less effective when unmodeled processes such as irrigation are the dominant source of systematic biases. In this article, we demonstrate an alternative approach, i.e. anomaly correction for overcoming this limitation. Unlike the rescaling approaches, the proposed method does not scale remote sensing soil moisture retrievals to the model climatology, but it extracts the temporal variability information from the retrievals. The study demonstrates this approach through the assimilation of soil moisture retrievals from the Soil Moisture Active Passive mission into the Noah land surface model. The results demonstrate that DA using the anomaly correction method can better capture the effect of irrigation on soil moisture in agricultural areas while providing comparable performance to the DA integrations using rescaling approaches in non-irrigated areas. These findings emphasize the need to reduce inconsistencies between remote sensing and the models so that assimilation methods can employ information from remote sensing more directly to develop representations of unmodeled processes such as irrigation.

soil moisture↗

Soft-Switching Solid-State Transformer With Reduced Conduction Loss

Solid-state transformers (SSTs) are a promising solution photovoltaic (PV), wind, traction, data center, battery energy storage system (BESS), and fast charging electric vehicle (EV) applications. The traditional SSTs are typically three-stage, i.e., hard-switching cascaded multilevel rectifiers and inverters with dual active bridge (DAB) converters, which leads to bulky passives, low efficiency, and high electromagnetic interference (EMI). This article proposes a new soft-switching solid-state transformer (S4T). The S4T has full-range zero-voltage switching (ZVS), electrolytic capacitor-less dc link, and controlled dv/dt, which reduces EMI. The S4T comprises two reverse-blocking current-source inverter (CSI) bridges, auxiliary branches for ZVS, and transformer magnetizing inductor as a reduced dc link with 60% ripple. Compared with the prior S4T, an effective change on the leakage inductance diode is made to reduce the number of the devices on the main power path by 20% for significant conduction loss saving and retain the same functionality of damping the resonance between the leakage and resonant capacitors and recycling trapped leakage energy. The conduction loss saving is crucial, being the dominating loss mechanism in SSTs. Importantly, the proposed single-stage SST not only holds the potential for high power density and high efficiency but also has full functionality, e.g., multiport dc loads integration, voltage regulation, and reactive power compensation, unlike the traditional single-stage matrix SST. The S4T can achieve single-stage isolated bidirectional dc–dc, ac–dc, dc–ac, or ac–ac conversion. It can also be configured input-series output-parallel (ISOP) in a modular way for medium-voltage (MV) grids. Hence, the S4T is a promising candidate for the SST. The full functionality, e.g., voltage buck–boost, multiport, etc., and the universality of the S4T for the dc–dc, dc–ac, and ac–ac conversion are verified through the simulations and experiments of two-port and three-port MV prototypes based on 3.3 kV SiC MOSFETs in dc–dc, dc–ac, and ac–ac modes at 2 kV.

14 SOLAR ENERGY↗

Microgrid Conceptual Design Guidebook | 2022

This guide is meant to assist communities – from residents to energy experts to decision makers – in developing a conceptual microgrid design that meets site-specific energy resilience goals. Using the framework described in this guidebook, stakeholders can come together and start to quantify site-specific vulnerabilities, identify the most significant risks to delivery of electricity, and establish electric outage tolerances across the community. In addition to establishing minimum service needs, this framework encourages communities to consider broader sustainability goals and policy constraints and begin to estimate up-front costs associated with the installation of alternative microgrid solutions. The framework guides a community through data collection and a high-level assessment of its needs, constraints, and priorities, prior to engaging engineers, vendors, and contractors. The first sections of this guidebook provide a high-level primer on electric systems. The latter sections include guidance for step-by-step data gathering and analysis of site conditions. The ultimate product resulting from the stepwise approach is a conceptual microgrid design. A conceptual design is defined as an initial design (10%-20% complete) that considers the specific threats, needs, limitations, and investment options for a given location. Going through this exercise and developing the conceptual microgrid design as a community ensures the same community members who will ultimately live with the solution are the developers of its foundational design. Often, these are also the very same people who understand system tolerances and needs the best and are therefore the ideal candidates for establishing these criteria. Especially when it comes to evaluating critical infrastructure, it is the community that best understands the most critical services. The framework is intended to facilitate a systematic approach to planning for resilience and provide a deeper understanding of how to use a framework to make decisions around microgrid solutions. Like many processes where tradeoffs need to be considered, this is often an iterative process. If this guide serves to help educate and empower communities who are beginning the process of deploying a microgrid, it has met the goal of its authors.

24 POWER TRANSMISSION AND DISTRIBUTION↗

High Water Contents in the Siberian Cratonic Mantle: An FTIR Study of Udachnaya Peridotite Xenoliths

Water is believed to be a key factor controlling the long-term stability of cratonic lithosphere, but mechanisms responsible for the water content distribution in the mantle remain poorly constrained. Water contents were obtained by FTIR in olivine, pyroxene and garnet for 20 well-characterized peridotite xenoliths from the Udachnaya kimberlite (central Siberian craton) and equilibrated at 2-7 GPa. Water contents in minerals do not appear to be related to interaction with the host kimberlite. Diffusion modeling indicates that the core of olivines preserved their original water contents. The Udachnaya peridotites show a broad range of water contents in olivine (6.5 +/- 1.1 to 323 +- 65 ppm H2O (2 sigma)), and garnet (0 - 23 +/- 6 ppm H2O). The water contents of olivine and garnet are positively correlated with modal clinopyroxene, garnet and FeO in olivine. Water-rich garnets are also rich in middle rare earth elements. This is interpreted as the result of interaction between residual peridotites and water rich-melts, consistent with modal and cryptic metasomatism evidenced in the Siberian cratonic mantle. The most water-rich Udachnaya minerals contain 2 to 3 times more water than those from the Kaapvaal craton, the only craton with an intact mantle root for which water data is available. The highest water contents in olivine and orthopyroxene in this study (>= 300 ppm) are found at the bottom of the lithosphere (> 6.5 GPa). This is in contrast with the Kaapvaal craton where the olivines of peridotites equilibrated at > 6.4 GPa have < 1 ppm H2O. The latter "dry" olivine may make the base of the Kaapvaal cratonic root strong and thus protects it from erosion by the convective mantle The calculated viscosity for water-rich Udachnaya peridotites at > 6 GPa is lower or similar (8.4× 10(exp 16) to 8.0× 10(exp 18) Pa./s) to that of the asthenosphere (<= 3.7x10(exp 18) Pa./s ). Such lithologies would not be able to resist delamination by the convecting asthenosphere. However, seismology studies as well as the high equilibration pressures of our samples indicate that the Udachnaya cratonic lithosphere is 220-250 km thick. Consequently, the water-rich peridotites are likely not representative of the overall Siberian cratonic lithosphere. Their composition is linked to spatially limited melt metasomatism in mantle regions above asthenospheric upwellings responsible for the kimberlite magmatism prior to their ascent and eruption.

Doucet, Luc S.↗

Electra: A Modular-Based Expansion of NASA's Supercomputing Capability

NASA has increasingly relied on high-performance computing (HPC) re- sources for computational modeling, simulation, and data analysis to meet the science and engineering goals of its missions in space exploration, aeronautics, and Earth and space science. The NASA Advanced Supercomputing (NAS) Division at Ames Research Center in Silicon Valley, Calif., hosts NASA’s premier supercomputing resources, integral to achieving and enhancing the success of the agency’s missions. NAS provides a balanced environment, funded under the High-End Computing Capability (HECC) project, comprised of world-class supercomputers, including its flagship distributed-memory cluster, Pleiades; high-speed networking; and massive data storage facilities, along with multi-disciplinary support teams for user support, code porting and optimization, and large-scale data analysis and scientific visualization. However, as scientists have increased the fidelity of their simulations and engineers are conducting larger parameter-space studies, the requirements for supercomputing resources have been growing by leaps and bounds. With the facility housing the HECC systems reaching its power and cooling capacity, NAS undertook a prototype project to investigate an alternative approach for housing supercomputers. Modular supercomputing, or container-based computing, is an innovative concept for expanding NASA’s HPC capabilities. With modular supercomputing, additional containers—similar to portable storage pods—can be connected together as needed to accommodate the agency’s ever-increasing demand for computing resources. In addition, taking advantage of the local weather permits the use of cooling technologies that would additionally save energy and reduce annual water usage. The first stage of NASA’s Modular Supercomputing Facility (MSF) prototype, which resulted in a 1,000 square-foot module on a concrete pad with room for 16 compute racks, was completed in Fall 2016 and an SGI (now HPE) computer system, named Electra, was deployed there in early 2017. Cooling is performed via an evaporative system built into the module, and preliminary experience shows a Power Usage Effectiveness (PUE) measurement of 1.03. Electra achieved over a petaflop on the LINPACK benchmark, sufficient to rank number 96 on the November 2016 TOP500 list [14]. The system consists of 1,152 InfiniBand-connected Intel Xeon Broadwell-based nodes. Its users access their files on a facility-wide file system shared by all HECC compute assets via Mellanox MetroX InfiniBand extenders, which connect the Electra fabric to Lustre routers in the primary facility over fiber-optic links about 900 feet long. The MSF prototype has exceeded expectations and is serving as a blueprint for future expansions. In the remainder of this chapter, we detail how modular data center technology can be used to expand an existing compute resource. We begin by describing NASA’s requirements for supercomputing and how resources were provided prior to the integration of the Electra module-based system.

Biswas, Rupak↗

HERO CarbonSAFE Phase 2 Project in the Columbia River Basalt Group

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Basalt CO2 storage has several advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. Associated decreases in pressure and CO2 migration risks, 3. Reduced long-term monitoring requirements with respect to plume tracking, 4. Widespread geographic distribution and, 5. Large storage potential due to thickness, porosity, and CO2 interactions with basalt. And for locations such as the Pacific Northwest, Hawaii, Iceland, India and Japan, basalts may offer the only economically feasible option for local CO2 storage. However, there are limited field-scale assessments of CO2 storage in basalt, and current carbon capture utilization and storage (CCUS) permitting and regulatory frameworks were developed for conventional saline reservoirs. HERO CarbonSAFE is designed to address research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts. The HERO CarbonSAFE project will assess feasibility of developing a commercial-scale (50+ million metric tons of CO2) geological storage complex within the Columbia River Basalt Group (CRBG), a layered continental flood basalt complex that underlies Calpine’s natural gas-fired Hermiston Power Project (HPP) in Hermiston, OR (Figure 1). Under this 2-year CarbonSAFE Phase II project, the HERO team will conduct a data acquisition campaign that includes drilling a stratigraphic well to a total depth of ~1,500 m into the thick layered basalts proximal to HPP. A comprehensive well logging and hydrologic testing program will be augmented with new core collected from flow zones and sealing units, and comprehensive laboratory testing to help refine the kinetic rates of mineralization. The newly acquired information will be integrated with existing data from regional wells to correlate basalt injection zone properties to develop storage hub/commercial-scale models. Using these models, the project team will evaluate injection scenarios to define the technical and economic potential for storing a minimum of 50 million metric tons of CO2 over a 30-year period, along with a robust sensitivity analysis on key parameters governing reservoir viability for sustainable injection over a commercial project lifetime. Specific technical objectives of HERO are: (1) assessing the reservoir response of a series of stacked layered reservoir flowtop sequences occurring in this area of the CRBG to commercial-scale injection volumes; (2) extending prior efforts by the project team to characterize the deep layered basalts encountered in regional studies, to leverage prior investments by U.S. Department of Energy’s (DOE) Carbon Storage program; (3) leveraging DOE’s mineralization characterization efforts to advance model parametrization for commercial scale injection of CO2 in basalts; (4) conducting risk assessments associated with scaling up to commercial storage hub injection goals, while validating DOE’s National Risk Assessment Partnership (NRAP) tools, to identify potential constraints that would prevent the CRBG from serving as a commercial-scale storage complex; (5) developing mitigation plans to address identified risks; (6) developing a commercial-scale injection and monitoring, verification and accounting (MVA) strategy; (7) utilizing computational models to define and minimize, if possible, the Area of Review (AoR) under Class VI regulations; and (8) developing a robust CO2 management strategy for CRBG that also considers a regional source/sink approach that is responsive to stakeholder needs and industrial demand. Specific institutional objectives are: (1) identifying and developing plans to mitigate the nontechnical challenges associated with the build-out of a commercial-scale storage complex within the CRBG with integrated CO2 sources; (2) implementing the community outreach plan; (3) conducting regulatory research, including a survey of issues related to pore space ownership, MVA and long-term assurance of mineralization-based storage, to support an eventual application for a UIC Class VI permit; (4) advancing the project’s plan for CO2 liability management; and (5) continuing to refine and update the project’s economic model. The final objective is the preparation of a comprehensive Site Characterization Plan that draws upon the technical and institutional feasibility assessments to prepare the project for future commercialization efforts.

58 GEOSCIENCES↗

Impact of Lightning and Convection on Reactive Nitrogen in the Tropical Free Troposphere

Latitudinal distributions of NO, NO(y), O3, CO, CH3I, and H2O mixing ratios at 8.9-12 km were obtained between 30deg N and 1deg S by DC-8 aircraft measurements made in February 1994 during Pacific Exploratory Mission-West B (PEM-West B). Very low NO(y), mixing ratios with a median value of 51 parts per trillion by volume (pptv) were observed at 9.5-12 km at 1deg N-14deg N during two flights made within 3 days. A very low median O3 mixing ratio of 19 parts per billion by volume (ppbv) and high mixing ratios of H2O and CH3I were simultaneously observed, suggesting that the low NO(y), values were probably due to the convective transport of air from the tropical marine boundary layer to this altitude. The median NO(y)/O3 ratio being a factor of 2 smaller than in the air masses in the tropical marine boundary layer might suggest the possibility that the heterogeneous removal of HNO3 during convective transport further reduced NO(y) levels. In addition to the measurements between 9.5 and 12 km, low values of NO(y) and O3 were observed between 4 and 12 km at 1deg N. Divergent wind fields at 200 and 1000 hPa and infrared (IR) cloud images show that there was large scale convection (greater than 1000 km x 1000 km) in the northeast of New Guinea Island centered around Odeg S and 150deg E as part of systematic convective activity of the Intertropical Convergence Zone (ITCZ) and the South Pacific Convergence Zone (SPCZ). This type of large scale convection could have transported air with low levels of NO(y) and O3 to the middle and upper troposphere over a wide area in the tropics. On the other hand, NO mixing ratios of 50-200 pptv and high NQ,/NOY ratios of 0.4-0.6 were observed at 9.5 km between 4deg S and 10deg S. High H2O Mixing ratios of 600-1200 parts per million by volume (ppmv) and low CO mixing ratios of 65 ppbv observed in the air mass indicated that the high NO values were probably due to NO production by lightning. Satellite observations showed relatively frequent lightning flashes over the New Guinea Island for 3 days prior to the aircraft measurements. These results are considered to be consistent with the idea that, in general, marine convection is not accompanied by lightning activity, whereas convection over land is. Because of the large areal extent of the influences from these processes, the convective transport of low NO(y) air and NO production by lightning should play critical roles in controlling the abundance of reactive nitrogen in the equatorial region.

Kawakami, S.↗

Physics-Informed Gaussian Process Regression for States Estimation and Forecasting in Power Grids

Real-time state estimation and forecasting are critical for the efficient operation of power grids. In this paper, a physics-informed Gaussian process regression (PhI-GPR) method is presented and used for forecasting and estimating the phase angle, angular speed, and wind mechanical power of a three-generator power grid system using sparse measurements. In standard data-driven Gaussian process regression (GPR), parameterized models for the prior statistics are fit by maximizing the marginal likelihood of observed data. In the PhI-GPR method, we propose to compute the prior statistics offline by solving stochastic differential equations (SDEs) governing the power grid dynamics. The short-term forecast of a power grid system dominated by wind generation is complicated by the stochastic nature of the wind and the resulting uncertainty in wind mechanical power. Here, we assume that the power grid dynamics are governed by swing equations, with the wind mechanical power fluctuating randomly in time. We solve these equations for the mean and covariances of the power grid states using the Monte Carlo simulation method. We demonstrate that the proposed PhI-GPR method can accurately forecast and estimate observed and unobserved states. For the considered problem, PhI-GPR has computational advantages over the ensemble Kalman filter (EnKF) method: In PhI-GPR, ensembles are computed offline and independently of the data acquisition process, whereas for EnFK, ensembles are computed online with data acquisition, rendering real-time forecast more challenging. We also demonstrate that the PhI-GPR forecast is more accurate than the EnKF forecast when the random mechanical wind power is non-Markovian. In contrast, the two methods produce similar forecasts for the Markovian mechanical wind power. For observed states, we show that PhI-GPR provides a forecast comparable to the standard data-driven GPR; both forecasts are significantly more accurate than the autoregressive integrated moving average (ARIMA) forecast. We also show that the ARIMA forecast is more sensitive to observation frequency and measurement errors than the PhI-GPR forecast.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Cast Glance Near Infrared Imaging Observations of the Space Shuttle During Hypersonic Re-Entry

High resolution calibrated infrared imagery of the Space Shuttle was obtained during hypervelocity atmospheric entries of the STS-119, STS-125 and STS128 missions and has provided information on the distribution of surface temperature and the state of the airflow over the windward surface of the Orbiter during descent. This data collect was initiated by NASA s Hypersonic Thermodynamic Infrared Measurements (HYTHIRM) team and incorporated the use of air- and land-based optical assets to image the Shuttle during atmospheric re-entry. The HYTHIRM objective is to develop and implement a set of mission planning tools designed to establish confidence in the ability of an existing optical asset to reliably acquire, track and return global quantitative surface temperatures of the Shuttle during entry. On Space Shuttle Discovery s STS-119 mission, NASA flew a specially modified thermal protection system tile and instrumentation package to monitor heating effects from boundary layer transition during re-entry. On STS-119, the windward airflow on the port wing was deliberately disrupted by a four-inch wide and quarter-inch tall protuberance built into the modified tile. In coordination with this flight experiment, a US Navy NP-3D Orion aircraft was flown 28 nautical miles below Discovery and remotely monitored surface temperature of the Orbiter at Mach 8.4 using a long-range infrared optical package referred to as Cast Glance. Approximately two months later, the same Navy Cast Glance aircraft successfully monitored the surface temperatures of the Orbiter Atlantis traveling at approximately Mach 14.3 during its return from the successful Hubble repair mission. In contrast to Discovery, Atlantis was not part of the Boundary Layer Transition (BLT) flight experiment, thus the vehicle was not configured with a protuberance on the port wing. In September 2009, Cast Glance was again successful in capturing infrared imagery and monitoring the surface temperatures on Discovery s next flight, STS-128. Again, NASA flew a specially modified thermal protection system tile and instrumentation package to monitor heating effects from boundary layer transition during re-entry. During this mission, Cast Glance was able to image laminar and turbulent flow phenomenology optimizing data collection for Mach 14.7. The purpose of this paper is to describe key elements associated with STS-119/125/128 mission planning and execution from the perspective of the Cast Glance flight crew that obtained the imagery. The paper will emphasize a human element of experience, expertise and adaptability seamlessly coupled with Cast Glance system and sensor technology required to manually collect the required imagery. Specific topics will include a near infrared (NIR) camera upgrade that was implemented just prior to the missions, how pre-flight radiance modeling was utilized to optimize the IR sensor configuration, communications, the development of aircraft test support positions based upon Shuttle trajectory information, support to contingencies such as Shuttle one orbit wave-offs/west coast diversions and then the Cast Glance perspective during an actual Shuttle imaging mission.

Tack, Steve↗

Exploring the Future Energy Value of Long-Duration Energy Storage

Long-duration energy storage is commonly viewed as a key technology for providing flexibility to the grid and broader energy systems over a multidecadal time frame. However, prior work has typically used present-day grid infrastructures to characterize the relationship between the duration and arbitrage value of storage in electricity markets. This study leverages established National Renewable Energy Laboratory grid planning and operations tools, analysis, and data to execute a price-taker model of an energy storage system for several 8760 h price series representative of current and future contiguous United States grid infrastructures with varying shares of variable renewable energy (VRE). We find that the total value of energy storage typically increases with VRE shares, but any increase in the relative value of longer storage durations over time depends on the region and grid mix. Some regions see incremental value increasing notably, up to 20–40 h in 2050, while others do not. The negative effect of lower roundtrip efficiency on value is also found to be scenario-dependent, with the energy value in higher VRE scenarios being less sensitive to roundtrip efficiency and more supportive of longer storage durations. Long-duration storage value and deployment potential are a function of evolving electricity sector infrastructure, markets, and policy, making it critical to consistently revisit potential long-duration storage contributions to the grid.

14 SOLAR ENERGY↗

Determination of Eros Physical Parameters for Near Earth Asteroid Rendezvous Orbit Phase Navigation

Navigation of the orbit phase of the Near Earth steroid Rendezvous (NEAR) mission will re,quire determination of certain physical parameters describing the size, shape, gravity field, attitude and inertial properties of Eros. Prior to launch, little was known about Eros except for its orbit which could be determined with high precision from ground based telescope observations. Radar bounce and light curve data provided a rough estimate of Eros shape and a fairly good estimate of the pole, prime meridian and spin rate. However, the determination of the NEAR spacecraft orbit requires a high precision model of Eros's physical parameters and the ground based data provides only marginal a priori information. Eros is the principal source of perturbations of the spacecraft's trajectory and the principal source of data for determining the orbit. The initial orbit determination strategy is therefore concerned with developing a precise model of Eros. The original plan for Eros orbital operations was to execute a series of rendezvous burns beginning on December 20,1998 and insert into a close Eros orbit in January 1999. As a result of an unplanned termination of the rendezvous burn on December 20, 1998, the NEAR spacecraft continued on its high velocity approach trajectory and passed within 3900 km of Eros on December 23, 1998. The planned rendezvous burn was delayed until January 3, 1999 which resulted in the spacecraft being placed on a trajectory that slowly returns to Eros with a subsequent delay of close Eros orbital operations until February 2001. The flyby of Eros provided a brief glimpse and allowed for a crude estimate of the pole, prime meridian and mass of Eros. More importantly for navigation, orbit determination software was executed in the landmark tracking mode to determine the spacecraft orbit and a preliminary shape and landmark data base has been obtained. The flyby also provided an opportunity to test orbit determination operational procedures that will be used in February of 2001. The initial attitude and spin rate of Eros, as well as estimates of reference landmark locations, are obtained from images of the asteroid. These initial estimates are used as a priori values for a more precise refinement of these parameters by the orbit determination software which combines optical measurements with Doppler tracking data to obtain solutions for the required parameters. As the spacecraft is maneuvered; closer to the asteroid, estimates of spacecraft state, asteroid attitude, solar pressure, landmark locations and Eros physical parameters including mass, moments of inertia and gravity harmonics are determined with increasing precision. The determination of the elements of the inertia tensor of the asteroid is critical to spacecraft orbit determination and prediction of the asteroid attitude. The moments of inertia about the principal axes are also of scientific interest since they provide some insight into the internal mass distribution. Determination of the principal axes moments of inertia will depend on observing free precession in the asteroid's attitude dynamics. Gravity harmonics are in themselves of interest to science. When compared with the asteroid shape, some insight may be obtained into Eros' internal structure. The location of the center of mass derived from the first degree harmonic coefficients give a direct indication of overall mass distribution. The second degree harmonic coefficients relate to the radial distribution of mass. Higher degree harmonics may be compared with surface features to gain additional insight into mass distribution. In this paper, estimates of Eros physical parameters obtained from the December 23,1998 flyby will be presented. This new knowledge will be applied to simplification of Eros orbital operations in February of 2001. The resulting revision to the orbit determination strategy will also be discussed.

Miller, J. K.↗

Keep it short: Exploring the impacts of configuration choices on the recent economics of solar-plus-battery and wind-plus-battery hybrid energy plants

Coupled renewable-battery powerplants differ from the traditional concept of independent siting of electricity resources within transmission networks. Prior research on the value proposition and cost savings from coupling did not consider the geographic constraint of co-location. This paper fills the gap by assessing how pricing volatility differences between nodes within electricity markets impact the system value of coupled renewable-battery projects as compared to independent VRE and battery installations. We use wholesale power market prices from 2012–2019 across the seven main U.S. independent system operators (ISOs) with a linear optimization program to compare the electricity market value of coupled projects to the value of the same underlying sub-components, deployed separately. We find that additional value from adding a 4-hour battery sized to 50% of renewable-plant nameplate capacity is $\$10$/MWh across ISOs on average. The highest boost occurs in California ($\$15$/MWh), where the value of adding storage to solar rises over time in tandem with increased solar penetration in the region. If renewables and batteries are deployed independently, we estimate that $\$12.5$/MWh of additional value could be achieved because of more flexibility on battery siting and operation. The $\$12.5$/MWh coupling penalty is reduced to $\$1.6$/MWh when considering alternative approaches to integrating battery storage. This result implies that renewable-battery power plants will play an increasing role in electricity systems if they can be built for $\$2$–$\$13$/MWh less than independent projects of comparable size. However, the wide regional variation in coupling penalties, along with the importance of conditions captured in our sensitivity cases, suggests the tradeoff between coupling penalties and savings will vary by situation. Therefore, roles exist for independent and coupled projects from a system optimization perspective.

14 SOLAR ENERGY↗

Are coupled renewable-battery power plants more valuable than independently sited installations?

Coupled renewable-battery powerplants differ from the traditional concept of independent siting of electricity resources within transmission networks. Prior research on the value proposition and cost savings from coupling did not consider the geographic constraint of co-location. This paper fills the gap by assessing how pricing volatility differences between nodes within electricity markets impact the system value of coupled renewable-battery projects as compared to independent VRE and battery installations. We use wholesale power market prices from 2012–2019 across the seven main U.S. independent system operators (ISOs) with a linear optimization program to compare the electricity market value of coupled projects to the value of the same underlying sub-components, deployed separately. We find that additional value from adding a 4-hour battery sized to 50% of renewable-plant nameplate capacity is $\$$10/MWh across ISOs on average. The highest boost occurs in California ($\$$15/MWh), where the value of adding storage to solar rises over time in tandem with increased solar penetration in the region. If renewables and batteries are deployed independently, we estimate that $\$$12.5/MWh of additional value could be achieved because of more flexibility on battery siting and operation. The $\$$12.5/MWh coupling penalty is reduced to $\$$1.6/MWh when considering alternative approaches to integrating battery storage. This result implies that renewable-battery power plants will play an increasing role in electricity systems if they can be built for $\$$2–$\$$13/MWh less than independent projects of comparable size. However, the wide regional variation in coupling penalties, along with the importance of conditions captured in our sensitivity cases, suggests the tradeoff between coupling penalties and savings will vary by situation. Therefore, roles exist for independent and coupled projects from a system optimization perspective.

14 SOLAR ENERGY↗

Planetary Protection Bioburden Analysis Program

This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware. This work was done by David Shores and Sharon P. Goza of Johnson Space Center; Cheyenne McKeegan, Rick Easley, Janet Way, and Shonn Everett of MEI Technologies; Mark Manning of PTI; and Mark Guerra, Ray Kraesig, and William Leu of Tietronix Software, Inc. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24211-1 Spitzer Telemetry Processing System NASA's Jet Propulsion Laboratory, Pasadena, California The Spitzer Telemetry Processing System (SirtfTlmProc) was designed to address objectives of JPL's Multi-mission Image Processing Lab (MIPL) in processing spacecraft telemetry and distributing the resulting data to the science community. To minimize costs and maximize operability, the software design focused on automated error recovery, performance, and information management. The system processes telemetry from the Spitzer spacecraft and delivers Level 0 products to the Spitzer Science Center. SirtfTlmProc is a unique system with automated error notification and recovery, with a real-time continuous service that can go quiescent after periods of inactivity. The software can process 2 GB of telemetry and deliver Level 0 science products to the end user in four hours. It provides analysis tools so the operator can manage the system and troubleshoot problems. It automates telemetry processing in order to reduce staffing costs. This work was done by Alice Stanboli, Elmain M. Martinez, and James M. McAuley of Caltech for NASA's Jet Propulsion Laboratory. For more information, contact iaoffice @jpl.nasa.gov. This software is available for commercial licensing. Please contact Dan Broderick at Daniel.F. Broderick@jpl.nasa.gov. Refer to NPO-47803. NASA Tech Briefs, September 2013 29 This rapid response computer program predicts Orbiter Wing Leading Edge (WLE) damage caused by ice or foam impact during a Space Shuttle launch (Program "IMPACT2"). The program was developed after the Columbia accident in order to assess quickly WLE damage due to ice, foam, or metal impact (if any) during a Shuttle launch. IMPACT2 simulates an impact event in a few minutes for foam impactors, and in seconds for ice and metal impactors. The damage criterion is derived from results obtained from one sophisticated commercial program, which requires hours to carry out simulations of the same impact events. The program was designed to run much faster than the commercial program with prediction of projectile threshold velocities within 10 to 15% of commercial-program values. The mathematical model involves coupling of Orbiter wing normal modes of vibration to nonlinear or linear springmass models. IMPACT2 solves nonlinear or linear impact problems using classical normal modes of vibration of a target, and nonlinear/ linear time-domain equations for the projectile. Impact loads and stresses developed in the target are computed as functions of time. This model is novel because of its speed of execution. A typical model of foam, or other projectile characterized by material nonlinearities, impacting an RCC panel is executed in minutes instead of hours needed by the commercial programs. Target damage due to impact can be assessed quickly, provided that target vibration modes and allowable stress are known. This work was done by Robert Clark, Jr., Paul Cotter, and Constantine Michalopoulos of The Boeing Company for Johnson Space Center. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24988-1 Wing Leading Edge RCC Rapid Response Damage Prediction Tool (IMPACT2) Lyndon B. Johnson Space Center, Houston, Texas (3) the MSL statistics for only the accountable assays. Other options on the main menu include a data editing form and utility programs that produce various reports requested by the microbiologists and the project, and tools to generate the groupings for the final analyses. The analyses can be carried out in three ways: Each assay can be treated separately, the assays can be collectively treated for the whole zone as a group, or the assays can be collected in groups designated by the JPL Planetary Protection Manager. The latter approach was used to generate the final report because assays on the same equipment or similar equipment can be assumed to have been exposed to the same environment and cleaning. Thus, the statistics are improved by having a larger population, thereby reducing the standard deviation by the square root of N. For each method mentioned above, three reports are available. The first is a detailed report including all the data. This version was very useful in verifying the calculations. The second is a brief report that is similar to the full detailed report, but does not print out the data. The third is a grand total and summary report in which each assay requires only one line. For the first and second reports, most of the calculations are performed in the report section itself. For the third, all the calculations are performed directly in the query bound to the report. All the numerical results were verified by comparing them with Excel templates, then exporting the data from the Planetary Protection Analysis program to Excel.

Beaudet, Robert A.↗

A Multi-Model Framework for Assessing Long- and Short-Term Climate Influences on the Electric Grid

Climate change influences many aspects of the electric grid, but prior work and industry practices often ignore the potential effects of changing climate, or they only consider a single effect or individual effects in isolation. Challenges vary with each grid and include adapting to long-term trends such as changing temperature and precipitation or shorter-term events such as drought or storms that could increase in frequency or intensity. Here we present a multi-model framework designed to analyze the effects of long and short-term climate impacts in combination. This framework couples capacity expansion and production cost models with hydrologic models and future climate scenario data to analyze alternative climate and energy futures at high spatial, temporal, and process resolutions. Furthermore, we constructed and evaluated the results of a suite of simulated scenarios exploring climate impacts on capacity investment and stress-tested the resulting future infrastructures using hourly dispatch modeling under alternative drought and load conditions. We demonstrate the approach through a case study of the U.S. Western Interconnection, where climate impacts depend on interactions between temperature-induced load, water availability for hydropower, technology competitiveness, and demand flexibility. Changes in 2038 generating capacity range from -8.5-16.6 GW, and changes in 2038 transmission capacity range from -1-2 GW. Capacity increases are driven by higher load from higher temperatures, while capacity reductions can be achieved in scenarios with higher future hydropower availability and increased demand flexibility. Scenarios requiring additional capacity cost an additional $\$5$-$\$17$ billion (discounted) from 2018 to 2038; however, scenarios with capacity reductions cost $\$1$-$\$18$ billion less. Stress tests on four 2038 infrastructures demonstrated that the identified systems were able to serve at least 99.999% of load and 99.96% of reserves. However, drought and unexpected high-load conditions can result in reduced capacity to respond to contingency events we did not model. Although these results are system and scenario specific, they highlight the importance of considering multiple climate change impacts simultaneously in long-term planning efforts and demonstrate a multi-model, multiscale approach that can be flexibly applied to any system and set of climate change concerns.

24 POWER TRANSMISSION AND DISTRIBUTION↗