Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational fluid dynamics simulation”

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.

At least 523 records · Page 29

Development of an Efficient CFD Model for Nuclear Thermal Thrust Chamber Assembly Design

The objective of this effort is to develop an efficient and accurate computational methodology to predict both detailed thermo-fluid environments and global characteristics of the internal ballistics for a hypothetical solid-core nuclear thermal thrust chamber assembly (NTTCA). Several numerical and multi-physics thermo-fluid models, such as real fluid, chemically reacting, turbulence, conjugate heat transfer, porosity, and power generation, were incorporated into an unstructured-grid, pressure-based computational fluid dynamics solver as the underlying computational methodology. The numerical simulations of detailed thermo-fluid environment of a single flow element provide a mechanism to estimate the thermal stress and possible occurrence of the mid-section corrosion of the solid core. In addition, the numerical results of the detailed simulation were employed to fine tune the porosity model mimic the pressure drop and thermal load of the coolant flow through a single flow element. The use of the tuned porosity model enables an efficient simulation of the entire NTTCA system, and evaluating its performance during the design cycle.

Cheng, Gary↗

Experimental, Theoretical, and Computational Investigation of Separated Nozzle Flows

A detailed experimental, theoretical, and computational study of separated nozzle flows has been conducted. Experimental testing was performed at the NASA Langley 16-Foot Transonic Tunnel Complex. As part of a comprehensive static performance investigation, force, moment, and pressure measurements were made and schlieren flow visualization was obtained for a sub-scale, non-axisymmetric, two-dimensional, convergent- divergent nozzle. In addition, two-dimensional numerical simulations were run using the computational fluid dynamics code PAB3D with two-equation turbulence closure and algebraic Reynolds stress modeling. For reference, experimental and computational results were compared with theoretical predictions based on one-dimensional gas dynamics and an approximate integral momentum boundary layer method. Experimental results from this study indicate that off-design overexpanded nozzle flow was dominated by shock induced boundary layer separation, which was divided into two distinct flow regimes; three- dimensional separation with partial reattachment, and fully detached two-dimensional separation. The test nozzle was observed to go through a marked transition in passing from one regime to the other. In all cases, separation provided a significant increase in static thrust efficiency compared to the ideal prediction. Results indicate that with controlled separation, the entire overexpanded range of nozzle performance would be within 10% of the peak thrust efficiency. By offering savings in weight and complexity over a conventional mechanical exhaust system, this may allow a fixed geometry nozzle to cover an entire flight envelope. The computational simulation was in excellent agreement with experimental data over most of the test range, and did a good job of modeling internal flow and thrust performance. An exception occurred at low nozzle pressure ratios, where the two-dimensional computational model was inconsistent with the three-dimensional separation observed in the experiment. In general, the computation captured the physics of the shock boundary layer interaction and shock induced boundary layer separation in the nozzle, though there were some differences in shock structure compared to experiment. Though minor, these differences could be important for studies involving flow control or thrust vectoring of separated nozzles. Combined with other observations, this indicates that more detailed, three-dimensional computational modeling needs to be conducted to more realistically simulate shock-separated nozzle flows.

Hunter, Craig A.↗

Development of a Computational Chemical Vapor Deposition Model: Applications to Indium Nitride and Dicyanovinylaniline

A computational chemical vapor deposition (CVD) model is presented, that couples chemical reaction mechanisms with fluid dynamic simulations for vapor deposition experiments. The chemical properties of the systems under investigation are evaluated using quantum, molecular and statistical mechanics models. The fluid dynamic computations are performed using the CFD-ACE program, which can simulate multispecies transport, heat and mass transfer, gas phase chemistry, chemistry of adsorbed species, pulsed reactant flow and variable gravity conditions. Two experimental setups are being studied, in order to fabricate films of: (a) indium nitride (InN) from the gas or surface phase reaction of trimethylindium and ammonia; and (b) 4-(1,1)dicyanovinyl-dimethylaminoaniline (DCVA) by vapor deposition. Modeling of these setups requires knowledge of three groups of properties: thermodynamic properties (heat capacity), transport properties (diffusion, viscosity, and thermal conductivity), and kinetic properties (rate constants for all possible elementary chemical reactions). These properties are evaluated using computational methods whenever experimental data is not available for the species or for the elementary reactions. The chemical vapor deposition model is applied to InN and DCVA. Several possible InN mechanisms are proposed and analyzed. The CVD model simulations of InN show that the deposition rate of InN is more efficient when pulsing chemistry is used under conditions of high pressure and microgravity. An analysis of the chemical properties of DCVA show that DCVA dimers may form under certain conditions of physical vapor transport. CVD simulations of the DCVA system suggest that deposition of the DCVA dimer may play a small role in the film and crystal growth processes.

