Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipation rate”

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 55 records · Page 3

Conditional analysis of temperature and strain rate effects on dissipation structure in turbulent non-premixed jet flames

Here, this work presents results from simultaneous high-resolution temperature and velocity measurements in a series of turbulent non-premixed jet flames. The filtered Rayleigh scattering (FRS)-based temperature measurements demonstrate sufficient signal-to-noise (SNR) and spatial resolution to estimate the smallest scalar length scales and accurately determine dissipation rate fields. A comprehensive set of conditional statistics are used to characterize the small-scale structure, including the dependence of dissipation layer widths on Reynolds number, temperature, and dissipation magnitude. In general, the dissipation layer thickness decrease with increasing Reynolds number and increase with increasing temperature. However, dissipation layer widths show two distinct behaviors with respect to dissipation magnitude. For small dissipation values, increases in magnitude results in broadening of the dissipation layer, while for larger magnitude values of dissipation, the layer widths are thinned, highlighting the complexity of small-scale turbulent mixing. Additionally, measured ratios of the dissipation layer width to the Batchelor length scale are consistent across all Reynolds numbers and agree with previous studies in non-reacting flows. The unique aspect about the current set of measurements is the ability to examine the interaction of dissipation structure with turbulent flow parameters for the first time in turbulent non-premixed flames. Particularly, the strain rate/dissipation relationship is examined and compared to previous studies in non-reacting flows. It is found that the dissipation layers tend to align normal to the principal compressive strain axis and this tendency increases with increasing Reynolds number. For the lowest Reynolds number case, no dependence of the dissipation layer width nor dissipation rate magnitude on strain rate is found. However, for higher Reynolds numbers, a strong dependence of the dissipation layer width and dissipation rate magnitude on the principal compressive strain rate is observed. These results indicate the direct role of the compressive strain rate field on small-scale mixing structure in reacting flows.

42 ENGINEERING↗

A near-wall turbulence model and its application to fully developed turbulent channel and pipe flows

A near wall turbulence model and its incorporation into a multiple-time-scale turbulence model are presented. In the method, the conservation of mass, momentum, and the turbulent kinetic energy equations are integrated up to the wall; and the energy transfer rate and the dissipation rate inside the near wall layer are obtained from algebraic equations. The algebraic equations for the energy transfer rate and the dissipation rate inside the near wall layer were obtained from a k-equation turbulence model and the near wall analysis. A fully developed turbulent channel flow and fully developed turbulent pipe flows were solved using a finite element method to test the predictive capability of the turbulence model. The computational results compared favorably with experimental data. It is also shown that the present turbulence model could resolve the over shoot phenomena of the turbulent kinetic energy and the dissipation rate in the region very close to the wall.

Kim, S.-W.↗

Energy dissipation associated with crack extension in an elastic-plastic material

Crack extension in elastic-plastic material involves energy dissipation through the creation of new crack surfaces and additional yielding around the crack front. An analytical procedure, using a two-dimensional elastic-plastic finite element method, was developed to calculate the energy dissipation components during a quasi-static crack extension. The fracture of an isotropic compact specimen was numerically simulated using the critical crack-tip-opening-displacement (CTOD) growth criterion. Two specimen sizes were analyzed for three values of critical CTOD. Results from the analyses showed that the total energy dissipation rate consisted of three components: the crack separation energy rate, the plastic energy dissipation rate, and the residual strain energy rate. All three energy dissipation components and the total energy dissipation rate initially increased with crack extension and finally reached constant values.

Shivakumar, K. N.↗

Energy dissipation associated with crack extension in an elastic-plastic material

Crack extension in elastic-plastic material involves energy dissipation through the creation of new crack surfaces and additional yielding around the crack front. An analytical procedure, using a two-dimensional elastic-plastic finite element method, was developed to calculate the energy dissipation components during a quasi-static crack extension. The fracture of an isotropic compact specimen was numerically simulated using the critical crack-tip-opening-displacement (CTOD) growth criterion. Two specimen sizes were analyzed for three values of critical CTOD. Results from the analysis showed that the total energy dissipation rate consisted of three components: the crack separation energy rate, the plastic energy dissipation rate, and the residual strain energy rate. All three energy dissipation components and the total energy dissipation rate initially increased with crack extension and finally reached constant values.

Shivakumar, K. N.↗

Analysis of direct numerical simulation data of a Mach 4.5 transitional boundary-layer flow

