SEARCH · Search NASA
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.
Early options studied in the satisfaction of LACIE implementation computational demands are described as well as the ultimate selection and development of an array processing solution to the problem. The economic justification, as a function of required LANDSAT analysis, is provided. The suitability of such processors for LACIE and other applications is discussed.
The paper discusses the design of a higher harmonic blade pitch control system for flight testing on an OH-6A helicopter. Alternative designs for both the mechanical and electronic subsystems are also presented. Among the recommendations set forth are: (1) use electronic analog methods instead of FFT software, delegating spectral analysis and self testing to an Electronic Control Unit; (2) use a digital rather than analog computer for increased flexibility in solution processing; and (3) locate HHC actuators in the nonrotating system and separate these actuators from the primary control system. It is concluded that a target weight for a prototype control system would be less than 1% of design gross weight and the production weight of the future HHC system may be 0.5 of the weight of the aircraft.
A modification of an implicit approximate-factorization finite-difference algorithm applied to partial differential equations is presented. This algorithm is applied to the two- and three-dimensional Euler equations in general curvilinear coordinates. The modification transforms the coupled system of equations into an uncoupled diagonal form that requires less computational work. For steady-state applications, the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. The diagonal algorithm reduces the storage requirement of the implicit solution process and therefore has an important effect on the application of implicit finite-difference schemes to vector processors. Results are presented for realistic two-dimensional transonic flow fields about airfoils. Computation costs are reduced to 24-34%.
Results are presented from a quasi-three-dimensional calculation of steady (relative), inviscid, adiabatic, subsonic/shock-free transonic flow on multiple hub-to-shroud stream surfaces through turbomachinery blade rows. The quasi-three-dimensional technique incorporates some three-dimensional effects while retaining much of the simplicity of two-dimensional computational methods. Three typical turbomachinery flowfield calculations are presented including an axial-flow compressor rotor, a turbine stator vane cascade, and a radial-inflow turbine rotor. The calculations were performed using quasi-three-dimensional extensions of existing two-dimensional methods. The current results represent an intermediate step in the complete quasi-three-dimensional solution process. However, the results demonstrate the usefulness of the quasi-three-dimensional technique in complementing and extending the applicability of the two-dimensional methods of turbomachinery flow analysis.
A computer code that generates shock-free transonic compressor cascade shapes while taking into account viscosity effects is developed. The mathematical model for the inviscid flow field is the full potential equation. The Kutta-Joukowski condition is satisfied by varying the free stream angle at downstream infinity. A boundary fitted computational grid of C-type is generated using a sequence of conformal mapping and nonorthogonal coordinate stretching and shearing transformations. The inviscid calculation is performed sequentially on up to four consecutively refined grids thereby accelerating the convergence of the solution process. The full potential equation is solved using a finite area technique and rotated, type-dependent finite differencing. Artificial viscosity of the first order is added in a fully conservative form. Shock-free cascade airfoil shapes are obtained using the fictitious gas concept of Sobieczky and the method of characteristics in the rheograph plane. Viscosity effects are incorporated via a boundary layer displacement thickness. The integral boundary layer code is based on Rotta's turbulence model and assumes transition region of zero length.
An engineering and software specification which was written for a computer program to calculate aeroelastic structural loads including the effects of nonlinear aerodynamics is presented. The procedure used in the program for an iterative aeroelastic solution (PIAS) is to alternately execute two computer codes: one to calculate aerodynamic loads for a specific wing shape, and another to calculate the deflected shape caused by this loading. A significant advantage to the design of PIAS is that the initial aerodynamic module can be replaced with others. The leading edge vortex (LEV) program is used as the aerodynamic module in PIAS. This provides the capability to calculate aeroelastic loads, including the effects of a separation induced leading edge vortex. The finite element method available in ATLAS Integrated structural analysis and design system is used to determine the deflected wing shape for the applied aerodynamics and inertia loads. The data management capabilities in ATLAS are used by the execution control monitor (ECM) of PIAS to control the solution process.
The steady state displacements of a rotating advanced turboprop are computed using the geometrically nonlinear capabilities of COSMIC NASTRAN Rigid Format 4 and MSC NASTRAN Solution 64. A description of the modified Newton-Raphson algorithm used by Solution 64 and the iterative scheme used by Rigid Format 4 is provided. A representative advanced turboprop, SR3, was used for the study. Displacements for SR3 are computed for rotational speeds up to 10,000 rpm. The results show Solution 64 to be superior for computating displacements of flexible rotating structures. This is attributed to its ability to update the displacement dependent centrifugal force during the solution process.
A linear formulation of rotor airframe coupling intended for vibration analysis in airframe structural design is presented. The airframe is represented by a finite element analysis model; the rotor is represented by a general set of linear differential equations with periodic coefficients; and the connections between the rotor and airframe are specified through general linear equations of constraint. Coupling equations are applied to the rotor and airframe equations to produce one set of linear differential equations governing vibrations of the combined rotor airframe system. These equations are solved by the harmonic balance method for the system steady state vibrations. A feature of the solution process is the representation of the airframe in terms of forced responses calculated at the rotor harmonics of interest. A method based on matrix partitioning is worked out for quick recalculations of vibrations in design studies when only relatively few airframe members are varied. All relations are presented in forms suitable for direct computer implementation.
In an attempt to provide a set of accurate standard test problems for computational code developers, a series of inviscid airfoil test cases were chosen by the AGARD Working Group 07, a subpanel of the AGARD Fluid Dynamics Panel. The cases include three different airfoils at transonic to supersonic conditions. A large number of international experts responded in this effort with computations that have been contrasted for accuracy and consistency. This paper is a summary of the authors' contribution to this study. In particular, the important aspects of the solution process that made it possible to obtain the high level of accuracy needed in this study are stressed.
The orbit of the Lageos satellite is well suited as a reference frame for studying the rotation of the earth and the relative motion of points on the earth's crust. The satellite laser measurements can determine the location of a set of tracking stations in an appropriate terrestrial coordinate system. The motion of the earth's rotation axis relative to this system can be studied on the basis of the established tracking station locations. The present investigation is concerned with an analysis of 7.7 years of Lageos laser ranging data. In the first solution considered, the entire data span was used to adjust a single set of station positions simultaneously with orbit and earth rotation parameters. Attention is given to the accuracy of earth rotation parameters which are determined as an inherent part of the solution process.
The first priority of the NSTS program is to make the Space Shuttle system fully operational and cost effective in providing routine access to space. In support of this priority an integrated logistics support system was planned, structured and is being implemented to support a launch-on-time goal of 95 percent. In achieving a 95 percent spares 'fill rate' in an environment of small fleet size, highly unique and high cost assets, significant spares cost can be incurred. A major portion of these costs are for the additional spares required when repair or acquisition times are lengthy. This paper provides a fundamental analysis of the costs and other factors involved in the spare and repair process and provides an optimized cost and process solution for the Space Shuttle program - a centralized repair depot located at KSC.
A general mathematical model and solution methodologies are being developed for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic, or static thermomechanical loads. Among the system responses, which were associated with these load conditions, were thermal buckling, creep buckling, and ratcheting. Thus, geometric as well as material-type nonlinearities (of high order) can be anticipated and must be considered in the development of the mathematical model. Furthermore, this must also be accommodated in the solution process.
An optimal control solution process was developed for a general class of nonlinear dynamical systems. The method combines control theory, perturbation methods, and Van Loan's recent matrix exponential results. A variety of applications support the practical utility of this method. Nonlinear rigid body optimal maneuvers are routinely solved. Flexible body dynamical systems of an order greater than 40 were solved. The method fails occasionally due to poor convergence of the perturbation expansion or numerical difficulties associated with computing the matrix exponential. The method is attractive because it appears to be a good candidate for semi-automation; no initial guess is required, and it usually converges at 2nd or 3rd order in minutes of machine time.
The weathering of planetary surfaces is treated. Both physical and chemical weathering (reactions between minerals or mineraloids and planetary volatiles through oxidation, hydration, carbonation, or solution processes) are discussed. Venus, earth, and Mars all possess permanent atmospheres such that weathering should be expected to significantly affect their respective surfaces. In contrast, Mercury and the moon lack permanent atmospheres but conceivably could experience surface weathering in response to transient atmospheres generated by volcanic or impact cratering events. Weathering processes can be postulated for other rocky objects including Io, Titan, asteroids, and comets.
A three-dimensional, quasi-steady, full potential flow solver was adapted to include realistic rotor wake influence for the aerodynamic analysis of helicopter rotors. The method uses an inner/outer domain technique to accommodate wake effects. Nonlinear flow is computed in the inner domain using a finite difference solution method. The wake is modeled using prescribed wake techniques to allow for the inclusion of realistic wake geometries. Portions of the wake passing inside the inner domain are treated using an embedded vortex technique. The procedure couples the wake influence with the inner domain solution in a consistent and efficient solution process. Correlation with measured lifting transonic data in hover and forward flight is shown which demonstrates the merits of the approach.
A 3-D, quasi-steady, full potential flow solver was adapted to include realistic wake influence for the aerodynamic analysis of helicopter rotors. The method is based on a finite difference solution of the full potential equation, using an inner and outer domain procedure for the blade flowfield to accommodate wake effects. The nonlinear flow is computed in the inner domain region using a finite difference solution method. The wake is modeled by a vortex lattice using prescribed geometry techniques to allow for the inclusion of realistic rotor wakes. The key feature of the analysis is that vortices contained within the finite difference mesh (inner domain) were treated with a vortex embedding technique while the influence of the remaining portion of the wake (in the outer domain) is impressed as a boundary condition on the outer surface of the finite difference mesh. The solution procedure couples the wake influence with the inner domain solution in a consistent and efficient solution process. The method has been applied to both hover and forward flight conditions. Correlation with subsonic and transonic hover airload data is shown which demonstrates the merits of the approach.
An analytic technique for accounting for the joint effects of Earth oblateness and atmospheric drag on close-Earth satellites is investigated. The technique is analytic in the sense that explicit solutions to the Lagrange planetary equations are given; consequently, no numerical integrations are required in the solution process. The atmospheric density in the technique described is represented by a rotating spherical exponential model with superposed effects of the oblate atmosphere and the diurnal variations. A computer program implementing the process is discussed and sample output is compared with output from program NSEP (Numerical Satellite Ephemeris Program). NSEP uses a numerical integration technique to account for atmospheric drag effects.
A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.