Digital Commons Network^™

Runge-Kutta Upwind Multigrid Multi-Block Three-Dimensional Thin Layer Navier-Stokes Solver, Frank E. Cannizzaro

Mechanical & Aerospace Engineering Theses & Dissertations

A state-of-the-art computer code has been developed that incorporates a modified Runge-Kutta time integration scheme, Upwind numerical techniques, Multigrid acceleration, and Multi-block capabilities (RUMM). A three-dimensional thin-layer formulation of the Navier-Stokes equations is employed. For turbulent flow cases, the Baldwin-Lomax algebraic turbulence model is used. Two different upwind techniques are available, van Leer's flux-vector splitting and Roe's flux-difference splitting. Full approximation multigrid plus implicit residual and corrector smoothing were implemented to enhance the rate of convergence. Multi-block capabilities were developed to provide geometric flexibility. This feature allows the developed computer code to accommodate any grid topology or grid configuration with …

A Computational Aerodynamic Design Optimization Method Using Sensitivity Analysis, Mohamed E. Eleshaky

Mechanical & Aerospace Engineering Theses & Dissertations

A new and efficient procedure for aerodynamic shape optimization is presented. The salient lineaments of this procedure are: (1) using a discrete sensitivity analysis approach to determine analytically the aerodynamic sensitivity coefficients; (2) obtaining the flowfield solution either by a computational fluid dynamics (CFD) analysis or, alternatively, by a flowfield extrapolation method which is based on a truncated Taylor's series; (3) defining the aerodynamic shape in such a way that it is not restricted to any class of surfaces and the optimizer automatically shapes the aerodynamic configuration to any arbitrary geometry; and (4) requiring no expertise other than that needed …

Shape Sensitivity Analysis And Optimization Of Skeletal Structures And Geometrically Nonlinear Solids, Ching-Hung Chuang

Mechanical & Aerospace Engineering Theses & Dissertations

Formulations and computational schemes for shape design sensitivity analysis and optimization have been developed for both skeletal structures and geometrically nonlinear elastic solids. The continuum approach, which is based on the weak variational form of the governing differential equation and the concept of the material derivative, plays a central role in such a development.

In the first part of this work, the eigenvalue and eigenvector sensitivity equations for skeletal structures are derived with respect to configuration variables of joint and support locations. This derivation is done by the domain method as well as the boundary method. The discrete approach for …

Control Of Low-Speed Turbulent Separated Flow Over A Backward-Facing Ramp, John C. Lin

Mechanical & Aerospace Engineering Theses & Dissertations

The relative performance and flow phenomena associated with several devices for controlling turbulent separated flow were investigated at low speeds. Relative performance of the devices was examined for flow over a curved, backward-facing ramp in a wind tunnel, and the flow phenomena were examined in a water tunnel using dye-flow visualization. Surface static pressure measurements and oil-flow visualization results from the wind tunnel tests indicated that transverse grooves, longitudinal grooves, submerged vortex generators, vortex generator jets (VGJ’s), Viets’ fluidic flappers, elongated arches +a + a (positive angle of attack), and large-eddy breakup devices (LEBU’s) +a + a placed near the …

Unsteady Euler And Navier-Stokes Computations Around Oscillating Delta Wing Including Dynamics, Ahmed Abd-El-Bar Ahmed Salman

Mechanical & Aerospace Engineering Theses & Dissertations

Unsteady flows around rigid or flexible delta wings with and without oscillating leading-edge flaps are considered. These unsteady flow problems are categorized under two classes of problems. In the first class, the wing motion is prescribed a priori and in the second class, the wing motion is obtained as a part of the solution. The formulation of the first class includes either the unsteady Euler or unsteady Navier-Stokes equations for the fluid dynamics and the unsteady linearized Navier-displacement (ND) equations for the grid deformation.

The problem of unsteady transonic flow past a bicircular-arc airfoil undergoing prescribed thickening-thinning oscillation is studied …

Explicit Multistage Schemes For The Solution Of The Three-Dimensional Compressible Euler/Navier-Stokes Equations, Alaa Ahmed Elmiligui

Mechanical & Aerospace Engineering Theses & Dissertations

