SEARCH · Search NASA

Results for “graph model”

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.

SourceResearch connectionFrom yearThrough year

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

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Transonic airfoil calculations including wind tunnel wall-interference effects

The results of Reynolds-averaged time-dependent inviscid and turbulent compressible Navier-Stokes computations using the implicit finite-difference approach of Steger (1978), modified by incorporating a pressure boundary condition, (PBC) to account for wall interference are compared with experimental data on a NACA 64A010 airfoil (Johnson and Bachalo, 1980) in graphs and briefly characterized. The computational approach is the same as that used by King and Johnson (1980), but a 137 x 50 mesh is used instead of a 97 x 35 mesh, and special care is taken in resolving the nose, shock, and trailing-edge regions. Imposition of PBC is shown to improve significantly the accuracy of the computations for the flowfield on the upper surface of the airfoil, shifting the shock forward to its experimentally measured position in the case of turbulent flow. The failure of the method, even with PBC, to match the experimental shock location in the case of a flow with a separation bubble is attributed to inadequacies in the algebraic turbulence model employed (Baldwin and Lomax, 1978).

King, L. S.↗

NASA NTRS2007

Predicting Rocket or Jet Noise in Real Time

A semi-empirical theoretical model and a C++ computer program that implements the model have been developed for use in predicting the noise generated by a rocket or jet engine. The computer program, entitled the Realtime Rocket and Jet Engine Noise Analysis and Prediction Software, is one of two main subsystems of the Acoustic Prediction/Measurement Tool, which comprises software, acoustic instrumentation, and electronic hardware combined to afford integrated capabilities for real-time prediction and measurement of noise emitted by rocket and jet engines. [The other main subsystem, consisting largely of acoustic instrumentation and electronic hardware, is described in Wireless Acoustic Measurement System, which appears elsewhere in this section.] The theoretical model was derived from the fundamental laws of fluid mechanics, as first was done by M. J. Lighthill in his now famous theory of aerodynamically generated sound. The far-field approximation of the Lighthill theory is incorporated into this model. Many other contributions from various researchers have also been introduced into the model. The model accounts for two noise components: shear noise and self noise. The final result of the model is expressed in terms of a volume integral of the acoustic intensities attributable to these two components, subject to various directivity coefficients. The computer program was written to solve the volume integral. The inputs required by the program are two data files from a computational fluid dynamics (CFD) simulation of the flow of interest: the computational-grid file and the solution file. The CFD solution should be one that has been obtained for conditions that closely approximate those of an experimental test that is yet to be performed. In the current state of development of the model and software, it is recommended that the observation points lie along a radius at an angle >60 from the jet axis. The software provides, and is driven via, a graphical user interface, which facilitates its use. Optionally, the program accepts additional input in the form of data on the measured sound pressure level as a function of frequency at a given far-field location, preferably at an angle of 90 from the jet axis. The user is prompted to use default empirical constants or to choose constants based the measurement data. The user can view the results and compare them with other computational or experimental data. Once satisfied with the results, the user can save a graph of the results in a file that can be imported into documents.

Frendi, Kader↗

NASA NTRS1993

Structure in gamma ray burst time profiles: Correlations with other observables

One of the current debates raging in the world of gamma-ray burst physics is whether the sources of these enigmatic bursts arise from a single or from multiple distributions. Several authors contend that the histograms of GRB observables imply the latter. The two most-likely candidate components are galactic and cosmological. For example, some researchers claim that a dip in the V/V max distribution is a result of such a two-component source distribution. Others have used a parameter called the 'burst variability' calculated by dividing the maximum count rate on the 64-msec timescale by that from the 1024-msec timescale to show that a correlation of this parameter with bursts brightness implies a two-component model. This method has met vigorous criticism. We have developed two parameters that measure the variability or structure in the time profiles of BATSE gamma-ray bursts. Both parameters ('structure' and 'spike height') are based on the statistics of 'runs up' and 'runs down.' In short, the structure parameter is the observed number of runs (at several lengths) minus the number expected in a chance distribution. The 'spike height' is the sum of all run heights minus the expected sum. These two are straight-forward to calculate, robust, and measure the variability over the complete profile--not just at the peak. We have applied this algorithm to the profiles of 156 GRB's. In this paper we present graphs of the two parameters as functions of the following: (1) burst duration, (2) burst hardness ratio, (3) V/V max, (4) source galactic longitude, and (5) source galactic latitude. We seek correlations as well as groupings in the data that might indicate a multi-component source distribution.

