SEARCH · Engineering Papers
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.
Implementation in the programming language Julia of the Restricted Alternating Anderson Picard (AAP) acceleration for (non)-linear fixed point iterations.
In this report we present two approaches to data remapping between one-dimensional meshes implemented with the c++ programming language. Our goal was to test the performance of two search algorithms, linear and binary, and verify the accuracy of our implementations of the two methods. We first introduce the concept of data remap and meshing components, as well as their various uses. We then delve into the differences between point-wise and conservative remap, the algorithms used in the implementations, and lastly confirm the implementations work as intended when given various inputs. We expect that, after profiling, the binary search algorithm will be more efficient than the linear algorithm for sorted sets of data, the point-wise remap implementation to accurately approximate the data transfer between two meshes, and the conservative remap implementation to conserve the area underneath the curve of two distinct meshes.
Analysis of linear elastic thin shells of revolution subjected to arbitrary temperatures and loads - computer program users manual
MTRAC program for transient analysis of circuits with square loop magnetic cores based on linearization of nonlinear circuit elements and Kirchhoff current law
Recommendations for improving the methodology of pogo suppression for the space shuttle include: Consideration of inter-pump location for accumulator or active device, inclusion of tank outflow effects in dynamic structural analysis, the use of simplified transfer functions in systems studies, three phase dynamic testing program for turbopump with development of dynamic flowmeters, and the use of a linearized mathematical model for engine physics studies.
Results are presented of research into the problem of finding an excitation of a given antenna such that the desired radiation pattern is approximated to within acceptable limits. This is to be done in such a fashion that boundary conditions involving hardware limitations may be inserted into the problem. The intended application is synthesis of multiple shaped beam antennas. Since this is perhaps the most difficult synthesis problem an antenna engineer is likely to encounter, the approach taken was to include as a by-product capability for synthesizing simpler patterns. The synthesis technique has been almost totally computerized. The class of antennas which may be synthesized with the computer program are those which may be represented as planar (continuous or discrete) current distributions. The technique is not limited in this sense and could indeed by extended to include, for example, the synthesis of conformal arrays or current distributions on the surface of reflectors. The antenna types which the program is set up to synthesize are: line source, rectangular aperture, circular aperture, linear array, rectangular array, and arbitrary planar array.
This paper discusses numerical aspects of computing maximum likelihood estimates for linear dynamical systems in state-vector form. Different gradient-based nonlinear programming methods are discussed in a unified framework and their applicability to maximum likelihood estimation is examined. The problems due to singular Hessian or singular information matrix that are common in practice are discussed in detail and methods for their solution are proposed. New results on the calculation of state sensitivity functions via reduced order models are given. Several methods for speeding convergence and reducing computation time are also discussed.
Several experiments with J- and blade-stiffened graphite/epoxy panels were conducted to obtain insight into how well experimental data could be correlated with analysis for the buckling behavior of open-section stiffened composite compression panels. Although some nonlinear behavior was observed during the experiments, adequate correlation with analysis was obtained to justify the use of linear, thin-plate buckling analysis in a minimum-weight design synthesis program for J- and blade-configurations. Results from two design studies using this program are presented. In the first study the minimum weights of Jand blade-configurations for two different material systems (graphite/epoxy and aluminum) are determined subject to buckling and strength constraints for a wide range of the compressive load index. In the second study the minimum weights required for graphite/epoxy blade-stiffened panels to satisfy additional stiffness constraints typical of medium-size commercial aircraft wing structures are determined. Both minimum-weight studies indicate that graphite/epoxy open-section stiffened panels can be designed so that weight savings of 30 to 50% are possible compared with the most efficient aluminum designs.
Computerization of aerodynamic theory has progressed to state where analysis of complete aircraft configurations can be performed in single program. Aerodynamic Preliminary Analysis System, APAS, is comprehensive aerodynamic analysis system, based on linearized potential theory. Three-dimensional configurations (with or without jet flaps) having multiple nonplanar surfaces of arbitrary planform and open or closed slender bodies of noncircular contour may be analyzed with APAS. As preliminary design aid, APAS allows designer to survey systematically large number of alternative configurations and component geometries economically.
A computer program and procedure for plotting the contour of the data transferred from the Marshall Space Flight Center solar magnetography is presented. The plotted data then can be easily compared with solar data from other sources, such as the Solar Maximum Mission. From the data file for circular polarization the longitudinal program plots the contours for filtered longitudinal plot and intensity plot by choosing the positive and negative contour levels, intensity levels, and also X,Y plotting ranges which need to be used. In a similar manner the transverse program generates the transverse contour plot, azimuth plot, and intensity plot from the linear polarization data files.
Existing integrated infrared detector array technology is being evaluated under low-background conditions to determine its applicability in orbiting astronomical applications where extended integration times and photometric accuracy are of interest. Preliminary performance results of a 1 x 20 elements InSb CCD array under simulated astronomical conditions are presented. Using the findings of these tests, improved linear- and area-array technology will be developed for use in NASA programs such as the Shuttle Infrared Telescope Facility. For wavelengths less than 30 microns, extrinsic silicon and intrinsic arrays with CCD readout will be evaluated and improved as required, while multiplexed arrays of Ge:Ga for wavelengths in the range 30 to 120 microns will be developed as fundamental understanding of this material improves. Future efforts will include development of improved drive and readout circuitry, and consideration of alternate multiplexing schemes.
An integrated system for the analysis of supersonic configurations consists of an executive driver and eight basic computer programs that build up force coefficients of an selected configuration. System employs modified linearized theory for calculation of surface pressures and employs supersonic-area-rule concepts in combination with linearized theory for calculation of aerodynamic force coefficients.
A general approach is required for describing matter of behavior when the failure is likely to involve growth and coalescence of a large number of fractures. Failures of this kind appear frequently in rapid dynamic processes, particularly in the formation of spall fragments. An approach to formulating constitutive relations that accounts for the opening, shear and growth of an ensemble of cracks is discussed. The approach accounts for plastic flow accompanying fragmentation. The resulting constitutive relations were incorporated into a Lagrangian computer program. A theoretical approach to coalescence is described. The simplest formulation uses a linear Liouville equation, with crack growth limited by the mean free path of cracks, assumed constant. This approach allows for an anisotropic distribution of cracks. An alternative approach in which the decrease of the mean free path with increasing crack size is accounted for, but the crack distribution is assumed isotropic is described. A reduction of the governing Liouville equation to an ordinary differential equation of third order is possible, and the result can be used to determine how mean free path decreases with increasing crack size.
The combined effects of blade torsion and dynamic inflow on the aeroelastic stability of an elastic rotor blade in forward flight are studied. The governing sets of equations of motion (fully nonlinear, linearized, and multiblade equations) used in this study are derived symbolically using a program written in FORTRAN. Stability results are presented for different structural models with and without dynamic inflow. A combination of symbolic and numerical programs at the proper stage in the derivation process makes the obtainment of final stability results an efficient and straightforward procedure.
Interactive controls analysis (INCA) program developed to provide user-friendly environment for design and analysis of linear control systems, primarily feedback control. Designed for use with both small- and large-order systems. Using interactive-graphics capability, INCA user quickly plots root locus, frequency response, or time response of either continuous-time system or sampled-data system. Configuration and parameters easily changed, allowing user to design compensation networks and perform sensitivity analyses in very convenient manner. Written in Pascal and FORTRAN.
Recent developments in the area of nonlinear control theory have shown how coordiante changes in the state and input spaces can be used with nonlinear feedback to transform certain nonlinear ordinary differential equations into equivalent linear equations. These feedback linearization techniques are applied to resolve two problems arising in the control of spacecraft equipped with control moment gyroscopes (CMGs). The first application involves the computation of rate commands for the gimbals that rotate the individual gyroscopes to produce commanded torques on the spacecraft. The second application is to the long-term management of stored momentum in the system of control moment gyroscopes using environmental torques acting on the vehicle. An approach to distributing control effort among a group of redundant actuators is described that uses feedback linearization techniques to parameterize sets of controls which influence a specified subsystem in a desired way. The approach is adapted for use in spacecraft control with double-gimballed gyroscopes to produce an algorithm that avoids problematic gimbal configurations by approximating sets of gimbal rates that drive CMG rotors into desirable configurations. The momentum management problem is stated as a trajectory optimization problem with a nonlinear dynamical constraint. Feedback linearization and collocation are used to transform this problem into an unconstrainted nonlinear program. The approach to trajectory optimization is fast and robust. A number of examples are presented showing applications to the proposed NASA space station.
The neighboring extremal feedback control law is developed for systems with a piecewise linear control for the case where the optimal control is obtained by nonlinear programming techniques. To develop the control perturbation for a given deviation from the nominal path, the second variation is minimized subject to the constraint that the final conditions be satisfied. This process leads to a feedback relationship between the control perturbation and the measured deviation from the nominal state. A simple example, the lunar launch problem, is used to demonstrate the validity of the guidance law. For model errors on the order of 5 percent, the results indicate that 5 percent errors occur in the final conditions.
The neighboring optimal feedback control law is developed for systems with a piecewise linear control for the case where the optimal control is obtained by nonlinear programming techniques. To develop the control perturbation for a given deviation from the nominal path, the second variation is minimized subject to the constraint that the final conditions be satisfied (neighboring suboptimal control). This process leads to a feedback relationship between the control perturbation and the measured deviation from the nominal state. Neighboring suboptimal control is applied to the lunar launch problem. Two approaches, single optimization and multiple optimization for calculating the gains are used, and the gains are tested in a guidance simulation with a mismatch in the acceleration of gravity. Both approaches give acceptable results, but multiple optimization keeps the perturbed path closer to the nominal path.