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 73 records · Page 4

Evaluation of FIDAP on some classical laminar and turbulent benchmarks

The numerical accuracy of the fluid-dynamics code FIDAP is investigated by means of test computations on two-dimensional flows. The theoretical basis of the algorithm is briefly reviewed, and results for laminar flow in a wall-driven cavity and for both laminar and turbulent flows over a backward-facing step are presented in extensive tables and graphs and characterized in detail. Good agreement with published experimental data and previous computations is obtained for laminar flows using a version of FIDAP without streamlined upwinding (STU). For turbulent flows, the addition of STU is found to be useful, although it causes separation-zone reattachment length to be underpredicted.

Sohn, J. L.↗

Velocity-pressure integrated versus penalty finite element methods for high Reynolds number flows

Velocity-pressure integrated and consistent penalty finite element computations of high Reynolds number laminar flows are presented. In both methods the pressure has been interpolated using linear shape functions for a triangular element which is contained inside the biquadratic flow element. It has been shown previously that the pressure interpolation method, when used in conjunction with the velocity-pressure integrated method, yields accurate computational results for high-Reynolds-number flows. It is shown in this paper that use of the same pressure interpolation method in the consistent penalty finite element method yields computational results which are comparable to those of the velocity-pressure integrated method for both the velocity and the pressure fields. Accuracy of the two finite element methods has been demonstrated by comparing the computational results with available experimental data and/or fine grid finite difference computational results. Advantages and disadvantages of the two finite element methods are discussed on the basis of accuracy and convergence nature. Example problems considered include a lid-driven cavity flow of Reynolds number 10000, a laminar backward-facing step flow and a laminar flow through a nest of cylinders.

Kim, S.-W.↗

Numerical simulation of pressure oscillations in a ramjet combustor

Large-eddy simulations of compressible cold flow in a ramjet combustor configuration were performed. The objectives were to investigate the mechanisms for vortex-acoustic interaction in such a device and to develop a simulation method that can be extended to cases with combustion heat release to study combustion instability. From the simulation, it was found that the separated shear layer rolls up into concentrated vortices that merge to form large coherent structures. These vortices interact with the choked nozzle downstream to produce an axial acoustic dipole. The spectrum of the pressure fluctuation at the base of the backward-facing step shows that there are two types of oscillations: an acoustic resonant mode and a vortex-acoustic coupled mode. Based on the flow physics observed in the simulations, a simple one-dimensional model for the vortex-acoustic coupled mode was proposed. The eigenvalue problem based on this model was solved to obtain the frequency of the coupled mode.

Jou, Wen-Huei↗

A least-squares finite element method for incompressible Navier-Stokes problems

A least-squares finite element method, based on the velocity-pressure-vorticity formulation, is developed for solving steady incompressible Navier-Stokes problems. This method leads to a minimization problem rather than to a saddle-point problem by the classic mixed method, and can thus accommodate equal-order interpolations. This method has no parameter to tune. The associated algebraic system is symmetric, and positive definite. Numerical results for the cavity flow at Reynolds number up to 10,000 and the backward-facing step flow at Reynolds number up to 900 are presented.

Jiang, Bo-Nan↗

Understanding transition and turbulence through direct simulations

Direct simulations consist in solving the full Navier-Stokes equations, without any turbulence model, and describing all the detailed features of the flow. Usually the flows are three-dimensional and time-dependent and contain both coarse and fine structures, which makes the numerical task very challenging in terms of both the algorithm and the computational effort. Most of the work until now has involved spectral methods, which are highly accurate but not very flexible in terms of geometry or complex equations. For that reason, future work will also rely on high-order finite-difference or other methods. Direct simulations complement experimental work, and both contribute to the theory and the empirical knowledge of turbulence. Once such a simulation has been shown to be accurate, the flow field is completely known in three dimensions and time, including the pressure, the vorticity and any other quantity. On the other hand, most simulations to date solved the incompressible equations in rather simple geometries, and direct simulations will always be limited to moderate Reynolds numbers. Extensive simulations have been conducted in homogeneous turbulence, channel flows, boundary layers, and mixing layers. Much effort is devoted to addressing flows with compressibility and chemical reactions, and to new geometries such as a backward-facing step.

