Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level Set Method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

The NASA Continuous Risk Management Process

As an intern this summer in the GRC Risk Management Office, I have become familiar with the NASA Continuous Risk Management Process. In this process, risk is considered in terms of the probability that an undesired event will occur and the impact of the event, should it occur (ref., NASA-NPG: 7120.5). Risk management belongs in every part of every project and should be ongoing from start to finish. Another key point is that a risk is not a problem until it has happened. With that in mind, there is a six step cycle for continuous risk management that prevents risks from becoming problems. The steps are: identify, analyze, plan, track, control, and communicate & document. Incorporated in the first step are several methods to identify risks such as brainstorming and using lessons learned. Once a risk is identified, a risk statement is made on a risk information sheet consisting of a single condition and one or more consequences. There can also be a context section where the risk is explained in more detail. Additionally there are three main goals of analyzing a risk, which are evaluate, classify, and prioritize. Here is where a value is given to the attributes of a risk &e., probability, impact, and timeframe) based on a multi-level classification system (e.g., low, medium, high). It is important to keep in mind that the definitions of these levels are probably different for each project. Furthermore the risks can be combined into groups. Then, the risks are prioritized to see what risk is necessary to mitigate first. After the risks are analyzed, a plan is made to mitigate as many risks as feasible. Each risk should be assigned to someone in the project with knowledge in the area of the risk. Then the possible approaches to choose from are: research, accept, watch, or mitigate. Next, all risks, mitigated or not, are tracked either individually or in groups. As the plan is executed, risks are re-evaluated, and the attribute values are adjusted as necessary. Metrics are established and monitored as tools for risk tracking. Also a trigger or threshold should be set on the metric data that indicates when an action is needed. Results of this tracking are usually evaluated and reported in a relevant format at weekly or monthly meetings. Choosing controls is the subsequent step, which involves the effects of the tracking. The three basic controls are: close, continue tracking, and re- plan. Finally communicate & document is the last step, but occurs throughout the process. It is vital that main risks, plans, changes, and progress are known by everyone in the project. A good way to keep everyone updated and inform other projects of common issues is by thoroughly documenting project risks. NASA sees value in risk management and believes that projects have greater probability or success by using the NASA Continuous Risk Management Process.

Pokorny, Frank M.↗

The SunPy Project: An Interoperable Ecosystem for Solar Data Analysis

The SunPy Project is a community of scientists and software developers creating an ecosystem of Python packages for solar physics. The project includes the sunpy core package as well as a set of affiliated packages. The sunpy core package provides general purpose tools to access data from different providers, read image and time series data, and transform between commonly used coordinate systems. Affiliated packages perform more specialized tasks that do not fall within the more general scope of the sunpy core package. In this article, we give a high-level overview of the SunPy Project, how it is broader than the sunpy core package, and how the project curates and fosters the affiliated package system. We demonstrate how components of the SunPy ecosystem, including sunpy and several affiliated packages, work together to enable multi-instrument data analysis workflows. We also describe members of the SunPy Project and how the project interacts with the wider solar physics and scientific Python communities. Finally, we discuss the future direction and priorities of the SunPy Project.

Solar physics↗

Ionospheric calibration for single frequency altimeter measurements

This report investigates the potential of using Global Positioning System (GPS) data and a model of the ionosphere to supply a measure of the sub-satellite Total Electron Current (TEC) of the required accuracy (10 TECU rms) for the purpose of calibrating single frequency radar altimeter measurements. Since climatological (monthly mean) models are known to be in error by as much as 50 percent, this work focused on the Parameterized Real-Time Ionospheric Specification Model (PRISM) which has the capability to improve model accuracy by ingesting (adjusting to) in situ ionospheric measurements. A set of globally distributed TEC measurements were generated using GPS data and were used as input to improve the accuracy of the PRISM model. The adjusted PRISM TEC values were compared to TOPEX dual frequency TEC measurements (which are considered truth) for a number of TOPEX sub-satellite tracks. The adjusted PRISM values generally compared to the TOPEX measurements within the 10 TECU accuracy requirements when the sub-satellite track passed within 300 to 400 km of the GPS TEC data or when the track passed through a night time ionosphere. However, when the sub-satellite points were greater than 300 to 400 km away from the GPS TEC data or when a local noon ionosphere was sampled, the adjusted PRISM values generally differed by greater than 10 TECU rms with data excursions from the TOPEX TEC measurements of as much as 40 TECU (an 8 cm path delay error at K band). Therefore, it can be concluded from this analysis that an unrealistically large number of GPS stations would be needed to predict sub-satellite TEC at the 10 TECU level in the day time ionosphere using a model such as PRISM. However, a technique currently being studied at the Jet Propulsion Laboratory (JPL) may provide a means of supplying adequate TEC data to meet the 10 TECU ionospheric correction accuracy when using a realistic number of ionospheric stations. This method involves using global GPS TEC data to estimate a global grid of vertical ionospheric TEC as a function of time (i.e. every one half hour) in a sun-fixed longitude frame. Working in a sun-fixed longitude frame, one is not limited by the spatial decorrelation distance of the ionosphere, but instead is limited more by the temporal correlations of the ionosphere in the sun-fixed frame which are a smaller effect. It is the opinion of the authors that using the global sun-fixed TEC grid data, in particular, ingesting it into PRISM, offers the best possibility of meeting the 10 TECU ionospheric correction accuracy requirement, and should be the subject of further study.

Schreiner, William S.↗

A compact solution to computational acoustics

