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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

High Order Filter Methods for the Non-ideal Compressible MHD Equations

The generalization of a class of low-dissipative high order filter finite difference methods for long time wave propagation of shock/turbulence/combustion compressible viscous gas dynamic flows to compressible MHD equations for structured curvilinear grids has been achieved. The new scheme is shown to provide a natural and efficient way for the minimization of the divergence of the magnetic field numerical error. Standard divergence cleaning is not required by the present filter approach. For certain non-ideal MHD test cases, divergence free preservation of the magnetic fields has been achieved.

Yee, H. C.↗
NASA NTRS2004

Adaptive Low Dissipative High Order Filter Methods for Multiscale MHD Flows

Adaptive low-dissipative high order filter finite difference methods for long time wave propagation of 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 [divergence of 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 NTRS2002

Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Sjoegreen, Bjoern↗
NASA NTRS2000

Progress in the Development of a Class of Efficient Low Dissipative High Order Shock-capturing Methods

In a series of papers, Olsson (1994, 1995), Olsson & Oliger (1994), Strand (1994), Gerritsen Olsson (1996), Yee et al. (1999a,b, 2000) and Sandham & Yee (2000), the issue of nonlinear stability of the compressible Euler and Navier-Stokes Equations, including physical boundaries, and the corresponding development of the discrete analogue of nonlinear stable high order schemes, including boundary schemes, were developed, extended and evaluated for various fluid flows. High order here refers to spatial schemes that are essentially fourth-order or higher away from shock and shear regions. The objective of this paper is to give an overview of the progress of the low dissipative high order shock-capturing schemes proposed by Yee et al. (1999a,b, 2000). This class of schemes consists of simple non-dissipative high order compact or non-compact central spatial differencings and adaptive nonlinear numerical dissipation operators to minimize the use of numerical dissipation. The amount of numerical dissipation is further minimized by applying the scheme to the entropy splitting form of the inviscid flux derivatives, and by rewriting the viscous terms to minimize odd-even decoupling before the application of the central scheme (Sandham & Yee). The efficiency and accuracy of these scheme are compared with spectral, TVD and fifth- order WENO schemes. A new approach of Sjogreen & Yee (2000) utilizing non-orthogonal multi-resolution wavelet basis functions as sensors to dynamically determine the appropriate amount of numerical dissipation to be added to the non-dissipative high order spatial scheme at each grid point will be discussed. Numerical experiments of long time integration of smooth flows, shock-turbulence interactions, direct numerical simulations of a 3-D compressible turbulent plane channel flow, and various mixing layer problems indicate that these schemes are especially suitable for practical complex problems in nonlinear aeroacoustics, rotorcraft dynamics, direct numerical simulation or large eddy simulation of compressible turbulent flows at various speeds including high-speed shock-turbulence interactions, and general long time wave propagation problems. These schemes, including entropy splitting, have also been extended to freestream preserving schemes on curvilinear moving grids for a thermally perfect gas (Vinokur & Yee 2000).

Yee, H. C.↗
NASA NTRS1987

Electron dynamics and potential jump across slow mode shocks

In the de Hoffmann-Teller reference frame, the cross-shock electric field is simply the thermoelectric field responsible for preserving charge neutrality. As such, it gives information regarding the heating and dissipation occurring within the shock. The total cross-shock potential can be determined by integrating a weighted electron pressure gradient through the shock, but this requires knowledge of the density and temperature profiles. Here, a recently proposed alternative approach relying on particle dynamics is exploited to provide an independent estimate of this potential. Both determinations are applied to slow mode shocks which form the plasma sheet boundary in the deep geomagnetic tail as observed by ISEE 3. The two methods correlate well. There is no indication of the expected transition from resistive to viscous shocks, although the highest Mach number shocks show the highest potentials. The implications of these results for the electron dissipation mechanisms and turbulence at the shock are discussed.

Schwartz, Steven J.↗
NASA NTRS1985

A design technique for determining actuator gains in spacecraft vibration control

A design procedure is described which determines the gains of a diagonal damping matrix to control the vibrations of a flexible structure with application to orbiting spacecraft. The procedure is based on minimizing the energy dissipated by control actuators using nonlinear mathematical programming. Each damping gain is assumed to be an active viscous damper and the design process is formulated so that the force or torque output of the actuator does not exceed a specified value. The response of the structure at some specified time after the termination of the disturbance is constrained to be less than some prescribed value based upon spacecraft mission performance requirements. A grillage example is used to demonstrate the design process for determining gains for two representative cases. Resulting designs are verified by a finite element analysis of the structure augmented by the control actuators.

Horner, G. C.↗
NASA NTRS1985

The effects of strain heating in lithospheric stretching models

The deformation by stretching of a continental type lithosphere has been formulated so that the problem can be solved by a continuum mechanical approach. The deformation, stress state, and temperature distribution are constrained to satisfy the physical laws of conservation of mass, energy, momentum, and an experimentally defined rheological response. The conservation of energy equation including a term of strain energy dissipation is given. The continental lithosphere is assumed to have the rheology of an isotropic, incompressible, nonlinear viscous, two layered solid.

Stanton, M.↗
NASA NTRS1999

Weakly Nonlinear Description of Parametric Instabilities in Vibrating Flows

This project focuses on the effects of weak dissipation on vibrational flows in microgravity and in particular on (a) the generation of mean flows through viscous effects and their reaction on the flows themselves, and (b) the effects of finite group velocity and dispersion on the resulting dynamics in large domains. The basic mechanism responsible for the generation of such flows is nonlinear and was identified by Schlichting [21] and Longuet-Higgins. However, only recently has it become possible to describe such flows self-consistently in terms of amplitude equations for the parametrically excited waves coupled to a mean flow equation. The derivation of these equations is nontrivial because the limit of zero viscosity is singular. This project focuses on various aspects of this singular problem (i.e., the limit C equivalent to (nu)((g)(h(exp 3)))exp -1/2 << 1,where nu is the kinematic viscosity and h is the liquid depth) in the weakly nonlinear regime. A number of distinct cases is identified depending on the values of the Bond number, the size of the nonlinear terms, distance above threshold and the length scales of interest. The theory provides a quantitative explanation of a number of experiments on the vibration modes of liquid bridges and related experiments on parametric excitation of capillary waves in containers of both small and large aspect ratio. The following is a summary of results obtained thus far.

Knobloch, E.↗
NASA NTRS2011

A Sub-Hertz, Low-Frequency Vibration Isolation Platform

One of the major technical problems deep-space optical communication (DSOC) systems need to solve is the isolation of the optical terminal from vibrations produced by the spacecraft navigational control system and by the moving parts of onboard instruments. Even under these vibration perturbations, the DSOC transceivers (telescopes) need to be pointed l000 fs of times more accurately than an RF communication system (parabolic antennas). Mechanical resonators have been extensively used to provide vibration isolation for groundbased, airborne, and spaceborne payloads. The effectiveness of these isolation systems is determined mainly by the ability of designing a mechanical oscillator with the lowest possible resonant frequency. The Low-Frequency Vibration Isolation Platform (LFVIP), developed during this effort, aims to reduce the resonant frequency of the mechanical oscillators into the sub-Hertz region in order to maximize the passive isolation afforded by the 40 dB/decade roll-off response of the resonator. The LFVIP also provides tip/tilt functionality for acquisition and tracking of a beacon signal. An active control system is used for platform positioning and for dampening of the mechanical oscillator. The basic idea in the design of the isolation platform is to use a passive isolation strut with an approximately equal to 100-mHz resonance frequency. This will extend the isolation range to lower frequencies. The harmonic oscillator is a second-order lowpass filter for mechanical disturbances. The resonance quality depends on the dissipation mechanisms, which are mainly hysteretic because of the low resonant frequency and the absence of any viscous medium. The LFVIP system is configured using the well-established Stewart Platform, which consists of a top platform connected to a base with six extensible struts (see figure). The struts are attached to the base and to the platform via universal joints, which permit the extension and contraction of the struts. The struts ends are connected in pairs to the base and to the platform, forming an octahedron. The six struts provide the vibration isolation due to the properties of mechanical oscillators that behave as second-order lowpass filters for frequencies above the resonance. At high frequency, the ideal second-order low-pass filter response is spoiled by the distributed mass and the internal modes of membrane and of the platform with its payload.

Ortiz, Gerardo, G.↗
NASA NTRS1974

Farfield viscous effects in nonlinear noise propagation

Discussion of the method of parametric differentiation in application to predictions of farfield noise propagation in both lossless and dissipative media. It is shown that, in the lossless medium, the governing equation, transformed to parameter space, reduces to a wave equation in the farfield. In the dissipative medium, the system of nonlinear partial differential equations, transformed to parameter space, reduces to a linear partial differential equation of the propagating type which contains a third-order derivative as well as the wave operator.

Harris, W. L., Sr.↗
DOE OSTI2024

Self-organization in collisionless, high- β turbulence

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, β $\gg$ 1, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-β, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (~5%) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

A diagonal algorithm for the method of pseudocompressibility

The method of pseudocompressibility has been shown to be an efficient method for obtaining a steady-state solution to the incompressible Navier-Stokes equations. Recent improvements to this method include the use of a diagonal scheme for the inversion of the equations at each iteration. The necessary transformations have been derived for the pseudocompressibility equations in generalized coordinates. The diagonal algorithm reduces the computing time necessary to obtain a steady-state solution by a factor of nearly three. Implicit viscous terms are maintained in the equations, and it has become possible to use fourth-order implicit dissipation. The steady-state solution is unchanged by the approximations resulting from the diagonalization of the equations. Computed results for flow over a two-dimensional backward-facing step and a three-dimensional cylinder mounted normal to a flat plate are presented for both the old and new algorithms. The accuracy and computing efficiency of these algorithms are compared.

Rogers, S. E.↗
NASA NTRS1987

A Diagonal Algorithm for the Method of Pseudocompressibility

The method of pseudocompressibility has been found to be an efficient method for obtaining a steady-state solution to the incompressible Navier-Stokes equations. Recent improvements to this method include the use of a diagonal scheme for the inversion of the equations equations at each iteration. The necessary transformations have been derived for the pseudocompressibility equations in generalized coordinates. The diagonal algorithm reduces the computing time necessary to obtain a steady-state solution by a factor of nearly three. Implicit viscous terms are maintained in the equations, and it has become possible to use fourth-order implicit dissipation. The steady-state solution is unchanged by the approximations resulting from the diagonalization of the equations. Computed results for flow over a two-dimensional backward-facing step and a three-dimensional cylinder mounted normal to a flat plate are presented for both the old and new algorithms. The computing efficiency of these algorithms are compared. Identical solutions are obtained from both algorithms which compare well with experimental results.

Rogers, Stuart E.↗
NASA NTRS2002

Multi-dimensional Upwind Fluctuation Splitting Scheme with Mesh Adaption for Hypersonic Viscous Flow

A multi-dimensional upwind fluctuation splitting scheme is developed and implemented for two-dimensional and axisymmetric formulations of the Navier-Stokes equations on unstructured meshes. Key features of the scheme are the compact stencil, full upwinding, and non-linear discretization which allow for second-order accuracy with enforced positivity. Throughout, the fluctuation splitting scheme is compared to a current state-of-the-art finite volume approach, a second-order, dual mesh upwind flux difference splitting scheme (DMFDSFV), and is shown to produce more accurate results using fewer computer resources for a wide range of test cases. A Blasius flat plate viscous validation case reveals a more accurate upsilon-velocity profile for fluctuation splitting, and the reduced artificial dissipation production is shown relative to DMFDSFV. Remarkably, the fluctuation splitting scheme shows grid converged skin friction coefficients with only five points in the boundary layer for this case. The second half of the report develops a local, compact, anisotropic unstructured mesh adaptation scheme in conjunction with the multi-dimensional upwind solver, exhibiting a characteristic alignment behavior for scalar problems. The adaptation strategy is extended to the two-dimensional and axisymmetric Navier-Stokes equations of motion through the concept of fluctuation minimization.

Wood, William A., III↗
NASA NTRS1984

Ion viscous effects on magnetosonic shock

Profiles of magnetic field and plasma bulk velocity associated with perpendicular shocks are investigated theoretically by solving magnetohydrodynamic equations. The main dissipation mechanism is assumed to be due to ion viscosity. Shock structures with zero temperatures and finite pressure are examined separately. The numerical solutions indicate that ion viscosity tends to smooth profiles of magnetic field and bulk velocity for high Mach number shocks. For smaller ion viscosity and lower Mach number, a damped wave train would appear downstream. The results suggest that kinetic effects such as ion reflections and microturbulence are necessary in order to explain magnetic field overshoots as well as oscillatory magnetic-field structure in high Mach number shocks.

Lin, C. S.↗
NASA NTRS1988

A numerical model for gravity wave dissipation in the thermosphere

Two simplified models have been developed for the internal gravity wave dissipation due to viscosity, thermal conduction, and ion-drag in a multilayered, isothermal thermosphere. Both models use the WKB approximation, ray theory, and the time-averaged equations of gravity wave energy conservation. One model uses all the equations appropriate to a dissipative atmosphere, while the other uses the dispersion equation and polarization relations applicable to a nondissipative atmosphere, neglecting the viscous and thermal conduction contributions to the energy flux. Results from these models are compared to each other and to the results obtained by Klostermeyer (1973), using a full-wave model.

Hickey, M. P.↗
NASA NTRS1980

A calculation procedure for viscous flow in turbomachines, volume 2

Turbulent flow within turbomachines having arbitrary blade geometries is examined. Effects of turbulence are modeled using two equations, one expressing the development of the turbulence kinetic energy and the other its dissipation rate. To account for complicated blade geometries, the flow equations are formulated in terms of a nonorthogonal boundary fitted coordinate system. The analysis is applied to a radial inflow turbine. The solution obtained indicates the severity of the complex interaction mechanism that occurs between the different flow regimes (i.e., boundary layers, recirculating eddies, separation zones, etc.). Comparison with nonviscous flow solutions tend to justify strongly the inadequacy of using the latter with standard boundary layer techniques to obtain viscous flow details within turbomachine rotors. Capabilities and limitations of the present method of analysis are discussed.

Khalil, J.↗
NASA NTRS2001

Application of OVERFLOW to Hypersonic Perfect Gas Flowfields

A validation of the matrix dissipation option in OVERFLOW was completed, and a single set of constants was found which yield good inviscid shock capturing, viscous heat and skin friction predictions, and work over a wide Mach and Reynolds number range. This scheme works well with multigrid, and produces fast convergence to steady state for 2 and 3 dimensional problems, with accuracy comparable to Roe upwinding, and without 'carbuncle' problems for the blunt body flowfields.

Olsen, M. E.↗