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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Lightning activity on Jupiter

Photographic observations of the nightside of Jupiter by the Voyager 1 spacecraft show the presence of extensive lightning activity. Detection of whistlers by the plasma wave analyzer confirms the optical observations and implies that many flashes were not recorded by the Voyager camera because the intensity of the flashes was below the threshold sensitivity of the camera. Measurements of the optical energy radiated per flash indicate that the observed flashes had energies similar to that for terrestrial superbolts. The best estimate of the lightning energy dissipation rate of 0.0004 W/sq m was derived from a consideration of the optical and radiofrequency measurements. The ratio of the energy dissipated by lightning compared to the convective energy flux is estimated to be between 0.000027 and 0.00005. The terrestrial value is 0.0001.

Borucki, W. J.↗
NASA NTRS1988

Tidal dissipation, surface heat flow, and figure of viscoelastic models of Io

The deformation of Io, the tidal dissipation rate, and its interior spatial distribution are investigated by means of numerical simulations based on (1) a three-layer model (with dissipation in the mantle) or (2) a four-layer model (with dissipation in the asthenosphere). The mathematical derivation of the models is outlined; the selection of the input-parameter values is explained; the results are presented in extensive graphs and contour maps; and the constraints imposed on the models by observational data on the hot-spot distribution, tidal deformation, and gravity field are discussed in detail. It is found that both dissipation mechanisms may play a role on Io: model (2) is better able to explain the concentration of hot spots near the equator, while the presence of a large hot spot near the south pole (if confirmed by observations) would favor model (1).

Segatz, M.↗
NASA NTRS1985

An implicit finite-difference code for a two-equation turbulence model for three-dimensional flows

An implicit finite difference code was developed which solves the transport equations for the turbulence kinetic energy and its dissipation rate in generalized coordinates in three dimensions. The finite difference equations are solved using the Beam-Warming algorithm. The kinetic energy-dissipation code, KEM, provides the closure; i.e., the turbulent viscosity for calculation of either compressible or incompressible flows. Turbulent internal flow over a backward-facing step has been calculated using the present code in conjunction with the Incompressible Navier-Stokes Code, INS3D. The results are in good agreement with experiments and two dimensional computations of other researchers.

Kaul, U. K.↗
NASA NTRS1998

Numerical Modeling Studies of Wake Vortex Transport and Evolution Within the Planetary Boundary Layer

In support of the wake vortex effect of the Terminal Area Productivity program, we have put forward four tasks to be accomplished in our proposal. The first task is validation of two-dimensional wake vortex-turbulence interaction. The second task is investigation of three-dimensional interaction between wake vortices and atmospheric boundary layer (ABL) turbulence. The third task is ABL studies. The, fourth task is addition of a Klemp-Durran condition at the top boundary for TASS model. The accomplishment of these tasks will increase our understanding of the dynamics of wake vortex and improve forecasting systems responsible for air safety and efficiency. The first two tasks include following three parts: (a) Determine significant length scale for vortex decay and transport, especially the length scales associated with the onset of Crow instability (Crow, 1970); (b) Study the effects of atmospheric turbulence on the decay of the wake vortices; and (c) Determine the relationships between decay rate, transport properties and atmospheric parameters based on large eddy simulation (LES) results and the observational data. These parameters may include turbulence kinetic energy, dissipation rate, wind shear and atmospheric stratification. The ABL studies cover LES modeling of turbulence structure within planetary boundary layer under transition and stable stratification conditions. Evidences have shown that the turbulence in the stable boundary layer can be highly intermittent and the length scales of eddies are very small compared to those in convective case. We proposed to develop a nesting grid mesh scheme and a modified Klemp-Durran conditions (Klemp and Wilhelmson, 1978) at the top boundary for TASS model to simulate planetary boundary layer under stable stratification conditions. During the past year, our group has made great efforts to carry out the above mentioned four tasks simultaneously. The work accomplished in the last year will be described in the next section.

Lin, Yuh-Lang↗
NASA NTRS1989