This paper demonstrates that the linearized, dimensional Euler equations for acoustic computation can be accurately solved as a set of decoupled first-order wave equations, and that if ordered properly, this system of simple waves has unambiguous, easily implemented boundary conditions, allowing waves of same group speeds to pass through numerical boundaries or comply with wall conditions. Thus, the task of designing a complex multi-dimensional scheme with approximate far-field boundary conditions reduces to the design of higher order schemes for the one-dimensional simple wave equation. A compact finite-difference scheme and a characteristically exact but numerically n(th) order accurate boundary condition are introduced for solving the first order wave equation. Spanning a three-point two-level stencil, this low-dispersion implicit scheme has a third order spatial accuracy when used on nonuniform meshes, fourth order accurate on uniform meshes, and a temporal accuracy of second order due to the choice of trapezoidal integration for algorithmic simplicity. The robustness and accuracy of the scheme are demonstrated through a series of numerical experiments and comparisons with published results. When tested on the one-dimensional wave equation on a uniform grid, this scheme allows a Gaussian wave packet to pass through any finite domain with low numerical dispersion characteristic of a spatially fourth-order scheme and reflections at numerical boundaries maintained below truncation error. On highly stretched and irregular grids, only mild dispersions are found in the solution while solutions by other methods fail or are severely distorted. Yet, this scheme is no more sophisticated to solve or implement than the Crank-Nicolson scheme. This scheme has been tested on four categories of the ICASE/LaRC benchmark problems, which include propagation of acoustic and convective waves in Cartesian and cylindrical domains, reflection of acoustic wave at stationary/moving boundaries, and sound generation by gust-blade interaction.

Fung, K.-Y.↗

Origin and Reactivity of the Martian Soil: A 2003 Micromission

The role of water in the development of the martian surface remains a fundamental scientific question. Did Mars have one or more "warm and wet" climatic episodes where liquid water was stable at the surface? If so, the mineral phases present in the soils should be consistent with a history of aqueous weathering. More generally, the formation of hydrated mineral phases on Mars is a strong indicator of past habitable surface environments. The primary purpose of this investigation is to help resolve the question of whether such aqueous indicators are present on Mars by probing the upper meter for diagnostic mineral species. According to Burns [1993], the formation of the ferric oxides responsible for the visible color of Mars are the result of dissolution of Fe (+2) phases from basalts followed by aqueous oxidation and precipitation of Fe" mineral assemblages. These precipitates likely included iron oxyhydroxides such as goethite (a-FeOOH) and lepidocrocite (g-FeOOH), but convincing evidence for these phases at the surface is still absent. The stability of these minerals is enhanced beneath the surface, and thus we propose a subsurface search for hydroxylated iron species as a test for a large-scale chemical weathering process based on interactions with liquid water. It is also possible that the ferric minerals on Mars are not aqueous alteration products of the rocks. A chemical study of the Pathfinder landing site concluded that the soils are not directly derived from the surrounding rocks and are enhanced in Mg and Fe. The additional source of these elements might be from other regions of Mars and transported by winds, or alternatively, from exogenic sources. Gibson [1970] proposed that the spectral reflectivity of Mars is consistent with oxidized meteoritic material. Yen and Murray [1998] further extend Gibson's idea and show, in the laboratory, that metallic iron can be readily oxidized to maghemite and hematite under present-day martian surface conditions (in the absence of liquid water). A test for a meteoritic component of the soil can be conducted, as described below, by searching for the presence of Ni at the martian surface. The average abundance of nickel in an Fe-Ni meteorite is about 7% and, if present at measurable levels in the soil, would be indicative of an exogenic contribution. In addition, it may be possible to directly search for mineral phases common in meteorites. An understanding of the formation and evolution of the martian soil would not be complete without addressing the unusual reactivity discovered by the Viking Landers The presence of an inorganic oxidant, possibly one produced as a results of photochemical processes, is the most widely accepted explanation of the Viking results. Are these chemical species simply adsorbed on soil grains, or have they reacted with the metal oxide substrates and altered the mineral structures? Could a completely different (non-photochemical) process be responsible for the soil reactivity? The various ideas for the nature of this putative oxidant could be constrained by a measurement of the change in reactivity with depth. Different compositions will have different lifetimes and mobilities and thus will have different vertical profiles. Because the oxidizing compounds are believed to actively destroy organic molecules, determination of the reactivity gradient also has significant implications for the search for life on Mars. A DS2-based microprobe system can be instrumented for a 2003 micromission to investigate the origin and reactivity of the martian soil. These measurements would provide invaluable information regarding the climate history and exobiological potential of the planet. The NMR, X ray and chemiresistor measurement approach described embodies a highly synergistic and general set of soil interrogation methods for elements, compounds, and crystal structures and can also be applied to other geologic questions of interest. For example, if the capability for precise targeting of the probes is available, then in-situ investigations of suspected evaporite and hydrothermal deposits would be possible with the same set of instruments. Additional information is contained in the original.

Yen, Albert S.↗

Extended Testability Analysis Tool

