Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reynolds-averaged Navier Stokes”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Time-accurate Navier-Stokes computations of self-excited two-dimensional unsteady cavity flows

Time-accurate Navier Stokes computations were performed to study self-excited oscillatory shear layer behavior over a two-dimensional rectangular cavity and to gain insight into the vortex-acoustic coupling that sustains these oscillations. Results were obtained solving the full, unsteady compressible Navier-Stokes equations for Reynolds-averaged turbulent air flow past a cavity with L/d = 2. Conditions duplicated test conditions in experiments by Rossiter at Mach numbers varied from 0.6 to 1.2. The computational scheme was the USA 2D series code developed by Rockwell's Science Center with an algebraic turbulence model for the shear layer and recirculating region. Two physical conditions were simulated: an isolated cavity situated (1) in an infinite plane wall and (2) in an infinite channel with an opposing parallel wall. Computed results for (1) showed excellent agreement with Rossiter's measured amplitudes and frequencies. A tuned acoustic duct resonance effect was observed in (2), demonstrating a controlled change in the acoustic feedback mechanism and resultant amplification in computed oscillations.

Dougherty, N. Sam↗
NASA NTRS1992

Structure of supersonic turbulent flow past a swept compression corner

The structure of the shock wave/turbulent boundary-layer interaction generated by a 3D swept compression corner has been investigated through a combined experimental and theoretical research program. The flowfield geometry is defined by the streamwise compression angle alpha and the sweep angle lambda of the corner. The present study examines two different configurations, namely (alpha, lambda) = (24 deg, 40 deg) and (24 deg, 60 deg) at Mach 3 and Re sigma infinity about 9 x 10 exp 5. The theoretical model is the 3D Reynolds-averaged compressible Navier-Stokes equations with turbulence incorporated using a turbulent eddy viscosity. The calculated flowfields display general agreement with experimental data for surface pressure and good agreement with experimental flowfield profiles of pitot pressure and yaw angle. The principal feature of the flowfield is a large vortical structure approximately aligned with the corner. The entrainment of incoming fluid into the vortical structure is strongly affected by the sweep angle lambda. Viscous (turbulent and molecular) effects appear to be important only in the immediate vicinity of the surface and in an isolated region within the interaction and near the corner.

Knight, Doyle D.↗
NASA NTRS1992

Quasiconical flowfield structure of the three-dimensional single fin interaction

A series of conical and three-dimensional computations have been performed for the swept oblique shock wave/turbulent boundary-layer interaction generated by a 20-deg sharp fin at Mach 4 and freestream Reynolds number of 2.18 x 10 exp 5 based on the incoming boundary-layer thickness. The Reynolds-averaged compressible Navier-Stokes equations are employed with turbulence incorporated using the Baldwin-Lomax and Jones-Launder models. The computed results are basically similar for both turbulence models and display general agreement with experimental data for surface pressure and surface flow direction, although underestimating the size of the primary vortex. The computed three-dimensional flowfield displays quasiconical behavior of the surface pressure, surface flow direction, and flowfield contours of static pressure, density, and Mach number over the extent of the computational domain except for an inception region near the fin leading edge. Certain features of the flowfield model are not observed in the computations, namely, a 'normal' shock near the attachment line, transonic shocklets in the expansion region, and secondary separation. The absence of these features in the computation is believed to be indirectly attributable to limitations in the turbulence models.

Knight, Doyle D.↗
NASA NTRS1995

Effect of turbulence models on criticality conditions in swirling flows

The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon, RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.

Gatski, Thomas B.↗
NASA NTRS1995

Validation of the NPARC code for nozzle afterbody flows at transonic speeds

The NPARC code, a Reynolds-averaged full Navier-Stokes code, was validated for nozzle afterbody (boatail) flow fields at transonic speeds. The flow fields about three geometries were studied: an axisymmetric nozzle with attached flow; an axisymmetric nozzle with separated flow: and a two-dimensional (rectangular) nozzle with separated flow. Three turbulence models, Baldwin-Lomax, Baldwin-Barth, and Chien k-epsilon, were used to determine the effect of turbulence model selection on the flow field solution. Static pressure distributions on the nozzle surfaces and pitot pressure measurements in the exhaust plume were examined. Results from the NPARC code compared very well with experimental data for all cases. For attached flow fields, the effect of the turbulence models showed no discernable differences. The Baldwin-Barth model yielded better results than either the Chien k-epsilon or the Baldwin-Lomax model for separated flow fields.

Debonis, James R.↗
NASA NTRS1996

Improved Models for the Porous Surface with Passive Control

