Search NASA⌕ Search

SEARCH · Search NASA

Results for “Navier Stokes 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 163 records · Page 9

Growth Characteristics Downstream of a Shallow Bump: Computation and Experiment

Measurements of the velocity field created by a shallow bump on a wall revealed that an energy peak in the spanwise spectrum associated with the driver decays and an initially small-amplitude secondary mode rapidly grows with distance downstream of the bump. Linear theories could not provide an explanation for this growing mode. The present Navier-Stokes simulation replicates and confirms the experimental results. Insight into the structure of the flow was obtained from a study of the results of the calculations and is presented.

Joslin, Ronald D.↗

Navier-Stokes Computations on Full Wing-Body Configuration with Oscillating Control Surfaces

Unsteady Navier-Stokes simulations have been performed for vortical flows over an "arrow-wing" configuration of a supersonic transport in the transonic regime. Computed steady pressures and integrated force coefficients with and without control surface deflection at a moderate angle of attack are compared with experiment. For unsteady cases, oscillating trailing-edge control surfaces are modeled by using moving grids. Response characteristics between symmetric and antisymmetric oscillatory motions of the control surfaces on the left and right wings are studied. The antisymmetric case produces higher lift than the steady case with no deflection and the unsteady symmetric case produces higher lift than the antisymmetric case. The detailed analysis of the wake structure revealed a strong interaction between the primary vortex and the wake vortex sheet from the flap region when the flap is deflected up.

Obayashi, Shigeru↗

Study of a Wake Recovery Mechanism in a High-Speed Axial Compressor Stage

This work addresses the significant differences in compressor rotor wake mixing loss which exist in a stage environment relative to a rotor in isolation. The wake decay for a rotor in isolation is due solely to viscous dissipation which is an irreversible process and thus leads to a loss in both total pressure and efficiency. Rotor wake decay in the stage environment is due to both viscous mixing and the inviscid strain imposed on the wake fluid particles by the stator velocity field. This straining process, referred to by Smith (1993) as recovery, is reversible and for a 2D rotor wake leads to an inviscid reduction of the velocity deficit of the wake. A model for the rotor wake decay process is developed and used to quantify the viscous dissipation effects relative to those of inviscid wake stretching. The model is verified using laser anemometer measurements acquired in the wake of a transonic rotor operated in isolation and in a stage configuration at near peak efficiency and near stall operating conditions. Additional insight is provided by a time-accurate 3D Navier-Stokes simulation of the compressor stator flow field at the corresponding stage loading levels. Results from the wake decay model exhibit good agreement with the experimental data. Data from the model, laser anemometer measurements, and numerical simulations indicate that for the rotor/stator spacing used in this work, which is typical of core compressors, rotor wake straining (stretching) is the primary decay process in the stator passage with viscous mixing playing only a minor role. The implications of these results on compressor stage design are discussed.

VanZante, Dale E.↗

Inlet Development for a Rocket Based Combined Cycle, Single Stage to Orbit Vehicle Using Computational Fluid Dynamics

Design and analysis of the inlet for a rocket based combined cycle engine is discussed. Computational fluid dynamics was used in both the design and subsequent analysis. Reynolds averaged Navier-Stokes simulations were performed using both perfect gas and real gas assumptions. An inlet design that operates over the required Mach number range from 0 to 12 was produced. Performance data for cycle analysis was post processed using a stream thrust averaging technique. A detailed performance database for cycle analysis is presented. The effect ot vehicle forebody compression on air capture is also examined.

DeBonis, J. R.↗

Effects of Hot Streak Shape on Rotor Heating in a High-Subsonic Single-Stage Turbine

Experimental data have shown that combustor temperature non-uniformities can lead to the excessive heating of first-stage rotor blades in turbines. This heating of the rotor blades can lead to thermal fatigue and degrade turbine performance. The results of recent studies have shown that variations in the circumferential location (clocking) of the hot streak relative to the first-stage vane airfoils can be used to minimize the adverse effects of the hot streak. The effects of the hot streak/airfoil count ratio on the heating patterns of turbine airfoils have also been evaluated. In the present investigation, three-dimensional unsteady Navier-Stokes simulations have been performed for a single-stage high-pressure turbine operating in high subsonic flow. In addition to a simulation of the baseline turbine, simulations have been performed for circular and elliptical hot streaks of varying sizes in an effort to represent different combustor designs. The predicted results for the baseline simulation show good agreement with the available experimental data. The results of the hot streak simulations indicate: that a) elliptical hot streaks mix more rapidly than circular hot streaks, b) for small hot streak surface area the average rotor temperature is not a strong function of hot streak temperature ratio or shape, and c) hot streaks with larger surface area interact with the secondary flows at the rotor hub endwall, generating an additional high temperature region.