The Extended Testability Analysis (ETA) Tool is a software application that supports fault management (FM) by performing testability analyses on the fault propagation model of a given system. Fault management includes the prevention of faults through robust design margins and quality assurance methods, or the mitigation of system failures. Fault management requires an understanding of the system design and operation, potential failure mechanisms within the system, and the propagation of those potential failures through the system. The purpose of the ETA Tool software is to process the testability analysis results from a commercial software program called TEAMS Designer in order to provide a detailed set of diagnostic assessment reports. The ETA Tool is a command-line process with several user-selectable report output options. The ETA Tool also extends the COTS testability analysis and enables variation studies with sensor sensitivity impacts on system diagnostics and component isolation using a single testability output. The ETA Tool can also provide extended analyses from a single set of testability output files. The following analysis reports are available to the user: (1) the Detectability Report provides a breakdown of how each tested failure mode was detected, (2) the Test Utilization Report identifies all the failure modes that each test detects, (3) the Failure Mode Isolation Report demonstrates the system s ability to discriminate between failure modes, (4) the Component Isolation Report demonstrates the system s ability to discriminate between failure modes relative to the components containing the failure modes, (5) the Sensor Sensor Sensitivity Analysis Report shows the diagnostic impact due to loss of sensor information, and (6) the Effect Mapping Report identifies failure modes that result in specified system-level effects.

Melcher, Kevin↗

Lead-Lag Control for Helicopter Vibration and Noise Reduction

As a helicopter transitions from hover to forward flight, the main rotor blades experience an asymmetry in flow field around the azimuth, with the blade section tangential velocities increasing on the advancing side and decreasing on the retreating side. To compensate for the reduced dynamic pressure on the retreating side, the blade pitch angles over this part of the rotor disk are increased. Eventually, a high enough forward speed is attained to produce compressibility effects on the advancing side of the rotor disk and stall on the retreating side. The onset of these two phenomena drastically increases the rotor vibratory loads and power requirements, thereby effectively establishing a limit on the maximum achievable forward speed. The alleviation of compressibility and stall (and the associated decrease in vibratory loads and power) would potentially result in an increased maximum forward speed. In the past, several methods have been examined and implemented to reduce the vibratory hub loads. Some of these methods are aimed specifically at alleviating vibration at very high flight speeds and increasing the maximum flight speed, while others focus on vibration reduction within the conventional flight envelope. Among the later are several types passive as well as active schemes. Passive schemes include a variety of vibration absorbers such as mechanical springs, pendulums, and bifilar absorbers. These mechanism are easy to design and maintain, but incur significant weight and drag penalties. Among the popular active control schemes in consideration are Higher Harmonic Control (HHC) and Individual Blade Control (IBC). HHC uses a conventional swash plate to generate a multi-cyclic pitch input to the blade. This requires actuators capable of sufficiently high power and bandwidth, increasing the cost and weight of the aircraft. IBC places actuators in the rotating reference frame, requiring the use of slip rings capable of transferring enough power to the actuators. Both schemes cause an increase in pitch link loads. Trailing Edge Flap (TEF) deployment can also used to generate unsteady aerodynamic forces and moments that counter the original vibratory loads, and thereby reduce rotor vibrations. While the vibrations absorbers, HHC, IBC, and TEF concepts discussed above attempt to reduce the vibratory loads, they do not specifically address the phenomena causing the vibrations at high advance ratios. One passive method that attempts to directly alleviate compressibility and stall, instead of reducing the ensuing vibrations, is the use of advanced tip designs. Taper, sweep, anhedral, and the manipulation of other geometric properties of the blade tips can reduce the severity of stall and compressibility effects , as well as reduce rotor power. A completely different approach to solve these problems is the tiltrotor configuration. As the forward velocity of the aircraft increases, the rotors, in this case, are tilted forward until they are perpendicular to the flow and act as propellers. This eliminates the edgewise flow encountered by conventional rotors and circumvents all the problems associated with flow asymmetry. However, the success involves a tremendous increase in cost and complexity of the aircraft. Another possible approach that has been proposed for the alleviation of vibratory loads at high forward flight speeds involves the use of controlled lead-lag motions to reduce the asymmetry in flow. A correctly phased 1/rev controlled lag motion could be introduced such that it produces a backward velocity on the advancing side and a forward velocity on the retreating side, to delay compressibility effects and stall to a higher advance ratio. Using a large enough lead-lag amplitude, the tip velocities could be reduced to levels encountered in hover. This concept was examined by two groups in the 1950's and early 1960's. In the United States, the Research Labs Division of United Aircraft developed a large lead-lag motion rotor, meant to achieve lag motion amplitudes up to 45 degrees. In order to reduce the required actuation force, the blade hinges were moved to 40% of the blade radius to increase the rotating lag frequency to approximately 1/rev. The blade hinges were redesigned to produce a flap-lag coupling so the large flapwise aerodynamic loads could be exploited to actuate the blades in the lag direction. A wind tunnel test of this rotor concept revealed actuation and blade motion scheduling problems. The project was eventually discontinued due to these problems and high blade stresses. Around the same time, at Boelkow in Germany, a similar lead-lag rotor program was conducted under the leadership of Hans Derschmidt. Here, too, the blade hinges were moved outboard to 34% radius to reduce the actuation loads. The main difference between this and the United Aircraft program was the use of a mechanical actuation scheme with maximum lead-lag motions of 400. This program was also discontinued for unclear reasons. The present study is directed toward conducting a comprehensive analytical examination to evaluate the effectiveness of controlled lead-lag motions in reducing vibratory hub loads and increasing maximum flight speed. Since both previous studies on this subject were purely experimental, only a limited data set and physical understanding of the problem was obtained. With the currently available analytical models and computational resources, the present effort is geared toward developing an in-depth physical understanding of the precise underlying mechanisms by which vibration reduction may be achieved. Additionally, in recognition of the fact that large amplitude lead-lag motions would - (i) be difficult to implement, and (ii) produce very large blade stresses; the present study examines the potential of only moderate-to-small lead-lag motions for reduction of vibratory hub loads. Using such an approach, the emphasis is not on eliminating the periodic variations in tangential velocity at the blade tip, but at best reducing these variations slightly so that compressibility and stall are delayed to slightly higher advance ratios. This study was conducted in two steps. In the first step, a hingeless helicopter rotor was modeled using rigid blades undergoing flap-lag-torsion rotations about spring restrained hinges and bearings. This model was then modified by separating the lead-lag degree of freedom into two components, a free and a prescribed motion. Using this model, a parametric study of the effect of phase and amplitude of a prescribed lead-lag motion on hub vibration was conducted. The data gathered was analyzed to obtain an understanding of the basic physics of the problem and show the capability of this method to reduce vibration and expand the flight envelope. In the second half of the study, the similar analysis was conducted using an elastic blade model to confirm the effects predicted by the simpler model.

