Search NASA⌕ Search

SEARCH · Search NASA

Results for “Problem structuring methods”

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 1,765 records · Page 98

Interplanetary Propagation of Coronal Mass Ejections

Although more than ten thousand coronal mass ejections (CMEs) are produced during each solar cycle at the Sun, only a small fraction hits the Earth. Only a small fraction of the Earth-directed CMEs ultimately arrive at Earth depending on their interaction with the solar wind and other large-scale structures such as coronal holes and CMEs. The interplanetary propagation is essentially controlled by the drag force because the propelling force and the solar gravity are significant only near the Sun. Combined remote-sensing and in situ observations have helped us estimate the influence of the solar wind on the propagation of CMEs. However, these measurements have severe limitations because the remote-sensed and in-situ observations correspond to different portions of the CME. Attempts to overcome this problem are made in two ways: the first is to model the CME and get the space speed of the CME, which can be compared with the in situ speed. The second method is to use stereoscopic observation so that the remote-sensed and in-situ observations make measurements on the Earth-arriving part of CMEs. The Solar Terrestrial Relations Observatory (STEREO) mission observed several such CMEs, which helped understand the interplanetary evolution of these CMEs and to test earlier model results. This paper discusses some of these issues and updates the CME/shock travel time estimates for a number of CMEs.

Gopalswamy, Nat↗

Spatio-Temporal Video Segmentation with Shape Growth or Shrinkage Constraint

We propose a new method for joint segmentation of monotonously growing or shrinking shapes in a time sequence of noisy images. The task of segmenting the image time series is expressed as an optimization problem using the spatio-temporal graph of pixels, in which we are able to impose the constraint of shape growth or of shrinkage by introducing monodirectional infinite links connecting pixels at the same spatial locations in successive image frames. The globally optimal solution is computed with a graph cut. The performance of the proposed method is validated on three applications: segmentation of melting sea ice floes and of growing burned areas from time series of 2D satellite images, and segmentation of a growing brain tumor from sequences of 3D medical scans. In the latter application, we impose an additional intersequences inclusion constraint by adding directed infinite links between pixels of dependent image structures.

Segmentation↗

Statistical analysis of the modal properties of large structural systems.

A theory is developed to predict eigenvalue and eigenvector uncertainty in large dynamic models. The uncertainty is based on physical property uncertainty and should not be confused with numerical roundoff, although the method can be extended to include the latter. The theory, when implemented on a computer, is used to analyze the uncertainties in frequencies and mode shapes based on uncertainties in mass, stiffness, modulus of elasticity, etc. The method incorporates a linear statistical model which is quite adequate for handling property uncertainties of 10% or more. The model is not limited to small systems but uses certain statistical assumptions as well as selective matrix manipulations to keep the size of all matrix operations to within the number of degrees of freedom of the system. Examples are given for two longitudinal vibration problems, and the results are supported by a Monte Carlo analysis.

Collins, J. D.↗

Scattering from dielectric structures above impedance surfaces and resistive sheets

Interest in understanding of electromagnetic interaction with rough surfaces has prompted the study of scattering from typical dielectric humps over impedance surfaces. It is shown that the Green's function of the problem for a resistive sheet resembles that of the impedance surface. Hence both problems are considered here. A numerical solution for the scattered field of a two-dimensional dielectric object, possibly inhomogeneous, with arbitrary cross section above the impedance surface or resistive sheet is sought. First the Green's function of the problem is derived based on the exact image theory. This form of the Green's function is amenable to numerical computation. Then the induced polarization currents are calculated by casting the integral equations into a matrix equation via the method of moments. Numerical problems in calculation of the Green's function when both source and observation points are close to the surface are discussed. Comparison of numerical results with a perturbation solution shows excellent agreement between the two methods.

Sarabandi, Kamal↗

The dynamics and control of large flexible space structures - 12, supplement 11

The rapid 2-D slewing and vibrational control of the unsymmetrical flexible SCOLE (Spacecraft Control Laboratory Experiment) with multi-bounded controls is considered. Pontryagin's Maximum Principle is applied to the nonlinear equations of the system to derive the necessary conditions for the optimal control. The resulting two point boundary value problem is then solved by using the quasilinearization technique, and the near minimum time is obtained by sequentially shortening the slewing time until the controls are near the bang-bang type. The tradeoff between the minimum time and the minimum flexible amplitude requirements is discussed. The numerical results show that the responses of the nonlinear system are significantly different from those of the linearized system for rapid slewing. The SCOLE station-keeping closed loop dynamics are re-examined by employing a slightly different method for developing the equations of motion in which higher order terms in the expressions for the mast modal shape functions are now included. A preliminary study on the effect of actuator mass on the closed loop dynamics of large space systems is conducted. A numerical example based on a coupled two-mass two-spring system illustrates the effect of changes caused in the mass and stiffness matrices on the closed loop system eigenvalues. In certain cases the need for redesigning control laws previously synthesized, but not accounting for actuator masses, is indicated.

