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 1,225 records · Page 68

Coupled CFD-PBE Predictions of Renal Stone Size Distributions in the Nephron in Microgravity

In this paper, a deterministic model is developed to assess the risk of critical renal stone formation for astronauts during space travel. A Population Balance Equation (PBE) model is used to compute the size distribution of a population of nucleating, growing and agglomerating renal calculi as they are transported through different sections of the nephron. The PBE model is coupled to a Computational Fluid Dynamics (CFD) model that solves for steady state flow of urine and transport of renal calculi along with the concentrations of ionic species, calcium and oxalate, in the nephron using an Eulerian two-phase mathematical framework. Parametric simulation are performed to study stone size enhancement and steady state volume fraction distributions in the four main sections of the nephron under weightlessness conditions. Contribution of agglomeration to the stone size distribution and effect of wall friction on the stone volume fraction distributions are carefully examined. Case studies using measured astronaut urinary calcium and oxalate concentrations in microgravity as input indicate that under nominal conditions the largest stone sizes developed in Space will be still considerably below the critical range for problematic stone development. However, results also indicate that the highest stone volume fraction occurs next to the tubule and duct walls. This suggests that there is an increased potential for wall adhesion with the possibility of evolution towards critical stone sizes.

Crystal Growth↗

The Isolated Synthetic Jet in Crossflow: A Benchmark for Flow Control Simulation

An overview of the data acquisition, reduction, and uncertainty of experimental measurements made of the flowfield created by the interaction of an isolated synthetic jet and a turbulent boundary layer is presented. The experimental measurements were undertaken to serve as the second of three computational fluid dynamics validation databases for Active Flow Control. The validation databases were presented at the NASA Langley Research Center Workshop on CFD Validation of Synthetic Jets and Turbulent Separation Control in March, 2004. Detailed measurements were made to document the boundary conditions for the flow and also for the phase-averaged flowfield itself. Three component Laser-Doppler Velocimetry, 2-D Particle Image Velocimetry, and Stereo Particle Image Velocimetry were utilized to document the phase-averaged velocity field and the turbulent stresses.

Schaeffler, Norman W.↗

Turbulence Model Validation through Joint Experimental / Computational Studies of Separated Flow Over a Three-Dimensional Tapered Bump: Part II - Computational Investigation

This report serves as documentation of activities related to Computational Fluid Dynamics (CFD) analysis applied to the Smooth Body Separation Experiment designed to produce validation-quality experimental data from wind-tunnel tests performed at the University of Notre Dame. All CFD analysis utilized the Boeing CFD (BCFD) flow solver. BCFD was used to analyze flow over the Boeing Speed Bump with second-order accuracy in space and time. Inviscid fluxes were computed using the HLLE scheme for steady-state simulations, and the Low-Dissipation Roe scheme for time-accurate simulations. The working fluid was modeled as calorically perfect air using Sutherland's law for viscosity. Thermal diffusivity was determined by assuming a constant molecular Prandtl number of 0.72, and a turbulent Prandtl number of 0.9. Reynolds Averaged Navier-Stokes (RANS) simulations were conducted using the one-equation Spalart-Allmaras turbulence model with rotation correction3 and the quadratic constitutive relation (SA-RC-QCR), as well as Menter's two-equation, Shear-Stress Transport (SST) model. Time-accurate computations were conducted using Delayed Detached eddy Simulation (DDES) with the one-equation RANS model, and Large Eddy Stress balance (LESb) with the two-equation RANS model. Implicit dual-time stepping was used to achieve second-order temporal accuracy. The experimental investigation conducted during the course of this study is described in detail in an accompanying report.

Flow separation↗

Time-Accurate Computational Fluid Dynamics Analysis of a Supersonic Turbine to Assess Structural Mode Excitation

Gas turbines can experience significant unsteady loading during operation which can cause damage such as blade cracking or bearing wear. While computational fluid dynamics (CFD) modeling of full three-dimensional gas turbine flow to verify meanline design performance is standard practice within the industry, the use of time-accurate CFD makes it possible to assess the unsteady loading which can be especially problematic at operating conditions where unsteady fluid loading aligns with structural modes. The Fluid Dynamics branch (ER42) at NASA Marshall Space Flight Center has developed a time-accurate CFD analysis methodology using the unstructured, density-based, Loci-CHEM solver which is capable of simulating the relative motion required for turbomachinery applications through sliding interface meshes. The highly-parallel unstructured solver allows for full three-dimensional, time-accurate simulation of turbines including upstream and downstream components such as manifolds, nozzles, and stators, any of which could excite a turbine structural mode. Recently, Loci-CHEM was used to analyze the FASTRAC supersonic turbine to assess the loading on the rotor. Overall performance and time-averaged loading on the rotor was computed, including the axial force over the entire rotor disk which is required for force balance calculations on the turbopump. The FASTRAC turbine has supersonic nozzles upstream of the rotor that result in significant unsteady loading on the rotor blades at multiples of the nozzle pass frequency. The unsteady blade pressure loading was decomposed into the frequency domain for direct use in structural forced response analysis. Utilizing this type of CFD analysis coupled with structural analysis has proven to be a valuable tool in identifying structural modes that are being excited by unsteady fluid loading.

