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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A reduced order observer based controller design for H-infinity optimization

In this paper the H-infinity control problem is investigated. It is well-known that for this problem, in general, controllers of the same dynamic order as the given system are needed. However, in the case that the standard assumptions on two direct feedthrough matrices are not satisfied, it is here shown that dynamic compensators of a lower dynamical order can be found. This result can be derived by using the standard reduced order observer based controllers in the case that one or more states are measured without noise.

Stoorvogel, Anton A.↗
NASA NTRS1994

A hydrodynamic approach to cosmology: The mixed dark matter cosmological scenario

We compute the evolution of spatially flat, mixed cold and hot dark matter models containing both baryonic matter and two kinds of dark matter. Hydrodynamics is treated with a highly developed Eulerian hydrodynamic code (see Cen 1992). A standard particle-mesh (PM) code is also used in parallel to calculate the motion of the dark matter components. We adopt the following parameters: h equivalent to (sub 0)/100 km/s Mpc(exp -1) = 0.5, OMEGA(sub C) = 0.3, and OMEGA(sub B) = 0.06, with amplitude of the perturbation spectrum fixed by the Cosmic Background Explorer Satellite (COBE) Dark Matter Radiation (DMR) measurements (Smoot et al. 1992) being sigma (sub 8) = 0.67. Four different boxes are simulated with box sizes of L = (64, 16, 4, 1) h(exp -1) Mpc, respectively, the two small boxes providing good resolution but little valid information due to absence of large-scale power. We use 128(exp 3) approximate 10(exp 6.3) baryonic cells, 128(exp .3) cold dark matter particles, and 2 x 128(exp 3) hot dark matter particles. In addition to the dark matter we follow separately six baryonic species (H, H(+), He, He(+), He(++), e(-)) with allowance for both (nonequilibrium) collisional and radiative ionization in every cell. The background radiation field is also followed in detail with allowance made for both continuum and line processes, to allow nonequilibrium heating and cooling processes to be followed in detail. The mean final Zeldovich-Sunyaev y parameter is estimated to be y Bar = (5.4 + or - 2.7) x 10(exp -7) below currently attainable observations, with a rms fluctuation of approximately delta bar y = (0.6 + or - 3.0) x 10(exp -7) on arcminute scales. The rate of galaxy formation peaks at an even later epoch (z approximate 0.3) than in the standard (OMEGA = 1, sigma sub 8 = 0.67) cold dark matter (CDM) model (z approximate 0.5) and, at a redshift of z = 4, is nearly a factor of 100 lower than for the CDM model with the same value of sigma sub 8. With regard to mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum has a broad peak in the vicinity of M(sub B) = 10(exp 9.5) solar mass with a reasonable fit to the Schechter luminosity function if the ratio of baryon mass to blue light is approximately 4. In addition, one very large PM simulation was made in a box with size (320 h(exp - 1) Mpc) containing 3 x 200(exp 3) = 10(exp 7.4) particles. Utilizing this simulation we find that the model yields a cluster mass function which is about a factor of 4 higher than observed, but a cluster-cluster correlation length marginally lower than observed, but that both are closer to observations than in the (COBE) normalized CDM model. The one-dimensional pairwise velocity dispersion is 605 + or - 8 km/s at 1/h separation, lower than that of the DCM model normalized to COBE, but still significant higher than observations (Davis & Peebles 1983). A plausible velocity bias b(sub v) = 0.8 + or - 0.1 on this scale will reduce but not remove the discrepancy. The velocity auto-correlat ion function has a coherence length of 40/h Mpc, which is somewhat lower than the observed counterpart. In all these respects the model would be improved by decreasing the cold fraction of the dark OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDB). But formation of galaxies and clusters of galaxies is much later in this model than in COBE-normalized CDM, perhaps too late. To improve on these constraints a larger ratio of OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDM)) is required than the value of 0.67 adopted here. It does not seem possible to find a value for this ratio which would satisfy all tests. Overall, the model is similar both on large and intermediate scales to the standard CDM model normalized to the same value of sigma(sub B), but the problem with regard to late formation of galaxies is more severe in this model than in that CDM model. Adding hot dark matter, significantly improves the ability of the COBE-normalized CDM scenario to fit existing observations, but the model is in fact not as good as the CDM model with the same sigma(sub 8) and is still probably unsatisfactory with regard to several critical tests.

Cen, Renyue↗
NASA NTRS1989

Multiplex gas chromatography: an alternative concept for gas chromatographic analysis of planetary atmospheres

