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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Aeroacoustic and aerodynamic applications of the theory of nonequilibrium thermodynamics

Recent developments in the field of nonequilibrium thermodynamics associated with viscous flows are examined and related to developments to the understanding of specific phenomena in aerodynamics and aeroacoustics. A key element of the nonequilibrium theory is the principle of minimum entropy production rate for steady dissipative processes near equilibrium, and variational calculus is used to apply this principle to several examples of viscous flow. A review of nonequilibrium thermodynamics and its role in fluid motion are presented. Several formulations are presented of the local entropy production rate and the local energy dissipation rate, two quantities that are of central importance to the theory. These expressions and the principle of minimum entropy production rate for steady viscous flows are used to identify parallel-wall channel flow and irrotational flow as having minimally dissipative velocity distributions. Features of irrotational, steady, viscous flow near an airfoil, such as the effect of trailing-edge radius on circulation, are also found to be compatible with the minimum principle. Finally, the minimum principle is used to interpret the stability of infinitesimal and finite amplitude disturbances in an initially laminar, parallel shear flow, with results that are consistent with experiment and linearized hydrodynamic stability theory. These results suggest that a thermodynamic approach may be useful in unifying the understanding of many diverse phenomena in aerodynamics and aeroacoustics.

Horne, W. Clifton↗
NASA NTRS2002

Extension of Low Dissipative High Order Hydrodynamics Schemes for MHD Equations

The objective of this paper is to extend our recently developed highly parallelizable nonlinear stable high order schemes for complex multiscale hydrodynamic applications to the viscous MHD (magnetohydrodynamic) equations. These schemes employed multiresolution wavelets as adaptive numerical dissipation controls to limit the amount and to aid the selection and/or blending of the appropriate types of dissipation to be used. The new scheme is formulated for both the conservative and non-conservative form of the MHD equations in curvi-linear grids. The three features of the present MHD scheme over existing schemes in the open literature are as follows. First, the scheme is constructed for long-time integrations of shock/turbulence/combustion magnetized flows. Available schemes are too diffusive for long-time integrations and/or turbulence/combustion problems. Second, unlike existing schemes for the conservative MHD equations which suffer from ill-conditioned eigen-decompositions, the present scheme makes use of a well-conditioned eigen-decomposition to solve the conservative form of the MHD equations. This is due to, partly. the fact that the divergence of the magnetic field condition is a different type of constraint from its incompressible Navier-Stokes cousin. Third, a new approach to minimize the numerical error of the divergence free magnetic condition for high order scheme is introduced.

Yee, H. C.↗
NASA NTRS1987

Artificial dissipation and central difference schemes for the Euler and Navier-Stokes equations

An artificial dissipation model, including boundary treatment, that is employed in many central difference schemes for solving the Euler and Navier-Stokes equations is discussed. Modifications of this model such as the eigenvalue scaling suggested by upwind differencing are examined. Multistage time stepping schemes with and without a multigrid method are used to investigate the effects of changes in the dissipation model on accuracy and convergence. Improved accuracy for inviscid and viscous airfoil flow is obtained with the modified eigenvalue scaling. Slower convergence rates are experienced with the multigrid method using such scaling. The rate of convergence is improved by applying a dissipation scaling function that depends on mesh cell aspect ratio.

Swanson, R. C.↗
NASA NTRS1987

Artificial dissipation and central difference schemes for the Euler and Navier-Stokes equations

An artificial dissipation model, including boundary treatment, that is employed in many central difference schemes for solving the Euler and Navier-Stokes equations is discussed. Modifications of this model such as the eigenvalue scaling suggested by upwind differencing are examined.Multistage time stepping schemes with and without a multigrid method are used to investigate the effects of changes in the dissipation model on accuracy and convergence. Improved accuracy for inviscid and viscous airfoil flows is obtained with the modified eigenvalue scaling. Slower convergence rates are experienced with the multigrid method using such scaling. The rate of convergence is improved by applying a dissipation scaling function that depends on mesh cell aspect ratio.

Swanson, R. C.↗
NASA NTRS1995

An assessment of viscous effects in computational simulation of benign and burst vortex flows on generic fighter wind-tunnel models using TEAM code

