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.
The impact of computational fluid dynamics on the design of fluid dynamics devices is discussed. An efficient numerical procedure for solving the Navier-Stokes equations in three dimensions is discussed and flow simulations made on an inexpensive desk top work station are presented.
An unsteady, three dimensional Navier-Stokes solution in rotating frame formulation for turbomachinery applications is presented. Casting the governing equations in a rotating frame enabled the freezing of grid motion and resulted in substantial savings in computer time. The turbine blade was computationally simulated and probabilistically evaluated in view of several uncertainties in the aerodynamic, structural, material and thermal variables that govern the turbine blade. The interconnection between the computational fluid dynamics code and finite element structural analysis code was necessary to couple the thermal profiles with the structural design. The stresses and their variations were evaluated at critical points on the Turbine blade. Cumulative distribution functions and sensitivity factors were computed for stress responses due to aerodynamic, geometric, mechanical and thermal random variables.
Under the Propulsion Discipline of NASA s Fundamental Aeronautics Program s Hypersonics Project, a test apparatus, for testing a scramjet isolator model, is being constructed at NASA's Langley Research Center. The test apparatus will incorporate a 1-inch by 2-inch by 15-inch-long scramjet isolator model supplied with 2.1 lbm/sec of unheated dry air through a Mach 2.5 converging-diverging nozzle. The planned research will incorporate progressively more challenging measurement techniques to characterize the flow field within the isolator, concluding with the application of the Laser-Induced Thermal Acoustic (LITA) measurement technique. The primary goal of this research is to use the data acquired to validate Computational Fluid Dynamics (CFD) models employed to characterize the complex flow field of a scramjet isolator. This paper describes the test apparatus being constructed, pre-test CFD simulations, and the LITA measurement technique.
Supersonic combustion performance of a bi-component gaseous hydrocarbon fuel mixture is one of the primary aspects under investigation in the HIFiRE Flight 2 experiment. In-flight instrumentation and post-test analyses will be two key elements used to determine the combustion performance. Pre-flight computational fluid dynamics (CFD) analyses provide valuable information that can be used to optimize the placement of a constrained set of wall pressure instrumentation in the experiment. The simulations also allow pre-flight assessments of performance sensitivities leading to estimates of overall uncertainty in the determination of combustion efficiency. Based on the pre-flight CFD results, 128 wall pressure sensors have been located throughout the isolator/combustor flowpath to minimize the error in determining the wall pressure force at Mach 8 flight conditions. Also, sensitivity analyses show that mass capture and combustor exit stream thrust are the two primary contributors to uncertainty in combustion efficiency.
The desire to reduce or eliminate the operational restrictions of supersonic aircraft over populated areas has led to extensive research at NASA. Restrictions are due to the disturbance of the sonic boom, caused by the coalescence of shock waves formed by the aircraft. A study has been performed focused on reducing the magnitude of the sonic boom N-wave generated by airplane components with a focus on shock waves caused by the exhaust nozzle plume. Testing was completed in the 1-foot by 1-foot supersonic wind tunnel to study the effects of an exhaust nozzle plume and shock wave interaction. The plume and shock interaction study was developed to collect data for computational fluid dynamics (CFD) validation of a nozzle plume passing through the shock generated from the wing or tail of a supersonic vehicle. The wing or tail was simulated with a wedgeshaped shock generator. This test entry was the first of two phases to collect schlieren images and off-body static pressure profiles. Three wedge configurations were tested consisting of strut-mounted wedges of 2.5- degrees and 5-degrees. Three propulsion configurations were tested simulating the propulsion pod and aft deck from a low boom vehicle concept, which also provided a trailing edge shock and plume interaction. Findings include how the interaction of the jet plume caused a thickening of the shock generated by the wedge (or aft deck) and demonstrate how the shock location moved with increasing nozzle pressure ratio.
In this article, we present a reduced order modeling approach suitable for active control of fluid dynamical systems based on proper orthogonal decomposition (POD). The rationale behind the reduced order modeling is that numerical simulation of Navier-Stokes equations is still too costly for the purpose of optimization and control of unsteady flows. We examine the possibility of obtaining reduced order models that reduce computational complexity associated with the Navier-Stokes equations while capturing the essential dynamics by using the POD. The POD allows extraction of certain optimal set of basis functions, perhaps few, from a computational or experimental data-base through an eigenvalue analysis. The solution is then obtained as a linear combination of these optimal set of basis functions by means of Galerkin projection. This makes it attractive for optimal control and estimation of systems governed by partial differential equations. We here use it in active control of fluid flows governed by the Navier-Stokes equations. We show that the resulting reduced order model can be very efficient for the computations of optimization and control problems in unsteady flows. Finally, implementational issues and numerical experiments are presented for simulations and optimal control of fluid flow through channels.
An acoustic prediction capability for supersonic axisymmetric jets was developed on the basis of OVERFLOW Navier-Stokes CFD (Computational Fluid Dynamics) code of NASA Langley Research Center. Reynolds-averaged turbulent stresses in the flow field are modeled with the aid of Spalart-Allmaras one-equation turbulence model. Appropriate acoustic and outflow boundary conditions were implemented to compute time-dependent acoustic pressure in the nonlinear source-field. Based on the specification of acoustic pressure, its temporal and normal derivatives on the Kirchhoff surface, the near-field and the far-field sound pressure levels are computed via Kirchhoff surface integral, with the Kirchhoff surface chosen to enclose the nonlinear sound source region described by the CFD code. The methods are validated by a comparison of the predictions of sound pressure levels with the available data for an axisymmetric turbulent supersonic (Mach 2) perfectly expanded jet.
A major component of any liquid propellant rocket is the propellant injection system. Issues of interest include the degree of liquid vaporization and its impact on the combustion process, the pressure and temperature fields in the combustion chamber, and the cooling of the injector face and chamber walls. The Finite Difference Navier-Stokes (FDNS) code is a primary computational tool used in the MSFC Computational Fluid Dynamics Branch. The branch has dedicated a significant amount of resources to development of this code for prediction of both liquid and solid fuel rocket performance. The FDNS code is currently being upgraded to include the capability to model liquid/gas multi-phase flows for fuel injection simulation. An important aspect of this effort is benchmarking the code capabilities to predict existing experimental injection data. The objective of this MSFC/ASEE Summer Faculty Fellowship term was to evaluate the capabilities of the modified FDNS code to predict flow fields with liquid injection. Comparisons were made between code predictions and existing experimental data. A significant portion of the effort included a search for appropriate validation data. Also, code simulation deficiencies were identified.
Electron beam freeform fabrication (EBF3) is a member of an emerging class of direct manufacturing processes known as solid freeform fabrication (SFF); another member of the class is the laser deposition process. Successful application of the EBF3 process requires precise control of a number of process parameters such as the EB power, speed, and metal feed rate in order to ensure thermal management; good fusion between the substrate and the first layer and between successive layers; minimize part distortion and residual stresses; and control the microstructure of the finished product. This is the only effort thus far that has addressed computer simulation of the EBF3 process. The models developed in this effort can assist in reducing the number of trials in the laboratory or on the shop floor while making high-quality parts. With some modifications, their use can be further extended to the simulation of laser, TIG (tungsten inert gas), and other deposition processes. A solid mechanics-based finite element code, ABAQUS, was chosen as the primary engine in developing these models whereas a computational fluid dynamics (CFD) code, Fluent, was used in a support role. Several innovative concepts were developed, some of which are highlighted below. These concepts were implemented in a number of new computer models either in the form of stand-alone programs or as user subroutines for ABAQUS and Fluent codes. A database of thermo-physical, mechanical, fluid, and metallurgical properties of stainless steel 304 was developed. Computing models for Gaussian and raster modes of the electron beam heat input were developed. Also, new schemes were devised to account for the heat sink effect during the deposition process. These innovations, and others, lead to improved models for thermal management and prediction of transient/residual stresses and distortions. Two approaches for the prediction of microstructure were pursued. The first was an empirical approach involving the computation of thermal gradient, solidification rate, and velocity (G,R,V) coupled with the use of a solidification map that should be known a priori. The second approach relies completely on computer simulation. For this purpose a criterion for the prediction of morphology was proposed, which was combined with three alternative models for the prediction of microstructure; one based on solidification kinetics, the second on phase diagram, and the third on differential scanning calorimetry data. The last was found to be the simplest and the most versatile; it can be used with multicomponent alloys and rapid solidification without any additional difficulty. For the purpose of (limited) experimental validation, finite element models developed in this effort were applied to three different shapes made of stainless steel 304 material, designed expressly for this effort with an increasing level of complexity. These finite element models require large computation time, especially when applied to deposits with multiple adjacent beads and layers. This problem can be overcome, to some extent, by the use of fast, multi-core computers. Also, due to their numerical nature coupled with the fact that solid mechanics- based models are being used to represent the material behavior in liquid and vapor phases as well, the models have some inherent approximations that become more pronounced when dealing with multi-bead and multi-layer deposits.
The objective of this research is to simulate steady and unsteady viscous flows, including rotor/stator interaction and tip clearance effects in turbomachinery. The numerical formulation for steady flow developed here includes an efficient grid generation scheme, particularly suited to computational grids for the analysis of turbulent turbomachinery flows and tip clearance flows, and a semi-implicit, pressure-based computational fluid dynamics scheme that directly includes artificial dissipation, and is applicable to both viscous and inviscid flows. The values of these artificial dissipation is optimized to achieve accuracy and convergency in the solution. The numerical model is used to investigate the structure of tip clearance flows in a turbine nozzle. The structure of leakage flow is captured accurately, including blade-to-blade variation of all three velocity components, pitch and yaw angles, losses and blade static pressures in the tip clearance region. The simulation also includes evaluation of such quantities of leakage mass flow, vortex strength, losses, dominant leakage flow regions and the spanwise extent affected by the leakage flow. It is demonstrated, through optimization of grid size and artificial dissipation, that the tip clearance flow field can be captured accurately. The above numerical formulation was modified to incorporate time accurate solutions. An inner loop iteration scheme is used at each time step to account for the non-linear effects. The computation of unsteady flow through a flat plate cascade subjected to a transverse gust reveals that the choice of grid spacing and the amount of artificial dissipation is critical for accurate prediction of unsteady phenomena. The rotor-stator interaction problem is simulated by starting the computation upstream of the stator, and the upstream rotor wake is specified from the experimental data. The results show that the stator potential effects have appreciable influence on the upstream rotor wake. The predicted unsteady wake profiles are compared with the available experimental data and the agreement is good. The numerical results are interpreted to draw conclusions on the unsteady wake transport mechanism in the blade passage.
This is a final report as far as our work at University of Minnesota is concerned. The report describes our research progress and accomplishments in development of high performance computing methods and tools for 3D finite element computation of aerodynamic characteristics and fluid-structure interactions (FSI) arising in airdrop systems, namely ram-air parachutes and round parachutes. This class of simulations involves complex geometries, flexible structural components, deforming fluid domains, and unsteady flow patterns. The key components of our simulation toolkit are a stabilized finite element flow solver, a nonlinear structural dynamics solver, an automatic mesh moving scheme, and an interface between the fluid and structural solvers; all of these have been developed within a parallel message-passing paradigm.
This paper discusses incompressible Navier-Stokes solution methods with an emphasis on the pseudocompressibility method. A steady-state flow solver based on the pseudocompressibility approach is then described. This flow solver code has been used to analyze the internal flow in the Space Shuttle main engine hot-gas manifold. Salient features associated with this three-dimensional realistic flow simulation are discussed. Numerical solutions relevant to the current engine analysis and the redesign effort are discussed along with experimental results. This example demonstrates the potential of computational fluid dynamics as a design tool for aerospace applications.
Incompressible Navier-Stokes solution methods are discussed with an emphasis on the pseudocompressibility method. A steady-state flow solver based on the pseudocompressibility approach is then described. This flow-solver code was used to analyze the internal flow in the Space Shuttle main engine hot-gas manifold. Salient features associated with this three-dimensional realistic flow simulation are discussed. Numerical solutions relevant to the current engine analysis and the redesign effort are discussed along with experimental results. This example demonstrates the potential of computational fluid dynamics as a design tool for aerospace applications.
A computational method to simulate unsteady flows involving moving rigid boundaries and interference has been developed. The method is used to solve inviscid equations governing the fluid flow and the dynamic equations governing the motion of rigid bodies. A second-order accurate, upwind-biased, and alterating-direction-implicit method is used to solve the governing equations of the flow. A kinematic domain decomposition (KDD) procedure is extended to treat 3D problems with a high degree of accuracy and generality. The method under consideration is applied to both transonic and supertransonic flows. Both cases involve flow past a cylinder which is forced to pitch sinusoidally near a vertical wall. Benefits of the proposed approach include accurate calculation of the flow around 3D moving multiple bodies with interference; reduction of a numerical error; in particular, the dispersion error which strongly affects wave propagation; and minimization of the phase error which is accumulated according to the time advance procedure.
A pitot-rake survey of the simulated exhaust of a half-span scramjet nozzle model was conducted in the Langley 20-Inch Mach 6 Tunnel to provide an additional data set for computational fluid dynamics (CFD) code comparisons. A wind-tunnel model was tested with a 26-tube pitot rake that could be manually positioned along the mid-semispan plane of the model. The model configuration had an external expansion surface of 20 degrees and an internal cowl expansion of 12 degrees; tests were also performed with a flow fence. Tests were conducted at a free-stream Reynolds number of approximately 6.5 x 10(exp 6) per foot and a model angle of attack of -0.75 degrees. The two exhaust gas mediums that were tested were air and a Freon 12-argon mixture. Each medium was tested at two jet total pressures at approximately 28 and 14 psia. This document presents the flow-field survey results in graphical as well as tabular form, and several observations concerning the results are discussed. The surveys reveal the major expected flow-field characteristics for each test configuration. For a 50-percent freon 12 and 50-percent argon mixture by volume (Fr-Ar), the exhaust jet pressures were slightly higher than those for air. The addition of a flow fence slightly raised the pitot pressure for the Fr-Ar mixture, but it produced little change for air. For the Fr-Ar exhaust, the plume was larger and the region between the shock wave and plume was smaller.
The Aeroassist Flight Experiment (AFE), scheduled to be performed in 1994, will serve as a precursor for aeroassisted space transfer vehicles (ASTV's) and is representative of entry concepts being considered for missions to Mars. Rationale for the AFE is reviewed briefly as are the various experiments carried aboard the vehicle. The approach used to determine hypersonic aerodynamic and aerothermodynamic characteristics over a wide range of simulation parameters in ground-based facilities is presented. Facilities, instrumentation and test procedures employed in the establishment of the data base are discussed. Measurements illustrating the effects of hypersonic simulation parameters, particularly normal-shock density ratio (an important parameter for hypersonic blunt bodies), and attitude on aerodynamic and aerothermodynamic characteristics are presented, and predictions from computational fluid dynamic (CFD) computer codes are compared with measurement.
COLD-SAT Dynamic Model (CSDM) computer code implements six-degree-of-freedom, rigid-body mathematical model for simulation of spacecraft in orbit around Earth. Investigates flow dynamics and thermodynamics of subcritical cryogenic fluids in microgravity. Consists of three parts: translation model, rotation model, and slosh model. Written in FORTRAN 77.
A combustor liner was computationally simulated and probabilistically evaluated in view of the several uncertainties in the aerodynamic, structural, material and thermal variables that govern the combustor liner. The interconnection between the computational fluid dynamics code and the finite element structural analysis codes was necessary to couple the thermal profiles with structural design. The stresses and their variations were evaluated at critical points on the liner. Cumulative distribution functions and sensitivity factors were computed for stress responses due to the aerodynamic, mechanical and thermal random variables. It was observed that the inlet and exit temperatures have a lot of influence on the hoop stress. For prescribed values of inlet and exit temperatures, the Reynolds number of the flow, coefficient of thermal expansion, gas emissivity and absorptivity and thermal conductivity of the material have about the same impact on the hoop stress. These results can be used to quickly identify the most critical design variables in order to optimize the design and make it cost effective.