A computational investigation of afterbody flow using a passive control method is conducted. The passive control method consists of a porous surface placed over a plenum. The purpose of the passive control method is to exploit the adverse pressure gradient present in afterbody flow in an attempt to reduce boundary layer separation and afterbody drag. Four different porous wall models are used to model the transpiration velocity in the region of passive control. A three-dimensional, time-dependent, Reynolds-averaged, simplified Navier-Stokes solver, PAB3D, is used to simulate afterbody flow with and without passive control. Three afterbody configurations with boat-tail angles of 10, 20, and 30 deg. are used to obtain two-dimensional solutions with a freestream Mach number of 0.6 and nozzle pressure ratio of 6. The region of passive control was initially placed from 20-60% of the nozzle length. The effect of the porous placement and porous extent is also studied. Baseline (no porosity) two-dimensional solutions are qualitatively similar to experimental data but under-predict the magnitude of the pressure recovery. Results for the subsonic solutions show losses in the pressure recovery for some cases with passive control. Three-dimensional effects are also investigated and seen to be very significant. Three-dimensional baseline solutions, for both sub- and super-critical freestream Mach numbers, compare very favorably with the experimental data in comparison to the two-dimensional solution. Future work is required to examine three-dimensional afterbody flows with passive porosity.

Kraushaar, Sandra L.↗
NASA NTRS1986

Thermodynamic evaluation of transonic compressor rotors using the finite volume approach

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite volume, time marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier Stokes equations. The entire calculation is performed in the physical domain. The method is currently limited to the calculation of attached flows. The features of the current method can be summarized as follows. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is now needed. Different time steps are used in the different governing equations to improve the convergence speed of the viscous calculations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes. A special discretization technique is also used to stabilize these calculations. A special formulation of the energy equation is used to provide improved transient behavior of solutions which use the full energy equation. The method is then compared with a wide variety of test cases. The freestream Mach numbers range from 0.075 to 2.8 in the calculations. Transonic viscous flow in a converging diverging nozzle is calculated with the method; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good. Essentially incompressible turbulent boundary layer flow in a adverse pressure gradient is calculated and the computed distribution of mean velocity and shear stress are in good agreement with the measurements. At the other end of the Mach number range, a flat plate turbulent boundary layer with a freestream Mach number of 2.8 is calculated using the full energy equation; the computed total temperature distribution and recovery factor agree well with the measurements when a variable Prandtl number is used through the boundary layer.

Nicholson, S.↗
NASA NTRS1989

Numerical simulation and comparison of symmetrical/supercritical airfoils for the near tip region of a helicopter in forward flight

Aerodynamic loads on a multi-bladed helicopter rotor in forward flight at transonic tip conditions are calculated. The unsteady, three-dimensional, time-accurate compressible Reynolds-averaged thin layer Navier-Stokes equations are solved in a rotating coordinate system on a body-conformed, curvilinear grid of C-H topology. Detailed boundary layer and global numerical comparisons of NACA-0012 symmetrical and CAST7-158 supercritical airfoils are made under identical forward flight conditions. The rotor wake effects are modeled by applying a correction to the geometric angle of attack of the blade. This correction is obtained by computing the local induced downwash velocity with a free wake analysis program. The calculations are performed on the Numerical Aerodynamic Simulation Cray 2 and the VPS32 (a derivative of a Cyber 205 at the Langley Research Center) for a model helicopter rotor in forward flight.

Badavi, F. F.↗
NASA NTRS1990

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three-dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form are iterated in time until a steady solution is reached. Evidence is given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at domain boundaries is proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects are accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. For the first, an investigation on the model behavior in case of multiple boundaries is performed. The flow in a developing S-duct is then solved in the laminar regime at Reynolds number (Re) 790 and in the turbulent regime at Re=40,000 using the Baldwin-Lomax model . The Stanitz elbow is then solved using an inviscid version of the same code at M(sub inlet)=0.4. Grid dependence and convergence rate are investigated showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re=2.5x10(exp 6) is solved with the Baldwin-Lomax and the q-omega models. Both approaches showed satisfactory agreement with experiments, although the q-omega model is slightly more accurate.

Michelassi, V.↗
NASA NTRS1991

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form were iterated in time until a steady solution was reached. Evidence was given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at the domain boundaries was proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects were accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. The flow in a developing S-duct was then solved in the laminar regime in a Reynolds number (Re) of 790 and in the turbulent regime at Re equals 40,000 by using the Baldwin-Lomax model. The Stanitz elbow was then solved by using an invicid version of the same code at M sub inlet equals 0.4. Grid dependence and convergence rate were investigated, showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re equals 2.5 times 10(exp 6) was solved with the Baldwin-Lomax and the q-omega models. Both approaches show satisfactory agreement with experiments, although the q-omega model was slightly more accurate.

Michelassi, V.↗
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 NTRS1999

Active Control of Flow Separation Over an Airfoil

