SEARCH · Engineering Papers

Results for “turbulent flows”

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,153 records · Page 64

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Large-eddy simulation of compressible turbulence

The increase in the range of length scales with increasing Reynolds number limits the direct simulation of turbulent flows to relatively simple geometries and low Reynolds numbers. However, since most flows of engineering interest occur at much higher Reynolds number than is currently within the capabilities of full simulation, prediction of these flow fields can only be obtained by solving some suitably-averaged set of governing equations. In the traditional Reynolds-averaged approach, the Navier-Stokes equations are averaged over time. This in turn yields correlations between various turbulence fluctuations. It is these terms, e.g. the Reynolds stresses, for which a turbulence model must be derived. Turbulence modeling of incompressible flows has received a great amount of attention in the literature. An area of research that has received comparatively less attention is the modeling of compressible turbulent flows. An approach to simulating compressible turbulence at high Reynolds numbers is through the use of Large-Eddy Simulation (LES). In LES the dependent variables are decomposed into a large-scale (resolved) component and a sub-grid scale component. It is the small-scale components of the velocity field which are presumably more homogeneous than the large scales and, therefore, more easily modeled. Thus, it seems plausible that simpler models, which should be more universal in character than those employed in second-order closure schemes, may be developed for LES of compressible turbulence. The objective of the present research, therefore, is to explore models for the Large-Eddy Simulation of compressible turbulent flows. Given the recent successes of Zeman in second order closure modeling of compressible turbulence, model development was guided by principals employed in second-order closures.

Squires, Kyle D.↗

NASA NTRS2004

Assessment of LAURA for Laminar Supersonic Shallow Cavities

The ability of the Laura flow solver to predict local heating augmentation factors for shallow cavities is assessed. This assessment is part of a larger e ort within the Space Shuttle return-to-flight program to develop technologies to support on-orbit tile repair decisions. The comparison is made against global phosphor thermography images taken in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel. The cavities are rectangular in shape, with lengths L/H of 14 20 and depths H/ of 1.1 5.2. The fully laminar results, for Re = 300, show good agreement between the data sets. For Re = 503, the wind tunnel data indicates boundary layer transition with turbulent flow both within and downstream of the cavity. The turbulent flow structures are significantly di erent from the laminar predictions, with order of magnitude increases in the heating augmentations. Because of the di erent flow structures, no simple bump factor can be used to correct the laminar calculations to account for the turbulent heating levels. A fine gradation in wind tunnel cases will be required to clearly delineate the laminar-to-turbulent transition point, and hence the limits of applicability of the laminar numerical approach.

Wood, William A.↗

NASA NTRS2010

Hypersonic Navier-Stokes Comparisons to Orbiter Flight Data

During the STS-119 flight of Space Shuttle Discovery, two sets of surface temperature measurements were made. Under the HYTHIRM program3 quantitative thermal images of the windward side of the Orbiter with a were taken. In addition, the Boundary Layer Transition Flight Experiment 4 made thermocouple measurements at discrete locations on the Orbiter wind side. Most of these measurements were made downstream of a surface protuberance designed to trip the boundary layer to turbulent flow. In this paper, we use the US3D computational fluid dynamics code to simulate the Orbiter flow field at conditions corresponding to the STS-119 re-entry. We employ a standard two-temperature, five-species finite-rate model for high-temperature air, and the surface catalysis model of Stewart.1 This work is similar to the analysis of Wood et al . 2 except that we use a different approach for modeling turbulent flow. We use the one-equation Spalart-Allmaras turbulence model8 with compressibility corrections 9 and an approach for tripping the boundary layer at discrete locations. In general, the comparison between the simulations and flight data is remarkably good

Candler, Graham V.↗

NASA NTRS1983

Modeling of turbulent separated flows for aerodynamic applications

Steady, high speed, compressible separated flows modeled through numerical simulations resulting from solutions of the mass-averaged Navier-Stokes equations are reviewed. Emphasis is placed on benchmark flows that represent simplified (but realistic) aerodynamic phenomena. These include impinging shock waves, compression corners, glancing shock waves, trailing edge regions, and supersonic high angle of attack flows. A critical assessment of modeling capabilities is provided by comparing the numerical simulations with experiment. The importance of combining experiment, numerical algorithm, grid, and turbulence model to effectively develop this potentially powerful simulation technique is stressed.

