Engineering topics

Rogers, M. M.

Publications and source records attributed to Rogers, M. M..

At least 19 records

NASA NTRS2002

Waves in Turbulent Stably-Stratified Shear Flow

Two approaches for the identification of internal gravity waves in sheared and unsheared homogeneous stratified turbulence are investigated. First, the phase angle between the vertical velocity and density fluctuations is considered. It is found, however, that a continuous distribution of the phase angle is present in weakly and strongly stratified flow. Second, a projection onto the solution of the linearized inviscid equations of motion of unsheared stratified flow is investigated. It is found that a solution of the fully nonlinear viscous Navier-Stokes equations can be represented by the linearized inviscid solution. The projection yields a decomposition into vertical wave modes and horizontal vortical modes.

Jacobitz, F. G.↗

NASA NTRS2002

Waves in Turbulent Stably Stratified Shear Flow

Two approaches for the identification of internal gravity waves in sheared and unsheared homogeneous stratified turbulence are investigated. First, the phase angle between the vertical velocity and density fluctuations is considered. It was found, however, that a continuous distribution of the phase angle is present in weakly and strongly stratified flow. Second, a projection onto the solution of the linearized inviscid equations of motion of unsheared stratified flow is investigated. It was found that a solution of the fully nonlinear viscous Navier-Stokes equations can be represented by the linearized inviscid solution. The projection yields a decomposition into vertical wave modes and horizontal vortical modes.

Jacobitz, F. G.↗

NASA NTRS1999

Aspects of Turbulent / Non-Turbulent Interfaces

A distinct boundary between turbulent and non-turbulent regions in a fluid of otherwise constant properties is found in many laboratory and engineering turbulent flows, including jets, mixing layers, boundary layers and wakes. Generally, the flow has mean shear in at least one direction within t he turbulent zone, but the non-turbulent zones have no shear (adjacent laminar shear is a different case, e.g. transition in a boundary layer). There may be purely passive differences between the turbulent and non-turbulent zones, e.g. small variations in temperature or scalar concentration, for which turbulent mixing is an important issue. The boundary has several major characteristics of interest for the present study. Firstly, the boundary advances into the non-turbulent fluid, or in other words, nonturbulent fluid is entrained. Secondly, the change in turbulence properties across the boundary is remarkably abrupt; strong turbulent motions come close to the nonturbulent fluid, promoting entrainment. Thirdly, the boundary is irregular with a continually changing convoluted shape, which produces statistical intermittency. Its shape is contorted at all scales of the turbulent motion.

Bisset, D. K.↗

NASA NTRS1999

Comparison of AGE and Spectral Methods for the Simulation of Far-Wakes

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.

Bisset, D. K.↗

NASA NTRS1996

On modeling pressure diffusion in non-homogeneous shear flows

New models are proposed for the 'slow and 'rapid' parts of the pressure diffusive transport based on the examination of DNS databases for plane mixing layers and wakes. The model for the 'slow' part is non-local, but requires the distribution of the triple-velocity correlation as a local source. The latter can be computed accurately for the normal component from standard gradient diffusion models, but such models are inadequate for the cross component. More work is required to remedy this situation.

Demuren, A. O.↗

NASA NTRS1996

A-priori testing of sub-grid models for chemically reacting nonpremixed turbulent shear flows

The beta-assumed-pdf approximation of (Cook & Riley 1994) is tested as a subgrid model for the LES computation of nonpremixed turbulent reacting flows, in the limit of cold infinitely fast chemistry, for two plane turbulent mixing layers with different degrees of intermittency. Excellent results are obtained for the computation of integrals properties such as product mass fraction, and the model is applied to other quantities such as powers of the temperature and the pdf of the scalar itself. Even in these cases the errors are small enough to be useful in practical applications. The analysis is extended to slightly out of equilibrium problems such as the generation of radicals, and formulated in terms of the pdf of the scalar gradients. It is shown that the conditional gradient distribution is universal in a wide range of cases whose limits are established. Within those limits, engineering approximations to the radical concentration are also possible. It is argued that the experiments in this paper are essentially in the limit of infinite Reynolds number.

Jimenez, J.↗

NASA NTRS1995

A similarity hypothesis for the two-point correlation tensor in a temporally evolving plane wake

