Search NASASearch

Engineering topics

Dowell, Earl H.

Publications and source records attributed to Dowell, Earl H..

At least 19 records

Recent Advances in Heliogyro Solar Sail Structural Dynamics, Stability, and Control Research

Results from recent NASA sponsored research on the structural dynamics, stability, and control characteristics of heliogyro solar sails are summarized. Specific areas under investigation include coupled nonlinear finite element analysis of heliogyro membrane blade with solar radiation pressure effects, system identification of spinning membrane structures, and solarelastic stability analysis of heliogyro solar sails, including stability during blade deployment. Recent results from terrestrial 1-g blade dynamics and control experiments on "rope ladder" membrane blade analogs, and small-scale in vacuo system identification experiments with hanging and spinning high-aspect ratio membranes will also be presented. A low-cost, rideshare payload heliogyro technology demonstration mission concept is used as a mission context for these heliogyro structural dynamics and solarelasticity investigations, and is also described. Blade torsional dynamic response and control are also shown to be significantly improved through the use of edge stiffening structural features or inclusion of modest tip masses to increase centrifugal stiffening of the blade structure. An output-only system identification procedure suitable for on-orbit blade dynamics investigations is also developed and validated using ground tests of spinning sub-scale heliogyro blade models. Overall, analytical and experimental investigations to date indicate no intractable stability or control issues for the heliogyro solar sail concept.

Wilkie, W. Keats

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part II: Nonlinear Theory and Extended Aerodynamics

Conical shell theory and a supersonic potential flow aerodynamic theory are used to study the nonlinear pressure buckling and aeroelastic limit cycle behavior of the thermal protection system for NASA's Hypersonic Inflatable Aerodynamic Decelerator. The structural model of the thermal protection system consists of an orthotropic conical shell of the Donnell type, resting on several circumferential elastic supports. Classical Piston Theory is used initially for the aerodynamic pressure, but was found to be insufficient at low supersonic Mach numbers. Transform methods are applied to the convected wave equation for potential flow, and a time-dependent aerodynamic pressure correction factor is obtained. The Lagrangian of the shell system is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the governing differential-algebraic equations of motion. Aeroelastic limit cycle oscillations and buckling deformations are calculated in the time domain using a Runge-Kutta method in MATLAB. Three conical shell geometries were considered in the present analysis: a 3-meter diameter 70 deg. cone, a 3.7-meter 70 deg. cone, and a 6-meter diameter 70 deg. cone. The 6-meter configuration was loaded statically and the results were compared with an experimental load test of a 6-meter HIAD. Though agreement between theoretical and experimental strains was poor, the circumferential wrinkling phenomena observed during the experiments was captured by the theory and axial deformations were qualitatively similar in shape. With Piston Theory aerodynamics, the nonlinear flutter dynamic pressures of the 3-meter configuration were in agreement with the values calculated using linear theory, and the limit cycle amplitudes were generally on the order of the shell thickness. The effect of axial tension was studied for this configuration, and increasing tension was found to decrease the limit cycle amplitudes when the circumferential elastic supports were neglected, but resulted in more complex behavior when the supports were included. The nominal flutter dynamic pressure of the 3.7-meter configuration was significantly lower than that of the 3-meter, and it was found that two sets of natural modes coalesce to flutter modes near the same dynamic pressure. This resulted in a significant drop in the limit cycle frequencies at higher dynamic pressures, where the flutter mode with the lower frequency becomes more critical. Pre-buckling pressure loads and the aerodynamic pressure correction factor were studied for all geometries, and these effects resulted in significantly lower flutter boundaries compared with Piston Theory alone. The maximum dynamic pressure predicted by aerodynamic simulations of a proposed 3.7-meter HIAD vehicle was still lower than any of the calculated flutter dynamic pressures, suggesting that aeroelastic effects for this vehicle are of little concern.

Goldman, Benjamin D.

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part I: Linear Theory

