Search NASA⌕ Search

SEARCH · Search NASA

Results for “viscous dissipation”

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 163 records · Page 9

Aspects of a high-resolution scheme for the Navier-Stokes equations

In this paper we emphasize the importance of the form of the numerical dissipation model in computing accurate viscous flow solutions. A high-resolution scheme for viscous flows based on three-point central differencing and a matrix dissipation is considered. The various components of this scheme, including 'entropy fix', limiter function, and boundary-point dissipation are discussed. By analyzing boundary-point dissipation stencils, we confirm that with the matrix dissipation model the normal numerical dissipation terms in the streamwise momentum equation are independent of the Reynolds number. Such independence is not achieved with a scalar dissipation form. The accuracy of the central-difference scheme, with and without matrix dissipation, and the flux-difference split scheme of Roe, which is classified as a high-resolution scheme, is compared. For this comparison, three high Reynolds number laminar flows are considered. Solutions of the Navier-Stokes equations are obtained for low-speed flow over a flat plate, transonic flow over an airfoil with transition near the leading edge, and hypersonic flow over a compression ramp. The emphasis of the comparison is primarily on the details of the viscous flows. The necessity of the high-resolution property is revealed.

Swanson, R. C.↗

A theoretical evaluation of rigid baffles in suppression of combustion instability

An analytical technique for the prediction of the effects of rigid baffles on the stability of liquid propellant combustors is presented. This analysis employs both two and three dimensional combustor models characterized by concentrated combustion sources at the chamber injector and a constant Mach number nozzle. An eigenfunction-matching method is used to solve the linearized partial differential equations describing the unsteady flow field for both models. Boundary layer corrections to this unsteady flow are used in a mechanical energy dissipation model to evaluate viscous and turbulence effects within the flow. An integral stability relationship is then employed to predict the decay rate of the oscillations. Results of this analysis agree qualitatively with experimental observations and show that sufficient dissipation exists to indicate that the proper mechanism of baffle damping is a fluid dynamic loss.

Combustion instability↗

Theoretical evaluation of rigid baffles in the suppression of combustion instability

An analytical technique for the prediction of the effects of rigid baffles on the stability of liquid propellant combustors is presented. This analysis employs both two and three dimensional combustor models characterized by concentrated combustion sources at the chamber injector and a constant Mach number nozzle. An eigenfunction-matching method is used to solve the linearized partial differential equations describing the unsteady flow field for both models. Boundary layer corrections to this unsteady flow are in a mechanical energy dissipation model to evaluate viscous and turbulence effects within the flow. An integral instability relationship is then employed to predict the decay rate of the oscillations. Results of this analysis agree qualitatively with experimental observations and show that sufficient dissipation exists to indicate that the proper mechanism of baffle damping is a fluid dynamic loss. The response of the dissipation model to varying baffle blade length, mean flow Mach number, oscillation amplitude, baffle configuration, and oscillation mode is examined.

Combustion instability↗

A theoretical evaluation of rigid baffles in suppression of combustion instability

An analytical technique for the prediction of the effects of rigid baffles on the stability of liquid propellant combustors is presented. This analysis employs a three dimensional combustor model characterized by a concentrated combustion source at the chamber injector and a constant Mach number nozzle. An eigenfunction-matching method is used to solve the linearized partial differential equations describing the unsteady flow field. Boundary layer corrections to this unsteady flow are used in a mechanical energy dissipation model to evaluate viscous and turbulence effects within the flow. An integral stability relationship is then employed to predict the decay rate of the oscillations. Results of this analysis agree qualitatively with experimental observations and show that sufficient dissipation exists to indicate that the proper mechanism of baffle damping is a fluid dynamic loss. The response of the dissipation model to varying baffle blade length, mean flow Mach number and oscillation amplitude is examined.

Combustion instability↗

Small-Scale Dissipation in Binary-Species Transitional Mixing Layers

