SEARCH · Search NASA

Results for “Mixing layer”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Simulation of terrain effects using a mesoscale mixed-layer model

The model discussed here is described in the study by Han et al. (1982) and is based on the shallow-water equations. It is noted that the mixed layer is another name for the earth's turbulent boundary layer when the latter is well stirred vertically, as under typical daytime conditions over land. The model abbreviates the vertical resolution by using only one and one-half layers; in this way, the computer power can be concentrated on the horizontal resolution of topographic effects. It is found that a steady flow pattern evolves over terrain when steady forcing occurs in the absence of surface heating or mixed-layer entrainment and that this simplification is removed when surface heating and entrainment do occur. While these findings are regarded as interesting, it is believed that they may be model dependent. The pressure adjustments upon the mixed layer caused by air movements above the mixed layer may be sufficient to preclude the evolution of a steady state on the mesoscale, even under the most ideal conditions.

Deardorff, J. W.↗

NASA NTRS1985

Coherent structures in a turbulent mixing layer - A comparison between direct numerical simulations and experiments

An eduction scheme has been developed in an attempt to determine the characteristics of large-scale vortical structures in a turbulent mixing layer. This analysis scheme has been applied to a set of experimental data taken in a new, larger mixing layer facility designed to minimize boundary and resonance effects. A similar scheme has been developed to apply to the results of a direct numerical simulation of a temporally growing mixing layer. A comparison of the two approaches shows important similarities in the coherent structures. The numerical simulations indicate that low levels of coherent forcing can dramatically change the evolution of the mixing layer. In the absence of such forcing, the numerical simulations and experiments show a lack of regularity in the transverse position, spacing, amplitude, shape and spanwise coherence of the large-scale vortical structures.

Metcalfe, R. W.↗

NASA NTRS1991

Three-dimensional simulations of large eddies in the compressible mixing layer

Consideration is given to the effect of compressibility on a plane mixing layer that is a prototype free shear layer, amenable to study by numerical simulation and experiment. The full time-dependent compressible Navier-Stokes equations are solved numerically for a temporally evolving mixing layer employing a mixed spectral and high-order finite difference method. Simulations with random initial conditions confirm the prediction of linear stability theory that at high Mach numbers oblique waves grow faster than two-dimensional waves. Simulations are presented of the nonlinear temporal evolution of the most rapidly amplified linear instability waves.

Sandham, N. D.↗

NASA NTRS1987

Comparison of experimental and computational techniques for plane mixing layers

In this paper, results from two experimental and two computational investigations of plane turbulent mixing layers are presented and compared. The experimental results include flow visualization data using the smoke laser technique and mean flow and turbulence measurements obtained with hot X-wires and a two-component LDV. Reasonably good agreement is found among these techniques, at least for turbulence quantities up to second order. Reynolds-averaged computations are successful at capturing the complete evolution of the mixing layer, including wake effects in the near field and approach to self-preservation in the far field. The two-dimensional vortex method shows excellent qualitative and quantitative agreement with measured data for the forced mixing layer. For the unforced layer, the results seem to indicate that a three-dimensional computation may be necessary.

Mehta, R. D.↗

NASA NTRS1995

Linear Stability Regime Analysis of the Compressible Reacting Mixing Layer

Previous investigations have shown that a compressible reacting mixing layer can develop two peaks in the mean density weighted vorticity profile. Linear stability analyses show that at these peaks two distinct 'outer' instability modes appear in addition to the more common central mode, which exists unaccompanied in incompressible nonreacting flows. The present study parametrically analyzes the effects of compressibility, heat release, stoichiometry, and density ratio on the amplification rate and obliquity of each stability mode. The mean profiles used in the spatial stability calculation are generated by self-similar solutions of the compressible boundary layer equations combined with the assumption of infinitely fast chemistry. It is shown that the influence of stoichiometry and density ratio on the peaks of the density weighted vorticity profile determines which modes will dominate. Of particular interest are the conditions where two modes are equally amplified, causing the mixing layer to develop into a 'colayer' structure.

Day, M. J.↗

NASA NTRS1994

Direct Simulation of a Self-Similar Turbulent Mixing Layer