Measurements of turbulence and its evolution and variability during MAP

The understanding of turbulence in the middle atmosphere has improved considerably during the MAP period. For a theoretical viewpoint, several advances were made including understanding the ways in which turbulence is generated, and the differences between the rates of diffusion of momentum and heat. Experimentally, a proper understanding of how radars can be used to measure turbulence has emerged, and turbulent energy dissipation rates in the middle atmosphere were measured with MF, HF, and VHF radars. New rocket techniques were developed which have enabled detailed studies of the fine structure of turbulence to be made. While some discrepancies between techniques still exist, these will undoubtedly be resolved soon, and these different techniques are already providing a great improvement in the understanding of turbulence on a global scale.

Hocking, W. K.↗
NASA NTRS1976

Critical examination of two-equation turbulence closure models

Comparison of the Jones-Launder, Ng-Spalding, Saffman-Wilcox, and Wilcox-Traci two-equation turbulence models has been conducted. It was shown that the Saffman-Wilcox and Wilcox-Traci dissipation-rate formulations admit straightforward integration through the viscous sublayer, whereas integration through the viscous sublayer is a more difficult issue with the Jones-Launder dissipation-function and the Ng-Spalding length-scale formulations. Numerical computations were conducted in which the models were applied to four equilibrium boundary layer flows including adverse, zero, and favorable pressure gradients. Computations of zero pressure gradient flow over a convex wall composed the final part of the comparison.

Chambers, T. L.↗
NASA NTRS1992

Kolmogorov Behavior of Near-Wall Turbulence and Its Application in Turbulence Modeling

The near-wall behavior of turbulence is re-examined in a way different from that proposed by Hanjalic and Launder and followers. It is shown that at a certain distance from the wall, all energetic large eddies will reduce to Kolmogorov eddies (the smallest eddies in turbulence). All the important wall parameters, such as friction velocity, viscous length scale, and mean strain rate at the wall, are characterized by Kolmogorov microscales. According to this Kolmogorov behavior of near-wall turbulence, the turbulence quantities, such as turbulent kinetic energy, dissipation rate, etc. at the location where the large eddies become Kolmogorov eddies, can be estimated by using both direct numerical simulation (DNS) data and asymptotic analysis of near-wall turbulence. This information will provide useful boundary conditions for the turbulent transport equations. As an example, the concept is incorporated in the standard k-epsilon model which is then applied to channel and boundary flows. Using appropriate boundary conditions (based on Kolmogorov behavior of near-wall turbulence), there is no need for any wall-modification to the k-epsilon equations (including model constants). Results compare very well with the DNS and experimental data.

Shih, Tsan-Hsing↗
NASA NTRS1992

Reynolds stress calculations of homogeneous turbulent shear flow with bounded energy states

Reynolds stress calculations of homogeneous turbulent shear flow are conducted with a second-order closure model modified to account for non-equilibrium vortex stretching in the dissipation rate transport equation, as recently proposed by Bernard and Speziale. As with the earlier reported k-epsilon model calculations incorporating this vortex stretching effect, a production-equals-dissipation equilibrium is obtained with bounded turbulent kinetic energy and dissipation. However, this equilibrium is not achieved until the dimensionless time greater than 60, an elapsed time that is at least twice as large as any of those considered in previous numerical and physical experiments on homogeneous shear flow. Direct quantitative comparisons between the model predictions and the results of experiments are quite favorable. In particular, it is shown that the inclusion of this non-equilibrium vortex stretching effect has the capability of explaining the significant range of production to dissipation ratios observed in experiments.

Speziale, Charles G.↗
NASA NTRS1993

Reynolds stress calculations of homogeneous turbulent shear flow with bounded energy states

