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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Double-core evolution. 5: Three-dimensional effects in the merger of a red giant with a dwarf companion

The evolution of the common envelope phase of a binary system consisting of a 4.67 solar mass red giant and a 0.94 solar mass dwarf is studied using smoothed particle hydrodynamics. We demonstrate that the three-dimensional effects associated with the gravitational tidal torques lead to a rapid decay of the orbit on timescales approximately less than 1 yr. The relative orbit of the two cores in the common envelope is initally eccentric and tends to circularize as the orbital separation of the two cores decreases. The angular momentum lost from the orbital motion is distributed throughout the common envelope, and the double core does not evolve to a state of co-rotation for the evolutionary time followed. The energy dissipated from the relative orbit and deposited in the common envelope results in the ejection of approximately 13% of the mass of the envelope. The mass is ejected in all directions, but there is a preference for mass ejection in the orbital plane of the binary system. For example, approximately 80% of the ejected mass lies within 30 deg of the binary orbital plane. Because gravitational forces are long range, most of the energy and angular momentum is imparted to a small fraction of the common envelope resulting in an efficiency of the mass ejection process of approximately 15%. The core of the red giant executes significant displacement with respect to the center of mass of the system and contributes nearly equally to the total energy dissipation rate during the latter phases of the evolution. The degree of departure from synchronism of the initial binary system can be an important property of the system which can affect the outcome of the common envelope phase.

Terman, James L.↗

NASA NTRS1994

Performance of a 2-micrometer coherent Doppler lidar for wind measurements

Measurements of boundary layer winds are presented using a 2-micrometer coherent Doppler lidar and the optimal performance of the maximum likelihood estimator. The systematic error for single-shot estimates was estimated as 3.6 cm/s using measurements from a stationary hard target. The estimation error for measurements of the radial component of the wind field was determined, as well as the fraction of the estimates that are randomly distributed over the velocity search space, when the signal power is low and speckle fading is important. The results from actual data are compared with the results from ideal simulations. The first direct estimation of the spatial structure function of the radial wind field and of the energy dissipation rate is presented for both horizontal and vertical directions of propagation. The rms estimation error of the velocity estimates is found to be within 30% of ideal performance based on simulation.

Frehlich, Rod↗

NASA NTRS1993

CUPRI observations of PMSE during Salvo B of NLC-91: Evidence of both partial relfection and turbulent scatter

During the first rocket sequence (called Salvo B) of the NLC-91 campaign, the Cornell University Portable Radar Interferometer (CUPRI) observed two simultaneously occurring layers of Polar Mesophere Summer Echoes (PMSE). during the time of the Turbo B flight, the high time-resolution CUPRI Doppler spectra exhibited sawtooth-like discontinuities in the lower layer which we interpret to be a distorted partial reflection layer which was advected across the radar beam. The upper layer, on the other hand, appeared to be caused by turbulent scatter and we estimate the turbulence energy dissipation rate in the upper layer at the time of the Turbo B flight to have been approximately 0.04 W/kg. Futhermore, a shift in the antenna beam direction from vertical to 8 deg off zenith revealed an aspect sensitivity of approximately 5 dB in the lower layer but none in the upper layer. We conclude that, at this particular time, turbulent scatter was responsible for the upper layer while some form of partial reflection was dominant in the lower layer.

Cho, John Y. N.↗

NASA NTRS1995

Development of methods for predicting large crack growth in elastic-plastic work-hardening materials in fully plastic conditions

The objects of the first, exploratory, stage of the project were listed as: (1) to make a detailed and critical review of the Boundary Element method as already published and with regard to elastic-plastic fracture mechanics, to assess its potential for handling present concepts in two-dimensional and three-dimensional cases. To this was subsequently added the Finite Volume method and certain aspects of the Finite Element method for comparative purposes; (2) to assess the further steps needed to apply the methods so far developed to the general field, covering a practical range of geometries, work hardening materials, and composites: to consider their application under higher temperature conditions; (3) to re-assess the present stage of development of the energy dissipation rate, crack tip opening angle and J-integral models in relation to the possibilities of producing a unified technology with the previous two items; and (4) to report on the feasibility and promise of this combined approach and, if appropriate, make recommendations for the second stage aimed at developing a generalized crack growth technology for its application to real-life problems.

