Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable speed limit”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Scalability analysis of heavy-duty gas turbines using data-driven machine learning

With the increasing integration of variable renewable energy sources into power systems, the role of flexible power generation technologies like gas turbines (GT) in rapid grid balancing remains crucial. This sustained importance underscores the need for scaled and precise modeling of GT to ensure effective integration within evolving energy frameworks. While physics-driven GT models integrate thermodynamics, fluid dynamics, and combustion principles, they often rely on approximate mathematical representations to accommodate scaling that may not capture the actual complex dynamics for GTs and inertial effects associated to GTs with different ratings. In this study, a data-driven model is proposed using machine learning (ML) techniques to conduct GT scalability analysis and performance evaluation with high accuracy. The ML model, trained on data from various operating conditions and performance parameters, aims to uncover intricate relationships and patterns, resembling GT characteristics at different scales (ratings). The model is developed to capture complex system interaction and to adapt to changing operational scenarios at different capacities, providing valuable insights of power system dynamics. In this study, the real-time digital simulator platform was employed to generate training data for the ML model and assess its dynamic characteristics. The ultimate objective was to develop a detailed modeling framework based on governing equations and data-driven ML capable of predicting key performance indicators, in thermal systems such as GTs, including power output, speed, fuel consumption, and exhaust temperature under diverse operating conditions at different scales. The developed ML framework demonstrated high accuracy, with mean relative errors for GT power prediction, reference speed, exhaust temperature, and compressor pressure ratio (CPR) parameters consistently below 0.1% across typical load fluctuation scenarios. Maximum deviations were limited to approximately 0.5 K for exhaust temperature and 0.009 for CPR, underscoring the model’s ability to replicating dynamic GT behavior with high precision. The adaptability of the ML model enables its application across diverse operational conditions and its extension to other thermal systems. By leveraging advanced ML techniques, this study presents a robust and scalable modeling framework that enhances GT simulation precision, facilitating improved integration into evolving power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1999

Physical Control of Biological Productivity Off the Coast of Peru During the 1997-1998 El Nino

Satellite observations and an ecosystem model are used to understand the variability in the planktonic ecosystem off Peru for the period January 1996 to May 1998. The objective of this study is to quantify the changes in the ecosystem components, carbon pathways, and available food for small pelagic fish that occur associated with the change in physical forcing due to El Nino. Two periods are distinguished based on the observed sea level anomaly: a La Nina (LaN) period (1996) in which sea level was below normal and El Nino (EN), the average conditions for December 1997, in which the sea level was anomalously high. There are three phytoplankton size classes (pico-, nano-, and net-phytoplankton) which compete for nutrients and are eaten by three zooplankton size classes. The ecosystem model is forced by alongshore wind speed measured by the NASA Scatterometer (NSCAT) and the European Remote-sensing Satellites (ERS-1 and ERS-2). Larger, slower growing organisms are more sensitive to physical disturbance than smaller organisms (Carr, 1998]. In the present simulation as well, the primary effect of the El Nino (reduced nutrient supply, and increased temperature) is to reduce the biomass of large cells (netphytoplankton) and consequently of the zooplankton that rely on large cells as food source. EN conditions are accompanied by a rearrangement of carbon pathways: comparable uptake goes into reduced biomass accumulation, increased losses to respiration, reduced carbon export, and much reduced carbon available to fish. The star indicates the remotely sensed biomass (assuming a constant carbon to chlorophyll ratio of 60) as measured by the Ocean Color and Temperature Sensor (Nov.-Dec. 1996) and the Sea-viewing Wide Field-of-view Sensor (Dec. 1997). The model, which assumes no light limitation, overestimates total phytoplankton biomass. Additional Information is contained in the original.

Carr, Mary-Elena↗
NASA NTRS2013

RUSHMAPS: Real-Time Uploadable Spherical Harmonic Moment Analysis for Particle Spectrometers

RUSHMAPS is a new onboard data reduction scheme that gives real-time access to key science parameters (e.g. moments) of a class of heliophysics science and/or solar system exploration investigation that includes plasma particle spectrometers (PPS), but requires moments reporting (density, bulk-velocity, temperature, pressure, etc.) of higher-level quality, and tolerates a lowpass (variable quality) spectral representation of the corresponding particle velocity distributions, such that telemetry use is minimized. The proposed methodology trades access to the full-resolution velocity distribution data, saving on telemetry, for real-time access to both the moments and an adjustable-quality (increasing quality increases volume) spectral representation of distribution functions. Traditional onboard data storage and downlink bandwidth constraints severely limit PPS system functionality and drive cost, which, as a consequence, drives a limited data collection and lower angular energy and time resolution. This prototypical system exploit, using high-performance processing technology at GSFC (Goddard Space Flight Center), uses a SpaceCube and/or Maestro-type platform for processing. These processing platforms are currently being used on the International Space Station as a technology demonstration, and work is currently ongoing in a new onboard computation system for the Earth Science missions, but they have never been implemented in heliospheric science or solar system exploration missions. Preliminary analysis confirms that the targeted processor platforms possess the processing resources required for realtime application of these algorithms to the spectrometer data. SpaceCube platforms demonstrate that the target architecture possesses the sort of compact, low-mass/power, radiation-tolerant characteristics needed for flight. These high-performing hybrid systems embed unprecedented amounts of onboard processing power in the CPU (central processing unit), FPGAs (field programmable gate arrays), and DSP (digital signal processing) elements. The fundamental computational algorithm de constructs 3D velocity distributions in terms of spherical harmonic spectral coefficients (which are analogous to a Fourier sine-cosine decomposition), but uses instead spherical harmonics Legendre polynomial orthogonal functions as a basis for the expansion, portraying each 2D angular distribution at every energy or, geometrically, spherical speed-shell swept by the particle spectrometer. Optionally, these spherical harmonic spectral coefficients may be telemetered to the ground. These will provide a smoothed description of the velocity distribution function whose quality will depend on the number of coefficients determined. Successfully implemented on the GSFC-developed processor, the capability to integrate the proposed methodology with both heritage and anticipated future plasma particle spectrometer designs is demonstrated (with sufficiently detailed design analysis to advance TRL) to show specific science relevancy with future HSD (Heliophysics Science Division) solar-interplanetary, planetary missions, sounding rockets and/or CubeSat missions.

Figueroa-Vinas, Adolfo↗
DOE OSTI2025