Gandhi, Farhan↗

The optical and near-infrared colors of galaxies, 1: The photometric data

We present optical and near-infrared photometry and spectroscopic redshifts of a well defined sample of 171 field galaxies selected from three high galactic latitude fields. This data set forms the basis for subsequent studies to characterize the trends, dispersion, and evolution of rest-frame colors and image structure. A subset of 143 galaxies constitutes a magnitude-limited sample to B approx. 19.9-20.75 (depending on field), with a median redshift of 0.14, and a maximum redshift of 0.54. This subset is statistically representative in its sampling of the apparent color distribution of galaxies. Thirty six galaxies were selected to have the reddest red-optical colors in two redshift intervals between 0.2 less than z less than 0.3. Photometric passbands are similar to U, B, V, I, and K, and sample galaxy spectral energy distributions between 0.37 and 2.2 micrometers in the observed frame, or down to 0.26 micrometers in the rest frame for the most distant galaxies. B and K images of the entire sample are assembled to form the first optical and near-infrared atlas of a statistically-representative sample of field galaxies. We discuss techniques for faint field-galaxy photometry, including a working definition of a total magnitude, and a method for matching magnitudes in different passbands and different seeing conditions to ensure reliable, integrated colors. Photographic saturation, which substantially affects the brightest 12% of our sample in the optical bands, is corrected with a model employing measured plate-density distributions for each galaxy, calibrated via similar measurements for stars as a function of known saturation level. Both the relative and absolute calibration of our photometry are demonstrated.

Bershady, Matthew A.↗

Video event trigger and tracking system using fuzzy comparators

A video observation method and apparatus, the apparatus having a frame storage mechanism, a dividing mechanism, a plurality of fuzzy comparators and a trigger signal mechanism. The frame storage mechanism stores at least one non-current video frame of a viewing field. The dividing mechanism divides a current video frame of the viewing field and the at least one non-current video frame into a plurality of corresponding trigger sections. The plurality of fuzzy comparators each compare and detect a fuzzy logic difference between one trigger section of the current video frame and the corresponding trigger sections of the at least one non-current video frame, the number of fuzzy comparators being selected so that every trigger section of the current video frame is compared. The trigger signal mechanism provides a trigger signal when a fuzzy logic difference is detected between any of the corresponding current and non-current trigger sections. A video observation mechanism and data reducing mechanism may be included with the above apparatus or alone with only a frame storage mechanism, a single generic comparator and a trigger signal mechanism. The video observation mechanism provides a video data stream, wherein each pixel of each frame of a viewing field is provided as multiple bits of data. The data reducing mechanism reduces each set of multiple bits of data which correspond to each pixel to one bit of binary data based on whether the pixel has a level of grey which is above or below a threshold level of grey.

Williams, Glenn L.↗

Effect of Free Stream Turbulence on the Flow-Induced Background Noise of In-Flow Microphones

When making noise measurements of sound sources in flow using microphones immersed in an air stream or wind tunnel, the factor limiting the dynamic range of the measurement is, in many cases, the noise caused by the flow over the microphone. To lower this self-noise, and to protect the microphone diaphragm, an aerodynamic microphone forebody is usually mounted on the tip of the omnidirectional microphone. The microphone probe is then pointed into the wind stream. Even with a microphone forebody, however, the self-noise persists, prompting further research in the area of microphone forebody design for flow-induced self-noise reduction. The magnitude and frequency characteristics of in-flow microphone probe self-noise is dependent upon the exterior shape of the probe and on the level of turbulence in the onset flow, among other things. Several recent studies present new designs for microphone forebodies, some showing the forbodies' self-noise characteristics when used in a given facility. However, these self-noise characteristics may change when the probes are used in different facilities. The present paper will present results of an experimental investigation to determine an empirical relationship between flow turbulence and self-noise levels for several microphone forebody shapes as a function of frequency. As a result, the microphone probe self-noise for these probes will be known as a function of freestream turbulence, and knowing the freestream turbulence spectra for a given facility, the probe self-noise can be predicted. Flow-induced microphone self-noise is believed to be related to the freestream. turbulence by three separate mechanisms. The first mechanism is produced by large scale, as compared to the probe size, turbulence which appears to the probe as a variation in the angle of attack of the freestream. flow. This apparent angle of attack variation causes the pressure along the probe surface to fluctuate, and at the location of the sensor orifice this fluctuating surface pressure is sensed by the diaphragm as noise. The second mechanism is caused by the convection of smaller sized turbulence, on the order of the probe cross-section, which passes nearby or strikes the probe giving rise to a fluctuating pressure at the sensor orifice. And, the third mechanism is related to fine scale turbulence through its effects on boundary layer growth and transition to a turbulent boundary layer. The method for relating the probe self-noise to the freestream turbulence will be based on the method of K. J. Young5 from Boeing, who developed the technique and presented flow noise results for a Bruel & Kjaer Type 0385, 1/4 inch (6.35 mm) nose cone. The experimental set-up used in the present experiment is similar to that of Young and is described in the present paper. Finally, flow noise predictions are made using the empirical correlations. These predictions are then compared with actual flow noise measurements made in the National Full-Scale Aerodynamics Complex 40- by 80-Foot Wind Tunnel at NASA Ames Research Center.