Three direct numerical simulations of incompressible turbulent plane mixing layers have been performed. All the simulations were initialized with the same two velocity fields obtained from a direct numerical simulation of a turbulent boundary layer with a momentum thickness Reynolds number of 300 computed by Spalart (J. Fluid Mech. 187, 61, 1988). In addition to a baseline case with no additional disturbances, two simulations were begun with two-dimensional disturbances of varying strength in addition to the boundary layer turbulence. After a development stage, the baseline case and the case with weaker additional two-dimensional disturbances evolve self-similarly, reaching visual thickness Reynolds numbers of up to 20 000. This self-similar period is characterized by a lack of large-scale organized pairings, a lack of streamwise vortices in the 'braid' regions, and scalar mixing that is characterized by 'marching' Probability Density Functions (PDFs). The case begun with strong additional two-dimensional disturbances only becomes approximately self-similar, but exhibits sustained organized large-scale pairings, clearly defined braid regions with streamwise vortices that span them, and scalar PDFs that are 'nonmarching.' It is also characterized by much more intense vertical velocity fluctuations than the other two cases. The statistics and structures in several experiments involving turbulent mixing layers are in better agreement with those of the simulations that do not exhibit organized pairings.

Rogers, Michael M.↗

NASA NTRS1987

Numerical simulation of excited jet mixing layers

A numerical simulation of unsteady flow in jet mixing layers, both with and without external excitation, has been performed by solving the time-dependent compressible Navier-Stokes equations. Computations were performed on a CRAY X-MP computer using MacCormick's explicit finite difference algorithm. Different excitation methods were investigated and were shown to be very effective in controlling the well organized periodic production, shedding and pairing of large scale vortex structures. It is found that pressure excitation was generally more effective than temperature excitation, and that grid refinement results in substantial improvement in the resolution of unsteady features. The location and orientation, in addition to the frequency, of the excitation source are shown to have a significant influence on the production and interaction of large scale vortex structures in the jet mixing layer.

Scott, J. N.↗

NASA NTRS2003

DNS and LES of a Shear-Free Mixing Layer

The purpose of this work is twofold. First, given the computational resources available today, it is possible to reach, using DNS, higher Reynolds numbers than in Briggs et al.. In the present study, the microscale Reynolds numbers reached in the low- and high-energy homogeneous regions are, respectively, 32 and 69. The results reported earlier can thus be complemented and their robustness in the presence of increased turbulence studied. The second aim of this work is to perform a detailed and documented LES of the shear-free mixing layer. In that respect, the creation of a DNS database at higher Reynolds number is necessary in order to make meaningful LES assessments. From the point of view of LES, the shear-free mixing-layer is interesting since it allows one to test how traditional LES models perform in the presence of an inhomogeneity without having to deal with difficult numerical issues. Indeed, as argued in Briggs et al., it is possible to use a spectral code to study the shear-free mixing layer and one can thus focus on the accuracy of the modelling while avoiding contamination of the results by commutation errors etc. This paper is organized as follows. First we detail the initialization procedure used in the simulation. Since the flow is not statistically stationary, this initialization procedure has a fairly strong influence on the evolution. Although we will focus here on the shear-free mixing layer, the method proposed in the present work can easily be used for other flows with one inhomogeneous direction. The next section of the article is devoted to the description of the DNS. All the relevant parameters are listed and comparison with the Veeravalli & Warhaft experiment is performed. The section on the LES of the shear-free mixing layer follows. A detailed comparison between the filtered DNS data and the LES predictions is presented. It is shown that simple eddy viscosity models perform very well for the present test case, most probably because the flow seems to be almost isotropic in the small-scale range that is not resolved by the LES.

Knaepen, B.↗

NASA NTRS1990

A study of the topology of dissipating motions in direct numerical simulations of time-developing compressible and incompressible mixing layers

