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.

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Control system design for flexible structures using data models

The dynamics and control of flexible aerospace structures exercises many of the engineering disciplines. In recent years there has been considerable research in the developing and tailoring of control system design techniques for these structures. This problem involves designing a control system for a multi-input, multi-output (MIMO) system that satisfies various performance criteria, such as vibration suppression, disturbance and noise rejection, attitude control and slewing control. Considerable progress has been made and demonstrated in control system design techniques for these structures. The key to designing control systems for these structures that meet stringent performance requirements is an accurate model. It has become apparent that theoretically and finite-element generated models do not provide the needed accuracy; almost all successful demonstrations of control system design techniques have involved using test results for fine-tuning a model or for extracting a model using system ID techniques. This paper describes past and ongoing efforts at Ohio University and NASA MSFC to design controllers using 'data models.' The basic philosophy of this approach is to start with a stabilizing controller and frequency response data that describes the plant; then, iteratively vary the free parameters of the controller so that performance measures become closer to satisfying design specifications. The frequency response data can be either experimentally derived or analytically derived. One 'design-with-data' algorithm presented in this paper is called the Compensator Improvement Program (CIP). The current CIP designs controllers for MIMO systems so that classical gain, phase, and attenuation margins are achieved. The center-piece of the CIP algorithm is the constraint improvement technique which is used to calculate a parameter change vector that guarantees an improvement in all unsatisfied, feasible performance metrics from iteration to iteration. The paper also presents a recently demonstrated CIP-type algorithm, called the Model and Data Oriented Computer-Aided Design System (MADCADS), developed for achieving H(sub infinity) type design specifications using data models. Control system design for the NASA/MSFC Single Structure Control Facility are demonstrated for both CIP and MADCADS. Advantages of design-with-data algorithms over techniques that require analytical plant models are also presented.

Irwin, R. Dennis↗

NASA NTRS1950

The Boundary Layers in Fluids with Little Friction

The vortices forming in flowing water behind solid bodies are not represented correctly by the solution of the potential theory nor by Helmholtz's jets. Potential theory is unable to satisfy the condition that the water adheres at the wetted bodies, and its solutions of the fundamental hydrodynamic equations are at variance with the observation that the flow separates from the body at a certain point and sends forth a highly turbulent boundary layer into the free flow. Helmholtz's theory attempts to imitate the latter effect in such a way that it joins two potential flows, jet and still water, nonanalytical along a stream curve. The admissibility of this method is based on the fact that, at zero pressure, which is to prevail at the cited stream curve, the connection of the fluid, and with it the effect of adjacent parts on each other, is canceled. In reality, however, the pressure at these boundaries is definitely not zero, but can even be varied arbitrarily. Besides, Helmholtz's theory with its potential flows does not satisfy the condition of adherence nor explain the origin of the vortices, for in all of these problems, the friction must be taken into account on principle, according to the vortex theorem.

Blasius, H.↗

NASA NTRS2014

Fuzzy Model-based Pitch Stabilization and Wing Vibration Suppression of Flexible Wing Aircraft.

This paper presents a fuzzy nonlinear controller to regulate the longitudinal dynamics of an aircraft and suppress the bending and torsional vibrations of its flexible wings. The fuzzy controller utilizes full-state feedback with input constraint. First, the Takagi-Sugeno fuzzy linear model is developed which approximates the coupled aeroelastic aircraft model. Then, based on the fuzzy linear model, a fuzzy controller is developed to utilize a full-state feedback and stabilize the system while it satisfies the control input constraint. Linear matrix inequality (LMI) techniques are employed to solve the fuzzy control problem. Finally, the performance of the proposed controller is demonstrated on the NASA Generic Transport Model (GTM).

Stability and Control↗

NASA NTRS2013

MESSENGER at Mercury: Early Orbital Operations

The MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft, launched in August 2004 under NASA's Discovery Program, was inserted into orbit about the planet Mercury in March 2011. MESSENGER's three flybys of Mercury in 2008-2009 marked the first spacecraft visits to the innermost planet since the Mariner 10 flybys in 1974-1975. The unprecedented orbital operations are yielding new insights into the nature and evolution of Mercury. The scientific questions that frame the MESSENGER mission led to the mission measurement objectives to be achieved by the seven payload instruments and the radio science experiment. Interweaving the full set of required orbital observations in a manner that maximizes the opportunity to satisfy all mission objectives and yet meet stringent spacecraft pointing and thermal constraints was a complex optimization problem that was solved with a software tool that simulates science observations and tracks progress toward meeting each objective. The final orbital observation plan, the outcome of that optimization process, meets all mission objectives. MESSENGER's Mercury Dual Imaging System is acquiring a global monochromatic image mosaic at better than 90% coverage and at least 250 m average resolution, a global color image mosaic at better than 90% coverage and at least 1 km average resolution, and global stereo imaging at better than 80% coverage and at least 250 m average resolution. Higher-resolution images are also being acquired of targeted areas. The elemental remote sensing instruments, including the Gamma-Ray and Neutron Spectrometer and the X-Ray Spectrometer, are being operated nearly continuously and will establish the average surface abundances of most major elements. The Visible and Infrared Spectrograph channel of MESSENGER's Mercury Atmospheric and Surface Composition Spectrometer is acquiring a global map of spectral reflectance from 300 to 1450 nm wavelength at a range of incidence and emission angles. Targeted areas have been selected for spectral coverage into the ultraviolet with the Ultraviolet and Visible Spectrometer (UVVS). MESSENGER's Mercury Laser Altimeter is acquiring topographic profiles when the slant range to Mercury's surface is less than 1800 km, encompassing latitudes from 20 deg. S to the north pole. Topography over the remainder of the southern hemisphere will be derived from stereo imaging, radio occultations, and limb profiles. MESSENGER's radio science experiment is determining Mercury's gravity field from Doppler signals acquired during frequent downlinks. MESSENGER's Magnetometer is measuring the vector magnetic field both within Mercury's magnetosphere and in Mercury's solar wind environment at an instrument sampling rate of up to 20 samples/s. The UVVS is determining the three-dimensional, time-dependent distribution of Mercury's exospheric neutral and ionic species via their emission lines. During each spacecraft orbit, the Energetic Particle Spectrometer measures energetic electrons and ions, and the Fast Imaging Plasma Spectrometer measures the energies and mass per charge of thermal plasma components, both within Mercury's magnetosphere and in Mercury's solar-wind environment. The primary mission observation sequence will continue for one Earth year, until March 2012. An extended mission, currently under discussion with NASA, would add a second year of orbital observations targeting a set of focused follow-on questions that build on observations to date and take advantage of the more active Sun expected during 2012-2013. MESSENGER's total primary mission cost, projected at $446 M in real-year dollars, is comparable to that of Mariner 10 after adjustment for inflation.

mercury↗

NASA NTRS2012

MESSENGER at Mercury: Early Orbital Operations

The MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft, launched in August 2004 under NASA's Discovery Program, was inserted into orbit about the planet Mercury in March 2011. MESSENGER's three flybys of Mercury in 2008-2009 marked the first spacecraft visits to the innermost planet since the Mariner 10 flybys in 1974-1975. The unprecedented orbital operations are yielding new insights into the nature and evolution of Mercury. The scientific questions that frame the MESSENGER mission led to the mission measurement objectives to be achieved by the seven payload instruments and the radio science experiment. Interweaving the full set of required orbital observations in a manner that maximizes the opportunity to satisfy all mission objectives and yet meet stringent spacecraft pointing and thermal constraints was a complex optimization problem that was solved with a software tool that simulates science observations and tracks progress toward meeting each objective. The final orbital observation plan, the outcome of that optimization process, meets all mission objectives. MESSENGER's Mercury Dual Imaging System is acquiring a global monochromatic image mosaic at better than 90%coverage and at least 250 m average resolution, a global color image mosaic at better than 90%coverage and at least 1 km average resolution, and global stereo imaging at better than 80%coverage and at least 250 m average resolution. Higher-resolution images are also being acquired of targeted areas. The elemental remote sensing instruments, including the Gamma-Ray and Neutron Spectrometer and the X-Ray Spectrometer, are being operated nearly continuously and will establish the average surface abundances of most major elements. The Visible and Infrared Spectrograph channel of MESSENGER's Mercury Atmospheric and Surface Composition Spectrometer is acquiring a global map of spectral reflectance from 300 to 1450 nm wavelength at a range of incidence and emission angles. Targeted areas have been selected for spectral coverage into the ultraviolet with the Ultraviolet and Visible Spectrometer (UVVS). MESSENGER's Mercury Laser Altimeter is acquiring topographic profiles when the slant range to Mercury's surface is less than 1800 km, encompassing latitudes from 201S to the north pole. Topography over the remainder of the southern hemisphere will be derived from stereo imaging, radio occultations, and limb profiles. MESSENGER's radio science experiment is determining Mercury's gravity field from Doppler signals acquired during frequent downlinks. MESSENGER's Magnetometer is measuring the vector magnetic field both within Mercury's magnetosphere and in Mercury's solar wind environment at an instrument sampling rate of up to 20 samples/s. The UVVS is determining the three-dimensional, time-dependent distribution of Mercury's exospheric neutral and ionic species via their emission lines. During each spacecraft orbit, the Energetic Particle Spectrometer measures energetic electrons and ions, and the Fast Imaging Plasma Spectrometer measures the energies and mass per charge of thermal plasma components, both within Mercury's magnetosphere and in Mercury's solar-wind environment. The primary mission observation sequence will continue for one Earth year, until March 2012. An extended mission, currently under discussion with NASA, would add a second year of orbital observations targeting a set of focused follow-on questions that build on observations to date and take advantage of the more active Sun expected during 2012-2013. MESSENGER's total primary mission cost, projected at $446 M in real-year dollars, is comparable to that of Mariner 10 after adjustment for inflation.

Planetary Science↗

NASA NTRS1984

Development of a dynamically adaptive grid method for multidimensional problems

An approach to solution adaptive grid generation for use with finite difference techniques, previously demonstrated on model problems in one space dimension, has been extended to multidimensional problems. The method is based on the popular elliptic steady grid generators, but is 'dynamically' adaptive in the sense that a grid is maintained at all times satisfying the steady grid law driven by a solution-dependent source term. Testing has been carried out on Burgers' equation in one and two space dimensions. Results appear encouraging both for inviscid wave propagation cases and viscous boundary layer cases, suggesting that application to practical flow problems is now possible. In the course of the work, obstacles relating to grid correction, smoothing of the solution, and elliptic equation solvers have been largely overcome. Concern remains, however, about grid skewness, boundary layer resolution and the need for implicit integration methods. Also, the method in 3-D is expected to be very demanding of computer resources.

Holcomb, J. E.↗

NASA NTRS2012

Creation of an Upper Stage Trajectory Capability Boundary to Enable Booster System Trade Space Exploration