Spalart, P. R.↗

Computations of three-dimensional steady and unsteady viscous incompressible flows

The INS3D family of computational fluid dynamics computer codes is presented. These codes are used to as tools in developing and assessing algorithms for solving the incompressible Navier-Stokes equations for steady-state and unsteady flow problems. This work involves applying the codes to real-world problems involving complex three-dimensional geometries. The algorithms utilized include the method of pseudocompressibility including both central and upwind differencing, several types of artificial dissipation schemes, approximate factorization, and an implicit line-relaxation scheme. These codes have been validated using a wide range of problems including flow over a backward-facing step, driven cavity flow, flow through various types of ducts, and steady and unsteady flow over a circular cylinder. Many diverse flow applications have been solved using these codes including parts of the Space Shuttle Main Engine, problems in naval hydrodynamics, low-speed aerodynamics, and biomedical fluid flows. The presentation details several of these, including the flow through a Space Shuttle Main Engine inducer, vortex shedding behind a circular cylinder, and flow through an artificial heart.

Kwak, Dochan↗

A steadying effect of acoustic excitation on transitory stall

The effect of acoustic excitation on a class of separated flows with a transitional boundary layer at the point of separation is considered. Experimental results on the flow over airfoils, a two-dimensional backward-facing step, and through large angle conical diffusers are presented. In all cases, the separated flow undergoes large amplitude fluctuations, much of the energy being concentrated at unusually low frequencies. In each case, an appropriate high frequency acoustic excitation is found to be effective in reducing the fluctuations substantially. The effective excitation frequency scales on the initial boundary layer thickness and the effect is apparently achieved through acoustic tripping of the separating boundary layer.

Zaman, K. B. M. Q.↗

The effects of streamline curvature and swirl on turbulent flows in curved ducts

A technique for improving the numerical predictions of turbulent flows with the effect of streamline curvature is developed. Separated flows, the flow in a curved duct, and swirling flows are examples of flow fields where streamline curvature plays a dominant role. A comprehensive literature review on the effect of streamline curvature was conducted. New algebraic formulations for the eddy viscosity incorporating the kappa-epsilon turbulence model are proposed to account for various effects of streamline curvature. The loci of flow reversal of the separated flows over various backward-facing steps are employed to test the capability of the proposed turbulence model in capturing the effect of local curvature. The inclusion of the effect of longitudinal curvature in the proposed turbulence model is validated by predicting the distributions of the static pressure coefficients in an S-bend duct and in 180 degree turn-around ducts. The proposed turbulence model embedded with transverse curvature modification is substantiated by predicting the decay of the axial velocities in the confined swirling flows. The numerical predictions of different curvature effects by the proposed turbulence models are also reported.

Cheng, Chih-Hsiung↗

Pressure algorithm for elliptic flow calculations with the PDF method

An algorithm to determine the mean pressure field for elliptic flow calculations with the probability density function (PDF) method is developed and applied. The PDF method is a most promising approach for the computation of turbulent reacting flows. Previous computations of elliptic flows with the method were in conjunction with conventional finite volume based calculations that provided the mean pressure field. The algorithm developed and described here permits the mean pressure field to be determined within the PDF calculations. The PDF method incorporating the pressure algorithm is applied to the flow past a backward-facing step. The results are in good agreement with data for the reattachment length, mean velocities, and turbulence quantities including triple correlations.

Anand, M. S.↗

Analysis of two-equation turbulence models for recirculating flows