The objective of this study was to develop a high-resolution-explicit-multi-block numerical algorithm, suitable for efficient computation of the three-dimensional, time-dependent Euler and Navier-Stokes equations. The resulting algorithm has employed a finite volume approach, using MUSCL-type differencing to obtain state variables at cell interface. Variable interpolations were written in the $\kappa$-scheme formulation. Inviscid fluxes were calculated via Roe's flux-difference splitting, and van Leer's flux-vector splitting techniques, which are considered state of the art. The viscous terms were discretized using a second-order, central-difference operator.

Two classes of explicit time integration has been investigated for solving the compressible inviscid/viscous flow problems--two-stage predictor-corrector schemes, …

Finite Element Frequency Domain Solution Of Nonlinear Panel Flutter With Temperature Effects And Fatigue Life Analysis, David Yongxiang Xue

Mechanical & Aerospace Engineering Theses & Dissertations

A frequency domain solution method for nonlinear panel flutter with thermal effects using a consistent finite element formulation has been developed. The von Karman nonlinear strain-displacement relation is used to account for large deflections, the quasi-steady first-order piston theory is employed for aerodynamic loading and the quasi-steady thermal stress theory is applied for the thermal stresses with a given change of the temperature distribution, ΔΤ (x, y, z). The equation of motion under a combined thermal-aerodynamic loading can be mathematically separated into two equations and then solved in sequence: (1) thermal-aerodynamic postbuckling and (2) limit-cycle oscillation. The Newton-Raphson iteration technique …

Integrated System Identification And Adaptive State Estimation For Control Of Flexible Space Structures, Chung-Wen Chen

Mechanical & Aerospace Engineering Theses & Dissertations

Accurate state information is crucial for control of flexible space structures in which the state feedback strategy is used. The performance of a state estimator relies on accurate knowledge about both the system and its disturbances, which are represented by system model and noise covariances respectively. For flexible space structures, due to their great flexibility, obtaining good models from ground testing is not possible. In addition, the characteristics of the systems in operation may vary due to temperature gradient, reorientation, and deterioration of material, etc. Moreover, the disturbances during operation are usually not known. Therefore, adaptive methods for system identification …

Large-Amplitude Finite Element Flutter Analysis Of Composite Panels In Hypersonic Flow, Carl E. Gray Jr.

Mechanical & Aerospace Engineering Theses & Dissertations

A finite-element approach is presented for determining the nonlinear flutter characteristics of two-dimensional isotropic and three-dimensional composite laminated thin panels using the third-order-piston, transverse loading, aerodynamic theory. The unsteady, hypersonic, aerodynamic theory and the von Karman large deflection plate theory are used to formulate the aeroelastic problem. Nonlinear flutter analyses are performed to assess the influence of the higher-order aerodynamic theory on the structure's limit-cycle amplitude and the dynamic pressure of the flow velocity. A solution procedure is presented to solve the nonlinear panel flutter and large-amplitude free vibration finite element equations. This procedure is a linearized updated mode with …

Two-Dimensional Heat Loss From A Building Slab Including Convective Effects In Saturated Soil, William W. Rust Iii

Mechanical & Aerospace Engineering Theses & Dissertations

The heat loss from a building slab was investigated. The continuity equation, Darcy's Law and the energy equation were formulated to include the temperature dependence of viscosity and density of water. The governing equations and appropriate boundary conditions were transformed into dimensionless variables. A finite difference numerical scheme was constructed based on the Gauss-Seidel method by lines and solved iteratively in alternating directions. A correlation between the geometrical characteristics of the domain, the convective surface heating parameters, and the total nondimensional slab heat loss in two dimensions was discovered. Furthermore, the correlation was extended to three-dimensional slabs and produced good …

Kinematic Synthesis Of Deployable-Foldable Truss Structures Using Graph Theory, Dirk B. Warnaar

Mechanical & Aerospace Engineering Theses & Dissertations

A graph theoretic approach is applied to the conceptual design of deployable truss structures. The characteristics that relate to the inter-connectivity of the elements of a deployable truss structure can be captured in a schematic representation, called a graph. A procedure is presented that enables the exhaustive generation of these graphs for structures of any given number of nodes and links and which are foldable onto a plane or onto a line.

A special type of truss structures, called truss modules, is presented. Graphs of this class of structures form a subset of the graphs of truss structures. Two procedures …

Prediction And Control Of Asymmetric Vortical Flows Around Slender Bodies Using Navier-Stokes Equations, Tin-Chee Wong

Mechanical & Aerospace Engineering Theses & Dissertations

