Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear programming”

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 811 records · Page 45

Parallel Conjugate Gradient: Effects of Ordering Strategies, Programming Paradigms, and Architectural Platforms

The Conjugate Gradient (CG) algorithm is perhaps the best-known iterative technique to solve sparse linear systems that are symmetric and positive definite. A sparse matrix-vector multiply (SPMV) usually accounts for most of the floating-point operations within a CG iteration. In this paper, we investigate the effects of various ordering and partitioning strategies on the performance of parallel CG and SPMV using different programming paradigms and architectures. Results show that for this class of applications, ordering significantly improves overall performance, that cache reuse may be more important than reducing communication, and that it is possible to achieve message passing performance using shared memory constructs through careful data ordering and distribution. However, a multi-threaded implementation of CG on the Tera MTA does not require special ordering or partitioning to obtain high efficiency and scalability.

Oliker, Leonid↗

Turbofan aft duct suppressor study program listing and user's guide

A description of the structure of the Annular Flow Duct Program (AFDP) for the calculation of acoustic suppression due to treatment in a finite length annular duct carrying sheared flow is presented. Although most appropriate for engine exhaust ducts, this program can be used to study sound propagation in any duct that maintains annular geometry over a considerable length of the duct. The program is based on the modal analysis of sound propagation in ducts with axial segments of different wall impedances. For specified duct geometry, wall impedance, flow and acoustic conditions in the duct (including mode amplitude distribution of the source) and duct termination reflection characteristics, the program calculates the suppression due to the treatment in the duct. The presence of forward and backward traveling modes in the duct due to the reflection and redistribution of modes at segment interfaces and duct end terminations are taken into account in the calculations. The effects of thin wall boundary layers (with a linear or mean flow velocity profile) on the acoustic propagation are also included in the program. A functional description of the major subroutines is included and a sample run is provided with an explanation of the output.

Joshi, M. C.↗

Comparison of two weighted integration models for the cueing task: linear and likelihood

In a task in which the observer must detect a signal at two locations, presenting a precue that predicts the location of a signal leads to improved performance with a valid cue (signal location matches the cue), compared to an invalid cue (signal location does not match the cue). The cue validity effect has often been explained with a limited capacity attentional mechanism improving the perceptual quality at the cued location. Alternatively, the cueing effect can also be explained by unlimited capacity models that assume a weighted combination of noisy responses across the two locations. We compare two weighted integration models, a linear model and a sum of weighted likelihoods model based on a Bayesian observer. While qualitatively these models are similar, quantitatively they predict different cue validity effects as the signal-to-noise ratios (SNR) increase. To test these models, 3 observers performed in a cued discrimination task of Gaussian targets with an 80% valid precue across a broad range of SNR's. Analysis of a limited capacity attentional switching model was also included and rejected. The sum of weighted likelihoods model best described the psychophysical results, suggesting that human observers approximate a weighted combination of likelihoods, and not a weighted linear combination.

NASA Program Biomedical Research and Countermeasur↗

ANISAP: A three-dimensional finite element program for laminated composites subjected to mechanical loading

ANISAP is a 3-D finite element FORTRAN 77 computer code for linear elastic, small strain, analysis of laminated composites with arbitrary geometry including free edges and holes. Individual layers may be isotropic or transversely isotropic in material principal coordinates; individual layers may be rotated off-axis about a global z-axis. The laminate may be a hybrid. Three different isoparametric elements, variable order of gaussian integration, calculation of stresses at element boundaries, and loading by either nodal displacement of forces are included in the program capability. Post processing capability includes failure analysis using the tensor polynominal failure criterion.

Burns, S. W.↗

Effects of Ordering Strategies and Programming Paradigms on Sparse Matrix Computations

The Conjugate Gradient (CG) algorithm is perhaps the best-known iterative technique to solve sparse linear systems that are symmetric and positive definite. For systems that are ill-conditioned, it is often necessary to use a preconditioning technique. In this paper, we investigate the effects of various ordering and partitioning strategies on the performance of parallel CG and ILU(O) preconditioned CG (PCG) using different programming paradigms and architectures. Results show that for this class of applications: ordering significantly improves overall performance on both distributed and distributed shared-memory systems, that cache reuse may be more important than reducing communication, that it is possible to achieve message-passing performance using shared-memory constructs through careful data ordering and distribution, and that a hybrid MPI+OpenMP paradigm increases programming complexity with little performance gains. A implementation of CG on the Cray MTA does not require special ordering or partitioning to obtain high efficiency and scalability, giving it a distinct advantage for adaptive applications; however, it shows limited scalability for PCG due to a lack of thread level parallelism.