Ford, Hugh↗

NASA NTRS1996

Reduced Stress Tensor and Dissipation and the Transport of Lamb Vector

We develop a methodology to ensure that the stress tensor, regardless of its number of independent components, can be reduced to an exactly equivalent one which has the same number of independent components as the surface force. It is applicable to the momentum balance if the shear viscosity is constant. A direct application of this method to the energy balance also leads to a reduction of the dissipation rate of kinetic energy. Following this procedure, significant saving in analysis and computation may be achieved. For turbulent flows, this strategy immediately implies that a given Reynolds stress model can always be replaced by a reduced one before putting it into computation. Furthermore, we show how the modeling of Reynolds stress tensor can be reduced to that of the mean turbulent Lamb vector alone, which is much simpler. As a first step of this alternative modeling development, we derive the governing equations for the Lamb vector and its square. These equations form a basis of new second-order closure schemes and, we believe, should be favorably compared to that of traditional Reynolds stress transport equation.

Wu, Jie-Zhi↗

NASA NTRS1995

A family of dynamic models for large-eddy simulation

Since its first application, the dynamic procedure has been recognized as an effective means to compute rather than prescribe the unknown coefficients that appear in a subgrid-scale model for Large-Eddy Simulation (LES). The dynamic procedure is usually used to determine the nondimensional coefficient in the Smagorinsky (1963) model. In reality the procedure is quite general and it is not limited to the Smagorinsky model by any theoretical or practical constraints. The purpose of this note is to consider a generalized family of dynamic eddy viscosity models that do not necessarily rely on the local equilibrium assumption built into the Smagorinsky model. By invoking an inertial range assumption, it will be shown that the coefficients in the new models need not be nondimensional. This additional degree of freedom allows the use of models that are scaled on traditionally unknown quantities such as the dissipation rate. In certain cases, the dynamic models with dimensional coefficients are simpler to implement, and allow for a 30% reduction in the number of required filtering operations.

Carati, D.↗

NASA NTRS1997

Modeling of Turbulent Swirling Flows

Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.

Shih, Tsan-Hsing↗

NASA NTRS1997

Dynamical System Analysis of Reynolds Stress Closure Equations

In this paper, we establish the causality between the model coefficients in the standard pressure-strain correlation model and the predicted equilibrium states for homogeneous turbulence. We accomplish this by performing a comprehensive fixed point analysis of the modeled Reynolds stress and dissipation rate equations. The results from this analysis will be very useful for developing improved pressure-strain correlation models to yield observed equilibrium behavior.

Girimaji, Sharath S.↗

NASA NTRS1998

Prediction of Transitional Flows in the Low Pressure Turbine

Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Huang, George↗

NASA NTRS1992

Computation of Turbulent Recirculating Flow in Channels, and for Impingement Cooling

Fully elliptic forms of the transport equations have been solved numerically for two flow configurations. The first is turbulent flow in a channel with transverse rectangular ribs, and the second is impingement cooling of a plane surface. Both flows are relevant to proposed designs for active cooling of hypersonic vehicles using supercritical hydrogen as the coolant. Flow downstream of an abrupt pipe expansion and of a backward-facing step were also solved with various near-wall turbulence models as benchmark problems. A simple form of periodicity boundary condition was used for the channel flow with transverse rectangular ribs. The effects of various parameters on heat transfer in channel flow with transverse ribs and in impingement cooling were investigated using the Yap modified Jones and Launder low Reynolds number k-epsilon turbulence model. For the channel flow, predictions were in adequate agreement with experiment for constant property flow, with the results for friction superior to those for heat transfer. For impingement cooling, the agreement with experiment was generally good, but the results suggest that improved modelling of the dissipation rate of turbulence kinetic energy is required in order to obtain improved heat transfer prediction, especially near the stagnation point. The k-epsilon turbulence model was used to predict the mean flow and heat transfer for constant and variable property flows. The effect of variable properties for channel flow was investigated using the same turbulence model, but comparison with experiment yielded no clear conclusions. Also, the wall function method was modified for use in the variable properties flow with a non-adiabatic surface, and an empirical model is suggested to correctly account for the behavior of the viscous sublayer with heating.

