SEARCH · Search NASA

Results for “computational fluid dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS2018

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation in a Single Injector Cooling Flow

NASA's Turbulent Heat Flux (THX) task under the Transformational Tools & Technologies project is focused on acquiring benchmark temperature and velocity data in turbulent shear layers for validation of computational fluid dynamics (CFD) codes. Standard CFD turbulence models lack the ability to accurately calculate a number of fundamental flow phenomena, including the turbulent transport of heat. In nearly all production class Reynolds-averaged Navier-Stokes (RANS) CFD codes, a gradient-diffusion approximation is used whereby a constant turbulent Prantdl number, Prt, is used to relate an eddy viscosity calculated for the momentum terms, to the turbulent heat flux terms. The default in most CFD codes is Prt = 0.9, however this is not representative of all flows. As discussed in Yoder et al. (2015) and Reynolds (1974), Prt = 0.7 is felt to be more appropriate for jets. Recent work, as discussed in Yoder (2016a), has explored more complex formulations including variable turbulent Prantdl number models. PARTICLE IMAGE VELOCIMETRY; TEMPERATURE MEASURING INSTRUMENTS; FILM COOLING; HIGH TEMPERATURE GASES; COMPUTATIONAL FLUID DYNAMICS; PREDICTION ANALYSIS TECHNIQUES; TURBULENCE MODELS; RAMAN SPECTRA; VELOCITY; MACH NUMBER; GAS TURBINE ENGINES; NOZZLE FLOW; THERMOCOUPLES; LASER INDUCED FLUORESCENCE; COOLING FLOWS (ASTROPHYSICS); RAMAN SPECTROSCOPY; HEAT FLUX

Wernet, Mark P.↗

DOE OSTI2024

Successive Procedure for Solution Verification Based on User Needs

This paper discusses a revised solution verification procedure for computational fluid dynamics simulations to estimate the uncertainties in the quantities of interest based on discretization error models. This proposed procedure builds upon current procedures described in ASME V&V 20 but provides more guidance in determining the necessary number of mesh levels to build reliable discretization error models. Such guidance is particularly useful for practicing engineers without prior experience in solution verification. The key features of this proposed solution verification procedure are the ability to determine the need for additional mesh levels iteratively and the seamless treatment for underdetermined, exact, and overdetermined solutions of the power series approximation to the discretization error models. This study applies the proposed procedure to a set of synthetic examples to demonstrate the revised procedure’s clarity in determining the number of mesh solutions required for a reliable estimate of the discretization error in computational fluid dynamics settings. Additionally, this proposed procedure prevents a potential pathway in the current procedure in ASME V&V 20 that may lead to unreasonably small discretization errors.

Weinmeister, Justin↗

NASA NTRS2014

Multi-Point Combustion System: Final Report

A low-NOx emission combustor concept has been developed for NASA's Environmentally Responsible Aircraft (ERA) program to meet N+2 emissions goals for a 70,000 lb thrust engine application. These goals include 75 percent reduction of LTO NOx from CAEP6 standards without increasing CO, UHC, or smoke from that of current state of the art. An additional key factor in this work is to improve lean combustion stability over that of previous work performed on similar technology in the early 2000s. The purpose of this paper is to present the final report for the NASA contract. This work included the design, analysis, and test of a multi-point combustion system. All design work was based on the results of Computational Fluid Dynamics modeling with the end results tested on a medium pressure combustion rig at the UC and a medium pressure combustion rig at GRC. The theories behind the designs, results of analysis, and experimental test data will be discussed in this report. The combustion system consists of five radially staged rows of injectors, where ten small scale injectors are used in place of a single traditional nozzle. Major accomplishments of the current work include the design of a Multipoint Lean Direct Injection (MLDI) array and associated air blast and pilot fuel injectors, which is expected to meet or exceed the goal of a 75 percent reduction in LTO NOx from CAEP6 standards. This design incorporates a reduced number of injectors over previous multipoint designs, simplified and lightweight components, and a very compact combustor section. Additional outcomes of the program are validation that the design of these combustion systems can be aided by the use of Computational Fluid Dynamics to predict and reduce emissions. Furthermore, the staging of fuel through the individually controlled radially staged injector rows successfully demonstrated improved low power operability as well as improvements in emissions over previous multipoint designs. Additional comparison between Jet- A fuel and a hydrotreated biofuel is made to determine viability of the technology for use with alternative fuels. Finally, the operability of the array and associated nozzles proved to be very stable without requiring additional active or passive control systems. A number of publications have been publish

