Search NASA⌕ Search

SEARCH · Search NASA

Results for “EDDY DIFFUSION”

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 361 records · Page 20

Plasma transport in the Io torus - The importance of microscopic diffusion

This paper considers the question of whether the distribution of mass in the Io plasma torus is consistent with the concept of interchange eddy transport. Specifically, the flux tube content exhibits a gradual decrease with increasing radial distance from the source near Io without any evidence for substantial density irregularity associated with the plasma source or loss. Using a simple one-dimensional numerical model to simulate macroscopic interchange eddy transport, it is demonstrated that this smooth equilibrium distribution of mass can occur but only with the inclusion of a minimal level of small scale microscopic mixing at a rate approaching Bohm diffusion. Otherwise, the system exhibits a chaotic appearance which never approaches an equilibrium distribution. Various physical mechanisms for the microscopic diffusion process which is required to provide a sufficiently rapid mixing of material between the macroscopic eddies are discussed.

Mei, YI↗

Three-dimensional computation of mixing of transverse injector in a ducted supersonic airstream

Numerical solutions of the three-dimensional mass-averaged elliptic Navier-Stokes equations, including species transport, are obtained for nonreacting, turbulent, mixing flow fields for the case of transverse sonic injection of a secondary gas into a supersonic airstream through a circular orifice injector. Results are presented for flow through a constant area duct and through a duct with a rearward-facing step upstream of the injector. The equations are numerically integrated using MacCormack's explicit method and turbulence is included using the Baldwin-Lomax algebraic eddy viscosity model. In the species transport and energy equations, diffusion coefficients based on Fick's law and an assumption of unit Lewis number are applied. The computations were performed on a CDC-VPS-32 (extended version of Cyber-205) using a grid consisting of approximately 200,000 points. The computed results are compared with experimentally observed penetration and spreading boundaries for an injected gas at two dynamic pressure ratios. Three-dimensional flow field structures are dipicted in terms of static pressure, mass fractions of species and velocity vectors.

Uenishi, K.↗

Turbulent solutions of the equations of fluid motion

Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. Initial three-dimensional cosine velocity fluctuations and periodic boundary conditions are used in most of the work considered. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations, the initially nonrandom flow develops into an apparently random turbulence. Thus randomness or turbulence can arise as a consequence of the structure of the Navier-Stokes equations. The cases considered include turbulence which is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A mean shear is present in some cases. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components, and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Deissler, R. G.↗

Turbulent solutions of equations of fluid motion

Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence, such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations the initially nonrandom flow develops into an apparently random turbulence. The cases considered include turbulence that is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Deissler, R. G.↗

Scalar rate correlation at a turbulent liquid free surface - A two-regime correlation for high Schmidt numbers

An experimental correlation is derived for gas absorption at a turbulent, shear-free liquid interface. The correlation is expressed in terms of the liquid-side turbulence intensity, liquid-side macroscale, and the properties of the diffusing gas and solvent. The transfer coefficient increases linearly with rms velocity up to a point where the eddy Reynolds number reaches a critical (Schmidt number dependent) value. At higher velocities, there is a more rapid linear rise. The slope of the lower Reynolds number region is proportional to the square root of the diffusivity; at Reynolds numbers much higher than that of the break point, the slope becomes independent of diffusivity.

