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.
Three dimensional motion analysis and vibration characteristics of dynamically scaled model of asymmetrically clustered launch vehicle
A method and an associated digital computer program for evaluating the vibrational characteristics of large liquid-filled rigid wall tanks of general shape are presented. A solution procedure was developed in which slosh modes and frequencies are computed for systems mathematically modeled as assemblages of liquid finite elements. To retain sparsity in the assembled system mass and stiffness matrices, a compressible liquid element formulation was incorporated in the program. The approach taken in the liquid finite element formulation is compatible with triangular and quadrilateral structural finite elements so that the analysis of liquid motion can be coupled with flexible tank wall motion at some future time. The liquid element repertoire developed during the course of this study consists of a two-dimensional triangular element and a three-dimensional tetrahedral element.
Methods are presented for representing unsteady aerodynamic loadings, valid for arbitrary motion, for the two-dimensional typical section with a trailing-edge control surface in incompressible flow and a three-dimensional lifting surface with leading- or trailing-edge control surfaces in subsonic compressible flow. Loads for the two-dimensional incompressible case are obtained via analytic continuation of the Theodorsen function into the complex plane, with guidance from a time-domain approximation of the Wagner indicial function. For the three-dimensional case, it is proposed that oscillatory generalized aerodynamic forces be approximated by Pade fractions, thereby permitting an obvious continuation into the complex plane. A theoretical justification for this procedure is briefly outlined, and some practical aspects of its implementation are discussed.
Particular attention is given to models of the upper atmosphere that deal with the diurnal variation and the latitudinal density dependence. The upper atmosphere is considered as essentially an ideal gas, or rather a weakly ionized plasma, where physical and chemical processes of considerable complexity take place. The final goal for the theoretical models of the upper atmosphere is to derive the time and space dependence of the physical quantities density, composition, temperature, pressure, and state of motion by theoretical considerations from basic physical assumptions. Questions of the dynamics of the upper atmosphere are discussed together with one- and two-dimensional models of the upper atmosphere, three-dimensional models, the geomagnetic effect, and the semiannual density variation.
At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three-dimensional problems.
Three-dimensional unsteady transonic flow through an axial turbomachine stage is described in terms of a pair of two-dimensional formulations pertaining to orthogonal surfaces, namely, a blade-to-blade surface and a hub-to-casing surface. The resulting systems of nonlinear, inviscid, compressible equations of motion are solved by an explicit finite-difference technique. The blade-to-blade program includes the periodic interaction between rotor and stator blade rows. Treatment of the boundary conditions and of the blade slipstream motion by a characteristic type procedure is discussed in detail. Harmonic analysis of the acoustic far field produced by the blade row interaction, including an arbitrary initial transient, is outlined. Results from the blade-to-blade program are compared with experimental measurements of the rotating pressure field at the tip of a high-speed fan. The hub-to-casing program determines circumferentially averaged flow properties on a meridional plane. Blade row interactions are neglected in this formulation, but the force distributions over the entire blade surface for both the rotor and stator are obtained. Results from the hub-to-casing program are compared with a relaxation method solution for a subsonic rotor. Results are also presented for a quiet fan stage which includes transonic flow in both the rotor and stator and a normal shock in the stator.
This investigation explores the effects of chordwise forces and deformations and steady-state deformation due to lift on the static and dynamic aeroelastic stability of a uniform cantilever wing. Results of this analysis are believed to have practical applications for high-performance sailplanes and certain RPV's. The airfoil cross section is assumed to be symmetric and camber bending is neglected. Motions in vertical bending, fore-and-aft bending, and torsion are considered. A differential equation model is developed, which included the nonlinear elastic bending-torsion coupling that accompanies even moderate deflections. A linearized expansion in small time-dependent deflections is made about a steady flight condition. The stability determinant of the linearized system then contains coefficients that depend on steady displacements. Loads derived from two-dimensional incompressible aerodynamic theory are used to obtain the majority of the results, but cases using three-dimensional subsonic compressible theory are also studied. The stability analysis is carried out in terms of the dynamically uncoupled natural modes of vibration of the uniform cantilever.
Vortical motion in the early morning sector of the plasma sheet was detected by means of two-dimensional and three-dimensional plasma measurements obtained by the LASL/MPI analyzers on the ISEE 1 and 2 satellites. The vortices, when present, are manifested as a recurrent or continuing rotation of the bulk flow vector in a plane sometimes moderately inclined with respect to the ecliptic plane. The preferred sense of the rotation is clockwise when viewed from above the ecliptic plane. The vortex rotation period ranges from 5 to 20 min., and several rotations can occur in a relatively uninterrupted sequence. The plasma vortices are estimated to be several earth radii in size. The possibility of a 'vortex street' convecting earthward along the axis of the magnetotail is considered.
Analytical and empirical studies of a finite difference method for the solution of the transonic flow about harmonically oscillating wings and airfoils are presented. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady equations for small disturbances. Since sinusoidal motion is assumed, the unsteady equation is independent of time. Three finite difference investigations are discussed including a new operator for mesh points with supersonic flow, the effects on relaxation solution convergence of adding a viscosity term to the original differential equation, and an alternate and relatively simple downstream boundary condition. A method is developed which uses a finite difference procedure over a limited inner region and an approximate analytical procedure for the remaining outer region. Two investigations concerned with three-dimensional flow are presented. The first is the development of an oblique coordinate system for swept and tapered wings. The second derives the additional terms required to make row relaxation solutions converge when mixed flow is present. A finite span flutter analysis procedure is described using the two-dimensional unsteady transonic program with a full three-dimensional steady velocity potential.
A detailed method of calculating the three-dimensional interaction of gas atoms with solid surfaces has been developed in order to provide a quantitative means of determining atmospheric evolution and global distribution on exospheric planetary bodies. The computational method involves accurate structural models of the surface and subsurface, and numerical solution of 60-dimensional differential equations of motion in a Monte Carlo calculation of the collision process. Preliminary results for the He-alpha-quartz system are described.
A theoretical investigation of three dimensional relativistic klystron action is described. The relativistic axisymmetric equations of motion are derived from the time-dependent Lagrangian function for a charged particle in electromagnetic fields. An analytical expression of the fringing RF electric and magnetic fields within and in the vicinity of the interaction gap and the space-charge forces between axially and radially elastic deformable rings of charges are both included in the formulation. This makes an accurate computation of electron motion through the tunnel of the cavities and the drift tube spaces possible. Method of analysis is based on Lagrangian formulation. Bunching is computed using a disk model of electron stream in which the electron stream is divided into axisymmetric disks of equal charge and each disk is assumed to consist of a number of concentric rings of equal charges. The Individual representative groups of electrons are followed through the interaction gaps and drift tube spaces. Induced currents and voltages in interacting cavities are calculated by invoking the Shockley-Ramo theorem.
A high power multicavity klystron amplifier was designed and a computation package containing all equations and procedures needed is presented. The rigorously derived three dimensional relativistic axisymmetric equations of motion are used to compute the bunched current and the induced RF gap voltage for all interaction cavities except the input and second cavities, where the linear space charge wave theory data are employed in order to reduce the computation time. Both distance step and time step integration methods are used to compute the Fourier coefficients of both the beam current and induced current.
The separation of three dimensional turbulent boundary layers from the lee of flight vehicles at high angles of attack is investigated. The separation results in dominant, large scale, coiled vortex motions that pass along the body in the general direction of the free stream. In all cases of three dimensional flow separation and reattachment, the assumption of continuous vector fields of skin friction lines and external flow streamlines, coupled with simple laws of topology, provides a flow grammar whose elemental constituents are the singular points: the nodes, spiral nodes (foci), and saddles. The phenomenon of three dimensional separation may be construed as either a local or a global event, depending on whether the skin friction line that becomes a line of separation originates at a node or a saddle point.
The two-dimensional steady and time dependent properties of plane Poiseuille and plane Couette flows are analyzed using iterative techniques and full numerical simulation of the Navier-Stokes equations. It is shown that the finite-amplitude two-dimensional states investigated are strongly unstable to very small three-dimensional perturbations. It is also shown, through full numerical simulation, that this explosive secondary instability can explain the subcritical transitions that occur in real flows. Finally, it is shown that the three-dimensional instability can be analyzed by a linear stability analysis of a two-dimensional flow consisting of the basic parallel flow and a steady (or quasi-steady) finite-amplitude two-dimensional cellular motion.
Rapid double-pulse holography was employed to obtain detailed, two-dimensional images of the shock forming during simulated flutter in a transonic flowfield. The experiment comprised a linear cascade of airfoils externally oscillated in torsion and viewed tangentially at the shock surface. Three biconvex airfoils were subjected to harmonic pitching motion about the midchord axis at a frequency of 0.53 while immersed in a Mach 0.81 flow. Failure to produce observable shocks led to use of choked flow with a Mach number near one, of which 50 holograms were taken. The images revealed a narrow shock surface with a spanwise variation in the shock properties. The method is concluded to be useful for examining transonic flowfield shocks in the presence of airfoil flutter.
The separation of three dimensional turbulent boundary layers from the lee of flight vehicles at high angles of attack is investigated. The separation results in dominant, large scale, coiled vortex motions that pass along the body in the general direction of the free stream. In all cases of three dimensional flow separation and reattachment, the assumption of continuous vector fields of skin friction lines and external flow streamlines, coupled with simple laws of topology, provides a flow grammar whose elemental constituents are the singular points: the nodes, spiral nodes (foci), and saddles. The phenomenon of three dimensional separation may be constrained as either a local or a global event, depending on whether the skin friction line that becomes a line of separation originates at a node or a saddle point.
This paper describes the experimental pressure and temperature results obtained when rotating a shaft in an eccentric lucite casing at velocities ranging from 209 to 628 rad/s (2000 to 6000 rpm). The results are presented in terms of three-dimensional plots and contour maps. Photographic evidence is presented to illustrate how the downstream and upstream regions of the cavity develop and evolve into the well-known finger patterns. A comparison between pressure and temperature profiles in air-saturated oil and carbon dioxide-saturated oil is presented; the origin and nature of the gases contained in the cavity are discussed. Three analytical models (Swift-Stieber, separation, Floberg) predicting the formation of the cavitation zone are presented and evaluated in light of the experimental results. A motion-picture supplement depicting these results is available upon request.
A numerical procedure is presented for computing the unsteady transonic flow field about three dimensional swept wings undergoing general time dependent motion. The outer inviscid portion of the flow is assumed to be governed by the modified unsteady transonic small disturbance potential equation which is integrated in the time domain by means of an efficient alternating direction implicit approximate factorization algorithm. Gross dominant effects of the shock boundary layer interaction are accounted for by a simple empirically defined model. Viscous flow regions adjacent to the wing surface and in the trailing wake are described by a set of integral equations appropriate for compressible turbulent shear layers. The two dimensional boundary layer equations are applied quasi-statically stripwise across the span. Coupling with the outer inviscid flow is implemented through use of the displacement thickness concept within the limitations of small disturbance theory. Validity of the assumptions underlying the method is established by comparison with experimental data for the flow about a high aspect ratio transport wing having an advanced airfoil section.