Bainum, Peter M.↗

Transient excitation and mechanical admittance test techniques for prediction of payload vibration environments

Transient excitation forces were applied separately to simple beam-and-mass launch vehicle and payload models to develop complex admittance functions for the interface and other appropriate points on the structures. These measured admittances were then analytically combined by a matrix representation to obtain a description of the coupled system dynamic characteristics. Response of the payload model to excitation of the launch vehicle model was predicted and compared with results measured on the combined models. These results are also compared with results of earlier work in which a similar procedure was employed except that steady-state sinusoidal excitation techniques were included. It is found that the method employing transient tests produces results that are better overall than the steady state methods. Furthermore, the transient method requires far less time to implement, and provides far better resolution in the data. However, the data acquisition and handling problem is more complex for this method. It is concluded that the transient test and admittance matrix prediction method can be a valuable tool for development of payload vibration tests.

Kana, D. D.↗

An algorithm for optimal structural design with frequency constraints

The paper presents a finite element method for minimum weight design of structures with lower-bound constraints on the natural frequencies, and upper and lower bounds on the design variables. The design algorithm is essentially an iterative solution of the Kuhn-Tucker optimality criterion. The three most important features of the algorithm are: (1) a small number of design iterations are needed to reach optimal or near-optimal design, (2) structural elements with a wide variety of size-stiffness may be used, the only significant restriction being the exclusion of curved beam and shell elements, and (3) the algorithm will work for multiple as well as single frequency constraints. The design procedure is illustrated with three simple problems.

Kiusalaas, J.↗

TAP 2: A finite element program for thermal analysis of convectively cooled structures

A finite element computer program (TAP 2) for steady-state and transient thermal analyses of convectively cooled structures is presented. The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature-dependent thermal parameters is performed using the Newton-Raphson iteration method. Transient analyses are performed using an implicit Crank-Nicolson time integration scheme with consistent or lumped capacitance matrices as an option. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗

Development of a locally mass flux conservative computer code for calculating 3-D viscous flow in turbomachines

The VANS successive approximation numerical method was extended to the computation of three dimensional, viscous, transonic flows in turbomachines. A cross-sectional computer code, which conserves mass flux at each point of the cross-sectional surface of computation was developed. In the VANS numerical method, the cross-sectional computation follows a blade-to-blade calculation. Numerical calculations were made for an axial annular turbine cascade and a transonic, centrifugal impeller with splitter vanes. The subsonic turbine cascade computation was generated in blade-to-blade surface to evaluate the accuracy of the blade-to-blade mode of marching. Calculated blade pressures at the hub, mid, and tip radii of the cascade agreed with corresponding measurements. The transonic impeller computation was conducted to test the newly developed locally mass flux conservative cross-sectional computer code. Both blade-to-blade and cross sectional modes of calculation were implemented for this problem. A triplet point shock structure was computed in the inducer region of the impeller. In addition, time-averaged shroud static pressures generally agreed with measured shroud pressures. It is concluded that the blade-to-blade computation produces a useful engineering flow field in regions of subsonic relative flow; and cross-sectional computation, with a locally mass flux conservative continuity equation, is required to compute the shock waves in regions of supersonic relative flow.

Walitt, L.↗

The fundamental structure function of oscillator noise models

Continuous-time models of oscillator phase noise x(t) usually have stationary nth differences, for some n. The covariance structure of such a model can be characterized in the time domain by the structure function: D sub n (t;gamma sub 1, gamma sub 2) = E delta (n) sub gamma sub 1 x(s+t) delta(n) sub gamma sub 2 x (s). Although formulas for the special case D sub 2 (0;gamma,gamma) (the Allan variance times 2 gamma(2)) exist for power-law spectral models, certain estimation problems require a more complete knowledge of (0). Exhibited is a much simpler function of one time variable, D(t), from which (0) can easily be obtained from the spectral density by uncomplicated integrations. Believing that D(t) is the simplest function of time that holds the same information as (0), D(t) is called the fundamental structure function. D(t) is computed for several power-law spectral models. Two examples are D(t) = K/t/(3) for random walk FM, D(t) = Kt(2) 1n/t/ for flicker FM. Then, to demonstrate its use, a BASIC program is given that computes means and variances of two Allan variance estimators, one of which incorporates a method of frequency drift estimation and removal.

