Search NASA⌕ Search

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.

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

At least 109 records · Page 6

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

This second of three reports on a computational study of a mixing layer laden with evaporating liquid drops presents the evaluation of Large Eddy Simulation (LES) models. The LES models were evaluated on an existing database that had been generated using Direct Numerical Simulation (DNS). The DNS method and the database are described in the first report of this series, Part 1 of a Computational Study of a Drop-Laden Mixing Layer (NPO-30719), NASA Tech Briefs, Vol. 28, No.7 (July 2004), page 59. The LES equations, which are derived by applying a spatial filter to the DNS set, govern the evolution of the larger scales of the flow and can therefore be solved on a coarser grid. Consistent with the reduction in grid points, the DNS drops would be represented by fewer drops, called computational drops in the LES context. The LES equations contain terms that cannot be directly computed on the coarser grid and that must instead be modeled. Two types of models are necessary: (1) those for the filtered source terms representing the effects of drops on the filtered flow field and (2) those for the sub-grid scale (SGS) fluxes arising from filtering the convective terms in the DNS equations. All of the filtered-sourceterm models that were developed were found to overestimate the filtered source terms. For modeling the SGS fluxes, constant-coefficient Smagorinsky, gradient, and scale-similarity models were assessed and calibrated on the DNS database. The Smagorinsky model correlated poorly with the SGS fluxes, whereas the gradient and scale-similarity models were well correlated with the SGS quantities that they represented.

Okongo, Nora↗

A numerical study of a class of TVD schemes for compressible mixing layers

At high Mach numbers the two-dimensional time-developing mixing layer develops shock waves, positioned around large-scale vortical structures. A suitable numerical method has to be able to capture the inherent instability of the flow, leading to the roll-up of vortices, and also must be able to capture shock waves when they develop. Standard schemes for low speed turbulent flows, for example spectral methods, rely on resolution of all flow-features and cannot handle shock waves, which become too thin at any realistic Reynolds number. The performance of a class of second-order explicit total variation diminishing (TVD) schemes on a compressible mixing layer problem was studied. The basic idea is to capture the physics of the flow correctly, by resolving down to the smallest turbulent length scales, without resorting to turbulence or sub-grid scale modeling, and at the same time capture shock waves without spurious oscillations. The present study indicates that TVD schemes can capture the shocks accurately when they form, but (without resorting to a finer grid) have poor accuracy in computing the vortex growth. The solution accuracy depends on the choice of limiter. However a larger number of grid points are in general required to resolve the correct vortex growth. The low accuracy in computing time-dependent problems containing shock waves as well as vortical structures is partly due to the inherent shock-capturing property of all TVD schemes. In order to capture shock waves without spurious oscillations these schemes reduce to first-order near extrema and indirectly produce clipping phenomena, leading to inaccuracy in the computation of vortex growth. Accurate simulation of unsteady turbulent fluid flows with shock waves will require further development of efficient, uniformly higher than second-order accurate, shock-capturing methods.

Sandham, N. D.↗

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

Aspects of Supercritical Turbulence: Direct Numerical Simulation of O2/H2 and C7H16/N2 Temporal Mixing Layers

Results from Direct Numerical Simulations of temporal, supercritical mixing layers for two species systems are analyzed to elucidate species-specific turbulence aspects. The two species systems, O2/H2 and C7HG16/N2, have different thermodynamic characteristics; thus, although the simulations are performed at similar reduced pressure (ratio of the pressure to the critical pressure), the former system is dose to mixture ideality and has a relatively high solubility with respect to the latter, which exhibits strong departures from mixture ideality Due to the specified, smaller initial density stratification, the C7H16/N2 layers display higher growth and increased global molecular mixing as well as larger turbulence levels. However, smaller density gradients at the transitional state for the O2/H2 system indicate that on a local basis, the layer exhibits an enhanced mixing, this being attributed to the increased solubility and to mixture ideality. These thermodynamic features are shown to affect the irreversible entropy production (i.e. the dissipation), which is larger for the O2/H2 layer and is primarily concentrated in high density-gradient magnitude regions that are distortions of the initial density stratification boundary. In contrast, the regions of largest dissipation in the C7H16/N2 layer are located in high density-gradient magnitude regions resulting from the mixing of the two fluids.