Marvin, J. G.↗

NASA NTRS1986

Modeling of turbulent separated flows for aerodynamic applications

Marvin, J. G.↗

NASA NTRS1985

Experimental and numerical investigation of supersonic flow through a square duct

Steady, developing, adiabatic supersonic flow in a square duct is investigated for an inlet Mach number of 3.91 and a unit Reynolds number of 1.8 x 10 to the 6th/m. The numerical results for laminar flow show that two secondary flow cells develop in the near vicinity of the corner which are centered about the corner bisector and distort the primary flow in this region. For turbulent flow, the experimental results indicate that two secondary flow cells also develop about the corner bisector, but are directed in an opposite sense to that observed for the laminar case. Numerical results based on the Baldwin-Lomax model show that this model is incapable of predicting turbulence-generated secondary flow cells. For a suitable choice of constants, the Gessner-Emery model is able to predict the strength of these cells, but is deficient with respect to predicting their positions in the flow and their distorting influence on the primary flow. These observations are based on comparisons made in this paper between predicted and measured total pressure contours, cross flow velocity profiles, and local wall shear stress distributions.

Davis, D. O.↗

NASA NTRS1982

Reentry heating analysis of space shuttle with comparison of flight data

Surface heating rates and surface temperatures for a space shuttle reentry profile were calculated for two wing cross sections and one fuselage cross section. Heating rates and temperatures at 12 locations on the wing and 6 locations on the fuselage are presented. The heating on the lower wing was most severe, with peak temperatures reaching values of 1240 C for turbulent flow and 900 C for laminar flow. For the fuselage, the most severe heating occured on the lower glove surface where peak temperatures of 910 C and 700 C were calculated for turbulent flow and laminar flow, respectively. Aluminum structural temperatures were calculated using a finite difference thermal analyzer computer program, and the predicted temperatures are compared to measured flight data. Skin temperatures measured on the lower surface of the wing and bay 1 of the upper surface of the wing agreed best with temperatures calculated assuming laminar flow. The measured temperatures at bays two and four on the upper surface of the wing were in quite good agreement with the temperatures calculated assuming separated flow. The measured temperatures on the lower forward spar cap of bay four were in good agreement with values predicted assuming laminar flow.

Gong, L.↗

NASA NTRS1992

Volcanic recycling of carbonate deposits on Mars

One question of great interest to those who study the evolution of the Martian atmosphere is: if there was an early, dense atmosphere that was removed, is there any mechanism that could restore it? In the case of an atmosphere removed largely by the formation of carbonates, the only obvious means of restoring it is by the thermal decomposition of the carbonates. Decomposition of carbonates under turbulently flowing lava holds great promise as a means of resupplying the atmosphere with CO2. Huppert and colleagues have modeled the emplacement of terrestrial komatiite flows and found that komatiites, even when flowing over previously emplaced and cooled komatiite flows, could melt and erode this rock to a significant depth. Based on this work, I have begun modeling the erosion of Martian carbonate deposits under turbulently flowing, komatiitic lava. Initial results from this modeling indicate that a high-volume lava flow, emerging at a temperature of, say, 1600 degrees, is capable of eroding several meters of carbonate deposits per day. If such a flow is active for a hundred days, several hundreds of meters of carbonate could be decomposed. If this process occurred over a large area, a bar or more of CO2 could be injected back into the atmosphere over an extremely short period of time. The implications of such an occurrence are intriguing. For instance, if a relatively late pulse of volcanism (such as is suggested by Frey) were to cause a large flow of lava over carbonate deposits in the northern lowlands, the resulting pulse of CO2 into the atmosphere could conceivably restore the climate to one in which liquid water could exist on the surface, or ice could flow.

Schaefer, M. W.↗

NASA NTRS1993

Computation of flows in a turn-around duct and a turbine cascade using advanced turbulence models