This paper describes the creation, by temporal direct numerical simulation and the analysis based on the Reynolds stress transport equations, of a high-quality data set that represents the laminar-turbulent transition of a high-speed boundary-layer flow. Following Pruett and Zang (1992), and with the help of algorithmic refinements, the evolution of an axial, Mach 4.5 boundary-layer flow along a hollow cylinder is simulated numerically. Favre-averaged Reynolds stress transport equations are derived in generalized curvilinear coordinates and are then specialized to the cylindrical geometry at hand. Reynolds stresses and various turbulence quantities, such as turbulent kinetic energy and turbulent Mach number, are calculated from the numerical data at various stages of the transition process. The kinetic energy 'budgets' are constructed from the transport equations. Various contributing terms for the evolution of kinetic energy, like the rates of production and dissipation, transport, and diffusion, are presented. The compressible dissipation rate is small in comparison with the solenoidal dissipation rate for all times. The pressure-dilatation term is of the same order of magnitude as the compressible dissipation rate.

Dinavahi, Surya P. G.↗

Near-wall modelling of compressible turbulent flows

Work was carried out to extend the near-wall models formulated for the incompressible Reynolds stress equations to compressible flows. The idea of splitting the compressible dissipation function into a solenoidal part that is not sensitive to changes of compressibility indicators and a compressible part that is directly affected by these changes is adopted. This means that all models involving the dissipation rate could be expressed in terms of the solenoidal dissipation rate and an equation governing its transport could be formulated to close the set of compressible Reynolds stress equations. The near-wall modelling of the dissipation rate equation is investigated and its behavior near a wall is studied in detail using k-epsilon closure. It is found that all existing modelled equations give the wrong behavior for the dissipation rate near a wall. Improvements are suggested and the resultant behavior is found to be in good agreement with near-wall data. Furthermore, the present modified k-epsilon closure is used too calculate a flat plate boundary layer and the results are compared with four existing k-epsilon closures. These comparisons show that all closures tested give essentially the same flow properties, except in a region very close to the wall. In this region, the present k-epsilon closure calculations are in better agreement with measurements and direct simulation data; in particular, the behavior of the dissipation rate.

So, Ronald M. C.↗

Version 2 - Tidal Energy Resource Characterization, Velocity and Turbulence Measurements, Cook Inlet, AK, 2021

Version 2 of MHKDR submission #475, which also contains the metadata for these measurements. Files here are stored in netCDF4 format, which can be read in Python using the "netcdf4" or "xarray" packages or in MATLAB using the "ncread()" function. Changes are as follows: MWM1 (mid-water mooring 1) ADV - Time length of data now ends when bottom track data is lost - Improved filtering and quality control of bottom track measurements - improved motion correction of velocity - Improved Doppler noise subtraction - improved velocity, turbulent kinetic energy (TKE), Reynolds stress, time interval (TI), and dissipation rate estimations MWM1 Down-Looking ADCP - Improved filtering and quality control of bottom track measurements - improved motion correction of velocity - Improved Doppler noise subtraction - improved TI and dissipation rate estimations - Removed TKE calculations - accuracy could not be verified MWM1 Up-Looking ADCP - Improved Doppler noise subtraction - improved TI and dissipation rate estimations - Removed TKE calculations - accuracy could not be verified - Velocity unchanged MWM2 ADV - Removes data after instrument was hit by debris - Improved filtering and quality control of bottom track measurements - improved motion correction of velocity - Improved Doppler noise subtraction - improved velocity, TKE, Reynolds stress, TI, and dissipation rate estimations MWM2 Down-Looking ADCP - Improved filtering and quality control of bottom track measurements - improved motion correction of velocity - Improved Doppler noise subtraction - improved TI and dissipation rate estimations - Removed TKE calculations - accuracy could not be verified MWM2 Up-Looking ADCP - Improved filtering and quality control of bottom track measurements - improved motion correction of velocity - Improved Doppler noise subtraction - improved TI and dissipation rate estimations - Removed TKE calculations - accuracy could not be verified THEOM (bottom lander) ADV - Improved Doppler noise subtraction - improved TKE, Reynolds stress, TI, and dissipation rate estimations - Velocity unchanged THEOM ADCP (up-looking) - Improved Doppler noise subtraction - improved TKE, TI, and dissipation rate estimations - Removed vertical TKE component (w'w' bar) - Velocity unchanged

16 TIDAL AND WAVE POWER↗

Protostellar disks and the primitive solar nebula

