Engineering Papers⌕ Search

Engineering topics

Mansour, N. N.

Publications and source records attributed to Mansour, N. N..

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Advancements in turbulence modeling

Direct simulation data are used to test newly proposed nonlinear models in the Reynolds stress anisotropy for the rapid pressure-strain terms. It is found that the nonlinear models will perform better than the linear models. In addition to testing the rapid pressure-strain models, a new k-epsilon model is proposed for near-wall flow simulation. The model is proposed for near-wall flow simulation. The model was developed by fitting low-Reynolds number data from a direct simulation of a channel flow. A new wall-correction for the eddy-viscosity as well as wall-corrections for the epsilon-equation are proposed. The model will give better predictions of the turbulent kinetic profile in the near-wall region than existing models.

Mansour, N. N.↗
NASA NTRS1989

Cross correlation and length scales in turbulent flows near surfaces

Two kinds of length scales are used in turbulent flows; 'functional length scales' such as mixing length, dissipation length L(sub epsilon), etc., and 'flow-field length scales' derived from cross correlations of velocity, pressure, etc. in the flow. Some connection between these scales are derived here. We first consider the cross correlation R(sub vv)(y,y(sub 1)) of the normal components u at two heights y, y(sub 1) above a rigid surface, normalized by the velocity y(sub 1) (greater than y). For shear-free boundary layers it is found theoretically, and in field and numerical experiments that R(sub vv) approximately equals y/y(sub 1). For shear layers it is also found that R(sub vv) approximately equals f(y/y(sub 1)) less than or equal to y,y(sub 1). This function f differs slightly between low Reynolds number numerical simulations and field experiments. The lateral structure defined by R(sub vv)(y,r(sub 3); y(sub 1),0) is also self similar and shows that the eddies centered at about y(sub 1) appear to have constant lateral width a(sub 3) above and below y(sub 1), where a(sub 3, sup +) approximately equals 7+1/(1.4dU(sup +)/dy(sup +)), when normalized on u(sub *) and v, where U is the mean velocity. Results for L(sub epsilon, sup -1) from direct numerical simulation are found to compare well with the formula L(sub epsilon, sup -1) = A(sub B)/y + A(sub S)dU/dy/v, for unidirectional and reversing turbulent boundary layers and channel flow, except near where dU/dy approximately equals 0. The conclusion is that the large-scale eddy structure and length scales in these flows are determined by a combination of shear and blocking, and that the vertical component of turbulence has a self-similar structure in both kinds of boundary layer.

Hunt, J. C. R.↗
NASA NTRS1989

Second-Order Modeling of Low-Reynolds-Number Turbulence Near Walls

This paper presents a set of second-order closure models for low-Reynolds-number turbulence near the wall. Existing closure models for the Reynolds-stress equations were modified to show proper near-wall behavior. A dissipation-rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The equations are solved for the mean velocity, the Reynolds stresses, and the dissipation rate of the turbulent kinetic energy. The calculations are compared with both direct numerical simulations and with measurements. It is shown that the present models perform well in predicting the behavior of the turbulence near a wall. Significant improvements over previous models in predicting the components of the Reynolds stress tensor are obtained in the present models.

Shih, T.-H.↗
NASA NTRS1988

Subharmonic resonance in a mixing layer

The subharmonic resonance phenomenon in a spatially evolving mixing layer is studied using direct simulations of the 2-D Navier-Stokes equations. The computational domain extends to + or - infinity in the cross-stream direction with U(1) = 1.25 and U(2) = 0.25 imposed at + infinity and - infinity respectively. The domain is finite in the streamwise direction with inflow and outflow boundary conditions imposed at x/delta sub omega = 0 and 100, respectively. A hyperbolic-tangent mean velocity profile is assumed at the inlet and the Reynolds number based on the inlet vorticity thickness and velocity difference is Re = 600. It is observed that the phase angle between the fundamental and its subharmonic plays a key role in the spatial development of these modes. Contour plots of vorticity show that varying the phase will have a dramatic effect on the dynamics of the vortices. Pairing or shredding is observed depending on the phase. Fourier decomposition of the time traces show that the fundamental grows, saturates and decays with the downstream distance. The subharmonic has a similar behavior. However, the level at which the modes will saturate is affected by the phase. At 0 deg phase, it was found that as the fundamental saturates, the growth rate of the subharmonic is enhanced. At 90 deg phase, it was found that as the fundamental saturates, the growth rate of the subharmonic is inhibited. In the later case, the growth rate of the subharmonic recovers after saturation of the fundamental. These results are in qualitative agreement with experimental data.

Mansour, N. N.↗
NASA NTRS1988

Reynolds-stress and dissipation-rate budgets in a turbulent channel flow

