Results for “flow simulation” — Page 75

NASA NTRS1993

Measurements and computational analysis of heat transfer and flow in a simulated turbine blade internal cooling passage

Visual and quantitative information was obtained on heat transfer and flow in a branched-duct test section that had several significant features of an internal cooling passage of a turbine blade. The objective of this study was to generate a set of experimental data that could be used to validate computer codes for internal cooling systems. Surface heat transfer coefficients and entrance flow conditions were measured at entrance Reynolds numbers of 45,000, 335,000, and 726,000. The heat transfer data were obtained using an Inconel heater sheet attached to the surface and coated with liquid crystals. Visual and quantitative flow field results using particle image velocimetry were also obtained for a plane at mid channel height for a Reynolds number of 45,000. The flow was seeded with polystyrene particles and illuminated by a laser light sheet. Computational results were determined for the same configurations and at matching Reynolds numbers; these surface heat transfer coefficients and flow velocities were computed with a commercially available code. The experimental and computational results were compared. Although some general trends did agree, there were inconsistencies in the temperature patterns as well as in the numerical results. These inconsistencies strongly suggest the need for further computational studies on complicated geometries such as the one studied.

Russell, Louis M.↗

NASA NTRS1993

Measurements and computational analysis of heat transfer and flow in a simulated turbine blade internal cooling passage

Visual and quantitative information was obtained on heat transfer and flow in a branched-duct test section that had several significant features of an internal cooling passage of a turbine blade. The objective of this study was to generate a set of experimental data that could be used to validate computer codes for internal cooling systems. Surface heat transfer coefficients and entrance flow conditions were measured at entrance Reynolds numbers of 45,000, 335,000, and 726,000. The heat transfer data were obtained using an Inconel heater sheet attached to the surface and coated with liquid crystals. Visual and quantitative flow field results using particle image velocimetry were also obtained for a plane at mid channel height for a Reynolds number of 45,000. The flow was seeded with polystyrene particles and illuminated by a laser light sheet. Computational results were determined for the same configurations and at matching Reynolds numbers; these surface heat transfer coefficients and flow velocities were computed with a commercially available code. The experimental and computational results were compared. Although some general trends did agree, there were inconsistencies in the temperature patterns as well as in the numerical results. These inconsistencies strongly suggest the need for further computational studies on complicated geometries such as the one studied.

Russell, Louis M.↗

NASA NTRS1989

Numerical simulation of unsteady flows generated by dissociating nitrogen diffusion

An explicit time accurate predictor-corrector scheme, similar to MacCormack's, has been used to simulate unsteady nonequilibrium laminar chemically-reacting diffusion-reaction generated flows. A zonal approach, along with nonstructured boundaries, allows computation of internal and external flow with the same computer code. Results for equilibrium flows generated by dissociating, vibrationally relaxing nitrogen in a rectangular chamber are presented. The unsteady flows were generated strictly by the chemical reactions and concentration gradients.

Marraffa, Lionel↗

NASA NTRS1992

Numerical simulation of transient hypervelocity flow in an expansion tube

Several numerical simulations of the transient flow of helium in an expansion tube are presented in an effort to identify some of the basic mechanisms which cause the noisy test flows seen in experiments. The calculations were performed with an axisymmetric Navier-Stokes code based on a finite volume formulation and upwinding techniques. Although laminar flow and ideal bursting of the diaphragms was assumed, the simulations showed some of the important features seen in experiments. In particular, the discontinuity in tube diameter of the primary diaphragm station introduced a transverse perturbation to the expanding driver gas and this perturbation was seen to propagate into the test gas under some flow conditions. The disturbances seen in the test flow can be characterized as either small amplitude, low frequency noise possibly introduced during shock compression or large amplitude, high frequency noise associated with the passage of the reflected head of the unsteady expansion.

Jacobs, P. A.↗

NASA NTRS1992

Numerical simulation of transient hypervelocity flow in an expansion tube

Several numerical simulations of the transient flow of helium in an expansion tube are presented. The aim of the exercise is to provide further information on the operational problems of the NASA Langley expansion tube. The calculations were performed with an axisymmetric Navier-Stokes code based on a finite-volume formulation and upwinding techniques. Although laminar flow and ideal bursting of the diaphragms was assumed, the simulations showed some of the important features seen in the experiments. In particular, the discontinuity in the tube diameter at the primary diaphragm station introduced a transverse perturbation to the expanding driver gas, and this perturbation was seen to propagate into the test gas under some flow conditions. The disturbances seen in the test flow can be characterized as either 'small-amplitude' noise possibly introduced during shock compression or 'large-amplitude' noise associated with the passage of the reflected head of the unsteady expansion.

