Search NASA⌕ Search

SEARCH · Search NASA

Results for “exponential advantage”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Prediction and causal reasoning in planning

Nonlinear planners are often touted as having an efficiency advantage over linear planners. The reason usually given is that nonlinear planners, unlike their linear counterparts, are not forced to make arbitrary commitments to the order in which actions are to be performed. This ability to delay commitment enables nonlinear planners to solve certain problems with far less effort than would be required of linear planners. Here, it is argued that this advantage is bought with a significant reduction in the ability of a nonlinear planner to accurately predict the consequences of actions. Unfortunately, the general problem of predicting the consequences of a partially ordered set of actions is intractable. In gaining the predictive power of linear planners, nonlinear planners sacrifice their efficiency advantage. There are, however, other advantages to nonlinear planning (e.g., the ability to reason about partial orders and incomplete information) that make it well worth the effort needed to extend nonlinear methods. A framework is supplied for causal inference that supports reasoning about partially ordered events and actions whose effects depend upon the context in which they are executed. As an alternative to a complete but potentially exponential-time algorithm, researchers provide a provably sound polynomial-time algorithm for predicting the consequences of partially ordered events.

Dean, T.↗
NASA NTRS2015

Pioneering a Biobased UAS

With the exponential growth of interest in unmanned aerial vehicles (UAVs) and their vast array of applications in both space exploration and terrestrial uses such as the delivery of medicine and monitoring the environment, the 2014 Stanford-Brown-Spelman iGEM team is pioneering the development of a fully biological UAV for scientific and humanitarian missions. The prospect of a biologically-produced UAV presents numerous advantages over the current manufacturing paradigm. First, a foundational architecture built by cells allows for construction or repair in locations where it would be difficult to bring traditional tools of production. Second, a major limitation of current research with UAVs is the size and high power consumption of analytical instruments, which require bulky electrical components and large fuselages to support their weight. By moving these functions into cells with biosensing capabilities - for example, a series of cells engineered to report GFP, green fluorescent protein, when conditions exceed a certain threshold concentration of a compound of interest, enabling their detection post-flight - these problems of scale can be avoided. To this end, we are working to engineer cells to synthesize cellulose acetate as a novel bioplastic, characterize biological methods of waterproofing the material, and program this material's systemic biodegradation. In addition, we aim to use an "amberless" system to prevent horizontal gene transfer from live cells on the material to microorganisms in the flight environment. So far, we have: successfully transformed Gluconacetobacter hansenii, a cellulose-producing bacterium, with a series of promoters to test transformation efficiency before adding the acetylation genes; isolated protein bands present in the wasp nest material; transformed the cellulose-degrading genes into Escherichia coli; and we have confirmed that the amberless construct prevents protein expression in wild-type cells. In addition, as part of our human outreach project, we have been in touch with leaders in the fields of UAVs, synthetic biology, and earth sciences, and it is clear that biodegradable UAVs could have a significant impact on the industry.

Escherichia↗
NASA NTRS2015

Pioneering a Biobased UAS

With the exponential growth of interest in unmanned aerial vehicles (UAVs) and their vast array of applications in both space exploration and terrestrial uses such as the delivery of medicine and monitoring the environment, the 2014 Stanford-Brown-Spelman iGEM team is pioneering the development of a fully biological UAV for scientific and humanitarian missions. The prospect of a biologically-produced UAV presents numerous advantages over the current manufacturing paradigm. First, a foundational architecture built by cells allows for construction or repair in locations where it would be difficult to bring traditional tools of production. Second, a major limitation of current research with UAVs is the size and high power consumption of analytical instruments, which require bulky electrical components and large fuselages to support their weight. By moving these functions into cells with biosensing capabilities – for example, a series of cells engineered to report GFP, green fluorescent protein, when conditions exceed a certain threshold concentration of a compound of interest, enabling their detection post-flight – these problems of scale can be avoided. To this end, we are working to engineer cells to synthesize cellulose acetate as a novel bioplastic, characterize biological methods of waterproofing the material, and program this material’s systemic biodegradation. In addition, we aim to use an “amberless” system to prevent horizontal gene transfer from live cells on the material to microorganisms in the flight environment. So far, we have: successfully transformed Gluconacetobacter hansenii, a cellulose-producing bacterium, with a series of promoters to test transformation efficiency before adding the acetylation genes; isolated protein bands present in the wasp nest material; transformed the cellulose-degrading genes into Escherichia coli; and we have confirmed that the amberless construct prevents protein expression in wild-type cells. In addition, as part of our human outreach project, we have been in touch with leaders in the fields of UAVs, synthetic biology, and earth sciences, and it is clear that biodegradable UAVs could have a significant impact on the industry.