Quality-Controlled Meteorological Data from the Flood Control District of Maricopa County (FCDMC) Network, Phoenix, Arizona (1987-2024)

This dataset contains 15- or 30-minute interval meteorological data from the Flood Control District of Maricopa County (FCDMC), Arizona, USA, covering eight key variables across multiple sensor stations between 1987 and 2024. Each variable is stored as a separate CSV file, containing time-series data that have undergone rigorous quality control (QC) procedures and, where appropriate, short-gap interpolation for consistency. The quality control (QC) pipeline consisted of four sequential tests: (1) a range test to ensure all values fall within physically realistic limits, (2) a step test to identify abrupt and implausible changes between consecutive records, (3) a proximity test that validates flagged values from step test using data from nearby stations and exceedance probability thresholds, and (4) a persistence test to detect and remove periods of unrealistically constant readings. These thresholds were calibrated to Arizona’s environmental conditions and sensor specifications. After QC, short gaps (≤2 hours) were linearly interpolated to ensure consistent temporal resolution, except for wind variables. Due to a major upgrade in FCDMC’s data transmission system, only ALERT-2 protocol data (2016–2024) for wind variables are included; earlier ALERT-1 data were excluded because of irregular sampling and high missing rates. This dataset supports regional climate and infrastructure resilience studies by providing standardized, high-resolution meteorological data for the greater Phoenix metropolitan area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Recurring Planetary Debris Transits and Circumstellar Gas around White Dwarf ZTF J0328–1219

We present follow-up photometry and spectroscopy of ZTF J0328-1219, strengthening its status as a white dwarf exhibiting transiting planetary debris. Using TESS and Zwicky Transient Facility photometry, along with follow-up high-speed photometry from various observatories, we find evidence for two significant periods of variability at 9.937 and 11.2 hr. We interpret these as most likely the orbital periods of different debris clumps. Changes in the detailed dip structures within the light curves are observed on nightly, weekly, and monthly timescales, reminiscent of the dynamic behavior observed in the first white dwarf discovered to harbor a disintegrating asteroid, WD 1145+017. We fit previously published spectroscopy along with broadband photometry to obtain new atmospheric parameters for the white dwarf, with M * = 0.731 ± 0.023 M⊙, T eff = 7630 ± 140 K, and [Ca/He] = - 9.55 ± 0.12. With new high-resolution spectroscopy, we detect prominent and narrow Na D absorption features likely of circumstellar origin, with velocities 21.4 ± 1.0 km s -1 blueshifted relative to atmospheric lines. We attribute the periodically modulated photometric signal to dusty effluents from small orbiting bodies such as asteroids or comets, but we are unable to identify the most likely material that is being sublimated, or otherwise ejected, as the environmental temperatures range from roughly 400 to 700 K.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

NNSS Soils Monitoring: Plutonium Valley (CAU 366) FY2019

Desert Research Institute (DRI) conducted a field assessment of the potential for radionuclide-contaminated soil to be transported from the Plutonium Valley Dispersion Sites Contamination Area (CA) because of both wind and storm water runoff. This activity supported U.S. Department of Energy (DOE) Environmental Management Nevada Program (EM NV) efforts to establish post-closure monitoring plans for the Plutonium Valley Dispersion Sites Corrective Action Unit (CAU) 366. The DRI task was intended to identify the likely mechanism(s) of transport and determine the meteorological conditions that might cause the movement of radionuclide-contaminated soils. The emphasis of the work was on collecting sediment transported by channelized storm runoff and measuring airborne dust concentrations and associated wind conditions. These data will facilitate an appropriate post-closure monitoring program. In 2011, DRI installed meteorological monitoring stations north and south of the Plutonium Valley CA—known as Pu Valley North and Pu Valley South, respectively—as well as a fluvial sediment sampling station within the CA. Since installation, temperature, wind speed, wind direction, relative humidity, precipitation, solar radiation, barometric pressure, soil temperature, volumetric soil moisture content, and airborne particulate concentrations have been collected at both meteorological stations. The maximum, minimum, and average or total (as appropriate) for each of these parameters were recorded for each 10-minute interval. The sediment sampling station included an automatically activated sampling pump with collection bottles for suspended sediment, which was activated when sufficient flow was present in the channel, and passive traps for bedload material that was transported down the channel during runoff events. This report presents data collected from these stations during fiscal year (FY) 2019. During the FY2019 (October 1, 2018, through September 30, 2019) reporting period, the warmest month was July and the coldest month was February. Monthly total solar radiation was highest in June or July (depending on the station) and lowest in December. Monthly mean wind speeds were highest during February. At the Pu Valley North station, winds were commonly northerly, northwesterly, and southerly. At the Pu Valley South station, the wind direction was more variable and winds frequently came from most directions, excluding westerly and northwesterly winds. Winds above 15.0 miles per hour (mph) (24.1 kilometers per hour [km/hr]) were frequently southerly at both stations, and less-frequent strong winds were northerly (at the Pu Valley South station) or northwesterly (at the Pu Valley North station). Monthly average relative humidity ranged from 18 percent in August to 62 percent in January and February. Monthly total precipitation ranged from zero at both stations in September to 2.23 inches (in) (56.6 millimeters [mm]) during February at the Pu Valley South station and 2.09 in (53.1 mm) during March at the Pu Valley North station. From October 1, 2018, through September 30, 2019, the total precipitation was 8.23 in (209 mm) and 8.15 in (207 mm) at the Pu Valley South and Pu Valley North stations, respectively. The largest daily precipitation totals of 1.46 in (37.1 mm) and 1.14 in (29.0 mm) occurred at the Pu Valley South and Pu Valley North stations, respectively, on March 6, 2019. However, the ISCO autosampler was not activated during any storm event in FY2019 because the channel water depth measurements were insufficient for suspended sediment sample collection. Additionally, because no significant flow through the channel was indicated, no bedload samples were collected. Light breezes of 0.0 mph to 5.0 mph (0.0 km/hr to 8.0 km/hr) occurred most frequently (approximately 55 percent to 57 percent of the time). The frequency of occurrence diminished exponentially as the wind speed increased such that winds in excess of 25.0 mph (40.2 km/hr) occurred less than 0.01 percent of the time. Winds in excess of 15.0 mph (24.1 km/hr) were most commonly southerly, which is likely controlled by the topography of the valley because topographic highs define the east and west sides of the valley and converge toward the north.The concentrations of PM2.5 (particulate matter with an aerodynamic diameter ≤2.5 micrometers [μm]) and PM10 (particulate matter with an aerodynamic diameter ≤10 μm) in the air generally increase as wind speed increases. However, the particle profiler at the Pu Valley South station regularly malfunctioned, which resulted in poor data quality and a limited period of record within FY2019. Therefore, discussion of the annual dust observations is mostly limited to the Pu Valley North station. At the Pu Valley North station, significant increases in windblown dust concentrations were observed when wind speeds exceeded 15.0 mph (24.1 km/hr), especially during southerly winds. When all wind speeds were considered, high dust concentrations at the Pu Valley North station were frequently associated with northerly, northwesterly, or southerly winds. The majority of PM10 transport at the Pu Valley North station occurred when winds were blowing northerly and northwesterly because of the high frequency of winds up to 15.0 mph (24.1 km/hr). The ratio of PM10 to PM2.5 is a qualitative indicator of the proximity of dust sources to the observation point and the values of the PM10 to PM2.5 ratio tend to be higher nearer to the source. The PM10 to PM2.5 ratio increased from 1.2 to 2.8 with higher wind speeds at the Pu Valley North station. This increase suggests that some locally sourced PM10 was present in the air when local winds exceeded 15.0 mph (24.1 km/hr). Local dust suspension is interpreted during major dust events when peaks in the PM10 to PM2.5 ratio occur simultaneously with the maximum wind speed and the peak PM10 concentration. Analysis of the five major dust events in FY2019 revealed that northerly winds were dominant in three events and southerly winds were dominant in two events. Local suspension of dust within Plutonium Valley is suspected to have occurred during two of the five major events (April 9, 2019, and September 28 to 29, 2019), and a nonlocal source of dust is presumed for the three other major dust events (October 14, 2018; May 16, 2019; and September 2, 2019). The three suspected regional dust transport events exhibited winds that typically did not exceed the 15.0 mph (24.1 km/hr) wind speed threshold for local dust transport and featured PM10 to PM2.5 ratios that indicate large-scale regional dust transport. Similar event characteristics (e.g., wind speed and PM10 concentration) were observed at air monitoring stations approximately 70 mi (113 km) to the northwest within the Tonopah Test Range (TTR). Therefore, these three regional events may not reflect transport from the Plutonium Valley CA. Routine monitoring of meteorological and hydrologic parameters at Plutonium Valley has been discontinued at the request of EM NV, effective at the conclusion of FY2019. Environmental data acquisition has been completed, and both Pu Valley North and Pu Valley South meteorological stations within Plutonium Valley as well as the ISCO autosampler station within the CA were removed on October 14, 2019.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Improving Computational Efficiency of Prognostics Algorithms in Resource-Constrained Settings