Oliker, Leonid↗

Improvements to the missile aerodynamic prediction code DEMON3

The computer program DEMON3 was developed for the aerodynamic analysis of nonconventional supersonic configurations comprising a body with noncircular cross section and up to two wing or fin sections. Within a wing or fin section, the lifting surfaces may be cruciform, triform, planar, or low profile layouts; the planforms of the lifting surfaces allow for breaks in sweep. The body and fin sections are modeled by triplet and constant u-velocity panels, respectively, accounting for mutual body-fin interference. Fin thickness effects are included for the use of supersonic planar source panels. One of the unique features of DEMON3 is the modeling of high angle of attack vortical effects associated with the lifting surfaces and the body. In addition, shock expansion and Newtonian pressure calculation methods can be optionally engaged. These two dimensional nonlinear methods are augmented by aerodynamic interference determined from the linear panel methods. Depending on the geometric details of the body, the DEMON3 program can be used to analyze nonconventional configurations at angles of attack up to 25 degrees for Mach numbers from 1.1 to 6. Calculative results and comparisons with experimental data demonstrate the capabilities of DEMON3. Limitations and deficiencies are listed.

Dillenius, Marnix F. E.↗

HYTESS-Hypothetical Turbofan-Engine Simplified Simulation

Simulated characteristics mimic those of advanced turbofan engines. Computer program developed to offer those interested in engine dynamics and controls research efficient, realistic, and easily-used engine simulation. Simulation developed from linearized operating-point models but still retains essential nonlinear engine effects. Representative of hypothetical, low-bypass-ratio, twin-spool, axialflow turbofan engine. Program written in FORTRAN IV.

Merrill, W.↗

qsp4pde v1.0

The software program is a collection of Python implementations of quantum circuits for solving linear partial differential equations using quantum signal processing. The circuits are implemented using the qiskit SDK.