The objective is to obtain quantitative information on the turbulent transport of mass, angular momentum, and energy under the conditions that characterize the solar nebula, by direct numerical calculations. These calculations were made possible by research conducted on supercomputers (Cray XMP and Cray 2) by the Ames Computational Fluid Dynamics Branch. Techniques were developed that permitted the accurate representation of turbulent flows over the full range of important eddy sizes. So far, these techniques were applied (and verified) primarily in mundane laboratory situations, but they have a strong potential for astrophysical applications. A sequence of numerical experiments were conducted to evaluate the Reynold's stress tensor, turbulent heat transfer rate, turbulent dissipation rate, and turbulent kinetic energy spectrum, as functions of position, for conditions relevant to the solar nebula. Emphasis is placed on the variation of these properties with appropriate nondimensional quantities, so that relations can be derived that will be useful for disk modeling under a variety of hypotheses and initial conditions.

Cassen, P. M.↗

Propagation and Breaking at High Altitudes of Gravity Waves Excited by Tropospheric Forcing

An anelastic approximation is used with a time-variable coordinate transformation to formulate a two-dimensional numerical model that describes the evolution of gravity waves. The model is solved using a semi-Lagrangian method with monotone (nonoscillatory) interpolation of all advected fields. The time-variable transformation is used to generate disturbances at the lower boundary that approximate the effect of a traveling line of thunderstorms (a squall line) or of flow over a broad topographic obstacle. The vertical propagation and breaking of the gravity wave field (under conditions typical of summer solstice) is illustrated for each of these cases. It is shown that the wave field at high altitudes is dominated by a single horizontal wavelength; which is not always related simply to the horizontal dimension of the source. The morphology of wave breaking depends on the horizontal wavelength; for sufficiently short waves, breaking involves roughly one half of the wavelength. In common with other studies, it is found that the breaking waves undergo "self-acceleration," such that the zonal-mean intrinsic frequency remains approximately constant in spite of large changes in the background wind. It is also shown that many of the features obtained in the calculations can be understood in terms of linear wave theory. In particular, linear theory provides insights into the wavelength of the waves that break at high altitudes, the onset and evolution of breaking. the horizontal extent of the breaking region and its position relative to the forcing, and the minimum and maximum altitudes where breaking occurs. Wave breaking ceases at the altitude where the background dissipation rate (which in our model is a proxy for molecular diffusion) becomes greater than the rate of dissipation due to wave breaking, This altitude, in effect, the model turbopause, is shown to depend on a relatively small number of parameters that characterize the waves and the background state.

Prusa, Joseph M.↗

An optical technique for remote sensing of near-surface turbulence

Measurements of turbulent kinetic energy dissipation rates (TKE) or temperature dissipation rates of the the near-surface boundary layer are needed to understand air-sea exchange processes and rates. The capability to accurately estimate these variables by means of a remote technique is relevant to a number of questions ranging from the air-sea transfer of heat and gas to the fate of pollutants.

turbulence optical instrumentation surface process↗

Progress in turbulence modeling for complex flow fields including effects of compressibility

Two second-order-closure turbulence models were devised that are suitable for predicting properties of complex turbulent flow fields in both incompressible and compressible fluids. One model is of the "two-equation" variety in which closure is accomplished by introducing an eddy viscosity which depends on both a turbulent mixing energy and a dissipation rate per unit energy, that is, a specific dissipation rate. The other model is a "Reynolds stress equation" (RSE) formulation in which all components of the Reynolds stress tensor and turbulent heat-flux vector are computed directly and are scaled by the specific dissipation rate. Computations based on these models are compared with measurements for the following flow fields: (a) low speed, high Reynolds number channel flows with plane strain or uniform shear; (b) equilibrium turbulent boundary layers with and without pressure gradients or effects of compressibility; and (c) flow over a convex surface with and without a pressure gradient.

Wilcox, D. C.↗

Investigation of the dilatational dissipation in compressible homogeneous shear flow

The dilatational dissipation rate within compressible homogeneous turbulent shear flow is studied using data from direct numerical simulations. It is found that the dilatational dissipation rate is mainly associated with large scale acoustic waves. Eddy shocklets are observed; however, they have little contribution to the average dissipation rate. A mechanism for the generation of eddy shocklets is shown to be the focusing of acoustic waves. Turbulence models for the dilatational dissipation rate are compared with data from the simulations. It is found that the formulation of Zeman used by Viegas & Rubesin to calculate a compressible mixing layer agrees well with the numerical simulation results for turbulent Mach numbers less than 0.3. However, it is also found that, for the Mach number range occurring in mixing layers, the model does not accurately represent the theory upon which it is based.

Blaisdell, G. A.↗

Turbulence Modeling Effects on the Prediction of Equilibrium States of Buoyant Shear Flows

