Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “range finding algorithm”

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 163 records · Page 9

Validation of Ozone Profile Retrievals Derived from the OMPS LP Version 2.5 Algorithm Against Correlative Satellite Measurements

The Limb Profiler (LP) is a part of the Ozone Mapping and Profiler Suite launched on board of the Suomi NPP satellite in October 2011. The LP measures solar radiation scattered from the atmospheric limb in ultraviolet and visible spectral ranges between the surface and 80 km. These measurements of scattered solar radiances allow for the retrieval of ozone profiles from cloud tops up to 55 km. The LP started operational observations in April 2012. In this study we evaluate more than 5.5 years of ozone profile measurements from the OMPS LP processed with the new NASA GSFC version 2.5 retrieval algorithm. We provide a brief description of the key changes that had been implemented in this new algorithm, including a pointing correction, new cloud height detection, explicit aerosol correction and a reduction of the number of wavelengths used in the retrievals. The OMPS LP ozone retrievals have been compared with independent satellite profile measurements obtained from the Aura Microwave Limb Sounder (MLS), Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) and Odin Optical Spectrograph and InfraRed Imaging System (OSIRIS). We document observed biases and seasonal differences and evaluate the stability of the version 2.5 ozone record over 5.5 years. Our analysis indicates that the mean differences between LP and correlative measurements are well within required +/- 10% between 18 and 42 km. In the upper stratosphere and lower mesosphere (>43 km) LP tends to have a negative bias. We find larger biases in the lower stratosphere and upper troposphere, but LP ozone retrievals have significantly improved in version 2.5 compared to version 2 due to the implemented aerosol correction. In the northern high latitudes we observe larger biases between 20 and 32 km due to the remaining thermal sensitivity issue. Our analysis shows that LP ozone retrievals agree well with the correlative satellite observations in characterizing vertical, spatial and temporal ozone distribution associated with natural processes, like the seasonal cycle and quasi-biennial oscillations. We found a small positive drift approx. 0.5%/yr in the LP ozone record against MLS and OSIRIS that is more pronounced at altitudes above 35 km. This pattern in the relative drift is consistent with a possible 100m drift in the LP sensor pointing detected by one of our altitude-resolving methods.

Ozone Mapping↗

Contact Multigraph Routing: Overview and Implementation

In Delay Tolerant Networking (DTN), the standard routing algorithm used to navigate time-varying networks has been Contact Graph Routing (CGR). In CGR, a globally distributed list of contacts, periods during which two DTN nodes may communicate, is used to construct a contact graph, in which contacts are vertices. A version of Dijkstra’s algorithm can then be used to find paths through this model of the timevarying network. However, since contact graphs may be large compared to the network, potentially growing with the square of the number of network nodes and linearly with the time interval represented, the resulting algorithm does not scale well with the size of the network or time. Any improvement to the routing algorithm will bring significant returns to scale. In a previous paper, we briefly introduced an alternative to the contact graph model for routing. This alternative model is based on a multigraph (a graph in which there may be multiple edges between a pair of vertices) where vertices represent network nodes instead of contacts. A version of Dijkstra’s algorithm in these multigraph models reduces the time needed to perform the same routing computations done in the existing CGR algorithm. Moreover, a modified version of Yen’s algorithm for multigraphs is included. Our variation of CGR, which we call Contact Multigraph Routing (CMR), provides an in-line replacement for the previously used pathfinding algorithms. This paper describes an implementation created based on the CMR approach, and experimental comparisons to traditional CGR are given. In addition, we explore some additional modifications to the routing pipeline traditionally assumed in CGR. These modifications range from the theoretical to the practical in terms of size and scope. We step forward our understanding of sheaftheoretic networking and describe how to model the routing pipeline using sheaves. We detail some enhanced route selection criteria that addresses some of the added complexity of DTNbased systems. We also include a future works section on future improvements and implementations that would be of service to the broader DTN community.

contact graph routing↗

Heuristics for Relevancy Ranking of Earth Dataset Search Results

As the Variety of Earth science datasets increases, science researchers find it more challenging to discover and select the datasets that best fit their needs. The most common way of search providers to address this problem is to rank the datasets returned for a query by their likely relevance to the user. Large web page search engines typically use text matching supplemented with reverse link counts, semantic annotations and user intent modeling. However, this produces uneven results when applied to dataset metadata records simply externalized as a web page. Fortunately, data and search provides have decades of experience in serving data user communities, allowing them to form heuristics that leverage the structure in the metadata together with knowledge about the user community. Some of these heuristics include specific ways of matching the user input to the essential measurements in the dataset and determining overlaps of time range and spatial areas. Heuristics based on the novelty of the datasets can prioritize later, better versions of data over similar predecessors. And knowledge of how different user types and communities use data can be brought to bear in cases where characteristics of the user (discipline, expertise) or their intent (applications, research) can be divined. The Earth Observing System Data and Information System has begun implementing some of these heuristics in the relevancy algorithm of its Common Metadata Repository search engine.

science data management↗

Relevancy Ranking of Satellite Dataset Search Results