Steady and unsteady vortex-dominated flows around slender bodies at high angles of attack are solved using the unsteady, compressible Navier-Stokes equations. An implicit upwind, finite-volume scheme is used for the numerical computations.

For supersonic flows past pointed bodies, the locally-conical flow assumption has been used. Asymmetric flows past five-degree semiapex cones using the thin-layer Navier-Stokes equations at different angles of attack, freestream Mach numbers, Reynolds numbers, grid fineness, computational domain size, sources of disturbances and cross-section shapes have been studied. The onset of flow asymmetry occurs when the relative incidence of pointed forebodies exceeds certain critical values. At these critical …

Numerical Modeling Of Flame Lift-Off Phenomenon And Calculation Of Thermal Loads On A Methane Fuel Injector With Complex Geometry, Taj O. Mohieldin

Mechanical & Aerospace Engineering Theses & Dissertations

A numerical study has been conducted to analyze a fuel injector with three in-line cylinder geometry that has been adopted as a model for investigating the combustion phenomenon in the 8-Foot High Temperature Tunnel (HTT) combustor at the NASA Langley Research Center. The primary objective here is to analyze the flame lift-off phenomenon in the three cylinder fuel injector geometry in two-dimensions. The fluid mechanics model used in the analysis includes time-averaged Navier-Stokes equations that are employed in conjunction with a two-equation k-$\epsilon$ model for predicting the effects of turbulence. Calculations were performed with three chemistry models, namely fast chemistry, …

Integrated Control Of Thermally Distorted Large Space Antennas, Robert H. Tolson

Mechanical & Aerospace Engineering Theses & Dissertations

Studies on controlling the thermal distortion of large space antennae have generally investigated a single orbital position and have optimized actuator locations based on minimizing the RMS surface deviation from the original parabolic shape. One study showed the benefits of directly using far zone electric field characteristics as the performance measure; but, this approach resulted in a nonlinear programming problem. The objective of the current study is to develop an approach to designing a control system that (1) recognizes the time dependence of the distortion and (2) controls variables that are directly related to far field performance in a quadratic …

Navier-Stokes Simulations Of Flows About Complex Configurations Using Domain Decomposition Techniques, Kamran Fouladi-Semnani

Mechanical & Aerospace Engineering Theses & Dissertations

An algorithm is developed to obtain numerical simulations of flows about complex configurations composed of multiple and nonsimilar components with arbitrary geometries. The algorithm uses a hybridization of the domain decomposition techniques for grid generation and to reduce the computer memory requirement. Three dimensional, Reynolds-averaged, unsteady, compressible, and complete Navier-Stokes equations are solved on each of the subdomains by a fully-vectorized, finite-volume, upwind-biased, approximately-factored, and multigrid method. The effect of Reynolds stresses is incorporated through an algebraic turbulence model with several modifications for interference flows. The present algorithm combines the advantages of an efficient, geometrically conservative, minimally and automatically dissipative …

An Experimental Study Of An Axisymmetric Turbulent Boundary Layer Disturbed By A Periodic Freestream, Chithrabhanu Koodalattupuram

Mechanical & Aerospace Engineering Theses & Dissertations

Behavior of an axisymmetric equilibrium turbulent boundary layer disturbed by a propeller wake in the freestream was investigated experimentally. Tests were conducted in a low speed wind tunnel and measurements of turbulence quantities were made using an X wire probe and a constant temperature anemometer. The boundary layer flow on a cylindrical body was characterized by measuring its gross parameters and comparing them with classical values. Propeller speed was measured using an electronic circuit whose output signal was also used to trigger hot wire probe. Gross boundary later characteristics of the disturbed boundary layer did not deviate appreciably from the …

Control And Mechanism Interaction For Ground-Based Testing Of Space Structures, Li-Farn Yang

Mechanical & Aerospace Engineering Theses & Dissertations

In the ground-based validation testing, the adverse effect of terrestrial conditions such as a gravitational force interferes with the dynamic behavior of space structures. A suspension system is developed to assess the structural characteristics in a simulated zero-gravity environment. Using a mechanisms approach, the synthesis of a noncircular disk with a torsional spring at its rotational axis is designed to counteract the gravitational force of test structures during the testing.

The multibody dynamics of a flexible steel beam carried on a rigid trolley has been investigated. The system is constructed in such a way that the rapid and large-angle slewing …

