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.
This report documents the research conducted to update computer models for dynamic fluid flow simulation of the E-1 test stand subsystems at te NASA John C. Stennis Space Center.Work also involved significant upgrades to the capabilities of EASY/ROCKETS library through the inclusion of the NIST-12 thermodynamic property database and development of new control system modules.
Three-dimensional numerical investigations on the nozzle start-up side load physics were performed. The objective of this study is to identify the three-dimensional side load physics and to compute the associated aerodynamic side load using an anchored computational methodology. The computational methodology is based on an unstructured-grid, and pressure-based computational fluid dynamics formulation, and a simulated inlet condition based on a system calculation. Finite-rate chemistry was used throughout the study so that combustion effect is always included, and the effect of wall cooling on side load physics is studied. The side load physics captured include the afterburning wave, transition from free- shock to restricted-shock separation, and lip Lambda shock oscillation. With the adiabatic nozzle, free-shock separation reappears after the transition from free-shock separation to restricted-shock separation, and the subsequent flow pattern of the simultaneous free-shock and restricted-shock separations creates a very asymmetric Mach disk flow. With the cooled nozzle, the more symmetric restricted-shock separation persisted throughout the start-up transient after the transition, leading to an overall lower side load than that of the adiabatic nozzle. The tepee structures corresponding to the maximum side load were addressed.
The results of a numerical study with flow and boundary conditions based on an experiment of a continuum sonic jet interacting with rarefied flow about a sharp leading edge flat plate at zero incidence are presented. Comparisons are made between computational fluid dynamics (CFD) and direct simulation Monte Carlo (DSMC) solutions which provide an assessment of applying each technique to the flow conditions of the experiment. An analysis of the CFD results revealed a correlation between the interaction interface of the jet continuum breakdown surface and a non-dimensional parameter derived from jet and free stream flow conditions. Using the breakdown surface from the correlation, the continuum jet was uncoupled from the interaction, thus allowing an uncoupled CFD-DSMC solution to be obtained. Also, a nearest neighbor collision algorithm, similar to the subcell technique, was implemented in the DSMC solution technique. The comparison between CFD and DSMC results shows good qualitative agreement in the interaction region and good quantitative agreement elsewhere.
In support of launch vehicle base heating and pressure prediction efforts using the Loci-CHEM Navier-Stokes computational fluid dynamics solver, 35 numerical simulations of the NASA TND-1093 wind tunnel test have been modeled and analyzed. This test article is composed of four JP-4/LOX 500 lbf rocket motors exhausting into a Mach 2 - 3.5 wind tunnel at various ambient pressure conditions. These water-cooled motors are attached to a base plate of a standard missile forebody. We explore the base heating profiles for fully coupled finite-rate chemistry simulations, one-way coupled RAMP (Reacting And Multiphase Program using Method of Characteristics)-BLIMPJ (Boundary Layer Integral Matrix Program - Jet Version) derived solutions and variable and constant specific heat ratio frozen flow simulations. Variations in turbulence models, temperature boundary conditions and thermodynamic properties of the plume have been investigated at two ambient pressure conditions: 255 lb/sq ft (simulated low altitude) and 35 lb/sq ft (simulated high altitude). It is observed that the convective base heat flux and base temperature are most sensitive to the nozzle inner wall thermal boundary layer profile which is dependent on the wall temperature, boundary layer s specific energy and chemical reactions. Recovery shock dynamics and afterburning significantly influences convective base heating. Turbulence models and external nozzle wall thermal boundary layer profiles show less sensitivity to base heating characteristics. Base heating rates are validated for the highest fidelity solutions which show an agreement within +/-10% with respect to test data.
No abstract available
A research program to provide physical insight into disk-gap-band parachute operation in the supersonic regime on Mars was conducted. The program included supersonic wind tunnel tests, computational fluid dynamics and fluid structure interaction simulations. Specifically, the nature and cause of the "area oscillation" phenomenon were investigated to determine the scale, aerodynamic, and aero-elastic dependence of the supersonic parachute collapse and re-inflation event. A variety of non-intrusive, temporally resolved, and high resolution diagnostic techniques were used to interrogate the flow and generate validation datasets. The results of flow visualization, particle image velocimetry, load measurements, and photogrammetric reconstruction will be presented. Implications to parachute design, use, and verification will also be discussed.
As part of comprehensive efforts to develop physics-based risk assessment techniques for space systems at NASA, coupled computational fluid and rigid body dynamic simulations were carried out to investigate the flow mechanisms that accelerate tank fragments in bursting pressurized vessels. Simulations of several configurations were compared to analyses based on the industry-standard Baker explosion model, and were used to formulate an improved version of the model. The standard model, which neglects an external fluid, was found to agree best with simulation results only in configurations where the internal-to-external pressure ratio is very high and fragment curvature is small. The improved model introduces terms that accommodate an external fluid and better account for variations based on circumferential fragment count. Physics-based analysis was critical in increasing the model's range of applicability. The improved tank burst model can be used to produce more accurate risk assessments of space vehicle failure modes that involve high-speed debris, such as exploding propellant tanks and bursting rocket engines.
The primary purpose of this study is to determine the extent to which the size and shape of the wind tunnel sting affect the accuracy of the base pressure corrections applied to measured axial force. The study also includes an assessment of the overall accuracy of the corrections. To accomplish these goals, Computational Fluid Dynamics is used to simulate a simplified version of the geometry of the Space Launch System Block 1B Cargo configuration, paired with a range of wind tunnel sting sizes, over a variety of ascent flight conditions. The base pressure correction method used in the wind tunnel is emulated on the base pressures from the simulated flows and results are compared to direct integration of the base pressures. Differences in results between the two methods provides an assessment of the accuracy of the base force correction method and how that accuracy is affected by sting size.
The Multirotor Test Bed (MTB) is a new capability for testing a wide array of advanced vertical take-off and landing (VTOL) rotor configurations, with a primary focus on testing in the U.S. Army 7- by 10-ft Wind Tunnel at NASA Ames Research Center. The MTB was designed to allow adjustment of the vertical, lateral, and longitudinal placement of each rotor, as well as allow tilt adjustment of each rotor and pitch adjustment of the whole assembly. Each rotor can tilt forward 90 deg and backwards 5 deg. In addition, the entire MTB can tilt forward 20 deg and backwards 10 deg. This flexibility allows the system to be tested in many different configurations. There is a six-axis load cell under each rotor assembly, to measure both the steady and dynamic loads produced by each rotor. The wind tunnel scales can measure loads on the full assembly. The overall goal of the MTB project is to help gain a better understanding of the performance, control, interactional aerodynamics, and acoustics of multirotor systems. A hybrid CFD tool called RotCFD (Rotorcraft Computational Fluid Dynamics) was used to simulate the MTB in several testing configurations. This paper explains the method of running the RotCFD simulations and explores the results from the simulations. The objective of this paper is to compare the RotCFD simulation results with the MTB wind tunnel test data, seeking to further validate RotCFD for multirotor systems and assess the influence of aerodynamic interactions on individual rotor performance.
Advanced multirotor vertical flight aircraft concepts are emerging faster than rigorous individualized tests can investigate their utility and performance. There are several analysis tools that predict multirotor performance and flow characteristics, but the accuracy of these predictions is still being debated due to lack of experimental data from multirotor tests that are needed to validate the analyses. The objective of this paper is to simulate multirotor configurations using two different mid-fidelity rotorcraft analysis tools, Comprehensive Hierarchical Aeromechanics Rotorcraft Model (CHARM) and Rotorcraft Computational Fluid Dynamics (RotCFD), and compare the simulation results to experimental data from a wind tunnel test of the Multirotor Test Bed (MTB). The MTB, developed by NASA Ames Research Center, is a new capability for testing a wide array of advanced vertical take-off and landing (VTOL) rotor configurations, with a primary focus on testing in the U.S. Army 7-by 10-Foot Wind Tunnel at NASA Ames Research Center. The MTB was designed to allow adjustment of the vertical, lateral, and longitudinal placement of up to six rotors, as well as allow tilt adjustment of each rotor and pitch adjustment of the whole assembly. The six-axis load cells under each rotor give the MTB the capability of measuring the rotor performance in a wide array of configurations. The overall goal of the MTB project is to help gain a better understanding of the performance, control, interactional aerodynamics, and acoustics of multirotor and tilting-rotor systems. For the work presented here, the MTB data were used to validate RotCFD and CHARM results for several multirotor test configurations. With confidence in both analyses established by the validation exercise, additional simulations were performed to explore quadrotor configurations that will be tested during the MTB’s second wind tunnel entry planned for 2022. This second tunnel entry will examine quadrotor configurations that represent published NASA reference designs for urban air mobility concept vehicles. Results from this paper confirm the ability of RotCFD and CHARM to simulate multirotor aerodynamic interactions on individual rotor performance under edgewise forward-flight conditions.
This paper summarizes findings from the sixth AIAA Propulsion Aerodynamics Workshop (PAW) nozzle test case. The focus of the workshop participants was turbulent supersonic round jets with variation in jet Mach number and jet heating. Three nozzles were utilized with on-design jet Mach numbers of 1.36, 1.63, and 2.00. A heated jet for each of these three nozzles at on-design conditions, a temperature matched on-design condition at Mach 1.63, and an off-design heated condition for the Mach 1.63 nozzle were the five cases examined. Workshop participants utilized Reynolds-averaged Navier-Stokes (RANS) and Large Eddy Simulation (LES) based techniques. Computational fluid dynamics (CFD) solutions were compared with non-intrusive measurements of velocities and temperatures, including both mean values and turbulence statistics. In general, the RANS approaches were unable to capture the trend of reduced mixing as the jet Mach number increased. Incorporation of compressibility corrections improved the agreement with experimental measurements in terms of trends with varying jet Mach number, but also increased the jet potential core lengths to be in worse agreement with data. RANS methods were able to reasonably capture trends of jet heating. The LES-based approaches as a whole were able to capture the effects of compressibility more accurately. Effects of jet heating were captured well by LES. None of the LES approaches introduced disturbances from the wall boundary layers of the nozzle into the jet mixing layer, which would affect the jet flow physics modeling in the initial part of the jet shear layer, but would be very computationally demanding.
Upgrades to the Mobile Launcher Rainbird nozzles were implemented for the Artemis II launch to increase water flow rates and improve the mobile launcher deck water coverage while minimizing the water throw distance to avoid impacting the Space Launch System (SLS) rocket nozzles. This analysis describes the use of a three-dimensional, multi-phase, free surface transient computational fluid dynamics (CFD) model to simulate the water flowing out through the Rainbird nozzle and ensure compliance with current Artemis requirements. Flow rates, pressures, velocities, nozzle throw distances and coverage areas were analyzed to ensure that system requirements are met. The completed results were used to make necessary modifications to the Ignition Overpressure Protection and Sound Suppression (IOP/SS) system Rainbird design by reshaping the nozzle used in the system.
Future cryogenic propulsion systems will require efficient methods with which to transfer cryogenic propellants from a depot storage tank to a customer receiver tank to minimize cost and maximize reusability. The Reduced Gravity Cryogenic Transfer project is currently developing advanced cryogenic fluid management technology and developing and validating new numerical models for three phases of transfer: line chilldown, tank chilldown, and tank fill. Additionally, multiple liquid nitrogen (LN 2 ) parabolic flight transfer rigs are being designed by universities and NASA to investigate the gravitational sensitivities that exist in these three technologies. In order to maximize the collection of low-g data during flights, it is required to extract as much (LN 2 as possible from the supply tank, despite variable gravity levels. The purpose of this paper is to present computational fluid dynamics (CFD) volume of fluid simulations of (LN 2 behavior in the supply tank onboard parabolic flights to validate the optimal design of a bi-directional propellant management device (PMD) using the commercial software FLOW-3D. A parametric study is conducted on the effects of gravity level, fill level, pore size, open area, thickness, and type of baffle on PMD performance. Based on results, the PMD as designed exceeds the targeted expulsion efficiency.
Outline - Introduction: X-57 CFD task overview; Motivation. Part I, Computational simulations without propulsion: Establishing CFD (Computational Fluid Dynamics) Best Practices - Grid generation - Mesh refinement study - Numerical methods - Wind tunnel validation study; Power-Off Aerodynamic Database Results. Part II, Computational simulations with propulsion: Cruise Power-On Database; High-Lift Power-On Database. Summary.
NASA has planned a supercomputer for computational fluid dynamics research since the mid-1970's. With the approval of the Numerical Aerodynamic Simulation Program as a FY 1984 new start, Congress requested an assessment of the program's objectives, projected short- and long-term uses, program design, computer architecture, user needs, and handling of proprietary and classified information. Specifically requested was an examination of the merits of proceeding with multiple high speed processor (HSP) systems contrasted with a single high speed processor system. The panel found NASA's objectives and projected uses sound and the projected distribution of users as realistic as possible at this stage. The multiple-HSP, whereby new, more powerful state-of-the-art HSP's would be integrated into a flexible network, was judged to present major advantages over any single HSP system.
Computational fluid dynamics (CFD) has the potential to improve the historical rocket injector design process by evaluating the sensitivity of performance and injector-driven thermal environments to the details of the injector geometry and key operational parameters. Methodical verification and validation efforts on a range of coaxial injector elements have shown the current production CFD capability must be improved in order to quantitatively impact the injector design process. This paper documents the status of a focused effort to compare and understand the predictive capabilities and computational requirements of a range of CFD methodologies on a set of single element injector model problems. The steady Reynolds-Average Navier-Stokes (RANS), unsteady Reynolds-Average Navier-Stokes (URANS) and three different approaches using the Large Eddy Simulation (LES) technique were used to simulate the initial model problem, a single element coaxial injector using gaseous oxygen and gaseous hydrogen propellants. While one high-fidelity LES result matches the experimental combustion chamber wall heat flux very well, there is no monotonic convergence to the data with increasing computational tool fidelity. Systematic evaluation of key flow field regions such as the flame zone, the head end recirculation zone and the downstream near wall zone has shed significant, though as of yet incomplete, light on the complex, underlying causes for the performance level of each technique. 1 Aerospace Engineer and Combustion CFD Team Leader, MS ER42, NASA MSFC, AL 35812, Senior Member, AIAA. 2 Professor and Director, Computational Combustion Laboratory, School of Aerospace Engineering, 270 Ferst Dr., Atlanta, GA 30332, Associate Fellow, AIAA. 3 Reilly Professor of Engineering, School of Mechanical Engineering, 585 Purdue Mall, West Lafayette, IN 47907, Fellow, AIAA. 4 Principal Member of Technical Staff, Combustion Research Facility, 7011 East Avenue, MS9051, Livermore, CA 94550, Associate Fellow, AIAA. 5 J. L. and G. H. McCain Endowed Chair, Mechanical Engineering, 104 Research Building East, University Park, PA 16802, Fellow, AIAA. American Institute of Aeronautics and Astronautics 1
The Wind-US computational fluid dynamics (CFD) flow solver was used to simulate dual-mode direct-connect ramjet/scramjet engine flowpath tests conducted in the University of Virginia (UVa) Supersonic Combustion Facility (SCF). The objective was to develop a computational capability within Wind-US to aid current hypersonic research and provide insight to flow as well as chemistry details that are not resolved by instruments available. Computational results are compared with experimental data to validate the accuracy of the numerical modeling. These results include two fuel-off non-reacting and eight fuel-on reacting cases with different equivalence ratios, split between one set with a clean (non-vitiated) air supply and the other set with a vitiated air supply (12 percent H2O vapor). The Peters and Rogg hydrogen-air chemical kinetics model was selected for the scramjet simulations. A limited sensitivity study was done to investigate the choice of turbulence model and inviscid flux scheme and led to the selection of the k-epsilon model and Harten, Lax and van Leer (for contact waves) (HLLC) scheme for general use. Simulation results show reasonably good agreement with experimental data and the overall vitiation effects were captured.
The thin-layer, Reynolds-averaged, Navier-Stokes equations are used to simulate the transonic viscous flow about the complete F-16A fighter aircraft. These computations demonstrate how computational fluid dynamics (CFD) can be used to simulate turbulent viscous flow about realistic aircraft geometries. A zonal grid approach is used to provide adequate viscous grid clustering on all aircraft surfaces. Zonal grids extend inside the F-16A inlet and up to the compressor face while power on conditions are modeled by employing a zonal grid extending from the exhaust nozzle to the far field. A simple solution adaptive grid procedure is used on the wing surface and good agreement with experimental data is obtained. Computations for the F-16A in side slip are also presented.