As the Variety of Earth science datasets increases, science researchers find it more challenging to discover and select the datasets that best fit their needs. The most common way of search providers to address this problem is to rank the datasets returned for a query by their likely relevance to the user. Large web page search engines typically use text matching supplemented with reverse link counts, semantic annotations and user intent modeling. However, this produces uneven results when applied to dataset metadata records simply externalized as a web page. Fortunately, data and search provides have decades of experience in serving data user communities, allowing them to form heuristics that leverage the structure in the metadata together with knowledge about the user community. Some of these heuristics include specific ways of matching the user input to the essential measurements in the dataset and determining overlaps of time range and spatial areas. Heuristics based on the novelty of the datasets can prioritize later, better versions of data over similar predecessors. And knowledge of how different user types and communities use data can be brought to bear in cases where characteristics of the user (discipline, expertise) or their intent (applications, research) can be divined. The Earth Observing System Data and Information System has begun implementing some of these heuristics in the relevancy algorithm of its Common Metadata Repository search engine.

science data management↗

Validation of Long-Term Global Aerosol Climatology Project Optical Thickness Retrievals Using AERONET and MODIS Data

A comprehensive set of monthly mean aerosol optical thickness (AOT) data from coastal and island AErosol RObotic NETwork (AERONET) stations is used to evaluate Global Aerosol Climatology Project (GACP) retrievals for the period 1995-2009 during which contemporaneous GACP and AERONET data were available. To put the GACP performance in broader perspective, we also compare AERONET and MODerate resolution Imaging Spectroradiometer (MODIS) Aqua level-2 data for 2003-2009 using the same methodology. We find that a large mismatch in geographic coverage exists between the satellite and ground-based datasets, with very limited AERONET coverage of open-ocean areas. This is especially true of GACP because of the smaller number of AERONET stations at the early stages of the network development. Monthly mean AOTs from the two over-the-ocean satellite datasets are well-correlated with the ground-based values, the correlation coefficients being 0.81-0.85 for GACP and 0.74-0.79 for MODIS. Regression analyses demonstrate that the GACP mean AOTs are approximately 17%-27% lower than the AERONET values on average, while the MODIS mean AOTs are 5%-25% higher. The regression coefficients are highly dependent on the weighting assumptions (e.g., on the measure of aerosol variability) as well as on the set of AERONET stations used for comparison. Comparison of over-the-land and over-the-ocean MODIS monthly mean AOTs in the vicinity of coastal AERONET stations reveals a significant bias. This may indicate that aerosol amounts in coastal locations can differ significantly from those in adjacent open-ocean areas. Furthermore, the color of coastal waters and peculiarities of coastline meteorological conditions may introduce biases in the GACP AOT retrievals. We conclude that the GACP and MODIS over-the-ocean retrieval algorithms show similar ranges of discrepancy when compared to available coastal and island AERONET stations. The factors mentioned above may limit the performance of the validation procedure and cause us to caution against a direct extrapolation of the presented validation results to the entirety of the GACP dataset.

MODIS (radiometry)↗

EHW Approach to Temperature Compensation of Electronics

Efforts are under way to apply the concept of evolvable hardware (EHW) to compensate for variations, with temperature, in the operational characteristics of electronic circuits. To maintain the required functionality of a given circuit at a temperature above or below the nominal operating temperature for which the circuit was originally designed, a new circuit would be evolved; moreover, to obtain the required functionality over a very wide temperature range, there would be evolved a number of circuits, each of which would satisfy the performance requirements over a small part of the total temperature range. The basic concepts and some specific implementations of EHW were described in a number of previous NASA Tech Briefs articles, namely, "Reconfigurable Arrays of Transistors for Evolvable Hardware" (NPO-20078), Vol. 25, No. 2 (February 2001), page 36; Evolutionary Automated Synthesis of Electronic Circuits (NPO- 20535), Vol. 26, No. 7 (July 2002), page 37; "Designing Reconfigurable Antennas Through Hardware Evolution" (NPO-20666), Vol. 26, No. 7 (July 2002), page 38; "Morphing in Evolutionary Synthesis of Electronic Circuits" (NPO-20837), Vol. 26, No. 8 (August 2002), page 31; "Mixtrinsic Evolutionary Synthesis of Electronic Circuits" (NPO-20773) Vol. 26, No. 8 (August 2002), page 32; and "Synthesis of Fuzzy-Logic Circuits in Evolvable Hardware" (NPO-21095) Vol. 26, No. 11 (November 2002), page 38. To recapitulate from the cited prior articles: EHW is characterized as evolutionary in a quasi-genetic sense. The essence of EHW is to construct and test a sequence of populations of circuits that function as incrementally better solutions of a given design problem through the selective, repetitive connection and/or disconnection of capacitors, transistors, amplifiers, inverters, and/or other circuit building blocks. The connection and disconnection can be effected by use of field-programmable transistor arrays (FPTAs). The evolution is guided by a search-andoptimization algorithm (in particular, a genetic algorithm) that operates in the space of possible circuits to find a circuit that exhibits an acceptably close approximation of the desired functionality. The evolved circuits can be tested by mathematical modeling (that is, computational simulation) only, tested in real hardware, or tested in combinations of computational simulation and real hardware.

Stoica, Adrian↗

VLSI neuroprocessors

