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 199 records · Page 11

Experimental and Theoretical Studies of Pulsating Turbulent Flow

The objective of this investigation was to study the effects of small amplitude sinusoidal pulsations on fully developed turbulent flow in a tube from both experimental and theoretical viewpoints. Theoretical models for the macroscopic behavior of pulsating turbulent tube flow were developed for the two cases of very low and very high pulsation frequencies. The models are based on assumptions of quasi-steady and frozen eddy viscosity flow behavior, respectively. The models successfully predict unsteady velocity profiles, thereby supporting the currently proposed definitions of frequency regimes in pulsating turbulent flow. Experimental measurements were made of the time-dependent pressure drop and velocity profiles over the range of frequency-to-Reynolds number ratios from 0.0095 to 0.24. The two macroscopic models developed in this study predict unsteady velocity profiles which are in moderately good agreement with the experiments in their respective frequency regimes, and a previously developed quasi-steady model is found to predict experimental velocity profiles well in both the quasisteady and the frozen eddy viscosity frequency regimes. The effect of flow pulsations on the dissipation of turbulence energy in the vicinity of the wall was measured in the lower transition frequency regime. The long-time averaged dissipation was observed to be unchanged from the steady flow dissipation, within the accuracy of the experiment. A theoretical model of the periodic viscous sublayer was also developed and applied to pulsating flow in a tube, in order to investigate the effects of flow pulsations on the rate of production of turbulence in the region of the wall. The periodic viscous sublayer model predicts sublayer growth periods in steady flow which agree with the published experimental data. When the model is applied to pulsating flow, the response of the sublayer growth period falls into three frequency regimes, the parameters of which are in approximate agreement with the frequency regimes which are defined on the basis of macroscopic flow behavior. The sublayer renewal cycle exhibits quasi-steady flow behavior when the sublayer growth period is much less than the pulsation period, transition behavior when these two periods are approximately equal, and frozen eddy viscosity behavior when the sublayer period is much longer than the pulsation period. The effect of the sublayer growth and renewal cycle on the level of turbulence was investigated by two methods. The velocity fluctuations seen by a point velocity probe located close to the wall were predicted from the model in one method and the rate of turbulence production was estimated from the frequency of sublayer renewal events in the other.

Kingston, G. C.↗

Method and apparatus for generating microshells of refractory materials

A system is described for forming accurately spherical and centered fluid-filled shells, especially of high melting temperature material. Material which is to form the shells is placed in a solid form in a container, and the material is rapidly heated to a molten temperature to avoid recrystallization and the possible generation of unwanted microbubbles in the melt. Immediately after the molten shells are formed, they drop through a drop tower whose upper end is heated along a distance of at least one foot to provide time for dissipation of surface waves on the shells while they cool to a highly viscous, or just above melting temperature so that the bubble within the shell will not rise and become off centered. The rest of the tower is cryogenically cooled to cool the shell to a solid state.

Lee, Mark C.↗

A Physical Model of Moulin Evolution on the Greenland Ice Sheet

Nearly all proglacial water discharge from the 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 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 both the timing and variability of meltwater inputs to the bed. This evolution can impact both the form of the subglacial hydrologic system and associated response of ice motion. Here, we develop a physical model of moulin formation and evolution to constrain the role of englacial processes in shaping the form and structure of the subglacial hydrologic system. Within this model, moulin geometry is controlled by a balance of viscous and elastic deformation and is dependent on that deformation, refreezing, and the dissipation of turbulent and sensible heat energy. All of which are dependent on the characteristics of the available supraglacial meltwater and the surrounding ice. We find moulin geometry is responsive to changes in these parameters over the course of hours to days, indicating that diurnal and multi-day variations in melt can substantially alter the geometry of a moulin and, consequently, the pressure-discharge relationship at the bed of the ice sheet. Therefore, there is no single moulin shape that can appropriately represent englacial storage across the Greenland Ice Sheet.

Andrews, Lauren C.↗

Development of computational methods for heavy lift launch vehicles

The research effort has been focused on the development of an advanced flow solver for complex viscous turbulent flows with shock waves. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. A new computer program named CENS3D has been developed for viscous turbulent flows with discontinuities. Details of the code are described in Appendix A and Appendix B. With the developments of the numerical algorithm and dissipation model, the simulation of three-dimensional viscous compressible flows has become more efficient and accurate. The results of the research are expected to yield a direct impact on the design process of future liquid fueled launch systems.

Yoon, Seokkwan↗

Grid Convergence of High Order Methods for Multiscale Complex Unsteady Viscous Compressible Flows

Grid convergence of several high order methods for the computation of rapidly developing complex unsteady viscous compressible flows with a wide range of physical scales is studied. The recently developed adaptive numerical dissipation control high order methods referred to as the ACM and wavelet filter schemes are compared with a fifth-order weighted ENO (WENO) scheme. The two 2-D compressible full Navier-Stokes models considered do not possess known analytical and experimental data. Fine grid solutions from a standard second-order TVD scheme and a MUSCL scheme with limiters are used as reference solutions. The first model is a 2-D viscous analogue of a shock tube problem which involves complex shock/shear/boundary-layer interactions. The second model is a supersonic reactive flow concerning fuel breakup. The fuel mixing involves circular hydrogen bubbles in air interacting with a planar moving shock wave. Both models contain fine scale structures and are stiff in the sense that even though the unsteadiness of the flows are rapidly developing, extreme grid refinement and time step restrictions are needed to resolve all the flow scales as well as the chemical reaction scales.

Sjoegreen, B.↗

Damping and vibration considerations for the design of optical systems in a launch/space environment

Engineering philosophies for the design of optical systems launched into space and operating in a vacuum or cryovacuum environment are reviewed. Particular attention is given to sources of energy dissipation which are usually lumped under a single modal parameter denoted as the equivalent viscous damping coefficient. Caging and/or damping system components or application of viscoelastic materials and/or dry friction devices are considered to be alternative methods for stabilizing instruments sensitive to motion.

Richard, Ralph M.↗

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows: Minimization of Delta * B Numerical Error - II

The generalization of a class of low-dissipative high order filter finite difference schemes for long time wave propagation of shock/turbulence/combustion compressible viscous gas dynamic flows to compressible MHD equations for structured curvilinear grids has been developed. The new scheme consists of a divergence free preserving high order spatial base scheme with a filter approach which can be divergence-free preserving depending on the type of filter operator being used, the method of applying the filter step, and the type of flow problem to be considered. 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 (Delta * B) numerical error in the sense that no standard divergence cleaning is required. Performance evaluation of these variants, and the key role that the proper treatment of their corresponding numerical boundary conditions can play will be illustrated. Many levels of grid refinement and detailed comparison with several commonly used compressible MHD shock-capturing schemes will be sought. For certain MHD 2-D test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Sjoegreen, Bjoern↗

Divergence Free High Order Filter Methods for the 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 diver- gence cleaning is not required by the present filter approach. For certain MHD test cases, divergence free preservation of the magnetic fields has been achieved.

Yea, H. C.↗

Divergence Free High Order Filter Methods for Multiscale Non-ideal MHD Flows

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

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

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

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↗

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

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

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

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

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

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