Goeke, Jerry↗

NASA NTRS1995

Cooperative Testing of Rocket Injectors That Use Gaseous Oxygen and Hydrogen

Gaseous oxygen and hydrogen propellants used in a special engine energy cycle called Full-Flow Staged Combustion are believed to significantly increase the lifetime of a rocket engine's pumps. The cycle can also reduce the operating temperatures of the engine. Improving the lifetime of the hardware reduces its overall maintenance and operations costs, and is critical to reducing costs for the joint NASA/industry Reusable Launch Vehicle (RLV). The work in this project will demonstrate the performance and lifetime of one-element and many-element combustors with gaseous O2/H2 injectors. This work supporting the RLV program is a cooperative venture of the NASA Lewis Research Center, the NASA Marshall Space Flight Center, Rocketdyne, and the Pennsylvania State University. Information about gas-gas rocket injector performance with O2/H2 is very limited. Because of this paucity of data, new testing is needed to improve the knowledge base for testing and designing new injectors for the RLV and to improve computer models that predict the combusting gas flows of new injector designs. Therefore, detailed observations and measurements of the combusting flow from many-element injectors in a rocket engine are being sought. These observations and measurements will be done with three different tools: schlieren photography, ultraviolet imaging, and Raman spectroscopy. The schlieren system will take photos of the density differences in combusting flow, the ultraviolet movies will determine the location of the hydroxyl (OH) radical in the combustion flow, and the Raman spectroscopic measurements will provide the combustion temperature and amount of water (H2O), hydrogen (H2), and oxygen (O2) in the combustor. Marshall is providing overall program management, design and computational fluid dynamics (CFD) analyses, as well as funding for the work at Penn State. An existing, windowed combustor and several injectors will be provided by Rocketdyne--two injectors for the initial screening tests and one with an optimized design based on the best design found in the screening tests. Lewis will provide a nozzle and several injectors for the screening test program. The configuration of the injectors will be based on a design chosen by all the participants, and their elements will be based on the coaxial and impinging flow. Lewis also will provide the instrumentation for the flow-field measurements: schlieren, ultraviolet imaging, and Raman spectroscopy. In addition, thermocouples will measure heat flow on the injector face. Other traditional measurements of rocket performance will be made as well: chamber pressure, mass flow of each propellant, purge flow, and the barrier cooling gas flow. Penn State will conduct single-element testing with the injector elements from both the Rocketdyne and the jointly designed injectors. A wide variety of traditional and nontraditional injector designs will be tested in this program. The results will be valuable in computational fluid dynamics code validation and overall rocket combustion efficiency measurements. Correlations between combustion efficiency, laser measurements of species, and ultraviolet and visible light photography will also be made. Thus far, several different single-element injectors have been tested at Penn State and Lewis. The experimental setup of a rocket engine with a viewing window is shown. The combusting flow is shown. The results are helping engineers design the many element injectors.

Source record↗

NASA NTRS2002

Tip-Clearance Vortex Characterized With Three-Dimensional Digital Particle Image Velocimetry