In engineering and aerospace applications, it is vital to operational success to have insight into the expected performance and health of physical systems. The field of prognostics and health management provides quantitative methods for monitoring, predicting, and managing system health. Prognostics algorithms can be employed to assess the current state of a system, propagate the system throughout time, and predict potential anomalies or failures that may occur. While they can provide accurate prediction results, effective prognostics algorithms can be challenging to use in resource-constrained settings due to computational limitations and high computational latency, leading to obsolete predictions. Thus, computationally efficient and accurate algorithms are necessary for future remaining useful life predictions. In this work, we implement new algorithmic approaches for prediction, quantitatively compare them via a battery degradation use-case, and provide recommendations of potential improvements to a prognostics framework. One approach to prediction is through sampling, whereby the current state of a physical system is sampled many times and each sample is propagated forward until failure is reached, resulting in a distribution of failure values. To improve the efficiency of this process, we implemented five new algorithmic approaches to prediction, including three distinct sampling methods (standard Monte Carlo, Quasi-Monte Carlo, and Latin Hypercube Sampling), a variable time step algorithm, and a variable sample size algorithm. To compare the algorithms, we employ a variety of metrics designed specifically to analyze both computational efficiency and model accuracy. Our metrics include accuracy to compare the average predicted value to ground truth, mean absolute deviation to illustrate dispersion, specific percentile error to describe accuracy within a user-defined risk tolerance, and code run-time. To quantitatively analyze our results, we employ a use-case of degradation of a Lithium-ion battery. We use an electrochemistry-based model to describe the current health state of the battery, and implement our prediction algorithms to propagate forward in time until end-of-discharge (EOD) is reached. Notably, through this work it was found that none of our sampling approaches had a significant impact on computational efficiency or model accuracy in predicting EOD of the battery. We find that while the sampling methods are unique, the distributions they generate are similar, ultimately producing final predictions that are nearly identical. In exploring the effect of the time step within the prediction algorithm, we found that prediction accuracy was highly dependent on the time step used, and that implementing a variable time step within a particular prediction may provide an increase in computational efficiency while also maintaining prediction accuracy. Finally, implementing a variable sample size also affected prediction, and our results show that tuning both the magnitude and timing of the sample size adjustment can result in improved computation speed and maintained prediction accuracy. Taken together, our findings highlight the challenge of performing prognostics in resource-constrained settings, and illustrate the potential of developing new prediction algorithms to improve computational efficiency.

prognostics↗
DOE OSTI2020

Macro-physical Properties of Shallow Cumulus from Integrated ARM Observations (Final Report)