The problem of trajectory optimization is important in all space missions. The solution of this problem enables one to specify the optimum thrust steering program which should be followed to achieve a specified mission objective, simultaneously satisfying the constraints.1 It is well known that whether or not the ascent trajectory is optimal can have a significant impact on propellant usage for a given payload, or on payload weight for the same gross vehicle weight.2 Consequently, ascent guidance commands are usually optimized in some fashion. Multi-stage vehicles add complexity to this analysis process as changes in vehicle properties in one stage propagate to the other stages through gear ratios and changes in the optimal trajectory. These effects can cause an increase in analysis time as more variables are added and convergence of the optimizer to system closure requires more analysis iterations. In this paper, an approach to simplifying this multi-stage problem through the creation of an upper stage capability boundary is presented. This work was completed as part of a larger study focused on trade space exploration for the advanced booster system that will eventually form a part of NASA s new Space Launch System.3 The approach developed leverages Design of Experiments and Surrogate Modeling4 techniques to create a predictive model of the SLS upper stage performance. The design of the SLS core stages is considered fixed for the purposes of this study, which results in trajectory parameters such as staging conditions being the only variables relevant to the upper stage. Through the creation of a surrogate model, which takes staging conditions as inputs and predicts the payload mass delivered by the SLS upper stage to a reference orbit as the response, it is possible to identify a "surface" of staging conditions which all satisfy the SLS requirement of placing 130 metric tons into low-Earth orbit (LEO).3 This identified surface represents the 130 metric ton capability boundary for the upper stage, such that if the combined first stage and boosters can achieve any one staging point on that surface, then the design is identified as feasible. With the surrogate model created, design and analysis of advanced booster concepts is streamlined, as optimization of the upper stage trajectory is no longer required in every design loop.

Walsh, Ptrick↗

DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

Triple galaxies and a hidden mass problem

The authors consider a homogeneous sample of 84 triple systems of galaxies with components brighter than m = 15.7, located in the northern sky and satisfying an isolation criterion with respect to neighboring galaxies in projection. The distributions of basic dynamical parameters for triplets have median values as follows: radial velocity dispersion 133 km/s, mean harmonic radius 63 kpc, absolute magnitude of galaxies M sub B equals -20.38, crossing time tau = 0.04 H(sup minus 1). For different ways of estimation the median mass-to-luminosity ratio is (20 - 30). A comparison of the last value with the ones for single and binary galaxies shows the presence of a virial mass excess for triplets by a factor 4. The mass-to-luminosity ratio is practically uncorrelated with linear size of triplets or with morphological types of their components. We note that a significant part of the virial excess may be explained by the presence of nonisolated triple configurations in the sample, which are produced by debris of more populous groups of galaxies.

Karachentsev, I. D.↗

DOE OSTI2025

Multi-Resolution UAV Path Replanning for Inspection of Tailings Dams

Autonomous inspection of large and complex structures with a commercial unmanned aerial vehicle (UAV) is a challenging problem that has been addressed in recent years. In this paper, we address the global motion planning problem of creating autonomous inspection missions for UAVs considering photogrammetry constraints. We focus on the inspection of large tailings dams, which are dam structures used to store waste byproducts of mining. Our method uses a prior sparse point cloud of the dam to generate a voxel grid, where paths satisfying photogrammetry constraints are tested for collisions. We then apply the A* algorithm as a local planner to avoid obstacles within the global mission. Moreover, we address the problem of changing routes online by using octree-based multi-resolution grids for efficient and fast pathfinding. Our results, obtained using tridimensional maps of an actual coal mine tailings dam, show that using octrees for multi-resolution motion planning is faster than using a fixed voxel grid in online missions while inspecting large structures.

42 ENGINEERING↗

NASA NTRS1973

On electromagnetic field problems in inhomogeneous media

Analysis of electromagnetic fields in inhomogeneous media is of practical interest in general scattering and propagation problems and in the study of lenses. For certain types of inhomogeneities, the fields may be represented in terms of two scalars. In a general orthogonal coordinate system, these potentials satisfy second order differential equations. Exact solutions of these equations are known only for a few particular cases and in general, an approximate or numerical technique must be employed. The present work reviews and generalizes some of the main methods of attack of the problem. The results are presented in a form appropriate for numerical computation.

Mohsen, A.↗

NASA NTRS2017

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗

DOE OSTI2021

Neural Lyapunov Control for Power System Transient Stability: A Deep Learning-Based Approach

