SEARCH · Search NASA

Results for “Computational fluid dynamics code”

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 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Investigation of the Flow Physics Driving Stall-Side Flutter in Advanced Forward Swept Fan Designs

Flutter-free operation of advanced transonic fan designs continues to be a challenging task for the designers of aircraft engines. In order to meet the demands of increased performance and lighter weight, these modern fan designs usually feature low-aspect ratio shroudless rotor blade designs that make the task of achieving adequate flutter margin even more challenging for the aeroelastician. This is especially true for advanced forward swept designs that encompass an entirely new design space compared to previous experience. Fortunately, advances in unsteady computational fluid dynamic (CFD) techniques over the past decade now provide an analysis capability that can be used to quantitatively assess the aeroelastic characteristics of these next generation fans during the design cycle. For aeroelastic applications, Mississippi State University and NASA Glenn Research Center have developed the CFD code TURBO-AE. This code is a time-accurate three-dimensional Euler/Navier-Stokes unsteady flow solver developed for axial-flow turbomachinery that can model multiple blade rows undergoing harmonic oscillations with arbitrary interblade phase angles, i.e., nodal diameter patterns. Details of the code can be found in Chen et al. (1993, 1994), Bakhle et al. (1997, 1998), and Srivastava et al. (1999). To assess aeroelastic stability, the work-per-cycle from TURBO-AE is converted to the critical damping ratio since this value is more physically meaningful, with both the unsteady normal pressure and viscous shear forces included in the work-per-cycle calculation. If the total damping (aerodynamic plus mechanical) is negative, then the blade is unstable since it extracts energy from the flow field over the vibration cycle. TURBO-AE is an integral part of an aeroelastic design system being developed at Honeywell Engines, Systems & Services for flutter and forced response predictions, with test cases from development rig and engine tests being used to validate its predictive capability. A recent experimental program (Sanders et al., 2002) was aimed at providing the necessary unsteady aerodynamic and vibratory response data needed to validate TURBO-AE for fan flutter predictions. A comparison of numerical TURBO-AE simulations with the benchmark flutter data is given in Sanders et al. (2003), with the data used to guide the validation of the code and define best practices for performing accurate unsteady simulations. The agreement between the analyses and the predictions was quite remarkable, demonstrating the ability of the analysis to accurately model the unsteady flow processes driving stall-side flutter.

Sanders, Albert J.↗

NASA NTRS2011

Sonic Boom Computations for a Mach 1.6 Cruise Low Boom Configuration and Comparisons with Wind Tunnel Data

Accurate analysis of sonic boom pressure signatures using computational fluid dynamics techniques remains quite challenging. Although CFD shows accurate predictions of flow around complex configurations, generating grids that can resolve the sonic boom signature far away from the body is a challenge. The test case chosen for this study corresponds to an experimental wind-tunnel test that was conducted to measure the sonic boom pressure signature of a low boom configuration designed by Gulfstream Aerospace Corporation. Two widely used NASA codes, USM3D and AERO, are examined for their ability to accurately capture sonic boom signature. Numerical simulations are conducted for a free-stream Mach number of 1.6, angle of attack of 0.3 and Reynolds number of 3.85x10(exp 6) based on model reference length. Flow around the low boom configuration in free air and inside the Langley Unitary plan wind tunnel are computed. Results from the numerical simulations are compared with wind tunnel data. The effects of viscous and turbulence modeling along with tunnel walls on the computed sonic boom signature are presented and discussed.

Elmiligui, Alaa A.↗

NASA NTRS1982

Design of supercritical swept wings

Computational fluid dynamics are used to discuss problems inherent to transonic three-dimensional flow past supercritical swept wings. The formulation for a boundary value problem for the flow past the wing is provided, including consideration of weak shock waves and the use of parabolic coordinates. A swept wing code is developed which requires a mesh of 152 x 10 x 12 points and 200 time cycles. A formula for wave drag is calculated, based on the idea that the conservation form of the momentum equation becomes an entropy inequality measuring the drag, expressible in terms of a small-disturbance equation for a potential function in two dimensions. The entropy inequality has been incorporated in a two-dimensional code for the analysis of transonic flow over airfoils. A method of artificial viscosity is explored for optimum pressure distributions with design, and involves a free boundary problem considering speed over only a portion of the wing.

Garabedian, P.↗

NASA NTRSDate not supplied

Aerothermal Analysis of the Dragonfly Titan Entry