Electronic and optoelectronic hardware implementations of highly parallel computing architectures address several ill-defined and/or computation-intensive problems not easily solved by conventional computing techniques. The concurrent processing architectures developed are derived from a variety of advanced computing paradigms including neural network models, fuzzy logic, and cellular automata. Hardware implementation technologies range from state-of-the-art digital/analog custom-VLSI to advanced optoelectronic devices such as computer-generated holograms and e-beam fabricated Dammann gratings. JPL's concurrent processing devices group has developed a broad technology base in hardware implementable parallel algorithms, low-power and high-speed VLSI designs and building block VLSI chips, leading to application-specific high-performance embeddable processors. Application areas include high throughput map-data classification using feedforward neural networks, terrain based tactical movement planner using cellular automata, resource optimization (weapon-target assignment) using a multidimensional feedback network with lateral inhibition, and classification of rocks using an inner-product scheme on thematic mapper data. In addition to addressing specific functional needs of DOD and NASA, the JPL-developed concurrent processing device technology is also being customized for a variety of commercial applications (in collaboration with industrial partners), and is being transferred to U.S. industries. This viewgraph p resentation focuses on two application-specific processors which solve the computation intensive tasks of resource allocation (weapon-target assignment) and terrain based tactical movement planning using two extremely different topologies. Resource allocation is implemented as an asynchronous analog competitive assignment architecture inspired by the Hopfield network. Hardware realization leads to a two to four order of magnitude speed-up over conventional techniques and enables multiple assignments, (many to many), not achievable with standard statistical approaches. Tactical movement planning (finding the best path from A to B) is accomplished with a digital two-dimensional concurrent processor array. By exploiting the natural parallel decomposition of the problem in silicon, a four order of magnitude speed-up over optimized software approaches has been demonstrated.

Kemeny, Sabrina E.↗

Parameterizing Spectral Surface Reflectance Relationships for the Dark Target Aerosol Algorithm Applied to A Geostationary Imager

Originally developed for the Moderate Resolution Imaging Spectroradiometer (MODIS) in polar, sun-synchronous low-earth orbit (LEO), the Dark Target (DT) aerosol retrieval algorithm relies on the assumption of a Surface Reflectance Parameterization (SRP) over land surfaces. Specifically for vegetated and dark-soiled surfaces, values of surface reflectance in blue and red visible-wavelength bands are assumed to be nearly linearly related to each other and to the value in a shortwave infrared (SWIR) wavelength band. This SRP also includes dependencies on scattering angle and a normalized difference vegetation index computed from two SWIR bands (NDVISWIR). As the DT retrieval algorithm is being ported to new sensors to continue and expand the aerosol data record, we assess whether the MODIS-assumed SRP can be used for these sensors. Here, we specifically assess SRP for the Advanced Baseline Imager (ABI) aboard, the Geostationary Operational Environmental Satellite (GOES)-16/East (ABIE). First, we find that using MODIS-based SRP leads to higher biases and artificial diurnal signatures in aerosol optical depth (AOD)retrievals from ABIE. The primary reason appears to be that geostationary orbit (GEO) encounters an entirely different set of observation geometry than does LEO, primarily with regards to solar angles coupled with fixed view angles. Therefore, we have developed a new SRP for GEO that draws the angular shape of the surface bidirectional reflectance. We also introduce modifications to the parametrization of both red-SWIR and blue-red spectral relationships to include additional information.The revised Red-SWIR SRP includes solar zenith angle, NDVISWIR, and land-type percentage from an ancillary database. The blue-red SRP adds dependencies on the scattering angle and NDVISWIR. The new SRPs improve the AOD retrieval of ABIE in terms of overall less bias and mitigation of the overestimation around local noon. The average bias of DT AOD compared to AERONET AOD shows a reduction from 0.08 to 0.03, while the bias of local solar noon decreases from 0.12 to 0.03.The agreement between DT and AERONET AOD is established through regression slope of 1.06 and y-intercept of 0.01with correlation coefficient is0.74. By using the new SRP, the percentage of data falling within the expected error range (±0.05 + 15%) is notably risen from 54% to 78%.

Mijin Kim↗

Launch Vehicle Propulsion Parameter Design Multiple Selection Criteria

The optimization tool described herein addresses and emphasizes the use of computer tools to model a system and focuses on a concept development approach for a liquid hydrogen/liquid oxygen single-stage-to-orbit system, but more particularly the development of the optimized system using new techniques. This methodology uses new and innovative tools to run Monte Carlo simulations, genetic algorithm solvers, and statistical models in order to optimize a design concept. The concept launch vehicle and propulsion system were modeled and optimized to determine the best design for weight and cost by varying design and technology parameters. Uncertainty levels were applied using Monte Carlo Simulations and the model output was compared to the National Aeronautics and Space Administration Space Shuttle Main Engine. Several key conclusions are summarized here for the model results. First, the Gross Liftoff Weight and Dry Weight were 67% higher for the design case for minimization of Design, Development, Test and Evaluation cost when compared to the weights determined by the minimization of Gross Liftoff Weight case. In turn, the Design, Development, Test and Evaluation cost was 53% higher for optimized Gross Liftoff Weight case when compared to the cost determined by case for minimization of Design, Development, Test and Evaluation cost. Therefore, a 53% increase in Design, Development, Test and Evaluation cost results in a 67% reduction in Gross Liftoff Weight. Secondly, the tool outputs define the sensitivity of propulsion parameters, technology and cost factors and how these parameters differ when cost and weight are optimized separately. A key finding was that for a Space Shuttle Main Engine thrust level the oxidizer/fuel ratio of 6.6 resulted in the lowest Gross Liftoff Weight rather than at 5.2 for the maximum specific impulse, demonstrating the relationships between specific impulse, engine weight, tank volume and tank weight. Lastly, the optimum chamber pressure for Gross Liftoff Weight minimization was 2713 pounds per square inch as compared to 3162 for the Design, Development, Test and Evaluation cost optimization case. This chamber pressure range is close to 3000 pounds per square inch for the Space Shuttle Main Engine.

Shelton, Joey Dewayne↗

Investigation of Magmatic Activities on Early Mars Using Igneous Mineral Chemistry in Gale Crater, Mars