We report that power system control and transient stability analysis play essential roles in secure system operation. Control of power systems typically involves highly nonlinear and complex dynamics. Most of the existing works address such problems with additional assumptions in system dynamics, leading to a requirement for a complete and general solution. This paper, therefore, proposes a novel control framework for various power system control and stability problems leveraging a learning-based approach. The proposed framework includes a two-module structure that iteratively and jointly learns the candidate Lyapunov function and control law via deep neural networks in a learning module. Meanwhile, it guides the learning procedure towards valid results satisfying Lyapunov conditions in a falsification module. The introduced termination criteria ensure provable system stability. This control framework is verified through several studies handling different types of power system control problems. The results show that the proposed framework is generalizable and can simplify the control design for complex power systems with the stability guarantee and enlarged region of attraction.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1972

Engine selection for transport and combat aircraft

The procedures that are used to select engines for transport and combat aircraft are discussed. In general, the problem is to select the engine parameters including engine size in such a way that all constraints are satisfied and airplane performance is maximized. This is done for four different classes of aircraft: (1) a long haul conventional takeoff and landing (CTOL) transport, (2) a short haul vertical takeoff and landing (VTOL) transport, (3) a long range supersonic transport (SST), and (4) a fighter aircraft. For the commercial airplanes the critical constraints have to do with noise while for the fighter, maneuverability requirements define the engine. Generally, the resultant airplane performance (range or payload) is far less than that achievable without these constraints and would suffer more if nonoptimum engines were selected.

Dugan, J. F., Jr.↗

NASA NTRS1975

The use of ERTS imagery in reservoir management and operation

The author has identified the following significant results. Real time data collection by orbiting satellite relay was found to be both reliable and feasible. ERTS imagery was assessed and it was shown that in most cases better spatial resolution and/or additional spectral bands would be required to satisfy NED's needs. A man-computer interactive system, using cathode ray tube display could solve the problem of an unwieldy mass of data for interpretation.

Cooper, S.↗

NASA NTRS1985

Introduction to cryogenic wind tunnels

The background to the evolution of the cryogenic wind tunnel is outlined, with particular reference to the late 60's/early 70's when efforts were begun to re-equip with larger wind tunnels. The problems of providing full scale Reynolds numbers in transonic testing were proving particularly intractible, when the notion of satisfying the needs with the cryogenic tunnel was proposed, and then adopted. The principles and advantages of the cryogenic tunnel are outlined, along with guidance on the coolant needs when this is liquid nitrogen, and with a note on energy recovery. Operational features of the tunnels are introduced with reference to a small low speed tunnel. Finally the outstanding contributions are highlighted of the 0.3-Meter Transonic Cryogenic Tunnel (TCT) at NASA Langley Research Center, and its personnel, to the furtherance of knowledge and confidence in the concept.

Goodyer, M. J.↗

NASA NTRS1987

Program VSAERO theory document: A computer program for calculating nonlinear aerodynamic characteristics of arbitrary configurations

The VSAERO low order panel method formulation is described for the calculation of subsonic aerodynamic characteristics of general configurations. The method is based on piecewise constant doublet and source singularities. Two forms of the internal Dirichlet boundary condition are discussed and the source distribution is determined by the external Neumann boundary condition. A number of basic test cases are examined. Calculations are compared with higher order solutions for a number of cases. It is demonstrated that for comparable density of control points where the boundary conditions are satisfied, the low order method gives comparable accuracy to the higher order solutions. It is also shown that problems associated with some earlier low order panel methods, e.g., leakage in internal flows and junctions and also poor trailing edge solutions, do not appear for the present method. Further, the application of the Kutta conditions is extremely simple; no extra equation or trailing edge velocity point is required. The method has very low computing costs and this has made it practical for application to nonlinear problems requiring iterative solutions for wake shape and surface boundary layer effects.

Maskew, Brian↗