A preliminary investigation of the geometry of flow patterns in numerically simulated compressible and incompressible mixing layers was carried out using 3-D critical point methodology. Motions characterized by high rates of kinetic energy dissipation and/or high enstrophy were of particular interest. In the approach the partial derivatives of the velocity field are determined at every point in the flow. These are used to construct the invariants of the velocity gradient tensor and the rate-of-strain tensor (P, Q, R, and P(sub s), Q(sub s), R(sub s) respectively). For incompressible flow the first invariant is zero. For the conditions of the compressible simulation, the first invariant is found to be everywhere small, relative to the second and third invariants, and so in both cases the local topology at a point is mainly determined by the second and third invariants. The data at every grid point is used to construct scatter plots of Q versus R and Q(sub s) versus R(sub s). Most points map to a cluster near the origin in Q-R space. However, fine scale motions, that is motions which are characterized by velocity derivatives which scale with the square root of R(sub delta), tend to map to regions which lie far from the origin. Definite trends are observed for motions characterized by high enstrophy and/or high dissipation. The observed trends suggest that, for these motions, the second and third invariants of the velocity gradient and rate-of-strain tensors are strongly correlated. The second and third invariants of the rate-of-strain tensor are related by K(-Q(sub s))(exp 3/2), which is consistent with the above scaling of velocity derivatives. The quantity K appears to depend on Reynolds number with an upper limit K = 2(the square root of 3)/9 corresponding to locally axisymmetric flow. For both the compressible and incompressible mixing layer, regions corresponding to high rates of dissipation are found to be characterized by comparable magnitudes of R(sub ij)R(sub ij) and S(sub ij)S(sub ij). For the incompressible mixing layer, regions characterized by the highest values of enstrophy are found to have relatively low strain rates.

Chen, J. H.↗

NASA NTRS2004

Part 1 of a Computational Study of a Drop-Laden Mixing Layer

This first of three reports on a computational study of a drop-laden temporal mixing layer presents the results of direct numerical simulations (DNS) of well-resolved flow fields and the derivation of the large-eddy simulation (LES) equations that would govern the larger scales of a turbulent flow field. The mixing layer consisted of two counterflowing gas streams, one of which was initially laden with evaporating liquid drops. The gas phase was composed of two perfect gas species, the carrier gas and the vapor emanating from the drops, and was computed in an Eulerian reference frame, whereas each drop was tracked individually in a Lagrangian manner. The flow perturbations that were initially imposed on the layer caused mixing and eventual transition to turbulence. The DNS database obtained included transitional states for layers with various liquid mass loadings. For the DNS, the gas-phase equations were the compressible Navier-Stokes equations for conservation of momentum and additional conservation equations for total energy and species mass. These equations included source terms representing the effect of the drops on the mass, momentum, and energy of the gas phase. From the DNS equations, the expression for the irreversible entropy production (dissipation) was derived and used to determine the dissipation due to the source terms. The LES equations were derived by spatially filtering the DNS set and the magnitudes of the terms were computed at transitional states, leading to a hierarchy of terms to guide simplification of the LES equations. It was concluded that effort should be devoted to the accurate modeling of both the subgridscale fluxes and the filtered source terms, which were the dominant unclosed terms appearing in the LES equations.

Okong'o, Nora A.↗

NASA NTRS1990

Linear stability of the reacting mixing layer

This work is aimed at understanding the stability of reacting mixing layers. Linear instability is analyzed for a wide variety of mixing layers to study the effects of the heat release. Analytic functions as well as laminar solutions of thin shear layer equations are used as the base flows. Both temporal and spatial developing layers are investigated. The laminar solutions are more consistent than the analytic functions. The stability properties are sensitive to the mean profiles. Including variable transport properties changes the mean profiles considerably. Chemical reaction during the instability is not very important; its primary effect is to change the laminar profiles. New modes are found in the outer part of the layer when the heat release is significant, and become dominant for large heat release. In general, heat release stabilizes the center mode but the outer mode is less affected.

Shin, D. S.↗

NASA NTRS2003

Irreversible Entropy Production in Two-Phase Mixing Layers