The Dragonfly mission will send a rotorcraft lander to the surface of Saturn’s moon Titan as part of the New-Frontiers program. This will be the first spacecraft to land on Titan since the Huygens probe’s descent in January of 2005 and only the second spacecraft to enter Titan’s atmosphere. The Dragonfly entry capsule is significantly larger than the Huygens probe and will experience higher aerothermal environments. This poster will provide an overview of the aerothermodynamic models used in the Dragonfly aeroshell design process. Titan provides a unique entry environment that has several fundamental differences from the environments at the more common entry destinations of Earth and Mars. The lack of atmospheric oxygen significantly reduces heatshield recession (expected to be negligible), which simplifies material response analysis but also reduces the efficiency of the thermal protection system. The atmos-phere is composed of mostly nitrogen and trace amounts of methane. During hypersonic entry this methane dissociates and leads to the formation of the molecule CN, which is known to radiate strongly in the shock layer. Current Dragonfly heating predictions estimate that radiation contributes about 50% of the total heat flux along the forebody and up to 90% on the aftbody, making radiation a key element of aerothermal analysis. To properly address the importance of radiation, a unique radiative-heating correlation was developed that incorporates normal-shock equilibrium chemistry predictions based on work from the Mars 2020 program. These new heating correlations agree well with Computational Fluid Dynamics (CFD) predictions and enable large-scale trajectory analyses by providing rapid heating predictions along candidate trajectories. Dragonfly aerothermal environments are generated with the Data Parallel Line Relaxation (DPLR) finite-volume Navier-Stokes solver alongside the Nonequilibrium Radiative Transport and Spectra (NEQAIR) radiation transport code. DPLR produces shock-aligned flowfields that incorporate the effects of both chemical and thermal nonequilibrium while NEQAIR solves the radiation transport equations along lines-of-sight through the flowfield solutions to predict the radiative heating on the aeroshell. Titan’s atmosphere produces complex thermal and chemical nonequilibria and so classical radiation approximations that extrapolate from a single line-of-sight, such as the Tangent-Slab method, were found to overpredict peak radiative heat fluxes by 10-20%. Therefore, all radiative heating is directly calculated using 3-D transport for both the forebody and aftbody to reduce unnecessary conservatism. A 21-species finite-rate chemistry model, including electrons and ionized species, is used in the CFD to ensure accurate vehicle heating predictions (the level of ionization is particularly for accurate radiative flux predictions). This poster will cover: 1) The models and assumptions used to simulate entry into Titan’s atmosphere, 2) Dragonfly aerothermal engineering correlations alongside statistical heating values, and 3) Driving aerothermal considerations for a Titan entry such as CN formation and radiation-flowfield coupling.

EDL↗

NASA NTRS1986

Advanced 3-D viscous SSME turbine rotor stator CFD algorithms

Current Space Shuttle Main Engine (SSME) high pressure fuel turbopump problems have generated a desire to analyze the flow field of rotating machinery. The status of Computational Fluid Dynamics (CFD) has reached the point that soon the capability to solve unsteady three-dimensional viscous flow fields will be at hand. The work presented involves upgrading the computational efficiency of an operational three-dimensional algorithm. The modifications include algorithm development, algorithm approximation and acceleration, and special coding optimizations. The overall result of these modifications has reduced processing time by nearly 80%.

Janus, J. Mark↗

NASA NTRS1993

Liquid propellant rocket engine combustion simulation with a time-accurate CFD method

Time-accurate computational fluid dynamics (CFD) algorithms are among the basic requirements as an engineering or research tool for realistic simulations of transient combustion phenomena, such as combustion instability, transient start-up, etc., inside the rocket engine combustion chamber. A time-accurate pressure based method is employed in the FDNS code for combustion model development. This is in connection with other program development activities such as spray combustion model development and efficient finite-rate chemistry solution method implementation. In the present study, a second-order time-accurate time-marching scheme is employed. For better spatial resolutions near discontinuities (e.g., shocks, contact discontinuities), a 3rd-order accurate TVD scheme for modeling the convection terms is implemented in the FDNS code. Necessary modification to the predictor/multi-corrector solution algorithm in order to maintain time-accurate wave propagation is also investigated. Benchmark 1-D and multidimensional test cases, which include the classical shock tube wave propagation problems, resonant pipe test case, unsteady flow development of a blast tube test case, and H2/O2 rocket engine chamber combustion start-up transient simulation, etc., are investigated to validate and demonstrate the accuracy and robustness of the present numerical scheme and solution algorithm.

Chen, Y. S.↗

NASA NTRS2006

CFD Simulation of the Space Shuttle Launch Vehicle with Booster Separation Motor and Reaction Control System Plumes

The primary objective of this paper is to demonstrate the capability of computational fluid dynamics (CFD) to simulate a very complicated flow field encountered during the space shuttle ascent. The flow field features nozzle plumes from booster separation motor (BSM) and reaction control system (RCS) jets with a supersonic incoming cross flow at speed of Mach 4. The overset Navier-Stokes code OVERFLOW, was used to simulate the flow field surrounding the entire space shuttle launch vehicle (SSLV) with high geometric fidelity. The variable gamma option was chosen due to the high temperature nature of nozzle flows and different plume species. CFD predicted Mach contours are in good agreement with the schlieren photos from wind tunnel test. Flow fields are discussed in detail and the results are used to support the debris analysis for the space shuttle Return To Flight (RTF) task.

