Search NASA⌕ Search

SEARCH · Search NASA

Results for “satisfiability problems”

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 793 records · Page 44

Failure detection and identification

Using the geometric concept of an unobservability subspace, a solution is given to the problem of detecting and identifying control system component failures in linear, time-invariant systems. Conditions are developed for the existence of a causal, linear, time-invariant processor that can detect and uniquely identify a component failure, first for the case where components can fail simultaneously, and then for the case where they fail only one at a time. Explicit design algorithms are provided when these conditions are satisfied. In addition to time-domain solvability conditions, frequency-domain interpretations of the results are given, and connections are drawn with results already available in the literature.

Massoumnia, Mohammad-Ali↗

Superresolution via sparsity constraints

The problem of recovering a measure mu supported on a lattice of span Delta is considered under the condition that measurements are only available concerning the Fourier Transform at frequencies of Omega or less. If Omega is much smaller than the Nyquist frequency pi/Delta and the measurements are noisy, then stable recovery of mu is generally impossible. It is shown here that if, in addition, it is known that mu satisfies certain sparsity constraints, then stable recovery is possible. This finding validates practical efforts in spectroscopy, seismic prospecting, and astronomy to provide superresolution by imposing support limitations in reconstruction.

Donoho, David L.↗

A Comparison of Structurally Connected and Multiple Spacecraft Interferometers

Structurally connected and multiple spacecraft interferometers are compared in an attempt to establish the maximum baseline (referred to as the "cross-over baseline") for which it is preferable to operate a single-structure interferometer in space rather than an interferometer composed of numerous, smaller spacecraft. This comparison is made using the total launched mass of each configuration as the comparison metric. A framework of study within which structurally connected and multiple spacecraft interferometers can be compared is presented in block diagram form. This methodology is then applied to twenty-two different combinations of trade space parameters to investigate the effects of different orbits, orientations, truss materials, propellants, attitude control actuators, onboard disturbance sources, and performance requirements on the cross-over baseline. Rotating interferometers and the potential advantages of adding active structural control to the connected truss of the structurally connected interferometer are also examined. The minimum mass design of the structurally connected interferometer that meets all performance-requirements and satisfies all imposed constraints is determined as a function of baseline. This minimum mass design is then compared to the design of the multiple spacecraft interferometer. It is discovered that the design of the minimum mass structurally connected interferometer that meets all performance requirements and constraints in solar orbit is limited by the minimum allowable aspect ratio, areal density, and gage of the struts. In the formulation of the problem used in this study, there is no advantage to adding active structural control to the truss for interferometers in solar orbit. The cross-over baseline for missions of practical duration (ranging from one week to thirty years) in solar orbit is approximately 400 m for non-rotating interferometers and 650 m for rotating interferometers.

Surka, Derek M.↗

A Novel Approach for Adaptive Signal Processing

Adaptive linear predictors have been used extensively in practice in a wide variety of forms. In the main, their theoretical development is based upon the assumption of stationarity of the signals involved, particularly with respect to the second order statistics. On this basis, the well-known normal equations can be formulated. If high- order statistical stationarity is assumed, then the equivalent normal equations involve high-order signal moments. In either case, the cross moments (second or higher) are needed. This renders the adaptive prediction procedure non-blind. A novel procedure for blind adaptive prediction has been proposed and considerable implementation has been made in our contributions in the past year. The approach is based upon a suitable interpretation of blind equalization methods that satisfy the constant modulus property and offers significant deviations from the standard prediction methods. These blind adaptive algorithms are derived by formulating Lagrange equivalents from mechanisms of constrained optimization. In this report, other new update algorithms are derived from the fundamental concepts of advanced system identification to carry out the proposed blind adaptive prediction. The results of the work can be extended to a number of control-related problems, such as disturbance identification. The basic principles are outlined in this report and differences from other existing methods are discussed. The applications implemented are speech processing, such as coding and synthesis. Simulations are included to verify the novel modelling method.

Chen, Ya-Chin↗

A nonlinear finite-element analysis of wings in steady incompressible flows with wake roll-up