An optical measurement technique known as Three-Dimensional Digital Particle Image Velocimetry (3-D DPIV) was used to characterize the tip clearance flow in NASA Glenn Research Center's low-speed axial compressor. 3-D DPIV is a technique in which a stereoscopic imaging system consisting of two cross-correlation cameras is used to record particles entrained in a flow as a laser light sheet is pulsed at two instances in time. Although 3-D DPIV has been used elsewhere, this is the first time it has been used to measure compressor tip clearance flows. In-house modifications of the DPIV system include the use of effective seeding technology and a novel system to perform a priori calibrations at all five measurement planes, greatly reducing facility run time. Computational fluid dynamics predictions, which are used to guide design changes toward improving the efficiency and operating range of turbomachinery, can be verified and improved by comparison with 3-D DPIV measurements of the actual tip clearance flow. This measurement campaign dealt with the characterization of the tip clearance vortex in the first stage of a four-stage axial compressor. The tip clearance vortex is formed in compressors operating with a clearance gap between the moving rotor blade tips and the stationary casing when a leakage flow, forced from the pressure side of the blade over the blade tip, forms a vortical structure on the suction side of the blade. 3-D DPIV is ideally suited to measure the clearance vortex for two reasons: (1) this technique captures the entire blade passage flow at one instant in time, so that wandering of the vortex during the measurement does not smear out velocity gradients in the flow field, and (2) the spanwise component of velocity changes sign across the vortex core, providing a more accurate measurement of the vortex location than was available with previous two-dimensional measurement approaches. These two attributes of the data will enable computational fluid dynamics researchers to validate their predictions to a level of accuracy not previously attainable. In an effort to understand the effects of this flow phenomenon on the operation of the low speed axial compressor, data were acquired at two mass flow coefficients: 0.395 (design operating point) and 0.35 (operating point just above stall). In order to identify the spatial extent, location, and magnitude of the tip clearance vortex for each mass flow condition, data were acquired at five equally spaced spanwise locations, from 90- to 100-percent span. The data presented shows a cross section of the tip clearance vortex at 94-percent span and a mass flow coefficient of 0.395, where the color contours represent the radial velocity component. The blue (flow towards the hub) and red (flow towards the casing) regions illustrate the extent of the tip clearance vortex. The interface between the blue and red regions delineates the location of the vortex core.

John, W. Trevor↗

NASA NTRS2009

A CFD Approach to Modeling Spacecraft Fuel Slosh

Energy dissipation and resonant coupling from sloshing fuel in spacecraft fuel tanks is a problem that occurs in the design of many spacecraft. In the case of a spin stabilized spacecraft, this energy dissipation can cause a growth in the spacecrafts' nutation (wobble) that may lead to disastrous consequences for the mission. Even in non-spinning spacecraft, coupling between the spacecraft or upper stage flight control system and an unanticipated slosh resonance can result in catastrophe. By using a Computational Fluid Dynamics (CFD) solver such as Fluent, a model for this fuel slosh can be created. The accuracy of the model must be tested by comparing its results to an experimental test case. Such a model will allow for the variation of many different parameters such as fluid viscosity and gravitational field, yielding a deeper understanding of spacecraft slosh dynamics. In order to gain a better understanding of the dynamics behind sloshing fluids, the Launch Services Program (LSP) at the NASA Kennedy Space Center (KSC) is interested in finding ways to better model this behavior. Thanks to past research, a state-of-the-art fuel slosh research facility was designed and fabricated at Embry Riddle Aeronautical University (ERAU). This test facility has produced interesting results and a fairly reliable parameter estimation process to predict the necessary values that accurately characterize a mechanical pendulum analog model. The current study at ERAU uses a different approach to model the free surface sloshing of liquid in a spherical tank using Computational Fluid Dynamics (CFD) methods. Using a software package called Fluent, a model was created to simulate the sloshing motion of the propellant. This finite volume program uses a technique called the Volume of Fluid (VOF) method to model the interaction between two fluids [4]. For the case of free surface slosh, the two fluids are the propellant and air. As the fuel sloshes around in the tank, it naturally displaces the air. Using the conservation of mass, momentum, and energy equations, as well as the VOF equations, one can predict the behavior of the sloshing fluid and calculate the forces, pressure gradients, and velocity field for the entire liquid as a function of time.

Marsell, Brandon↗

DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗

DOE OSTI2024

Accelerating transients with NekRS: GPU overlapping domain implementation and multi-rate timestepping