Numerical investigation has been carried out to evaluate the capability of the Algebraic Reynolds Stress Model (ARSM) and the Nonlinear Stress Model (NLSM) to predict strongly curved turbulent flow in a turn-around duct (TAD). The ARSM includes the near-wall damping term of pressure-strain correlation phi(sub ij,w), which enables accurate prediction of individual Reynolds stress components in wall flows. The TAD mean flow quantities are reasonably well predicted by various turbulence models. The ARSM yields better predictions for both the mean flow and the turbulence quantities than the NLSM and the k-epsilon (k = turbulent kinetic energy, epsilon = dissipation rate of k) model. The NLSM also shows slight improvement over the k-epsilon model. However, all the models fail to capture the recovery of the flow from strong curvature effects. The formulation for phi(sub ij,w) appears to be incorrect near the concave surface. The hybrid k-epsilon/ARSM, Chien's k-epsilon model, and Coakley's q-omega (q = the square root of k, omega = epsilon/k) model have also been employed to compute the aerodynamics and heat transfer of a transonic turbine cascade. The surface pressure distributions and the wake profiles are predicted well by all the models. The k-epsilon model and the k-epsilon/ARSM model provide better predictions of heat transfer than the q-omega model. The k-epsilon/ARSM solutions show significant differences in the predicted skin friction coefficients, heat transfer rates and the cascade performance parameters, as compared to the k-epsilon model. The k-epsilon/ARSM model appears to capture, qualitatively, the anisotropy associated with by-pass transition.

Lakshminarayana, B.↗

NASA NTRS1982

On the use of a global hydrogen-air combustion model in the calculation of turbulent reacting flows

A two-step global model describing the combustion of hydrogen in air at one atmosphere pressure is developed by comparing the temperature histories obtained from a 28-reaction H-O mechanism. Using criteria discussed herein, good agreement is obtained in the range of initial mixture temperatures from 1000K to 2000K and in the range of equivalence ratios from 0.2 to 2.0. The two-step global model is compared with the results obtained using an eight-step, H-O reaction mechanism in a computer program describing the turbulent diffusion of hydrogen in supersonic, axisymmetric and two-dimensional reacting flows. Comparisons of profiles of temperature, fuel mass flux, and pitot pressure are presented. The global model is judged to be adequate in flows which are not dominated by long ignition delay times.

Rogers, R. C.↗

NASA NTRS1991

Heat transfer in a compact heat exchanger containing rectangular channels and using helium gas

Development of a National Aerospace Plane (NASP), which will fly at hypersonic speeds, require novel cooling techniques to manage the anticipated high heat fluxes on various components. A compact heat exchanger was constructed consisting of 12 parallel, rectangular channels in a flat piece of commercially pure nickel. The channel specimen was radiatively heated on the top side at heat fluxes of up to 77 W/sq cm, insulated on the back side, and cooled with helium gas flowing in the channels at 3.5 to 7.0 MPa and Reynolds numbers of 1400 to 28,000. The measured friction factor was lower than that of the accepted correlation for fully developed turbulent flow, although the uncertainty was high due to uncertainty in the channel height and a high ratio of dynamic pressure to pressure drop. The measured Nusselt number, when modified to account for differences in fluid properties between the wall and the cooling fluid, agreed with past correlations for fully developed turbulent flow in channels. Flow nonuniformity from channel-to-channel was as high as 12 pct above and 19 pct below the mean flow.

Olson, D. A.↗

NASA NTRS1980

The influence of a high pressure gradient on unsteady velocity perturbations in the case of a turbulent supersonic flow

The amplification or reduction of unsteady velocity perturbations under the influence of strong flow acceleration or deceleration was studied. Supersonic flows with large velocity, pressure gradients, and the conditions in which the velocity fluctuations depend on the action of the average gradients of pressure and velocity rather than turbulence, are described. Results are analyzed statistically and interpreted as a return to laminar process. It is shown that this return to laminar implies negative values in the turbulence production terms for kinetic energy. A simple geometrical representation of the Reynolds stress production is given.

Dussauge, J. P.↗

NASA NTRS2001

The Summer 1997 Eruption at Pillan Patera on Io: Implications for Ultrabasic Lava Flow Emplacement