M.V. Fischels↗

Time-Accurate Computational Fluid Dynamics Analysis of a Supersonic Turbine to Assess Structural Mode Excitation

Gas turbines can experience significant unsteady loading during operation which can cause damage such as blade cracking or bearing wear. While computational fluid dynamics (CFD) modeling of full three-dimensional gas turbine flow to verify meanline design performance is standard practice within the industry, the use of time-accurate CFD makes it possible to assess the unsteady loading which can be especially problematic at operating conditions where unsteady fluid loading aligns with structural modes. The Fluid Dynamics branch (ER42) at NASA Marshall Space Flight Center has developed a time-accurate CFD analysis methodology using the unstructured, density-based, Loci-CHEM solver which is capable of simulating the relative motion required for turbomachinery applications through sliding interface meshes. The highly-parallel unstructured solver allows for full three-dimensional, time-accurate simulation of turbines including upstream and downstream components such as manifolds, nozzles, and stators, any of which could excite a turbine structural mode. Recently, Loci-CHEM was used to analyze the FASTRAC supersonic turbine to assess the loading on the rotor. Overall performance and time-averaged loading on the rotor was computed, including the axial force over the entire rotor disk which is required for force balance calculations on the turbopump. The FASTRAC turbine has supersonic nozzles upstream of the rotor that result in significant unsteady loading on the rotor blades at multiples of the nozzle pass frequency. The unsteady blade pressure loading was decomposed into the frequency domain for direct use in structural forced response analysis. Utilizing this type of CFD analysis coupled with structural analysis has proven to be a valuable tool in identifying structural modes that are being excited by unsteady fluid loading.

M V Fischels↗

CFD simulation of coaxial injectors

The development of improved performance models for the Space Shuttle Main Engine (SSME) is an important, ongoing program at NASA MSFC. These models allow prediction of overall system performance, as well as analysis of run-time anomalies which might adversely affect engine performance or safety. Due to the complexity of the flow fields associated with the SSME, NASA has increasingly turned to Computational Fluid Dynamics (CFD) techniques as modeling tools. An important component of the SSME system is the fuel preburner, which consists of a cylindrical chamber with a plate containing 264 coaxial injector elements at one end. A fuel rich mixture of gaseous hydrogen and liquid oxygen is injected and combusted in the chamber. This process preheats the hydrogen fuel before it enters the main combustion chamber, powers the hydrogen turbo-pump, and provides a heat dump for nozzle cooling. Issues of interest include the temperature and pressure fields at the turbine inlet and the thermal compatibility between the preburner chamber and injector plate. Performance anomalies can occur due to incomplete combustion, blocked injector ports, etc. The performance model should include the capability to simulate the effects of these anomalies. The current approach to the numerical simulation of the SSME fuel preburner flow field is to use a global model based on the MSFC sponsored FNDS code. This code does not have the capabilities of modeling several aspects of the problem such as detailed modeling of the coaxial injectors. Therefore, an effort has been initiated to develop a detailed simulation of the preburner coaxial injectors and provide gas phase boundary conditions just downstream of the injector face as input to the FDNS code. This simulation should include three-dimensional geometric effects such as proximity of injectors to baffles and chamber walls and interaction between injectors. This report describes an investigation into the numerical simulation of GH2/LOX coaxial injectors. The following sections will discuss the physical aspects of injectors, the CFD code employed, and preliminary results of a simulation of a single coaxial injector for which experimental data is available. It is hoped that this work will lay the foundation for the development of a unique and useful tool to support the SSME program.

Landrum, D. Brian↗

Down-to-earth supercomputing at NASA Ames