Kim, Hyeongjin [Lawrence Berkeley National Laborat↗

Study of dynamic response of elastic space stations

Analytical procedure and the requisite computer programs compute the dynamic responses of two large elastic space stations. The linearized equations of motion are derived from lagranges equations. Then the normal modes of free vibration of the nonrotating space station are used to define the elastic degrees of freedom.

Kamrath, P.↗

Numerical analysis of stiffened shells of revolution. Volume 1: Theory manual for STARS-2S, 2B, 2V digital computer programs

The theoretical analysis background for the STARS-2 (shell theory automated for rotational structures) program is presented. The theory involved in the axisymmetric nonlinear and unsymmetric linear static analyses, and the stability and vibrations (including critical rotation speed) analyses involving axisymmetric prestress are discussed. The theory for nonlinear static, stability, and vibrations analyses, involving shells with unsymmetric loadings are included.

Svalbonas, V.↗

Structural design for dynamic response reduction

A computer program for redesigning structural modes to reduce response has been initiated. The linear regulator approach in modal coordinates has been implemented. It is noted that the transformation of solution to physical structure is a major problem. It is concluded that the solution of stiffness equations and damping equations can be done separately as NXN set of (matrix Riccati) equations.

Hanks, B. R.↗

Time-independent Anisotropic Plastic Behavior by Mechanical Subelement Models

The paper describes a procedure for modelling the anisotropic elastic-plastic behavior of metals in plane stress state by the mechanical sub-layer model. In this model the stress-strain curves along the longitudinal and transverse directions are represented by short smooth segments which are considered as piecewise linear for simplicity. The model is incorporated in a finite element analysis program which is based on the assumed stress hybrid element and the iscoplasticity-theory.

Pian, T. H. H.↗

Prediction of LDEF ionizing radiation environment

The Long Duration Exposure Facility (LDEF) spacecraft flew in a 28.5 deg inclination circular orbit with an altitude in the range from 172 to 258.5 nautical miles. For this orbital altitude and inclination two components contribute most of the penetrating charge particle radiation encountered - the galactic cosmic rays and the geomagnetically trapped Van Allen protons. Where shielding is less than 1.0 g/sq cm geomagnetically trapped electrons make a significant contribution. The 'Vette' models together with the associated magnetic filed models were used to obtain the trapped electron and proton fluences. The mission proton doses were obtained from the fluence using the Burrell proton dose program. For the electron and bremsstrahlung dose we used the Marshall Space Flight Center (MSFC) electron dose program. The predicted doses were in general agreement with those measured with on-board thermoluminescent detector (TLD) dosimeters. The NRL package of programs, Cosmic Ray Effects on MicroElectronics (CREME), was used to calculate the linear energy transfer (LET) spectrum due to galactic cosmic rays (GCR) and trapped protons for comparison with LDEF measurements.

Watts, John W.↗

GRO: Black hole models for gamma-ray bursts

The Burst and Transient Source Experiment (BATSE) on board the Compton Gamma Ray Observatory (CGRO) has established that the distribution of gamma-ray bursts (GRB's) is isotropic but is bound radially. This finding suggests that the bursts are either cosmological or they originate from an extended Galactic halo. The implied luminosities and the observed variability of the GRB's on time scales as short as one millisecond suggest that they originate from compact objects. We are presently studying black hole models for GRB's. Any such model must produce a non-thermal photon spectrum to agree with the observed properties. For a wide range of burst parameters the assumed bursting source consists of a non-thermal electron-positron-photon plasma of very high density. It seems possible to produce such a plasma in accretion onto black holes. In our on-going work, we are developing the kinetic theory for a non-equilibrium pair plasma. The main new features of our work are as follows: (1) We do not assume the presence of a thermal electron bath. (2) Non-thermal, high-energy pairs are allowed to have an arbitrary concentration and energy distribution. (3) There is no soft photon source in our model; initially all the photons in the plasma are either energetic X-rays or gamma-rays. (4) The initial energy distribution of the pairs as well as photons is arbitrary. (5) We collect the analytical expressions for the kinetic kernels for all relevant processes. And (6) we present a different approach to finding the time-evolution of pair and photon spectra, which is a combination of the kinetic-theory and the non-linear Monte-Carlo schemes. We have developed many Monte-Carlo programs to model various process, to take into account the time evolution, and to incorporate various physical effects which are unique to non-thermal plasmas. The hydrodynamics of fireballs in GRB's was studied before. Applying results from kinetic theory will improve our understanding of these systems.

Ruderman, Malvin↗

Application of Hybrid Optimization-Expert System for Optimal Power Management on Board Space Power Station

The space power system has two sources of energy: photo-voltaic blankets and batteries. The optimal power management problem on-board has two broad operations: off-line power scheduling to determine the load allocation schedule of the next several hours based on the forecast of load and solar power availability. The nature of this study puts less emphasis on speed requirement for computation and more importance on the optimality of the solution. The second category problem, on-line power rescheduling, is needed in the event of occurrence of a contingency to optimally reschedule the loads to minimize the 'unused' or 'wasted' energy while keeping the priority on certain type of load and minimum disturbance of the original optimal schedule determined in the first-stage off-line study. The computational performance of the on-line 'rescheduler' is an important criterion and plays a critical role in the selection of the appropriate tool. The Howard University Center for Energy Systems and Control has developed a hybrid optimization-expert systems based power management program. The pre-scheduler has been developed using a non-linear multi-objective optimization technique called the Outer Approximation method and implemented using the General Algebraic Modeling System (GAMS). The optimization model has the capability of dealing with multiple conflicting objectives viz. maximizing energy utilization, minimizing the variation of load over a day, etc. and incorporates several complex interaction between the loads in a space system. The rescheduling is performed using an expert system developed in PROLOG which utilizes a rule-base for reallocation of the loads in an emergency condition viz. shortage of power due to solar array failure, increase of base load, addition of new activity, repetition of old activity etc. Both the modules handle decision making on battery charging and discharging and allocation of loads over a time-horizon of a day divided into intervals of 10 minutes. The models have been extensively tested using a case study for the Space Station Freedom and the results for the case study will be presented. Several future enhancements of the pre-scheduler and the 'rescheduler' have been outlined which include graphic analyzer for the on-line module, incorporating probabilistic considerations, including spatial location of the loads and the connectivity using a direct current (DC) load flow model.

Momoh, James↗

Extraction of Lateral-Directional Stability and Control Derivatives for the Basic F-18 Aircraft at High Angles of Attack

The results of parameter identification to determine the lateral-directional stability and control derivatives of an F-18 research aircraft in its basic hardware and software configuration are presented. The derivatives are estimated from dynamic flight data using a specialized identification program developed at NASA Dryden Flight Research Center. The formulation uses the linearized aircraft equations of motions in their continuous/discrete form and a maximum likelihood estimator that accounts for both state and measurement noise. State noise is used to model the uncommanded forcing function caused by unsteady aerodynamics, such as separated and vortical flows, over the aircraft. The derivatives are plotted as functions of angle of attack between 3 deg and 47 deg and compared with wind-tunnel predictions. The quality of the derivative estimates obtained by parameter identification is somewhat degraded because the maneuvers were flown with the aircraft's control augmentation system engaged, which introduced relatively high correlations between the control variables and response variables as a result of control motions from the feedback control system.

Iliff, Kenneth W.↗