Numerical Study Of Three-Dimensional Separation And Flow Control At A Wing/Body Junction, Balakrishnan Lakshmanan

Mechanical & Aerospace Engineering Theses & Dissertations

The problem of three-dimensional separation and flow control at a wing/body junction has been investigated numerically using a three-dimensional Navier-Stokes code. The numerical code employs an algebraic grid generation technique for generating the grid for unmodified junction and an elliptic grid generation technique for filleted fin junction. The results for laminar flow past a blunt fin/flat plate junction demonstrate that after grid refinement, the computations agree with experiment and reveal a strong dependency of the number of vortices at the junction on Mach number and Reynolds number. The numerical results for pressure distribution, particle paths and limiting streamlines for turbulent …

Low Speed Flowfield Characterization By Infrared Measurements Of Surface Temperatures, Ehud Gartenberg

Mechanical & Aerospace Engineering Theses & Dissertations

An experimental program was aimed at identifying areas in low speed aerodynamic research where infrared imaging systems can make significant contributions. Implementing a new technique, a long electrically heated wire was placed across a laminar jet. By measuring the temperature distribution along the wire with the IR imaging camera, the flow behavior was identified. Furthermore, using Nusselt number correlations, the velocity distribution could be deduced. The same approach was used to survey wakes behind cylinders in a wind-tunnel. This method is suited to investigate flows with position dependent velocities, e.g., boundary layers, confined flows, jets, wakes and shear layers. It …

An Adaptive Remeshing Finite Element Method For High-Speed Compressible Flows Using Quadrilateral And Triangular Elements, Gururaja R. Vemaganti

Mechanical & Aerospace Engineering Theses & Dissertations

In this study a new adaptive remeshing method for high speed compressible flow analysis is presented. The method uses quadrilateral elements where possible, and triangles are introduced as needed. The primary goal of this study is to develop a remeshing method which uses both the concepts of unstructured and structured meshes for the finite element analysis to predict accurate aerodynamic heating in problems related to high speed viscous flows. The remeshing method uses a solution based on an old mesh to create a new mesh based on an advancing front technique. In the present implementation, a structured mesh of quadrilaterals …

Effects Of Nose Bluntness And Shock-Shock Interactions On Blunt Bodies In Viscous Hypersonic Flows, Dal J. Singh

Mechanical & Aerospace Engineering Theses & Dissertations

A numerical study has been conducted to investigate the effects of blunt leading edges on the viscous flow field around a hypersonic vehicle such as the proposed National Aero-Space Plane. Attention is focused on two specific regions of the flow field. Analysis of these flow regions is required to accurately predict the overall flow field as well as to get necessary information on localized zones of high pressure and intense heating.

The forebody is modeled by slender cones and ogives with spherically blunted nose. A combination of Navier-Stokes and parabolized Navier-Stokes equations is used to compute the flow field. The …

Numerical Simulation Of Turbulent Flows Past Three-Dimensional Cavities, Shivakumar Srinivasan

Mechanical & Aerospace Engineering Theses & Dissertations

Computations have been performed to simulate turbulent supersonic, transonic, and subsonic flows past three-dimensional deep, transitional, and shallow cavities. Simulation of these self sustained oscillatory flows has been generated through time accurate solutions of Reynolds averaged full Navier-Stokes equations using the explicit MacCormack scheme. The Reynolds stresses have been modeled, using the Baldwin-Lomax algebraic turbulence model with certain modifications. The computational results include instantaneous and time averaged flow properties everywhere in the computational domain. Time series analyses have been performed for the instantaneous pressure values on the cavity floor. Comparison with experimental data is made in terms of the mean …

Viscous Shock Layer Analysis Of Hypersonic Flows Over Long Slender Vehicles, Kam-Pui Lee

Mechanical & Aerospace Engineering Theses & Dissertations

A method for solving the viscous shock-layer equations for hypersonic flows over long slender bodies is presented. The governing equations are solved by employing a spatial-marching implicit finite-difference technique. The two first-order equations, continuity and normal momentum, are solved simultaneously as a coupled set. This method yields a simple and computationally efficient technique.

Flows past hyperboloids and sphere cones with body half angles of five to 35 degrees are considered. The flow conditions included are from high Reynolds numbers at low altitudes to low Reynolds numbers at high altitudes. Detailed comparisons have been made with other predictions and experimental data …

Application Of Advanced Grid Generation Techniques For Flow Field Computations About Complex Configurations, Monchai Kathong