The requirements for computers for aerospace vehicle design problems, in particular fluid dynamics, are examined. The levels of approximation to Navier-Stokes equations are: (1) linearized inviscid form, (2) nonlinear inviscid equation, (3) Reynolds-averaged Navier-Stokes equation, (4) large eddy simulation with models for subgrid scale turbulence, and (5) full Navier-Stokes equations to calculate laminar/turbulent transition and turbulence dissipation. The Cray 2 with 256 megaword memory and 4.1 nsec clock and the Cray X-MP/48 with a 128 megaword memory and a 9.5 nsec clock of the Ames facility are described, and an example revealing their applicability is presented. The use of zonal methods to simulate aerodynamic flow fields numerically is discussed. Consideration is given to the development of a program for simulation of aerodynamic flow fields, the compilers, symbolic debuggers, the scheduling of multiple CPUs, and the analysis and storage of the data.

Stevens, Kenneth G., Jr.↗

Design and Analysis of Turbomachinery for Space Applications

This presentation provides an overview of CORSAIR, a three dimensional computational fluid dynamics software code for the analysis of turbomachinery components available from NASA, and discusses its potential use in the design of these parts. Topics covered include: time-dependent equations of motion, grid topology, turbulence models, boundary conditions, parallel simulations and miscellaneous capabilities.

Dorney, D.↗

Computational Investigation of Powered Descent for Human-Scale Mars Landers

NASA human-scale Mars entry, descent, and landing vehicles currently under consideration have system requirements beyond the current state of the art. Specifically, heritage approaches for deceleration using parachutes and initiating engines at subsonic speeds are not adequate to land these large vehicles (50-60 metric ton entry mass) precisely (within 50 meters of a surface target). Therefore, the present approach for landing humans on Mars utilizes entry vehicles with large aerodynamic surface areas that transition to the descent and landing phase by initiating retropropulsion at supersonic speeds. Vehicle performance, including stability and control, depends on accurate characterization of the aerodynamic-propulsive interference effects, especially at engine initiation, when the vehicles are transitioning from an entry attitude to that suitable for descent and landing. This paper summarizes the human Mars entry vehicle configurations, retropropulsion integration, performance assumptions, and expected flight environment at engine initiation. Cases spanning supersonic, transonic, and subsonic conditions during powered descent are analyzed using four computational fluid dynamics codes, all of which have been previously applied to supersonic retropropulsion at varying scales and conditions. The results of the analysis are used to develop an aerodynamic-propulsive interference model that has since been applied in trajectory simulations to augment the fidelity of human-scale Mars entry vehicle performance.

Ashley M Korzun↗

X-57 “Maxwell” High-Lift Propeller Testing and Model Development

NASA’s X-57 “Maxwell” distributed electric propulsion flight demonstrator has a high-lift system that includes 12 fixed-pitch high-lift propellers located upstream of the wing leading edge for lift augmentation at low speeds. These high-lift propellers are not required at higher speeds and are folded conformally along the nacelles to reduce drag when not in operation. Aircraft performance models and flight simulations incorporate propeller performance to predict thrust, moments, and power consumption and propeller model accuracy is important in identifying safe operating regimes for this aircraft. Current high-lift propeller performance models have been verified and calibrated against numerous computational fluid dynamics analyses under a variety of flight conditions. We performed a series of full-scale wind tunnel tests at the NASA Langley Research Center Low Speed Aeroacoustic Wind Tunnel to further validate these models, to identify any adverse operating conditions for these propellers, and to assess the accuracy of the facility’s new Propeller Test Stand. The results indicate that the models accurately predicted performance and that the right-handed propeller showed lower torque and thrust for a given propeller speed compared to the left-handed version. Acoustic data support these measurements, showing higher tonal noise for the left versus right propeller. We suspect this is due to slight differences in the material properties between the early production (right) and more recent (left) blade sets leading to greater detwisting in the right-handed blades under load. Both propellers demonstrated very stable operation throughout the test including during deployment, stowing, and windmilling. Furthermore, high-speed imagery of the spinning propeller indicated no blade oscillations or instabilities for any test case. We also leveraged these images to analyze the effects of propeller loading on blade detwisting and tip position. Unfortunately, the thrust measurements varied significantly during the test due to cross-loading of the Propeller Test Stand load cell, though trends in the thrust data were still observable. The torque measurement, however, remained very precise for all test points and was used for performance comparisons. As a result of this test, we validated and improved the X-57 high-lift propeller performance models, obtained detailed propeller acoustic data and high-speed imagery, observed the propeller operating safely during a variety of dynamic events, and identified two potentially adverse, low-propeller-speed conditions which will be avoided during normal aircraft operation.

X-57↗

Structure and Dynamics of Bleed-Controlled Impinging Shock/Turbulent-Boundary-Layer Interactions