Lestrade, John Patrick↗

NASA NTRS2011

Superpixel-Augmented Endmember Detection for Hyperspectral Images

Superpixels are homogeneous image regions comprised of several contiguous pixels. They are produced by shattering the image into contiguous, homogeneous regions that each cover between 20 and 100 image pixels. The segmentation aims for a many-to-one mapping from superpixels to image features; each image feature could contain several superpixels, but each superpixel occupies no more than one image feature. This conservative segmentation is relatively easy to automate in a robust fashion. Superpixel processing is related to the more general idea of improving hyperspectral analysis through spatial constraints, which can recognize subtle features at or below the level of noise by exploiting the fact that their spectral signatures are found in neighboring pixels. Recent work has explored spatial constraints for endmember extraction, showing significant advantages over techniques that ignore pixels relative positions. Methods such as AMEE (automated morphological endmember extraction) express spatial influence using fixed isometric relationships a local square window or Euclidean distance in pixel coordinates. In other words, two pixels covariances are based on their spatial proximity, but are independent of their absolute location in the scene. These isometric spatial constraints are most appropriate when spectral variation is smooth and constant over the image. Superpixels are simple to implement, efficient to compute, and are empirically effective. They can be used as a preprocessing step with any desired endmember extraction technique. Superpixels also have a solid theoretical basis in the hyperspectral linear mixing model, making them a principled approach for improving endmember extraction. Unlike existing approaches, superpixels can accommodate non-isometric covariance between image pixels (characteristic of discrete image features separated by step discontinuities). These kinds of image features are common in natural scenes. Analysts can substitute superpixels for image pixels during endmember analysis that leverages the spatial contiguity of scene features to enhance subtle spectral features. Superpixels define populations of image pixels that are independent samples from each image feature, permitting robust estimation of spectral properties, and reducing measurement noise in proportion to the area of the superpixel. This permits improved endmember extraction, and enables automated search for novel and constituent minerals in very noisy, hyperspatial images. This innovation begins with a graph-based segmentation based on the work of Felzenszwalb et al., but then expands their approach to the hyperspectral image domain with a Euclidean distance metric. Then, the mean spectrum of each segment is computed, and the resulting data cloud is used as input into sequential maximum angle convex cone (SMACC) endmember extraction.

Thompson, David R.↗

NASA NTRS2001

Pulse Detonation Engine Modeled