Fair-weather shallow cumuli (ShCu) play an important role in many climate-related processes. Irregular geometry of ShCu and their strong temporal and spatial variability make it challenging to observe ShCu holistically and to represent them correctly in climate models. To improve ShCu parameterizations, information on both vertically and horizontally resolved cloud properties is required. Commonly, the vertically resolved cloud properties are provided by zenith pointing lidar-radar observations with a very narrow field of view (FOV). Thus, these “pencil-beam” properties may not be representative of a larger surrounding area. Limited number of areal-averaged cloud properties, such as fractional sky cover (FSC), are offered typically by wide-FOV observations. The main goal of our project was to integrate advantages of the narrow-FOV (vertical structure of clouds) and wide-FOV (spatial arrangement of clouds) observations for an improved characterization of single-layer ShCu observed at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site for an 18-yr period (2000-2017). There are four major accomplishments of our project, First, an updated operational cloud classification for days with ShCu has been suggested and evaluated through a detailed comparison with the manually curated records. Our classification extends successfully the latest ARM cloud type Value Added Product (VAP) based on the Active Remote Sensing of Clouds (ARSCL) cloud product by incorporating both cloud fraction (CF) provided by narrow-FOV ceilometer data and FSC from wide-FOV images offered by a Total Sky Imager (TSI). Moreover, our classification allows one to identify impact of instrumentation changes at the SGP site, namely the transition to KAZRARSCL with the updated cloud radar, on the identification of periods with single-layer ShCu. Second, a new approach that resolves cloud area distributions for a given region (up to 4x4 km 2 ) has been suggested and cloud equivalent diameters (CEDs) have been estimated for the first time. These estimations have been performed over a wide range of cloud sizes (about 0.01–3.5 km) with high temporal resolution (30s) using wide-FOV TSI images and cloud base height (CBH) provided by complementary narrow-FOV lidar measurements. Our simple and computationally inexpensive approach offers a previously unavailable dataset for process studies in the convective boundary layer and evaluation of ShCu parameterizations in cloud-resolving models. Third, a long-term integrated record of ShCu macrophysical properties has been developed. The developed record represents the longest available compilation of events with ShCu and includes (i) a novel visualization of the spatial variability in cloud cover both along- and across-wind directions, (ii) updated estimates of narrow-FOV CF and wide-FOV FSC, (iii) updated narrow-FOV CBH, and (iv) complementary data, such as wind speed and direction from the 915-MHz Radar Wind Profiler (RWP) data. The developed record has been used successfully to assess conventional observational estimates of cloud cover and their sensitivity to the following two factors: (i) instrument-dependent cloud detection and data merging criteria and (ii) FOV configuration. Fourth, co-variability of the ShCu macrophysical properties and environmental parameters has been analyzed for a 3-yr period (2016-2018). Our initial analysis includes diurnal changes of FSCs obtained for clouds with small, moderate and large CEDs and several environmental parameters, such as lifted condensation level (LCL) and mixed layer height (zi). Preliminary results of our analysis suggest that the horizontal extent of ShCu is controlled substantially by the sign and magnitude of difference between these two parameters (zi-LCL): the CED tends to grow with increase of this difference (zi exceeds LCL). We have initiated relationships between the ShCu and key atmospheric parameters that control both the development and evolution of ShCu using our new data product, which combines effectively the advantages of narrow-FOV data offered by zenith pointing cloud radars and lidars and wide-FOV TSI images. While the latest instrumentation at the ARM sites may address these challenging relationships in the future, we believe that the historical ARM data at the SGP site has not yet been fully utilized. Overall, our data product can be used by researchers working on a wide range of climate-related projects. These projects may include (i) a comprehensive evaluation of outputs from the Large-Eddy Simulation (LES) and single-column models for their future improvement, (ii) the representativeness of “short-period” results obtained from the previous model and observational studies and (iii) the planning of future field campaigns with focus on improved understanding of the diurnal cycle of cumulus convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Retrieval of the sea spray aerosol mode from submicron particle size distributions and supermicron scattering during LASIC

Abstract. Improved quantification of sea spray aerosol concentration and size is important for determining aerosol effects on clouds and the climate, though attempts to accurately capture the size distribution of the sea spray mode remain limited by the availability of supermicron size distributions. In this work, we introduce a new approach to retrieving lognormal mode fit parameters for a sea spray aerosol mode by combining submicron size distributions with supermicron scattering measurements using a Mie inversion. Submicron size distributions were measured by an ultra-high-sensitivity aerosol spectrometer (UHSAS), and supermicron scattering was taken as the difference between <10 µm and <1 µm three-wavelength integrating nephelometer measurements (NEPH). This UHSAS-NEPH method was applied during background marine periods of the Department of Energy Atmospheric Radiation Measurement Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign on Ascension Island (November 2016–May 2017), when the contribution of sea spray aerosol was expected to represent a large fraction of the aerosol mass and total scattering. Lognormal sea spray modal parameters were retrieved from comparisons between nephelometer measurements and a lookup table of Mie theory-simulated scattering coefficients for low-error solutions that minimized the 0.4–1 µm residual in the UHSAS size distribution. We evaluated the UHSAS-NEPH method with a set of clean marine measurements in the North Atlantic that included supermicron size and chemical measurements, showing that measured supermicron size distributions are needed to constrain the sea spray number concentration but that mass concentration was reasonably characterized using supermicron scattering. For LASIC, the UHSAS-NEPH method retrieved sea spray mode properties for approximately 88 % of the background marine times when the scattering variability and total particle concentration were low (<± 5 Mm−1 and <400 cm−3, respectively), with mass mean diameter ranging from 0.6 to 1.9 µm (1.47 ± 0.17 µm), modal width ranging from 1.1 to 3.97 (2.4±0.3), and mass concentration ranging from 0.18 to 23.0 µg m−3 (8.37. ± 4.1 µg m−3). The measured nephelometer scattering at three wavelengths was found to constrain the mode width marginally at the largest particle sizes in the absence of additional size and chemical measurements for defining parameters for the Mie solutions. Comparing UHSAS-NEPH retrievals to those of a fitting algorithm applied only to the submicron UHSAS number size distribution showed that correlations between retrieved mass concentration and the available mass-based sea spray tracers (coarse scattering, wind speed, and chloride) are low when supermicron measurements are not considered. This work demonstrates the added value of supermicron scattering measurements for retrieving reasonable sea spray mass concentrations, providing the best-available observationally constrained estimate of the sea spray mode properties when supermicron size distribution measurements are not available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Wound Field and Hybrid Synchronous Machines for EV Traction with Brushless Capacitive Rotor Field Excitation (Final Report)