Porous boundary-layer bleed is often employed to alleviate the detrimental effects of shock/turbulent-boundary-layer interactions (STBLI) in high-speed vehicles. However, outstanding questions remain on the underlying mechanisms through which bleed influences STBLI dynamics. This work addresses this knowledge gap by performing a scale-resolved Large-Eddy Simulation of an impinging STBLI controlled by discrete bleed holes. The simulation has a freestream Mach number of 2.5, an oblique shock generator angle of 8 degrees, and a Reynolds number based on the incoming momentum thickness of 2800. Bleed substantially reduces the separation size and upstream influence of the STBLI compared to the canonical, uncontrolled case. Consequently, the characteristic, low-frequency motion of STBLI-induced separation is mitigated, as revealed by wall-pressure spectra. Energy is shifted from the low frequency separation behavior to a higher-frequency signature caused by the bleed shocks and expansions. The scale-resolved simulations also reveal that these small, localized structures severely limit the turbulence production traditionally seen in the uncontrolled interaction.

computational fluid dynamics↗

Evaluation of CFD as a Surrogate for Mach 2.4 to 4.6 Wind-Tunnel Testing – Project Overview

The debate over when wind-tunnel testing will be replaced by Computational Fluid Dynamics (CFD) comes and goes. More recently the debate has subsided with a more collaborative spirit between practitioners of these two disciplines resulting in significant improvements in the outcomes of both. There may come a time, however, when CFD has sufficient accuracy to supplant WTT as the dominant or perhaps only tool for aerodynamic simulation. If and/or when that happens, financial pressures favor efforts to close or severely limit the operations of wind tunnels. Presumably additional resources will go toward CFD to generate aerodynamic databases, load environments, and new aero/fluid-dynamic knowledge. It is therefore important to develop appropriate processes by which wind-tunnel closure decisions are made to ensure that facilities critical to industry and government research and development aren’t closed prematurely without proof that the available CFD tools have sufficient accuracy and low-enough cost (and enough experts and computational facilities) to take on the traditional role of wind tunnels. This paper will describe a project intended to answer the specific question of whether CFD can replace wind-tunnel testing for the limited Mach-number range 2.4 to 4.6. The project involves wind-tunnel testing and coordinated CFD for a variety of vehicle and flow-physics types in the high-speed leg of the Unitary Plan Wind Tunnel at NASA’s Langley Research Center.

James C Ross↗

Evaluation of CFD as a Surrogate for Wind-Tunnel Testing for Mach 2.4 to 4.6 - Project Overview

The debate over when wind-tunnel testing (WTT) will be replaced by Computational Fluid Dynamics (CFD) comes and goes. More recently the debate has subsided with a more collaborative spirit between practitioners of these two disciplines resulting in significant improvements in the outcomes of both. There may come a time, however, when CFD has sufficient accuracy to supplant WTT as the dominant or perhaps only tool for aerodynamic simulation. If and/or when that happens, financial pressure result in efforts to close or severely limit the operations of wind tunnels. Presumably additional resources will go toward CFD in order to generate aerodynamic databases, load environments, and new aero/fluid-dynamic knowledge. It is therefore important to develop appropriate processes by which wind-tunnel closure decisions are made to ensure that facilities critical to industry and government research and development aren’t closed prematurely without ensuring that the existing CFD tools have sufficient accuracy and low-enough cost (and enough experts and computational facilities) to take on the traditional role of wind tunnels. This paper will describe a project intended to answer the specific question of whether CFD can replace WTT for the limited Mach-number range 2.4 to 4.6. The wind tunnel being examined in this context is the high-speed leg of the Unitary Plan Wind Tunnel at NASA’s Langley Research Center (LaRC UPWT).

James C. Ross↗

Low-Dissipation Advection Schemes Designed for Large Eddy Simulations of Hypersonic Propulsion Systems

The 2nd-order upwind inviscid flux scheme implemented in the multi-block, structured grid, cell centered, finite volume, high-speed reacting flow code VULCAN has been modified to reduce numerical dissipation. This modification was motivated by the desire to improve the codes ability to perform large eddy simulations. The reduction in dissipation was accomplished through a hybridization of non-dissipative and dissipative discontinuity-capturing advection schemes that reduces numerical dissipation while maintaining the ability to capture shocks. A methodology for constructing hybrid-advection schemes that blends nondissipative fluxes consisting of linear combinations of divergence and product rule forms discretized using 4th-order symmetric operators, with dissipative, 3rd or 4th-order reconstruction based upwind flux schemes was developed and implemented. A series of benchmark problems with increasing spatial and fluid dynamical complexity were utilized to examine the ability of the candidate schemes to resolve and propagate structures typical of turbulent flow, their discontinuity capturing capability and their robustness. A realistic geometry typical of a high-speed propulsion system flowpath was computed using the most promising of the examined schemes and was compared with available experimental data to demonstrate simulation fidelity.