This report presents a study of dissipation (irreversible production of entropy) in three-dimensional, temporal mixing layers laden with evaporating liquid drops. The purpose of the study is to examine the effects of evaporating drops on the development of turbulent features in flows. Direct numerical simulations were performed to analyze transitional states of three mixing layers: one without drops, and two that included drops at different initial mass loadings. Without drops, the dissipation is essentially due to viscous effects. It was found that in the presence of drops, the largest contribution to dissipation was made by heating and evaporation of the drops, and that at large length scales, this contribution is positive (signifying that the drops reduce turbulence), while at small scales, this contribution is negative (the drops increase turbulence). The second largest contribution to dissipation was found to be associated with the chemical potential, which leads to an increase in turbulence at large scales and a decrease in turbulence at small scales. The next smaller contribution was found to be that of viscosity. The fact that viscosity effects are only third in order of magnitude in the dissipation is in sharp contrast to the situation for the mixing layer without the drops. The next smaller contribution - that of the drag and momentum of the vapor from the drops - was found to be negative at lower mass loading but to become positive at higher mass loading.

Okongo, Nora↗

NASA NTRS1992

The three-dimensional evolution of a plane mixing layer - The Kelvin-Helmholtz rollup

The hydrodynamic evolution of an incompressible plane mixing layer is addressed to elucidate scalar mixing in free shear flows. A detailed description of the onset of three-dimensionality in a mixing layer before or in the absence of pairing is presented. Various simulations were performed to investigate the sensitivity of these results to variations in initial conditions. These variations included changes in amplitude, wavelength, functional form, and relative phasing of the initial low-wavenumber disturbances. Pierrehumber and Widnall's (1982) translative instability eigenfunctions are found to include rib vortices in the braid region and oppositely signed streamwise vorticity in the roller core. The translative instability is an instability of the late-time oversaturated flow. Three-dimensional perturbation growth similar to that of the translative instability can occur whenever spanwise vorticity is present in the braid region. The nonlinear effects that occur when the initial rib circulation is sufficiently high are discussed.

Rogers, Michael M.↗

NASA NTRS1992

Sensitivity of the structure of untripped mixing layers to small changes in initial conditions

An experimental study was conducted concerning the influence of small changes in initial conditions on the near- and far-field evolution of the three-dimensional structure of a plan mixing layer. A two-stream mixing layer with a velocity ratio of 0.6 was generated with the initial boundary layers on the splitter plate laminar and was nominally two-dimensional. The initial conditions were changed slightly by interchanging the high- and low-speed sides of the wind tunnel, while maintaining the same velocities, and hence velocity ratio. This resulted in small changes in the initial boundary layer properties, and the perturbations present in the boundary layers were interchanged between the high- and low-speed sides for the two cases. The results indicate that, even with this relatively minor change in initial conditions, the near-field regions of the two cases differ significantly. The peak Reynolds stress levels in the near-field differ by up to 100 percent, and this is attributed to a difference in the location of the initial spanwise vortex roll-up. In addition, the positions and shapes of the individual streamwise vortical structures differ for the two cases, although the overall structures differ for the two cases, although the overall qualitative description of these structures is comparable. The subsequent reorganization and decay of the streamwise vortical structures is very similar for the two cases. As a result, in the far field, both mixing layers achieve similar structure, yielding comparable growth rates, Reynolds stress, distribution, and spectral content.

Plesniak, M. W.↗

NASA NTRS1990

An experimental study of secondary vortex structure in mixing layers

This report covers the first eight months of an experimental research project on the secondary vortex structure in plane mixing layers. The aim of the project is to obtain quantitative data on the behavior of the secondary structure in a turbulent mixing layer at reasonable reynolds numbers (Re(sub delta(sub w)) approx. 50,000). In particular, we hope to resolve the questions of how the scale of the secondary vortex structure changes with the scale of the mixing layer, and whether the structures are fixed in space, or whether they 'meander' in the spanwise direction.

Bell, J. H.↗

NASA NTRS2001

Development of a Hybrid RANS/LES Method for Turbulent Mixing Layers