Galileo data and numerical modeling were used to investigate the summer 1997 eruption at Pillan Patera on Io. This event, now defined as 'Pillanian' eruption style, included a high-temperature (greater than 1600 C), possibly ultrabasic, 140-km-high plume eruption that deposited dark, orthopyroxene-rich pyroclastic material over greater than 125,000 sq km, followed by emplacement of dark flow-like material over greater than 3100 sq km to the north of the caldera. We estimate that the high-temperature, energetic episode of this eruption had a duration of 52- 167 days between May and September 1997, with peak eruption temperatures around June 28, 1997. Galileo 20 m/pixel images of part of the Pillan flow field show a widespread, rough, pitted surface that is unlike any flow surface we have seen before. We suggest that th.s surface may have resulted from (1) a fractured lava crust formed during rapid, low-viscosity lava surging, perhaps including turbulent flow emplacement; (2) disruption of the lava flow by explosive interaction with a volatile-rich substrate: or (3) a combination of 1 and 2 with or without accumulation of pyroclastic materials on the surface. Well-developed flow lobes are observed, suggesting that this is a relatively distal part of the flow field. Shadow measurements at flow margins indicate a thickness of approx. 8-10 m. We have modeled the emplacement of putative ultrabasic flows from the summer 1997 Pillan eruption using constraints from new Galileo data. Results suggest that either laminar sheet flows or turbulent channelized flows could have traveled 50-150 km on a flat. unobstructed surface, which is consistent with the estimated length of the Pillan flow field (approx. 60 km). Our modeling suggests low thermal erosion rates (less than 0.1 m/d), and that the formation of deep (greater than 20 m) erosion channels was unlikely, especially distal to the source. We calculate a volumetric flow rate of approx. 2-7 x l0(exp 3) cu m/s, which is greater than those for typical Mauna Loa/Kilauea flows but comparable to those for the (1783) Laki eruption and the inferred flow rates of the Roza flows in the Columbia River flood basalts. The differences in ultrabasic eruption styles on Earth and Io appear to be controlled by the different eruption environments: Plumes at sites of ultrabasic eruptions on Io suggest strong magma-volatile: interactions on a low-gravity body lacking an atmosphere, whereas the geology at sites of komatiite eruptions on Earth suggest mostly submarine emplacement of thick flows with a pronounced lack of subaerial explosive activity.

Williams, David A.↗

NASA NTRS2001

The Summer 1997 Eruption at Pillan Patera on Io: Implications for Ultrabasic Lava Flow Emplacement

Galileo data and numerical modeling were used to investigate the summer 1977 eruption at Pillan Patera on Io. This event, now defined as "Pillanian" eruption style, included a high-temperature (greater than 1600 C), possible ultrabasic , 140-km-high plume eruption that deposited dark, orthopyroxene-rich pyroclastic material over greater than 125,000 sq km, followed by emplacement of dark flow-like material over greater than 3100 sq km to the north of the caldera. We estimate that the high-temperature, energetic episode of this eruption had a duration of 52 - 167 days between May and September 1997, with peak eruption temperatures around June 28, 1997. Galileo 20 m/pixel images of part of the Pillan flow field show a wide-spread, rough, pitted surface that is unlike any flow surface we have seen before. We suggest that this surface may have resulted from: 1. A fractured lava crust formed during rapid, low-viscosity lava surging, perhaps including turbulent flow emplacement. 2. Disruption of the lava flow by explosive interaction with a volatile-rich substrate. or 3. A combination of 1 and 2 with or without accumulation of pyroclastic material on the surface. Well-developed flow lobes are observed, suggesting that this is a relatively distant part of the flow field.Shadow measurements at flow margins indicate a thickness of-8 - 10 m. We have modeled the emplacement of putative ultrabasic flow from the summer 1997 Pillan eruption using constraints from new Galileo data. Results suggest that either laminar sheet flows or turbulent channelized flows could have traveled 50 - 150 km on a flat, unobstructed surface, which is consistent with the estimated length of the Pillan flow field (approx. 60 km). Our modeling suggests low thermal erosion rates (less than 4.1 m/d), and that the formation of deep (greater than 20 m) erosion channels was unlikely, especially distal to the source. We calculate a volumetric flow rate of approx. 2 - 7 x 10(exp 3)cu m/s, which is greater than those for typical Mauna Loa/Kilaueaq flows but comparable to those for the (1783) Laki eruption and the inferred flow rates of the Roza flows in the Columbia River flood basalts. The differences in ultrabasic eruption styles on Earth and Io appear to be controlled by the different eruption environments; Plumes at sites of ultrabasic eruptions on Io suggest strong magma-volatile interactions on a low-gravity body lacking an atmosphere, whereas the geology at sites of komatiite eruptions on Earth suggest mostly submarine emplacement of thick flows with a pronounced lack of subaerial explosive activity.

