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 325 records · Page 18

Navier-Stokes computations of prolate spheroids at angle of attack

Three-dimensional viscous flow calculations are made for a 6:1 prolate spheroid at conditions for which detailed experimental data are available. The computations are made with two finite-volume algorithms for the compressible Navier-Stokes equations, one using central differencing for the convective and pressure terms and the other using an upwind-biased flux-difference-splitting approach. The effects of grid density and artificial dissipation on the accuracy of the numerical results are included. Generally good agreement of the computations with the experimental results is obtained over a range of Reynolds numbers and angles-of-attack, up to 30 deg, although the results at lower Reynolds numbers are sensitive to the assumed transition location.

Vatsa, V. N.↗

High-Order Polynomial Expansions (HOPE) for flux-vector splitting

The Van Leer flux splitting is known to produce excessive numerical dissipation for Navier-Stokes calculations. Researchers attempt to remedy this deficiency by introducing a higher order polynomial expansion (HOPE) for the mass flux. In addition to Van Leer's splitting, a term is introduced so that the mass diffusion error vanishes at M equals 0. Several splittings for pressure are proposed and examined. The effectiveness of the HOPE scheme is illustrated for 1-D hypersonic conical viscous flow and 2-D supersonic shock-wave boundary layer interactions. Also, the authors give the weakness of the scheme and suggest areas for further investigation.

Liou, Meng-Sing↗

Dissipation due to obliquity-induced tides in a Kelvin-Voigt viscoelastic body

A calculation is conducted regarding the local and total dissipation rate of mechanical energy for an incompressible, Kelvin-Voigt body in a circular orbit with finite obliquity. The considered material is characterized by a shear modulus and a viscosity. Of the two parameters, only the shear modulus is known. An investigation is conducted regarding the variation of the total dissipation rate with viscosity. For small values of the viscosity the Kelvin-Voigt material approximates the behavior of an elastic solid. In this range, results for the earth-moon system are similar to those obtained by Peale and Cassen (1978) using a modified form of the elastic solution. For large values of the viscosity the results are qualitatively different from those in the 'elastic' range, as viscous effects dominate the material response.

Vanarsdale, W. E.↗

The connection between Venus' free obliquity and its CMB oblateness

The most striking feature of Venus rotational state is its slow retrograde rotation which is apparently maintained by a balance between solid tidal friction and thermal tidal torques. Solid tides tend to drive the spin toward synchronous rotation while thermal tides drive it away. A balance is achieved at a specific rate because of the inverse frequency dependence of the thermal tide to the semi-diurnal heating. Atmospheric models have been constructed to estimate the thermal tidal torque based on ground heating. The solid friction dissipation factor Q approximately equal to 50 can be deduced assuming rotation has achieved steady state. The most perplexing feature of Venus orientation is its non-zero free obliquity epsilon approximately equal to 1.5 deg relative to its orbit. Although solid tides and perhaps atmospheric tides tend to increase the free obliquity on a time 1/K(sub t) approximately equal to 1 x 10(exp 8) yr, viscous friction (CMF) at a core-mantle boundary (CMB) resulting from the differential angular orientation Delta-epsilon of the core and mantle spin axes should have damped the free obliquity on a time scale as short as 10(exp 6) yr. One means of achieving a balance similar to that controlling rotation is to introduce a comparatively large CMB ellipticity e(sub c) to reduce Delta-epsilon such that there is a balance between solid-thermal tides and CMF. The balance depends not so much on the potential frequency dependence on the tides as on the quadratic dependence of CMF on Delta-epsilon if the layer is turbulent.

Yoder, C. F.↗

Navier-Stokes solutions for transonic flow over a wing mounted in a tunnel

Three-dimensional viscous flow calculations are performed for a swept, NACA 0012 wing mounted inside a wind tunnel for which detailed experimental data is available. A Runge-Kutta time-stepping scheme is used for obtaining steady-state solutions to the thin-layer Navier-Stokes equations. Free-air comptuations are also performed to assess the wall-interference effects. The effects of grid density and artificial dissipation on the accuracy of numerical results are included. The effect of the wind-tunnel sidewall boundary layer on the flow pattern over the wing surface, particularly in the vicinity of the wing/wall juncture is found to be significant.

Vatsa, V. N.↗

Vortex methods for separated flows

The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented in an elementary fashion and includes the relationship with traditional point-vortex studies, the convergence to smooth solutions of the Euler equations, and the essential differences between two- and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. The overlap with the excellent review articles available is kept to a minimum and more emphasis is placed on the area of expertise, namely two-dimensional flows around bluff bodies. When solid walls are present, complete mathematical models are not available and a more heuristic attitude must be adopted. The imposition of inviscid and viscous boundary conditions without conformal mappings or image vortices and the creation of vorticity along solid walls are examined in detail. Methods for boundary-layer treatment and the question of the Kutta condition are discussed. Practical aspects and tips helpful in creating a method that really works are explained. The topics include the robustness of the method and the assessment of accuracy, vortex-core profiles, timemarching schemes, numerical dissipation, and efficient programming. Calculations of flows past streamlined or bluff bodies are used as examples when appropriate.