Dorney, Daniel J.↗

Computational Analysis for Rocket-Based Combined-Cycle Systems During Rocket-Only Operation

A series of Reynolds-averaged Navier-Stokes calculations were employed to study the performance of rocket-based combined-cycle systems operating in an all-rocket mode. This parametric series of calculations were executed within a statistical framework, commonly known as design of experiments. The parametric design space included four geometric and two flowfield variables set at three levels each, for a total of 729 possible combinations. A D-optimal design strategy was selected. It required that only 36 separate computational fluid dynamics (CFD) solutions be performed to develop a full response surface model, which quantified the linear, bilinear, and curvilinear effects of the six experimental variables. The axisymmetric, Reynolds-averaged Navier-Stokes simulations were executed with the NPARC v3.0 code. The response used in the statistical analysis was created from Isp efficiency data integrated from the 36 CFD simulations. The influence of turbulence modeling was analyzed by using both one- and two-equation models. Careful attention was also given to quantify the influence of mesh dependence, iterative convergence, and artificial viscosity upon the resulting statistical model. Thirteen statistically significant effects were observed to have an influence on rocket-based combined-cycle nozzle performance. It was apparent that the free-expansion process, directly downstream of the rocket nozzle, can influence the Isp efficiency. Numerical schlieren images and particle traces have been used to further understand the physical phenomena behind several of the statistically significant results.

Steffen, C. J., Jr.↗

Geometric Perturbations and Asymmetric Vortex Shedding about Slender Pointed Bodies

The flow about slender, pointed bodies can be characterized by different states with angle of attack. At moderate-to-high angles of attack (alpha approx. = 40 deg.), a steady, asymmetric vortex pattern develops along the body, leading to a net lateral force. At higher angles of attack (alpha approx. = 60 deg.), the aft-end of the body develops an unsteady von Karman shedding. As the angle of attack approaches 90 deg., the entire body length exhibits a time-dependent vortex shedding pattern. The current work uses three-dimensional, thin-layer Navier-Stokes simulations to investigate the physical mechanisms of asymmetric vortex shedding at alpha = 40 deg. and alpha = 60 deg. The development of an asymmetric vortex pattern via a convective instability mechanism is investigated using tip bumps, surface roughness, and tip curvature. It's found that surface roughness simulations can incite an asymmetric vortex state at alpha = 60 deg. which is consistent with the application of a tip bump, and the experimentally observed flowfield. The unsteady von Karman vortex shedding on the aft. portion of the body is also well resolved. The use of surface roughness did not incite a flow asymmetry at alpha = 40 deg., and it was necessary to simulate tip curvature at this angle of attack in order to generate an asymmetric vortex state.

Murman, Scott M.↗

Geometric Perturbations and Asymmetric vortex shedding About Slender Pointed Bodies

The flow about slender, pointed bodies can be characterized by different states with angle of attack. At moderate-to-high angles of attack (alpha approximately equals 40deg), a steady, asymmetric vortex pattern develops along the body, leading to a net lateral force. At higher angles of attack (alpha approximately equals 60deg), the aft-end of the body develops an unsteady von Karman shedding. As the angle of attack approaches 90deg, the entire body length exhibits a time-dependent vortex shedding pattern. The current work uses three-dimensional, thin-layer Navier-Stokes simulations to investigate the physical mechanisms of asymmetric vortex shedding at alpha = 40deg and alpha = 60deg. The development of an asymmetric vortex pattern via a convective instability mechanism is investigated using tip bumps, surface roughness, and tip curvature. It's found that surface roughness simulations can incite an asymmetric vortex state at alpha = 60deg which is consistent with the application of a tip bumps, and the experimentally observed flowfield. The unsteady von Karman vortex shedding on the aft portion of the body is also well resolved. The use of surface roughness did not incite a flow asymmetry at alpha = 40deg, and it was necessary to simulate tip curvature at this angle of attack in order to generate an asymmetric vortex state.

Murman, Scott M.↗

