Search NASA⌕ Search

SEARCH · Search NASA

Results for “Backward-facing step”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗
NASA NTRS1993

Calculations of turbulent separated flows

A numerical study of incompressible turbulent separated flows is carried out by using two-equation turbulence models of the K-epsilon type. On the basis of realizability analysis, a new formulation of the eddy-viscosity is proposed which ensures the positiveness of turbulent normal stresses - a realizability condition that most existing two-equation turbulence models are unable to satisfy. The present model is applied to calculate two backward-facing step flows. Calculations with the standard K-epsilon model and a recently developed RNG-based K-epsilon model are also made for comparison. The calculations are performed with a finite-volume method. A second-order accurate differencing scheme and sufficiently fine grids are used to ensure the numerical accuracy of solutions. The calculated results are compared with the experimental data for both mean and turbulent quantities. The comparison shows that the present model performs quite well for separated flows.

Zhu, J.↗
NASA NTRS1993

A coupled multi-block solution procedure for spray combustion in complex geometries

Turbulent spray-combusting flow in complex geometries is presently treated by a coupled implicit procedure that employs finite-rate chemistry and real gas properties for combustion, as well as the stochastic separated model for spray and a multiblock treatment for complex geometries. Illustrative numerical tests conducted encompass a steady-state nonreacting backward-facing step flow, a premixed single-phase combustion flow, and spray combustion flow in a gas turbine combustor.

Chen, Kuo-Huey↗
NASA NTRS1993

Preconditioned Conjugate Gradient methods for low speed flow calculations

An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations are integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and the convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the lower-upper (L-U)-successive symmetric over-relaxation iterative scheme is more efficient than a preconditioner based on incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional line Gauss-Seidel relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.

Ajmani, Kumud↗
NASA NTRS1993

Selected computations of transonic cavity flows

An efficient diagonal scheme implemented in an overset mesh framework has permitted the analysis of geometrically complex cavity flows via the Reynolds averaged Navier-Stokes equations. Use of rapid hyperbolic and algebraic grid methods has allowed simple specification of critical turbulent regions with an algebraic turbulence model. Comparisons between numerical and experimental results are made in two dimensions for the following problems: a backward-facing step; a resonating cavity; and two quieted cavity configurations. In three-dimensions the flow about three early concepts of the stratospheric Observatory For Infrared Astronomy (SOFIA) are compared to wind-tunnel data. Shedding frequencies of resolved shear layer structures are compared against experiment for the quieted cavities. The results demonstrate the progress of computational assessment of configuration safety and performance.

Atwood, Christopher A.↗
NASA NTRS1994

A New Reynolds Stress Algebraic Equation Model