The analysis demonstrated that the governing equations for the two-point velocity correlation tensor in the temporally evolving wake admit similarity solutions, which include the similarity solutions for the single-point moment as a special case. The resulting equations for the similarity solutions include two constants, beta and Re(sub sigma), that are ratios of three characteristic time scales of processes in the flow: a viscous time scale, a time scale characteristic of the spread rate of the flow, and a characteristic time scale of the mean strain rate. The values of these ratios depend on the initial conditions of the flow and are most likely measures of the coherent structures in the initial conditions. The occurrences of these constants in the governing equations for the similarity solutions indicates that these solutions, in general, will only be the same for two flows if these two constants are equal (and hence the coherent structures in the flows are related). The comparisons between the predictions of the similarity hypothesis and the data presented here and elsewhere indicate that the similarity solutions for the two-point correlation tensors provide a good approximation of the measures of those motions that are not significantly affected by the boundary conditions caused by the finite extent of real flows. Thus, the two-point similarity hypothesis provides a useful tool for both numerical and physical experimentalist that can be used to examine how the finite extent of real flows affect the evolution of the different scales of motion in the flow.

Ewing, D. W.↗

NASA NTRS1994

Modeling the two-point correlation of the vector stream function

A new model for the two-point vector stream function correlation has been developed using tensor invariant arguments and evaluated by the comparison of model predictions with DNS data for incompressible homogeneous turbulent shear flow. This two-point vector stream function model correlation can then be used to calculate the two-point velocity correlation function and other quantities useful in turbulence modeling. The model assumes that the two-point vector stream function correlation can be written in terms of the separation vector and a new tensor function that depends only on the magnitude of the separation vector. The model has a single free model coefficient, which has been chosen by comparison with the DNS data. The relative error of the model predictions of the two-point vector stream function correlation is only a few percent for a broad range of the model coefficient. Predictions of the derivatives of this correlation, which are of interest in turbulence modeling, may not be this accurate.

Oberlack, M.↗

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 NTRS1994

The Effect of Two-Dimensional Forcing on Self-Similar Turbulent Plane Wakes

Three direct numerical simulations of time-evolving turbulent plane wakes with velocity deficit Reynolds numbers of about 2,000 have been simulated using a spectral numerical method with up to 600 x 260 x 160 modes. The initial conditions for the simulations are generated from direct numerical simulations of a turbulent boundary layer (momentum thickness Reynolds number of 670), and varying amounts of additional two- dimensional, forcing. In order to preserve the self-similar flow evolution, the forcing is implemented by multiplying all the two-dimensional modes in the initial condition by a constant factor. In the "natural" case no additional forcing is used; in the "forced" and "heavily forced" cases this factor is 5 and 20, respectively. The wake spreading rate Is increased by factors of 1.7 and 7.1 for the two forced cases. The Reynolds stresses are also increased by a similar or even larger factor. These results indicate that the plane wake is much more sensitive to initial forcing than the plane mixing layer. As in the plane mixing layer, two-dimensional forcing promotes more organized large-scale vortical flow structures and these structures axe sometimes separated by "braid regions" containing streamwise "rib" vortices, unlike in the unforced wake.

Rogers, M. M.↗

NASA NTRS1992

Stochastic modeling of turbulent reacting flows

Direct numerical simulations of a single-step irreversible chemical reaction with non-premixed reactants in forced isotropic turbulence at R(sub lambda) = 63, Da = 4.0, and Sc = 0.7 were made using 128 Fourier modes to obtain joint probability density functions (pdfs) and other statistical information to parameterize and test a Fokker-Planck turbulent mixing model. Preliminary results indicate that the modeled gradient stretching term for an inert scalar is independent of the initial conditions of the scalar field. The conditional pdf of scalar gradient magnitudes is found to be a function of the scalar until the reaction is largely completed. Alignment of concentration gradients with local strain rate and other features of the flow were also investigated.

Fox, R. O.↗

NASA NTRS1991

Mixing transition and the cascade to small scales in a plane mixing layer

Direct numerical simulations of time-developing plane mixing layers have been performed using a variety of initial conditions. Up to two pairings of the dominant spanwise rollers have been simulated. When the flow is sufficiently three-dimensional, a pairing can cause the mixing layer to undergo a transition to small-scale turbulence. This small-scale transition is accompanied by increased mixing of a passive scalar, similar to observations of the mixing transition in experiments. As part of the transition process, thin vortex sheets are generated by vortex stretching and roll up in a two-dimensional mixing layer. This higher-order rollup is part of the cascade to small-scale turbulence. Estimates are obtained for the degree of three-dimensionality required for the pairing to initiate the transition.