The problem of lifting surfaces and complex aircraft configurations in steady incompressible flow is considered. For lifting surfaces the problem is formulated in terms of an integral equation relating the potential discontinuity on wing and wake to the normal derivative of the potential on the lifting surface. For complex configurations the problem is formulated in terms of an integral equation relating the potential to its normal derivative on the surface of the aircraft. The integral equation is approximated by a system of linear algebraic equations obtained by dividing the surfaces into small quadrilateral elements and by assuming the potential (or the potential discontinuity) and its normal derivative to be constant within each element. The wake geometry is obtained by iteration by satisfying the condition that the velocity be tangent to the surface of the wake and that the potential discontinuity be constant along the streamlines.

Suciu, E. O.↗

On Global Optimal Sailplane Flight Strategy

The derivation and interpretation of the necessary conditions that a sailplane cross-country flight has to satisfy to achieve the maximum global flight speed is considered. Simple rules are obtained for two specific meteorological models. The first one uses concentrated lifts of various strengths and unequal distance. The second one takes into account finite, nonuniform space amplitudes for the lifts and allows, therefore, for dolphin style flight. In both models, altitude constraints consisting of upper and lower limits are shown to be essential to model realistic problems. Numerical examples illustrate the difference with existing techniques based on local optimality conditions.

Sander, G. J.↗

Future payload isolation and pointing system technology

An overview is presented of the technology associated with the precision pointing of future spaceborne science instruments. High-level architectural options for dealing with the precision pointing problem are considered, and a representative NASA mission set for the 1990s is given. Pointing accuracy and stability requirements are extracted from the mission requirements, with the sub-0.1 arcsec stability requirement emerging as the primary driver. The state of the art of current technology is assessed, including an evaluation of gimbal systems, suspension systems, and actuator and sensor component technology. Areas where the technology needs to be pushed to satisfy future requirements are identified, and some promising design options are proposed.

Laskin, R. A.↗

Avionics Tether Operations Control

The activities described in this Final Report were authorized and performed under Purchase Order Number H32835D, issued as part of NASA contract number NAS8-00114. The period of performance of this PO was from March 1 to September 30, 2001. The primary work activity was the continued development and updating of the tether dynamic simulation tools GTOSS (Generalized Tethered Object System Simulation) and TSSIM (Tethered Satellite System) and use of these and other tools in the analysis of various tether dynamics problems. Several updated versions of GTOSS were delivered during the period of performance by the author of the simulation, Lang Associates' David Lang. These updates had mainly to do with updated documentation and an updated coordinate system definition to the J2000 standards. This Final Report is organized by the months in which the activities described were performed. The following sections review the Statement of Work (SOW) and activities performed to satisfy it.

Glaese, John R.↗

ECLEIRS: Exact conservation law embedded identification of reduced states for parameterized nonlinear conservation laws from sparse and noisy data

Multi-query applications such as parameter estimation, uncertainty quantification and design optimization for parameterized partial differential equation (PDE) systems are expensive. While reduced/latent state dynamics approaches for parameterized PDEs offer a viable alternative, these approaches rely on high-quality data and struggle with highly sparse spatiotemporal noisy measurements typically obtained from experiments. Furthermore, there is no guarantee that these models satisfy governing physical conservation laws. In this article, we propose a reduced state dynamics approach, referred to as ECLEIRS, that embeds exact conservation in the solution and flux representation by utilizing a space-time divergence-free neural network formulation. We compare ECLEIRS with other reduced state dynamics approaches, those that do not enforce any physical constraints and those with physics-informed loss functions, for three shock-propagation problems: 1-D advection, 1-D Burgers and 2-D Euler equations. In conclusion, the numerical experiments conducted in this study demonstrate that ECLEIRS provides the most accurate prediction of dynamics for unseen parameters even in the presence of highly sparse and noisy data.

97 MATHEMATICS AND COMPUTING↗

Planetary entry aerothermodynamics. II - Computational analyses and flight experience

Limitations concerning the possibility to simulate all the significant flow and thermal phenomena occurring during the entry of a space vehicle into a planetary atmosphere make it necessary to rely on computational analyses to obtain the required data for the design of the spacecraft needed for the NASA missions planned for the next two decades. 'Benchmark' computer programs concerned with complete, detailed, and accurate computational solutions of entry problems are considered along with programs representing engineering approximations for cases in which the accuracy provided by the benchmark programs is not needed. The information obtainable by computational analysis has to be supplemented by actual flight experience in order to meet the goals of the NASA entry-technology program. The individual space missions planned for the coming years are examined together with the possibilities for obtaining the data needed to satisfy the entry requirements in each case.