The simulation of nuclear transients using Computational Fluid Dynamics (CFD) presents significant computational challenges due to the inherent complexity and the wide separation in temporal scales between various flow physical phenomena. These disparities lead to high computational costs, often making the simulation of transients impractical without advanced techniques. Consequently, multiple research initiatives are being pursued by the NEAMS thermal-hydraulic area, some driven by academic institutions and some by national laboratories. Overall, they are exploring novel methods to make transient simulations more feasible and efficient. This report delves into recent advancements within the CFD code NekRS, specifically those achieved in Fiscal Year 2024 under the CONNECT effort, aimed at improving the performance and feasibility of transient simulations. The first major advancement involves the porting of NekRS to Aurora, one of the Department of Energy’s (DOE) most powerful supercomputers. Additionally, the report discusses the implementation of an overlapping domain capability within NekRS. This novel GPU-accelerated capability allows different spatial regions of the domain to be solved independently, enhancing the code’s efficiency, particularly when running large-scale simulations in complex domains. The scalability of this approach is demonstrated, highlighting its potential to transform how transients are approached in CFD simulations. Lastly, the report focuses on how this overlapping domain capability specifically accelerates transient simulations through multi-rate timestepping. By decoupling different regions and facilitating faster computations, this method offers a promising pathway to making nuclear transient simulations more computationally feasible, addressing one of the critical bottlenecks in the field. Together, these advancements represent a significant leap forward in transient simulation technology, bringing closer the possibility of handling highly complex nuclear scenarios with greater efficiency and accuracy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Computational Fluid Dynamics-Based Modeling of Methane Flows Around Oil and Gas Equipment

Recent studies estimate that emissions from oil and gas production facilities contribute between 20 and 50% of the total methane ( CH 4 ) emitted in the US; therefore, quantifying and reducing these emissions are crucial for achieving climate goals. Methane quantification depends on both measuring methane concentrations and converting them to emissions through a modeling framework. Currently, simple atmospheric dispersion models are primarily used to quantify emissions and concentrations, but these estimates are highly uncertain when quantifying emissions from complex aerodynamic sources, such as oil and gas facilities. This investigation used a CFD modeling approach, which can account for aerodynamic complexity but has hitherto not been used to model methane concentrations downwind of a methane release of a known rate, and compared it against in situ measurements. High-time-resolution (1 Hz) methane concentration and meteorological data were measured during experiments conducted at the METEC on 21 March and 11 July 2024. The METEC site configuration, measured wind data, and controlled emission rates were used as input for the CONVERGE CFD model to model downwind CH 4 concentration. The modeling was carried out between 20 and 70 m, from two different points of release in two separate controlled-release experiments, one from a separator and another from a wellhead. In these experiments, we found that the CFD model could predict the CH 4 concentrations downwind of the release to a good degree. The model was evaluated on multiple metrics to assess its performance in estimating methane concentrations at typical fence line distances (∼30 m). These results help us to understand external flows and the ability of CFD models to predict downwind concentrations in aerodynamically complex environments.

03 NATURAL GAS↗

NASA NTRS1979

Multi-level adaptive computations in fluid dynamics

The multi-level adaptive technique (MLAT) is a general strategy of solving continuous problems by cycling between coarser and finer levels of discretization. It provides very fast solvers together with adaptive, nearly optimal discretization schemes to general boundary-value problems in general domains. Here the state of the art is surveyed, emphasizing steady-state fluid dynamics applications, from slow viscous flows to transonic ones. Various new techniques are briefly discussed, including distributive relaxation schemes, the treatment of evolution problems, the combined use of upstream and central differencing, local truncation extrapolations, and other 'super-solver' techniques.

Brandt, A.↗

NASA NTRS1986

Perturbation Method for Computational Fluid-Dynamical Equations

Perturbation technique yields accurate flow solutions using as few as one-fourth number of grid points required by finite-difference methods. Technique originally developed to solve Euler equations of two-dimensional, steady, inviscid transonic flow about airfoils, applicable to arbitrary equation sets and higher dimensions. New perturbations scheme used in design cycle where potential solutions generated routinely; Euler perturbation method used in second-cut analysis. Method also used to couple other equation sets.

Chow, L. J.↗

NASA NTRS1986

The Space Shuttle, fluid dynamics, and computer graphics