Vortex flows on a twin-tail and a single-tail modular transonic vortex interaction (MTVI) model, representative of a generic fighter configuration, are computationally simulated in this study using the Three-dimensional Euler/Navier-Stokes Aerodynamic Method (TEAM). The primary objective is to provide an assessment of viscous effects on benign (10 deg angle of attack) and burst (35 deg angle of attack) vortex flow solutions. This study was conducted in support of a NASA project aimed at assessing the viability of using Euler technology to predict aerodynamic characteristics of aircraft configurations at moderate-to-high angles of attack in a preliminary design environment. The TEAM code solves the Euler and Reynolds-average Navier-Stokes equations on patched multiblock structured grids. Its algorithm is based on a cell-centered finite-volume formulation with multistage time-stepping scheme. Viscous effects are assessed by comparing the computed inviscid and viscous solutions with each other and experimental data. Also, results of Euler solution sensitivity to grid density and numerical dissipation are presented for the twin-tail model. The results show that proper accounting of viscous effects is necessary for detailed design and optimization but Euler solutions can provide meaningful guidelines for preliminary design of flight vehicles which exhibit vortex flows in parts of their flight envelope.

Kinard, Tim A.↗
NASA NTRS1985

The Formation of Preplanetary Disks from the Collapse of Rotating Molecular Cloud Cores

Solutions that describe the collapse of a molecular cloud core that is initially in unstable equilibrium, embedded within an envelope of uniform density, and rotating at the same rate as the envelope are given. Hydrodynamic equations, including self gravity, are deduced to a set of ordinary differential equations, which are solved by the method of matched asymptotic expansions. Results of these calculations are: (1) the range of stellar masses derived seems to correspond to realistic ranges of observed stellar masses and interstellar cloud parameters, (2) the proper measure of dissipation rate is the ratio of accretion time to viscous diffusion time, and (3) the pressure distribution on the surface of an accreting protostar is nonuniform in a way that favors the channeling of a stellar wind into a bipolar flow directed along the rotation axis.

Cassen, P.↗
NASA NTRS1993

Carrying the mass flux terms exactly in the first and second moment equations of compressible turbulence

In compressible turbulence models, it is assumed that the Favre-mean velocities are suitable approximations to the Reynolds-mean velocities in order to close unknown terms. This neglects, in the mean momentum and energy equations, the contribution to the stress and work terms by the mean of the fluctuating Favre velocity, a quantity proportional to the turbulent mass flux. As the stress and work terms do not introduce any new unknown correlations requiring closure in either k-epsilon or Reynolds stress closures and because the exact form of the terms can, with little additional work, be carried there is no need to make any modeling assumptions. In the Reynolds stress equations the viscous terms appear naturally in Reynolds variables while the problem is posed in Favre variables. In the process of splitting the viscous terms into the viscous transport terms, carried in Favre variables, and the dissipation terms, carried in Reynolds variables, important contributions from the mass flux appear. The accurate accounting of these terms is important for any consistent near wall modeling and the retention of the mass flux terms is important in complex compressible turbulent flows.

Ristorcelli, J. R.↗
NASA NTRS1992

Effect of mean velocity shear on the dissipation rate of turbulent kinetic energy

The dissipation rate of turbulent kinetic energy in incompressible turbulence is investigated using a two-scale DIA. The dissipation rate is shown to consist of two parts; one corresponds to the dissipation rate used in the current turbulence models of eddy-viscosity type, and another comes from the viscous effect that is closely connected with mean velocity shear. This result can elucidate the physical meaning of the dissipation rate used in the current turbulence models and explain part of the discrepancy in the near-wall dissipation rates between the current turbulence models and direct numerical simulation of the Navier-Stokes equation.

Yoshizawa, Akira↗
NASA NTRS2003

An Approximate Dissipation Function for Large Strain Rubber Thermo-Mechanical Analyses

Mechanically induced viscoelastic dissipation is difficult to compute. When the constitutive model is defined by history integrals, the formula for dissipation is a double convolution integral. Since double convolution integrals are difficult to approximate, coupled thermo-mechanical analyses of highly viscous rubber-like materials cannot be made with most commercial finite element software. In this study, we present a method to approximate the dissipation for history integral constitutive models that represent Maxwell-like materials without approximating the double convolution integral. The method requires that the total stress can be separated into elastic and viscous components, and that the relaxation form of the constitutive law is defined with a Prony series. Numerical data is provided to demonstrate the limitations of this approximate method for determining dissipation. Rubber cylinders with imbedded steel disks and with an imbedded steel ball are dynamically loaded, and the nonuniform heating within the cylinders is computed.