Olstad, W. B.↗

Multi-point inverse design of an infinite cascade of airfoils

This paper describes a practical method for the design of an infinite cascade in incompressible flow. The method is based on conformal mapping, and as a result it allows for multi-point design. The cascade blade to be determined is divided into a desired number of segments. Over each segment, the velocity distribution is prescribed together with an inlet or outlet flow angle at which this velocity distribution is to be achieved. In this way multi-point design requirements can be met. It is necessary to satisfy several conditions that arise to guarantee compatibility with the inlet and outlet flow as well as closure of the cascade blade. Satisfaction of these conditions does not necessarily result in a cascade with all the desired characteristics. For example, the cascade blades may be bulbous or crossed. Through Newton iteration, however, the desired characteristics may be prescribed by allowing for the adjustment of the design parameters that define the mathematical problem through conformal mapping. Several examples will be illustrated to demonstrate the capability of the method. It will be shown that the method is limited to the design of cascades with solidities of up to one.

Selig, Michael S.↗

Use of source distributions for evaluating theoretical aerodynamics of thin finite wings at supersonic speeds

A series of publications on the source-distribution methods for evaluating the aerodynamics of thin wings at supersonic speeds is summarized, extended, and unified. Included in the first part are the deviations of: (a) the linearized partial-differential equation for unsteady flow at a substantially constant Mach number. b) The source-distribution solution for the perturbation-velocity potential that satisfies the boundary conditions of tangential flow at the surface and in the plane of the wing; and (c) the integral equation for determining the strength and the location of sources to describe the interaction effects (as represented by upwash) of the bottom and top wing surfaces through the region between the finite wing boundary and the foremost Mach wave. The second part deals with steady-state thin-wing problems. The third part of the report approximates the integral equation for unsteady upwash and includes a solution of approximate equation. Expressions are then derived to evaluate the load distributions for time-dependent finite-wing motions.

Evvard, John C↗

Primal interface debonding formulation for finite strain isotropic plasticity

In this work, a framework is developed for modeling ductile damage of nonlinear materials whose plastic deformation is characterized using rate independent classical plasticity. This method relies on the assumption that the free energy can be decomposed into elastic, plastic and damage parts. A thermodynamically consistent method is derived which satisfies the second law of thermodynamics in the Clausius–Duhem inequality form. The dissipation associated with plasticity takes place in the domain only, while damage dissipation is localized to the interface. The method is developed using Variational Multiscale ideas to obtain definitions of the interface fluxes within a primal formulation analogous to the Discontinuous Galerkin method, which ensures weakly vanishing interface gap prior to reaching a damage initiation criterion. The local nonlinear problem to calculate both plastic deformation gradient and damage variable follows an incremental approach similar to classical plasticity return mapping algorithm. This elastoplastic damage formulation is developed for material undergoing finite strain, and it naturally accommodates a trapezoidal traction separation law (TSL) whose shape can be varied to model either ductile interface behavior or brittle interface behavior. The formulation's performance is assessed through modeling a patch test and a compact tension specimen.

42 ENGINEERING↗

A multilevel approach for minimum weight structural design including local and system buckling constraints

A rational multilevel approach for minimum weight structural design of truss and wing structures including local and system buckling constraints is presented. Overall proportioning of the structure is achieved at the system level subject to strength, displacement and system buckling constraints, while the detailed component designs are carried out separately at the component level satisfying local buckling constraints. Total structural weight is taken to be the objective function at the system level while employing the change in the equivalent system stiffness of the component as the component level objective function. Finite element analysis is used to predict static response while system buckling behavior is handled by incorporating a geometric stiffness matrix capability. Buckling load factors and the corresponding mode shapes are obtained by solving the eigenvalue problem associated with the assembled elastic stiffness and geometric stiffness matrices for the structural system. At the component level various local buckling failure modes are guarded against using semi-empirical formulas. Mathematical programming techniques are employed at both the system and component level.

Schmit, L. A., Jr.↗