Khoo, Boo-Cheong↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s), and passive scalars K(sub c). First, we took K(sub c) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub c) = K(sub h). Indeed, the turbulence model yelds a salt diffusivity that depends on R(sub i) and R(sub rho) = ((beta)(delta)S/(alpha)(delta)T/(delta)z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in General Circulation Models (GCM'S). The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s). and passive scalars K(sub c). First, we took K(sub s) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derived using the latest model, the KPP model by the NCAR group. Presently, we are running the same Geophysical Fluid Dynamics Laboratory (GFDL) code relaxing the assumptive K(sub s) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R rho (= Beta delta S/delta z/alpha delta T/delta z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of FLuids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusivities K for moment K(sub m), temperature K(sub h), salt K(sub g) and passive scalars K(sub c). First, we took K(sub g) = K(sub h) and tested the model using the Geophysical Fluid Dynamics Laboratory (GFDL) ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub g) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R Rho (= Beta Sigma S/Sigma z/Alpha Sigma T/Sigma z) in such a way that K(sub g) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, Vittorio↗

Experimental investigation of flow through an asymmetric plane diffuser

There is a need for experimental measurements in complex turbulent flows that originate from very well-defined initial conditions. Testing of large-eddy simulations and other higher-order computation schemes requires inlet boundary condition data that are not normally measured. The use of fully developed upstream conditions offers a solution to this dilemma so that the upstream conditions can be adequately computed at any level of sophistication. The plane diffuser experiment by Obi et al. (1993) has received a lot of attention because it has fully-developed inlet conditions and it includes separation from a smooth wall, subsequent reattachment and redevelopment of the downstream boundary layer. The objective of this study is to provide careful qualification and detailed measurements in a recreation of the Obi experiment. The work will include extensive documentation of the flow two-dimensionality and detailed measurements required for testing of flow computations.

Buice, Carl U.↗

An experimental investigation of the separating/reattaching flow over a backstep

Reynolds number effects on the evolution of the turbulence structure in the recovery region of a turbulent separated flow generated by a backward-facing step were studied experimentally. Five different flow conditions and four Reynolds numbers were investigated. All single point measurements were conducted on the plane of symmetry of the wind tunnel starting at approximately x/h = 10 for each flow condition. The upper limit of the measuring domain depends upon the physical size of the tunnel and ranged from 50 to 114 step heights. It is apparent that most of the boundary layer recovery was captured by the experiment. Normal stresses, shear stress, mixing length scales and eddy viscosity are shown, while the higher order velocity products related to the diffusion of the turbulence energy are planned for latter stages of the data analysis.

Jovic, Srboljub↗

An Evaluation of Size-Resolved Cloud Microphysics Scheme Numerics for use with Radar Observations Part II: Condensation and Evaporation

Accurate numerical modeling of clouds and precipitation is essential for weather forecasting and climate change research. While size-resolved (bin) cloud microphysics models predict particle size distributions without imposing shapes, results are subject to artificial size distribution broadening owing to numerical diffusion associated with various processes. Whereas Part 1 addressed collision-coalescence, here we investigate numerical diffusion that occurs in solving condensation and evaporation. In a parcel model framework, all of the numerical schemes examined converge to one solution of condensation and evaporation as the mass grid is refined, and the advection-based schemes are recommended over the reassigning schemes. Including Eulerian vertical advection in a column limits the convergence to some extent, but that limitation occurs at a sufficiently fine mass grid, and the number of iterations in solving vertical advection should be minimized to reduce numerical diffusion. Insubstantial numerical diffusion in solving condensation can be amplified if collision-coalescence is also active, which in turn can be substantially diminished if turbulence effects on collision are incorporated. Large-eddy simulations of a drizzling stratocumulus field reveal that changes in moments of Doppler spectra obtained using different mass grids are consistent with those obtained from the simplified framework, and that spectral moments obtained using a mass grid designed to effectively reduce numerical diffusion are generally closer to observations. Notable differences between the simulations and observations still exist, and our results suggest a need to consider whether factors other than numerical diffusion in the fundamental process schemes employed can cause such differences.

Cloud microphysics↗

On the modeling of the scalar correlations necessary to construct a second-order closure description of turbulent reacting flows

A tentative scheme is proposed for constructing second-order closure models of turbulent reacting flows, with particular reference to mixing processes in turbulent diffusion flows. The model is based on the concepts of typical structure and typical eddy, and treats various hypothetical turbulent flame situations, including the time dependence of species mass fractions passing a point in the flame, the most typical occurrence of species mass fractions, and the case of a rather bizarre typical eddy which shows the effects of nonmixedness.

Donaldson, C. DUP.↗

A two-dimensional model with coupled dynamics, radiation, and photochemistry. 1: Simulation of the middle atmosphere

Spatial and seasonal distribution of ozone and other trace gases are simulated using a interactive two-dimensional model of the stratosphere updated to include full gas phase chemistry. The model consists of a primitive equation dynamics module, a full radiative transfer scheme, and a comprehensive gas phase chemistry module. The circulation is derived from heating rates in the stratosphere that are calculated using model-generated ozone. In the troposphere, parameterized heating rates are adopted. The eddy momentum flux divergence in the zonal mean momentum equation is given by the eddy fluxes of potential vorticity. Eddy fluxes of potential vorticity and tracers are parameterized using a set of predetermined diffusion coefficients. The adopted values for K(sub yy), show a hemispherical asymmetry in that the values in the lower stratosphere are consistently smaller in the southern hemisphere. The asymmetry in K(sub yy) and in the parameterization of the tropospheric heating rate results in an asymmetry in the circulation giving rise to unique signatures in the trace gas distributions. The model successfully simulates the observed asymmetry in the column abundance of the springtime ozone maxima between northern and southern hemisphere. Results for other trace gases are in agreement with the gross observed features although specific differences exist.

Ko, Malcolm K. W.↗

Cloud Condensation in Titan's Lower Stratosphere

A 1-D condensation model is developed for the purpose of reproducing ice clouds in Titan's lower stratosphere observed by the Composite Infrared Spectrometer (CIRS) onboard Cassini. Hydrogen cyanide (HCN), cyanoacetylene (HC3N), and ethane (C2H6) vapors are treated as chemically inert gas species that flow from an upper boundary at 500 km to a condensation sink near Titan's tropopause (-45 km). Gas vertical profiles are determined from eddy mixing and a downward flux at the upper boundary. The condensation sink is based upon diffusive growth of the cloud particles and is proportional to the degree of supersaturation in the cloud formation regIOn. Observations of the vapor phase abundances above the condensation levels and the locations and properties of the ice clouds provide constraints on the free parameters in the model. Vapor phase abundances are determined from CIRS mid-IR observations, whereas cloud particle sizes, altitudes, and latitudinal distributions are derived from analyses of CIRS far-IR observations of Titan. Specific cloud constraints include: I) mean particle radii of2-3 J.lm inferred from the V6 506 cm- band of HC3N, 2) latitudinal abundance distributions of condensed nitriles, inferred from a composite emission feature that peaks at 160/cm , and 3) a possible hydrocarbon cloud layer at high latitudes, located near an altitude of 60 km, which peaks between 60 and 80 cm l . Nitrile abundances appear to diminish substantially at high northern latitudes over the time period 2005 to 2010 (northern mid winter to early spring). Use of multiple gas species provides a consistency check on the eddy mixing coefficient profile. The flux at the upper boundary is the net column chemical production from the upper atmosphere and provides a constraint on chemical pathways leading to the production of these compounds. Comparison of the differing lifetimes, vapor phase transport, vapor phase loss rate, and particle sedimentation, sheds light on temporal stability of the clouds.