Jacobs, P. A.↗

NASA NTRS2019

Effect of Spatial Filtering in Implicit Large-Eddy Simulations of Separated Flows

The relatively high Reynolds number of turbulent flows encountered in various applications puts these problems well beyond the reach of direct numerical simulation (DNS) at present. Meanwhile, lower-fidelity Reynolds-averaged Navier-Stokes (RANS) calculations are known to be not accurate enough in complex problems, such as smooth-body flow separation and other flows involving highly-unsteady phenomena. Hence, given the current infeasibility of DNS and the unsatisfactory performance of RANS, intermediate techniques such as large-eddy simulation (LES) and hybrid RANS-LES, whose fidelity lie between RANS and DNS, have received much attention for application to various problems of practical importance. Modeling of the effect of missing scales on resolved scales, also known as subgrid-scale (SGS) modeling, is an important subject for LES. SGS models can be broadly categorized as explicit or implicit approaches. The explicit approach is based on an SGS model that explicitly appears in the governing equations expressed in the form of so-called “filtered Navier-Stokes equations”, which describe the evolution of the turbulence scales resolved by the LES grid. The effect of the scales unresolved by the grid is represented by the SGS model. The implicit modeling approach, on the other hand, does not employ an explicit model but instead treats the intrinsic dissipation of the numerical discretization scheme as an implicit SGS model. An LES without an explicit SGS model is commonly termed as an implicit LES (ILES). The relative merits of one SGS modeling approach over another is a subject of ongoing debate. We have opted to employ an ILES methodology, based on high-order compact finite-difference and spatial filtering schemes, in our recent investigations of separated flow problems [1, 2]. Further discussion of our preference of ILES over explicit LES is provided in Uzun and Malik [2]. The spatial filtering operation, described in the next section, is treated as an implicit SGS model for the ILES. Some observations made during the course of our recent investigations, which pointed out to excessive numerical dissipation in certain parts of the flowfield, prompted us to take a closer look at the potential effect of the spatial filter on ILES predictions. This technical note is therefore devoted to spatial filter effects in the context of a high Reynolds number, transonic shock-induced separated flow.

Ali Uzun↗

NASA NTRS1993

The nonlinear Galerkin method: A multi-scale method applied to the simulation of homogeneous turbulent flows

Using results of Direct Numerical Simulation (DNS) in the case of two-dimensional homogeneous isotropic flows, the behavior of the small and large scales of Kolmogorov like flows at moderate Reynolds numbers are first analyzed in detail. Several estimates on the time variations of the small eddies and the nonlinear interaction terms were derived; those terms play the role of the Reynolds stress tensor in the case of LES. Since the time step of a numerical scheme is determined as a function of the energy-containing eddies of the flow, the variations of the small scales and of the nonlinear interaction terms over one iteration can become negligible by comparison with the accuracy of the computation. Based on this remark, a multilevel scheme which treats differently the small and the large eddies was proposed. Using mathematical developments, estimates of all the parameters involved in the algorithm, which then becomes a completely self-adaptive procedure were derived. Finally, realistic simulations of (Kolmorov like) flows over several eddy-turnover times were performed. The results are analyzed in detail and a parametric study of the nonlinear Galerkin method is performed.

Debussche, A.↗

NASA NTRS1989

On the structure of pressure fluctuations in simulated turbulent channel flow

Pressure fluctuations in a turbulent channel flow are investigated by analyzing a database obtained from a direct numerical simulation. Detailed statistics associated with the pressure fluctuations are presented. Characteristics associated with the rapid (linear) and slow (nonlinear) pressure are discussed. It is found that the slow pressure fluctuations are larger than the rapid pressure fluctuations throughout the channel except very near the wall, where they are about the same magnitude. This is contrary to the common belief that the nonlinear source terms are negligible compared to the linear source terms. Probability density distributions, power spectra, and two-point correlations are examined to reveal the characteristics of the pressure fluctuations. The global dependence of the pressure fluctuations and pressure-strain correlations are also examined by evaluating the integral associated with Green's function representations of them. In the wall region where the pressure-strain terms are large, most contributions to the pressure-strain terms are from the wall region (i.e., local), whereas away from the wall where the pressure-strain terms are small, contributions are global. Structures of instantaneous pressure and pressure gradients at the wall and the corresponding vorticity field are examined.