Reconstruction of atmospheric pollutant concentrations from remote sensing data - An application of distributed parameter observer theory

The reconstruction of a concentration distribution from spatially averaged and noise-corrupted data is a central problem in processing atmospheric remote sensing data. Distributed parameter observer theory is used to develop reconstructibility conditions for distributed parameter systems having measurements typical of those in remote sensing. The relation of the reconstructibility condition to the stability of the distributed parameter observer is demonstrated. The theory is applied to a variety of remote sensing situations, and it is found that those in which concentrations are measured as a function of altitude satisfy the conditions of distributed state reconstructibility.

Koda, M.↗

Consistent boundary conditions for Reduced Navier-Stokes (RNS) scheme applied to 3-dimensional internal viscous flows

A consistent and efficient set of boundary conditions is developed for the multi-sweep space-marching pressure-elliptic Reduced Navier-Stokes (RNS) scheme as applied for 3-D internal viscous flow problems. No-slip boundary conditions are directly imposed on the solid walls. There is no iteration procedure required in the cross plane to ensure mass conservation across each marching plane. The finite difference equations forming the coefficient matrix are ordered such that the surface normal velocity is specified on all the solid walls; unlike external flows, a pressure boundary condition in the cross plane is not required. Since continuity is directly satisfied at all points in the flow domain, the first order momentum equations can be solved directly for the pressure without the need for a Poisson pressure correction equation. The procedure developed herein can also be applied with periodic boundary conditions. The analysis is given for general compressible flows. Incompressible flow solutions are obtained, for straight and curved ducts of square cross section, to validate the procedure. These solutions are used to demonstrate the applicability of the RNS scheme, with the improved boundary conditions for internal flows with strong interaction, as would be encountered in ducts and turbomachinery geometries.

Reddy, D. R.↗

Optimizing fluvial flood mitigation strategies: A multi-objective approach for cost-effective and socially-aware infrastructure feasibility analysis

Effective levee planning must balance capital cost, risk reduction, and community priorities. These objectives are rarely optimized together. This study presents a feasibility phase, simulationin-the-loop framework that couples terrain-based flood modeling with a socially aware multiobjective optimizer. Flood risk is measured as Expected Annual Exposed Population (EAEP), obtained by integrating exposure over Annual Exceedance Probability (AEP) nodes, mirroring the Hydrologic Engineering Center's Flood Damage Reduction Analysis (HEC-FDA) expected-annual formulation but with people rather than dollars. Exposure per scenario is computed by overlaying binary inundation masks with a population surface at the tract level. Distributional fairness is encoded through a Group Benefit Share (GBS) constraint that requires high-SVI tracts to receive at least a baseline share of annualized benefits. Capital cost is represented by a height-dependent unit-cost model suitable for screening. This study addresses the two-objective problem, minimize cost and expected annual exposure subject to the GBS constraint, using Non-Dominated Sorting Genetic Algorithm II (NSGA-II) and leveraging Pareto front for feasibility phase decision making. Implemented with terrain-based flood modeling, GeoFlood, for rapid scenario evaluation, the framework is demonstrated in Southeast Texas. The results reveal clear trade-offs among cost, risk, and social benefits and identify non-dominated levee height configurations that satisfy the benefit-share floor. The contributions are a scalable decision support method that operationalizes expected annual population-based risk, embeds enforceable benefit-sharing guarantees, and uses lightweight simulation to explore large design spaces before higher fidelity design stages.

Flood mitigation↗

Numerical modelling of instantaneous plate tectonics

Assuming lithospheric plates to be rigid, 68 spreading rates, 62 fracture zones trends, and 106 earthquake slip vectors are systematically inverted to obtain a self-consistent model of instantaneous relative motions for eleven major plates. The inverse problem is linearized and solved iteratively by a maximum-likelihood procedure. Because the uncertainties in the data are small, Gaussian statistics are shown to be adequate. The use of a linear theory permits (1) the calculation of the uncertainties in the various angular velocity vectors caused by uncertainties in the data, and (2) quantitative examination of the distribution of information within the data set. The existence of a self-consistent model satisfying all the data is strong justification of the rigid plate assumption. Slow movement between North and South America is shown to be resolvable.

Minster, J. B.↗