Romani, Paul N.↗

Three-Dimensional, Space-Dependent Mesoscale Diffusivity: Derivation and Implications

Recently, we presented a parameterization of an arbitrary tracer 3D mesoscale flux that describes both diabatic and adiabatic regimes without using arbitrary tapering functions. However, we did not parameterize the mesoscale diffusivity, which is the subject of this work. A key difference between the present and previous diffusivity parameterizations is that in the latter, the two main ingredients, mesoscale drift velocity and eddy kinetic energy, were not parameterized but determined using present data, which deprives the models of predictive power. Since winds, stratification, etc., are predicted to change in the future, use of these parameterizations to study future climate scenarios becomes questionable. In this work, we parameterize drift velocity and eddy kinetic energy (vertical–horizontal components), which we first assess with data [WOCE, TOPEX/Poseidon (T/P), and North Atlantic Tracer Release Experiment (NATRE)] and then use in a coarse-resolution stand-alone ocean code under Coordinated Ocean-Ice Reference Experiment I (CORE-I) forcing. We present results for the global ocean temperature and salinity, Atlantic overturning circulation, meridional heat transport, and Drake Passage transport, which we compare with several previous studies. The temperature drift is less than that of five of seven previous OGCMs, and the salinity drift is among the smallest in those studies. The predicted winter Antarctic Circumpolar Current mixed layer depths (MLDs) are in good agreement with the data. Predicting the correct MLD is important in climate studies since models that predict very deep mixed layers transfer more of the radiative perturbation to the deep ocean, reducing surface warming (and vice versa).