Motivated by large eddy simulation (LES) modeling of supercritical turbulent flows, transitional states of databases obtained from direct numerical simulations (DNS) of binary-species supercritical temporal mixing layers were examined to understand the subgrid-scale dissipation, and its variation with filter size. Examination of the DSN-scale domain- averaged dissipation confirms previous findings that, out of the three modes of viscous, temperature and species-mass dissipation, the species-mass dissipation is the main contributor to the total dissipation. The results revealed that the percentage of species-mass by total dissipation is nearly invariant across species systems and initial conditions. This dominance of the species-mass dissipation is due to high-density-gradient magnitude (HDGM) regions populating the flow under the supercritical conditions of the simulations; such regions have also been observed in fully turbulent supercritical flows. The domain average being the result of both the local values and the extent of the HDGM regions, the expectations were that the response to filtering would vary with these flow characteristics. All filtering here is performed in the dissipation range of the Kolmogorov spectrum, at filter sizes from 4 to 16 times the DNS grid spacing. The small-scale (subgrid scale, SGS) dissipation was found by subtracting the filtered-field dissipation from the DNS-field dissipation. In contrast to the DNS dissipation, the SGS dissipation is not necessarily positive; negative values indicate backscatter. Backscatter was shown to be spatially widespread in all modes of dissipation and in the total dissipation (25 to 60 percent of the domain). The maximum magnitude of the negative subgrid- scale dissipation was as much as 17 percent of the maximum positive subgrid- scale dissipation, indicating that, not only is backscatter spatially widespread in these flows, but it is considerable in magnitude and cannot be ignored for the purposes of LES modeling. The Smagorinsky model, for example, is unsuited for modeling SGS fluxes in the LES because it cannot render backscatter. With increased filter size, there is only a modest decrease in the spatial extent of backscatter. The implication is that even at large LES grid spacing, the issue of backscatter and related SGS-flux modeling decisions are unavoidable. As a fraction of the total dissipation, the small-scale dissipation is between 10 and 30 percent of the total dissipation for a filter size that is four times the DNS grid spacing, with all OH cases bunched at 10 percent, and the HN cases spanning 24 30 percent. A scale similarity was found in that the domain-average proportion of each small-scale dissipation mode, with respect to the total small-scale dissipation, is very similar to equivalent results at the DNS scale. With increasing filter size, the proportion of the small-scale dissipation in the dissipation increases substantially, although not quite proportionally. When the filter size increases by four-fold, 52 percent for all OH runs, and 70 percent for HN runs, of the dissipation is contained in the subgrid-scale portion with virtually no dependence on the initial conditions of the DNS. The indications from the dissipation analysis are that modeling efforts in LES of thermodynamically supercritical flows should be focused primarily on mass-flux effects, with temperature and viscous effects being secondary. The analysis also reveals a physical justification for scale-similarity type models, although the suitability of these will need to be confirmed in a posteriori studies.

Bellan, Josette↗

Irreversible Entropy Production in Two-Phase Mixing Layers

This report presents a study of dissipation (irreversible production of entropy) in three-dimensional, temporal mixing layers laden with evaporating liquid drops. The purpose of the study is to examine the effects of evaporating drops on the development of turbulent features in flows. Direct numerical simulations were performed to analyze transitional states of three mixing layers: one without drops, and two that included drops at different initial mass loadings. Without drops, the dissipation is essentially due to viscous effects. It was found that in the presence of drops, the largest contribution to dissipation was made by heating and evaporation of the drops, and that at large length scales, this contribution is positive (signifying that the drops reduce turbulence), while at small scales, this contribution is negative (the drops increase turbulence). The second largest contribution to dissipation was found to be associated with the chemical potential, which leads to an increase in turbulence at large scales and a decrease in turbulence at small scales. The next smaller contribution was found to be that of viscosity. The fact that viscosity effects are only third in order of magnitude in the dissipation is in sharp contrast to the situation for the mixing layer without the drops. The next smaller contribution - that of the drag and momentum of the vapor from the drops - was found to be negative at lower mass loading but to become positive at higher mass loading.

Okongo, Nora↗

Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations

The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing nu- merical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten (1978) but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch to the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992) used by Gerritsen and Olsson (1996) in an adaptive mesh refinement method, to determine regions where re nement should be done. The other is the modification of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand alone options for numerical algorithm other than the Yee et al. scheme.

Sjoegreen, B.↗

Propagation of acoustic modes in the transitional ionosphere

Continuum theory is unsatisfactory for describing wave propagation in the transitional ionosphere, since the mean free path is of the same order as the characteristic length (i.e., Knudsen number in this region is of the order of unity). Therefore a transition model in which isotropic electrons behave as a fluid while the ions are governed by the kinetic equation modified by ion-electron collision effects is proposed to study acoustic wave propagation in the upper atmosphere in the altitude region from 500 to 2000 km. The results show that the dissipation of acoustic waves by electrons is through the viscous and thermal conduction effects, and the dissipation of acoustic waves by ions is through Landau damping. Comparisons of the characteristics of acoustic mode propagation and wave dissipation mechanisms in the collision-dominated, transitional from collisional to collisionless, and collisionless media are also discussed.

Hung, R. J.↗

Internal structure of Io and the global distribution of its topography