Gas chromatography (GC) is a powerful technique for analyzing gaseous mixtures. Applied to the earth's atmosphere, GC can be used to determine the permanent gases--such as carbon dioxide, nitrogen, and oxygen--and to analyze organic pollutants in air. The U.S. National Aeronautics and Space Administration (NASA) has used GC in spacecraft missions to Mars (the Viking Biology Gas Exchange Experiment [GEX] and the Viking Gas Chromatograph-Mass Spectrometer [GC-MS]) and to Venus (the Pioneer Venus Gas Chromatograph [PVGC] on board the Pioneer Venus sounder probe) for determining the atmospheric constituents of these two planets. Even though conventional GC was very useful in the Viking and Pioneer missions, spacecraft constraints and limitations intrinsic to the technique prevented the collection of more samples. With the Venus probe, for instance, each measurement took a relatively long time to complete (10 min), and successive samples could not be introduced until the previous samples had left the column. Therefore, while the probe descended through the Venusian atmosphere, only three samples were acquired at widely separated altitudes. With the Viking mission, the sampling rate was not a serious problem because samples were acquired over a period of one year. However, the detection limit was a major disadvantage. The GC-MS could not detect simple hydrocarbons and simple alcohols below 0.1 ppm, and the GEX could not detect them below 1 ppm. For more complex molecules, the detection limits were at the parts-per-billion level for both instruments. Finally, in both the Viking and Pioneer missions, the relatively slow rate of data acquisition limited the number of analyses, and consequently, the amount of information returned. Similar constraints are expected in future NASA missions. For instance, gas chromatographic instrumentation is being developed to collect and analyze organic gases and aerosols in the atmosphere of Titan (one of Saturn's satellites). The Titan-Cassini entry probe, which is being jointly planned by NASA and the European Space Agency (ESA), might be launched as early as 1994. As in the Pioneer mission, limited time--perhaps only 3-4 h--will be available for the completion of all analyses while the probe descends through the atmosphere. A conventional GC or GC-MS system would be able to analyze no more than two aerosol and two gas samples during the probe's descent. Conventional GC also is limited by the sensitivity of the detector and by the sample volume. For the Titan mission, the sensitivity problems will be worse because the atmospheric pressure at the time of instrument deployment is expected to be < 3 torr. Consequently, the sample volume might not be large enough to satisfy the detector sensitivity requirements. Because of such limitations, alternative GC analysis techniques have been investigated for future NASA missions. Multiplex gas chromatography has been investigated as a possible candidate for chemical analysis within a spacecraft or other restricted environment, and chemical modulators have been developed and used when needed with this technique to reduce the size and weight of the instrumentation. Also, several new multiplex techniques have been developed for use in specific applications.

NASA Center ARC↗
NASA NTRS1978

Determination of launch conditions for Spacelab satisfying the constraints of an absorption spectroscopy project

The conditions for absorption spectroscopy of the homosphere between 20 and 100 km from a space platform such as Spacelab are considered. Using nominal orbital elements, and taking launch data and time as parameters of the problem, formula were established that give the data necessary for determining the flight plan. The geographical coordinates of the tangential points of the solar illumination boundaries were also calculated. It is shown that the latitude coverage of the observation is closely related to launch conditions. Expressions for the pointing angles referred to an instantaneous terrestrial trihedron are given. These angular elements are needed to determine the real pointing angles in a reference system linked to Spacelab.

Vercheval, J.↗
NASA NTRS1976

Solution of the minimum time-to-climb problem by matched asymptotic expansions

Application of singular perturbation techniques to trajectory optimization problems of flight mechanics is discussed. The method of matched asymptotic expansions is used to obtain an approximate solution to the aircraft minimum time-to-climb problem. Outer, boundary-layer, and composite solutions are obtained to zeroth and first orders. A stability criterion is derived for the zeroth-order boundary-layer solutions (the theory requires a form of boundary-layer stability). A numerical example is considered for which it is shown that the stability criterion is satisfied and a useful numerical solution is obtained. The zeroth-order solution proves to be a poor approximation, but the first-order solution gives a good approximation for both the trajectory and the minimum time-to-climb. The computational cost of the singular perturbation solution is considerably less than that of a steepest descent solution. Thus singular perturbation methods appear to be promising for the solution of optimal control problems.

Ardema, M. D.↗
DOE OSTI2023

Bell inequalities in 2-2 scattering