Greenhall, C. A.↗

Mathematical modeling of the flow field and particle motion in a rotating bioreactor at unit gravity and microgravity

The biotechnology group at NASA Johnson Space Center is developing systems for culturing mammalian cells that stimulate some aspect of microgravity and provide a low shear environment for microgravity-based studies on suspension and anchorage dependent cells. The design of these vessels for culturing cells is based on the need to suspend cells and aggregates of cells and microcarrier beads continually in the culturing medium. The design must also provide sufficient circulation for adequate mass transfer of nutrients to the cells and minimize the total force on the cells. Forces, resulting from sources such as hydrodynamic fluid shear and collisions of cells and walls of the vessels, may damage delicate cells and degrade the formation of three dimensional structures. This study examines one particular design in both unit gravity and microgravity based on two concentric cylinders rotating in the same direction at different speeds to create a Couette flow between them. A numerical simulation for the flow field and the trajectories of particles in the vessel. The flow field for the circulation of the culturing medium is modeled by the Navier-Stokes equations. The forces on a particle are assumed to be drag from the fluid's circulation, buoyancy from the gravitational force and centrifugal force from the rotation of the vessel. The problem requires first solving the system of partial differential equations for the fluid flow by a finite difference method and then solving the system of ordinary differential equations for the trajectories by Gear's stiff method. Results of the study indicate that the trajectories in unit gravity and microgravity are very similar except for small spatial deviations on the fast time scale in unit gravity. The total force per unit cross sectional area on a particle in microgravity, however, is significantly smaller than the corresponding value in unit gravity, which is also smaller than anticipated. Hence, this study indicates that this design for a bioreactor with optimal rates of rotation can provide a good environment for culturing cells in microgravity with adequate circulation and minimal force on the cells.

Boyd, Ernest J.↗

Toward a More Robust Pruning Procedure for MLP Networks

Choosing a proper neural network architecture is a problem of great practical importance. Smaller models mean not only simpler designs but also lower variance for parameter estimation and network prediction. The widespread utilization of neural networks in modeling highlights an issue in human factors. The procedure of building neural models should find an appropriate level of model complexity in a more or less automatic fashion to make it less prone to human subjectivity. In this paper we present a Singular Value Decomposition based node elimination technique and enhanced implementation of the Optimal Brain Surgeon algorithm. Combining both methods creates a powerful pruning engine that can be used for tuning feedforward connectionist models. The performance of the proposed method is demonstrated by adjusting the structure of a multi-input multi-output model used to calibrate a six-component wind tunnel strain gage.

Stepniewski, Slawomir W.↗

Binary-Signal Recovery

Binary communication through long cables, opto-isolators, isolating transformers, or repeaters can become distorted in characteristic ways. The usual solution is to slow the communication rate, change to a different method, or improve the communication media. It would help if the characteristic distortions could be accommodated at the receiving end to ease the communication problem. The distortions come from loss of the high-frequency content, which adds slopes to the transitions from ones to zeroes and zeroes to ones. This weakens the definition of the ones and zeroes in the time domain. The other major distortion is the reduction of low frequency, which causes the voltage that defines the ones or zeroes to drift out of recognizable range. This development describes a method for recovering a binary data stream from a signal that has been subjected to a loss of both higher-frequency content and low-frequency content that is essential to define the difference between ones and zeroes. The method makes use of the frequency structure of the waveform created by the data stream, and then enhances the characteristics related to the data to reconstruct the binary switching pattern. A major issue is simplicity. The approach taken here is to take the first derivative of the signal and then feed it to a hysteresis switch. This is equivalent in practice to using a non-resonant band pass filter feeding a Schmitt trigger. Obviously, the derivative signal needs to be offset to halfway between the thresholds of the hysteresis switch, and amplified so that the derivatives reliably exceed the thresholds. A transition from a zero to a one is the most substantial, fastest plus movement of voltage, and therefore will create the largest plus first derivative pulse. Since the quiet state of the derivative is sitting between the hysteresis thresholds, the plus pulse exceeds the plus threshold, switching the hysteresis switch plus, which re-establishes the data zero to one transition, except now at the logic levels of the receiving circuit. Similarly, a transition from a one to a zero will be the most substantial and fastest minus movement of voltage and therefore will create the largest minus first derivative pulse. The minus pulse exceeds the minus threshold, switching the hysteresis switch minus, which re-establishes the data one to zero transition. This innovation has a large increase in tolerance for the degradation of the binary pattern of ones and zeroes, and can reject the introduction of noise in the form of low frequencies that can cause the voltage pattern to drift up or down, and also higher frequencies that are beyond the recognizable content in the binary transitions.