Mechanical & Aerospace Engineering Theses & Dissertations

In the computation of flow fields about complex configurations, it is very difficult to construct a boundary-fitted coordinate system. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the "multiple grids" or "zonal grids" approach, and its applications are investigated in this study. The method is a conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-stage Runge-Kutta time integration. The code is vectorized and programmed …

A Finite Element Formulation For The Large Deflection Random Response Of Thermally Buckled Structures, James Eugene Locke

Mechanical & Aerospace Engineering Theses & Dissertations

The effects of temperature and acoustic loading are included in a theoretical finite element large deflection formulation for thin, isotropic plate and beam type structures. Thermal loads are applied as steady-state temperature distributions, and acoustic loads are taken to be stationary and Gaussian with zero mean and uniform magnitude and phase over the surface of the structure. Material properties are considered to be independent of temperature. Also, inplane and rotary inertia terms are assumed to be neglegible, and all inplane edge conditions are taken to be immovable. For the random vibration analysis, cross correlation terms are included.

The nature of …

Integral Equation Solution Of The Full Potential Equation For Three-Dimensional, Steady, Transonic Wing Flows, Li-Chuan Chu

Mechanical & Aerospace Engineering Theses & Dissertations

An integral equation method for solving the full potential equation has been developed for three dimensional transonic vortex-wing flows. This method is capable of capturing shocks using the Murman-Cole type of finite difference scheme and is capable of predicting accurate and force-free wake shape as well.

Reading the full potential equation as Poisson's equation, the solution for the velocity field has been expressed in terms of an integral equation using Green's theorem. The solution consists of a surface integral of vorticity distribution on the wing and its free-vortex sheets and a volume integral of source distribution within a computational region …

Investigation Of Supersonic Chemically Reacting And Radiating Channel Flow, Mortaza Mani

Mechanical & Aerospace Engineering Theses & Dissertations

The two-dimensional time dependent Navier-Stokes equations are used to investigate supersonic flows undergoing finite rate chemical reaction and radiation interaction for a hydrogen-air system. The explicit multi-stage finite volume technique of Jameson is used to advance the governing equations in time until convergence is achieved. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The multi-dimensional radiative transfer equations for a nongray model are provided for a general configuration, and then reduced for a planer geometry. Both pseudo-gray and nongray models are used to represent the absorption-emission characteristics of the …

A Study Of Flow Separation In Transonic Flow Using Inviscid And Viscous Cfd Schemes, James Andrew Rhodes

Mechanical & Aerospace Engineering Theses & Dissertations

A comparison of flow separation in transonic flows is made using various computational schemes which solve the Euler and the Navier-Stokes equations of fluid mechanics. The flows examined are computed using several simple two-dimensional configurations including a backward facing step and a bump in a channel. Comparison of the results obtained using shock fitting and flux vector splitting methods are presented and the results obtained using the Euler codes are compared to results on the same configurations using a code which solves the Navier-Stokes equations.

Topology And Grid Adaption For High-Speed Flow Computations, Jamshid Samareh Abolhassani

Mechanical & Aerospace Engineering Theses & Dissertations

This study investigates the effects of grid topology and grid adaption on numerical solutions of the Navier-Stokes equations. In the first part of this study, a general procedure is presented for computation of high-speed flow over complex three-dimensional configurations. This includes the grid generation and solution algorithm for Navier-Stokes equations in a general three-dimensional curvilinear coordinate system. The flow field is simulated on the surface of a Butler wing in a uniform stream. Results are presented for Mach number 3.5 and a Reynolds number of 2,000,000. The O-type and H-type grids have been used for this study, and the results …

A Numerical Study Of Three-Dimensional Vortex Breakdown, Robert Edward Spall

Mechanical & Aerospace Engineering Theses & Dissertations

A numerical simulation of vortex breakdown using the time-dependent Navier-Stokes equations was performed. Unlike previous studies, the numerical algorithm, formulated in terms of the velocity and vorticity, was not restricted to axisymmetric flows. Prototype vortices were parameterized in terms of their Reynolds number and Rossby number. Two cases were studied extensively. In one case, a vortex was imbedded in a uniform free stream. In the other case, the axial velocity of the free stream was decelerated in the streamwise direction. For both cases breakdown occurred when the Rossby number was below a critical value. Vortex lines, particle traces and velocity …