Reynolds stress calculations of homogeneous turbulent shear flow are conducted with a second-order closure model modified to account for nonequilibrium vortex stretching in the dissipation rate transport equation as recently proposed by Bernard and Speziale (1992). As with the earlier reported K-epsilon model calculations incorporating this vortex stretching effect, a production-equals-dissipation equilibrium is obtained with bounded turbulent kinetic energy and dissipation. However, this equilibrium is now not achieved until the dimensionless time St greater than 60 - an elapsed time that is at least twice as large as any of those considered in previous numerical and physical experiments on homogeneous shear flow. Direct quantitative comparisons between the model predictions and the results of experiments are quite favorable. In particular, it is shown that the inclusion of this nonequilibrium vortex stretching effect has the capability of explaining the significant range of production to dissipation ratios observed in experiments.

Speziale, Charles G.↗
DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗
DOE OSTI2023

Adjustments to the law of the wall above an Amazon forest explained by a spectral link

Modification to the law of the wall represented by a dimensionless correction function ϕ RSL (z/h) is derived using atmospheric turbulence measurements collected at two sites in the Amazon in near-neutral stratification, where z is the distance from the forest floor and h is the mean canopy height. The sites are the Amazon Tall Tower Observatory for z/h∈ [1,2.3] and the Green Ocean Amazon (GoAmazon) site for z/h∈ [1,1.4]. Here, a link between the vertical velocity spectrum E ww (k) (k is the longitudinal wavenumber) and ϕ RSL is then established using a co-spectral budget (CSB) model interpreted by the moving-equilibrium hypothesis. The key finding is that ϕ RSL is determined by the ratio of two turbulent viscosities and is given as ν t,BL /ν t,RSL, where ν t,RSL = (1/A)∫$^{∞}_{0}$ τ(k)E ww (k)dk, ν t,BL = k v (z−d)u * , τ(k) is a scale-dependent decorrelation time scale between velocity components, A = C R /(1−C I ) = 4.5 is predicted from the Rotta constant C R = 1.8, and the isotropization of production constant C I = 3/5 given by rapid distortion theory, k v is the von Kármán constant, u * is the friction velocity at the canopy top, and d is the zero-plane displacement. Because the transfer of energy across scales is conserved in E ww (k) and is determined by the turbulent kinetic energy dissipation rate (ε), the CSB model also predicts that ϕ RSL scales with L BL /L d , where L BL is the length scale of attached eddies to z = d, and L d = u$_{*}^{3}$/ε is a macro-scale dissipation length.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

On Boundary-Value Problems for RANS Equations and Two-Equation Turbulence Models

Currently, in engineering computations for high Reynolds number turbulent flows, turbulence modeling continues to be the most frequently used approach to represent the effects of turbulence. Such models generally rely on solving either one or two transport equations along with the Reynolds-Averaged Navier–Stokes (RANS) equations. The solution of the boundary-value problem of any system of partial differential equations requires the complete delineation of the equations and the boundary conditions, including any special restrictions and conditions. In the literature, such a description is often incomplete, neglecting important details related to the boundary conditions and possible restrictive conditions, such as how to ensure satisfying prescribed values of the dependent variables of the transport equations in the far field of a finite domain. In this article, we discuss the possible influence of boundary values, as well as near-field and far-field behavior, on the solution of the RANS equations coupled with transport equations for turbulence modeling. In so doing, we defne the concept of a welldefined boundary-value problem. Additionally, a three-dimensional, rather than a simpler one-dimensional analysis is performed to analyze the near-wall and far-field behavior of the turbulence model variables. This allows an assessment of the decay rate of these variables required to realize the boundary conditions in the far field. This paper also addresses the impact of various transformations of two-equation models (e.g., the model of Wilcox) to remove the singular behavior of the dissipation rate (ω) at the surface boundary. Finally, the issue of well-posedness regarding the governing equations is considered. A compelling argument (although not a proof) for ill-posedness is made for both direct and inverse problems.

Navier–Stokes equations↗
DOE OSTI2024

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Theoretical and computational study of rotating isotropic turbulence