Gea, L. M.↗

NASA NTRS2006

CFD Simulation of the Space Shuttle Launch Vehicle with Booster Separation Motor and Reaction Control Plumes

Gea, L. M.↗

NASA NTRS1991

Euler solutions to nonlinear acoustics of non-lifting hovering rotor blades

For the first time a computational fluid dynamics (CFD) method is used to calculate directly the high-speed impulsive (HSI) noise of a non-lifting hovering rotor blade out to a distance of over three rotor radii. In order to accurately propagate the acoustic wave in a stable and efficient manner, an implicit upwind-biased Euler method is solved on a grid with points clustered along the line of propagation. A detailed validation of the code is performed for a rectangular rotor blade at tip Mach numbers ranging from 0.88 to 0.92. The agreement with experiment is excellent at both the sonic cylinder and at 2.18 rotor radii. The agreement at 3.09 rotor radii is still very good, showing improvements over the results from the best previous method. Grid sensitivity studies indicate that with special attention to the location of the boundaries a grid with approximately 60,000 points is adequate. This results in a computational time of approximately 40 minutes on a Cray-XMP. The practicality of the method to calculate HSI noise is demonstrated by expanding the scope of the investigation to examine the rectangular blade as well as a highly swept and tapered blade over a tip Mach number range of 0.80 to 0.95. Comparisons with experimental data are excellent and the advantages of planform modifications are clearly evident. New insight is gained into the mechanisms of nonlinear propagation and the minimum distance at which a valid comparison of different rotors can be made: approximately two rotor radii from the center of rotation.

Baeder, J. D.↗

NASA NTRS1991

Euler solutions to nonlinear acoustics of non-lifting rotor blades

Baeder, J. D.↗

NASA NTRS1992

Computational fluid dynamics combustion analysis evaluation

This study involves the development of numerical modelling in spray combustion. These modelling efforts are mainly motivated to improve the computational efficiency in the stochastic particle tracking method as well as to incorporate the physical submodels of turbulence, combustion, vaporization, and dense spray effects. The present mathematical formulation and numerical methodologies can be casted in any time-marching pressure correction methodologies (PCM) such as FDNS code and MAST code. A sequence of validation cases involving steady burning sprays and transient evaporating sprays will be included.

Kim, Y. M.↗

NASA NTRS1982

Tenth NASTRAN User's Colloquium

The development of the NASTRAN computer program, a general purpose finite element computer code for structural analysis, was discussed. The application and development of NASTRAN is presented in the following topics: improvements and enhancements; developments of pre and postprocessors; interactive review system; the use of harmonic expansions in magnetic field problems; improving a dynamic model with test data using Linwood; solution of axisymmetric fluid structure interaction problems; large displacements and stability analysis of nonlinear propeller structures; prediction of bead area contact load at the tire wheel interface; elastic plastic analysis of an overloaded breech ring; finite element solution of torsion and other 2-D Poisson equations; new capability for elastic aircraft airloads; usage of substructuring analysis in the get away special program; solving symmetric structures with nonsymmetric loads; evaluation and reduction of errors induced by Guyan transformation.

Source record↗

NASA NTRS1991

Computations of unsteady multistage compressor flows in a workstation environment

High-end graphics workstations are becoming a necessary tool in the Computational Fluid Dynamics (CFD) environment. In addition to their graphics capabilities, the latest generation of workstations have powerful floating point operation capabilities. As workstations become common, they could provide valuable computing time for applications, such as turbomachinery flow calculations. This paper discusses the issues involved in implementing an unsteady, viscous multistage turbomachinery code (STAGE-2) on workstations. The workstation version of STAGE-2 has then been used to study the effects of axial-gap spacing on the time-averaged and unsteady flow within a 2 1/2-stage compressor. Results include force polar plots, time-averaged pressure contours, standard deviation of pressure contours, time-averaged surface pressures and pressure amplitudes.

Gundy-Burlet, Karen L.↗

NASA NTRS2003