Unsteady Flow in a Supersonic Turbine with Variable Specific Heats

Modern high-work turbines can be compact, transonic, supersonic, counter-rotating, or use a dense drive gas. The vast majority of modern rocket turbine designs fall into these Categories. These turbines usually have large temperature variations across a given stage, and are characterized by large amounts of flow unsteadiness. The flow unsteadiness can have a major impact on the turbine performance and durability. For example, the Space Transportation Main Engine (STME) fuel turbine, a high work, transonic design, was found to have an unsteady inter-row shock which reduced efficiency by 2 points and increased dynamic loading by 24 percent. The Revolutionary Reusable Technology Turbopump (RRTT), which uses full flow oxygen for its drive gas, was found to shed vortices with such energy as to raise serious blade durability concerns. In both cases, the sources of the problems were uncovered (before turbopump testing) with the application of validated, unsteady computational fluid dynamics (CFD) to the designs. In the case of the RRTT and the Alternate Turbopump Development (ATD) turbines, the unsteady CFD codes have been used not just to identify problems, but to guide designs which mitigate problems due to unsteadiness. Using unsteady flow analyses as a part of the design process has led to turbine designs with higher performance (which affects temperature and mass flow rate) and fewer dynamics problems. One of the many assumptions made during the design and analysis of supersonic turbine stages is that the values of the specific heats are constant. In some analyses the value is based on an average of the expected upstream and downstream temperatures. In stages where the temperature can vary by 300 to 500 K, however, the assumption of constant fluid properties may lead to erroneous performance and durability predictions. In this study the suitability of assuming constant specific heats has been investigated by performing three-dimensional unsteady Navier-Stokes simulations for a supersonic turbine stage.

Dorney, Daniel J.↗

Design Evolution and Performance Characterization of the GTX Air-Breathing Launch Vehicle Inlet

The design and analysis of a second version of the inlet for the GTX rocket-based combine-cycle launch vehicle is discussed. The previous design did not achieve its predicted performance levels due to excessive turning of low-momentum comer flows and local over-contraction due to asymmetric end-walls. This design attempts to remove these problems by reducing the spike half-angle to 10- from 12-degrees and by implementing true plane of symmetry end-walls. Axisymmetric Reynolds-Averaged Navier-Stokes simulations using both perfect gas and real gas, finite rate chemistry, assumptions were performed to aid in the design process and to create a comprehensive database of inlet performance. The inlet design, which operates over the entire air-breathing Mach number range from 0 to 12, and the performance database are presented. The performance database, for use in cycle analysis, includes predictions of mass capture, pressure recovery, throat Mach number, drag force, and heat load, for the entire Mach range. Results of the computations are compared with experimental data to validate the performance database.

DeBonis, J. R.↗

Effects of Endwall Geometry and Stacking on Two-Stage Supersonic Turbine Performance

The drive towards high-work turbines has led to designs which can be compact, transonic, supersonic, counter rotating, or use a dense drive gas. These aggressive designs can lead to strong secondary flows and airfoil flow separation. In many cases the secondary and separated flows can be minimized by contouring the hub/shroud endwalls and/or modifying the airfoil stacking. In this study, three-dimensional unsteady Navier-Stokes simulations were performed to study three different endwall shapes between the first-stage vanes and rotors, as well as two different stackings for the first-stage vanes. The predicted results indicate that changing the stacking of the first-stage vanes can significantly impact endwall separation (and turbine performance) in regions where the endwall profile changes.

Dorney, Daniel J.↗

Effects of Impeller-Diffuser Interaction on Centrifugal Compressor Performance