Johnson, Arthur R.↗
NASA NTRS1982

Nutation and polar motion of the Earth

Effects of nonconservative (dissipative) processes in the Earth were studied. Specifically, the magnetic, viscous, and topographic coupling at the core-mantle boundary and the effects of physical dispersion due to anelasticity in the Earth's rheology were investigated and are described.

Smith, M. L.↗
NASA NTRS1987

Critical Mach numbers in classical magnetohydrodynamics

Stationary point analysis is used to compute generalized critical Mach numbers for finite-amplitude fast and slow shocks in classical MHD fluids. Particular attention is paid to the case where the resistive and thermal conduction dissipation scale lengths are comparable and much larger than the viscous scale lengths. With both resistivity and thermal conduction, the critical Mach number at which viscosity must be invoked is determined by the condition that the downstream flow speed equals the isothermal sound speed. It is also shown that resistivity and thermal conduction can provide convergent stationary point solutions for nearly all slow shocks, except perhaps switch-off-shocks.

Kennel, Charles F.↗
NASA NTRS2003

High Order Filter Methods for Shock/Turbulence MHD Flows

Low-dissipative high order filter finite difference methods for shock/turbulence/combustion compressible viscous MHD flows has been constructed. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (del (raided dot) B) numerical error in the sense that no standard divergence cleaning is required. For certain 2-D MHD test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Yee, H. C.↗
NASA NTRS2019

A Physical Model of Moulin Formation and Evolution

Nearly all proglacial water discharge from the present-day Greenland Ice Sheet is routed englacially via moulins. Identification of these moulins in high-resolution imagery is a frequent topic of study, but the processes controlling how and where moulins form, including on past ice sheets for which remote-sensing data are not available, remain poorly understood. Because moulins may reasonably compose approximately 10-15% of the englacial-subglacial hydrologic system, the evolution and shape of moulins can alter the timing of meltwater inputs to the bed. This evolution can impact both the form of the subglacial hydrologic system and the structure of associated geomorphological structures. Here, we develop a physical model of moulin formation and evolution to constrain the role of englacial processes in controlling the form and structure of the subglacial hydrologic system. Ice deformation within and around a moulin is both viscous and elastic, with the rate of turbulent and heat dissipation from water circulation in the moulin controlling both moulin wall melting and warming of the surrounding ice. We find moulin geometry is responsive to changes in these parameters over hours to days, indicating that diurnal and multi-day variations in surface melt can substantially alter the geometry of a moulin and the pressure-discharge relationship at the bed of the ice sheet. These results should be considered carefully when determining surface water inputs for subglacial hydrologic models. In the future, a parameter space study of these results will be combined with an analytic model to create a predictive, stochastic model of moulin and crevasse locations. This future model will be applicable to constraining the potential for surface-to-bed connections in regions where the exact ice-sheet surface morphology is not known, including ice sheets under future warming atmospheric conditions, and paleo ice sheets, where moulins created modern landforms.

Andrews, Lauren C.↗
NASA NTRS1992

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 4: Advanced fan section aerodynamic analysis

The purpose of this study is the development of a three-dimensional Euler/Navier-Stokes flow analysis for fan section/engine geometries containing multiple blade rows and multiple spanwise flow splitters. An existing procedure developed by Dr. J. J. Adamczyk and associates and the NASA Lewis Research Center was modified to accept multiple spanwise splitter geometries and simulate engine core conditions. The procedure was also modified to allow coarse parallelization of the solution algorithm. This document is a final report outlining the development and techniques used in the procedure. The numerical solution is based upon a finite volume technique with a four stage Runge-Kutta time marching procedure. Numerical dissipation is used to gain solution stability but is reduced in viscous dominated flow regions. Local time stepping and implicit residual smoothing are used to increase the rate of convergence. Multiple blade row solutions are based upon the average-passage system of equations. The numerical solutions are performed on an H-type grid system, with meshes being generated by the system (TIGG3D) developed earlier under this contract. The grid generation scheme meets the average-passage requirement of maintaining a common axisymmetric mesh for each blade row grid. The analysis was run on several geometry configurations ranging from one to five blade rows and from one to four radial flow splitters. Pure internal flow solutions were obtained as well as solutions with flow about the cowl/nacelle and various engine core flow conditions. The efficiency of the solution procedure was shown to be the same as the original analysis.