Escherichia↗
NASA NTRS1994

Kriging in the Shadows: Geostatistical Interpolation for Remote Sensing

It is often useful to estimate obscured or missing remotely sensed data. Traditional interpolation methods, such as nearest-neighbor or bilinear resampling, do not take full advantage of the spatial information in the image. An alternative method, a geostatistical technique known as indicator kriging, is described and demonstrated using a Landsat Thematic Mapper image in southern Chiapas, Mexico. The image was first classified into pasture and nonpasture land cover. For each pixel that was obscured by cloud or cloud shadow, the probability that it was pasture was assigned by the algorithm. An exponential omnidirectional variogram model was used to characterize the spatial continuity of the image for use in the kriging algorithm. Assuming a cutoff probability level of 50%, the error was shown to be 17% with no obvious spatial bias but with some tendency to categorize nonpasture as pasture (overestimation). While this is a promising result, the method's practical application in other missing data problems for remotely sensed images will depend on the amount and spatial pattern of the unobscured pixels and missing pixels and the success of the spatial continuity model used.

Rossi, Richard E.↗
NASA NTRS2003

New Active Optical Technique Developed for Measuring Low-Earth-Orbit Atomic Oxygen Erosion of Polymers

Polymers such as polyimide Kapton (DuPont) and Teflon FEP (DuPont, fluorinated ethylene propylene) are commonly used spacecraft materials because of desirable properties such as flexibility, low density, and in the case of FEP, a low solar absorptance and high thermal emittance. Polymers on the exterior of spacecraft in the low-Earth-orbit (LEO) environment are exposed to energetic atomic oxygen. Atomic oxygen reaction with polymers causes erosion, which is a threat to spacecraft performance and durability. It is, therefore, important to understand the atomic oxygen erosion yield E (the volume loss per incident oxygen atom) of polymers being considered in spacecraft design. The most common technique for determining E is a passive technique based on mass-loss measurements of samples exposed to LEO atomic oxygen during a space flight experiment. There are certain disadvantages to this technique. First, because it is passive, data are not obtained until after the flight is completed. Also, obtaining the preflight and postflight mass measurements is complicated by the fact that many polymers absorb water and, therefore, the mass change due to water absorption can affect the E data. This is particularly true for experiments that receive low atomic oxygen exposures or for samples that have a very low E. An active atomic oxygen erosion technique based on optical measurements has been developed that has certain advantages over the mass-loss technique. This in situ technique can simultaneously provide the erosion yield data on orbit and the atomic oxygen exposure fluence, which is needed for erosion yield determination. In the optical technique, either sunlight or artificial light can be used to measure the erosion of semitransparent or opaque polymers as a result of atomic oxygen attack. The technique is simple and adaptable to a rather wide range of polymers, providing that they have a sufficiently high optical absorption coefficient. If one covers a photodiode with a uniformly thick sheet of semitransparent polymer such as Kapton H polyimide, then as atomic oxygen erodes the polymer, the short-circuit current from the photodiode will increase in an exponential manner with fluence. This nonlinear response with fluence results in a lack of sensitivity for measuring low atomic oxygen fluences. However, if one uses a variable-thickness polymer or carbon sample, which is configured as shown in the preceding figure, then a linear response can be achieved for opaque materials using a parabolic well for a circular geometry detector or a V-shaped well for a rectangular-geometry detector. Variable-thickness samples can be fabricated using many thin polymer layers. For semitransparent polymers such as Kapton H polyimide, there is an initial short-circuit current that is greater than zero. This current has a slightly nonlinear dependence on atomic oxygen fluence in comparison to opaque materials such as black Kapton as shown in the graph. For this graph figure, the total thickness of Kapton H was assumed to be 0.03 cm. The photodiode short-circuit current shown in the graph was generated on the basis of preliminary measurements-a total reflectance rho of 0.0424 and an optical absorption coefficient a of 146.5 cm(sup -1). In addition to obtaining on-orbit data, the advantage of this active erosion and erosion yield measurement technique is its simplicity and reliance upon well-characterized fluence witness materials as well as a nearly linear photodiode short-circuit current dependence upon atomic oxygen fluence. The optical technique is useful for measuring either atomic oxygen fluence or erosion, depending on the information desired. To measure the atomic oxygen erosion yield of a test material, one would need to have two photodiode sensors, one for the test material and one that uses a known erosion yield material (such as Kapton) to measure the atomic oxygen fluence.