We consider Bell inequalities in 2-2 scattering of photons, gravitons, fermions, and pions. We choose measurement settings that give maximum Bell violation for maximally entangled states and calculate the relevant Bell inequalities for these processes. For photon scattering at low energies, quantum electrodynamics exhibits Bell violation for all scattering angles except for a small transverse region. This leads to a fine-tuning problem. Incorporating a light axion/axionlike particle removes the fine-tuning problem and constrains the axion-coupling–axion-mass parameters. Allowing for graviton exchange and demanding Bell violation in photon scattering, we find that the weak-gravity conjecture is satisfied. Quantum gravity effects on axion coupling is discussed. For 2-2 graviton scattering, we find that CEMZ bounds allow for, at most, small Bell violations. Restriction on the Weinberg angle is found by demanding Bell violation in Bhabha scattering. We use recent $S$-matrix bootstrap data for pions and photons to study the Bell parameter in the space of allowed $S$-matrices. In the photon case, we study the Bell parameters as a function of energy and find support for the effective field theory observations. We discuss Bell parameter for pion $S$-matrices, which are qutrits. For pions, we find that there is a minimization of a suitable Bell parameter for $S$-matrices which exhibit Regge behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1986

Calculations of the stability of some axisymmetric flows proposed as a model of vortex breakdown

The term vortex breakdown refers to the abrupt and drastic changes of structure that can sometimes occur in swirling flows. It was conjectured that the bubble type of breakdown can be viewed as an axisymmetric wave traveling upstream in a primarily columnar vortex flow. In this scenario the wave's upstream progress is impeded only when it reaches a critical amplitude and it loses stability to some nonaxisymmetric disturbance. The stability of some axisymmetric wavy flows to three dimensional disturbances, viewing the amplitude of the wave as a bifurcation parameter is examined. The stability of a set of related columnar vortex flows, constructed by taking the two dimensional flow at a single axial location and extending it throughout the domain without variation, is investigated. The method used will be to expand the perturbation velocity in a series of divergence free vectors which ensures that the continuity equation for the incompressible fluid is satisfied exactly by the computed velocity field. Projections of the stability equation onto the space of inviscid vector fields eliminated the pressure term from the equation and reduces the differential eigen problem to a generalized matrix eigen problem. Results are presented both for the one dimensional, columnar vortex flows and also for the wavy bubble flow.

Mhuiris, N. M. G.↗
NASA NTRS1984

Applications of CONMIN to wing design optimization with vortex flow effect

Slender wings on supersonic cruise configurations are expected to be thin and highly swept. As a result, edge-separated vortex flow is inevitable and must be accounted for in aerodynamic analysis and design. The present method is based on the method of suction analogy to calculate the total aerodynamic characteristics. The method requires the solution of the attached flow problem, the latter being solved by a low-order panel method in subsonic and supersonic flow. In essence, the lifting pressure is calculated by using a pressure-doublet distribution satisfying the Prandtl-Glauert equation. From the pressure distribution, the leading-edge suction is calculated. The latter is assumed to be the vortex lift through the method of suction analogy. For a cambered wing, the location of vortex-lift action point is important in predicting the aerodynamic characteristics. It is also seen that the effect of camber shape appears nonlinearly in all aerodynamic expressions. To design the camber shape, the camber slope is represented by a cosine Fourier series at each of several spanwise stations. The Fourier coefficients are the design variables. To design a leading-edge flap in the vortex flow (i.e., a vortex flap), the coordinates of corner points and the deflection angle are the design variables. The process of wing design is to determine the camber shape and twist distribution such that an objective function, typically the drag, is minimized, subject to various constraints.

Lan, C. E.↗
NASA NTRS1987

Optimization of orbital assignment and specification of service areas in satellite communications

The mathematical nature of the orbital and frequency assignment problem for communications satellites is explored, and it is shown that choosing the correct permutations of the orbit locations and frequency assignments is an important step in arriving at values which satisfy the signal-quality requirements. Two methods are proposed to achieve better spectrum/orbit utilization. The first, called the delta S concept, leads to orbital assignment solutions via either mixed-integer or restricted basis entry linear programming techniques; the method guarantees good single-entry carrier-to-interference ratio results. In the second, a basis for specifying service areas is proposed for the Fixed Satellite Service. It is suggested that service areas should be specified according to the communications-demand density in conjunction with the delta S concept in order to enable the system planner to specify more satellites and provide more communications supply.

Wang, Cou-Way↗
NASA NTRS1990

Generalized gradient algorithm for trajectory optimization

The generalized gradient algorithm presented and verified as a basis for the solution of trajectory optimization problems improves the performance index while reducing path equality constraints, and terminal equality constraints. The algorithm is conveniently divided into two phases, of which the first, 'feasibility' phase yields a solution satisfying both path and terminal constraints, while the second, 'optimization' phase uses the results of the first phase as initial guesses.

Zhao, Yiyuan↗
NASA NTRS1991

A game theoretic approach to a finite-time disturbance attenuation problem