Griebeler, Elmer L.↗

A Stability Analysis of Cylindrical Panels Using a Finite Element Formulation

A cylindrical finite element suitable for the linear stability analysis of cylindrical shells is developed. Energy principles and variational methods lead to a problem formulation which lends itself to physical interpretations of the governing matrices of the finite element. By properly grouping the terms which result from taking the second variation of the potential energy of the element, it is possible to identify three distinct types of matrices. The first matrix is the conventional stiffness matrix; the second is an initial stress stiffness matrix; and the third is an initial displacement stiffness matrix. With the assumption of linearity, the buckling problem is stated in terms of the classical linear real eigenvalue equation. This problem formulation was programmed on the CDC 6600 series computer. The computer program is used to analyze the buckling of a variety of structures. Columns, arches, flat plates and curved panels with and without cutouts are considered.

Snyder, R. E.↗

Review - Axial compressor stall phenomena

Stall in compressors can be associated with the initiation of several types of fluid dynamic instabilities. These instabilities and the different phenomena, surge and rotating stall, which result from them, are discussed in this paper. Assessment is made of the various methods of predicting the onset of compressor and/or compression system instability, such as empirical correlations, linearized stability analyses, and numerical unsteady flow calculation procedures. Factors which affect the compressor stall point, in particular inlet flow distortion, are reviewed, and the techniques which are used to predict the loss in stall margin due to these factors are described. The influence of rotor casing treatment (grooves) on increasing compressor flow range is examined. Compressor and compression system behavior subsequent to the onset of stall is surveyed, with particular reference to the problem of engine recovery from a stalled condition. The distinction between surge and rotating stall is emphasized because of the very different consequences on recoverability. The structure of the compressor flow field during rotating stall is examined, and the prediction of compressor performance in rotating stall, including stall/unstall hysteresis, is described.

Greitzer, E. M.↗

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).

Rodriguez, G.↗

Stochastic Control Synthesis of Systems with Structured Uncertainty

This paper presents a study on the design of robust controllers by using random variables to model structured uncertainty for both SISO and MIMO feedback systems. Once the parameter uncertainty is prescribed with probability density functions, its effects are propagated through the analysis leading to stochastic metrics for the system's output. Control designs that aim for satisfactory performances while guaranteeing robust closed loop stability are attained by solving constrained non-linear optimization problems in the frequency domain. This approach permits not only to quantify the probability of having unstable and unfavorable responses for a particular control design but also to search for controls while favoring the values of the parameters with higher chance of occurrence. In this manner, robust optimality is achieved while the characteristic conservatism of conventional robust control methods is eliminated. Examples that admit closed form expressions for the probabilistic metrics of the output are used to elucidate the nature of the problem at hand and validate the proposed formulations.

Padula, Sharon L.↗

Design and Analysis Tools for Concurrent Blackboard Systems

A blackboard system consists of a set of knowledge sources, a blackboard data structure, and a control strategy used to activate the knowledge sources. The blackboard model of problem solving is best described by Dr. H. Penny Nii of the Stanford University AI Laboratory: "A Blackboard System can be viewed as a collection of intelligent agents who are gathered around a blackboard, looking at pieces of information written on it, thinking about the current state of the solution, and writing their conclusions on the blackboard as they generate them. " The blackboard is a centralized global data structure, often partitioned in a hierarchical manner, used to represent the problem domain. The blackboard is also used to allow inter-knowledge source communication and acts as a shared memory visible to all of the knowledge sources. A knowledge source is a highly specialized, highly independent process that takes inputs from the blackboard data structure, performs a computation, and places the results of the computation in the blackboard data structure. This design allows for an opportunistic control strategy. The opportunistic problem-solving technique allows a knowledge source to contribute towards the solution of the current problem without knowing which of the other knowledge sources will use the information. The use of opportunistic problem-solving allows the data transfers on the blackboard to determine which processes are active at a given time. Designing and developing blackboard systems is a difficult process. The designer is trying to balance several conflicting goals and achieve a high degree of concurrent knowledge source execution while maintaining both knowledge and semantic consistency on the blackboard. Blackboard systems have not attained their apparent potential because there are no established tools or methods to guide in their construction or analyze their performance.

McManus, John W.↗