Banks, Bruce A.↗
NASA NTRS1999

Simulating the Formation of Lunar Floor-Fracture Craters Using Elastoviscoplastic Relaxation

Lunar floor-fracture craters formed during the height of mare basalt emplacement. Due to a general temporal and spatial relation with the maria, these craters, numbering some 200, may be diagnostic of the thermal structure of the crust during this time. As the name suggests, these craters exhibit brittle failure, generally limited to the central floor region. That, and a shallower depth than fresh lunar craters, has led to two main theories as to their formation: laccolith emplacement under the crater and viscous relaxation. The implications of each model for the state of the Moon's crust during this time are quite different, so the viability of each model must be checked. Laccolith emplacement has been treated elsewhere. However, previous attempts to study the relaxation of the craters have assumed only a uniform, Newtonian viscous response of the near surface to the topographic driving forces, and simply postulated that the fractures resulted from tensile stresses associated with floor uplift. Here, we use a more sophisticated rheological model that includes not only non-Newtonian viscous behavior (i.e., the viscosity is stress-dependent), but also incorporates elastic behavior and a plastic component to the rheology to directly simulate the formation of the floor fractures. The results of our simulations show that while elastoviscoplastic relaxation is potentially viable for larger floor-fracture craters, it is not viable for craters with diameters < or = 60 km, the size of the majority of floor-fracture craters. We employ the finite element method, a numerical technique well suited for boundary-value problems, via the commercially available MARC software package. To test the viability of topographic relaxation, our goal is to prepare the simulations as to maximize the amount of relaxation. We take advantage of the natural axisymmetry of craters, simulating one radial plane. Initial shapes are based on data for fresh craters from Pike. To simplify implementation, a fourth order polynomial is used for the basin, while a third order inverse function is used for the rim. This form closely approximates the long-wavelength behavior of complex craters, while ignoring higher-frequency topography, save the rim. This approximation is appropriate because crater relaxation is strongly controlled by long-wavelength topography. Loading is accomplished assuming a uniform gravity field (1.62 m/s-square) and a uniform density of 2900 kg/cubic m. The initial stress state is set to be hydrostatic, with an additional pressure term to account for any overlying topography. This additional pressure term is tapered exponentially with depth. While the simulations quickly settle on a preferred stress state, and while the final solution is fairly insensitive to the choice of the e-folding depth of the taper, selecting an e-folding depth close to the diameter of the crater sets the initial stress state near the preferred state. We assume a diuranally averaged surface temperature of -20 C, and allow temperature to increase with depth at a rate of 50 K//km. Assuming a thermal conductivity of 2 W/in/K, this gradient translates to a heat flow of 100mW/square m, an extremely high value for the Moon. Temperature, of course, will not increase without bound. To maximize relaxation, we allow our temperature profile to increase linearly until it reaches the solidus (assumed to be 1200C) at a depth of 24.4 km, at which point it is kept constant. The presence of melt will drop the bulk viscosity; however, we have no rheological control for partial melts. Therefore, we make no attempt to simulate this situation. Elastoviscoplastic rheological model. In general, geologic materials can behave in three main ways: elastically, viscously (via solid-state creep), and brittly (plasticity is a continuum approach to simulate this phenomenon). We combine these three deformation mechanisms in an extended Maxwell solid, where the total strain can be broken down into a simple summation of the elastic, creep, and plastic strains. In relaxation phenomena in general, the system takes advantage of any means possible to eliminate deviatoric stresses by relaxing away the topography. Previous analyses have only modeled the viscous response. Comparatively, the elastic response in our model can augment the relaxation, to a point. This effect decreases as the elastic response becomes stiffer; indeed, in the limit of infinite elastic Young's modulus (and with no plasticity), the solution converges on the purely viscous solution. Igneous rocks common to the lunar near-surface have Young's modulii in the range of 10-100 GPa. To maximize relaxation, we use a Young's modulus of 10 GPa. (There is negligible sensitivity to the other elastic modulus, the Poisson's ratio; we use 0.25.) For the viscous response, we use a flow law for steady-state creep in thoroughly dried Columbia diabase, because the high plagioclase (about 70 vol%) and orthopyroxene (about 17 vol%) content is similar to the composition of the lunar highland crust as described by remote sensing and sample studies: noritic anorthosite. This flow law is highly non-Newtonian, i.e., the viscosity is highly stress dependent. That, and the variability with temperature, stands in strong contrast to previous examinations of lunar floor-fracture crater relaxation. To model discrete, brittle faulting, we assume "Byerlee's rule," a standard geodynamical technique. We implement this "rule" with an-angle of internal friction of about 40 deg, and a higher-than-normal cohesion of about 3.2 MPa (to approximate the breaking of unfractured rock). The actual behavior of geologic materials is more complex than in our rheological model, so the uncertainties in the plasticity do not represent the state-of-the-art error. Additional information is contained in the original.

Dombard, A. J.↗
NASA NTRS2011

Model Adaptation for Prognostics in a Particle Filtering Framework

One of the key motivating factors for using particle filters for prognostics is the ability to include model parameters as part of the state vector to be estimated. This performs model adaptation in conjunction with state tracking, and thus, produces a tuned model that can used for long term predictions. This feature of particle filters works in most part due to the fact that they are not subject to the "curse of dimensionality", i.e. the exponential growth of computational complexity with state dimension. However, in practice, this property holds for "well-designed" particle filters only as dimensionality increases. This paper explores the notion of wellness of design in the context of predicting remaining useful life for individual discharge cycles of Li-ion batteries. Prognostic metrics are used to analyze the tradeoff between different model designs and prediction performance. Results demonstrate how sensitivity analysis may be used to arrive at a well-designed prognostic model that can take advantage of the model adaptation properties of a particle filter.

Saha, Bhaskar↗
NASA NTRS1987

Temperature characteristics of silicon space solar cells and underlying parameters

Silicon space cells, 2 cm x 2 cm, with 10 ohm-cm p-base resistivity, 8-mil base thickness, and no back-surface fields have been investigated over the temperature range from 301 to 223 K by measurements of dark forward and reverse current-voltage characteristics and current-voltage relations under illumination. From dark forward bias data, the first and second diode saturation currents, I01 and I02, are determined and hence the base diffusion length and lifetime of minority carriers as functions of temperature. Lifetime increases exponentially with temperature and is explained by a Shockley-Read-Hall model with deep recombination levels 0.245 eV above the valence band. The I02 variation with temperature follows the Sah-Noyce-Shockley-Choo model except at low temperature where extra transitions raise the value above the predicted level. Reverse bias current at low voltage is a thermally assisted tunneling process via deep levels which are observed in base recombination at higher temperatures. The tunneling effects tend to become independent of temperature in the low-temperature region. These results demonstrate the ability to deduce basic parameters such as lifetime from simple measurements and show that back-surface fields offer no advantage at temperatures below 230 K. The analysis also explains the fall in lifetimes observed as the base conductivity increases, attributing it to native defects (perhaps carbon-oxygen-vacancy complexes) rather than the concentration of base dopant.

Anspaugh, B. E.↗
NASA NTRS2005

Constraint-based Temporal Reasoning with Preferences

Often we need to work in scenarios where events happen over time and preferences are associated to event distances and durations. Soft temporal constraints allow one to describe in a natural way problems arising in such scenarios. In general, solving soft temporal problems require exponential time in the worst case, but there are interesting subclasses of problems which are polynomially solvable. In this paper we identify one of such subclasses giving tractability results. Moreover, we describe two solvers for this class of soft temporal problems, and we show some experimental results. The random generator used to build the problems on which tests are performed is also described. We also compare the two solvers highlighting the tradeoff between performance and robustness. Sometimes, however, temporal local preferences are difficult to set, and it may be easier instead to associate preferences to some complete solutions of the problem. To model everything in a uniform way via local preferences only, and also to take advantage of the existing constraint solvers which exploit only local preferences, we show that machine learning techniques can be useful in this respect. In particular, we present a learning module based on a gradient descent technique which induces local temporal preferences from global ones. We also show the behavior of the learning module on randomly-generated examples.

Khatib, Lina↗
NASA NTRS1986

Statistical performance of the intensified charged coupled device

The first analysis of the statistical performance of the intensified charge coupled device (CCD) is presented. The Poisson statistics of the arriving photons, the exponential pulse height distribution of the microchannel plate in the image intensifier, and the effects of discrete digitization of the charge in reading the CCD are included. A range of intensifier gains from the high levels that permit pulse counting to the low levels that can be used with the accumulation of charge on the CCD is considered. As a tool for this analysis, an analytical model of the performance of the detector system that is generally useful for predicting many aspects of system performance under a wide range of conditions is developed. The statistical penalty for operating the intensified CCD in the charge accumulation mode is not severe, and in some applications it may be desirable to accept this penalty to gain the advantages of increased dynamic range through charge accumulation. This investigation has revealed no pathological aspects of the operation of the detector system.

Sandel, B. R.↗
NASA NTRS1990

Low NO(x) potential of gas turbine engines

The purpose is to correlate emission levels of gas turbine engines. The predictions of NO(x) emissions are based on a review of the literature of previous low NO(x) combustor programs and analytical chemical kinetic calculations. Concepts included in the literature review consisted of lean-premixed-prevaporized (LPP), rich burn/quick quench/lean burn (RQL), and direct injection. The NO(x) emissions were found to be an exponential function of adiabatic combustion temperature over a wide range of inlet temperatures, pressures and (lean) fuel-air ratios. A simple correlation of NO(x) formation with time was not found. The LPP and direct injection (using gaseous fuels) concepts have the lowest NO(x) emissions of the three concepts. The RQL data has higher values of NO(x) than the LPP concept, probably due to the stoichiometric temperatures and NO(x) production that occur during the quench step. Improvements in the quick quench step could reduce the NO(x) emissions to the LPP levels. The low NO(x) potential of LPP is offset by the operational disadvantages of its narrow stability limits and its susceptibility to autoignition/flashback. The Rich-Burn/Quick-Quench/Lean-Burn (RQL) and the direct injection concepts have the advantage of wider stability limits comparable to conventional combustors.

Tacina, Robert R.↗
NASA NTRS1995

The correlated k-distribution technique as applied to the AVHRR channels

Correlated k-distributions have been created to account for the molecular absorption found in the spectral ranges of the five Advanced Very High Resolution Radiometer (AVHRR) satellite channels. The production of the k-distributions was based upon an exponential-sum fitting of transmissions (ESFT) technique which was applied to reference line-by-line absorptance calculations. To account for the overlap of spectral features from different molecular species, the present routines made use of the multiplication transmissivity property which allows for considerable flexibility, especially when altering relative mixing ratios of the various molecular species. To determine the accuracy of the correlated k-distribution technique as compared to the line-by-line procedure, atmospheric flux and heating rate calculations were run for a wide variety of atmospheric conditions. For the atmospheric conditions taken into consideration, the correlated k-distribution technique has yielded results within about 0.5% for both the cases where the satellite spectral response functions were applied and where they were not. The correlated k-distribution's principal advantages is that it can be incorporated directly into multiple scattering routines that consider scattering as well as absorption by clouds and aerosol particles.

Kratz, David P.↗
NASA NTRS2008

Recent Insights into the Nature of Turbulence in the Solar Wind

During the past several years, studies of solar wind turbulence using data from Cluster and other spacecraft, and results from new numerical simulations, have revealed new aspects of solar wind turbulence. I will try to highlight some of that research. At the shortest length scales and highest frequencies, there is renewed interest in determining how the turbulence dissipates, e.g., whether by kinetic Alfven waves or whistler turbulence. Finding observational evidence for exponential damping of solar wind fluctuations has proven challenging. New studies using a combination of flux gate and search coil magnetometer data from Cluster have extended this search (in the spacecraft frame of reference) to more than 10 Hertz. New models and simulations are also being used to study the dissipation. A detailed study of fluctuations in the magnetosheath suggests that turbulent dissipation could be occurring at very thin current sheets as had been suggested by two-dimensional MHD simulations more than 20 years ago. Data from the four Cluster spacecraft, now at their maximum separation of 10,000 km provide new opportunities to investigate the symmetry properties, scale lengths, and the relative proportion of magnetic energy in parallel and perpendicular wave numbers of solar wind turbulence. By utilizing well-calibrated electron data, it has been possible to take advantage of the tetrahedral separation of Cluster in the solar wind near apogee to measure directly the compressibility and vorticity of the solar wind plasma.

Goldstein, Melvun L.↗
NASA NTRS2005

Developing a Framework for Effective Network Capacity Planning

As Internet traffic continues to grow exponentially, developing a clearer understanding of, and appropriately measuring, network's performance is becoming ever more critical. An important challenge faced by the Information Resources Directorate (IRD) at the Johnson Space Center in this context remains not only monitoring and maintaining a secure network, but also better understanding the capacity and future growth potential boundaries of its network. This requires capacity planning which involves modeling and simulating different network alternatives, and incorporating changes in design as technologies, components, configurations, and applications change, to determine optimal solutions in light of IRD's goals, objectives and strategies. My primary task this summer was to address this need. I evaluated network-modeling tools from OPNET Technologies Inc. and Compuware Corporation. I generated a baseline model for Building 45 using both tools by importing "real" topology/traffic information using IRD's various network management tools. I compared each tool against the other in terms of the advantages and disadvantages of both tools to accomplish IRD's goals. I also prepared step-by-step "how to design a baseline model" tutorial for both OPNET and Compuware products.

Yaprak, Ece↗
NASA NTRS2007

Generative Representations for Automated Design of Robots

A method of automated design of complex, modular robots involves an evolutionary process in which generative representations of designs are used. The term generative representations as used here signifies, loosely, representations that consist of or include algorithms, computer programs, and the like, wherein encoded designs can reuse elements of their encoding and thereby evolve toward greater complexity. Automated design of robots through synthetic evolutionary processes has already been demonstrated, but it is not clear whether genetically inspired search algorithms can yield designs that are sufficiently complex for practical engineering. The ultimate success of such algorithms as tools for automation of design depends on the scaling properties of representations of designs. A nongenerative representation (one in which each element of the encoded design is used at most once in translating to the design) scales linearly with the number of elements. Search algorithms that use nongenerative representations quickly become intractable (search times vary approximately exponentially with numbers of design elements), and thus are not amenable to scaling to complex designs. Generative representations are compact representations and were devised as means to circumvent the above-mentioned fundamental restriction on scalability. In the present method, a robot is defined by a compact programmatic form (its generative representation) and the evolutionary variation takes place on this form. The evolutionary process is an iterative one, wherein each cycle consists of the following steps: 1. Generative representations are generated in an evolutionary subprocess. 2. Each generative representation is a program that, when compiled, produces an assembly procedure. 3. In a computational simulation, a constructor executes an assembly procedure to generate a robot. 4. A physical-simulation program tests the performance of a simulated constructed robot, evaluating the performance according to a fitness criterion to yield a figure of merit that is fed back into the evolutionary subprocess of the next iteration. In comparison with prior approaches to automated evolutionary design of robots, the use of generative representations offers two advantages: First, a generative representation enables the reuse of components in regular and hierarchical ways and thereby serves a systematic means of creating more complex modules out of simpler ones. Second, the evolved generative representation may capture intrinsic properties of the design problem, so that variations in the representations move through the design space more effectively than do equivalent variations in a nongenerative representation. This method has been demonstrated by using it to design some robots that move, variously, by walking, rolling, or sliding. Some of the robots were built (see figure). Although these robots are very simple, in comparison with robots designed by humans, their structures are more regular, modular, hierarchical, and complex than are those of evolved designs of comparable functionality synthesized by use of nongenerative representations.

Homby, Gregory S.↗
NASA NTRS2014

Phase Transitions in Planning Problems: Design and Analysis of Parameterized Families of Hard Planning Problems

There are two common ways to evaluate algorithms: performance on benchmark problems derived from real applications and analysis of performance on parametrized families of problems. The two approaches complement each other, each having its advantages and disadvantages. The planning community has concentrated on the first approach, with few ways of generating parametrized families of hard problems known prior to this work. Our group's main interest is in comparing approaches to solving planning problems using a novel type of computational device - a quantum annealer - to existing state-of-the-art planning algorithms. Because only small-scale quantum annealers are available, we must compare on small problem sizes. Small problems are primarily useful for comparison only if they are instances of parametrized families of problems for which scaling analysis can be done. In this technical report, we discuss our approach to the generation of hard planning problems from classes of well-studied NP-complete problems that map naturally to planning problems or to aspects of planning problems that many practical planning problems share. These problem classes exhibit a phase transition between easy-to-solve and easy-to-show-unsolvable planning problems. The parametrized families of hard planning problems lie at the phase transition. The exponential scaling of hardness with problem size is apparent in these families even at very small problem sizes, thus enabling us to characterize even very small problems as hard. The families we developed will prove generally useful to the planning community in analyzing the performance of planning algorithms, providing a complementary approach to existing evaluation methods. We illustrate the hardness of these problems and their scaling with results on four state-of-the-art planners, observing significant differences between these planners on these problem families. Finally, we describe two general, and quite different, mappings of planning problems to QUBOs, the form of input required for a quantum annealing machine such as the D-Wave II.

Problems↗
NASA NTRS2010

Hidden Statistics Approach to Quantum Simulations

Recent advances in quantum information theory have inspired an explosion of interest in new quantum algorithms for solving hard computational (quantum and non-quantum) problems. The basic principle of quantum computation is that the quantum properties can be used to represent structure data, and that quantum mechanisms can be devised and built to perform operations with this data. Three basic non-classical properties of quantum mechanics superposition, entanglement, and direct-product decomposability were main reasons for optimism about capabilities of quantum computers that promised simultaneous processing of large massifs of highly correlated data. Unfortunately, these advantages of quantum mechanics came with a high price. One major problem is keeping the components of the computer in a coherent state, as the slightest interaction with the external world would cause the system to decohere. That is why the hardware implementation of a quantum computer is still unsolved. The basic idea of this work is to create a new kind of dynamical system that would preserve the main three properties of quantum physics superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. In other words, such a system would reinforce the advantages and minimize limitations of both quantum and classical aspects. Based upon a concept of hidden statistics, a new kind of dynamical system for simulation of Schroedinger equation is proposed. The system represents a modified Madelung version of Schroedinger equation. It preserves superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. Such an optimal combination of characteristics is a perfect match for simulating quantum systems. The model includes a transitional component of quantum potential (that has been overlooked in previous treatment of the Madelung equation). The role of the transitional potential is to provide a jump from a deterministic state to a random state with prescribed probability density. This jump is triggered by blowup instability due to violation of Lipschitz condition generated by the quantum potential. As a result, the dynamics attains quantum properties on a classical scale. The model can be implemented physically as an analog VLSI-based (very-large-scale integration-based) computer, or numerically on a digital computer. This work opens a way of developing fundamentally new algorithms for quantum simulations of exponentially complex problems that expand NASA capabilities in conducting space activities. It has been illustrated that the complexity of simulations of particle interaction can be reduced from an exponential one to a polynomial one.

Zak, Michail↗
NASA NTRS2008

Spiking Neurons for Analysis of Patterns

Artificial neural networks comprising spiking neurons of a novel type have been conceived as improved pattern-analysis and pattern-recognition computational systems. These neurons are represented by a mathematical model denoted the state-variable model (SVM), which among other things, exploits a computational parallelism inherent in spiking-neuron geometry. Networks of SVM neurons offer advantages of speed and computational efficiency, relative to traditional artificial neural networks. The SVM also overcomes some of the limitations of prior spiking-neuron models. There are numerous potential pattern-recognition, tracking, and data-reduction (data preprocessing) applications for these SVM neural networks on Earth and in exploration of remote planets. Spiking neurons imitate biological neurons more closely than do the neurons of traditional artificial neural networks. A spiking neuron includes a central cell body (soma) surrounded by a tree-like interconnection network (dendrites). Spiking neurons are so named because they generate trains of output pulses (spikes) in response to inputs received from sensors or from other neurons. They gain their speed advantage over traditional neural networks by using the timing of individual spikes for computation, whereas traditional artificial neurons use averages of activity levels over time. Moreover, spiking neurons use the delays inherent in dendritic processing in order to efficiently encode the information content of incoming signals. Because traditional artificial neurons fail to capture this encoding, they have less processing capability, and so it is necessary to use more gates when implementing traditional artificial neurons in electronic circuitry. Such higher-order functions as dynamic tasking are effected by use of pools (collections) of spiking neurons interconnected by spike-transmitting fibers. The SVM includes adaptive thresholds and submodels of transport of ions (in imitation of such transport in biological neurons). These features enable the neurons to adapt their responses to high-rate inputs from sensors, and to adapt their firing thresholds to mitigate noise or effects of potential sensor failure. The mathematical derivation of the SVM starts from a prior model, known in the art as the point soma model, which captures all of the salient properties of neuronal response while keeping the computational cost low. The point-soma latency time is modified to be an exponentially decaying function of the strength of the applied potential. Choosing computational efficiency over biological fidelity, the dendrites surrounding a neuron are represented by simplified compartmental submodels and there are no dendritic spines. Updates to the dendritic potential, calcium-ion concentrations and conductances, and potassium-ion conductances are done by use of equations similar to those of the point soma. Diffusion processes in dendrites are modeled by averaging among nearest-neighbor compartments. Inputs to each of the dendritic compartments come from sensors. Alternatively or in addition, when an affected neuron is part of a pool, inputs can come from other spiking neurons. At present, SVM neural networks are implemented by computational simulation, using algorithms that encode the SVM and its submodels. However, it should be possible to implement these neural networks in hardware: The differential equations for the dendritic and cellular processes in the SVM model of spiking neurons map to equivalent circuits that can be implemented directly in analog very-large-scale integrated (VLSI) circuits.

Huntsberger, Terrance↗