A disturbance attenuation problem over a finite-time interval is considered by a game theoretic approach where the control, restricted to a function of the measurement history, plays against adversaries composed of the process and measurement disturbances, and the initial state. A zero-sum game, formulated as a quadratic cost criterion subject to linear time-varying dynamics and measurements, is solved by a calculus of variation technique. By first maximizing the quadratic cost criterion with respect to the process disturbance and initial state, a full information game between the control and the measurement residual subject to the estimator dynamics results. The resulting solution produces an n-dimensional compensator which expresses the controller as a linear combination of the measurement history. A disturbance attenuation problem is solved based on the results of the game problem. For time-invariant systems it is shown that under certain conditions the time-varying controller becomes time-invariant on the infinite-time interval. The resulting controller satisfies an H(infinity) norm bound.

Rhee, Ihnseok↗
NASA NTRS1992

Observations on computational methodologies for use in large-scale, gradient-based, multidisciplinary design incorporating advanced CFD codes

How a combination of various computational methodologies could reduce the enormous computational costs envisioned in using advanced CFD codes in gradient based optimized multidisciplinary design (MdD) procedures is briefly outlined. Implications of these MdD requirements upon advanced CFD codes are somewhat different than those imposed by a single discipline design. A means for satisfying these MdD requirements for gradient information is presented which appear to permit: (1) some leeway in the CFD solution algorithms which can be used; (2) an extension to 3-D problems; and (3) straightforward use of other computational methodologies. Many of these observations have previously been discussed as possibilities for doing parts of the problem more efficiently; the contribution here is observing how they fit together in a mutually beneficial way.

Newman, P. A.↗
NASA NTRS1994

Coupled 2-dimensional cascade theory for noise and unsteady aerodynamics of blade row interaction in turbofans. Volume 1: Theory development and parametric studies

Typical analytical models for interaction between rotor and stator in a turbofan analyze the effect of wakes from the rotor impinging on the stator, producing unsteady loading, and thereby generating noise. Reflection/transmission characteristics of the rotor are sometimes added in a separate calculation. In those models, there is a one-to-one relationship between wake harmonics and noise harmonics; that is, the BPF (blade passing frequency) wake harmonic causes only the BPF noise harmonic, etc. This report presents a more complete model in which flow tangency boundary conditions are satisfied on two cascades in relative motion for several harmonics simultaneously. By an extension of S.N. Smith's code for two dimensional flat plate cascades, the noise generation/frequency scattering/blade row reflection problem is solved in a single matrix inversion. It is found that the BPF harmonic excitation of the stator scatters considerable energy in the higher BPF harmonics due to relative motion between the blade rows. Furthermore, when swirl between the rotor and stator is modeled, a 'mode trapping' effect occurs which explains observations on fans operating at rotational speeds below BFP cuton: the BPF mode amplifies between blade rows by multiple reflections but cannot escape to the inlet and exit ducts. However, energy scattered into higher harmonics does propagate and dominates the spectrum at two and three times BPF. This report presents the complete derivation of the theory, comparison with a previous (more limited) coupled rotor/stator interaction theory due to Kaji and Okazaki, exploration of the mode trapping phenomenon, and parametric studies showing the effects of vane/blade ratio and rotor/stator interaction. For generality, the analysis applies to stages where the rotor is either upstream or downstream of the stator and to counter rotation stages. The theory has been coded in a FORTRAN program called CUP2D, documented in Volume 2 of this report. It is concluded that the new features of this analysis - unsteady coupling, frequency scattering, and flow turning between rotor and stator - have a profound effect on noise generation caused by rotor/stator interaction. Treating rotors and stators as isolated cascades is not adequate for noise analysis and prediction.

Hanson, Donald B.↗
NASA NTRS1994

Microgravity isolation system design: A case study

Many acceleration-sensitive, microgravity science experiments will require active vibration isolation from manned orbiters on which they will be mounted. The isolation problem, especially in the case of a tethered payload, is a complex three-dimensional one that is best suited to modern-control design methods. In this paper, extended H(sub 2) synthesis is used to design an active isolator (i.e., controller) for a realistic single-input-multiple-output (SIMO) microgravity vibration isolation problem. Complex mu-analysis methods are used to analyze the isolation system with respect to sensor, actuator, and umbilical uncertainties. The paper fully discusses the design process employed and the insights gained. This design case study provides a practical approach for isolation problems of greater complexity. Issues addressed include a physically intuitive state-space description of the system, disturbance and noise filters, filters for frequency weighting, and uncertainty models. The controlled system satisfies all the performance specifications and is robust with respect to model uncertainties.

Hampton, R. D.↗
NASA NTRS2001

