SEARCH · Search NASA
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.
Results have been compared from Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) of a temporal mixing layer laden with evaporating drops, to assess the ability of LES to reproduce detailed characteristics of DNS.
This paper presents a computational study of real-gas effects on the mean flow and temporal stability of heptane/nitrogen and oxygen/hydrogen mixing layers at supercritical pressures. These layers consist of two counterflowing free streams of different composition, temperature, and density. As in related prior studies reported in NASA Tech Briefs, the governing conservation equations were the Navier-Stokes equations of compressible flow plus equations for the conservation of total energy and of chemical- species masses. In these equations, the expressions for heat fluxes and chemical-species mass fluxes were derived from fluctuation-dissipation theory and incorporate Soret and Dufour effects. Similarity equations for the streamwise velocity, temperature, and mass fractions were derived as approximations to the governing equations. Similarity profiles showed important real-gas, non-ideal-mixture effects, particularly for temperature, in departing from the error-function profile, which is the similarity solution for incompressible flow. The temperature behavior was attributed to real-gas thermodynamics and variations in Schmidt and Prandtl numbers. Temporal linear inviscid stability analyses were performed using the similarity and error-function profiles as the mean flow. For the similarity profiles, the growth rates were found to be larger and the wavelengths of highest instability shorter, relative to those of the errorfunction profiles and to those obtained from incompressible-flow stability analysis. The range of unstable wavelengths was found to be larger for the similarity profiles than for the error-function profiles
The vortex-scalar element method, a scheme which utilizes vortex elements to discretize the region of high vorticity and scalar elements to represent species or temperature fields, is utilized in the numerical simulations of a two-dimensional reacting mixing layer. Computations are performed for a diffusion flame at high Reynolds and Peclet numbers without resorting to turbulence models. In the nonreacting flow, the mean and fluctuation profiles of a conserved scalar show good agreement with experimental measurements. Results for the reacting flow indicate that for temperature independent kinetics, the chemical reaction begins immediately downstream of the splitter plate where mixing starts. Results for the reacting flow with Arrhenius kinetics show an ignition delay, which depends on reactant temperature, before significant chemical reaction occurs. Harmonic forcing changes the structure of the layer, and concomitantly the rates of mixing and reaction, in accordance with experimental results. Strong stretch within the braids in the nonequilibrium kinetics case causes local flame quenching due to the temperature drop associated with the large convective fluxes.
A methodology for simulation of molecular mixing and the resulting velocity and temperature fields has been developed. The ideas are applied to the flow conditions present in the NASA Lewis Planar Reacting Shear Layer (PRSL) facility, and results compared to experimental data. A gaussian transverse turbulent velocity distribution is used in conjunction with a linearly increasing time scale to describe the mixing of different regions of the flow. Equilibrium reaction calculations are then performed on the mix to arrive at a new species composition and temperature. Velocities are determined through summation of momentum contributions. The analysis indicates a combustion efficiency of the order of 80 percent for the reacting mixing layer, and a turbulent Schmidt number of 2/3. The success of the model is attributed to the simulation of large-scale transport of fluid. The favorable comparison shows that a relatively quick and simple PC calculation is capable of simulating the basic flow structure in the reacting and non-reacting shear layer present in the facility given basic assumptions about turbulence properties.
A methodology for simulation of molecular mixing, and the resulting velocity and temperature fields has been developed. The ideas are applied to the flow conditions present in the NASA Lewis Research Center Planar Reacting Shear Layer (PRSL) facility, and results compared to experimental data. A gaussian transverse turbulent velocity distribution is used in conjunction with a linearly increasing time scale to describe the mixing of different regions of the flow. Equilibrium reaction calculations are then performed on the mix to arrive at a new species composition and temperature. Velocities are determined through summation of momentum contributions. The analysis indicates a combustion efficiency of the order of 80 percent for the reacting mixing layer, and a turbulent Schmidt number of 2/3. The success of the model is attributed to the simulation of large-scale transport of fluid. The favorable comparison shows that a relatively quick and simple PC calculation is capable of simulating the basic flow structure in the reacting and nonreacting shear layer present in the facility given basic assumptions about turbulence properties.
The present study is concerned with the optimization of conditional sampling techniques for eduction of naturally occurring large-scale structures in an axisymmetric mixing layer, and a comparison of the natural structure with that induced by means of controlled excitation. Questions regarding the sensitivity of the educed structure to the excitation amplitude are also explored. Aspects regarding the eduction with the high-speed-side reference signal are considered, taking into account structure identification, the smearing effect and a comparison with the structure under excitation, the effects of the threshold level, and the effect of a trigger-level window. Attention is given to structure characteristics and their dependence on Reynolds number and initial condition.
The structure and behavior of a streamwise rib vortex in a direct numerical simulation of a time-developing three-dimensional incompressible plane mixing layer is examined. Where the rib vortex is being stretched, the vorticity vector is primarily directed in the vortex axial direction and the radial and azimuthal velocity distribution is similar to that of a Burger's vortex. In the region where the vortex stretching is negative, there is a change in the local topology of the vortex. The axial flow is decelerated and a negative azimuthal component of vorticity is induced. These features are characteristic of vortex breakdown. The temporal evolution of the rib vortex is similar to the evolution of an axisymmetric vortex in the early stages of vortex breakdown. The effect of vortex breakdown on other parts of the flow is, however, not as significant as the interaction between the rib vortex and other vortices.
A formulation has been developed which combines the simplicity of an experimentally-derived well-established correclation for describing the development of a mixing layer, and a rigorous approach for the description of the dynamics and evaporation of dense or dilute clusters of drops in large coherent vortices.
Direct numerical simulations are performed of a confined three-dimentional, temporally developing, initially isothermal gas mixing layer with one stream laden with as many as 7.3 x 10***sup5*** evaporating hydrocarbon droplets, at moderate gas temperature and subsonic Mach number.
Coherent Vortex Simulation (CVS) filtering has been applied to Direct Numerical Simulation (DNS) data of forced and unforced time-developing turbulent mixing layers. CVS filtering splits the turbulent flow into two orthogonal parts, one corresponding to coherent vortices and the other to incoherent background flow. We have shown that the coherent vortices can be represented by few wavelet modes and that these modes are sufficient to reproduce the vorticity probability distribution function (PDF) and the energy spectrum over the entire inertial range. The remaining incoherent background flow is homogeneous, has small amplitude, and is uncorrelated. These results are compared with those obtained for the same compression rate using large eddy simulation (LES) filtering. In contrast to the incoherent background flow of CVS filtering, the LES subgrid scales have a much larger amplitude and are correlated, which makes their statistical modeling more difficult.
Several studies have shown that sub-mesoscales (SM 1km horizontal scale) play an important role in mixed layer dynamics. In particular, high resolution simulations have shown that in the case of strong down-front wind, the re-stratification induced by the SM is of the same order of the de-stratification induced by small scale turbulence, as well as of that induced by the Ekman velocity. These studies have further concluded that it has become necessary to include SM in ocean global circulation models (OGCMs), especially those used in climate studies. The goal of our work is to derive and assess an analytic parameterization of the vertical tracer flux under baroclinic instabilities and wind of arbitrary directions and strength. To achieve this goal, we have divided the problem into two parts: first, in this work we derive and assess a parameterization of the SM vertical flux of an arbitrary tracer for ocean codes that resolve mesoscales, M, but not sub-mesoscales, SM. In Part 2, presented elsewhere, we have used the results of this work to derive a parameterization of SM fluxes for ocean codes that do not resolve either M or SM. To carry out the first part of our work, we solve the SM dynamic equations including the non-linear terms for which we employ a closure developed and assessed in previous work. We present a detailed analysis for down-front and up-front winds with the following results: (a) down-front wind (blowing in the direction of the surface geostrophic velocity) is the most favorable condition for generating vigorous SM eddies; the de-stratifying effect of the mean flow and re-stratifying effect of SM almost cancel each other out,
The quantum yield of photosynthesis (mol C/mol photons) was calculated at six depths for the waters of the Marine Light-Mixed Layer (MLML) cruise of May 1991. As there were photosynthetically available radiation (PAR) but no spectral irradiance measurements for the primary production incubations, three ways are presented here for the calculation of the absorbed photons (AP) by phytoplankton for the purpose of calculating phi. The first is based on a simple, nonspectral model; the second is based on a nonlinear regression using measured PAR values with depth; and the third is derived through remote sensing measurements. We show that the results of phi calculated using the nonlinear regreesion method and those using remote sensing are in good agreement with each other, and are consistent with the reported values of other studies. In deep waters, however, the simple nonspectral model may cause quantum yield values much higher than theoretically possible.
Motivated by large eddy simulation (LES) modeling of supercritical turbulent flows, transitional states of databases obtained from direct numerical simulations (DNS) of binary-species supercritical temporal mixing layers were examined to understand the subgrid-scale dissipation, and its variation with filter size. Examination of the DSN-scale domain- averaged dissipation confirms previous findings that, out of the three modes of viscous, temperature and species-mass dissipation, the species-mass dissipation is the main contributor to the total dissipation. The results revealed that the percentage of species-mass by total dissipation is nearly invariant across species systems and initial conditions. This dominance of the species-mass dissipation is due to high-density-gradient magnitude (HDGM) regions populating the flow under the supercritical conditions of the simulations; such regions have also been observed in fully turbulent supercritical flows. The domain average being the result of both the local values and the extent of the HDGM regions, the expectations were that the response to filtering would vary with these flow characteristics. All filtering here is performed in the dissipation range of the Kolmogorov spectrum, at filter sizes from 4 to 16 times the DNS grid spacing. The small-scale (subgrid scale, SGS) dissipation was found by subtracting the filtered-field dissipation from the DNS-field dissipation. In contrast to the DNS dissipation, the SGS dissipation is not necessarily positive; negative values indicate backscatter. Backscatter was shown to be spatially widespread in all modes of dissipation and in the total dissipation (25 to 60 percent of the domain). The maximum magnitude of the negative subgrid- scale dissipation was as much as 17 percent of the maximum positive subgrid- scale dissipation, indicating that, not only is backscatter spatially widespread in these flows, but it is considerable in magnitude and cannot be ignored for the purposes of LES modeling. The Smagorinsky model, for example, is unsuited for modeling SGS fluxes in the LES because it cannot render backscatter. With increased filter size, there is only a modest decrease in the spatial extent of backscatter. The implication is that even at large LES grid spacing, the issue of backscatter and related SGS-flux modeling decisions are unavoidable. As a fraction of the total dissipation, the small-scale dissipation is between 10 and 30 percent of the total dissipation for a filter size that is four times the DNS grid spacing, with all OH cases bunched at 10 percent, and the HN cases spanning 24 30 percent. A scale similarity was found in that the domain-average proportion of each small-scale dissipation mode, with respect to the total small-scale dissipation, is very similar to equivalent results at the DNS scale. With increasing filter size, the proportion of the small-scale dissipation in the dissipation increases substantially, although not quite proportionally. When the filter size increases by four-fold, 52 percent for all OH runs, and 70 percent for HN runs, of the dissipation is contained in the subgrid-scale portion with virtually no dependence on the initial conditions of the DNS. The indications from the dissipation analysis are that modeling efforts in LES of thermodynamically supercritical flows should be focused primarily on mass-flux effects, with temperature and viscous effects being secondary. The analysis also reveals a physical justification for scale-similarity type models, although the suitability of these will need to be confirmed in a posteriori studies.
Results are presented of an experimental study tracing the effects of mild streamwise curvature (delta/R is less than 5 percent) on the 3D structure of a two-stream mixing layer at high Reynolds number (about 27,000). Measurements of the mean and fluctuating velocities were obtained on fine cross-plane grids at 11 streamwise locations using a rotatable cross-wide probe. Well-organized streamwise vorticity was generated in both cases through the braid instability. Although the vortex structures decayed with streamwise distance in both cases, the rate of decay for the unstable case was lower. As a result, the unstable layer exhibited noticeable spanwise variations in the mean velocity and Reynolds stresses. Both cases achieved linear growth, but the rate of growth for the unstable case was higher. The far-field spanwise-averaged peak Reynolds stresses were significantly higher for the destabilized case compared to the stabilized case, which exhibited levels comparable to those of a straight case.
It is well known that the growth rate of instability waves of a two-dimensional free shear layer is reduced greatly at supersonic convective Mach numbers. In previous works, it has been shown that new wave modes exist when the shear layers are bounded by a channel due to the coupling effect between the acoustic wave modes and the motion of the mixing layer. The present work studies the simultaneous propagation of multiple stability waves using numerical simulation. It is shown here that the coexistence of two wave modes in the flow field can lead to an oscillatory growth of disturbance energy with each individual wave mode propagating linearly. This is particularly important when the growth rates of the unstable waves are small. It is also shown here that the propagation of two neutrally stable wave modes can lead to a stationary periodic structure of rms fluctuations. In the numerical simulations presented here the forced wave modes are propagating at same frequency, but with different phase velocities. In order to track the growth of each wave mode as it propagates downstream, a numerical method that can effectively detect and separate the contribution of the individual wave is given. It is demonstrated that by a least square fitting of the disturbance field with eigenfunctions the amplitude of each wave mode can be found. Satisfactory results as compared to linear theory are obtained.
This paper summarizes our progress and experience in the development of a Computational-Fluid-Dynamics code on parallel computers to simulate three-dimensional spatially-developing mixing layers. In this initial study, the three-dimensional time-dependent Euler equations are solved using a finite-volume explicit time-marching algorithm. The code was first programmed in Fortran 77 for sequential computers. The code was then converted for use on parallel computers using the conventional message-passing technique, while we have not been able to compile the code with the present version of HPF compilers.
The construction and development of the multi-component traversing system and associated control hardware and software are presented. A hydrogen bubble/laser sheet flow visualization technique was developed to visually study the characteristics of the mixing layers. With this technique large-scale rollers arising from the Taylor-Gortler instability and its interaction with the primary Kelvin-Helmholtz structures can be studied.
Experimental data on the turbulent structure of free shear layers is examined in incompressible air streams. Both single-point and two-point data are used to examine the large coherent structure of the mixing layer. The use of spatial scales to define this large structure is explored. Application of these scales to modeling of free shear layers to show the similarity of such flows is demonstrated. Comparison of the statistical structure with observations from shadowgraph displays of plane shear layer flows shows good qualitative and quantitative agreement.