Allen, Christopher S.↗

Diurnal Reflectance Changes in Vegetation Observed with AVIRIS

Among the most important short-term dynamic biological processes are diurnal changes in canopy water relations. Plant regulation of water transport through stomatal openings affects other gaseous transport processes, often dramatically decreasing photosynthetic fixation of carbon dioxide during periods of water stress. Water stress reduces stomatal conductance of water vapor through the leaf surface and alters the diurnal timing of stomatal opening. Under non-water stressed conditions, stomates typically open soon after dawn and transpire water vapor throughout the daylight period. During stress periods, stomates may close for part of the day, generally near mid-day. Under prolonged stress conditions, stomatal closure shifts to earlier times during the day; stomates may close by mid-morning and remain closed until the following morning - or remain closed entirely. Under these conditions the relationship between canopy greenness (e.g., measured with a vegetation index or by spectral mixture analysis) and photosynthetic fixation of carbon is lost and the remotely sensed vegetation metric is a poor predictor of gas exchange. Prediction of stomatal regulation and exchange of water and trace gases is critical for ecosystem and climate models to correctly estimate budgets of these gases and understand or predict other processes like gross and net ecosystem primary production. Plant gas exchange has been extensively studied by physiologists at the leaf and whole plant level and by biometeorologists at somewhat larger scales. While these energy driven processes follow a predictable if somewhat asymmetric diurnal cycle dependent on soil water availability and the constraints imposed by the solar energy budget, they are nonetheless difficult to measure at the tree and stand levels using conventional methods. Ecologists have long been interested in the potential of remote sensing for monitoring physiological changes using multi-temporal images. Much of this research has focused on day-to-day changes in water use, especially for agricultural applications. Ustin et al. showed seasonal changes in canopy water content in chaparral shrub could be estimated using optical methods. Vanderbilt et al. followed asymmetric diurnal changes in the reflectance of a walnut orchard, but could not attribute specific reflectance changes to specific changes in canopy architecture or physiology. Forests and shrub lands in California experience prolonged periods of drought, sometimes extending six months without precipitation. The conifer and evergreen chaparral communities common to the foothill region around the central valley of California retain their foliage throughout the summer and have low transpiration rates despite high net radiation and temperature conditions. In contrast, grasslands and drought resistant deciduous species in the same habitat are seasonally dormant in summer. Because of differences in the mechanisms of drought tolerance, rooting depth and physiology between different plant communities in the region, it is likely that they display differences in diurnal water relations. The presence of diverse plant communities provides an opportunity to investigate possible diurnal landscape patterns in water relations that could be observed by an airborne hyperspectral scanner. This investigation of AVIRIS data collected over forest and shrub land represents the continuation of a prior investigation involving spectral mixture analysis of diurnal effects in the same AVIRIS data set.

Vanderbilt, V. C.↗

Autonomous Ocean World Exploration: Advancement of a Virtual Testbed