One objective of rover missions is exploring the geological context of the surroundings. Over the years, igneous petrology and sedimentology have been disconnected, the first investigating magmatic processes and volcanic activities, and the second seeking environmental conditions in the past and assessing the habitability of the planet. Although different, one is related to the other: igneous rocks are altered and broken down, leading to the formation of sedimentary rocks, which can in turn be used to back out the nature of their magmatic source. The Curiosity rover that landed in the 3.7 Gyr old impact crater Gale is traveling through sedimentary rocks. About fifty float rocks have been observed, and several of them with ambiguous texture and composition have been classified as igneous or sedimentary depending on studies such as Jake_M. The composition of several unambiguous igneous rocks has been analyzed [4- 6] but their heterogeneity at a larger instrumental (measurement size < 2 cm) scale prevents the measurement of a bulk composition as performed on Earth. An original approach avoiding these two last issues is to consider igneous mineral chemistry analyzed within igneous and sedimentary rocks to assess magmatic processes that could have formed them. Most Curiosity data are used to explore ancient environmental conditions, and a significant number of compositional analyses are under-explored for constraining magmatic activities. We will present how we can make use of sedimentary data for investigating igneous processes in the vicinity of Gale crater. Geological Context: We focus on the first 750 martian days, corresponding to measurements in a coherent lacustrine sedimentary unit called Bradbury, because all sedimentary rocks were sourced from the same watershed and appear to have a consistent source with minimal alteration [2-3]. Igneous detrital minerals including feldspar and pyroxene, are observed in sedimentary rocks. Monte Carlo models showed that minimal cation loss is observed based on the composition of all Bradbury rocks, implying negligible weathering [3]. Although clay minerals are detected in few rocks [7], chemical compositions of rocks can be explained by a mixture of primary igneous minerals [3]. Variation of composition within Bradbury rocks can be explained by mineral sorting and one distinct source component. While a common magmatic source is suggested, Bradbury sediments likely come from several volcanic eruptions from a single magmatic chamber [9- 10]. The occurrence of alkali minerals like sanidine and K-rich rocks throughout Bradbury supports the presence of a potassic component, likely trachytic, while plagioclase and a mafic composition suggest a basaltic component [8-9]. Instruments: Mineral chemistry can be estimated by three instruments onboard Curiosity. The CheMin instrument enables detection of mineral assemblages using X-ray diffraction (XRD). Using Rietveld refinement, each mineral is identified according to their 1D XRD pattern [11]. Note that distinction between pyroxene minerals is challenging with the CheMin instrument due to overlapping peaks on XRD patterns and low angular resolution of the instrument [12]. Then, using least square regression and optimization algorithms based on unit-cell parameters, mineral chemistry has been estimated by [11]. Plagioclase composition has been estimated using the NaAlSi3O8- CaAl2Si2O8 system and alkali feldspar is based on the NaAlSi3O8-KAlSi3O8 system (stars in Fig.1). Two mudstone samples (John Klein and Cumberland) and one sandstone sample (Windjana) were analyzed by CheMin at Yellowknife Bay and Kimberley, respectively. Figure 1. Ternary diagrams of feldspars (top) and pyroxene (bottom) quadrilateral. Stars correspond to CheMin composition and the gray patches to ChemCam composition. The colored dots are the composition of feldspar and pyroxene that crystallized during fractional crystallization at FMQ+1 of a melt extracted at distinct melting degree during the adiabatic ascent of a primitive mantle composition, without any water (left panels) and with 0.5 wt.% of water (right panels) at distinct pressure. The ChemCam instrument enables the analysis of the chemical compositions of rocks at hundreds of micrometer scale (350-550 μm) using laser induced breakdown spectroscopy (LIBS), which may provide the composition of minerals when they are larger than the beam spot (>550 μm) [13]. Within >5000 LIBS points, we performed a typical stoichiometric filtering allowing us to distinguish 56 feldspar and 10 pyroxene mineral compositions (grey patches in Fig. 1). Finally, the Alpha Particle X-ray Spectrometer (APXS) analyzes the composition of rocks with a 1.6 cm diameter spot size. Monte Carlo mass balance modeling allowed [3] to decipher a feldspar range varying between An30 and An40 (Fig. 1). Discussion: Although there could be a more complex history and other ways to form the whole compositional range of igneous minerals analyzed within the Bradbury formation, we are presenting here simple magmatic pathways commonly occurring on Earth using the thermodynamical softwares pMELTS and rhyoliteMELTS [14]. The objective is to find reasonable igneous processes that produce minerals that parallel the compositions of feldspar and pyroxene analyzed by the Curiosity rover. As commonly observed for mid-ocean ridge basalts, the adiabatic ascent of a primitive mantle composition [15] partially melting at 2 GPa has been modeled, followed by the extraction of a liquid at distinct degrees of partial melting, which undergoes fractional crystallization at an oxygen fugacity +1 log unit above the fayalite-magnetite-quartz (FMQ) buffer within the crust (0.02-0.4 GPa) with H2O = 0- 0.5 wt. %. These latter conditions correspond to those recorded within igneous clasts from the Noachian martian breccia NWA 7034 and paired and within Gale igneous rocks (colored dots in Fig. 1) [16-17]. To check the reliability of these 2-step models, we also tested fractional crystallization at similar conditions (FMQ+1; P=0.02-0.4 GPa; H2O = 0-0.5 wt. %) of starting compositions corresponding to that of magmas with distinct melting degrees obtained from isobaric experiments at 2 GPa [18]. Mineral compositions obtained from both models are similar. As shown on Fig. 1, the whole range of observed feldspar compositions cannot be reproduced by fractionation of one magma only. Indeed, while alkali feldspar and Na-plagioclase likely crystallized from fractional crystallization of a low-degree melt (here <15%), plagioclase and pyroxene can only be formed by fractional crystallization of a higher degree melt (here >19%). The corresponding liquid descent lines are broadly in agreement with compositions estimated by ChemCam corresponding to float igneous rocks (Fig. 2) [4-6]. Figure 2. Silica versus alkali content. Lines show the liquid lines of descent from magmas with distinct degrees of melting. Gray patches represents the composition of Gale igneous rocks [4-6]. Trachytic to rhyolitic magmas crystallize alkali feldspar, and andesite to dacite magmas likely form plagioclase. Therefore, at least two starting magmas at distinct melting degrees, which could easily come from a single mantle source, are necessary to explain the whole compositional range of feldspar and pyroxene analyzed within Bradbury rocks. Conclusion: Because rocks from the Bradbury formation are likely originating from the same magmatic source with minimal weathering as supported by several studies using different approaches, igneous mineral chemistry analyzed by CheMin and ChemCam allows us to back out reasonable magmatic pathways that could have crystallized them. Fractional crystallization of at least two starting magmas originating from distinct melting degrees of a single mantle source can explain the whole range of feldspar and pyroxene composition. Both alkaline and sub-alkaline liquids can be produced, with compositions corresponding to those of the igneous rocks analyzed by ChemCam within the Bradbury formation, highlighting the complexity of Mars magmatism.