Cardelino, Carlos↗

Design of a variable area diffuser for a 15-inch Mach 6 open-jet tunnel

The Langley 15-inch Mach 6 High Temperature Tunnel was recently converted from a Mach 10 Hypersonic Flow Apparatus. This conversion was effected to improve the capability of testing in Mach 6 air at relatively high reservoir temperatures not previously possible at Langley. Elevated temperatures allow the matching of the Mach numbers, Reynolds numbers, and ratio of wall-to-adiabatic-wall temperatures (TW/Taw) between this and the Langley 20-inch Mach 6 CF4 Tunnel. This ratio is also matched for Langley's 31-inch Mach 10 Tunnel and is an important parameter useful in the simulation of slender bodies such as National Aerospace Plane (NASP) configurations currently being studied. Having established the nozzle's operating characteristics, the decision was made to install another test section to provide model injection capability. This test section is an open-jet type, with an injection system capable of injecting a model from retracted position to nozzle centerline between 0.5 and 2 seconds. Preliminary calibrations with the new test section resulted in Tunnel blockage. This blockage phenomenon was eliminated when the conical center body in the diffuser was replaced. The issue then, is to provide a new and more efficient variable area diffuser configuration with the capability to withstand testing of larger models without sending the Tunnel into an unstart condition. Use of the 1-dimensional steady flow equation with due regard to friction and heat transfer was employed to estimate the required area ratios (exit area / throat area) in a variable area diffuser. Correlations between diffuser exit Mach number and area ratios, relative to the stagnation pressure ratios and diffuser inlet Mach number were derived. From these correlations, one can set upper and lower operating pressures and temperatures for a given diffuser throat area. In addition, they will provide appropriate input conditions for the full 3-dimensional computational fluid dynamics (CFD) code for further simulation studies.

Loney, Norman W.↗

Computational fluid dynamics - Current capabilities and directions for the future

Computational fluid dynamics (CFD) has made great strides in the detailed simulation of complex fluid flows, including some of those not before understood. It is now being routinely applied to some rather complicated problems and starting to affect the design cycle of aerospace flight vehicles and their components. It is being used to complement, and is being complemented by, experimental studies. Several examples are presented in the paper to illustrate the current state of the art. Included is a discussion of the barriers to accomplishing the basic objective of numerical simulation. In addition, the directions for the future in the discipline of computational fluid dynamics are addressed.

Kutler, Paul↗

A Comparison of Rotor Disk Modeling and Blade-Resolved CFD Simulations for NASA's Tiltwing Air Taxi

A multi-fidelity computational fluid dynamics analysis is carried out for NASA’s tiltwing air taxi concept operating in airplane and helicopter mode. High-fidelity simulations are computationally expensive due to individual rotor blade modeling in a time-dependent computational domain with rotating grids. The mid-fidelity rotor disk option, in its source term implementation, is explored as a more affordable alternative. Computations are performed with NASA’s OVERFLOW flow solver loosely-coupled with the comprehensive code CAMRAD II for appropriate rotor trim. Detailed comparisons are shown for the trim solution, airloads, wake geometry, and rotor performance. While the rotor disk model is able to capture the flow field with satisfactory agreement in airplane mode, it faces difficulties in helicopter mode due to the three-dimensional effects of the wake. Although this study is limited to a specific vehicle geometry, it is expected that the results are somewhat generalizable to the analysis of multi-rotor configurations.