Bellan, J.↗

Effects of Mesoscale Eddies in the Active Mixed Layer: Test of the Parametrisation in Eddy Resolving Simulations

In eddy resolving simulations, we test a mixed layer mesoscale parametrisation, developed recently by Canuto and Dubovikov [Ocean Model., 2011, 39, 200-207]. With no adjustable parameters, the parametrisation yields the horizontal and vertical mesoscale fluxes in terms of coarse-resolution fields and eddy kinetic energy (EKE). We compare terms of the parametrisation diagnosed from coarse-grained fields with the eddy mesoscale fluxes diagnosed directly from the high resolution model. An expression for the EKE in terms of mean fields has also been found to get a closed parametrisation in terms of the mean fields only. In 40 numerical experiments we simulated two types of flows: idealised flows driven by baroclinic instabilities only, and more realistic flows, driven by wind and surface fluxes as well as by inflow-outflow. The diagnosed quasi-instantaneous horizontal and vertical mesoscale buoyancy fluxes (averaged over 1-2 degrees and 10 days) demonstrate a strong scatter typical for turbulent flows, however, the fluxes are positively correlated with the parametrisation with higher (0.5-0.74) correlations at the experiments with larger baroclinic radius Rossby. After being averaged over 3-4 months, diffusivities diagnosed from the eddy resolving simulations are consistent with the parametrisation for a broad range of parameters. Diagnosed vertical mesoscale fluxes restratify mixed layer and are in a good agreement with the parametrisation unless vertical turbulent mixing in the upper layer becomes strong enough in comparison with mesoscale advection. In the latter case, numerical simulations demonstrate that the deviation of the fluxes from the parametrisation is controlled by dimensionless parameter estimating the ratio of vertical turbulent mixing term to mesoscale advection. An analysis using a modified omega-equation reveals that the effects of the vertical mixing of vorticity is responsible for the two-three fold amplification of vertical mesoscale flux. Possible physical mechanisms, responsible for the amplification of vertical mesoscale flux are discussed.

simulation↗

Growth Of Instabilities In Two Types Of Mixing Layers

Report discusses issues involved in computation of growth of instabilities and vorticity in mixing (shear) layer between two parallel flows having different speeds and separated initially by flat surface. Study focusses upon degree to which these growing flow structures hasten or retard their own further growth, and upon related question of differences between flows computed according to mathematical models of S-layer and T-layer mixing. Transformation method used offers suitable tool for study of such flows.

Kaul, Upender K.↗

Transitional mixing in a binary species supercritical, temporal mixing layer

The turbulent mixing of fluids at high pressure is a topic of much interest as it is relavant both to the natural phenomena and to technical applications. In the realm of technical applications, liquid rocket combustion presents a particular challenge as the operating conditions are supercritical with respect to both fuel and oxidizer, making mandatory the understanding of supercritical fluid behavior for a potentially explosive mixture. In this situation, numerical simulations with validated models can contribute information that would be otherwise impossible to obtain.

supercritical↗

Model of Mixing Layer With Multicomponent Evaporating Drops