Crook, Andrew J.↗
NASA NTRS2002

Modeling Disturbance Dynamics in Transitional and Turbulent Boundary Layers

The dynamics of an ensemble of linear disturbances in boundary-layer flows at various Reynolds numbers is studied through an analysis of the transport equations for the mean disturbance kinetic energy and energy dissipation rate. Effects of adverse and favorable pressure-gradients on the disturbance dynamics are also included in the analysis. Unlike the fully turbulent regime where nonlinear phase scrambling of the fluctuations affects the flow field even in proximity to the wall, the early stage transition regime fluctuations studied here are influenced across the boundary layer by the solid boundary. In addition, the dominating dynamics in the disturbance kinetic energy equation is governed by the energy production, pressure-transport and viscous diffusion - also in contrast to the fully turbulent regime. For the disturbance dissipation rate, a dynamic balance exists between the destruction and diffusion of dissipation.

Grosch, C. E.↗
NASA NTRS2004

Adaptive Numerical Dissipative Control in High Order Schemes for Multi-D Non-Ideal MHD

The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free of numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multi-resolution wavelets (WAV) (for the above types of flow feature). These filter approaches also provide a natural and efficient way for the minimization of Div(B) numerical error. The filter scheme consists of spatially sixth order or higher non-dissipative spatial difference operators as the base scheme for the inviscid flux derivatives. If necessary, a small amount of high order linear dissipation is used to remove spurious high frequency oscillations. For example, an eighth-order centered linear dissipation (AD8) might be included in conjunction with a spatially sixth-order base scheme. The inviscid difference operator is applied twice for the viscous flux derivatives. After the completion of a full time step of the base scheme step, the solution is adaptively filtered by the product of a 'flow detector' and the 'nonlinear dissipative portion' of a high-resolution shock-capturing scheme. In addition, the scheme independent wavelet flow detector can be used in conjunction with spatially compact, spectral or spectral element type of base schemes. The ACM and wavelet filter schemes using the dissipative portion of a second-order shock-capturing scheme with sixth-order spatial central base scheme for both the inviscid and viscous MHD flux derivatives and a fourth-order Runge-Kutta method are denoted.

Yee, H. C.↗
NASA NTRS2001

Low Dissipative High Order Numerical Simulations of Supersonic Reactive Flows

The objective of this paper is to evaluate the performance of a newly developed low dissipative sixth-order spatial and fourth-order temporal scheme for viscous reactive flows interacting with shock waves that contain fine scale flow structures. The accuracy and efficiency of the scheme, and to what degree the scheme can capture the correct physical wave speeds of stiff reactive flows will be included.

Sjoegreen, B.↗
NASA NTRS2022

Controls on Greenland Moulin Geometry and Evolution From the Moulin Shape Model

Nearly all meltwater from glaciers and ice sheets is routed englacially through moulins, which collectively comprise approximately 10–14% of the efficient englacial–subglacial hydrologic system. Therefore, the geometry and evolution of moulins has the potential to influence subglacial water pressure variations, ice motion, and the runoff hydrograph delivered to the ocean. We develop the Moulin Shape (MouSh) model, a time-evolving model of moulin geometry. MouSh models ice deformation around a moulin using both viscous and elastic rheologies and models melting within the moulin through heat dissipation from turbulent water flow, both above and below the water line. We force MouSh with idealized and realistic surface melt inputs. Our results show that variations in surface melt change the geometry of a moulin by approximately 30% daily and by over 100% seasonally. These size variations cause observable differences in moulin water storage capacity, moulin water levels, and subglacial channel size compared to a static, cylindrical moulin. Our results suggest that moulins are significant storage reservoirs for meltwater, with storage capacity and water levels varying over multiple timescales. Representing moulin geometry within subglacial hydrologic models would therefore improve their accuracy, especially over seasonal periods or in regions where overburden pressures are high.

Moulin↗