ARMD↗

LeRC reduced gravity fluid management technology program

A survey of the reduced gravity fluid management technology program is presented. Information on reduced gravity fluid behavior, techniques for thermal control of cryogenic tankage, and design for fluid management systems are discussed. The development of Spacelab experiments, propellant management systems for orbit transfer vehicles, and computer techniques for simulating reduced gravity fluid dynamic processes is reported.

Aydelott, J. C.↗

Geometric methods in computational fluid dynamics

General methods for the construction of geometric computational fluid dynamic algorithms are presented which simulate a variety of flow fields in various nontrivial regions. Included are: basic developments with tensors; various forms for the equations of motion; generalized numerical methods and boundary conditions; and methods for mesh generation to meet the strong geometric constraints of turbomachines. Coordinate generation is shown generally to yield mesh descriptions from one or more transformations that are smoothly joined together to form a composite mesh.

Eiseman, P. R.↗

Transonic aeroelastic analysis of launch vehicle configurations

A numerical study of the aeroelastic stability of typical launch vehicle configurations in transonic flight is performed. Recent computational fluid dynamics techniques are used to simulate the transonic aerodynamic flow fields, as opposed to relying on experimental data for the unsteady aerodynamic pressures. The flow solver is coupled to an appropriate structural representation of the vehicle. The aerodynamic formulation is based on the thin layer approximation to the Reynolds-Averaged Navier-Stokes equations, where the account for turbulent mixing is done by the two-layer Baldwin and Lomax algebraic eddy viscosity model. The structural-dynamic equations are developed considering free-free flexural vibration of an elongated beam with variable properties and are cast in modal form. Aeroelastic analyses are performed by integrating simultaneously in the two sets of equations. By tracing the growth or decay of a perturbed oscillation, the aeroelastic stability of a given constant configuration can be ascertained. The method is described in detail, and results that indicate its application are presented. Applications include some validation cases for the algorithm developed, as well as the study of configurations known to have presented flutter programs in the past.

Filgueirasdeazevedo, Joao Luiz↗

Implementation of a 3D mixing layer code on parallel computers

This paper summarizes our progress and experience in the development of a Computational-Fluid-Dynamics code on parallel computers to simulate three-dimensional spatially-developing mixing layers. In this initial study, the three-dimensional time-dependent Euler equations are solved using a finite-volume explicit time-marching algorithm. The code was first programmed in Fortran 77 for sequential computers. The code was then converted for use on parallel computers using the conventional message-passing technique, while we have not been able to compile the code with the present version of HPF compilers.

Roe, K.↗

Parallel ALLSPD-3D: Speeding Up Combustor Analysis Via Parallel Processing

The ALLSPD-3D Computational Fluid Dynamics code for reacting flow simulation was run on a set of benchmark test cases to determine its parallel efficiency. These test cases included non-reacting and reacting flow simulations with varying numbers of processors. Also, the tests explored the effects of scaling the simulation with the number of processors in addition to distributing a constant size problem over an increasing number of processors. The test cases were run on a cluster of IBM RS/6000 Model 590 workstations with ethernet and ATM networking plus a shared memory SGI Power Challenge L workstation. The results indicate that the network capabilities significantly influence the parallel efficiency, i.e., a shared memory machine is fastest and ATM networking provides acceptable performance. The limitations of ethernet greatly hamper the rapid calculation of flows using ALLSPD-3D.

Fricker, David M.↗

RSRM and ETM03 Internal Flow Simulations and Comparisons