Spalart, Philippe R.↗

Protoearth mass shedding and the origin of the moon

Darwin's (1980) theory of lunar formation from the earth by means of a rotationally driven dynamic fission instability is presently considered in view of viscous shear's maintenance of solid body rotation throughout the protoearth's accretion phase. Assuming the appropriateness of a polytropic account of the protoearth, it is unlikely that dynamic fission could have occurred; instantaneous spin-up following a giant impact would instead have led to mass shedding. The dynamical phenomenon of mass shedding is here explored on the basis of numerical models for a self-gravitating, axisymmetric, polytropic and dissipative protoearth. It is concluded that mass shedding from the protoearth mantle after a giant impact and explosion could have contributed substantial matter to a lunar disk.

Boss, A. P.↗

A spin-orbit constraint on the viscosity of a Mercurian liquid core

The escape of Mercury from the stable spin-orbit resonance in which the spin angular velocity is twice the orbital mean motion (2n) requires that the kinematic viscosity of a molten core with a laminar boundary layer be comparable to that of water (0.01 sq cm/s) and the tidal Q be less than about 100. If the boundary layer is turbulent, escape from the resonance is only consistent with a liquid core of low viscosity if the critical Reynolds number for the onset of turbulence is above about 500, the moment difference (B - A)/C is below about 0.00001, and the tidal dissipation factor Q is less than about 40. These conclusions depend on the assumptions that Mercury's obliquity was near 0 deg at the time of resonance passage, that the liquid core was not stably stratified at the time at which Mercury passed through the resonance, that a turbulent boundary layer can be characterized by a turbulent or eddy viscosity coefficient, and that the most important coupling between core and mantle is a viscous coupling at a smooth spherical boundary.

Peale, S. J.↗

Weak collisionless shocks mediated by ion gyroviscosity

Collisionless shocks are ubiquitous in space and astrophysical plasmas, and they are essential dynamical features of these systems. Lacking Coulomb collisions, these shocks are mediated by the anomalous dissipation provided by nonlinear plasma instabilities. By numerically resolving the structure of a steady-state, ion gyroviscous shock, we show that ion gyroviscosity, alone, can produce weak (M≲1.1, where M is the sonic Mach number) shocks in a collisionless, magnetized plasma. We emphasize that this effect does not require an appeal to plasma microturbulence. Moreover, while most collisionless systems may be unsuitable to support purely gyroviscous shocks, we argue that gyro-viscous heating may be an overlooked mechanism, generally; and it may be a key driver within magnetohydrodynamic shocks at large. In conclusion, representative examples include the plasma environments produced on the plasma liner experiment and the magnetized liner inertial fusion platforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computation of hypersonic flow through a narrow expansion slot

The compressible Navier-Stokes equations are numerically solved for hypersonic flow over a three-dimensional ramp with a narrow expansion slot. In a two-dimensional test case, it is shown that a higher order scheme is needed to avoid too much numerical dissipation and hence too much total-pressure loss. The calculation demonstrates the role of viscosity in the expansion process of hypersonic flow through the narrow slot. Cases with various wall temperatures and slot widths are studied. Calculations show that wall cooling reduces the thickness of boundary layer, and hence increases the flow expansion substantially. In the lower portion of the slot, inviscidly, the flow is dominated by a highly expanded low density fluid, and viscously, by a viscous layer. As a direct consequence, as soon as the wedge angle is large enough, the mass flux through the slot is almost constant.

Hung, Ching-Mao↗

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre↗

Emergent viscous hydrodynamics from a single quantum particle

We investigate an explicit example of how spatial decoherence can lead to hydrodynamic behavior in the late-time, long-wavelength regime of open quantum systems. We focus on the case of a single nonrelativistic quantum particle linearly coupled to a thermal bath of noninteracting harmonic oscillators at temperature T , a la Caldeira and Leggett. Taking advantage of decoherence in the position representation, we expand the reduced density matrix in powers of the off-diagonal spatial components, so that high-order terms are suppressed at late times. Truncating the resulting power series at second order leads to a set of dissipative transient hydrodynamic equations similar to the nonrelativistic limit of equations widely used in simulations of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions. Transport coefficients are directly determined by the damping constant γ , which quantifies the influence of the environment. The asymptotic limit of our hydrodynamic equations reduces to the celebrated Navier-Stokes equations for a compressible fluid in the presence of a drag force. Furthermore, our results shed new light on the onset of hydrodynamic behavior in open quantum systems where a system with few degrees of freedom is coupled to a large thermal environment.