Conical shell theory and piston theory aerodynamics are used to study the aeroelastic stability of the thermal protection system (TPS) on the NASA Hypersonic Inflatable Aerodynamic Decelerator (HIAD). Structural models of the TPS consist of single or multiple orthotropic conical shell systems resting on several circumferential linear elastic supports. The shells in each model may have pinned (simply-supported) or elastically-supported edges. The Lagrangian is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the equations of motion. The natural modes of vibration and aeroelastic stability boundaries are found by calculating the eigenvalues and eigenvectors of a large coefficient matrix. When the in-flight configuration of the TPS is approximated as a single shell without elastic supports, asymmetric flutter in many circumferential waves is observed. When the elastic supports are included, the shell flutters symmetrically in zero circumferential waves. Structural damping is found to be important in this case. Aeroelastic models that consider the individual TPS layers as separate shells tend to flutter asymmetrically at high dynamic pressures relative to the single shell models. Several parameter studies also examine the effects of tension, orthotropicity, and elastic support stiffness.

Goldman, Benjamin D.

Heliogyro Solar Sail Research at NASA

The recent successful flight of the JAXA IKAROS solar sail has renewed interest within NASA in spinning solar sail concepts for high-performance solar sailing. The heliogyro solar sail, in particular, is being re-examined as a potential game-changing architecture for future solar sailing missions. In this paper, we present an overview of ongoing heliogyro technology development and feasibility assessment activities within NASA. In particular, a small-scale heliogyro solar sail technology demonstration concept will be described. We will also discuss ongoing analytical and experimental heliogyro structural dynamics and controls investigations and provide an outline of future heliogyro development work directed toward enabling a low cost heliogyro technology demonstration mission ca. 2020.

Wilkie, W. Keats

Flutter Analysis of the Thermal Protection Layer on the NASA HIAD

A combination of classical plate theory and a supersonic aerodynamic model is used to study the aeroelastic flutter behavior of a proposed thermal protection system (TPS) for the NASA HIAD. The analysis pertains to the rectangular configurations currently being tested in a NASA wind-tunnel facility, and may explain why oscillations of the articles could be observed. An analysis using a linear flat plate model indicated that flutter was possible well within the supersonic flow regime of the wind tunnel tests. A more complex nonlinear analysis of the TPS, taking into account any material curvature present due to the restraint system or substructure, indicated that significantly greater aerodynamic forcing is required for the onset of flutter. Chaotic and periodic limit cycle oscillations (LCOs) of the TPS are possible depending on how the curvature is imposed. When the pressure from the base substructure on the bottom of the TPS is used as the source of curvature, the flutter boundary increases rapidly and chaotic behavior is eliminated.

Goldman, Benjamin D.

Aeroelastic Studies of a Rectangular Wing with a Hole: Correlation of Theory and Experiment

Two rectangular wing models with a hole have been designed and tested in the Duke University wind tunnel to better understand the effects of damage. A rectangular hole is used to simulate damage. The wing with a hole is modeled structurally as a thin elastic plate using the finite element method. The unsteady aerodynamics of the plate-like wing with a hole is modeled using the doublet lattice method. The aeroelastic equations of motion are derived using Lagrange's equation. The flutter boundary is found using the V-g method. The hole's location effects the wing's mass, stiffness, aerodynamics and therefore the aeroelastic behavior. Linear theoretical models were shown to be capable of predicting the critical flutter velocity and frequency as verified by wind tunnel tests.

Conyers, Howard J.

Aerodynamic and Aeroelastic Insights using Eigenanalysis

This paper presents novel analytical results for eigenvalues and eigenvectors produced using discrete time aerodynamic and aeroelastic models. An unsteady, incompressible vortex lattice aerodynamic model is formulated in discrete time; the importance of several modeling parameters is examined. A detailed study is made of the behavior of the aerodynamic eigenvalues both in discrete and continuous time. The aerodynamic model is then incorporated into aeroelastic equations of motion. Eigenanalyses of the coupled equations produce stability results and modal characteristics which are valid for critical and non-critical velocities. Insight into the modeling and physics associated with aeroelastic system behavior is gained by examining both the eigenvalues and the eigenvectors. Potential pitfalls in discrete time model construction and analysis are examined.

Heeg, Jennifer

System Identification and POD Method Applied to Unsteady Aerodynamics