Chang, Byong Hoon↗

NASA NTRS1999

Prediction of Heat and Mass Transfer in a Rotating Ribbed Coolant Passage With a 180 Degree Turn

Numerical results are presented for flow in a rotating internal passage with a 180 degree turn and ribbed walls. Reynolds numbers ranging from 5200 to 7900, and Rotation numbers of 0.0 and 0.24 were considered. The straight sections of the channel have a square cross section, with square ribs spaced one hydraulic diameter (D) apart on two opposite sides. The ribs have a height of 0.1D and are not staggered from one side to the other. The full three dimensional Reynolds Averaged Navier-Stokes equations are solved combined with the Wilcox k-omega turbulence model. By solving an additional equation for mass transfer, it is possible to isolate the effect of buoyancy in the presence of rotation. That is, heat transfer induced buoyancy effects can be eliminated as in naphthalene sublimation experiments. Heat transfer, mass transfer and flow field results are presented with favorable agreement with available experimental data. It is shown that numerically predicting the reattachment between ribs is essential to achieving an accurate prediction of heat/mass transfer. For the low Reynolds numbers considered, the standard turbulence model did not produce reattachment between ribs. By modifying the wall boundary condition on omega, the turbulent specific dissipation rate, much better agreement with the flow structure and heat/ mass transfer was achieved. It is beyond the scope of the present work to make a general recommendation on the omega wall boundary condition. However, the present results suggest that the omega boundary condition should take into account the proximity to abrupt changes in geometry.

Rigby, David L.↗

NASA NTRS1999

The Impact of Buoyancy and Flame Structure on Soot, Radiation and NOx Emissions from a Turbulent Diffusion Flame

It is hypothesized that the spatial structure of a turbulent diffusion flame plays an important role in determining the emissions of radiative energy, soot and NO, from a combustor. This structure, manifested in the two point statistics, is influenced by buoyancy. Radiation, soot and NOx emissions are the cumulative result of processes that occur throughout a flame. For example, radiation fluxes along a line of sight can be found from summing up the contributions from sources in individual pockets of hot soot that emit, and from sinks in cold soot that absorb. Soot and NOx are both the results of slow chemistry and are not equilibrium products. The time that is available for production and burnout is crucial in determining the eventual emissions of these pollutants. Turbulence models generally rely on a single point closure of the appropriate time averaged equations. Hence, spatial information is lost and needs to be modeled using solution variables such as turbulence kinetic energy and dissipation rate, often with the assumption of isotropy. However, buoyancy can affect the physical structure of turbulent flames and can change the spatial extent of soot bearing regions. Theoretical comparisons with models are best done in the limit of infinite Froude number because the inclusion of buoyancy in flow models introduces significant uncertainties. Hence, LII measurements of soot, measurements of radiation fluxes from soot, Particle Imaging Velocimetry (PIV) of the flow field and measurements of post flame NOX will be carried out on the NASA Lewis 2.2 sec drop tower and eventually on the parabolic flight aircraft. The drop rig will be a modified version of a unit that has been successfully used at Lewis in the past.

Kennedy, I. M.↗

NASA NTRS1998

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗

NASA NTRS1999

Asymptotic Behavior of Ensemble-Averaged Linear Disturbances in Homogeneous Shear Flow