Payre, V.↗

Convective Weather Avoidance with Uncertain Weather Forecasts

Convective weather events have a disruptive impact on air traffic both in terminal area and in en-route airspaces. In order to make sure that the national air transportation system is safe and efficient, it is essential to respond to convective weather events effectively. Traffic flow control initiatives in response to convective weather include ground delay, airborne delay, miles-in-trail restrictions as well as tactical and strategic rerouting. The rerouting initiatives can potentially increase traffic density and complexity in regions neighboring the convective weather activity. There is a need to perform rerouting in an intelligent and efficient way such that the disruptive effects of rerouting are minimized. An important area of research is to study the interaction of in-flight rerouting with traffic congestion or complexity and developing methods that quantitatively measure this interaction. Furthermore, it is necessary to find rerouting solutions that account for uncertainties in weather forecasts. These are important steps toward managing complexity during rerouting operations, and the paper is motivated by these research questions. An automated system is developed for rerouting air traffic in order to avoid convective weather regions during the 20- minute - 2-hour time horizon. Such a system is envisioned to work in concert with separation assurance (0 - 20-minute time horizon), and longer term air traffic management (2-hours and beyond) to provide a more comprehensive solution to complexity and safety management. In this study, weather is dynamic and uncertain; it is represented as regions of airspace that pilots are likely to avoid. Algorithms are implemented in an air traffic simulation environment to support the research study. The algorithms used are deterministic but periodically revise reroutes to account for weather forecast updates. In contrast to previous studies, in this study convective weather is represented as regions of airspace that pilots are likely to avoid. The automated system periodically updates forecasts and reassesses rerouting decisions in order to account for changing weather predictions. The main objectives are to reroute flights to avoid convective weather regions and determine the resulting complexity due to rerouting. The eventual goal is to control and reduce complexity while rerouting flights during the 20 minute - 2 hour planning period. A three-hour simulation is conducted using 4800 flights in the national airspace. The study compares several metrics against a baseline scenario using the same traffic and weather but with rerouting disabled. The results show that rerouting can have a negative impact on congestion in some sectors, as expected. The rerouting system provides accurate measurements of the resulting complexity in the congested sectors. Furthermore, although rerouting is performed only in the 20-minute - 2-hour range, it results in a 30% reduction in encounters with nowcast weather polygons (100% being the ideal for perfectly predictable and accurate weather). In the simulations, rerouting was performed for the 20-minute - 2-hour flight time horizon, and for the en-route segment of air traffic. The implementation uses CWAM, a set of polygons that represent probabilities of pilot deviation around weather. The algorithms were implemented in a software-based air traffic simulation system. Initial results of the system's performance and effectiveness were encouraging. Simulation results showed that when flights were rerouted in the 20-minute - 2-hour flight time horizon of air traffic, there were fewer weather encounters in the first 20 minutes than for flights that were not rerouted. Some preliminary results were also obtained that showed that rerouting will also increase complexity. More simulations will be conducted in order to report conclusive results on the effects of rerouting on complexity. Thus, the use of the 20-minute - 2-hour flight time horizon weather avoidance teniques performed in the simulation is expected to provide benefits for short-term weather avoidance.

Karahan, Sinan↗

Understanding the Scalability of Bayesian Network Inference using Clique Tree Growth Curves

Bayesian networks (BNs) are used to represent and efficiently compute with multi-variate probability distributions in a wide range of disciplines. One of the main approaches to perform computation in BNs is clique tree clustering and propagation. In this approach, BN computation consists of propagation in a clique tree compiled from a Bayesian network. There is a lack of understanding of how clique tree computation time, and BN computation time in more general, depends on variations in BN size and structure. On the one hand, complexity results tell us that many interesting BN queries are NP-hard or worse to answer, and it is not hard to find application BNs where the clique tree approach in practice cannot be used. On the other hand, it is well-known that tree-structured BNs can be used to answer probabilistic queries in polynomial time. In this article, we develop an approach to characterizing clique tree growth as a function of parameters that can be computed in polynomial time from BNs, specifically: (i) the ratio of the number of a BN's non-root nodes to the number of root nodes, or (ii) the expected number of moral edges in their moral graphs. Our approach is based on combining analytical and experimental results. Analytically, we partition the set of cliques in a clique tree into different sets, and introduce a growth curve for each set. For the special case of bipartite BNs, we consequently have two growth curves, a mixed clique growth curve and a root clique growth curve. In experiments, we systematically increase the degree of the root nodes in bipartite Bayesian networks, and find that root clique growth is well-approximated by Gompertz growth curves. It is believed that this research improves the understanding of the scaling behavior of clique tree clustering, provides a foundation for benchmarking and developing improved BN inference and machine learning algorithms, and presents an aid for analytical trade-off studies of clique tree clustering using growth curves.