White, Jeffrey A.↗

Determination of solid/porous wall boundary conditions from wind tunnel data for computational fluid dynamics codes

A computational and experimental study has been undertaken to investigate methods of modeling solid and porous wall boundary conditions in CFD codes. The procedure utilizes experimental measurements at the walls to develop a flow field solution based on the method of singularities. This flow field solution is then imposed as a boundary condition in a CFD simulation of the internal flow field. The effectiveness of this method in describing the boundary conditions at the wind tunnel walls using only sparse experimental measurements has been investigated. Verification of the approach using computational studies has been carried out using an incompressible flow solver. The current work demonstrates this technique for low speed flows and compares the result with experimental data obtained from a heavily instrumented variable porosity test section. Position and refinement of experimental measurements required to describe porous wall boundary conditions has also been considered for application to other porous wall wind tunnels.

Beutner, Thomas J.↗

Quasi 1D Modeling of Mixed Compression Supersonic Inlets

The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of the 2-dimensional bifurcated mixed compression supersonic inlet is being developed. The model utilizes computational fluid dynamics for both the supersonic and subsonic diffusers. The oblique shocks are modeled utilizing compressible flow equations. This model also implements variable geometry required to control the normal shock position. The model is flexible and can also be utilized to simulate other mixed compression supersonic inlet designs. The model was validated both in time and in the frequency domain against the legacy LArge Perturbation INlet code, which has been previously verified using test data. This legacy code written in FORTRAN is quite extensive and complex in terms of the amount of software and number of subroutines. Further, the legacy code is not suitable for closed loop feedback controls design, and the simulation environment is not amenable to systems integration. Therefore, a solution is to develop an innovative, more simplified, mixed compression inlet model with the same steady state and dynamic performance as the legacy code that also can be used for controls design. The new nonlinear dynamic model is implemented in MATLAB Simulink. This environment allows easier development of linear models for controls design for shock positioning. The new model is also well suited for integration with a propulsion system model to study inlet/propulsion system performance, and integration with an aero-servo-elastic system model to study integrated vehicle ride quality, vehicle stability, and efficiency.

Kopasakis, George↗

Calorimeter Heat Flux Trends in NASA’s Subscale Rotating Detonation Rocket Engine

A calorimeter-style rotating detonation rocket engine (RDRE) outer chamber body was designed, additively manufactured in GRCop-42, and hot-fire tested in test cell 105 at NASA Marshall Space Flight Center (MSFC) in Huntsville, AL in August 2024. Steady-state, spatially resolved heat flux and chamber pressure were measured at 13 and 9 locations from the injector face, respectively, for a variety of test conditions and hardware configurations with gaseous methane and oxygen propellants. The 2-in. diameter highly instrumented design, enabled by metal additive manufacturing (AM), allowed for simultaneous measurement of heat flux and chamber pressure through 44 long duration hot-fire tests without failure. Of the 44 tests, only 3 of the tests featured chamber pressures below 100 psia and are discussed in this work. Heat flux trends for this subscale hardware are explored with respect to local propellant mass flux for several different injectors, chamber lengths, and chamber “subsonic” area ratios (contraction ratios). Results show that for increasing chamber pressure, a reduced amount of energy per kilogram of propellant is transferred to the wall. Switching injectors affects the near-face heat flux and chamber pressure profiles. Chamber length does not noticeably affect the heat flux profile but increases bulk heat load. A subsonic area ratio greater than 1 reduces heat losses to the wall compared to a straight annulus case. Equivalence ratio shift did not significantly affect heat flux. Computational fluid dynamics (CFD) and conjugate heat transfer (CHT) analyses were performed to delineate 3-D heat transfer and coolant mass flow maldistribution effects and create a calorimeter transfer function to transform heat flux data and mitigate the profile distortions. Additional CFD/CHT simulations were performed using experimental data as boundary conditions to assess local nucleate boiling propensity during testing and to assess projected hot wall temperatures for future fatigue assessment. Vapor fractions were observed locally for several cases. It is anticipated that this effort will provide a common framework for coolant system design/analysis in calorimeter systems for rocket-based and air-breathing combustors, ensuring greater repeatability in calorimetry measurements across the research community.

Joseph Hernandez-McCloskey↗