Computational Study of Axisymmetric Off-Design Nozzle Flows

Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles operating off-design at transonic Mach numbers have been completed. These computations span the very difficult transonic flight regime with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined, including the Explicit Algebraic Stress model. Computations have been completed at freestream Mach numbers of 0.9 and 1.2, and nozzle pressure ratios (NPR) of 4 and 6. Calculations completed with variable time-stepping (steady-state) did not converge to a true steady-state solution. Calculations obtained using constant timestepping (timeaccurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was the result of using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show reasonable agreement with experimental data. The SST turbulence model demonstrates the best overall agreement with experimental data.

DalBello, Teryn↗

NASA NTRS2014

Long Hole Film Cooling Dataset for CFD Development - Flow and Film Effectiveness

An experiment investigating flow and heat transfer of long (length to diameter ratio of 18) cylindrical film cooling holes has been completed. In this paper, the thermal field in the flow and on the surface of the film cooled flat plate is presented for nominal freestream turbulence intensities of 1.5 and 8 percent. The holes are inclined at 30 deg above the downstream direction, injecting chilled air of density ratio 1.0 onto the surface of a flat plate. The diameter of the hole is 0.75 in. (approx. 0.02 m) with center to center spacing (pitch) of 3 hole diameters. Coolant was injected into the mainstream flow at nominal blowing ratios of 0.5, 1.0, 1.5, and 2.0. The Reynolds number of the freestream was approximately 11,000 based on hole diameter. Thermocouple surveys were used to characterize the thermal field. Infrared thermography was used to determine the adiabatic film effectiveness on the plate. Hotwire anemometry was used to provide flowfield physics and turbulence measurements. The results are compared to existing data in the literature. The aim of this work is to produce a benchmark dataset for Computational Fluid Dynamics (CFD) development to eliminate the effects of hole length to diameter ratio and to improve resolution in the near-hole region. In this report, a Time Filtered Navier Stokes (TFNS), also known as Partially Resolved Navier Stokes (PRNS), method that was implemented in the Glenn-HT code is used to model coolant-mainstream interaction. This method is a high fidelity unsteady method that aims to represent large scale flow features and mixing more accurately.

Shyam, Vikram↗

NASA NTRSDate not supplied

Propagation Analysis of the 3rd Sonic Boom Prediction Workshop Cases using sBOOM

The 3rd Sonic Boom Prediction Workshop was held on January 4-5, 2020 in conjunction with the AIAA SciTech conference. The workshop had 23 participants who represented 6 different countries and 18 unique organizations that included government, industry, and academia. The workshop was attended by over 50 individuals. The motivation for the workshop stems from the goal of obtaining supersonic commercial overland flight. In order to replace the current prohibition with a certification standard, an international effort is required to quantify the accuracy and reliability of prediction methods. The workshops also identify deficiencies in existing methods where further research should be focused. The workshop consisted of two days with the first day focused on near-field computational fluid dynamics (CFD) and the second day on propagation techniques. For the atmospheric propagation portion of the workshop, two required cases and one optional case were prescribed for the participants to exercise their propagation implementations. This paper presents and discusses the results that were obtained using the sonic boom propagation code sBOOM. Results show that sonic boom carpets could be much different under arbitrary atmospheric conditions as opposed to standard atmospheric conditions. Additionally, several unique features when using sBOOM are presented.

Sonic Boom↗

NASA NTRSDate not supplied

Towards a Third-Order Accurate, Second-Derivative-Free, Shock-Capturing Finite-Volume Method for Hypersonic Flows on Tetrahedral Grids

In this paper, we report progress in the development of a third-order accurate, second-derivative-free, shock-capturing finite-volume solver for three-dimensional unstructured grids. The method is economical in the sense that the computation and storage of second derivatives are not required for third-order accuracy. It is based on point-valued numerical solutions stored at cells, gradients computed and stored at nodes, and an efficient projected-derivative formula that eliminates the need for second derivatives in a quadratic solution interpolation. The projected-derivative formula is also used to eliminate second derivatives from a high-order flux quadrature formula, so that it can be implemented conveniently in the form of a numerical flux at a face center plus a correction term. Similarly, a high-order source quadrature formula can also be implemented in the form of a cell-center point evaluation plus a similar correction term. These features make it relatively straightforward to extend an existing second-order finite-volume code to third-order. This paper reports progress of implementing the method in the NASA VULCAN-CFD code and discusses the implementation of a high-order accurate limiter for shock capturing.

Weighted Least-Squares↗

NASA NTRSDate not supplied

Energy of Fluid (EOF) Simulation of Large-Scale Densification and Solidification of Hydrogen

The current research develops a finite volume-based Computation Fluid Dynamic (CFD) model utilizing an Energy of Fluid (EOF) approach to generate a simulation of the densification of liquid hydrogen (LH2) in the Integrated Refrigeration and Storage (IRAS) experimental tank for the Ground Operations Demonstration Unit for Liquid Hydrogen Project (GODU-LH2) at the Kennedy Space Center (KSC). The computational code will integrate a commercial software pressure-based model with User Defined Functions (UDF) generated and added to the model. The modified model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. This method will decrease computational time when compared to the built-in commercially available Enthalpy solvers. Currently the research focused on utilizing User Defined Scalars (UDS) compared to the enthalpy model and a lumped node analysis to show the written UDS is as accurate as the built-in solvers, and the models will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

C. Mahony↗