Mengshoel, Ole Jakob↗

Automated Antenna Design with Evolutionary Algorithms

Current methods of designing and optimizing antennas by hand are time and labor intensive, and limit complexity. Evolutionary design techniques can overcome these limitations by searching the design space and automatically finding effective solutions. In recent years, evolutionary algorithms have shown great promise in finding practical solutions in large, poorly understood design spaces. In particular, spacecraft antenna design has proven tractable to evolutionary design techniques. Researchers have been investigating evolutionary antenna design and optimization since the early 1990s, and the field has grown in recent years as computer speed has increased and electromagnetic simulators have improved. Two requirements-compliant antennas, one for ST5 and another for TDRS-C, have been automatically designed by evolutionary algorithms. The ST5 antenna is slated to fly this year, and a TDRS-C phased array element has been fabricated and tested. Such automated evolutionary design is enabled by medium-to-high quality simulators and fast modern computers to evaluate computer-generated designs. Evolutionary algorithms automate cut-and-try engineering, substituting automated search though millions of potential designs for intelligent search by engineers through a much smaller number of designs. For evolutionary design, the engineer chooses the evolutionary technique, parameters and the basic form of the antenna, e.g., single wire for ST5 and crossed-element Yagi for TDRS-C. Evolutionary algorithms then search for optimal configurations in the space defined by the engineer. NASA's Space Technology 5 (ST5) mission will launch three small spacecraft to test innovative concepts and technologies. Advanced evolutionary algorithms were used to automatically design antennas for ST5. The combination of wide beamwidth for a circularly-polarized wave and wide impedance bandwidth made for a challenging antenna design problem. From past experience in designing wire antennas, we chose to constrain the evolutionary design to a monopole wire antenna. The results of the runs produced requirements-compliant antennas that were subsequently fabricated and tested. The evolved antenna has a number of advantages with regard to power consumption, fabrication time and complexity, and performance. Lower power requirements result from achieving high gain across a wider range of elevation angles, thus allowing a broader range of angles over which maximum data throughput can be achieved. Since the evolved antenna does not require a phasing circuit, less design and fabrication work is required. In terms of overall work, the evolved antenna required approximately three person-months to design and fabricate whereas the conventional antenna required about five. Furthermore, when the mission was modified and new orbital parameters selected, a redesign of the antenna to new requirements was required. The evolutionary system was rapidly modified and a new antenna evolved in a few weeks. The evolved antenna was shown to be compliant to the ST5 mission requirements. It has an unusual organic looking structure, one that expert antenna designers would not likely produce. This antenna has been tested, baselined and is scheduled to fly this year. In addition to the ST5 antenna, our laboratory has evolved an S-band phased array antenna element design that meets the requirements for NASA's TDRS-C communications satellite scheduled for launch early next decade. A combination of fairly broad bandwidth, high efficiency and circular polarization at high gain made for another challenging design problem. We chose to constrain the evolutionary design to a crossed-element Yagi antenna. The specification called for two types of elements, one for receive only and one for transmit/receive. We were able to evolve a single element design that meets both specifications thereby simplifying the antenna and reducing testing and integration costs. The highest performance antenna found using a getic algorithm and stochastic hill-climbing has been fabricated and tested. Laboratory results correspond well with simulation. Aerospace component design is an expensive and important step in space development. Evolutionary design can make a significant contribution wherever sufficiently fast, accurate and capable software simulators are available. We have demonstrated successful real-world design in the spacecraft antenna domain; and there is good reason to believe that these results could be replicated in other design spaces.

Hornby, Gregory S.↗

Integrated Ultra-Wideband Tracking and Carbon Dioxide Sensing System Design for International Space Station Applications

A three-dimensional (3D) Ultra-Wideband (UWB) Time-of-Arrival (TOA) tracking system has been studied at NASA Johnson Space Center (JSC) to provide the tracking capability inside the International Space Station (ISS) modules for various applications. One of applications is to locate and report the location where crew experienced possible high level of carbon-dioxide (CO2) and felt upset. Recent findings indicate that frequent, short-term crew exposure to elevated CO2 levels combined with other physiological impacts of microgravity may lead to a number of detrimental effects, including loss of vision. To evaluate the risks associated with transient elevated CO2 levels and design effective countermeasures, doctors must have access to frequent CO2 measurements in the immediate vicinity of individual crew members along with simultaneous measurements of their location in the space environment. To achieve this goal, a small, low-power, wearable system that integrates an accurate CO2 sensor with an ultra-wideband (UWB) radio capable of real-time location estimation and data communication is proposed. This system would be worn by crew members or mounted on a free-flyer and would automatically gather and transmit sampled sensor data tagged with real-time, high-resolution location information. Under the current proposed effort, a breadboard prototype of such a system has been developed. Although the initial effort is targeted to CO2 monitoring, the concept is applicable to other types of sensors. For the initial effort, a micro-controller is leveraged to integrate a low-power CO2 sensor with a commercially available UWB radio system with ranging capability. In order to accurately locate those places in a multipath intensive environment like ISS modules, it requires a robust real-time location system (RTLS) which can provide the required accuracy and update rate. A 3D UWB TOA tracking system with two-way ranging has been proposed and studied. The designed system will be tested in the Wireless Habitat Testbed which simulates the ISS module environment. This report describes the research and development effort for this prototype integrated UWB tracking and CO2 sensing system. The remainder of the report is organized as follows. In Section II, the TOA tracking methodology is introduced and the 3D tracking algorithm is derived. The simulation results are discussed in Section III. In Section VI, prototype system design and field tests are discussed. Some concluding remarks and future works are presented in Section V.