The effects of turbulence modeling on the prediction of equilibrium states of turbulent buoyant shear flows were investigated. The velocity field models used include a two-equation closure, a Reynolds-stress closure assuming two different pressure-strain models and three different dissipation rate tensor models. As for the thermal field closure models, two different pressure-scrambling models and nine different temperature variance dissipation rate, Epsilon(0) equations were considered. The emphasis of this paper is focused on the effects of the Epsilon(0)-equation, of the dissipation rate models, of the pressure-strain models and of the pressure-scrambling models on the prediction of the approach to equilibrium turbulence. Equilibrium turbulence is defined by the time rate (if change of the scaled Reynolds stress anisotropic tensor and heat flux vector becoming zero. These conditions lead to the equilibrium state parameters. Calculations show that the Epsilon(0)-equation has a significant effect on the prediction of the approach to equilibrium turbulence. For a particular Epsilon(0)-equation, all velocity closure models considered give an equilibrium state if anisotropic dissipation is accounted for in one form or another in the dissipation rate tensor or in the Epsilon(0)-equation. It is further found that the models considered for the pressure-strain tensor and the pressure-scrambling vector have little or no effect on the prediction of the approach to equilibrium turbulence.

Zhao, C. Y.↗

Computing dissipation for molecular-level turbulence simulations

A major difficulty in the analysis of molecular-level simulations is that macroscopic flow quantities are inherently noisy due to molecular fluctuations. An important example for turbulent flows is the kinetic energy dissipation rate. Traditionally, this quantity is calculated from gradients of the macroscopic velocity field, which exacerbates the noise problem. The inability to accurately compute the dissipation rate makes meaningful comparison of molecular-level and continuum simulation results a serious challenge. Herein, we extend previously developed coarse-graining theories to derive an exact molecular-level expression for the dissipation rate, which would circumvent the need to compute gradients of noisy fields. Although the exact expression cannot feasibly be implemented in Sandia’s direct simulation Monte Carlo (DSMC) code SPARTA, we utilize an approximate “hybrid” approach and compare it to the conventional gradient-based approach for planar Couette flow and the two-dimensional Taylor-Green vortex, demonstrating that the hybrid approach is significantly more accurate. Finally, we explore the possibility of adopting a Lagrangian approach to calculate the energy dissipation rate.

74 ATOMIC AND MOLECULAR PHYSICS↗

Single point modeling of rotating turbulent flows

A model for the effects of rotation on turbulence is proposed and tested. These effects which influence mainly the rate of turbulence decay are modeled in a modified turbulent energy dissipation rate equation that has explicit dependence on the mean rotation rate. An appropriate definition of the rotation rate derived from critical point theory and based on the invariants of the deformation tensor is proposed. The modeled dissipation rate equation is numerically well behaved and can be used in conjunction with any level of turbulence closure. The model is applied to the two-equation kappa-epsilon turbulence model and is used to compute separated flows in a backward-facing step and an axisymmetric swirling coaxial jets into a sudden expansion. In general, the rotation modified dissipation rate model shows some improvements over the standard kappa-epsilon model.

Hadid, A. H.↗

A Self-Contained Mapping Closure Approximation for Scalar Mixing

Scalar turbulence exhibits interplays of coherent structures and random fluctuations over a broad range of spatial and temporal scales. This feature necessitates a probabilistic description of the scalar dynamics, which can be achieved comprehensively by using probability density functions (PDFs). Therefore, the challenge is to obtain the scalar PDFs (Lundgren 1967; Dopazo 1979). Generally, the evolution of a scalar is governed by three dynamical processes: advection, diffusion and reaction. In a PDF approach (Pope 1985), the advection and reaction can be treated exactly but the effect of molecular diffusion has to be modeled. It has been shown (Pope 1985) that the effect of molecular diffusion can be expressed as conditional dissipation rates or conditional diffusions. The currently used models for the conditional dissipation rates and conditional diffusions (Pope 1991) have resisted deduction from the fundamental equations and are unable to yield satisfactory results for the basic test cases of decaying scalars in isotropic turbulence, although they have achieved some success in a variety of individual cases. The recently developed mapping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989; Kraichnan 1990; Klimenko & Pope 2003) provides a deductive method for conditional dissipation rates and conditional di usions, and the models obtained can successfully describe the shape relaxation of the scalar PDF from an initial double delta distribution to a Gaussian one. However, the mapping closure approach is not able to provide the rate at which the scalar evolves. The evolution rate has to be modeled. Therefore, the mapping closure approach is not closed. In this letter, we will address this problem.

He, Guo-Wei↗