Engineering Papers⌕ Search

Engineering topics

Zeman, O.

Publications and source records attributed to Zeman, O..

Effects of turbulence compressibility and unsteadiness in compression corner flow

The structure of the separated flow region over a 20 degree compression corner at a free-stream Mach number of 2.84 is investigated computationally using a Reynolds averaged Navier Stokes (R.A.N.S.) solver and kappa-epsilon model. At this Mach number and ramp angle, a steady-state recirculation region of order delta(sub o) is observed, with onset of a 'plateau' in the wall pressure distribution near the corner. At lower ramp angles, separation is negligible, while at an angle of 24 degrees, separation regions of length 2 delta(sub o) are expected. Of interest here is the response of the mathematical model to inclusion of the pressure dilatation term for turbulent kinetic energy. Compared with the experimental data of Smits and Muck (1987), steady-state computations show improvement when the pressure dilatation term is included. Unsteady computations, using both unforced and then forced inlet conditions, did not predict the oscillatory motion of the separation bubble as observed in laboratory experiments. An analysis of the separation bubble oscillation and the turbulent boundary layer (T.B.L.) frequencies for this flow suggests that the bubble oscillations are of nearly the same order as the turbulent frequencies, and therefore difficult for the model to separate and resolve.

Brankovic, A.↗

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

Second-order closure PBL model with new third-order moments: Comparison with LES data

This paper contains two parts. In the first part, a new set of diagnostic equations is derived for the third-order moments for a buoyancy-driven flow, by exact inversion of the prognostic equations for the third-order moment equations in the stationary case. The third-order moments exhibit a universal structure: they all are a linear combination of the derivatives of all the second-order moments, bar-w(exp 2), bar-w theta, bar-theta(exp 2), and bar-q(exp 2). Each term of the sum contains a turbulent diffusivity D(sub t), which also exhibits a universal structure of the form D(sub t) = a nu(sub t) + b bar-w theta. Since the sign of the convective flux changes depending on stable or unstable stratification, D(sub t) varies according to the type of stratification. Here nu(sub t) approximately equal to wl (l is a mixing length and w is an rms velocity) represents the 'mechanical' part, while the 'buoyancy' part is represented by the convective flux bar-w theta. The quantities a and b are functions of the variable N(sub tau)(exp 2), where N(exp 2) = g alpha derivative of Theta with respect to z and tau is the turbulence time scale. The new expressions for the third-order moments generalize those of Zeman and Lumley, which were subsequently adopted by Sun and Ogura, Chen and Cotton, and Finger and Schmidt in their treatments of the convective boundary layer. In the second part, the new expressions for the third-order moments are used to solve the ensemble average equations describing a purely convective boundary laye r heated from below at a constant rate. The computed second- and third-order moments are then compared with the corresponding Large Eddy Simulation (LES) results, most of which are obtained by running a new LES code, and part of which are taken from published results. The ensemble average results compare favorably with the LES data.

Canuto, V. M.↗

Toward modeling wingtip vortices

Wingtip vortices are generated by lifting airfoils; their salient features are compactness and relatively slow rate of decay. The principal motivation for studying the far field evolution of wingtip vortices is the need to understand and predict the extent of the vortex influence during aircraft take-off or landing. On submarines a wingtip vortex ingested into a propeller can be a source of undesirable noise. The main objectives of this research are (1) to establish theoretical understanding of the principal mechanisms that govern the later (diffusive) stages of a turbulent vortex, (2) to develop a turbulence closure model representing the basic physical mechanisms that control the vortex diffusive stage, and further (3) to investigate coupling between the near and far field evolutions; in other words, to study the effect of initial conditions on the vortex lifetime and the ultimate state.

Zeman, O.↗

Investigation of the dilatational dissipation in compressible homogeneous shear flow

The dilatational dissipation rate within compressible homogeneous turbulent shear flow is studied using data from direct numerical simulations. It is found that the dilatational dissipation rate is mainly associated with large scale acoustic waves. Eddy shocklets are observed; however, they have little contribution to the average dissipation rate. A mechanism for the generation of eddy shocklets is shown to be the focusing of acoustic waves. Turbulence models for the dilatational dissipation rate are compared with data from the simulations. It is found that the formulation of Zeman used by Viegas & Rubesin to calculate a compressible mixing layer agrees well with the numerical simulation results for turbulent Mach numbers less than 0.3. However, it is also found that, for the Mach number range occurring in mixing layers, the model does not accurately represent the theory upon which it is based.

Blaisdell, G. A.↗

Response of a supersonic boundary layer to a compression corner

On the basis of direct numerical simulations of rapidly compressed turbulence, Zeman and Coleman have developed a model to represent rapid directional compression contribution to the pressure dilatation term in the turbulent kinetic energy equation. The model has been implemented in the CFD code for simulation of supersonic compression corner flow with an extended separated region. The computational results have shown a significant improvement with respect to the baseline solution given by the standard k- epsilon turbulence model which does not contain any compressibility corrections.

Vandromme, D.↗

Modeling of compressible turbulence

The work is directed toward understanding and modeling compressibility effects in turbulent flows. The primary objective is to investigate how compressibility influences the basic turbulence processes such as turbulence production and dissipation, length scale modification, spectral energy transfer, etc. The second objective is to develop parameterization schemes and models to incorporate compressibility into one-point closure models and into the subgrid scale models for LES techniques. The ultimate purpose of this research is to develop compressible turbulence models which are capable of handling the hypersonic regime. Specifically, the emphasis will be on the model capability of predicting: (1) heat and momentum transfer in hypersonic boundary layers; and (2) mixing and growth of high Mach number shear layers, jet and wakes.

Zeman, O.↗

A model for dispersion from area sources in convective turbulence

Four independent estimates of the vertical distribution of the eddy coefficient for dispersion of a passive contaminant from an extensive area source in a convective layer have been presented. The estimates were based on the following methods: (1) a second-order closure prediction, (2) field data of pollutant concentrations over Los Angeles, (3) lab measurements of particle dispersion, and (4) assumption of equality between momentum and mass transfer coefficients in the free convective limit. It is suggested that K-values estimated both from second-order closure theory and from Los Angeles measurements are systematically underestimated.

Crane, G.↗