This project focused on the development of wound field synchronous machines (WFSMs) and hybrid excitation synchronous machines (HESMs) with brushless capacitive power transfer for the field excitation. The target application for the machines developed is the main traction motor in electric vehicles. The magnetization in these types of machines is provided by a field winding on the rotor which is excited with DC current. The magnetization level in the machine can be varied by changing the magnitude of the field current. The variable magnetization or field is one of the key features of WFSMs. WFSMs are commonly used as generators however they have several attractive features for automotive traction applications. 1) No use of permanent magnet: Rare earth permanent magnets are primarily mined and processed in China. They have been subject to large price and supply variations and their export may be restricted during times of geopolitical tension. 2) Easy field weakening: Wound field synchronous machines have complete control of their field excitation. With proper design, this type of machine can electromagnetically have an infinite constant power speed range. 3) High power factor: With proper choice of the field excitation, WFSMs may be operated with high or even unity power factor. This potentially allows for the inverter connected to the stator winding to be downsized. In comparison induction machines and interior permanent magnet synchronous machines must supply reactive power to the stator increasing the kVA rating and cost of the inverter. 4) Reduced iron losses at high speed: By reducing the field excitation, the iron loss in the stator can be reduced. This is in comparison to interior permanent magnet synchronous machines which must use stator current to buck or reduce the flux produced by the permanent magnets. Generally, WFSMs have their highest efficiency at high speed. 5) Torque output at high temperatures: The magnetization provided by the field winding only depends on the field current and not on the field winding temperature. This is in contrast to permanent magnet machines where the permanent magnet flux decreases as the magnet temperature increases. Historically a number of approaches have been developed to provide DC current to the rotating field winding in WFSM’s including brushes and slip rings, low frequency brushless exciters, and high frequency rotary transformers and rectifiers. In this project, a different approach was used: brushless capacitive power transfer. Brushless capacitive power transfer uses two sets of rotating capacitors or electrodes in which an AC electric field is established by a high frequency inverter. A displacement current can flow through the airgap in the rotating capacitors which is rectified on the rotor using a diode bridge. The potential advantage of capacitive power transfer is that there is no need for heavy iron to guide magnetic flux. The electric flux lines terminate on the charges on the rotary capacitor surfaces. This should also limit the electric field outside the rotary capacitor airgaps. The main challenge with capacitive power transfer is that the capacitance and surface area of the rotating capacitors is small. Because capacitive power transfer systems are essentially a dual of a magnetic system, an Ampere per Hertz relationship is characteristic versus a Volts per Hertz relationship. A very high frequency power inverter must be used to provide sufficient excitation to the field winding. The concept of using capacitive power transfer to excite a high performance WFSM was initially developed in a previous U.S.A. Dept. of Energy project, DE-EE0006829. This project focused on increasing the power density of the WFSMs and reducing the cost and manufacturing complexity of the capacitive power transfer system. This project has demonstrated that WFSMs with brushless capacitive field power transfer can provide a high-power density and low-cost automotive powertrain technology.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1994

Technology Benefit Estimator (T/BEST): User's Manual

The Technology Benefit Estimator (T/BEST) system is a formal method to assess advanced technologies and quantify the benefit contributions for prioritization. T/BEST may be used to provide guidelines to identify and prioritize high payoff research areas, help manage research and limited resources, show the link between advanced concepts and the bottom line, i.e., accrued benefit and value, and to communicate credibly the benefits of research. The T/BEST software computer program is specifically designed to estimating benefits, and benefit sensitivities, of introducing new technologies into existing propulsion systems. Key engine cycle, structural, fluid, mission and cost analysis modules are used to provide a framework for interfacing with advanced technologies. An open-ended, modular approach is used to allow for modification and addition of both key and advanced technology modules. T/BEST has a hierarchical framework that yields varying levels of benefit estimation accuracy that are dependent on the degree of input detail available. This hierarchical feature permits rapid estimation of technology benefits even when the technology is at the conceptual stage. As knowledge of the technology details increases the accuracy of the benefit analysis increases. Included in T/BEST's framework are correlations developed from a statistical data base that is relied upon if there is insufficient information given in a particular area, e.g., fuel capacity or aircraft landing weight. Statistical predictions are not required if these data are specified in the mission requirements. The engine cycle, structural fluid, cost, noise, and emissions analyses interact with the default or user material and component libraries to yield estimates of specific global benefits: range, speed, thrust, capacity, component life, noise, emissions, specific fuel consumption, component and engine weights, pre-certification test, mission performance engine cost, direct operating cost, life cycle cost, manufacturing cost, development cost, risk, and development time. Currently, T/BEST operates on stand-alone or networked workstations, and uses a UNIX shell or script to control the operation of interfaced FORTRAN based analyses. T/BEST's interface structure works equally well with non-FORTRAN or mixed software analysis. This interface structure is designed to maintain the integrity of the expert's analyses by interfacing with expert's existing input and output files. Parameter input and output data (e.g., number of blades, hub diameters, etc.) are passed via T/BEST's neutral file, while copious data (e.g., finite element models, profiles, etc.) are passed via file pointers that point to the expert's analyses output files. In order to make the communications between the T/BEST's neutral file and attached analyses codes simple, only two software commands, PUT and GET, are required. This simplicity permits easy access to all input and output variables contained within the neutral file. Both public domain and proprietary analyses codes may be attached with a minimal amount of effort, while maintaining full data and analysis integrity, and security. T/BESt's sotware framework, status, beginner-to-expert operation, interface architecture, analysis module addition, and key analysis modules are discussed. Representative examples of T/BEST benefit analyses are shown.

Generazio, Edward R.↗
NASA NTRS2012

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.↗
NASA NTRS2019

Unmanned Aircraft Systems (UAS) Integration in the National Airspace System (NAS) Project: ACAS Xu HITL Overview

This presentation provides an overview of a recently completed human-in-the-loop simulation, conducted as part of the Unmanned Aircraft Systems (UAS) Integration in the National Airspace System (NAS) Project. This study examined how to present resolution advisories (RAs) issued by Airborne Collision Avoidance System (ACAS) Xu, which can be in the vertical dimension, the horizontal dimension, or both (i.e., "blended"). The study varied the location of the ACAS Xu traffic information - it was either presented within an 'integrated' display or a 'standalone' display, where the traffic information was separated from the vehicle control interfaces and navigational information. The presentation covers the test set-up, scenario design, and facilities utilized in this experiment. The primary dependent variables are listed, but no detailed results are included in this outbrief since data analysis is still underway. Goal: assess ACAS Xu Run 5 in a HITL setting to measure pilot & system performance – RWC alerting and guidance – RA alerting and guidance • Horizontal-only, vertical-only, and blended RAs • Strengthening RAs, reversals and added/removed RAs during blended maneuvers – Incorporate realistic sensor noise into the simulation environment – Compare pilot and DAA system performance back to NASA’s Phase 1 DAA MOPS V&V HITL • We leveraged the findings of an engineering analysis conducted in March to help inform how to implement the Xu logic – Part 1 focused on how to display horizontal-only, vertical-only and blended RAs – Part 2 focused on how to display automated RA responses Xu Engineering Analysis Results Part 1 objectives: – Characterize pilot responses to (canned) ACAS Xu RAs in a variety of display configurations • With vs. Without text accompaniment • Simple vs. ‘Advanced’ aural alerting • Results – Pilots struggled to meet 5 seconds initial response requirement • Particularly against horizontal and blended RAs – Pilots failed to respond more quickly to secondary RAs • Expected response time to subsequent RAs = 2.5 seconds No clear effect of the different alerting conditions on response times – Response times primarily driven by RA type • Pilots demonstrated high level of compliance with RAs and self-reported Xu alerting and guidance as being acceptable – 1/5 reported text as being necessary – 4/5 wanted the aural alert to retain the original RA sequence (i.e., issue a follow-on RA second rather than first) – 4/5 wanted to retain the “Maintain Heading/Vertical Speed” aural alert in the event that the pilot reached their target response at the time of secondary RA Based on these results, we decided on the following for the HITL experimental design: – Incorporated an ‘auto-fill’ feature in Vigilant Spirit that removes the need for pilots to manually enter a heading or altitude for RAs • Pilots only have to click ‘Send’ to upload the RA target heading/altitude – Did not include an RA ‘text box’ • The auto-fill feature largely replaces the purpose of the text box – Using a combination of the ‘Basic’ and ‘Advanced’ aural alerting • In case of blended maneuvers will issue “Maintain Heading/Vertical Rate” if pilot has reached 1st RA target at the time the 2nd RA is issued • Integration and testing with Xu Run 5 also resulted in us adding display logic to modify how horizontal RAs are presented – Target heading was shown to update at approx. 1 Hz making it difficult for pilots to implement – Display logic limited horizontal RA strengthening to once every 5 seconds • Did not impact timing of reversals, new RAs (i.e., blended), or CoC Engineering Analysis Results