The budgets for the Reynolds stresses and for the dissipation rate of the turbulence kinetic energy are computed using direct simulation data of a turbulent channel flow. The budget data reveal that all the terms in the budget become important close to the wall. For inhomogeneous pressure boundary conditions, the pressure-strain term is split into a return term, a rapid term, and a Stokes term. The Stokes term is important close to the wall. The rapid and return terms play different roles depending on the component of the term. A split of the velocity pressure-gradient term into a redistributive term and a diffusion term is proposed, which should be simpler to model. The budget data is used to test existing closure models for the pressure-strain term, the dissipation rate, and the transport rate. In general, further work is needed to improve the models.

Mansour, N. N.↗
NASA NTRS1988

Oscillations of drops in zero gravity with weak viscous effects

Nonlinear oscillations and other motions of large axially symmetric liquid drops in zero gravity are studied numerically by a boundary-integral method. The effect of small viscosity is included in the computations by retaining first-order viscous terms in the normal stress boundary condition. This is accomplished by making use of a partial solution of the boundary-layer equations which describe the weak vortical surface layer. Small viscosity is found to have a relatively large effect on resonant mode coupling phenomena.

Lundgren, T. S.↗
NASA NTRS1987

Reynolds-stress and dissipation rate budgets in a turbulent channel flow

The budgets for the Reynolds stresses and for the dissipation rate of the turbulence kinetic energy are computed using direct simulation data of a turbulent channel flow. The budget data reveal that all the terms in the budget become important close to the wall. For inhomogeneous pressure boundary conditions, the pressure-strain term is split into a return term, a rapid term, and a Stokes term. The Stokes term is important close to the wall. The rapid and return terms play different roles depending on the component of the term. A split of the velocity pressure-gradient term into a redistributive term and a diffusion term is proposed, which should be simpler to model. The budget data is used to test existing closure models for the pressure-strain term, the dissipation rate, and the transport rate. In general, further work is needed to improve the models.

Mansour, N. N.↗
NASA NTRS1987

On local approximations of the pressure-strain term in turbulence models

The results of numerical simulations of turbulent channel flows were used to examine the validity of the local approximation of the pressure-strain term in the Reynolds stress transport equation. Outside of the viscous sublayer the local approximation compares very well with the exact pressure strain. This agreement is due, at least in part, to the high correlation between the rapid pressure and its Laplacian, which suggests that only the near parts of the flow contribute to the rapid pressure at a point. In the viscous sublayer the distance over which the mean shear can be considered constant is comparable to the length scale in the normal direction of the correlations of velocity gradients, leading to failure of the local approximation.

Bradshaw, P.↗
NASA NTRS1987

A general form for the dissipation length scale in turbulent shear flows

It has been found that, for a wide range of turbulent wall-bounded shear flows with mean velocity profile U(y), the scale length L-epsilon determining the dissipation is approximately described in terms of distance from the wall (y), the mean shear (dU/dy), and the variance of the normal component of turbulence (bar u-squared), by the formula: L-epsilon to the -1 power is approximately equal to (A sub B / y) + A sub S (dU/dy / square root of bar u-squared) where L-epsilon = epsilon/(bar u-squared) to the 3/2 power. To match with shear-free boundary layers, A sub B is approximately 0.27, and to match with the log layer, A sub S is approximately 0.46. The shear flows tested here were: boundary layers over a flat plate, sink flow, oscillatory flow, and channel flow. The use of the square root of bar u-squared as a velocity scale minimizes the effects of Reynolds number. However, the formula fails within a distance of order L-epsilon for the regions where dU/dy = 0.

Hunt, J. C. R.↗
NASA NTRS1987

Reynolds stress models of homogeneous turbulence

Existing and new models for the rapid and the return terms in the Reynolds stress equations were tested in two ways. One, by direct comparison of the model with simulation data. The other, by simulating the flows using the models and comparing the predicted Reynolds stresses with the data. It was found that existing linear models can be improved and that nonlinear models are in better agreement with the simulation data for a wide variety of flows.

Shih, T.-H.↗
NASA NTRS1987

The decay of isotropic turbulence in a rapidly rotating frame