A global topography is presently calculated for two multilayer Io models in which dissipation occurs in a viscous asthenosphere and a solid mantle: (1) a 'thermal swell' model, in which topography and heat flow are positively correlated, and (2) a 'differentiated lithosphere' model, in which topography and heat flow are negatively correlated. Both the polar topography and the hypsometric distribution of elevations in the differentiated lithosphere model are better matched with observations than the thermal swell model. The shift of the equatorial basin-swell pattern indicates a recent zonal rotation of about 25 deg for Io's lithosphere.

Ross, M. N.↗

Quasi-Static Viscoelastic Finite Element Model of an Aircraft Tire

An elastic large displacement thick-shell mixed finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at the element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses is included in the mixed variational functional. The nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the nonlinear elastic equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to computationally simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Johnson, Arthur R.↗

A Viscoelastic Hybrid Shell Finite Element

An elastic large displacement thick-shell hybrid finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at he element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses in included in the mixed variational functional. Nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to numerically simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Johnson, Arthur↗

Navier-Stokes computation of compressible turbulent flows with a second order closure

The objective was the development of a complete second order closure for wall bounded flows, including all components of the dissipation rate tensor and a numerical solution procedure for the resulting system of equations. The main topics discussed are the closure of the pressure correlations and the viscous destruction terms in the dissipation rate equations and the numerical solution scheme based on a block-tridiagonal solver for the nine equations required for the prediction of plane or axisymmetric flows.

Dingus, C.↗

Studies of unsteady viscous flows using a two-equation model of turbulence

A two equation model of turbulence, based on the turbulent kinetic energy and energy dissipation, suitable for prediction of unsteady viscous flows, was developed. Also, the performance of the two equation model was compared with simpler algebraic models such as the Baldwin-Lomax two layer eddy viscosity model, and a model by Johnson and King which accounts for upstream history of the turbulent kinetic energy. A brief discussion of this study is given.

Sankar, L. N.↗

Coronal heating by the resonant absorption of Alfven waves: The effect of viscous stress tensor

The time-dependent linearized magnetohydrodynamics (MHD) equations for a fully compressible, low-beta, viscoresistive plasma are solved numerically using an implicit integration scheme. The full viscosity stress tensor (Braginskii 1965) is included with the five parameters eta(sub i) i = 0 to 4. In agreement with previous studies, the numerical simulations demonstrate that the dissipation on inhomogeneities in the background Alfven speed occurs in a narrow resonant layer. For an active region in the solar corona the values of eta(sub i) are eta(sub o) = 0.65 g/cm/s, eta(sub 1) = 3.7 x 10(exp -12) g/cm/s, eta(sub 2) = 4 eta(sub 1), eta(sub 3) = 1.4 x 10(exp -6) g/cm/s, eta(sub 4) = 2 eta(sub 3), with n = 10(exp 10)/cu cm, T = 2 x 10(exp 6) K, and B = 100 G. When the Lundquist number S = 10(exp 4) and R(sub 1) much greater than S (where R(sub 1) is the dimensionless shear viscous number) the width of the resistive dissipation layer d(sub r) is 0.22a (where a is the density gradient length scale) and d(sub r) approximately S(exp -1/3). When S much greater than R(sub 1) the shear viscous dissipation layer width d(sub r) scales as R(sub 1)(exp -1/3). The shear viscous and the resistive dissipation occurs in an overlapping narrow region, and the total heating rate is independent of the value of the dissipation parameters in agreement with previous studies. Consequently, the maximum values of the perpendicular velocity and perpendicular magnetic field scale as R(sub 1)(exp -1/3). It is evident from the simulations that for solar parameters the heating due to the compressive viscosity (R(sub 0) = 560) is negligible compared to the resistive and the shear viscous (R(sub 1)) dissipation and it occurs in a broad layer of order a in width. In the solar corona with S approximately equals 10(exp 4) and R(sub 1) approximately equals 10(exp 14) (as calculated from the Braginskii expressions), the shear viscous resonant heating is of comparable magnitude to the resistive resonant heating.

Ofman, L.↗

Nearly sonic and transsonic convective motions in the solar atmosphere related to the solar wind origin