Designing an aircraft without conventional control surfaces is of interest to aerospace community. In this direction, smart actuator devices such as synthetic jets have been proposed to provide aircraft maneuverability instead of control surfaces. In this article, a numerical study is performed to investigate the effects of unsteady suction and blowing on airfoils. The unsteady suction and blowing is introduced at the leading edge of the airfoil in the form of tangential jet. Numerical solutions are obtained using Reynolds-Averaged viscous compressible Navier-Stokes equations. Unsteady suction and blowing is investigated as a means of separation control to obtain lift on airfoils. The effect of blowing coefficients on lift and drag is investigated. The numerical simulations are compared with experiments from the Tel-Aviv University (TAU). These results indicate that unsteady suction and blowing can be used as a means of separation control to generate lift on airfoils.

Ravindran, S. S.↗
NASA NTRS1992

Navier-Stokes computations of cavity aeroacoustics with suppression devices

Effectiveness of two devices to suppress the cavity acoustics was computationally investigated. Two dimensional, computational simulations were performed for the transonic, turbulent flows past a cavity, which was first equipped with a rear face ramp and then with a spoiler. The Reynolds-averaged, unsteady, compressible, full Navier-Stokes equations were solved time accurately by a second-order accurate, implicit, upwind, finite-volume method. The effect of turbulence was included through a Baldwin-Lomax model with modifications for the multiple-wall effects and for the highly vortical flow with a shear layer. The results included instantaneous and time-averaged flow properties, and time-series analyses of the pressure inside the cavity, which compared favorably with the available experimental data. These results were also contrasted with the computed aeroacoustics of the same cavity (length-to-depth ratio of 4.5), but without a device, to demonstrate the suppression effectiveness.

Baysal, Oktay↗
NASA NTRS1986

Unsteady separated boundary layer in a transonic diffuser flow with self-excited oscillations

A numerical investigation of two-dimensional unsteady boundary layer in a transonic diffuser flow with self-excited oscillations and strong flow separation by solving the compressible, Reynolds-averaged, thin-layer Navier-Stokes equations with two-equations turbulence model is described. Three different meshes with constant streamwise mesh distribution and varying vertical mesh distribution were used. Results obtained indicate that a refinement of mesh studied here has minimal effect on the mean boundary layer flow but significantly increases the amplitude of oscillation of all flow variables. Comparisons of unsteady wall pressure, velocity profile, terminal shock, and separation pocket among computations and with experiment are presented.

Hsieh, T.↗
NASA NTRS1986

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗
NASA NTRS1986

Explicit finite-volume time-marching calculations of total temperature distributions in turbulent flow

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. This paper investigates the introduction of a new formulation of the energy equation which gives improved transient behavior as the calculation converges. The effect of variable Prandtl number on the temperature distribution through the boundary layer is also investigated. A turbulent boundary layer in an adverse pressure gradient (M = 0.55) is used to demonstrate the improved transient temperature distribution obtained when the new formulation of the energy equation is used. A flat plate turbulent boundary layer with a supersonic free-stream Mach number of 2.8 is used to investigate the effect of Prandtl number on the distribution of properties through the boundary layer. The computed total temperature distribution and recovery factor agree well with the measurements when a variable Prandtl number is used through the boundary layer.

Nicholson, Stephen↗
NASA NTRS1987

Downstream boundary effects on the frequency of self-excited oscillations in transonic diffuser flows

An investigation of downstream boundary effects on the frequency of self-excited oscillations in two-dimensional, separated transonic diffuser flows has been conducted numerically by solving the compressible, Reynolds-averaged, thin-layer Navier-Stokes equation with a two-equation turbulence model. It was found that the unsteady diffuser flowfields are very sensitive to the location of the downstream boundary. Extension of the diffuser downstream boundary significantly reduces the frequency and amplitude of oscillations for pressure, velocity and shock. Computational results suggest that the mechanism causing the self-excited oscillation changes from viscous convective wave dominated oscillations to inviscid acoustic wave dominated oscillations when the location of downstream boundary varies from 8.66 to 134.7 throat height. The existence of a suction slot in the experimental setup obscures the physical downstream boundary and, therefore, presents a difficulty for quantitative comparisons between computation and experiment.

Hsieh, T.↗
NASA NTRS1987

Numerical simulation of viscous flows over transonic aircraft configurations

Numerical simulation of compressible viscous flow fields is performed for a transonic transport configuration. A single structured grid system is constructed using analytical transformations such as conformal mapping, shearing/twisting/rotating/clustering/stretching transformations. The Reynolds-averaged, thin-layer Navier-Stokes equations are solved on a supercomputer, FACOM VP-400, using the LU-ADI factorization method.

Takanashi, Susumu↗