ETM03 (Engineering Test Motor-03) is an extended length RSRM (Reusable Solid Rocket Motor) designed to increase motor performance and create more severe internal environments compared with the standard four-segment RSRM motor configuration. This is achieved primarily through three unique design features. First is the incorporation of an additional RSRM center segment, second is a slight increase in throat diameter, and third is the use of an Extended Aft Exit Cone (EAEC). As a result of these design features, parameters such as web time, action time, head end pressure, web time average pressure, maximum thrust, mass flow rate, centerline Mach number, pressure and thrust integrals have all increased compared with nominal RSRM values. In some cases these increases are substantial. The primary objective of the ETM03 test program is to provide a platform for RSRM component margin testing. Test results will not only provide direct data concerning component performance under more adverse conditions, but serve as a second design data point for developing, validating and enhancing component analytical modeling techniques. To help component designers assess how the changes in motor environment will affect performance, internal flow simulations for both the nominal RSRM and ETM03 motor designs were completed to obtain comparisons of aero-thermal boundary conditions and system loads. Full geometries for both motors were characterized with two-dimensional axi-symmetric models at burn times of 1, 20, 54, 67 and 80-seconds. A sixth set considered burn times of 110 and 117-seconds for RSRM and ETM03, respectively. The simulations were performed using the computational fluid dynamics (CFD) commercial code FLUENT (trademark). Of particular interest were any differences between the two motor environments that could lead to a significant increase in system loads, or in internal insulation and/or nozzle component charring and erosion in ETM03 compared with RSRM. Based on these comparative analyses conducted in this study, the objective of ETM03 will be achieved by providing a more adverse operating environment for motor components than the nominal RSRM environment. For example: Higher chamber pressure drop in ETM03 than in RSRM; higher centerline Mach numbers approaching the nozzle in ETM03 than in RSRM; higher heat transfer rates for the internal insulation and nozzle components in ETM03 than in RSRM; and higher levels of droplet impingement and slag accumulation in ETM03 than in the RSRM.

Ahmad, R. A.↗

Combustion Model of Supersonic Rocket Exhausts in an Entrained Flow Enclosure

This paper describes the Computation Fluid Dynamics (CFD) model developed to simulate the supersonic rocket exhaust in an entrained flow cylinder. The model can be used to study the plume-induced environment due to static firing test of the Taurus II launch vehicle. The finite rate chemistry is used to model the combustion process involving rocket propellant (RP 1) and liquid oxidizer (LOX). A similar chemical reacting model is also used to simulate the mixing of rocket plume and ambient air. The model provides detailed information on the gas concentration and other flow parameters within the enclosed region thus allowing different operating scenarios to be examined in an efficient manner. It is shown that the real gas influence is significant and yields better agreement with the theory.

Vu, Bruce↗

Combustion Model of Supersonic Rocket Exhausts in an Entrained Flow Enclosure

This paper describes the Computational Fluid Dynamics (CFD) model developed to simulate the supersonic rocket exhaust in an entrained flow cylinder. The model can be used to study the plume-induced environment due to static firing tests of the Taurus-II launch vehicle. The finite-rate chemistry is used to model the combustion process involving rocket propellant (RP-1) and liquid oxidizer (LOX). A similar chemical reacting model is also used to simulate the mixing of rocket plume and ambient air. The model provides detailed information on the gas concentration and other flow parameters within the enclosed region, thus allowing different operating scenarios to be examined in an efficient manner. It is shown that the real gas influence is significant and yields better agreement with the theory.

Vu, Bruce T.↗

Extension of the Viscous Collision Limiting Direct Simulation Monte Carlo Technique to Multiple Species

There are many flows fields that span a wide range of length scales where regions of both rarefied and continuum flow exist and neither direct simulation Monte Carlo (DSMC) nor computational fluid dynamics (CFD) provide the appropriate solution everywhere. Recently, a new viscous collision limited (VCL) DSMC technique was proposed to incorporate effects of physical diffusion into collision limiter calculations to make the low Knudsen number regime normally limited to CFD more tractable for an all-particle technique. This original work had been derived for a single species gas. The current work extends the VCL-DSMC technique to gases with multiple species. Similar derivations were performed to equate numerical and physical transport coefficients. However, a more rigorous treatment of determining the mixture viscosity is applied. In the original work, consideration was given to internal energy non-equilibrium, and this is also extended in the current work to chemical non-equilibrium.

Liechty, Derek S.↗

Predicting Orion Launch Abort Acoustics