Pulse Detonation Engine Technology is currently being investigated at Glenn for both airbreathing and rocket propulsion applications. The potential for both mechanical simplicity and high efficiency due to the inherent near-constant-volume combustion process, may make Pulse Detonation Engines (PDE's) well suited for a number of mission profiles. Assessment of PDE cycles requires a simulation capability that is both fast and accurate. It should capture the essential physics of the system, yet run at speeds that allow parametric analysis. A quasi-one-dimensional, computational-fluid-dynamics-based simulation has been developed that may meet these requirements. The Euler equations of mass, momentum, and energy have been used along with a single reactive species transport equation, and submodels to account for dominant loss mechanisms (e.g., viscous losses, heat transfer, and valving) to successfully simulate PDE cycles. A high-resolution numerical integration scheme was chosen to capture the discontinuities associated with detonation, and robust boundary condition procedures were incorporated to accommodate flow reversals that may arise during a given cycle. The accompanying graphs compare experimentally measured and computed performance over a range of operating conditions for a particular PDE. Experimental data were supplied by Fred Schauer and Jeff Stutrud from the Air Force Research Laboratory at Wright-Patterson AFB and by Royce Bradley from Innovative Scientific Solutions, Inc. The left graph shows thrust and specific impulse, Isp, as functions of equivalence ratio for a PDE cycle in which the tube is completely filled with a detonable hydrogen/air mixture. The right graph shows thrust and specific impulse as functions of the fraction of the tube that is filled with a stoichiometric mixture of hydrogen and air. For both figures, the operating frequency was 16 Hz. The agreement between measured and computed values is quite good, both in terms of trend and magnitude. The error is under 10 percent everywhere except for the thrust value at an equivalence ratio of 0.8 in the left figure, where it is 14 percent. The simulation results shown were made using 200 numerical cells. Each cycle of the engine, approximately 0.06 sec, required 2.0 min of CPU time on a Sun Ultra2. The simulation is currently being used to analyze existing experiments, design new experiments, and predict performance in propulsion concepts where the PDE is a component (e.g., hybrid engines and combined cycles).

Paxson, Daniel E.↗

NASA NTRS1999

A Study of Direct and Cloud-Mediated Radiative Forcing of Climate Due to Aerosols

The Intergovernmental Panel on Climate Change (IPCC) has reported that in the southeastern US and eastern China, the general greenhouse warming due to anthropogenic gaseous emissions is dominated by the cooling effect of anthropogenic aerosols. To verify this model prediction in eastern China and southeastern US, we analyzed regional patterns of climate changes at 72 stations in eastern China during 1951- 94 (44 years), and at 52 stations in the southeastern US during 1949-94 (46 years) to detect the fingerprint of aerosol radiative forcing. It was found that the mean rates of change of annual mean daily, maximum, minimum temperatures and diurnal temperature range (DTR) in eastern China were 0.8, -0.2, 1.8, and -2.0 C/100 years respectively, while the mean rates of change of annual mean daily, maximum, minimum temperatures and DTR in the southeastern US were -0.2, -0.6, 0.2, and -0.8 C/100 years, respectively. This indicates that the high rate of increase in annual mean minimum temperature in eastern China results in a slightly warming trend of daily temperature, while the high rate of decrease in annual mean maximum temperature and low rate of increase in annual mean minimum temperature lead to the cooling trend of daily temperature in the southeastern US. We found that the warming from the longwave forcing due to both greenhouse gases and aerosols was completely counteracted by the shortwave aerosol forcing in the southeastern US in the past 46 years. A slightly overall warming trend in eastern China is evident; winters have become milder. This finding is explained by hypothesizing that increasing energy usage during the past 44 years has resulted in more coal and biomass burning, thus increasing the emission of absorbing soot and organic aerosols in eastern China. Such emissions, in addition to well-known Asia dust and greenhouse gases, may be responsible for the winter warming trend in eastern China that we have reported here. The sensitivity of aerosol radiative properties to aerosol composition, size distribution, relative humidity (RH) is examined for the following aerosol systems: inorganic and organic ions (Cl-, Br-, NO3 -, SO4 2-, Na+, NH4 +, K+, Ca2+, Mg2+, HCOO-, CH3COO-, CH3CH2COO-, CH3COCOO-, OOCCOO2-, MSA-1); water-insoluble inorganic and organic compounds (elemental carbon, n-alkanes, SiO2, Al2O3, Fe2O3 and other organic compounds). The partial molar refraction method was used to calculate the real part of the refractive index. It was found that the asymmetry factor increased by approximately 48% with the real part varying from 1.40 to 1.65, and the single scattering albedo decreased by 24% with the imaginary part varying from -0.005 to -0.1. The asymmetry factor increased by 5.4 times with the geometric standard deviation varying from 1.2 to 3.0. The radiation transmission is very sensitive to the change in size distribution; other factors are not as significant. To determine the aerosol direct radiative forcing (ADRF), the aerosol optical depth (AOD) values at the three operational wavelengths (415, 500 and 673 nm) were determined at a regionally representative site, namely, Mt. Gibbs (35.78 deg N, 82.29 deg W, elevation 2006 m) in Mt. Mitchell State Park, NC, and a site located in an adjacent valley (Black Mountain, 35.66 deg N, 82.38 deg W, elevation 951 m) in the southeastern US. The two sites are separated horizontally by 10 km and vertically by 1 km. It was found that the representative total AOD values at 500 nm at the valley site for highly polluted (HP), marine (M) and continental (C) air masses were 0.68 +/- 0.33, 0.29 +/- 0.19 and 0.10 +/- 0.04, respectively. A search-graph method was used to retrieve the columnar size distribution (number concentration N, effective radius reff and geometric standard deviation=?g) from the optical depth observations at three operational wavelengths. The ground albedo, single scattering albedo and imaginary part of the refractive index were calculated using a mathematically unique procedure involving a Mie code and a radiative transfer code in conjunction with the retrieved aerosol size distribution, AOD, and diffuse-direct irradiance ratio. It was found that N, r(eff) and sigma(g) were in the ranges of 10 to 1.7 x 10(exp 4)/cubic cm, 0.09 to 0.68 micrometers and 1.12 to 2.95, respectively. The asymmetry factor and single scattering albedo were in the ranges of 0.63 to 0.75 and 0.74 to 0.97 respectively. The ground albedo for the forested terrain and imaginary part of refractive index were found to be in the ranges of 0.06 to 0.29 and 0.005 to 0.051 respectively. On the basis of these aerosol radiative properties obtained at the research sites and computations using the Column Radiation Model (CRM) of National Center of Atmospheric Research (NCAR) Community Climate Model (CCM3), it was found that the average cloud-free 24-hour ADRF values were -13 +/- 8, -8 +/- 3, -33 +/- 16 W/square m for marine, continental, and polluted air masses, respectively. On the assumption that the fractional coverage of clouds is 0.61, it was estimated that the annual mean ADRF was 7 +/- 2 W/square m in the southeastern US. The review with respect to the current knowledge of organic acids shows that aerosol formate and acetate concentrations range from 0.02 to 5.3 nmol/cubic m and from 0.03 to 12.4 nmol/cubic m respectively, and that between 34% to 77% of formate and between 21% to 66% of acetate are present in the fine fraction of aerosols. It was found that although most (98-99%) of these volatile organic acids were present in the gas phase, their concentrations in the aerosol particles were sufficient to make them a good candidate for cloud condensation nuclei (CCN). It is hypothesized that organic acids are at least one of the primary sources of CCN in the atmosphere due to their ubiquitous presence in the troposphere, especially over the continental forested areas. The results of our measurements at Palmer Station, Antarctica show that the daily average CCN concentrations at 0.3% and 1% supersaturations ranged from 0.3 to 160/cubic cm and from 4 to 168/cubic cm, respectively, during the period from 17 January to 26 February, 1994. New evidence for substantial and definitive CCN enhancement near and within cloud has been observed at Mt. Mitchell, North Carolina. The results show that the average monthly CCN concentrations were 460 +/- 217, 386 +/- 286, 429 +/- 228 and 238 +/- 134/cubic cm for in-cloud, overcast, clear and rainy conditions, respectively. The typical CCN spectra show that there were a lot of small CCN produced and the ion concentrations (especially H+ and SO4 2-) were very high during the CCN enhancement period. The significantly positive correlation between black carbon (BC) and CCN at 1% supersaturation indicates that a percentage of the BC measured at the site may be in the form of an internal mixture and participated in the formation of CCN.

Yu, Shao-Cai↗

NASA NTRS2002

Reducing False Positives in Runtime Analysis of Deadlocks

This paper presents an improvement of a standard algorithm for detecting dead-lock potentials in multi-threaded programs, in that it reduces the number of false positives. The standard algorithm works as follows. The multi-threaded program under observation is executed, while lock and unlock events are observed. A graph of locks is built, with edges between locks symbolizing locking orders. Any cycle in the graph signifies a potential for a deadlock. The typical standard example is the group of dining philosophers sharing forks. The algorithm is interesting because it can catch deadlock potentials even though no deadlocks occur in the examined trace, and at the same time it scales very well in contrast t o more formal approaches to deadlock detection. The algorithm, however, can yield false positives (as well as false negatives). The extension of the algorithm described in this paper reduces the amount of false positives for three particular cases: when a gate lock protects a cycle, when a single thread introduces a cycle, and when the code segments in different threads that cause the cycle can actually not execute in parallel. The paper formalizes a theory for dynamic deadlock detection and compares it to model checking and static analysis techniques. It furthermore describes an implementation for analyzing Java programs and its application to two case studies: a planetary rover and a space craft altitude control system.

Bensalem, Saddek↗

NASA NTRS1990

BATSE software for the analysis of the gamma ray burst spatial distribution

The Burst and Transient Source Experiment (BATSE) on the Gamma Ray Observatory (GRO) is designed to study astronomical gamma ray sources and to provide better positional, spectral, and time resolution about these objects than has previously been possible from one experiment. The procedure to be used in the analysis of the gamma ray burst spatial distribution is presented. Data is input from BATSE via the Gamma Ray Burst Catalog (listing individual burst positions, flux values, and associated errors) and the Sky Sensitivity Map (which summarizes observational selection effects in table format). A FORTRAN program generates Monte Carlo burst catalogs, which are models to be compared to the actual distribution. The Monte Carlo models are then filtered through the Sky Sensitivity Map so that they suffer from the same selection effects as the actual catalog data. Additionally, each burst position is converted into a probability distribution to mimic BATSE positional sensitivity. The Burst Catalog, Monte Carlo burst catalog, and Sky Sensitivity Map are then passed onto an IDL program that compares the catalogs for statistical significance. The Sky Sensitivity Map is used to estimate how often each sky area is observed above the minimum flux level in question. Each burst found in this sky area is then weighted according to the frequency with which this sky area is observed. The catalogs are then compared via tests of homogeneity (based on their radial distributions) and isotropy (based upon their angular distributions). The results of the statistical comparisons along with graphs and charts of the summaries, are output from the IDL program for study.

Hakkila, Jon↗

NASA NTRS2008

Potential Use of a Bayesian Network for Discriminating Flash Type from Future GOES-R Geostationary Lightning Mapper (GLM) data

Continuous monitoring of the ratio of cloud flashes to ground flashes may provide a better understanding of thunderstorm dynamics, intensification, and evolution, and it may be useful in severe weather warning. The National Lighting Detection Network TM (NLDN) senses ground flashes with exceptional detection efficiency and accuracy over most of the continental United States. A proposed Geostationary Lightning Mapper (GLM) aboard the Geostationary Operational Environmental Satellite (GOES-R) will look at the western hemisphere, and among the lightning data products to be made available will be the fundamental optical flash parameters for both cloud and ground flashes: radiance, area, duration, number of optical groups, and number of optical events. Previous studies have demonstrated that the optical flash parameter statistics of ground and cloud lightning, which are observable from space, are significantly different. This study investigates a Bayesian network methodology for discriminating lightning flash type (ground or cloud) using the lightning optical data and ancillary GOES-R data. A Directed Acyclic Graph (DAG) is set up with lightning as a "root" and data observed by GLM as the "leaves." This allows for a direct calculation of the joint probability distribution function for the lighting type and radiance, area, etc. Initially, the conditional probabilities that will be required can be estimated from the Lightning Imaging Sensor (LIS) and the Optical Transient Detector (OTD) together with NLDN data. Directly manipulating the joint distribution will yield the conditional probability that a lightning flash is a ground flash given the evidence, which consists of the observed lightning optical data [and possibly cloud data retrieved from the GOES-R Advanced Baseline Imager (ABI) in a more mature Bayesian network configuration]. Later, actual GLM and NLDN data can be used to refine the estimates of the conditional probabilities used in the model; i.e., the Bayesian network is a learning network. Methods for efficient calculation of the conditional probabilities (e.g., an algorithm using junction trees), finding data conflicts, goodness of fit, and dealing with missing data will also be addressed.

Solakiewiz, Richard↗

NASA NTRS1992

Internal reversing flow in a tailpipe offtake configuration for SSTOVL aircraft

A generic one-third scale model of a tailpipe offtake system for a supersonic short takeoff vertical landing (SSTOVL) aircraft was tested at LeRC Powered Lift Facility. The model consisted of a tailpipe with twin elbows, offtake ducts, and flow control nozzles, plus a small ventral nozzle and a blind flange to simulate a blocked cruise nozzle. The offtake flow turned through a total angle of 177 degrees relative to the tailpipe inlet axis. The flow split was 45 percent to each offtake and 10 percent to the ventral nozzle. The main test objective was to collect data for comparison to the performance of the same configuration predicted by a computational fluid dynamics (CFD) analysis. Only the experimental results are given - the analytical results are published in a separate paper. Performance tests were made with unheated air at tailpipe-to-ambient pressure ratios up to 5. The total pressure loss through the offtakes was as high as 15.5 percent. All test results are shown as graphs, contour plots, and wall pressure distributions. The complex flow patterns in the tailpipe and elbows at the offtake openings are described with traversing flow angle probe and paint streak flow visualization data.

Mcardle, Jack G.↗

NASA NTRS1992

Internal reversing flow in a tailpipe offtake configuration for SSTOVL aircraft

Mcardle, Jack G.↗

NASA NTRS1994

Separated and Recovering Turbulent Boundary Layer Flow Behind a Backward Facing Step For Different Reynolds Numbers

Experimental results for a two-dimensional separated turbulent boundary layer behind a backward facing step for five different Reynolds numbers are reported. Results are presented in the form of tables, graphs and a floppy disk for an easy access of the data. Reynolds number based on the step height was varied by changing the reference velocity upstream of the step, U(sub o), and the step height, h. Hot-wire measurement techniques were used to measure three Reynolds stresses and four triple-velocity correlations. In addition, surface pressure and skin friction coefficients were measured. All hot-wire measurements were acquired in a measuring domain which excluded recirculating flow region due to the directional insensitivity of hot-wires. The downstream extent of the domain from the step was 51 h for the largest and I 14h for the smallest step height. This significant downstream length permitted extensive study of the flow recovery. Prediction of perturbed flows and their recovery is particularly attractive for popular turbulence models since variations of turbulence length and time scales and flow interactions in different regions are generally inadequately predicted. The data indicate that the flow in the free shear layer region behaves like the plane mixing layer up to about 2/3 of the mean reattachment length when the flow interaction with the wall commences the flow recovery to that of an ordinary turbulent boundary layer structure. These changes of the flow do not occur abruptly with the change of boundary conditions. A reattachment region represents a transitional region where the flow undergoes the most dramatic adjustments to the new boundary conditions. Large eddies, created in the upstream free-shear layer region, are being torn, recirculated, reentrained back into the main stream interacting with the incoming flow structure. It is foreseeable that it is quite difficult to describe the physics of this region in a rational and quantitative manner other than statistical. Downstream of the reattachment point the flow recovers at different rates near the wall, in the newly developing internal boundary layer, and in the outer part of the flow. It appears that Reynolds stresses do not fully recover up to the longest recovery length of 114 h.

Jovic, Srba↗

NASA NTRS1977

A life cycle cost economics model for projects with uniformly varying operating costs

A mathematical model is developed for calculating the life cycle costs for a project where the operating costs increase or decrease in a linear manner with time. The life cycle cost is shown to be a function of the investment costs, initial operating costs, operating cost gradient, project life time, interest rate for capital and salvage value. The results show that the life cycle cost for a project can be grossly underestimated (or overestimated) if the operating costs increase (or decrease) uniformly over time rather than being constant as is often assumed in project economic evaluations. The following range of variables is examined: (1) project life from 2 to 30 years; (2) interest rate from 0 to 15 percent per year; and (3) operating cost gradient from 5 to 90 percent of the initial operating costs. A numerical example plus tables and graphs is given to help calculate project life cycle costs over a wide range of variables.

Remer, D. S.↗

NASA NTRS1974

A study of remote sensing as applied to regional and small watersheds. Volume 2: Supporting technical details

The Stanford Watershed Model, the Kentucky Watershed Model and OPSET program, and the NASA-IBM system for simulation and analysis of watersheds are described in terms of their applications to the study of remote sensing of water resources. Specific calibration processes and input and output parameters that are instrumental in the simulations are explained for the following kinds of data: (1) hourly precipitation data; (2) daily discharge data; (3) flood hydrographs; (4) temperature and evaporation data; and (5) snowmelt data arrays. The Sensitivity Analysis Task, which provides a method for evaluation of any of the separate simulation runs in the form of performance indices, is also reported. The method is defined and a summary of results is given which indicates the values obtained in the simulation runs performed for Town Creek, Alabama; Alamosa Creek, Colorado; and Pearl River, Louisiana. The results are shown in tabular and plot graph form. For Vol. 1, see N74-27813.

Source record↗

NASA NTRS2002

Contingency Planning for Planetary Rovers

There has been considerable work in AI on planning under uncertainty. But this work generally assumes an extremely simple model of action that does not consider continuous time and resources. These assumptions are not reasonable for a Mars rover, which must cope with uncertainty about the duration of tasks, the power required, the data storage necessary, along with its position and orientation. In this paper, we outline an approach to generating contingency plans when the sources of uncertainty involve continuous quantities such as time and resources. The approach involves first constructing a "seed" plan, and then incrementally adding contingent branches to this plan in order to improve utility. The challenge is to figure out the best places to insert contingency branches. This requires an estimate of how much utility could be gained by building a contingent branch at any given place in the seed plan. Computing this utility exactly is intractable, but we outline an approximation method that back propagates utility distributions through a graph structure similar to that of a plan graph.

Dearden, Richard↗

NASA NTRS2003

Incremental Contingency Planning

There has been considerable work in AI on planning under uncertainty. However, this work generally assumes an extremely simple model of action that does not consider continuous time and resources. These assumptions are not reasonable for a Mars rover, which must cope with uncertainty about the duration of tasks, the energy required, the data storage necessary, and its current position and orientation. In this paper, we outline an approach to generating contingency plans when the sources of uncertainty involve continuous quantities such as time and resources. The approach involves first constructing a "seed" plan, and then incrementally adding contingent branches to this plan in order to improve utility. The challenge is to figure out the best places to insert contingency branches. This requires an estimate of how much utility could be gained by building a contingent branch at any given place in the seed plan. Computing this utility exactly is intractable, but we outline an approximation method that back propagates utility distributions through a graph structure similar to that of a plan graph.

Dearden, Richard↗

NASA NTRSDate not supplied

Aerial Vehicle Routing and Scheduling for UAS Traffic Management: A Hybrid Monte Carlo Tree Search Approach

We present the Multi-Route Weighted Package Delivery Problem (MRWPDP) and a scalable solution methodology as a major step towards enabling an airspace deconfliction service for drone delivery operations. The problem is motivated by Strategic deconfliction under the FAA’s “Unmanned Aircraft Systems Traffic Management” Concept of Operations. MRWPDP falls under a class of vehicle routing and scheduling problems, and as such is NP-Hard. In MRWPDP, a graph network is given which consists of depots, drop-off sites, and multiple routes connecting the two. In addition, routes are weighted by the associated ground risk and total travel distance for package delivery. The goal is to optimally schedule the departure time and assign routes to a known set of vehicles at the depot. We propose a heuristic solution to the problem by borrowing techniques from Mixed Integer Linear Programming (MILP), Constraint Programming, and Monte Carlo Tree Search (MCTS). The resulting hybrid framework is MCTS with Bound-and-Prune (BP) and rapid simulated updates (U), or MCTS-BP-U. This approach is able to quickly provide a feasible solution for MRWPDP, even for large problem instances up to 1000 vehicles. We provide a MILP formulation of MRWPDP and compare its performance against MCTS-BP-U in terms of solution quality. An agent-based model simulation is conducted as a final step to validate the efficacy of our approach.

air traffic scheduling↗

NASA NTRSDate not supplied

Aerial Vehicle Routing and Scheduling for UAS Traffic Management: A Hybrid Monte Carlo Tree Search Approach

air traffic scheduling↗