A mathematical model of a three-dimensional mixing layer laden with evaporating fuel drops composed of many chemical species has been derived. The study is motivated by the fact that typical real petroleum fuels contain hundreds of chemical species. Previously, for the sake of computational efficiency, spray studies were performed using either models based on a single representative species or models based on surrogate fuels of at most 15 species. The present multicomponent model makes it possible to perform more realistic simulations by accounting for hundreds of chemical species in a computationally efficient manner. The model is used to perform Direct Numerical Simulations in continuing studies directed toward understanding the behavior of liquid petroleum fuel sprays. The model includes governing equations formulated in an Eulerian and a Lagrangian reference frame for the gas and the drops, respectively. This representation is consistent with the expected volumetrically small loading of the drops in gas (of the order of 10 3), although the mass loading can be substantial because of the high ratio (of the order of 103) between the densities of liquid and gas. The drops are treated as point sources of mass, momentum, and energy; this representation is consistent with the drop size being smaller than the Kolmogorov scale. Unsteady drag, added-mass effects, Basset history forces, and collisions between the drops are neglected, and the gas is assumed calorically perfect. The model incorporates the concept of continuous thermodynamics, according to which the chemical composition of a fuel is described probabilistically, by use of a distribution function. Distribution functions generally depend on many parameters. However, for mixtures of homologous species, the distribution can be approximated with acceptable accuracy as a sole function of the molecular weight. The mixing layer is initially laden with drops in its lower stream, and the drops are colder than the gas. Drop evaporation leads to a change in the gas-phase composition, which, like the composition of the drops, is described in a probabilistic manner

Bellan, Josette↗

Direct simulations of chemically reacting turbulent mixing layers

The report presents the results of direct numerical simulations of chemically reacting turbulent mixing layers. The work consists of two parts: (1) the development and testing of a spectral numerical computer code that treats the diffusion reaction equations; and (2) the simulation of a series of cases of chemical reactions occurring on mixing layers. The reaction considered is a binary, irreversible reaction with no heat release. The reacting species are nonpremixed. The results of the numerical tests indicate that the high accuracy of the spectral methods observed for rigid body rotation are also obtained when diffusion, reaction, and more complex flows are considered. In the simulations, the effects of vortex rollup and smaller scale turbulence on the overall reaction rates are investigated. The simulation results are found to be in approximate agreement with similarity theory. Comparisons of simulation results with certain modeling hypotheses indicate limitations in these hypotheses. The nondimensional product thickness computed from the simulations is compared with laboratory values and is found to be in reasonable agreement, especially since there are no adjustable constants in the method.

Riley, J. J.↗

A statistical model of turbulence in two-dimensional mixing layers

A statistical model of turbulence in fully developed two-dimensional incompressible turbulent mixing layers is proposed. The development of this model is largely motivated by the recent experimental observations of Brown and Roshko (1974). The model is based on the proposition that the turbulence of a fully developed two-dimensional incompressible mixing layer is in a state of quasi-equilibrium. The model is used to predict the second-order turbulence statistics of the flow including single-point turbulent Reynolds stress distribution, intensity of turbulent velocity components, rms turbulent pressure fluctuations, power spectra and two-point space-time correlation functions. It is shown that numerical results compare favorably with available experimental measurements.

Tam, C. K. W.↗

Effects of streamwise vorticity injection on turbulent mixing layer development

This paper investigates the effect of injected strong streamwise vorticity on the structure and development of a plane mixing layer originating from tripped boundary layers. The experiments were conducted in a mixing layer wind tunnel consisting of two separate legs which were driven individually by centrifugal blowers. It was found that, while the vorticity injection increased the growth rate in the near-field, the asymptotic growth rate was reduced by a factor of about two, together with the peak Reynolds stress levels. The result is attributed to the effect of the relatively strong and short wavelength streamwise vorticity in making the spanwise structures more three-dimensional and slowing down their pairing process, thus reducing entrainment, and hence growth.

Bell, James H.↗

Computation of turbulent high speed mixing layers using a two-equation turbulence model

A two-equation turbulence model was extended to be applicable for compressible flows. A compressibility correction based on modelling the dilational terms in the Reynolds stress equations were included in the model. The model is used in conjunction with the SPARK code for the computation of high speed mixing layers. The observed trend of decreasing growth rate with increasing convective Mach number in compressible mixing layers is well predicted by the model. The predictions agree well with the experimental data and the results from a compressible Reynolds stress model. The present model appears to be well suited for the study of compressible free shear flows. Preliminary results obtained for the reacting mixing layers are included.