The two-equation kappa-epsilon model is used to analyze turbulent separated flow past a backward-facing step. It is shown that if the model constraints are modified to be consistent with the accepted energy decay rate for isotropic turbulence, the dominant features of the flow field, namely the size of the separation bubble and the streamwise component of the mean velocity, can be accurately predicted. In addition, except in the vicinity of the step, very good predictions for the turbulent shear stress, the wall pressure, and the wall shear stress are obtained. The model is also shown to provide good predictions for the turbulence intensity in the region downstream of the reattachment point. Estimated long time growth rates for the turbulent kinetic energy and dissipation rate of homogeneous shear flow are utilized to develop an optimal set of constants for the two equation kappa-epsilon model. The physical implications of the model performance are also discussed.

Thangam, S.↗

A steadying effect of acoustic excitation on transitory stall

The effect of acoustic excitation on a class of separated flows with a transitional boundary layer at the point of separation is considered. Experimental results on the flow over airfoils, a two-dimensional backward-facing step, and through large angle conical diffusers are presented. In all cases, the separated flow undergoes large amplitude fluctuations, much of the energy being concentrated at unusually low frequencies. In each case, an appropriate high frequency acoustic excitation is found to be effective in reducing the fluctuations substantially. The effective excitation frequency scales on the initial boundary layer thickness and the effect is apparently achieved through acoustic tripping of the separating boundary layer.

Zaman, K. B. M. Q.↗

An experimental investigation of the separating/reattaching flow over a backstep

Reynolds number effects on the evolution of the turbulence structure in the recovery region of a turbulent separated flow generated by a backward-facing step were studied experimentally. Five different flow conditions and four Reynolds numbers were investigated. All single point measurements were conducted on the plane of symmetry of the wind tunnel starting at approximately x/h = 10 for each flow condition. The upper limit of the measuring domain depends upon the physical size of the tunnel and ranged from 50 to 114 step heights. It is apparent that most of the boundary layer recovery was captured by the experiment. Normal stresses, shear stress, mixing length scales and eddy viscosity are shown, while the higher order velocity products related to the diffusion of the turbulence energy are planned for latter stages of the data analysis.

Jovic, Srboljub↗

On turbulent flows dominated by curvature effects

A technique for improving the numerical predictions of turbulent flows with the effect of streamline curvature is developed. Separated flows and the flow in a curved duct are examples of flowfields where streamline curvature plays a dominant role. New algebraic formulations for the eddy viscosity incorporating the k-epsilon turbulence model are proposed to account for various effects of streamline curvature. The loci of flow reversal of the separated flows over various backward-facing steps are employed to test the capability of the proposed turbulence model in capturing the effect of local curvature.

Cheng, G. C.↗

A least-squares finite element method for incompressible Navier-Stokes problems

A least-squares finite element method, based on the velocity-pressure-vorticity formulation, is developed for solving steady incompressible Navier-Stokes problems. This method leads to a minimization problem rather than to a saddle-point problem by the classic mixed method and can thus accommodate equal-order interpolations. This method has no parameter to tune. The associated algebraic system is symmetric, and positive definite. Numerical results for the cavity flow at Reynolds number up to 10,000 and the backward-facing step flow at Reynolds number up to 900 are presented.

Jiang, Bo-Nan↗

LDV measurements on a rectangular wing with a simulated glaze ice accretion

LDV measurement results are presented for the upper surface of a rectangular semispan wing with and without simulated glaze ice accretion. Inspection of the model centerline flow field indicates that a large region of reverse flow exists aft of the ice horn on the iced model. At alpha = 0 deg, this region extends to 7 percent chord, while at alpha = 4 deg the bubble grows to more than 12 percent chord. At alpha = 8 deg, the time-averaged separation bubble is measured well beyond 50 percent chord. Experimental and computational flow visualization support these findings. The flow in the vicinity of the ice shape contains many of the features of flow over a backward-facing step.

Khodadoust, A.↗

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 is 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 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 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↗

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.↗

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↗