Kim, John↗

NASA NTRS1989

On the structure of pressure fluctuations in simulated turbulent channel flow

Pressure fluctuations in a turbulent channel flow are investigated by analyzing a database obtained from a direct numerical simulation. Detailed statistics associated with the pressure fluctuations are presented. Characteristics associated with the rapid (linear) and slow (nonlinear) pressure are discussed. It is found that the slow pressure fluctuations are larger than the rapid pressure fluctuations throughout the channel except very near the wall, where they are about the same magnitude. This is contrary to the common belief that the nonlinear source terms are negligible compared to the linear source terms. Probability density distributions, power spectra, and two-point correlations are examined to reveal the characteristics of the pressure fluctuations. The global dependence of the pressure fluctuations and pressure-strain correlations are also examined by evaluating the integral associated with Green's function representations of them. In the wall region where the pressure-strain terms are large, most contributions to the pressure-strain terms are from the wall region (i.e., local), whereas away from the wall where the pressure-strain terms are small, contributions are global. Structures of instantaneous pressure and pressure gradients at the wall and the corresponding vorticity field are examined.

Kim, John↗

DOE OSTI2024

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗

NASA NTRS2004

Detached-Eddy Simulations of Separated Flow Around Wings With Ice Accretions: Year One Report

A computational investigation was performed to assess the effectiveness of Detached-Eddy Simulation (DES) as a tool for predicting icing effects. The AVUS code, a public domain flow solver, was employed to compute solutions for an iced wing configuration using DES and steady Reynolds Averaged Navier-Stokes (RANS) equation methodologies. The configuration was an extruded GLC305/944-ice shape section with a rectangular planform. The model was mounted between two walls so no tip effects were considered. The numerical results were validated by comparison with experimental data for the same configuration. The time-averaged DES computations showed some improvement in lift and drag results near stall when compared to steady RANS results. However, comparisons of the flow field details did not show the level of agreement suggested by the integrated quantities. Based on our results, we believe that DES may prove useful in a limited sense to provide analysis of iced wing configurations when there is significant flow separation, e.g., near stall, where steady RANS computations are demonstrably ineffective. However, more validation is needed to determine what role DES can play as part of an overall icing effects prediction strategy. We conclude the report with an assessment of existing computational tools for application to the iced wing problem and a discussion of issues that merit further study.

Choo, Yung K.↗

NASA NTRS2018

An Approach for Utilizing Power Flow Modeling for Simulations of Hybrid Electric Propulsion Systems

This paper describes an approach to creating simulations of the electric components for a hybrid electric propulsionsystem. The proposed modeling technique is based on power/load flow modeling and is designed to provide a modularframework that includes buses, lines, and other electrical components that can be connected together to form the electrical distribution system. The purpose of this paper is to detail an electric distribution system modeling technique and to demonstrate how these models may be integrated with turbomachinery simulations. These general modeling techniques were created to be utilized for system and control design studies. Additionally, steady-state and dynamic performance for a proposed model example is compared with data from a hardware in the loop simulation.

Power Flow↗

NASA NTRS1985

Description and calibration of the Langley Hypersonic CF4 tunnel: A facility for simulating low gamma flow as occurs for a real gas

The Langley Hypersonic CF4 Tunnel is a Mach 6 facility which simulates an important aspect of dissociative real-gas phenomena associated with the reentry of blunt vehicles, i.e., the decrease in the ratio of specific heats (gamma) that occurs within the shock layer of the vehicle. A general description of this facility is presented along with a discussion of the basic components, instrumentation, and operating procedure. Pitot-pressure surveys were made at the nozzle exit and downstream of the exit for reservoir temperatures from 1020 to 1495 R and reservoir pressures from 1000 to 2550 psia. A uniform test core having a diameter of circa 11 in. (0.55 times the nozzle-exit diameter) exists at the maximum value of reservoir pressure and temperature. The corresponding free-stream Mach number is 5.9, the unit Reynolds number is 4 x 10 to the 5th power per foot, the ratio of specific heats immediately behind a normal shock is 1.10, and the normal-shock density ratio is 12.6. When the facility is operated at reservoir temperatures below 1440 R, irregularities occur in the pitot-pressure profile within a small region about the nozzle centerline. These variations in pitot pressure indicate the existence of flow distrubances originating in the upstream region of the nozzle. This necessitates testing models off centerline in the uniform flow between the centerline region and either the nozzle boundary layer or the lip shock originating at the nozzle exit. Samples of data obtained in this facility with various models are presented to illustrate the effect of gamma on flow conditions about the model and the importance of knowing the magnitude of this effect.