Eddies; Mesoscale processes; Mixing; Ocean models;↗

Reynolds-Stress Budgets in an Impinging Shock Wave/Boundary-Layer Interaction

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Comparisons with experimental data showed a sensitivity of the current prediction to the modeling of the sidewalls. This was found to be common among various computational studies in the literature where periodic boundary conditions were used in the spanwise direction, as was the case in the present work. Thus, although the experiment was quasi-two-dimensional, the present simulation was determined to be two-dimensional. Quantities present in the exact equation of the Reynolds-stress transport, i.e., production, molecular diffusion, turbulent transport, pressure diffusion, pressure strain, dissipation, and turbulent mass flux were calculated. Reynolds-stress budgets were compared with past large-eddy simulation and direct numerical simulation datasets in the undisturbed portion of the turbulent boundary layer to validate the current approach. The budgets in SBLI showed the growth in the production term for the primary normal stress and energy transfer mechanism was led by the pressure strain term in the secondary normal stresses. The pressure diffusion term, commonly assumed as negligible by turbulence model developers, was shown to be small but non-zero in the normal stress budgets, however it played a key role in the primary shear stress budget.

shock wave/boundary-layer interactio↗

Enhancement of Mixing and Reaction in High-speed Combustor Flowfields

Research has been underway for a number of years, both in the United States and abroad, to develop advanced aerospace propulsion systems for use late in this century and beyond. One program is now underway at the NASA Langley Research Research has been underway for a number of years, both in the United States and abroad, to develop advanced aerospace propulsion systems for use late in this century and beyond. One program is now underway at the NASA Langley Research Center to develop a hydrogen-fueled supersonic combustion ramjet (scramjet) that is capable of propelling a vehicle at hypersonic speeds in the atmosphere. A part of that research has been directed toward the optimization of the scramjet combustor and, in particular, the efficiency of fuel-air mixing and reaction taking place in the engine. In the high-speed vehicle configurations currently being considered, achieving a high combustor efficiency becomes particularly difficult. With increasing combustor Mach number, the degree of fuel-air mixing that can be achieved through natural convective and diffusive processes is reduced leading to an overall decrease in combustion efficiency and thrust. Compressible shear layers and jets provide a good model for studying the physical processes occurring in high speed mixing and combustion in a scramjet. Mixing layers are characterized by large scale eddies that form due to the high shear that is present between the fuel and the air streams. These eddies entrain the fuel and air into the mixing region between the fluids leading to increased surface area and locally steep concentration gradients. Molecular diffusion then occurs across the strained interfaces. In an early study of high-speed mixing, Brown and Roshko show that the spreading rate of a supersonic mixing layer decreases with increasing Mach number, exhibiting a factor of three decrease in spread rate as compared with an incompressible mixing layer with the same density ratio. They conclude that the reduced spread rate is primarily due to compressibility. Papamoschou and Roshko also observe that the spreading rate of compressible mixing layers is significantly reduced over that of incompressible layers. To characterize the structure of the flow quantitatively, they define a convective Mach number. The reduction in mixing layer spreading rate is shown in these experiments to correlate well with increasing convective Mach number.

Drummond, J. Philip↗