human systems integration↗
NASA NTRSDate not supplied

Urban Air Mobility Noise: Current Practice, Gaps, and Recommendations

AN Air Mobility (UAM) is an opportunity for aviation to improve transportation systems across the world. Representative UAM vehicle attributes include electrical vertical takeoff and landing (eVTOL) vehicles that can accommodate up to 6 passengers (or equivalent cargo), are possibly autonomous, perform missions of up to 100 nautical miles at altitudes up to 3000 ft. above ground level, have flight speeds up to 200 knots, and weigh between 800 and 8000 pounds. Along with the many anticipated benefits, there will be noise issues that need to be addressed. In 2018, NASA formed an Urban Air Mobility Noise Working Group (UNWG) to assemble noise experts from industry, universities and government agencies to identify, discuss, and address UAM noise issues. This oral presentation summarizes technology gaps and goals associated with four areas of interest: Tools & Technologies, Ground & Flight Testing, Human Response & Metrics, and Regulation & Policy, and is drawn from a draft white paper [1] by the same title. Tools & Technologies include noise prediction tools and noise reduction technologies that have been developed for conventional rotorcraft and fixed-wing vehicles that may be applicable or need to be modified for UAM. Prediction tools need to be able to account for variable speed rotors and other temporal variation effects that impact community noise. A reprioritization of noise sources needs to be done since UAM vehicles include multiple rotors/propellers, often in proximity to one another and/or the airframe, with dynamic transition, and new noise sources such as electric motors or hybrid-electric propulsion. Scattering and propagation methods need to be developed that include the vehicle components and surfaces near a receiver such as buildings and vertiports. Validation databases are needed to quantify prediction uncertainties. Prediction tools used to evaluate community noise will need source models appropriate for a wide range of UAM vehicles. Existing noise reduction technologies need to be evaluated and new noise reduction technologies should be developed in anticipation of future noise requirements. Although the prediction and treatment of interior cabin noise is a secondary goal, it is recognized that new tools and methods may be needed due to the uniqueness of the vehicle design and the presence of both acoustic and structure-borne loads. Ground & Flight Testing has been a critical part of validating noise reduction technologies and verifying that an air vehicle is ready for certification. UAM vehicles introduce new challenges for test procedures such as different source noise directivity, unsteady sources due to maneuvers, and a variety of takeoff and approach trajectories. The operating environment will be more complex than current aircraft with the introduction of vertiports in populated areas with “urban canyons” making reflections an important part of noise prediction and annoyance. It is expected that new test procedures and measurement methods will be necessary. Consideration will need to be given for both piloted and autonomous operations. Human Response & Metrics may be very different for UAM noise compared to current experience with airport noise. Current metrics used to certify rotorcraft and fixed-wing aircraft may not be as useful for evaluating UAM noise. Operations at lower altitudes may influence annoyance. Psychoacoustic and community testing will be needed to quantify annoyance and assess appropriate metrics. In addition to conventional noise level metrics, considerations such as audibility and temporal variation of the sound may be required. Differences between indoor or outdoor exposure will have an impact with dependence on urban and residential flight paths. Aircraft noise is currently regulated at a national level and typically involves partnerships with the industry to establish regulations. Regulators realize that current policies and procedures may not be appropriate for some of the emerging air vehicles and new procedures may be needed to address UAM noise. Development of new policies and procedures are needed so that local communities do not hastily attempt to establish their own restrictions that will both limit growth of the market and create an inconsistent and confusing regulatory environment. To expedite this development, it is crucial early measurement data are shared through partnership arrangements to support both noise certification and noise modeling/noise assessment. At the same time, an effective engagement strategy should be developed to address local community noise issues associated with UAM vehicles and flight operations as they arise.

urban air mobility↗
NASA NTRS1996

Response of the neuromuscular unit to spaceflight: what has been learned from the rat model