Midden, Raymond E.↗

NASA NTRS1982

A Numerical simulation of transition in plane channel flow

A numerical simulation of the final stages of transition to turbulence in plane channel flow at a Reynolds number of 7500 is described. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two and three dimensional finite amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32 by 33 by 32 grid. Solutions indicate the existence of structures similar to those observed in the laboratory and which are characteristic of various stages of transition that lead to final breakdown. Details of the resulting flow field after breakdown indicate the evolution of streak-like formations found in turbulent flows. Although the flow field does approach a steady state (turbulent channel flow), implementation of subgrid-scale terms are necessary to obtain proper turbulent statistics.

Goglia, G.↗

DOE OSTI2016

Taylor flow in intermediate diameter channels: Simulation and hydrodynamic models

The two-phase Taylor flow pattern has been studied extensively. However, limited information is available for flows in intermediate diameter channels (5 ≲ Bo ≲ 40, or 6mm ≲ D ≲ 27mm for ambient gas–water flows), as found in air-lift and bubble pumps. Previous investigations have primarily evaluated Taylor flow models in terms of overall pressure drop, which incorporates hydrostatic and multiple hydrodynamic components. Thus, individual sub-models and sources of error could not be directly assessed. Here in this investigation, volume-of-fluid (VOF) based Taylor flow simulations are performed over a wide range of laminar and turbulent conditions in the intermediate Bond number regime (5 < Bo < 20, 250 < N f < 1000, and 20 < Re j < 8100). Results are applied to individually evaluate hydrodynamic sub-models for bubble-region frictional pressure drop gradient (∇p f,b ), slug frictional pressure drop gradient (∇p f,s ), and flow transition pressure loss (Δp trans ). Based on these results, recommendations are provided for selection of hydrodynamic sub-models. These hydrodynamic closure models are integrated with kinematic flow models to yield a complete intermediate Bond number Taylor flow formulation for which all submodels were independently validated. The resulting model achieves improved accuracy for predicting experimental liquid flow rates compared with previous Taylor flow models (81% of cases within 50% of measured flows rates).

42 ENGINEERING↗

NASA NTRS2000

SUPREM-DSMC: A New Scalable, Parallel, Reacting, Multidimensional Direct Simulation Monte Carlo Flow Code

An AFRL/NRL team has recently been selected to develop a scalable, parallel, reacting, multidimensional (SUPREM) Direct Simulation Monte Carlo (DSMC) code for the DoD user community under the High Performance Computing Modernization Office (HPCMO) Common High Performance Computing Software Support Initiative (CHSSI). This paper will introduce the JANNAF Exhaust Plume community to this three-year development effort and present the overall goals, schedule, and current status of this new code.

Campbell, David↗

NASA NTRS2007

SUPG Finite Element Simulations of Compressible Flows for Aerothermodynamic Applications

This viewgraph presentation reviews the Streamline-Upwind Petrov-Galerkin (SUPG) Finite Element Simulation. It covers the background, governing equations, weak formulation, shock capturing, inviscid flux discretization, time discretization, linearization, and implicit solution strategies. It also reviews some applications such as Type IV Shock Interaction, Forward-Facing Cavity and AEDC Sharp Double Cone.

Kirk, Benjamin S.↗

DOE OSTI2025

Modeling and simulation of multiphase flows

Multiphase flow can be found in a wide range of natural and industrial systems, such as sediment flow in rivers and coasts to pipeline transport, and reactors of all types. Computational fluid dynamics (CFD) modeling provides a rigorous physics-based method to visualize and quantify these flows and to assess their influence on their surroundings. This presentation will highlight computational research on two varied technologies in the energy sector: carbon capture and wellbore acoustics.

acoustic analysis↗