MHD equations are considered for the solar atmosphere. 15 different simplest MHD regimes are indicated for the momentum transport equation depending on the mutual binary interplay between 6 possible and locally dominant terms: non-stationarity and inhomogeneity of the flow, gas pressure, magnetic tensions, viscous and gravity forces. These regimes are delimited by five physically independent dimensionless parameters, for example, Strouhal, sonic Mach, alfvenic Mach, Reynolds and Froude numbers or their combinations. Another partially overlapping classification of the simplest regimes may be introduced based on the energy conservation equation. There are also 15 independent binary combinations between nonstationary and inhomogeneous convective terms, dissipative energy sinks and sources (viscous, heat-conductive, Joule and radiative ones) in the energy equation. More complicated regimes are considered with multiple dominated terms. All these MHD regimes play their important role somewhere in the solar atmosphere complicated by the tensor transport coefficients in the magnetically dominated regions of the upper atmosphere. Nearly sonic and transsonic nonstationary convective motions with ascending and descending flows are observed in the solar chromosphere. the transition region and the lower pans of the solar corona together with related horizontal velocity components. This convection represents a kind of the 'cocoonery' manufacturing nonstationary vortices generated here and partially connected to the photosphere and to the solar wind. The solar wind originates from this powerful transsonic muddle in the solar atmosphere as a tiny fraction of the streamlines which are temporarily getting detached from the 'cocoons; and going to the infinity. The topologically complicated instantaneous 'runaway surface' around the Sun, i.e., the surface which separates outgoing to the infinity streams from other finite flows in the solar atmosphere was not described in the literature and needs additional investigation. We conclude that a simple one-connected smooth and quasistationary 'critical surface' which is often supposed to be placed somewhere in the solar corona (say, at 4 solar radii) in reality does not exist. Observational and theoretical arguments favor instead of this a highly structured, disordered, patched and time-dependent sonic transition surface in the solar atmosphere permanently fluctuating at positional dependent heights starting sometimes from photospheric (and maybe subphotospheric) levels and more frequently from the chromosphenc network up to several solar radii in the solar corona.

Veselovsky, I. S.↗

Numerical simulation of steady and unsteady viscous flow in turbomachinery using pressure based algorithm

The objective of this research is to simulate steady and unsteady viscous flows, including rotor/stator interaction and tip clearance effects in turbomachinery. The numerical formulation for steady flow developed here includes an efficient grid generation scheme, particularly suited to computational grids for the analysis of turbulent turbomachinery flows and tip clearance flows, and a semi-implicit, pressure-based computational fluid dynamics scheme that directly includes artificial dissipation, and is applicable to both viscous and inviscid flows. The values of these artificial dissipation is optimized to achieve accuracy and convergency in the solution. The numerical model is used to investigate the structure of tip clearance flows in a turbine nozzle. The structure of leakage flow is captured accurately, including blade-to-blade variation of all three velocity components, pitch and yaw angles, losses and blade static pressures in the tip clearance region. The simulation also includes evaluation of such quantities of leakage mass flow, vortex strength, losses, dominant leakage flow regions and the spanwise extent affected by the leakage flow. It is demonstrated, through optimization of grid size and artificial dissipation, that the tip clearance flow field can be captured accurately. The above numerical formulation was modified to incorporate time accurate solutions. An inner loop iteration scheme is used at each time step to account for the non-linear effects. The computation of unsteady flow through a flat plate cascade subjected to a transverse gust reveals that the choice of grid spacing and the amount of artificial dissipation is critical for accurate prediction of unsteady phenomena. The rotor-stator interaction problem is simulated by starting the computation upstream of the stator, and the upstream rotor wake is specified from the experimental data. The results show that the stator potential effects have appreciable influence on the upstream rotor wake. The predicted unsteady wake profiles are compared with the available experimental data and the agreement is good. The numerical results are interpreted to draw conclusions on the unsteady wake transport mechanism in the blade passage.

Lakshminarayana, B.↗

A 3-dimensional Navier-Stokes-Euler code for blunt-body flow computations

The shock-layer flowfield is obtained with or without viscous and heat-conducting dissipations from the conservative laws of fluid dynamics equations using a shock-fitting implicity finite-difference technique. The governing equations are cast in curvilinear-orthogonal coordinates and transformed to the domain between the shock and the body. Another set of equations is used for the singular coordinate axis, which, together with a cone generator away from the stagnation point, encloses the computation domain. A time-dependent alternating direction implicit factorization technique is applied to integrate the equations with local-time increment until a steady solution is reached. The shock location is updated after the flowfield computation, but the wall conditions are implemented into the implicit procedure. Innovative procedures are introduced to define the initial flowfield, to treat both perfect and equilibrium gases, to advance the solution on a coarse-to-fine grid sequence, and to start viscous flow computations from their corresponding inviscid solutions. The results are obtained from a grid no greater than 28 by 18 by 7 and converged within 300 integration steps. They are of sufficient accuracy to start parabolized Navier-Stokes or Euler calculations beyond the nose region, to compare with flight and wind-tunnel data, and to evaluate conceptual designs of reentry spacecraft.

Li, C. P.↗

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.↗