The search for life (extinct or extant) and potentially habitable bodies in our solar system and beyond is one of the 12 priority science questions outlined in the National Acadamies’ 2022 decadal survey [5]. Extraterrestrial destinations containing liquid water present an opportunity to search for life as we know it, and in recent years an increasing number of such locations have been discovered within our solar system. Several Jovian moons—Europa, Ganymede, and Callisto [10]—and the Saturnian moons Enceladus [8] and Titan [9] are known or suspected to harbor massive subsurface oceans. Of these "ocean worlds", Europa is the focus of at least one planned NASA orbiter mission, Europa Clipper [4], and an early lander mission concept, the Europa Lander [2, 3]. Whereas most robotic missions to the Moon and Mars (e.g. orbiters, rovers, landers) to date have had ground controllers on Earth tightly involved in mission operations, missions to more distant worlds will require a high degree of onboard autonomy due to long communication lags and blackouts, harsh environments (radiation, cold), and more limited battery and hardware life. The past decade has seen great advances in both AI technologies and computing scalability and performance that offer promising solutions for spacecraft autonomy and motivate the software system and research programs described in this paper. The Ocean Worlds Autonomy Testbed for Exploration, Research, and Simulation (OceanWATERS) [1], which has been in development at the NASA Ames Research Center since 2018, is a virtual environment for testing lander autonomy solutions. It is built on the Robot Operating System (ROS), runs on consumer-grade Linux workstations, and was released as open source in 2020. OceanWATERS provides a physical and visual simulation of a prototypical lander in a Europa-like environment (Figure 1). The lander was modeled after requirements and specifications made in JPL’s Europa Lander Study of 2016 [3]. Simulated lander systems include stereo cameras and spotlights mounted on an antenna mast that pans and tilts, a 6 degrees of freedom (DoF) robotic arm with a force-torque sensor and two interchangeable end effectors, and a battery pack power system. The environment consists of multiple terrain models including a highly detailed model sourced from the FROST dataset [11], simulation of surrounding planetary bodies based on an ephemeris model, and lighting from the sun with associated surface illumination, reflectance, and shadows. Operations supported by OceanWATERS include panoramic and directed imaging of the environment and lander workspace, Cartesian and joint-level arm commanding, grinding of the terrain surface (e.g. digging a trench), and scooping of ground material (Figure 2) which can be discarded or collected as science samples in a receptacle that can be emptied (science operations themselves are not simulated). These operations are realized as ROS Actions and are complimented by a wide selection of telemetry that is continually produced by each lander subsystem. The power system model is driven by the open-source Generic Software Architecture for Prognostics (GSAP) [11] that predicts the battery’s remaining useful life and other characteristics. As a testbed for high-level autonomy, OceanWATERS provides an execution framework based on PLEXIL [12], an open-source plan specification language and execution engine developed largely at Ames. NASA's initial development of OceanWATERS, as well the Ocean Worlds Lander Autonomy Testbed (OWLAT) [6], a complimentary physical testbed developed at JPL, was the first step in a plan for realizing candidate onboard autonomy solutions for such planetary landers. In 2020 NASA solicited applications for its Autonomous Robotics Research for Ocean Worlds (ARROW) program, and in 2021 the similar Concepts for Ocean worlds Life Detection Technology (COLDTech) program. Collectively six research teams, based in universities and companies across the United States, were awarded grants to develop and demonstrate autonomy solutions on OceanWATERS and OWLAT. These 1–2-year projects have now finished or are nearing completion, and a wide variety of autonomy challenges in ocean world surface missions were addressed. Prototyped and demonstrated solutions have included autonomous discovery, response and adaptation to system faults and unexpected environmental events, world model synthesis through perception, plan synthesis using learned models, methods to optimize sample target selection and prioritize science data transmission, extension of PLEXIL for stochastic decision-making, and an integration of a model of JPL’s mission-ready COLDArm [7]. Technologies used in these projects include many forms of machine learning, causal reasoning, automated planning, Markov decision processes, formal methods, and other advanced techniques. A more detailed summary of the ARROW and COLDTech projects is given herein. OceanWATERS has had significant enhancements since its open-source release in 2020. Many of its new features were driven or shaped by feedback from the ARROW and COLDTech teams and requirements of their projects. In support of enabling autonomous adaptation to spacecraft faults (a specific capability solicited by both programs), a fault injection and detection framework was developed that supports a wide and growing range of fault types such as locked joints, image loss, and battery failures. The power system model was completed and integrated into the simulator, starting as a single-cell battery model and later upgraded to a multi-cell model with associated faults such as cell disconnection. Arm/terrain interaction was improved by adding a force-torque sensor and associated faults, and an analytic dig force model based on the Balovnev bucket force equations. Environment fidelity was increased by modeling terrain deformation resulting from digging and scooping; visual improvements were made in textures, lighting, and shadows. To facilitate interoperation with OWLAT, a unified command and telemetry interface between the testbeds was developed at the ROS level, along with a PLEXIL interface. The number of lander operations was greatly expanded (e.g. with Cartesian-based arm and antenna movement), and a framework was designed for users to build their own lander actions. A GUI for PLEXIL plan selection was created (Figure 3), and an expansive set of plans were added, such as those that illustrate patterns for fault handling. This paper provides a self-contained high-level description of OceanWATERS, focusing on more detailed coverage of the aforementioned enhancements. It provides a high-level summary of the projects undertaken by participants in the ARROW and COLDTech programs and how these efforts have helped shape OceanWATERS. Finally, potential future work and directions for the testbed are listed, as likely informed by the recent planetary science decadal survey [5].

K Michael Dalal↗

A quick look method of detecting water of hydration in small solar system bodies

The action of aqueous alteration of the near-subsurface material on asteroids and probably some satellites in the Solar System constitutes part of the formation history of the Solar System. The C-class asteroids (and subclasses B, G, and F) were initially believed to have undergone aqueous alteration based on their low albedos and neutral broadband visible and near-infrared colors. These spectra exhibit a sharp drop at wavelengths shorter than 0.55 microns due to a strong ferric oxide intervalence change transfer transition. This IVCT comprises multiple absorptions that are not uniquely indicative of phyllosilicates, but rather are present in the spectrum of any object containing Fe(2+) and Fe(3+) in its surface material. A definitive indication of aqueous alteration came when the broad IR absorption feature having a minimum near 3.0 microns indicative of structural hydroxyl (OH) and interlayer and and absorbed water in phyllosilicates was identified in the IR photometry of many C-class asteroids. Additional mineralogical compositional evidence of aqueously-altered asteroids has come as the result of high SNR narrowband spectrophotometry in the visible and near-infrared spectral regions taken using a CCD/spectrograph combination. An absorption feature centered at 0.7 microns indicative of an Fe(2+) - Fe(3+) charge transfer transition in oxidized iron in phyllosilicates in spectra of some low-albedo asteroids, especially C and G class, and CM2 carbonaceous chondrite meteorites was identified in the CCD spectra. The correlation and covariance of the 0.7-micron and 3.0-micron features were examined by comparing observations of asteroids common to both the CCD reflectance spectra and the 3.0-micron multicolor photometry data sets. Thirty-one pairs of observations were included in this training group. The results of the statistical study indicate that with a 95% confidence level, 84% of the objects observed either having or not having the 0.7-micron feature will correspondingly have the 3.0-micron water of hydration feature.

Vilas, F.↗

Characterization of zinc selenide single crystals