In order to expand the predictive capability of single-point turbulence closure models too account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. The calibration is based on homogeneous shear flow where disturbances can be described by rapid distort,ion theory (RDT). The relationship between RDT and linear stability theory is exploit,c d in order to obtain a closed set, of modeled equations. The linear disturbance equations are solved directly so that, the numerical simulation yields a database from which the closure coefficient,s in the ensemble-averaged disturbance equations can he determined.

Thacker, W. D.↗

NASA NTRS1993

Computation of Heat Transfer from Impinging Turbulent Jets

An axisymmetric jet impinging on a flat plate was simulated using a low-Reynolds number k-epsilon model. The accuracy or the numerical calculations was validated by comparing the predicted flowfields with experimental data for various nozzle-plate distance to jet diameter ratios (H/d) and Reynolds numbers. Heat transfer predictions for various H/d showed good agreement with experimental data in the wall jet region, and the complex behavior of heat transfer for small H/d is predicted. However, predictions in the neighborhood of the stagnation point were poor for small H/d The results suggest that improved modeling in the equation for the dissipation rate of the turbulence kinetic energy is required.

Chang, B. H.↗

NASA NTRS2000

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

The fundamental objective of this research is study behavior of aircraft wake vortices within atmospheric boundary layer (ABL) in support of developing the system, Aircraft VOrtex Spacing System (AVOSS), under NASA's Terminal Area Productivity (TAR) program that will control aircraft spacing within the narrow approach corridors of airports. The purpose of the AVOSS system is to increase airport capacity by providing a safe reduction in separation of aircraft compared to the now-existing flight rules. In our first funding period (7 January 19994 - 6 April 1997), we have accomplished extensive model development and validation of ABL simulations. Using the validated model, in our second funding period (7 April 1997 - 6 April 2000) we have investigated the effects of ambient atmospheric turbulence on vortex decay and descent, Crow instability, and wake vortex interaction with the ground. Recognizing the crucial influence of ABL turbulence on wake vortex behavior, we have also developed a software generating vertical profiles of turbulent kinetic energy (TKE) or energy dissipation rate (EDR), which are, in turn, used as input data in the AVOSS prediction algorithms.

Lin, Yuh-Lang↗

NASA NTRS2000

Turbulence Model Predictions of Strongly Curved Flow in a U-Duct

The ability of three types of turbulence models to accurately predict the effects of curvature on the flow in a U-duct is studied. An explicit algebraic stress model performs slightly better than one- or two-equation linear eddy viscosity models, although it is necessary to fully account for the variation of the production-to-dissipation-rate ratio in the algebraic stress model formulation. In their original formulations, none of these turbulence models fully captures the suppressed turbulence near the convex wall, whereas a full Reynolds stress model does. Some of the underlying assumptions used in the development of algebraic stress models are investigated and compared with the computed flowfield from the full Reynolds stress model. Through this analysis, the assumption of Reynolds stress anisotropy equilibrium used in the algebraic stress model formulation is found to be incorrect in regions of strong curvature. By the accounting for the local variation of the principal axes of the strain rate tensor, the explicit algebraic stress model correctly predicts the suppressed turbulence in the outer part of the boundary layer near the convex wall.

Rumsey, Christopher L.↗

NASA NTRS2000

Airborne Coherent Lidar for Advanced In-Flight Measurements (ACLAIM) Flight Testing of the Lidar Sensor

The purpose of the ACLAIM program is ultimately to establish the viability of light detection and ranging (lidar) as a forward-looking sensor for turbulence. The goals of this flight test are to: 1) demonstrate that the ACLAIM lidar system operates reliably in a flight test environment, 2) measure the performance of the lidar as a function of the aerosol backscatter coefficient (beta), 3) use the lidar system to measure atmospheric turbulence and compare these measurements to onboard gust measurements, and 4) make measurements of the aerosol backscatter coefficient, its probability distribution and spatial distribution. The scope of this paper is to briefly describe the ACLAIM system and present examples of ACLAIM operation in flight, including comparisons with independent measurements of wind gusts, gust-induced normal acceleration, and the derived eddy dissipation rate.

Soreide, David C.↗