Despite the inherent limitations placed on spaceflight investigations, much has been learned about the adaptations of the neuromuscular system to weightlessness from studies of rats flown for relatively short periods (approximately 4-22 days). Below is a summary of the major effects of spaceflight observed in muscles of rats that are not in their rapid growth stage: 1. Skeletal muscles atrophy rapidly during spaceflight; significant atrophy is observed as early as after 4 days of flight. 2. The atrophic response appears to be related to the primary function of the muscle. In the hindlimb, the relative amount of atrophy can be characterized as slow extensors > fast extensors > fast flexors. This pattern of relative atrophy does not appear to be occurring in the forelimb; however, not enough data are available to draw any definitive conclusions at this time. 3. Both slow and fast fibers atrophy during spaceflight, with the largest fibers within an individual muscle generally showing the greatest atrophic response. Interestingly, the amount of fiber atrophy appears to reach a plateau after about 14 days of flight. 4. Adaptations have been observed in the concentration and content of all muscle proteins pools, with the protein pools in slow muscles the most affected. 5. Some slow and fast fibers in predominantly slow and fast muscles show shifts in their histochemical and biochemical properties, toward those observed in a "faster" phenotype. 6. Some fibers, presumably expressing slow MHC isoforms before flight, begin to express fast MHC isoforms during flight. 7. The oxidative capacity of the muscles or fibers is relatively unaffected by spaceflight, particularly in the slow muscles. Any change in whole-body fatigability associated with spaceflight most likely reflects the loss in muscle and fiber mass. 8. The glycolytic capacity of the muscles and muscle fibers is enhanced after spaceflight. This metabolic adaptation seems to be related to the shift in the contractile proteins towards "faster" isoforms. 9. The vascularity of muscles appears to be maintained after flight, based, at least, on histological observations of capillarity. 10. The force capabilities of the muscles and fibers appear to decrease in parallel with the decreases in size, i.e., the specific tension is not significantly affected after flight. 11. Changes in the speed-related properties of the slow muscles are consistent with the adaptations in the myosin molecule, i.e., the slow muscles and some fibers in the slow muscles become "faster." 12. Some muscle fiber and neuromuscular junction damage has been observed after flight, particularly in the slow muscles. The extent of damage may be related to the amount of time that the muscles are allowed to reload before removal, i.e., in general, shorter intervals result in less fiber damage. 13. Adaptations in the motor (ventral horn) and sensory (spinal ganglia) neurons have been quite variable, but this may be related to the amount of time that the muscles are allowed to reload before removal. Morphological adaptations after relatively short periods of reloading may reflect a decrease in the activation of the neural elements during flight.

manned↗
NASA NTRSDate not supplied

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗
NASA NTRSDate not supplied

Global Reference Atmospheric Model (GRAM) Advancements and Additions

Introduction:The Global Reference Atmospheric Model (GRAM) is one of the most widely used engineering models of the atmosphere. GRAM development and maintenance has been led by NASA Marshall Space Flight Center (MSFC). The NASA Science Mission Directorate (SMD) has provided funding support to upgrade the GRAMs since Fiscal Year 2018. NASA Langley Research Center has been working with MSFC on the upgrades.This presentation will provide details regarding the upgrades that have been made to the existing GRAMs, the development of new GRAMs, as well as the ongoing objectives, tasks, and milestones re-lated to the GRAM upgrades funded by NASA SMD. GRAM: The GRAMs are engineering-oriented atmospheric models that estimate mean values and statistical variations of the atmospheric properties for numerous planetary destinations.They provide mean values and variability for any point in the atmosphere as well as seasonal, geographic, and altitude variations. GRAM outputs include atmospheric density, temperature, pressure, winds, and chemical composition along a user-defined path.Theyhave been widely used by the engineering community because of their ability to create realistic dispersions. GRAMs have been integrated into high fidelity flight dynamic simulations of launch, entry, descent and landing (EDL), aerobraking and aerocapture. MSFC has been developing and updating GRAMs since 1974; GRAMs are currently available for Earth, Mars, Venus, Neptune, and Titan. GRAM Upgrade Status: Code Moderization.The planetary GRAMs are being rearchitected from Fortran to a common object-ori-entedC++ framework called the GRAM Suite. This new architecture creates a common GRAM library of data models and utilities. The first C++ releases of the existing planetary GRAMs (Mars, Venus, Neptune, and Titan-GRAM) in the GRAM Suite are straight conversions from the latest Fortran version. Model Upgrades.The focus of the model upgrade task is to improve the atmosphere models in the existing GRAMs and to establish a foundation for developing GRAMs for additional destinations. The GRAM ephemerishas been upgraded to the NASA Navigation and Ancillary Information Facility (NAIF) SPICE toolkit (version N0066). The calculation of the speed of sound has also been improved in the GRAMs. The GRAM team has received updated Mars General Circulation Model (MGCM) datasets from NASA Ames Research Center. Mars Global Ionosphere-Thermosphere Model (M-GITM) data is being obtained toreplace the Mars Thermospheric General Circulation Model (MTGCM) data in legacy Mars-GRAM. M-GITM and updated MGCM data will be incorporated into a future GRAM Suite release.Twoprojects that will improve the atmospheric model data in the GRAMs have been funded by the GRAM team since Fiscal Year 2020. Sanjay Limaye and Patrick Fry at the University of Wisconsin are reanalyzing the Venus Express radio occultation observations and analyzing the Akatsukiradio occultation observations.This will lead tothe calculation of number density, temperature, and pressure profiles for the 40-90 km altitude range. Kunio Sayanagi, Justin Garland, and Ryan McCabeat Hampton University are developing empirical global models for Venus,Jupiter, Saturn, Uranus, Neptune, and Titan that incorporates the latest data available for each of these planetary destinations. Upgraded GRAM Releases. GRAM Suite Version 1.0 was released in May 2020 and contains the rearchitected Neptune-GRAM, including the common GRAM framework and planet–specific code. GRAM Suite Version 1.1 was released in September 2020 and addedthe rearchitected Titan-GRAM to the GRAM Suite. A User Guide and Programmer’s Manualarereleased with all GRAMs. The rearchitected Mars and Venus-GRAMs will be released in upcomingversions of the GRAM Suite. New GRAM Releases. New GRAMs have been de-veloped for Uranus and Jupiter. Uranus-GRAM is based on an individual profile generated by Gary Allen (ARC) from Voyager 2 occultation data and will be released in GRAM Suite Version 1.2. Jupiter-GRAM is based on individual profile produced from Al Seiff’s Ju-piter model[1].Jupiter-GRAM will be released in GRAM Suite Version 1.3. Saturn-GRAM is currently under development and will be released in a future version of the GRAM Suite. Conclusions: GRAMs are vital and frequently used toolsets. Releases of the GRAM Suite, upgrades of the existing planetary GRAMs, and development of new planetary GRAMs are ongoing. NASA SMD funding has been essential to addressing current limitations and accomplishing GRAM developmental goals. Continua-tion ofSMD fundingwill ensure the development, up-grades, and maintenance of the GRAMs. References: [1] Seiff, A., et al. (1998) JGR, 103, 22,857-22,889. Acknowledgments: The authors gratefully acknowledge support from the NASA SMD

atmospheric models↗
NASA NTRS2011

Solar Dynamics Observatory Guidance, Navigation, and Control System Overview