ZnSe single crystals of high quality and low impurity levels are desired for use as substrates in optoelectronic devices. This is especially true when the device requires the formation of homoepitaxial layers. While ZnSe is commercially available, it is at present extremely expensive due to the difficulty of growing single crystal boules with low impurity content and the resultant low yields. Many researchers have found it necessary to heat treat the crystals in liquid Zn in order to remove the impurities, lower the resistivity and activate the photoluminescence at room temperature. The physical vapor transport method (PVT) has been successfully used at MSFC to grow many single crystals of II-VI semiconducting materials including ZnSe. The main goal at NASA has been to try to establish the effect of gravity on the growth parameters. To this effect, crystals have been grown vertically upwards or horizontally. Both (111) and (110) oriented ZnSe crystals have been obtained via unseeded PVT growth. Preliminary characterization of the horizontally grown crystals has revealed that Cu is a major impurity and that the low temperature photoluminescence spectra is dominated by the copper peak. The ratio of the copper peak to the free exciton peak is being used to determine variations in composition throughout the crystal. It was the intent of this project to map the copper composition of various crystals via photoluminescence first, then measure their electrical resistivity and capacitance as a function of frequency before proceeding with a heat treatment designed to remove the copper impurities. However, equipment difficulties with the photoluminescence set up, having to establish a procedure for measuring the electrical properties of the as-grown crystals and time limitations made us re-evaluate the project goals. Vertically grown samples designated as ZnSe-25 were chosen to be measured electrically since they were not expected to show as much variation in their composition through their cross-section as the horizontally grown samples.

Gerhardt, Rosario A.↗

A Scenario-Based Approach to Assess Continuity Gaps in Earth Observations

Decision analysis processes outlined in systems engineering references, such as the Systems Engineering Body of Knowledge [1] or the NASA Systems Engineering Handbook [2], recommend following a series of steps to support decision-making for engineering applications. These steps are consistent across the literature and typically involve defining objectives, defining relevant criteria against which candidate alternatives can be assessed, selecting an evaluation method, assessing alternatives, and making a recommendation based on this assessment. These approaches are well-suited to assess systems for which sets of common objectives and constraints can be identified. A vast body of literature describes their application to numerous engineering problems, and many methods have been developed since the middle of the 20th century to support the assessment of candidate alternatives for such systems [3, 4, 5, 6]. These approaches, however, show some limitations when decisions pertain to systems of systems. These collections of individual systems typically do not share a common set of objectives and constraints, tend to be highly complex, and the definition of assessable, high-value candidate alternatives poses a challenge due to interdependencies. The formulation, development, operation, and funding of Earth observing missions are such that gaps may occur between missions, therefore impacting the continuity of measurements. These Earth observing mission architectures constitute systems of systems for which a common set of objectives and constraints can be challenging to define. When tasked with assessing gaps that may occur for a series of spaceborne missions that address similar science parameters, the authors therefore proposed to depart from the traditional objective- and criteria-based approach and instead adopted a scenario-based approach. This approach does not require a prioritized set of common objectives; rather, it assesses the impact of possible decisions on the system-of-systems and characterizes its possible future states. It also offers the flexibility to adjust assumptions, update inputs, and refine supporting models over time, while enabling the rapid, early identification of challenges and key decision points. By simulating the impact of potential decisions on the entirety of the system-of-systems, the approach enables the synchronization of multiple decisions that typically occur at the level of the individual system. This paper provides an overview of the approach that was developed to assess continuity gaps for a series of operational and planned spaceborne Earth observing missions, a description of its application to the problem at hand, and a discussion of currently known limitations.

scenario earth observations continuity gaps↗

A Scenario-Based Approach to Assess Continuity Gaps in Earth Observations

Decision analysis processes outlined in systems engineering references, such as the Systems Engineering Body of Knowledge [1] or the NASA Systems Engineering Handbook [2], recommend following a series of steps to support decision-making for engineering applications. These steps are consistent across the literature and typically involve defining objectives, defining relevant criteria against which candidate alternatives can be assessed, selecting an evaluation method, assessing alternatives, and making a recommendation based on this assessment. These approaches are well-suited to assess systems for which sets of common objectives and constraints can be identified. A vast body of literature describes their application to numerous engineering problems, and many methods have been developed since the middle of the 20th century to support the assessment of candidate alternatives for such systems [3, 4, 5, 6]. These approaches, however, show some limitations when decisions pertain to systems of systems. These collections of individual systems typically do not share a common set of objectives and constraints, tend to be highly complex, and the definition of assessable, high-value candidate alternatives poses a challenge due to interdependencies. The formulation, development, operation, and funding of Earth observing missions are such that gaps may occur between missions, therefore impacting the continuity of measurements. These Earth observing mission architectures constitute systems of systems for which a common set of objectives and constraints can be challenging to define. When tasked with assessing gaps that may occur for a series of spaceborne missions that address similar science parameters, the authors therefore proposed to depart from the traditional objective- and criteria-based approach and instead adopted a scenario-based approach. This approach does not require a prioritized set of common objectives; rather, it assesses the impact of possible decisions on the system-of-systems and characterizes its possible future states. It also offers the flexibility to adjust assumptions, update inputs, and refine supporting models over time, while enabling the rapid, early identification of challenges and key decision points. By simulating the impact of potential decisions on the entirety of the system-of-systems, the approach enables the synchronization of multiple decisions that typically occur at the level of the individual system. This paper provides an overview of the approach that was developed to assess continuity gaps for a series of operational and planned spaceborne Earth observing missions, a description of its application to the problem at hand, and a discussion of currently known limitations.