Ni, Jianjun (David)↗

Drop Size Distribution Measurements Supporting the NASA Global Precipitation Measurement Mission: Infrastructure and Preliminary Results

Global Precipitation Measurement Mission (GPM) retrieval algorithm validation requires datasets that characterize the 4-D structure, variability, and correlation properties of hydrometeor particle size distributions (PSD) and accumulations over satellite fields of view (5 -- 50 km). Key to this process is the combined use of disdrometer and polarimetric radar platforms. Here the disdrometer measurements serve as a reference for up-scaling dual-polarimetric radar observations of the PSD to the much larger volumetric sampling domain of the radar. The PSD observations thus derived provide a much larger data set for assessing DSD variability, and satellite-based precipitation retrieval algorithm assumptions, in all three spatial dimensions for a range of storm types and seasons. As one component of this effort, the GPM Ground Validation program recently acquired five 3rd generation 2D Video disdrometers as part of its Disdrometer and Radar Observations of Precipitation Facility (DROP), currently hosted in northern Alabama by the NASA Marshall Space Flight Center and the University of Alabama in Huntsville. These next-generation 2DVDs were operated and evaluated in different phases of data collection under the scanning domain of the UAH ARMOR C-band dual-polarimetric radar. During this period approximately 7500 minutes of PSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSDs were then used as input to T-matrix code for computation of polarimetric radar moments at C-band. The combined dataset was then analyzed with two basic objectives in mind: 1) the investigation of seasonal variability in the rain PSD parameters as observed by the 2DVDs; 2) the use of combined polarimetric moments and observed gamma distribution parameters in a functional form to retrieve PSD parameters in 4-D using the ARMOR radar for precipitation occurring in different seasons and for different rain system types. Preliminary results suggest that seasonal variations in the DSD parameters do occur, but are most pronounced when comparing tropical PSDs to either winter or summer convective precipitation. For example the previously documented shift to relatively smaller drop diameters in higher number concentrations for equivalent rain rate bins was observed in tropical storm rainbands occurring over Huntsville. On a more inter seasonal basis empirical fits between parameters such as D0 and ZDR do not appear to exhibit robust seasonal biases- i.e., one fit seems to work for all seasons within acceptable standard error (O[10%]) for estimates of D0. In polarimetric retrievals of the vertical variability in PSD (rain layer) for a tropical rainband we find that the Do varies with height when partitioned by specified precipitation categories (e.g., convective or stratiform, heavy and light stratiform etc.) but this variation is of order 10-20% and is smaller than the difference in D0 observed between the basic delineation of convective and stratiform precipitation types. Currently we are expanding our analysis of the vertical structure of the PSD to include several seasonally and/or dynamically-different storm system types (e.g., winter convection and stratiform events; summer mid-latitude convective etc.) sampled by ARMOR. The study will present the results of our combined analyses.

Petersen, Walter A.↗

Enhancements to the SSME transfer function modeling code

This report details the results of a one year effort by Ohio University to apply the transfer function modeling and analysis tools developed under NASA Grant NAG8-167 (Irwin, 1992), (Bartholomew, 1992) to attempt the generation of Space Shuttle Main Engine High Pressure Turbopump transfer functions from time domain data. In addition, new enhancements to the transfer function modeling codes which enhance the code functionality are presented, along with some ideas for improved modeling methods and future work. Section 2 contains a review of the analytical background used to generate transfer functions with the SSME transfer function modeling software. Section 2.1 presents the 'ratio method' developed for obtaining models of systems that are subject to single unmeasured excitation sources and have two or more measured output signals. Since most of the models developed during the investigation use the Eigensystem Realization Algorithm (ERA) for model generation, Section 2.2 presents an introduction of ERA, and Section 2.3 describes how it can be used to model spectral quantities. Section 2.4 details the Residue Identification Algorithm (RID) including the use of Constrained Least Squares (CLS) and Total Least Squares (TLS). Most of this information can be found in the report (and is repeated for convenience). Section 3 chronicles the effort of applying the SSME transfer function modeling codes to the a51p394.dat and a51p1294.dat time data files to generate transfer functions from the unmeasured input to the 129.4 degree sensor output. Included are transfer function modeling attempts using five methods. The first method is a direct application of the SSME codes to the data files and the second method uses the underlying trends in the spectral density estimates to form transfer function models with less clustering of poles and zeros than the models obtained by the direct method. In the third approach, the time data is low pass filtered prior to the modeling process in an effort to filter out high frequency characteristics. The fourth method removes the presumed system excitation and its harmonics in order to investigate the effects of the excitation on the modeling process. The fifth method is an attempt to apply constrained RID to obtain better transfer functions through more accurate modeling over certain frequency ranges. Section 4 presents some new C main files which were created to round out the functionality of the existing SSME transfer function modeling code. It is now possible to go from time data to transfer function models using only the C codes; it is not necessary to rely on external software. The new C main files and instructions for their use are included. Section 5 presents current and future enhancements to the XPLOT graphics program which was delivered with the initial software. Several new features which have been added to the program are detailed in the first part of this section. The remainder of Section 5 then lists some possible features which may be added in the future. Section 6 contains the conclusion section of this report. Section 6.1 is an overview of the work including a summary and observations relating to finding transfer functions with the SSME code. Section 6.2 contains information relating to future work on the project.