A general turbulent constitutive relation is directly applied to propose a new Reynolds stress algebraic equation model. In the development of this model, the constraints based on rapid distortion theory and realizability (i.e. the positivity of the normal Reynolds stresses and the Schwarz' inequality between turbulent velocity correlations) are imposed. Model coefficients are calibrated using well-studied basic flows such as homogeneous shear flow and the surface flow in the inertial sublayer. The performance of this model is then tested in complex turbulent flows including the separated flow over a backward-facing step and the flow in a confined jet. The calculation results are encouraging and point to the success of the present model in modeling turbulent flows with complex geometries.

Shih, Tsan-Hsing↗
NASA NTRS1994

Stability of separating subsonic boundary layers

The primary and subharmonic instabilities of separating compressible subsonic two-dimensional boundary layers in the presence of a two-dimensional roughness element on a flat plate are investigated. The roughness elements considered are humps and forward- and backward-facing steps. The use of cooling and suction to control these instabilities is studied. The similarities and differences between the instability characteristics of separating boundary layers and those of the boundary layer over a flat plate with a zero pressure gradient are pointed out and discussed. The theoretical results agree qualitatively and quantitatively with the experimental data of Dovgal and Kozlov. Cooling and suction decrease the growth rates of primary and subharmonic waves in the attached-flow regions but increase them in the separated-flow regions.

Masad, Jamal A.↗
NASA NTRS1994

Single point modeling of rotating turbulent flows

A model for the effects of rotation on turbulence is proposed and tested. These effects which influence mainly the rate of turbulence decay are modeled in a modified turbulent energy dissipation rate equation that has explicit dependence on the mean rotation rate. An appropriate definition of the rotation rate derived from critical point theory and based on the invariants of the deformation tensor is proposed. The modeled dissipation rate equation is numerically well behaved and can be used in conjunction with any level of turbulence closure. The model is applied to the two-equation kappa-epsilon turbulence model and is used to compute separated flows in a backward-facing step and an axisymmetric swirling coaxial jets into a sudden expansion. In general, the rotation modified dissipation rate model shows some improvements over the standard kappa-epsilon model.

Hadid, A. H.↗
NASA NTRS1994

Large-eddy simulation of a plane wake

Previously the theoretical development leading to the dynamic localization model (DLM) for large-eddy simulation (LES) was presented. The method has been successfully applied to isotropic turbulence, channel flow, and the flow over a backward-facing step. Here we apply the model to the computation of the temporally developing place wake. The two main objectives of this project are: (1) Use the model to perform an LES of a time developing plane wake and compare the results with direction numerical simulation (DNS) data to see if important statistical measures can be readily predicted, and to provide a relative evaluation of the several versions of the model in terms of predictive capability and cost; and (2) If the tests in (1) show that the model generates reliable predictions, then use the LES to study various aspects of the physics of turbulent wakes and mixing layers.

Ghosal, Sandip↗
NASA NTRS1995

Upwind-Differencing Computations Of Incompressible Flows

Report discusses application of upwind-differencing numerical-integration scheme in conjunction with method of pseudocompressibility to computations of flows of incompressible fluids, with emphasis on steady-state flows. Three test cases presented: driven flow in square cavity, flow over backward-facing step, and flow around circular cylinder.

Rogers, Stuart↗
NASA NTRS1996

Annual Research Briefs - 1996, Center for Turbulence Research

Partial contents include: (1) Simulation of jet impingement heat transfer with the k-e-V2 model; (2) A new methodology for turbulence modelers using Direct Numerical Simulation (DNS) database analysis; (3) A particle representation model for the deformation of homogeneous turbulence; (4) Modeling of inhomogeneous compressible turbulence using a two-scale statistical theory; (5) Lewis number effects on partially premixed flames; (6) Conditional moment closure of mixing and reaction in turbulent non-premixed combustion; (7) Dynamic models for LES of turbulent front propagation with a spectral method; (8) Jet noise using large eddy simulation; (9) Towards numerical simulations of trailing-edge aeroacoustics; (10) A model for the onset of vortex breakdown; (11) Compressing an elliptic vortex: transition to turbulence by tumble breakdown; (12) A new method for the adaptive control of vortex-wall interactions; (13) Symmetries in turbulent boundary layer flows; (14) Near-wall models in large eddy simulations of flow behind a backward-facing step; (15) Assessment of the subgrid-scale models at low and high Reynolds numbers; (16) Large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids; (17) Progress on LES of flow past a circular cylinder; (18) Experimental investigation of flow through an asymmetric plane diffuser; (19) Progress in the large-eddy simulation of an asymmetric plane diffuser; (20) A numerical study of self-similarity in a turbulent plane wake using large-eddy simulation; (21) On the generation of turbulent inflow conditions for boundary layer simulations, and (22) A new class of finite difference schemes.

COMPUTATIONAL FLUID DYNAMICS↗
NASA NTRS1995

Comparison of Implicit Schemes for the Incompressible Navier-Stokes Equations

For a computational flow simulation tool to be useful in a design environment, it must be very robust and efficient. To develop such a tool for incompressible flow applications, a number of different implicit schemes are compared for several two-dimensional flow problems in the current study. The schemes include Point-Jacobi relaxation, Gauss-Seidel line relaxation, incomplete lower-upper decomposition, and the generalized minimum residual method preconditioned with each of the three other schemes. The efficiency of the schemes is measured in terms of the computing time required to obtain a steady-state solution for the laminar flow over a backward-facing step, the flow over a NACA 4412 airfoil, and the flow over a three-element airfoil using overset grids. The flow solver used in the study is the INS2D code that solves the incompressible Navier-Stokes equations using the method of artificial compressibility and upwind differencing of the convective terms. The results show that the generalized minimum residual method preconditioned with the incomplete lower-upper factorization outperforms all other methods by at least a factor of 2.

Rogers, Stuart E.↗
NASA NTRS1992

Computation of Turbulent Recirculating Flow in Channels, and for Impingement Cooling

Fully elliptic forms of the transport equations have been solved numerically for two flow configurations. The first is turbulent flow in a channel with transverse rectangular ribs, and the second is impingement cooling of a plane surface. Both flows are relevant to proposed designs for active cooling of hypersonic vehicles using supercritical hydrogen as the coolant. Flow downstream of an abrupt pipe expansion and of a backward-facing step were also solved with various near-wall turbulence models as benchmark problems. A simple form of periodicity boundary condition was used for the channel flow with transverse rectangular ribs. The effects of various parameters on heat transfer in channel flow with transverse ribs and in impingement cooling were investigated using the Yap modified Jones and Launder low Reynolds number k-epsilon turbulence model. For the channel flow, predictions were in adequate agreement with experiment for constant property flow, with the results for friction superior to those for heat transfer. For impingement cooling, the agreement with experiment was generally good, but the results suggest that improved modelling of the dissipation rate of turbulence kinetic energy is required in order to obtain improved heat transfer prediction, especially near the stagnation point. The k-epsilon turbulence model was used to predict the mean flow and heat transfer for constant and variable property flows. The effect of variable properties for channel flow was investigated using the same turbulence model, but comparison with experiment yielded no clear conclusions. Also, the wall function method was modified for use in the variable properties flow with a non-adiabatic surface, and an empirical model is suggested to correctly account for the behavior of the viscous sublayer with heating.

Chang, Byong Hoon↗
NASA NTRS1987

A Diagonal Algorithm for the Method of Pseudocompressibility

The method of pseudocompressibility has been found to be an efficient method for obtaining a steady-state solution to the incompressible Navier-Stokes equations. Recent improvements to this method include the use of a diagonal scheme for the inversion of the equations equations at each iteration. The necessary transformations have been derived for the pseudocompressibility equations in generalized coordinates. The diagonal algorithm reduces the computing time necessary to obtain a steady-state solution by a factor of nearly three. Implicit viscous terms are maintained in the equations, and it has become possible to use fourth-order implicit dissipation. The steady-state solution is unchanged by the approximations resulting from the diagonalization of the equations. Computed results for flow over a two-dimensional backward-facing step and a three-dimensional cylinder mounted normal to a flat plate are presented for both the old and new algorithms. The computing efficiency of these algorithms are compared. Identical solutions are obtained from both algorithms which compare well with experimental results.

Rogers, Stuart E.↗
NASA NTRS1994

Turbulent Heat Transfer Characteristics in the Shear Layer of a Separated Flow

Experiments were performed to study the evolution of the heat transfer structure in a separated free shear layer region of an incompressible separated turbulent boundary layer flow behind a backward-facing step. While there is an abundance of velocity field measurements of separated flows, heat transfer measurements are rather scarce, thus limiting assessment of the heat transfer physics and its accurate modeling. The purpose of the paper is twofold: to improve an understanding of effects of flow separation on heat transfer characteristics, and to provide data for turbulence modeling and computation. The boundary layer upstream of the step was turbulent and fully developed. A constant temperature surface boundary condition was imposed upstream and downstream of the step for the heat transfer study. An internal mixing-layer like flow forms and grows from the step lip within the original boundary layer. The turbulent structure of the flow evolving downstream, however, does not switch immediately to that of a mixing layer over the entire shear layer thickness. Measurements of mean and fluctuating velocity and temperature fields indicate that the internal layer spreads gradually in the transverse direction while the outer part of the original boundary layer is effectively unperturbed. The results in this paper have not been previously reported.

Jovic, S.↗
NASA NTRS1994

Incompressible Navier-Stokes Computations with Heat Transfer

The existing pseudocompressibility method for the system of incompressible Navier-Stokes equations is extended to heat transfer problems by including the energy equation. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. Both forced and natural convection problems are examined. Numerical results from turbulent reattaching flow behind a backward-facing step will be compared against experimental measurements for the forced convection case. The validity of Boussinesq approximation to simplify the buoyancy force term will be investigated. The natural convective flow structure generated by heat transfer in a vertical rectangular cavity will be studied. The numerical results will be compared by experimental measurements by Morrison and Tran.

Kiris, Cetin↗
NASA NTRS2017

High-Spatial-Resolution OH and CH2O PLIF Visualization in a Dual-Mode Scramjet Combustor

A high-spatial-resolution planar laser-induced fluorescence (PLIF) imaging system was constructed and used to image a cavity-stabilized, premixed ethylene-air flame. The flame was created within a continuous flow, electrically-heated supersonic combustion facility consisting of a Mach 2 nozzle, an isolator with flush-wall fuel injectors, a combustor with a cavity flameholder of height 9 mm and optical access, and an extender. Tests were conducted at total temperature 1200 K, total pressure 300 kPa, equivalence ratio near 0.4 in the combustor, and Mach number near 0.6 in the combustor. A frequency-doubled Nd:YAG laser pumped a dye laser, which produced light at 283.55 nm. The beam was shaped into a light sheet with full width half-maximum 25 microns, which illuminated a streamwise plane that bisected the cavity. An intensified camera system imaged OH in this plane with a square 6.67 mm field of view and in-plane resolution 39 microns. Images were taken between the backward-facing step and 120 mm downstream of the step. OH structures as small as 110 microns were observed. CH2O was excited using 352.48 nm light; the smallest observed CH2O structures were approximately 200 microns wide. Approximately 15,000 images per species were processed and used to compute composite images.

Geipel, Clayton M.↗
NASA NTRS2017

High-Resolution OH and CH2O Visualization in a Premixed Cavity-Anchored Ethylene-Air Flame in a M = 0.6 Flowfield

OH and CH2O were imaged in a premixed, cavity-anchored, ethylene-air turbulent flame using a high resolution planar laser-induced fluorescence (PLIF) system. The electrically-heated, continuous flow facility (UVa Supersonic Combustion Facility, Configuration E) consisted of a Mach 2 nozzle, an isolator with fuel injectors, a test section with a cavity flame holder and optical access, and an extender. Standard test conditions comprised total temperature 1200 K, total pressure 300 kPa, local equivalence ratio near 0.4, and local Mach number near 0.6. OH PLIF data was also collected for a case with reduced total temperature and another with reduced equivalence ratio. OH and CH2O were excited in separate experiments with light sheets at 283.55 nm and 352.48 nm, respectively. A light sheet of approximate thickness 25 ?m illuminated the stream-wise midplane. This plane was imaged for 120 mm downstream of the backward-facing step. The intensified camera system imaged OH with magnification 1.97, a square 6.67 mm field of view, and in-plane resolution of 39 ?m. The smallest observed OH structures observed were approximately 100 ?m wide. The CH2O PLIF image signal was much weaker; the smallest observed structures were approximately 200 ?m wide. Composite fluorescence images were computed for the observed area.

Geipel, Clayton M.↗