scenario earth observations continuity gaps↗

Estimates of Ground Temperature and Atmospheric Moisture from CERES Observations

A method is developed to retrieve surface ground temperature (Tg) and atmospheric moisture using clear sky fluxes (CSF) from CERES-TRMM observations. In general, the clear sky outgoing long-wave radiation (CLR) is sensitive to upper level moisture (q(sub h)) over wet regions and Tg over dry regions The clear sky window flux from 800 to 1200 /cm (RadWn) is sensitive to low level moisture (q(sub j)) and Tg. Combining these two measurements (CLR and RadWn), Tg and q(sub h) can be estimated over land, while q(sub h) and q(sub t) can be estimated over the oceans. The approach capitalizes on the availability of satellite estimates of CLR and RadWn and other auxiliary satellite data. The basic methodology employs off-line forward radiative transfer calculations to generate synthetic CSF data from two different global 4-dimensional data assimilation products. Simple linear regression is used to relate discrepancies in CSF to discrepancies in Tg, q(sub h) and q(sub t). The slopes of the regression lines define sensitivity parameters that can be exploited to help interpret mismatches between satellite observations and model-based estimates of CSF. For illustration, we analyze the discrepancies in the CSF between an early implementation of the Goddard Earth Observing System Data Assimilation System (GEOS-DAS) and a recent operational version of the European Center for Medium-Range Weather Prediction data assimilation system. In particular, our analysis of synthetic total and window region SCF differences (computed from two different assimilated data sets) shows that simple linear regression employing (Delta)Tg and broad layer (Delta)q(sub l) from 500 hPa to surface and (Delta)q(sub h) from 200 to 500 hPa provides a good approximation to the full radiative transfer calculations, typically explaining more than 90% of the 6-hourly variance in the flux differences. These simple regression relations can be inverted to "retrieve" the errors in the geophysical parameters. Uncertainties (normalized by standard deviation) in the monthly mean retrieved parameters range from 7% for (Delta)T to about 20% for (Delta)q(sub t). Our initial application of the methodology employed an early CERES-TRMM data set (CLR and Radwn) to assess the quality of the GEOS2 data. The results showed that over the tropical and subtropical oceans GEOS2 is, in general, too wet in the upper troposphere (mean bias of 0.99 mm) and too dry in the lower troposphere (mean bias of -4.7 mm). We note that these errors, as well as a cold bias in the Tg, have largely been corrected in the current version of GEOS-2 with the introduction of a land surface model, a moist turbulence scheme and the assimilation of SSTM/I total precipitable water.

Wu, Man Li C.↗

Virtual Flight Demonstration of the Stratospheric Dual-Aircraft Platform

A baseline configuration for the dual-aircraft platform (DAP) concept is described and evaluated in a physics-based flight dynamics simulations for two month-long missions as a communications relay in the lower stratosphere above central Florida, within 150-miles of downtown Orlando.The DAP configuration features two large glider-like (130 ft wing span) unmanned aerial vehicles connected via a long adjustable cable (total extendible length of 3000 ft) which effectively sail without propulsion using available wind shear. Use of onboard LiDAR wind profilers to forecast wind distributions are found to be necessary to enable the platform to efficiently adjust flight conditions to remain sailing by finding sufficient wind shear across the platform. The aircraft derive power from solar cells, like a conventional solar aircraft, but also extract wind power using the propeller as a turbine when there is an excess of wind shear available.Month-long atmospheric profiles (at 3-5 min intervals) in the vicinity of 60,000-ft are derived from archived data measured by the 50-Mhz Doppler Radar Wind Profiler at Cape Canaveral and used in the DAP flight simulations. A cursory evaluation of these datasets show that sufficient wind shear for DAP sailing is persistent, suggesting that DAP could potentially sail over 90% of the month-long durations even when limited by modest ascent/descent rates.DAP's novel guidance software uses a non-linear constrained optimization technique to define waypoints such that sailing mode of flight is maintained where possible, and minimal thrust is required where sailing is not practical. A set of constraints are identified which result in waypoints that enable efficient flight (i.e., minimal use of propulsion) over the two month-long flight simulations. Waypoint solutions may need to be tabulated for a wide range of potential atmospheric conditions and stored onboard for quick retrieval on a real DAP.DAP's flight control software uses an unconventional mixture of spacecraft and aircraft control techniques. Flight simulations confirms that this controls approach enables the platform to consistently reach successive waypoints over the month-long flight simulations.The ability of DAP to transition between the sailing mode (i.e., cable tension is high) and standard formation flight (i.e., cable tension is low) is a vital capability (e.g., to enable intermittent turns while stationkeeping). A new method to perform these transitions has been identified and characterized with flight simulation which requires special aircraft modifications.The energy-usage of the DAP configuration during two month-long stationkeeping missions over central Florida (i.e., stationkeeping over Orlando) is evaluated and compared to that of a pure solar aircraft of the same weight and aerodynamic performance. DAP is shown to consistently reduce net propulsion usage while simultaneously increasing solar energy capture.A baseline 700 GHz communications system is described and its performance evaluated for the proposed mission over central Florida. It is found that the variable roll orientation of the aircraft would increase the power required to maintain coverage over the stationkeeping radius of 150 miles (e.g., by as much as 100% when DAP is 150 miles from Orlando), compared to level flight. This effect can be mitigated via additional antenna design complexity or a more restricted stationkeeping radius.

Demonstrations↗