Hydrodynamics↗

Transonic viscous flow calculations for a turbine cascade with a two equation turbulence model

A numerical method for the study of steady, transonic, turbulent viscous flow through plane turbine cascades is presented. The governing equations are written in Favre-averaged form and closed with a first order model. The turbulent quantities are expressed according to a two-equation kappa-epsilon model where low Reynolds number and compressibility effects are included. The solution is obtained by using a pseudo-unsteady method with improved perturbation propagation properties. The equations are discretized in space by using a finite volume formulation. An explicit multistage dissipative Runge-Kutta algorithm is then used to advance the flow equations in the pseudo-time. First results of calculations compare fairly well with experimental data.

Boretti, A. A.↗

Black hole accretion disks with coronae

Observations suggest the existence of both hot and cold dark matter in the centers of active galactic nuclei. Recent spectral models require a major fraction of power to be dissipated in the hot matter. We study the case when the hot matter forms a corona around a standard cold alpha-disk. In particular, we investigate the case when a major fraction, f, of the power released when the cold matter accretes is transported to and dissipated in the corona. This has major effects on the cold disk, making it colder, more geometrically thin, denser, and having larger optical depths. One important consequence is the disappearance of the effectively optically thin zone as well as of the radiation pressure dominated zone for values of f sufficiently closed to unity. The disappearance of the radiation pressure dominated zone will result in a cold disk with only a gas pressure dominated zone that is stable against thermal and viscous instabilities. We also show that the pressure ( and the radiation) from the corona will only affect the surface layers of the cold disk. Our results disagree with those of other recent work on accretion disks with coronae. We find those works to be based on unphysical assumptions.

Svensson, Roland↗

Numerical Experiment with Time and Spatial Accuracy of Navier-Stokes Computation For Helicopter Problems

Helicopter flowfields are highly unsteady, nonlinear and three-dimensional. In forward flight and in hover, the rotor blades interact with the tip vortex and wake sheet developed by either itself or the other blades. This interaction, known as blade-vortex interactions (BVI), results in unsteady loading of the blades and can cause a distinctive acoustic signature. Accurate and cost-effective computational fluid dynamic solutions that capture blade-vortex interactions can help rotor designers and engineers to predict rotor performance and to develop designs for low acoustic signature. Such a predictive method must preserve a blade's shed vortex for several blade revolutions before being dissipated. A number of researchers have explored the requirements for this task. This paper will outline some new capabilities that have been added to the NASA Ames' OVERFLOW code to improve its overall accuracy for both vortex capturing and unsteady flows. To highlight these improvements, a number of case studies will be presented. These case studies consist of free convection of a 2-dimensional vortex, dynamically pitching 2-D airfoil including light-stall, and a full 3-D unsteady viscous solution of a helicopter rotor in forward flight In this study both central and upwind difference schemes are modified to be more accurate. Central difference scheme is chosen for this simulation because the flowfield is not dominated by strong shocks. The feature of shock-vortex interaction in such a flow is less important than the dominant blade-vortex interaction. The scheme is second-order accurate in time and solves the thin-layer Navier-Stokes equations in fully-implicit manner at each time-step. The spatial accuracy is either second and fourth-order central difference or third-order upwind difference using Roe-flux and MUSCLE scheme. This paper will highlight and demonstrate the methods for several sample cases and for a helicopter rotor. Preliminary computations on a rotor were performed by using this method and are in the process of documentation.

Ahmad, Jasim↗

Three dimensional flow field inside compressor rotor, including blade boundary layers

The flow in a turbomachinery blade passage has a predominant flow direction. The viscous diffusion in the streamwise direction is usually small and the elliptic influence is transmitted upstream through the pressure field. Starting with a guessed pressure field, it is possible to converge on the full elliptic solution by iterating between a parabolic solution and an iteration of the pressure field. The main steps of the calculation are given. The blade boundary layers which are three dimensional with laminar, transitional, turbulent, and separation zones are investigated. The kinetic energy is analyzed, and the dissipation equation is presented. Measurements were made of the three dimensional flow inside an axial flow compressor passage.

Pouagare, M.↗

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible; (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred; (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used; and (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for nonlinear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows. The minimization of employing very fine grids to overcome the production of spurious numerical solutions and/or instability due to under-resolved grids is also sought. The incremental studies to illustrate the performance of the approach are summarized. Extensive testing and full implementation of the approach is forthcoming. The results shown so far are very encouraging.

Yee, Helen C.↗

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible. (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred. (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used. (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for non-linear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows.

Yee, H. C.↗