This research program focuses on characterizing the effect of impeller-diffuser interactions in a centrifugal compressor stage on its performance using unsteady threedimensional Reynolds-averaged Navier-Stokes simulations. The computed results show that the interaction between the downstream diffuser pressure field and the impeller tip clearance flow can account for performance changes in the impeller. The magnitude of performance change due to this interaction was examined for an impeller with varying tip clearance followed by a vaned or vaneless diffuser. The impact of unsteady impeller-diffuser interaction, primarily through the impeller tip clearance flow, is reflected through a time-averaged change in impeller loss, blockage and slip. The results show that there exists a tip clearance where the beneficial effect of the impeller-diffuser interaction on the impeller performance is at a maximum. A flow feature that consists of tip flow back leakage was shown to occur at design speed for the centrifugal compressor stage. This flow phenomenon is described as tip flow that originates in one passage, flows downstream of the impeller trailing edge and then returns to upstream of the impeller trailing edge of a neighboring passage. Such a flow feature is a source of loss in the impeller. A hypothesis is put forth to show that changing the diffuser vane count and changing impeller-diffuser gap has an analogous effect on the impeller performance. The centrifugal compressor stage was analyzed using diffusers of different vane counts, producing an impeller performance trend similar to that when the impeller-diffuser gap was varied, thus supporting the hypothesis made. This has the implication that the effect impeller performance associated with changing the impeller-diffuser gap and changing diffuser vane count can be described by the non-dimensional ratio of impeller-diffuser gap to diffuser vane pitch. A procedure is proposed and developed for isolating impeller passage blockage change without the need to define the region of blockage generation (which may incur a certain degree of arbitrariness). This method has been assessed for its applicability and utility.

Tan, Choon S.↗

Rotor-Stator Interaction Performance Effects

Decreased axial spacing between blade rows in an axial compressor stage is thought to increase stage performance because of an unsteady process that occurs in the downstream blade row and acts on the upstream blade row wakes. This process results in the "recovery" of part of the wake energy before all of this energy is irreversibly lost due to viscous diffusion. To study the wake-blade interaction mechanism, researchers at the NASA Lewis Research Center acquired two-component Laser Fringe Anemometer measurements of the rotor wake in the single-stage transonic compressor at two stage loading levels. The detailed measurements were acquired for one stator pitch in circumference at axial positions from the rotor trailing edge to 20 percent of the stator axial chord, at the exit of the stator passage, and downstream of the stator row including the stator wake. These data show that the changes in wake energy that occur inside the stator passage are not due to viscous dissipation alone, and thus the data provide evidence that "wake recovery" is occurring. A time-accurate, three-dimensional Navier Stokes simulation of the compressor stator was done at the corresponding stage loading levels. The measurements and simulations are being used in combination to show the effects of stator blade loading, quantify the effects of viscosity, and quantify the stage efficiency gain due to the wake recovery process. The accuracy of simple models of the wake recovery process is being evaluated in an effort to include the effects of wake recovery in the NASA-developed Average Passage code for multistage turbomachinery simulations.

VanZante, Dale E.↗

Three-Dimensional Navier-Stokes Calculations Using the Modified Space-Time CESE Method

The space-time conservation element solution element (CESE) method is modified to address the robustness issues of high-aspect-ratio, viscous, near-wall meshes. In this new approach, the dependent variable gradients are evaluated using element edges and the corresponding neighboring solution elements while keeping the original flux integration procedure intact. As such, the excellent flux conservation property is retained and the new edge-based gradients evaluation significantly improves the robustness for high-aspect ratio meshes frequently encountered in three-dimensional, Navier-Stokes calculations. The order of accuracy of the proposed method is demonstrated for oblique acoustic wave propagation, shock-wave interaction, and hypersonic flows over a blunt body. The confirmed second-order convergence along with the enhanced robustness in handling hypersonic blunt body flow calculations makes the proposed approach a very competitive CFD framework for 3D Navier-Stokes simulations.

Chang, Chau-lyan↗

Gaseous Surrogate Hydrocarbons for a Hifire Scramjet that Mimic Opposed Jet Extinction Limits for Cracked JP Fuels