The representation of unsteady aerodynamic flow fields in terms of global aerodynamic modes has proven to be a useful method for reducing the size of the aerodynamic model over those representations that use local variables at discrete grid points in the flow field. Eigenmodes and Proper Orthogonal Decomposition (POD) modes have been used for this purpose with good effect. This suggests that system identification models may also be used to represent the aerodynamic flow field. Implicit in the use of a systems identification technique is the notion that a relative small state space model can be useful in describing a dynamical system. The POD model is first used to show that indeed a reduced order model can be obtained from a much larger numerical aerodynamical model (the vortex lattice method is used for illustrative purposes) and the results from the POD and the system identification methods are then compared. For the example considered, the two methods are shown to give comparable results in terms of accuracy and reduced model size. The advantages and limitations of each approach are briefly discussed. Both appear promising and complementary in their characteristics.

Tang, Deman

System Identification of a Vortex Lattice Aerodynamic Model

The state-space presentation of an aerodynamic vortex model is considered from a classical and system identification perspective. Using an aerodynamic vortex model as a numerical simulator of a wing tunnel experiment, both full state and limited state data or measurements are considered. Two possible approaches for system identification are presented and modal controllability and observability are also considered. The theory then is applied to the system identification of a flow over an aerodynamic delta wing and typical results are presented.

Juang, Jer-Nan

Characterizing Aeroelastic Systems Using Eigenanalysis, Explicitly Retaining The Aerodynamic Degrees of Freedom

Discrete time aeroelastic models with explicitly retained aerodynamic modes have been generated employing a time marching vortex lattice aerodynamic model. This paper presents analytical results from eigenanalysis of these models. The potential of these models to calculate the behavior of modes that represent damped system motion (noncritical modes) in addition to the simple harmonic modes is explored. A typical section with only structural freedom in pitch is examined. The eigenvalues are examined and compared to experimental data. Issues regarding the convergence of the solution with regard to refining the aerodynamic discretization are investigated. Eigenvector behavior is examined; the eigenvector associated with a particular eigenvalue can be viewed as the set of modal participation factors for that particular mode. For the present formulation of the equations of motion, the vorticity for each aerodynamic element appears explicitly as an element of each eigenvector in addition to the structural dynamic generalized coordinates. Thus, modal participation of the aerodynamic degrees of freedom can be assessed in M addition to participation of structural degrees of freedom.

Heeg, Jennifer

Aerodynamic and Aeroelastic Insights using Eigenanalysis

This paper presents novel analytical results for eigenvalues and eigenvectors produced using discrete time aerodynamic and aeroelastic models. An unsteady, incompressible vortex lattice aerodynamic model is formulated in discrete time; the importance of several modeling parameters is examined. A detailed study is made of the behavior of the aerodynamic eigenvalues both in discrete and continuous time. The aerodynamic model is then incorporated into aeroelastic equations of motion. Eigenanalyses of the coupled equations produce stability results and modal characteristics which are valid for critical and non-critical velocities. Insight into the modeling and physics associated with aeroelastic system behavior is gained by examining both the eigenvalues and the eigenvectors. Potential pitfalls in discrete time model construction and analysis are examined.

Heeg, Jennifer

Reduced-Order Aerodynamic Model and Its Application to a Nonlinear Aeroelastic System

Starting from a finite state model for a two-dimensional aerodynamic flow over an airfoil, the eigenmodes of the aerodynamic flow are determined. Using a small number of these aerodynamic eigenmodes, ie., a reduced-order model, the aeroelastic model is formed by coupling them to a typical section structural model with a trailing-edge flap. A free-play nonlinearity is modeled. Results are shown from the finite state model, the reduced-order model, and previous theoretical and experimental work. All results are in good agreement.

Tang, Deman

An efficient procedure for cascade aeroelastic stability determination using nonlinear, time-marching aerodynamic solvers

A numerical eigenvalue problem formulation and a practical calculation procedure for exact eigenvalues and corresponding eigenvectors are developed and applied to a nonlinear, two-dimensional, time-marching full potential solver for cascade aeroelastic stability analysis. This procedure is based on the Lanczos recursive method and it directly calculates stability information about a nonlinear steady state. It is compared to conventional approaches in the frequency and time domains developed earlier and is found to be 100-10.000 times more computationally efficient. Eigenvalue constellations and the flutter results for flow through a cascade SR5 propfan airfoil are presented.