Williams, David A.↗

NASA NTRS1979

Computation of supersonic viscous flows over ogive-cylinders at angle of attack

The parabolic Navier-Stokes (PNS) marching finite-difference method is applied to 3-D viscous flow over pointed ogive-cylinders, and to turbulent flow over a cone. Ogive computations were performed using the new technique recently reported by Vigneron, Rakich, and Tannehill. Comparison is made with experiment and inviscid computations. The present results show that this method, which neglects part of the pressure gradient in the x-momentum equation, is nevertheless valid for flows with a strong favorable pressure gradient. In addition, turbulent separated flow over a cone has been computed using the older PNS code due to Lubard and Helliwell. It is found that one must freeze the turbulent eddy-viscosity model upstream of 3-D separation to get agreement with experiment.

Rakich, J. V.↗

NASA NTRS1996

Numerical Simulation of Flow in a Whirling Annular Seal and Comparison With Experiments

The turbulent flow field inside a whirling annular seal was simulated by using SCISEAL, a three-dimensional computational fluid dynamics code. The rotor center described a circular synchronous whirl. A rotating frame transformation was used to make the problem quasi-steady. The flow field at an axial Reynolds number of 24,000 and a Taylor number of 6600 was simulated. The standard kappa-epsilon model with wall functions and the low-Reynolds-number model were used to treat turbulence. An experimentally measured velocity field was used at the inlet boundary. Numerical predictions of the velocities and the stator wall pressures compared well with experimental data. Both turbulence models yielded nearly the same results. The capability of the SCISEAL code to analyze this complex flow field was demonstrated; the isotropic turbulence models performed adequately on this nonisotropic turbulence flow.

Athavale, M. M.↗

NASA NTRS1998

Current Status on the use of Parallel Computing in Turbulent Reacting Flow Computations Involving Sprays, Monte Carlo PDF and Unstructured Grids

The state of the art in multidimensional combustor modeling as evidenced by the level of sophistication employed in terms of modeling and numerical accuracy considerations, is also dictated by the available computer memory and turnaround times afforded by present-day computers. With the aim of advancing the current multi-dimensional computational tools used in the design of advanced technology combustors, a solution procedure is developed that combines the novelty of the coupled CFD/spray/scalar Monte Carlo PDF (Probability Density Function) computations on unstructured grids with the ability to run on parallel architectures. In this approach, the mean gas-phase velocity and turbulence fields are determined from a standard turbulence model, the joint composition of species and enthalpy from the solution of a modeled PDF transport equation, and a Lagrangian-based dilute spray model is used for the liquid-phase representation. The gas-turbine combustor flows are often characterized by a complex interaction between various physical processes associated with the interaction between the liquid and gas phases, droplet vaporization, turbulent mixing, heat release associated with chemical kinetics, radiative heat transfer associated with highly absorbing and radiating species, among others. The rate controlling processes often interact with each other at various disparate time 1 and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and liquid phase evaporation in many practical combustion devices.

Raju, M. S.↗

NASA NTRS1972

Diffusion flame in homologous turbulent shear flows.

The simplifed statistical theory developed previously is employed to analyze the equilibrium and near-equilibrium combustion of initially unmixed reactants. It is found that the flame zone in the limit of large Damkohler numbers is very thick and is of the order of the local integral scale of turbulence. This is in contrast to the existing phenomenological theories which predict the infinitesimally thin flame sheet, in the same limit, as it is with the laminar diffusion flame. Qualitative agreements with the available experimental results are shown. It is found that singularities exist at the edges of the flame which are removed as Damkohler number is reduced. Also, it is found that the heat transfer may take place against the local mean temperature gradient in certain regions within the flame.

Chung, P. M.↗