Computational fluid dynamics (CFD) scale-resolving simulations were performed to help assess the Orion launch abort system’s vibro-acoustic environment in collaboration with the Orion Loads and Dynamics team at the Johnson Space Center (JSC). The primary objectives are to establish best practices and to validate CFD predictions for surface fluctuating pressure spectra with ground, wind tunnel, and flight test data. Excellent agreement is obtained between CFD predictions by the Cartesian adaptive mesh refinement (AMR) module of the Launch, Ascent, and Vehicle Aerodynamics (LAVA) solver and measurements from a ground test representative of the start of a pad abort scenario, and with measurements from a wind tunnel test case corresponding to a transonic, high total angle of attack ascent abort scenario. A simplified (static and at rest) simulation of a pad abort flight test shows good agreement with the flight data, whereas a more costly simulation where the vehicle is accelerating relative to the ground from rest to Mach 0.4 shows consistent over-prediction of the spectra. Simulation results for a supersonic ascent abort flight test are in reasonable agreement with the measurements, but show larger excursions in regions where non-linear interactions between the plume turbulence, shocks, and turbulent boundary layer are strongest – likely due to CFD modeling assumptions, the effect of which are discussed. Overall, the CFD results demonstrate that scale-resolving simulation methodologies like the one presented can help characterize full-scale launch abort vehicle near-field noise and can complement wind tunnel, ground and flight tests by providing trends and reducing uncertainty for scenarios that cannot or are not planned to be tested.

Orion↗

PlanetarY Telemetric Helicopter for Investigation and Analysis (PYTHIA): A Rotorcraft for Martian Lava Tube Exploration

The PlanetarY Telemetric Helicopter for Investigation and Analysis (PYTHIA) project presents an early feasibility study into future Mars rotorcraft. With the success of Ingenuity and the current development of the Mars Sample Retrieval Helicopters, there is motivation to explore additional vehicle concepts for Mars exploration. This work presents an early conceptual design of a lava tube exploring quadrotor. The nominal mission for the PYTHIA quadrotor includes a two-phase in-depth exploration of one or mutliple lava tubes in the area of interest, Arsia Mons. A sizing analysis was completed using the NASA Design and Analysis of Rotorcraft (NDARC) tool. An eight-bladed quadrotor was selected based on the NDARC sizing sweeps. With preliminary vehicle and blade sizing completed, a flow visualization study was conducted. Two simulations were completed in Rotorcraft Computational Fluid Dynamics (RotCFD): one of the quadrotor alone and another in a lava tube. This paper introduces a baseline reference mission, preliminary vehicle design, and initial sizing analysis of a rotorcraft to explore Martian lava tubes.

PlanetarY↗

Propulsion Simulations with the Unstructured-Grid CFD Tool TetrUSS

A computational investigation has been completed to assess the capability of the NASA Tetrahedral Unstructured Software System (TetrUSS) for simulation of exhaust nozzle flows. Three configurations were chosen for this study: (1) a fluidic jet effects model, (2) an isolated nacelle with a supersonic cruise nozzle, and (3) a fluidic pitchthrust- vectoring nozzle. These configurations were chosen because existing data provided a means for measuring the ability of the TetrUSS flow solver USM3D for simulating complex nozzle flows. Fluidic jet effects model simulations were compared with structured-grid CFD (computational fluid dynamics) data at Mach numbers from 0.3 to 1.2 at nozzle pressure ratios up to 7.2. Simulations of an isolated nacelle with a supersonic cruise nozzle were compared with wind tunnel experimental data and structured-grid CFD data at Mach numbers of 0.9 and 1.2, with a nozzle pressure ratio of 5. Fluidic pitch-thrust-vectoring nozzle simulations were compared with static experimental data and structured-grid CFD data at static freestream conditions and nozzle pressure ratios from 3 to 10. A fluidic injection case was computed with the third configuration at a nozzle pressure ratio of 4.6 and a secondary pressure ratio of 0.7. Results indicate that USM3D with the S-A turbulence model provides accurate exhaust nozzle simulations at on-design conditions, but does not predict internal shock location at overexpanded conditions or pressure recovery along a boattail at transonic conditions.

Deere, Karen A.↗