Narayan, J. R.↗

Development of a two-stream mixing layer from tripped and untripped boundary layers

The effects of the state of the initial boundary layers on the development of a two-stream, plane mixing layer, with a velocity ratio of 0.6 are experimentally investigated. Spanwise-average profiles are compared for the first time. The results indicate that both the near and far-field growth rates for the untripped case are significantly higher than the tripped case. The maximum Reynolds stresses and higher-order products for the two cases behave very differently in the near-field, but asymptote to approximately the same constant levels far downstream. The mean velocity and turbulence profiles in this region also collapse adequately for the two cases when plotted in similarity coordinates. The distance required to achieve self-similarity is distinctly shorter for the tripped case, in contrast to previous observations. The higher growth rate for the untripped case is attributed to the presence of streamwise vortices which result in additional entrainment by the mixing layer.

Bell, James H.↗

Spanwise averaging of plane mixing layer properties

An attempt is made to demonstrate the effects of streamwise vorticity on such mixing layer global properties as layer thickness and peak Reynolds stresses. Attention is given to differences between results obtained from spanwise averaging over a number of velocity profiles and those given by the more conventional approach of analyzing a single profile measured along the test section centerline. It is found that the growth rate and streamwise development of peak Reynolds stresses in the near field can be significantly affected by spanwise variations in the mixing layer.

Bell, James H.↗

Mixing and non-equilibrium chemical reaction in a compressible mixing layer

The effects of compressibility, chemical reaction exothermicity, and non-equilibrium chemical modeling in a reacting plane mixing layer were investigated by means of two dimensional direct numerical simulations. The chemical reaction was irreversible and second order of the type A + B yields Products + Heat. The general governing fluid equations of a compressible reacting flow field were solved by means of high order finite difference methods. Physical effects were then determined by examining the response of the mixing layer to variation of the relevant non-dimensionalized parameters. The simulations show that increased compressibility generally results in a suppressed mixing, and consequently a reduced chemical reaction conversion rate. Reaction heat release was found to enhance mixing at the initial stages of the layer growth, but had a stabilizing effect at later times. The increased stability manifested itself in the suppression or delay of the formation of large coherent structures within the flow. Calculations were performed for a constant rate chemical kinetics model and an Arrhenius type kinetic prototype. The choice of the model was shown to have an effect on the development of the flow. The Arrhenius model caused a greater temperature increase due to reaction than the constant kinetic model. This had the same effect as increasing the exothermicity of the reaction. Localized flame quenching was also observed when the Zeldovich number was relatively large.

Steinberger, Craig J.↗

Stability of mixing layers

The research program for the first year of this project (see the original research proposal) consists of developing an explicit marching scheme for solving the parabolized stability equations (PSE). Performing mathematical analysis of the computational algorithm including numerical stability analysis and the determination of the proper boundary conditions needed at the boundary of the computation domain are implicit in the task. Before one can solve the parabolized stability equations for high-speed mixing layers, the mean flow must first be found. In the past, instability analysis of high-speed mixing layer has mostly been performed on mean flow profiles calculated by the boundary layer equations. In carrying out this project, it is believed that the boundary layer equations might not give an accurate enough nonparallel, nonlinear mean flow needed for parabolized stability analysis. A more accurate mean flow can, however, be found by solving the parabolized Navier-Stokes equations. The advantage of the parabolized Navier-Stokes equations is that its accuracy is consistent with the PSE method. Furthermore, the method of solution is similar. Hence, the major part of the effort of the work of this year has been devoted to the development of an explicit numerical marching scheme for the solution of the Parabolized Navier-Stokes equation as applied to the high-seed mixing layer problem.

Tam, Christopher↗