Irwin, R. Dennis↗

Application of neural network to remote sensing of soil moisture using theoretical polarimetric backscattering coefficients

This paper investigates the potential application of neural network to inversion of soil moisture using polarimetric remote sensing data. The neural network used for the inversion of soil parameters is multi-layer perceptron trained with the back-propagation algorithm. The training data include the polarimetric backscattering coefficients obtained from theoretical surface scattering models together with an assumed nominal range of soil parameters which are comprised of the soil permittivity and surface roughness parameters. Soil permittivity is calculated from the soil moisture and the assumed soil texture based on an empirical formula at C-, L-, and P-bands. The rough surface parameters for the soil surface, which is described by the Gaussian random process, are the root-mean-square (rms) height and correlation length. For the rough surface scattering, small perturbation method is used for the L-band frequency, and Kirchhoff approximation is used for the C-band frequency to obtain the corresponding backscattering coefficients. During the training, the backscattering coefficients are the inputs to the neural net and the output from the net are compared with the desired soil parameters to adjust the interconnecting weights. The process is repeated for each input-output data entry and then for the entire training data until convergence is reached. After training, the backscattering coefficients are applied to the trained neural net to retrieve the soil parameters which are compared with the desired soil parameters to verify the effectiveness of this technique. Several cases are examined. First, for simplicity, the correlation length and rms height of the soil surface are fixed while soil moisture is varied. Soil moisture obtained using the neural networks with either L-band or C-band backscattering coefficients for the HH and VV polarizations as inputs is in good agreement with the desired soil moisture. The neural net output matches the desired output for the soil moisture range of 16 to 60 percent for the C-band case. The next case investigated is to vary both soil moisture and rms height while keeping the correlation length fixed. For this case, C-band backscattering coefficients are not sufficient for retrieving two parameters because the Kirchhoff approximation gives the same HH and VV backscattering coefficients. Therefore, the backscattering coefficients at two different frequency bands are necessary to find both the soil moisture and rms height. Finally, the neural nets are also applied to simultaneously invert soil moisture, rms height, and correlation length. Overall, the soil moisture retrieved from the neural network agrees very well with the desired soil moisture. This suggests that the neural network shows potential for retrieval of soil parameters from remote sensing data.

Wang, L.↗

Simulation Technology at NASA

NASA Ames Research Center is home to several high-fidelity research flight and air-traffic control simulation facilities which, together with an experienced workforce, produce high-quality research data and findings that have proven to be applicable in the real world. These assets include the Vertical Motion Simulator (VMS), Crew Vehicle Systems Research Facility (CVSRF), Future Flight Central (FFC) air traffic control tower simulator, and several air-traffic control (ATC) simulators. The VMS combines a high-fidelity simulation capability with an adaptable simulation environment, enabling customization for numerous human-in-the-loop research applications. The distinctive feature of the VMS is its unparalleled large amplitude, high-fidelity motion capability. In over 30 years of continuous operation, the VMS has contributed significantly to the body of knowledge in a range of disciplines directly benefiting several aerospace programs and flight safety, including the design and development of flight control systems for the Joint Strike Fighter, Space Shuttle Orbiter, and rotorcraft. It continues to be used for researching new vehicle configurations, vehicle control and safety, transfer-of-training, etc., by NASA, other government agencies, and Industry.The CVSRF includes two motion-based flight simulators: a Boeing 747-400 full flight simulator and the reconfigurable Advanced Concepts Flight Simulator (ACFS). These simulators are primarily used to research air-traffic management concepts and procedures, advanced navigation and avionics concepts, and cockpit human factors. FFC is a full-sized control tower simulator with a 360-degree external field-of-view display system and reconfigurable system architecture. FFC and the ATC simulators are used for testing air-traffic management automation and decision support tools and demonstrate their feasibility in a realistic environment prior to technology transfer for implementation in the National Airspace System (NAS).To support integrated simulations and flight-tests for NASA's Unmanned Aircraft Systems (UAS) in the National Airspace System (NAS) Project, NASA developed a distributed test environment incorporating Live, Virtual, Constructive, (LVC) concepts. Development of the software enabling the LVC is conducted primarily at the Distributed Simulation Research Lab (DSRL) at NASA Ames. The LVC components provide the core infrastructure supporting simulation of UAS operations by integrating live and virtual aircraft in a realistic air traffic environment. This provides the ability to conduct tests more efficiently by promoting the use of existing distributed assets. The LVC infrastructure was used in several human-in-the-loop simulations to evaluate acceptance of Detect and Avoid (DAA) advisories used by UAS pilots to maintain well clear of other virtual traffic and to negotiate maneuvers with air traffic control. It is currently being used to support testing of self-separation algorithms between unmanned and manned aircraft in live flight. Further simulations with more comprehensive air traffic scenarios mixing live and virtual aircraft is planned.In the current fiscal environment, maintaining and upgrading these high-fidelity simulation assets and retaining the skilled workforce necessary to meet future research needs is the primary non-technical challenge.

CVSRF↗