The various flow regimes that occur in rotating isotropic turbulence and their associated turbulence structure are studied. Direct numerical simulations are conducted for moderate and rapid rotation rates and comparisons are made with the predictions of generalized Eddy Damped Quasi-Normal Markovian (EDQNM) approximations. It is shown that at high rotation rates the nonlinear transfer terms remain small and that the development of the spectrum is through pure viscous decay. This shows that the effect of rapid rotation on the turbulence kinetic energy is through the shutting off of the production term in the dissipation rate equation. At moderate and relatively rapid rates, the rotation causes a discernible anisotropy to develop in the integral length scale with a mild trend toward a two-dimensionalization of the flow. As the flow decays, the Rossby number decreases, leading again to the shutting off of the nonlinear transfer terms. For extremely rapid rotation rates, the anisotropies in the integral length scales are small and direct numerical simulation results are in good agreement with rapid distortion theory.

Mansour, N. N.↗
NASA NTRS1990

Bridge feedback for active damping augmentation

A method is described for broadband damping augmentation of a structural system in which the active members (with feedback control) were developed such that their mechanical input impedance can be electrically adjusted to maximize the energy dissipation rate in the structural system. The active member consists of sensors, an actuator, and a control scheme. A mechanical/electrical analogy is described to model the passive structures and the active members in terms of their impedance representation. As a result, the problem of maximizing dissipative power is analogous to the problem of impedance matching in the electrical network. Closed-loop performance was demonstrated for single- and multiple-active-member controlled truss structure.

Chen, G.-S.↗
NASA NTRS1981

The Eulerian- and Lagrangian-mean flows induced by stationary, dissipating planetary waves

The Eulerian- and the Lagrangian-mean flows induced by stationary, dissipating planetary waves are discussed by employing a simple channel model on a beta-plane. It is assumed that the wave is excited by the bottom undulation and dissipated by Newtonian cooling with relaxation time alpha and by Rayleigh friction with (lambda)(alpha), lambda being constant. Three cases where lambda is equal to one are discussed: (1) the basic zonal wind U sub 0 and the dissipation rate alpha are both constant; (2) U sub 0 varies with height while alpha is constant; and (3) U sub 0 and alpha both vary with height. In case (1), the Eulerian- and the Lagrangian-mean fields are shown to depend on the difference between the dissipation scale-height and the density scale-height. In case (2) and case (3), it is shown that the results for case (1) are modified under slightly more realistic situations.

Takahashi, M.↗
NASA NTRS1978

Analysis of turbulent free jet hydrogen-air diffusion flames with finite chemical reaction rates

The nonequilibrium flow field resulting from the turbulent mixing and combustion of a supersonic axisymmetric hydrogen jet in a supersonic parallel coflowing air stream is analyzed. Effective turbulent transport properties are determined using the (K-epsilon) model. The finite-rate chemistry model considers eight reactions between six chemical species, H, O, H2O, OH, O2, and H2. The governing set of nonlinear partial differential equations is solved by an implicit finite-difference procedure. Radial distributions are obtained at two downstream locations of variables such as turbulent kinetic energy, turbulent dissipation rate, turbulent scale length, and viscosity. The results show that these variables attain peak values at the axis of symmetry. Computed distributions of velocity, temperature, and mass fraction are also given. A direct analytical approach to account for the effect of species concentration fluctuations on the mean production rate of species (the phenomenon of unmixedness) is also presented. However, the use of the method does not seem justified in view of the excessive computer time required to solve the resulting system of equations.

Sislian, J. P.↗
NASA NTRS1991

Source function approach to turbulence bulk properties

A turbulence model to compute bulk properties is presented in which a prescribed source function is used to provide the rate of energy input into the turbulent flow, and the EDQNM model is used to treat the nonlinear transfer in the Navier-Stokes equations. The predictions of the model are tested against (1) the measured Nusselt number versus Rayleigh number relation in turbulent laboratory convection, and (2) the measured bulk kinetic energies and dissipation rates in turbulent channel flow for Reynolds numbers Re = 12,300 and 30,800. In addition, a sensitivity study is performed with respect to the choice of Kolmogorov constant. Generally, the predictions of the model are in reasonable accord with the available experimental data.

Canuto, V. M.↗