Moser, R. D.↗

NASA NTRS1988

Phase decorrelation, streamwise vortices and acoustic radiation in mixing layers

Several direct numerical simulations were performed and analyzed to study various aspects of the early development of mixing layers. Included are the phase jitter of the large-scale eddies, which was studied using a 2-D spatially-evolving mixing layer simulation; the response of a time developing mixing layer to various spanwise disturbances; and the sound radiation from a 2-D compressible time developing mixing layer.

Ho, C. M.↗

NASA NTRS1988

A numerical study of bifurcations in a barotropic shear flow

In the last few years, more and more evidence has emerged suggesting that transition to turbulence may be viewed as a succession of bifurcations to deterministic chaos. Most experimental and numerical observations have been restricted to Rayleigh-Benard convection and Taylor-Couette flow between concentric cylinders. An attempt is made to accurately describe the bifurcation sequence leading to chaos in a 2-D temporal free shear layer on the beta-plane. The beta-plane is a locally Cartesian reduction of the equations describing the dynamicss of a shallow layer of fluid on a rotating spherical planet. It is a valid model for large scale flows of interest in meteorology and oceanography.

Huerre, P.↗

NASA NTRS1988

Dynamics of coherent structures in a plane mixing layer

An incompressible, time developing 3-D mixing layer with idealized initial conditions was simulated numerically. Consistent with the suggestions from experimental measurements, the braid region between the dominant spanwise vortices or rolls develops longitudinal vortices or ribs, which are aligned upstream and downstream of a roll and produce spanwise distortion of the rolls. The process by which this distortion occurs is explained by studying a variety of quantities of dynamic importance (e.g., production of enstrophy, vortex stretching). Other quantities of interest (dissipation, helicity density) are also computed and discussed. The currently available simulation only allows the study of the early evolution (before pairing) of the mixing layer. New simulations in progress will relieve this restriction.

Hussain, Fazle↗

NASA NTRS1988

Sensitivity of mixing layers to three-dimensional forcing

It is well known that turbulent mixing layers are dominated by large scale, fairly coherent structures, and that these structures are related to the stability characteristics of the flow. These facts have led researchers to attempt controlling such flows by selectively forcing certain unstable modes, which can in addition have the effect of suppressing other modes. Much of the work on controlling the mixing layer has relied on forcing 2-D instabilities. The results of forcing 3-D instabilities are addressed. The objectives of the work are twofold: to understand how a mixing layer responds to 3-D perturbations, and to test the validity of an amplitude expansion in predicting the mixing layer development. The amplitude expansion could be very useful in understanding and predicting the 3-D response of the flow to a variety of initial conditions.

Riley, J. J.↗

NASA NTRS1988

Advances in turbulent transport modeling based on direct simulations of turbulence

Recent work using numerical simulations to study the modeling of scalar transport is outlined. A model developed by reference to simulations for homogeneous shear flow is presented. This model predicts cross-gradient transports that can give rise to the illusion of counter-gradient transport. The model is found to work remarkably well when compared to simulations for channel flow, but fails to predict important features of the scalar fluxes in mixing layers. Reasons for this failure are discussed and suggestions are made for future work.

Reynolds, W. C.↗

NASA NTRS1987

Direct numerical simulation of buoyantly driven turbulence

Numerical simulations of homogeneous turbulence subject to buoyant forcing were performed. The presence of a mean temperature gradient combined with a gravitational field results in a forcing term in the momentum equations. The development of the turbulence was studied and compared to the decay of similar fields in the absence of gravity. In the buoyantly driven field, the vorticity is preferentially aligned with the intermediate eigenvector of the strain-rate tensor and the local temperature gradient is more likely to be aligned with the most compressive eigenvector. These relationships are qualitatively similar to those observed in previous shear flow results studied by Ashurst (1987). A tensor diffusivity model for passive scalar transport developed from shear flow results in Rogers, Moin, and Reynolds (1986) also predicts this buoyant scalar transport, indicating that the relationship between the scalar flux and the Reynolds stress is similar in both flows.

Ashurst, William T.↗