Mahajan, Aparajit J.

Asymptotic modal analysis and statistical energy analysis

Asymptotic Modal Analysis (AMA) is a method which is used to model linear dynamical systems with many participating modes. The AMA method was originally developed to show the relationship between statistical energy analysis (SEA) and classical modal analysis (CMA). In the limit of a large number of modes of a vibrating system, the classical modal analysis result can be shown to be equivalent to the statistical energy analysis result. As the CMA result evolves into the SEA result, a number of systematic assumptions are made. Most of these assumptions are based upon the supposition that the number of modes approaches infinity. It is for this reason that the term 'asymptotic' is used. AMA is the asymptotic result of taking the limit of CMA as the number of modes approaches infinity. AMA refers to any of the intermediate results between CMA and SEA, as well as the SEA result which is derived from CMA. The main advantage of the AMA method is that individual modal characteristics are not required in the model or computations. By contrast, CMA requires that each modal parameter be evaluated at each frequency. In the latter, contributions from each mode are computed and the final answer is obtained by summing over all the modes in the particular band of interest. AMA evaluates modal parameters only at their center frequency and does not sum the individual contributions from each mode in order to obtain a final result. The method is similar to SEA in this respect. However, SEA is only capable of obtaining spatial averages or means, as it is a statistical method. Since AMA is systematically derived from CMA, it can obtain local spatial information as well.

Dowell, Earl H.

Asymptotic modal analysis of a rectangular acoustic cavity excited by wall vibration

Asymptotic modal analysis, a method that has recently been developed for structural dynamical systems, has been applied to a rectangular acoustic cavity. The cavity had a flexible vibrating portion on one wall, and the other five walls were rigid. Banded white noise was transmitted through the flexible portion (plate) only. Both the location along the wall and the size of the plate were varied. The mean square pressure levels of the cavity interior were computed as a ratio of the result obtained from classical modal analysis to that obtained from asymptotic modal analysis for the various plate configurations. In general, this ratio converged to 1.0 as the number of responding modes increased. Intensification effects were found due to both the excitation location and the response location. The asymptotic modal analysis method was both efficient and accurate in solving the given problem. The method has advantages over the traditional methods that are used for solving dynamics problems with a large number of responding modes.

Peretti, Linda F.

Study of intensification zones in a rectangular acoustic cavity

The interior acoustic field of a rectangular acoustic cavity, which is excited by the structural vibration of one of its walls, or a portion of the wall, has been studied. Particularly, the spatial variations of sound pressure levels from the peak levels at the boundaries (intensification zones) to the uniform interior are considered. Analytical expressions, which describe the intensification zones, are obtained using the methodology of asymptotic modal analysis. These results agree well with results computed by a discrete summation over all of the modes. The intensification zones were also modeled as a set of oblique waves incident upon a surface. The result for a rigid surface agrees with the asymptotic modal analysis result. In the presence of an absorptive surface, the character of the intensification zone is dramatically changed. The behavior of the acoustic field near an absorptive wall is described by an expression containing the rigid wall result plus additional terms containing impedance information. The important parameter in the intensification zone analysis is the bandwidth to center frequency ratio. The effect of bandwidth is separated from that of center frequency by expanding the expression about the center frequency wave number. The contribution from the bandwidth is second order in bandwidth to center frequency ratio.

Peretti, Linda F.

Eigenvalue calculation procedure for an Euler/Navier-Stokes solver with application to flows over airfoils

A Lanczos procedure is presently applied to a Navier-Stokes (N-S) solver for eigenvalues and eigenvectors associated with the small-perturbation analysis of the N-S equations' finite-difference representation for airfoil flows; the matrix used is very large, sparse, real, and nonsymmetric. The Lanczos procedure is shown to furnish complete spectral information for the eigenvalues, as required for transient-stability analysis of N-S solvers.

Mahajan, Aparajit J.

Role of artificial viscosity in Euler and Navier-Stokes solvers

A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small time-dependent perturbation about a steady inviscid flow over an NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 deg. The method has been applied to inviscid flows here, but as discussed is also applicable to viscous flows. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundaries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time-marching scheme are also presented. The eigenvalue procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of artificial viscosity.

Mahajan, Aparajit J.