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.
Report presents a new method for solving linear equations developed by Bromwich, which is suited to determine the motion for any particular gust, when the machine started from equilibrium, without the trouble of determining the constants of integration in the complementary function. (author)
A method is derived for calculating the damping coefficients in pitch and roll for a series of triangular wings and a restricted series of sweptback wings at supersonic speeds. The elementary "supersonic source" solution of the linearized equation of motion is used to find the potential function of a line of doublets, and the flows are obtained by surface distributions of these doublet lines. The damping derivatives for triangular wings are found to be a function of the ratio of the tangent of the apex angle to the tangent of the Mach angle. As this ratio becomes equal to and greater than 1.0 for triangular wings, the damping derivatives, in pitch and in roll, become constant. The damping derivative in roll becomes equal to one-half the value calculated for an infinite rectangular wing, and the damping derivative in pitch for pitching about the apex becomes equal to 3.375 times that of an infinite rectangular wing.
A graphite/bismaleimide laminate was prepared without the usual fiber treatment and was tested over a wide range of stress states to measure its ply cracking strength. These tests were conducted using off-axis flexure specimens and produced fiber-matrix interface failure data over a correspondingly wide range of interface stress states. The absence of fiber treatment, weakened the fiber-matrix interfaces and allowed these tests to be conducted at load levels that did not yield the matrix. An elastic micromechanics computer code was used to calculate the fiber-matrix interface stresses at failure. Two different fiber-array models (square and diamond) were used in these calculations to analyze the effects of fiber arrangement as well as stress state on the critical interface stresses at failure. This study showed that both fiber-array models were needed to analyze interface stresses over the range of stress states. A linear equation provided a close fit to these critical stress combinations and, thereby, provided a fiber-matrix interface failure criterion. These results suggest that prediction procedures for laminate ply cracking can be based on micromechanics stress analyses and appropriate fiber-matrix interface failure criteria. However, typical structural laminates may require elastoplastic stress analysis procedures that account for matrix yielding, especially for shear-dominated ply stress states.
In this report, we provide some examples of French, Russian, Chinese, and Japanese launch vehicles that have utilized fins in their designs. Next, the aerodynamic design of the fins is considered in Section 3. Some comments on basic static stability and control theory are followed by a brief description of an aerodynamic characteristics prediction code that was used to estimate the characteristics of a modified NLS 1.5 Stage vehicle. Alternative fin designs are proposed and some estimated aerodynamic characteristics presented and discussed. Also included in Section 3 is a discussion of possible methods of enhancement of the aerodynamic efficiency of fins, such as vortex generators and jet flaps. We consider the construction of fins for launch vehicles in Section 4 and offer an assessment of the state-of-the-art in the use of composites for aerodynamic control surfaces on high speed vehicles. We also comment on the use of smart materials for launch vehicle fins. The dynamic stability and control of a launch vehicle that utilizes both thrust vector control (engine nozzle gimballing) and movable fins is the subject addressed in Section 5. We give a short derivation of equations of motion for a launch vehicle moving in a vertical plane above a spherical earth, discuss the use of a gravity-turn nominal trajectory, and give the form of the period equations linearized about such a nominal. We then consider feedback control of vehicle attitude using both engine gimballing and fin deflection. Conclusions are stated and recommendations made in Section 6. An appendix contains aerodynamic data in tabular and graphical formats.
A model for signal delay computation in RC networks is presented. The strength of the paradigm is its generality and simplicity. The definition of delay is applicable to RC meshes with potential resistive attenuating paths to ground. The algorithms can also be applied to undriven circuits (static charge sharing) and circuits with initial charge. To compute the delays, each node in the network is explored locally to derive a system of sparse linear equations. The solutions of the system are delay values based on the Elmore time constant at each point in the circuit.
A dense linear system solver package recently developed at the University of Texas at Austin for distributed memory machine (e.g. Intel Paragon) has been reviewed and analyzed. The package contains about 45 software routines, some written in FORTRAN, and some in C-language, and forms the basis for parallel/distributed solutions of systems of linear equations encountered in many problems of scientific and engineering nature. The package, being studied by the Computer Applications Branch of the Analysis and Computation Division, may provide a significant computational resource for NASA scientists and engineers in parallel/distributed computing. Since the package is new and not well tested or documented, many of its underlying concepts and implementations were unclear; our task was to review, analyze, and critique the package as a step in the process that will enable scientists and engineers to apply it to the solution of their problems. All routines in the package were reviewed and analyzed. Underlying theory or concepts which exist in the form of published papers or technical reports, or memos, were either obtained from the author, or from the scientific literature; and general algorithms, explanations, examples, and critiques have been provided to explain the workings of these programs. Wherever the things were still unclear, communications were made with the developer (author), either by telephone or by electronic mail, to understand the workings of the routines. Whenever possible, tests were made to verify the concepts and logic employed in their implementations. A detailed report is being separately documented to explain the workings of these routines.
A mathematical model containing the essential features embodied in the noise suppression of lined ejectors is presented. Although some simplification of the physics is necessary to render the model mathematically tractable, the current model is the most versatile and technologically advanced at the current time. A system of linearized equations and the boundary conditions governing the sound field are derived starting from the equations of fluid dynamics. A nonreflecting boundary condition is developed. In view of the complex nature of the equations, a parametric study requires the use of numerical techniques and modern computers. A finite element algorithm that solves the differential equations coupled with the boundary condition is then introduced. The numerical method results in a matrix equation with several hundred thousand degrees of freedom that is solved efficiently on a supercomputer. The model is validated by comparing results either with exact solutions or with approximate solutions from other works. In each case, excellent correlations are obtained. The usefulness of the model as an optimization tool and the importance of variable impedance liners as a mechanism for achieving broadband suppression within a lined ejector are demonstrated.
The stability has been investigated of the unsteady flow past an infinite flat plate when it is moved impulsively from rest, in its own plane. For small times the instantaneous stability of the flow depends on the linearized equations of motion which reduce in this problem to the Orr-Sommerfeld equation. It is known that the flow for certain values of Reynolds number, frequency and wave number is unstable to Tollmien-Schlichting waves, as in the case of the Blasius boundary layer flow past a flat plate. With increase in time, the unstable waves only undergo growth for a finite time interval, and this growth rate is itself a function of time. The influence of finite amplitude effects is studied by solving the full Navier-Stokes equations. It is found that the stability characteristics are markedly changed both by the consideration of the time evolution of the flow, and by the introduction of finite amplitude effects.
The scope of this project dealt with the investigation of the requirements to support distributed computing of scientific computations over a cluster of cooperative workstations. Various experiments on computations for the solution of simultaneous linear equations were performed in the early phase of the project to gain experience in the general nature and requirements of scientific applications. A specification of a distributed integrated computing environment, DICE, based on a distributed shared memory communication paradigm has been developed and evaluated. The distributed shared memory model facilitates porting existing parallel algorithms that have been designed for shared memory multiprocessor systems to the new environment. The potential of this new environment is to provide supercomputing capability through the utilization of the aggregate power of workstations cooperating in a cluster interconnected via a local area network. Workstations, generally, do not have the computing power to tackle complex scientific applications, making them primarily useful for visualization, data reduction, and filtering as far as complex scientific applications are concerned. There is a tremendous amount of computing power that is left unused in a network of workstations. Very often a workstation is simply sitting idle on a desk. A set of tools can be developed to take advantage of this potential computing power to create a platform suitable for large scientific computations. The integration of several workstations into a logical cluster of distributed, cooperative, computing stations presents an alternative to shared memory multiprocessor systems. In this project we designed and evaluated such a system.
A coupled multiple-method analysis procedure for use late in conceptual design or early in preliminary design of aircraft structures is described. Using this method, aircraft wing structures are represented with equivalent plate models, and structural details such as engine/pylon structure, landing gear, or a 'stick' model of a fuselage are represented with beam finite element models. These two analysis methods are implemented in an integrated multiple-method formulation that involves the assembly and solution of a combined set of linear equations. The corresponding solution vector contains coefficients of the polynomials that describe the deflection of the wing and also the components of translations and rotations at the joints of the beam members. Two alternative approaches for coupling the methods are investigated; one using transition finite elements and the other using Lagrange multipliers. The coupled formulation is applied to the static analysis and vibration analysis of a conceptual design model of a fighter aircraft. The results from the coupled method are compared with corresponding results from an analysis in which the entire model is composed of finite elements.
The stability of ion distributions with velocity space holes caused by chaotic ion acceleration in the vicinity of the field reversal region in the Earth's magnetotail is investigated. It is shown that such distributions are unstable to electrostatic waves with wave vector directed along the Sun-Earth line. The solution of the quasi-linear equation shows that the saturation of the instability occurs due to velocity space trapping of the low-energy background ions and results in a flattening of the reduced ion distribution function in the region of the original velocity space hole. The energy level of the excited waves in the vicinity of the Earth's midplane is estimated as a function of the initial size of the empty velocity space region.
A discrete sensitivity analysis algorithm had previously been developed and applied to two-dimensional aerodynamic optimization problems, where the computational domains were discretized by using single grids. The sparse, unsymmetric systems of linear equations resulting from this algorithm were solved by a direct matrix inversion matrix. However, for large two-dimensional problems and, practically, all three-dimensional problems, direct inversion methods become inapplicable, primarily due to the prohibitive computer storage needed. In an attempt to alleviate such hindrances, the sensitivity analysis with domain decomposition (SADD) scheme was developed. This scheme divides the computational domain into smaller and nonoverlapping subdomains (multiblock grids) that are solved separately. Then, the final solution is constructed from the subdomain solutions. As the number of grid points in the interface boundaries of the subdomains becomes large, the computer memory required to store the effective coefficient matrix of these interface points starts to increase. Presented in this Technical Note is a preconditioned iterative procedure to overcome this particular problem.
Three methods to optimize rotorcraft aeromechanical behavior for those cases where the rotorcraft plant can be adequately represented by a linear model system matrix were identified and implemented in a stand-alone code. These methods determine the optimal control vector which minimizes the vibration metric subject to constraints at discrete time points, and differ from the commonly used non-optimal constraint penalty methods such as those employed by conventional controllers in that the constraints are handled as actual constraints to an optimization problem rather than as just additional terms in the performance index. The first method is to use a Non-linear Programming algorithm to solve the problem directly. The second method is to solve the full set of non-linear equations which define the necessary conditions for optimality. The third method is to solve each of the possible reduced sets of equations defining the necessary conditions for optimality when the constraints are pre-selected to be either active or inactive, and then to simply select the best solution. The effects of maneuvers and aeroelasticity on the systems matrix are modelled by using a pseudo-random pseudo-row-dependency scheme to define the systems matrix. Cases run to date indicate that the first method of solution is reliable, robust, and easiest to use, and that it was superior to the conventional controllers which were considered.
The stability calculations made for a shock-free supersonic jet using the model based on parabolized stability equations are presented. In this analysis the large scale structures, which play a dominant role in the mixing as well as the noise radiated, are modeled as instability waves. This model takes into consideration non-parallel flow effects and also nonlinear interaction of the instability waves. The stability calculations have been performed for different frequencies and mode numbers over a range of jet operating temperatures. Comparisons are made, where appropriate, with the solutions to Rayleigh's equation (linear, inviscid analysis with the assumption of parallel flow). The comparison of the solutions obtained using the two approaches show very good agreement.
Multigrid methods are optimal; that is, their rate of convergence is independent of the number of grid points, because they use a nested sequence of coarse grids to represent different scales of the solution. This nesting does, however, usually lead to certain restrictions of the permissible size of the discretised problem. In cases where the modeler is free to specify the whole problem, such constraints are of little importance because they can be taken into consideration from the outset. We consider the situation in which there are other competing constraints on the resolution. These restrictions may stem from the physical problem (e.g., if the discretised operator contains experimental data measured on a fixed grid) or from the need to avoid limitations set by the hardware. In this paper we discuss a modification to the cell-centered multigrid algorithm, so that it can be used br problems with any resolution. We discuss in particular a coarsening strategy and choice of intergrid transfer operators that can handle grids with both an even or odd number of cells. The method is described and applied to linear equations obtained by discretization of two- and three-dimensional second-order elliptic PDEs.
This paper presents stability calculations made for a shock-free supersonic jet using the model based on parabolized stability equations. In this analysis the large-scale structures, which play a dominant role in the mixing as well as the noise radiated, are modeled as instability waves. This model takes into consideration non-parallel flow effects and also nonlinear interaction of the instability waves. The stability calculations have been performed for different frequencies and mode numbers over a range of jet operating temperatures. Comparisons are made, where appropriate, with the solutions to Rayleigh's equation (linear, inviscid analysis with the assumption of parallel flow). The comparison of the solutions obtained using the two approaches show very good agreement.
In this paper we describe a block Lanczos algorithm that is used as the key building block of a software package for the extraction of eigenvalues and eigenvectors of large sparse symmetric generalized eigenproblems. The software package comprises: a version of the block Lanczos algorithm specialized for spectrally transformed eigenproblems; an adaptive strategy for choosing shifts, and efficient codes for factoring large sparse symmetric indefinite matrices. This paper describes the algorithmic details of our block Lanczos recurrence. This uses a novel combination of block generalizations of several features that have only been investigated independently in the past. In particular new forms of partial reorthogonalization, selective reorthogonalization and local reorthogonalization are used, as is a new algorithm for obtaining the M-orthogonal factorization of a matrix. The heuristic shifting strategy, the integration with sparse linear equation solvers and numerical experience with the code are described in a companion paper.
This report characterizes the sets of small disturbances possible in cylindrical and annular ducts with mean flow whose axial and tangential components vary arbitrarily with radius. The linearized equations of motion are presented and discussed, and then exponential forms for the axial, circumferential, and time dependencies of any unsteady disturbances are assumed. The resultant equations form a generalized eigenvalue problem, the solution of which yields the axial wavenumbers and radial mode shapes of the unsteady disturbances. Two numerical discretizations are applied to the system of equations: (1) a spectral collocation technique based on Chebyshev polynomial expansions on the Gauss-Lobatto points, and (2) second and fourth order finite differences on uniform grids. The discretized equations are solved using a standard eigensystem package employing the QR algorithm. The eigenvalues fall into two primary categories: a discrete set (analogous to the acoustic modes found in uniform mean flows) and a continuous band (analogous to convected disturbances in uniform mean flows) where the phase velocities of the disturbances correspond to the local mean flow velocities. Sample mode shapes and eigensystem distributions are presented for both sheared axial and swirling flows. The physics of swirling flows is examined with reference to hydrodynamic stability and completeness of the eigensystem expansions. The effect of assuming exponential dependence in the axial direction is discussed.