Several samples are presented of CFD-generated flow ribbon images of fuel flow in the SSME. The images result from simulations of H2 flow through a design iteration of the main injector, the hot gas manifold and the fuel bowl section. The color-coded ribbons provide three-dimensional perspectives which indicate the direction and orientation of the flow and the velocity magnitude in intervals along the flow path.

Belie, R. G.↗

NASA NTRS1989

Experimental and computational flow-field results for an all-body hypersonic aircraft

A comprehensive test program is defined which is being implemented in the NASA/Ames 3.5 foot Hypersonic Wind Tunnel for obtaining data on a generic all-body hypersonic vehicle for computational fluid dynamics (CFD) code validation. Computational methods (approximate inviscid methods and an upwind parabolized Navier-Stokes code) currently being applied to the all-body model are outlined. Experimental and computational results on surface pressure distributions and Pitot-pressure surveys for the basic sharp-nose model (without control surfaces) at a free-stream Mach number of 7 are presented.

Cleary, Joseph W.↗

NASA NTRS1990

The calculation of rotor/fuselage interaction for two-dimensional bodies

Unsteady rotor wake interactions with the empennage, tail boom, and other aerodynamic surfaces have a significant influence on the aerodynamic performance of the helicopter, ride quality, and vibration. A Computational Fluid Dynamic (CFD) method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary 2-D bodies was developed to address this helicopter problem. The vorticity and flow field velocities are calculated on a body-fitted computational mesh using an uncoupled iterative solution. The interacting vortex wake is represented by an array of discrete vortices which, in turn, are represented by a finite core model. The evolution of the interacting vortex wake is calculated by Lagrangian techniques. The flow around circular and elliptic cylinders in the absence of an interacting vortex wake was calculated. These results compare very well with other numerical results and with results obtained from experiment and thereby demonstrate the accuracy of the viscous solution. The interaction of a simulated rotor wake with the flow about 2-D bodies, representing cross sections of fuselage components, was calculated to address the vortex interaction problem. The vortex interaction was calculated for the flow about a circular and an elliptic cylinder at 45 and 90 degrees incidence. The results demonstrate the significant variation in lift and drag on the 2-D bodies during the vortex interaction.

Stremel, Paul M.↗

NASA NTRS1990

A finite element solver for 3-D compressible viscous flows

Computation of the flow field inside a space shuttle main engine (SSME) requires the application of state of the art computational fluid dynamic (CFD) technology. Several computer codes are under development to solve 3-D flow through the hot gas manifold. Some algorithms were designed to solve the unsteady compressible Navier-Stokes equations, either by implicit or explicit factorization methods, using several hundred or thousands of time steps to reach a steady state solution. A new iterative algorithm is being developed for the solution of the implicit finite element equations without assembling global matrices. It is an efficient iteration scheme based on a modified nonlinear Gauss-Seidel iteration with symmetric sweeps. The algorithm is analyzed for a model equation and is shown to be unconditionally stable. Results from a series of test problems are presented. The finite element code was tested for couette flow, which is flow under a pressure gradient between two parallel plates in relative motion. Another problem that was solved is viscous laminar flow over a flat plate. The general 3-D finite element code was used to compute the flow in an axisymmetric turnaround duct at low Mach numbers.

Reddy, K. C.↗

NASA NTRS1995

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.↗

NASA NTRS1994

CFD Optimization on Network-Based Parallel Computer System

Combining multiple engineering workstations into a network-based heterogeneous parallel computer allows application of aerodynamic optimization with advance computational fluid dynamics codes, which is computationally expensive in mainframe supercomputer. This paper introduces a nonlinear quasi-Newton optimizer designed for this network-based heterogeneous parallel computer on a software called Parallel Virtual Machine. This paper will introduce the methodology behind coupling a Parabolized Navier-Stokes flow solver to the nonlinear optimizer. This parallel optimization package has been applied to reduce the wave drag of a body of revolution and a wing/body configuration with results of 5% to 6% drag reduction.

Cheung, Samson H.↗

NASA NTRS1994

CFD Optimization on Network-Based Parallel Computer System

Cheung, Samson H.↗