The Solar Dynamics Observatory (SDO) was designed and built at the Goddard Space Flight Center, launched from Cape Canaveral on February 11, 2010, and reached its final geosynchronous science orbit on March 16, 2010. The purpose of SDO is to observe the Sun and continuously relay data to a dedicated ground station. SDO remains Sun-pointing throughout most of its mission for the instruments to take measurements of the Sun. The SDO attitude control system (ACS) is a single-fault tolerant design. Its fully redundant attitude sensor complement includes sixteen coarse Sun sensors (CSSs), a digital Sun sensor (DSS), three two-axis inertial reference units (IRUs), and two star trackers (STs). The ACS also makes use of the four guide telescopes included as a part of one of the science instruments. Attitude actuation is performed using four reaction wheels assemblies (RWAs) and eight thrusters, with a single main engine used to provide velocity-change thrust for orbit raising. The attitude control software has five nominal control modes, three wheel-based modes and two thruster-based modes. A wheel-based Safehold running in the attitude control electronics box improves the robustness of the system as a whole. All six modes are designed on the same basic proportional-integral-derivative attitude error structure, with more robust modes setting their integral gains to zero. This paper details the final overall design of the SDO guidance, navigation, and control (GN&C) system and how it was used in practice during SDO launch, commissioning, and nominal operations. This overview will include the ACS control modes, attitude determination and sensor calibration, the high gain antenna (HGA) calibration, and jitter mitigation operation. The Solar Dynamics Observatory mission is part of the NASA Living With a Star program, which seeks to understand the changing Sun and its effects on the Solar System, life, and society. To this end, the SDO spacecraft carries three Sun-observing instruments: Helioseismic and Magnetic Imager (HMI), led by Stanford University; Atmospheric Imaging Assembly (AIA), led by Lockheed Martin Space and Astrophysics Laboratory; and Extreme Ultraviolet Variability Experiment (EVE), led by the University of Colorado. The basic mission is to observe the Sun for a very high percentage of the 5-year mission (10-year goal) with long stretches of uninterrupted observations and with constant, high-data-rate transmission to a dedicated ground station to be located in White Sands, New Mexico. These goals guided the design of the spacecraft bus that will carry and service the three-instrument payload. Overarching design goals for the bus are geosynchronous orbit, near-constant Sun observations with the ability to fly through eclipses, and constant HGA contact with the dedicated ground station. A three-axis stabilized ACS is needed both to point at the Sun accurately and to keep the roll about the Sun vector correctly positioned with respect to the solar north pole. This roll control is especially important for the magnetic field imaging of HM I. The mission requirements have several general impacts on the ACS design. Both the AIA and HMI instruments are very sensitive to the blurring caused by jitter. Each has an image stabilization system (ISS) with some ability to filter out high frequency motion, but below the bandwidth of the ISS the control system must compensate for disturbances within the ACS bandwidth or avoid exciting jitter at higher frequencies. Within the ACS bandwidth, the control requirement imposed by AIA is to place the center of the solar disk no more than 2 arc sec, 3 , from a body-defined target based on one of the GTs that accompany the instrument. This body-defined target, called the science reference boresight (SRB), was determined from the postlaunch orientation of the GTs by averaging the bounding telescope boresights for pitch to get a pitch SRB coordinate, and by averaging the bounding boresights for yaw toet the yaw SRB coordinate. The location of this SRB in the 0.5-deg field-of-view for each GT then becomes the central target for each telescope; one GT is selected for use as the ACS controlling guide telescope (CGT) at any given time. Fine Sun-pointing is effected based on this SRB for all three instruments when the Sun is within the linear range of the CGT. In addition to limiting jitter, HMI science requires averaging several observations, making the instrument sensitive to low frequency motion that induces differential motion between each observation. This requires the spacecraft attitude to be stable about the roll axis to approximately 10 arcsec over a ten-minute period. Instrument calibrations require that the spacecraft point the SRB up to 2.5 degrees in pitch and yaw away from the center of the Sun, placing the Sun outside the field-of-view of the guide telescopes. In such instances, when the GTs cannot provide the definitive target for the ACS, on-board attitude determination combined with ephemeris prediction of the Sun direction must provide the definitive target. EVE is capable of observing the Sun with less dependence on attitude control. However, the ground data processing needs for calibrations result in the most strict attitude knowledge requirements for the mission: [35,70,70] arcsec, 3 , of knowledge with respect to the center of the solar disk. In addition to driving the ACS sensor selection, the knowledge requirements, which have their effect primarily during Inertial mode calibrations, drive the accuracy requirements for the solar ephemeris. The need to achieve and maintain geosynchronous orbit (GEO) drove the need for high-efficiency propulsive systems and appropriate attitude control. The main engine provided high specific impulse for the maneuvers to attain GEO, while the smaller ACS thrusters managed the disturbance torques of the larger engine and provided the capability for much smaller adjustment burns on orbit. SDO s large solar profile means that solar radiation pressure is a large torque disturbance, and the momentum buildup from this disturbance and the GEO altitude drives the ACS to use thrusters to manage vehicle momentum. The demanding data capture budget for the mission, however, requires SDO to avoid frequent thruster maneuvers, while concerns about on-orbit jitter restrict the maximum desired wheel speeds desired from the RWAs. The plan for on-orbit wheel speed and momentum management will be discussed as well as what is now being done in operation after the jitter environment was characterized. The SDO ACS hardware complement is single-fault tolerant. Two main processors carry virtually identical copies of the command and data handling and ACS software, and two identical attitude control electronics (ACE) boxes carry Coldfire processors with contingency ACS software and other hardware interface cards; the ACE structure allows reaction wheels to be commanded by the Sun-pointing Safehold independent of the Mil Std 1553 data bus. The sixteen Adcole CSSs are grouped into primary and backup sets of eight sensors, each set providing the ability to calculate a sun vector. Each set of eight eyes provides full 4 -steradian coverage. The Adcole DSS comprises an optics head and a separate electronics box providing a 1553 data interface. The electronics box is mounted inside the Faraday cage created by the spacecraft bus module. The DSS head with its 32- deg square FOV is mounted on the instrument module with its boresight along the spacecraft X axis, nearly aligned with the Sun during observations. Adcole has designed the DSS calibration parameters so that the accuracy is 0.24 arcminutes within 10 deg of the boresight, and diminishes to 3 arcminutes as the Sun moves towards the edges of its FOV . This DSS calibration scheme provides higher accuracy attitude determination over the range of the instrument calibration maneuvers.

Morgenstern, Wendy M.↗