This paper describes, first, the top-down methodology used to define simple gaseous surrogate hydrocarbon (HC) fuel mixtures for a hypersonic scramjet combustion subtask of the HiFIRE program. It then presents new and updated Opposed Jet Burner (OJB) extinction-limit Flame Strength (FS) data obtained from laminar non-premixed HC vs. air counterflow diffusion flames at 1-atm, which follow from earlier investigations. FS represents a strain-induced extinction limit based on cross-section-average air jet velocity, U(sub air), that sustains combustion of a counter jet of gaseous fuel just before extinction. FS uniquely characterizes a kinetically limited fuel combustion rate. More generally, Applied Stress Rates (ASRs) at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or t) can directly be compared with extinction limits determined numerically using either a 1-D or (preferably) a 2-D Navier Stokes simulation with detailed transport and finite rate chemistry. The FS results help to characterize and define three candidate surrogate HC fuel mixtures that exhibit a common FS 70% greater than for vaporized JP-7 fuel. These include a binary fuel mixture of 64% ethylene + 36% methane, which is our primary recommendation. It is intended to mimic the critical flameholding limit of a thermally- or catalytically-cracked JP-7 like fuel in HiFIRE scramjet combustion tests. Our supporting experimental results include: (1) An idealized kinetically-limited ASR reactivity scale, which represents maximum strength non-premixed flames for several gaseous and vaporized liquid HCs; (2) FS characterizations of Colket and Spadaccini s suggested ternary surrogate, of 60% ethylene + 30% methane + 10% n-heptane, which matches the ignition delay of a typical cracked JP fuel; (3) Data showing how our recommended binary surrogate, of 64% ethylene + 36% methane, has an identical FS; (4) Data that characterize an alternate surrogate of 44% ethylene + 56% ethane with identical FS and nearly equal molecular weights; this could be useful when systematically varying the fuel composition. However, the mixture liquefies at much lower pressure, which limits on-board storage of gaseous fuel; (5) Dynamic Flame Weakening results that show how oscillations in OJB input flow (and composition) can weaken (extinguish) surrogate flames up to 200 Hz, but the weakening is 2.5x smaller compared to pure methane; and finally, (6) FS limits at 1-atm that compare with three published 1-D numerical OJB extinction results using four chemical kinetic models. The methane kinetics generally agree closely at 1-atm, whereas, the various ethylene models predict extinction limits that average ~ 45% high, which represents a significant problem for numerical simulation of surrogate-based flameholding in a scramjet cavity. Finally, we continue advocating the FS approach as more direct and fundamental for assessing idealized scramjet flameholding potentials than measurements of "unstrained" premixed laminar burning velocity or blowout in a Perfectly Stirred Reactor.

Pellett, Gerald L.↗

Computational Analysis of Arc-Jet Wedge Tests Including Ablation and Shape Change

Coupled fluid-material response analyses of arc-jet wedge ablation tests conducted in a NASA Ames arc-jet facility are considered. These tests were conducted using blunt wedge models placed in a free jet downstream of the 6-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. The fluid analysis includes computational Navier-Stokes simulations of the nonequilibrium flowfield in the facility nozzle and test box as well as the flowfield over the models. The material response analysis includes simulation of two-dimensional surface ablation and internal heat conduction, thermal decomposition, and pyrolysis gas flow. For ablating test articles undergoing shape change, the material response and fluid analyses are coupled in order to calculate the time dependent surface heating and pressure distributions that result from shape change. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator. Effects of the test article shape change on fluid and material response simulations are demonstrated, and computational predictions of surface recession, shape change, and in-depth temperatures are compared with the experimental measurements.

Goekcen, Tahir↗

Rotary-Wing Relevant Compressor Aero Research and Technology Development Activities at Glenn Research Center

Technical challenges of compressors for future rotorcraft engines are driven by engine-level and component-level requirements. Cycle analyses are used to highlight the engine-level challenges for 3000, 7500, and 12000 SHP-class engines, which include retention of performance and stability margin at low corrected flows, and matching compressor type, axial-flow or centrifugal, to the low corrected flows and high temperatures in the aft stages. At the component level: power-to-weight and efficiency requirements impel designs with lower inherent aerodynamic stability margin; and, optimum engine overall pressure ratios lead to small blade heights and the associated challenges of scale, particularly increased clearance-to-span ratios. The technical challenges associated with the aerodynamics of low corrected flows and stability management impel the compressor aero research and development efforts reviewed herein. These activities include development of simple models for clearance sensitivities to improve cycle calculations, full-annulus, unsteady Navier-Stokes simulations used to elucidate stall, its inception, and the physics of stall control by discrete tip-injection, development of an actuator-duct-based model for rapid simulation of nonaxisymmetric flow fields (e.g., due inlet circumferential distortion), advanced centrifugal compressor stage development and experimentation, and application of stall control in a T700 engine.

Welch, Gerard E.↗

Visualization and Quantification of Rotor Tip Vortices in Helicopter Flows

This paper presents an automated approach for effective extraction, visualization, and quantification of vortex core radii from the Navier-Stokes simulations of a UH-60A rotor in forward flight. We adopt a scaled Q-criterion to determine vortex regions and then perform vortex core profiling in these regions to calculate vortex core radii. This method provides an efficient way of visualizing and quantifying the blade tip vortices. Moreover, the vortices radii are displayed graphically in a plane.

Quantification↗