Significant research has been underway for several years in NASA Glenn Research Center's nozzle branch to develop advanced computational methods for simulating turbulent flows in exhaust nozzles. The primary efforts of this research have concentrated on improving our ability to calculate the turbulent mixing layers that dominate flows both in the exhaust systems of modern-day aircraft and in those of hypersonic vehicles under development. As part of these efforts, a hybrid numerical method was recently developed to simulate such turbulent mixing layers. The method developed here is intended for configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. Interest in Large Eddy Simulation (LES) methods have increased in recent years, but applying an LES method to calculate the wide range of turbulent scales from small eddies in the wall-bounded regions to large eddies in the mixing region is not yet possible with current computers. As a result, the hybrid method developed here uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall-bounded regions entering a mixing section and uses a LES procedure to calculate the mixing-dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. With this technique, closure for the RANS equations is obtained by using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The LES equations are closed using the Smagorinsky subgrid scale model. Although the function of the Cebeci-Smith model to replace all of the turbulent stresses is quite different from that of the Smagorinsky subgrid model, which only replaces the small subgrid turbulent stresses, both are eddy viscosity models and both are derived at least in part from mixing-length theory. The similar formulation of these two models enables the RANS and LES equations to be solved with a single solution scheme and computational grid. The hybrid RANS-LES method has been applied to a benchmark compressible mixing layer experiment in which two isolated supersonic streams, separated by a splitter plate, provide the flows to a constant-area mixing section. Although the configuration is largely two dimensional in nature, three-dimensional calculations were found to be necessary to enable disturbances to develop in three spatial directions and to transition to turbulence. The flow in the initial part of the mixing section consists of a periodic vortex shedding downstream of the splitter plate trailing edge. This organized vortex shedding then rapidly transitions to a turbulent structure, which is very similar to the flow development observed in the experiments. Although the qualitative nature of the large-scale turbulent development in the entire mixing section is captured well by the LES part of the current hybrid method, further efforts are planned to directly calculate a greater portion of the turbulence spectrum and to limit the subgrid scale modeling to only the very small scales. This will be accomplished by the use of higher accuracy solution schemes and more powerful computers, measured both in speed and memory capabilities.

Georgiadis, Nicholas J.↗

NASA NTRS1996

The structure of variable property, compressible mixing layers in binary gas mixtures

We present the results of a study of the structure of a parallel compressible mixing layer in a binary mixture of gases. The gases included in this study are hydrogen (H2), helium (He), nitrogen (N2), oxygen (02), neon (Ne) and argon (Ar). Profiles of the variation of the Lewis and Prandtl numbers across the mixing layer for all thirty combinations of gases are given. It is shown that the Lewis number can vary by as much as a factor of eight and the Prandtl number by a factor of two across the mixing layer. Thus assuming constant values for the Lewis and Prandtl numbers of a binary gas mixture in the shear layer, as is done in many theoretical studies, is a poor approximation. We also present profiles of the velocity, mass fraction, temperature and density for representative binary gas mixtures at zero and supersonic Mach numbers. We show that the shape of these profiles is strongly dependent on which gases are in the mixture as well as on whether the denser gas is in the fast stream or the slow stream.

Kozusko, F.↗

DOE OSTI2023

On the evolution of scalar iso-surface area density in a turbulent mixing layer

In this paper, the transport equation for scalar iso-surface area density (Σ) in a turbulent, temporally developing mixing layer is examined. Exploring the spatial and temporal evolution of the terms in the Σ transport equation is vital to improving our understanding of turbulent flows characterized by distinct interfaces, e.g. the flame surface or the turbulent/non-turbulent interface. Previous work reported by the authors identified that Σ exhibits self-similar behaviour consistent with the development of the temporal mixing layer. Accordingly, each of the terms in the Σ transport equation is found to behave in a self-similar manner, though there are notable differences in the self-similar behaviours for each term. Based on the results presented herein, it is suggested that the rate of change of Σ (∂Σ/∂t) and the advection term scale with hλ Φ /ΔU, where h is the width of the mixing layer, λ Φ is the scalar Taylor length scale and ΔU is the velocity difference. The production and destruction terms are found to scale with an additional factor (ReSc) 1/2 . In contrast, the molecular diffusion term is found to scale with a factor (ReSc) –1/2 compared to ∂Σ/∂t. Importantly, it is found that the difference between the production and destruction terms, or net surface ‘stretch’, scales with the same factor as ∂Σ/∂t and the advection term, which may have a significant impact on how the evolution of Σ is understood and modelled in turbulent flows.

42 ENGINEERING↗