A direct numerical simulation of the decay of initially isotropic turbulence in a rapidly rotating frame was conducted. This 128 x 128 x 128 simulation was completed for a Reynolds number Re sub lambda = 15.3 and a Rossby number Ro sub lambda = 0.07 based on the initial turbulent kinetic energy and Taylor microscale. The numerical results indicate that the turbulence remains essentially isotropic during the major part of the decay (i.e., beyond the point where the turbulent kinetic energy has decayed to less than 10 percent of its initial value). The rapid rotation has the primary effect of shutting off the energy transfer so that the turbulence dissipation (and hence the rate of decay of the turbulent kinetic energy) is substantially reduced. Consequently, the anisotropy tensor remains essentially unchanged while the energy spectrum undergoes a nearly linear viscous decay (the same results that are predicted by Rapid Distortion Theory which is only formally valid for much shorter elapsed times. Surprisingly, no Taylor-Proudman reorganization of the flow to a two-dimensional state is observed. The implications that these results have on turbulence modeling are discussed briefly along with prospective future research.

Speziale, C. G.↗
NASA NTRS1987

Near-wall k-epsilon turbulence modeling

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

Mansour, N. N.↗
NASA NTRS1987

Passive scalar entrainment and mixing in a forced, spatially-developing mixing layer

Numerical simulations are performed for the forced, spatially-developing plane mixing layer in two and three dimensions. Transport of a passive scalar field is included in the computation. This, together with the allowance for spatial development in the simulations, affords the opportunity for study of the asymmetric entrainment of irrotational fluid into the layer. The inclusion of a passive scalar field provides a means for simulating the effect of this entrainment asymmetry on the generation of 'products' from a 'fast' chemical reaction. Further, the three-dimensional simulations provide useful insight into the effect of streamwise structures on these entrainment and 'fast' reaction processes. Results from a two-dimensional simulation indicate 1.22 parts high-speed fluid are entrained for every one part low-speed fluid. Inclusion of streamwise vortices at the inlet plane of a three-dimensional simulation indicate a further increase in asymmetric entrainment - 1.44:1. Results from a final three-dimensional simulation are presented. In this case, a random velocity perturbation is imposed at the inlet plane. The results indicate the 'natural' development of the large spanwise structures characteristic of the mixing layer.

Lowery, P. S.↗
NASA NTRS1987

Propagation of a premixed flame in a divided-chamber combustor

The propagation of premixed ethylene-air mixtures (of 0.5, 0.525, 0.55, and 0.65 equivalence ratios) in a divided-chamber combustor was investigated. The vessel, divided by a small cylindrical prechamber, had optical access (for laser-schlieren videography) and was instrumented by a pressure transducer. For the Reynolds numbers of 1870, 2300, and 2830, the observed spatial development of the laminar flames showed that the flame position and shape could be scaled by a characteristic time, based on the burned gas flame speed and the length of the prechamber. Above a Reynolds number of 4330, this scaling breaks down the appearance of Kelvin-Helmholtz instabilities. The observed flame propagation was compared with predictions obtained with a numerical model of flame propagation. The calculated spatial and temporal development of the flame in the main combustion chamber agreed with the experimental observations only for the lowest Reynolds number (1870).

Cattolica, R. J.↗
NASA NTRS1987

Near-wall k-epsilon turbulence modeling

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

Mansour, N. N.↗
NASA NTRS1985

Simulation of turbulent mixing layers

A hybrid numerical scheme, combined with a scalar convection scheme, is used to compute turbulent mixing layers. The hybrid scheme is used to solve the vorticity equation by dividing the domain into two regions, applying a Lagrangian method in one region and an Eulerian method in the other. The scalar convection scheme is used to track an interface which separates a dyed fluid from a dye-free fluid, thereby distinguishing the fluid that starts on the high speed side from the fluid that starts on the low speed side. A two-dimensional mixing layer has been simulated using this new method; the roll up of a vortex sheet into structures that interact and pair is observed. The growth of the momentum thickness and the development of the mean turbulence statistics are found to be consistent with experimental observations. Forcing of the layer at moderate and low frequencies shows a dramatic effect on the growth of the momentum thickness and the entrainment ratio. Low frequency forcing enhances the growth of the layer; moderate frequency forcing inhibits the growth.

Mansour, N. N.↗
NASA NTRS1984

Numerical simulation of the tip vortex off a low-aspect-ratio wing at transonic speed

The viscous transonic flow around a low aspect ratio wing was computed by an implicit, three dimensional, thin-layer Navier-Stokes solver. The grid around the geometry of interest is obtained numerically as a solution to a Dirichlet problem for the cube. A low aspect ratio wing with large sweep, twist, taper, and camber is the chosen geometry. The topology chosen to wrap the mesh around the wing with good tip resolution is a C-O type mesh. The flow around the wing was computed for a free stream Mach number of 0.82 at an angle of attack of 5 deg. At this Mach number, an oblique shock forms on the upper surface of the wing, and a tip vortex and three dimensional flow separation off the wind surface are observed. Particle path lines indicate that the three dimensional flow separation on the wing surface is part of the roots of the tip vortex formation. The lifting of the tip vortex before the wing trailing edge is observed by following the trajectory of particles release around the wing tip.

Mansour, N. N.↗