Theoretical Investigation of Thermo-Mechanical Behavior of Carbon Nanotube-Based Composites Using the Integral Transform Method

This project uses the integral transform technique to model the problem of nanotube behavior as an axially symmetric system of shells. Assuming that the nanotube behavior can be described by the equations of elasticity, we seek a stress function x which satisfies the biharmonic equation: del(exp 4) chi = [partial deriv(r(exp 2)) + partial deriv(r) + partial deriv(z(exp 2))] chi = 0. The method of integral transformations is used to transform the differential equation. The symmetry with respect to the z-axis indicates that we only need to consider the sine transform of the stress function: X(bar)(r,zeta) = integral(from 0 to infinity) chi(r,z)sin(zeta,z) dz.

Pawloski, Janice S.↗
NASA NTRS1992

Least-squares finite element solution of 3D incompressible Navier-Stokes problems

Although significant progress has been made in the finite element solution of incompressible viscous flow problems. Development of more efficient methods is still needed before large-scale computation of 3D problems becomes feasible. This paper presents such a development. The most popular finite element method for the solution of incompressible Navier-Stokes equations is the classic Galerkin mixed method based on the velocity-pressure formulation. The mixed method requires the use of different elements to interpolate the velocity and the pressure in order to satisfy the Ladyzhenskaya-Babuska-Brezzi (LBB) condition for the existence of the solution. On the other hand, due to the lack of symmetry and positive definiteness of the linear equations arising from the mixed method, iterative methods for the solution of linear systems have been hard to come by. Therefore, direct Gaussian elimination has been considered the only viable method for solving the systems. But, for three-dimensional problems, the computer resources required by a direct method become prohibitively large. In order to overcome these difficulties, a least-squares finite element method (LSFEM) has been developed. This method is based on the first-order velocity-pressure-vorticity formulation. In this paper the LSFEM is extended for the solution of three-dimensional incompressible Navier-Stokes equations written in the following first-order quasi-linear velocity-pressure-vorticity formulation.

Jiang, Bo-Nan↗
NASA NTRS2005

SOFIA'S Challenge: Scheduling Airborne Astronomy Observations

The Stratospheric Observatory for Infrared Astronomy (SOFIA) is NASA's next generation airborne astronomical observatory, and will commence operations in 2005. The facility consists of a 747-SP modified to accommodate a 2.5 meter telescope. SOFIA is expected to fly an average of 140 science flights per year over its 20 year lifetime. Depending on the nature of the instrument used during flight, 5-15 observations per flight are expected. The SOFIA telescope is mounted aft of the wings on the port side of the aircraft and is articulated through a range of 20deg to 60deg of elevation. The telescope has minimal lateral flexibility; thus, the aircraft must turn constantly to maintain the telescope's focus on an object during observations. A significant problem in future SOFIA operations is that of scheduling flights in support of observations. Investigators are expected to propose small numbers of observations, and many observations must be grouped together to make up single flights. Flight planning for the previous generation airborne observatory, the Kuiper Airborne Observatory (KAO), was done by hand; planners had to choose takeoff time, observations to perform, and decide on setup-actions (called "dead-legs") to position the aircraft prior to observing. This task frequently required between 6-8 hours to plan one flight The scope of the flight planning problem for supporting GI observations with the anticipated flight rate for SOFIA makes the manual approach for flight planning daunting. In response, we have designed an Automated Flight Planner (AFP) that accepts as input a set of requested observations, designated flight days, weather predictions and fuel limitations, and searches automatically for high-quality flight plans that satisfy all relevant aircraft and astronomer specified constraints. The AFP can generate one candidate flight plan in 5-10 minutes, of computation time, a feat beyond the capabilities of human flight planners. The rate at which the AFP can generate flights enables humans to assess and analyze complex tradeoffs between fuel consumption, estimated science quality and the percentage of scheduled observations. Due to the changing nature of SOFIA scheduling problems, this functionality will play a crucial role in optimizing science and minimizing costs during operations. In the full paper, we will summarize the technical challenges that have been met in order to build this system. These include: design of the search algorithm, design of appropriate heuristics and approximations, and reduction in the size of the search space. We will also describe technical challenges that are currently being addressed, including the extension of the existing approach to handle new solution criteria. Finally, we will describe a variety of cultural challenges that the astronomical community must address in order to successfully use SOFIA, and describe how the AFT can be used to address some of these challenges. Specifically, many of the intended science users are accustomed to using ground-based or space-based observatories; we will identify some differences that arise due to the nature of airborne observatories, and how the AFT can be extended to provide useful services to ease these cultural differences.

Frank, Jeremy↗
DOE OSTI2023

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