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 109 records · Page 6

A compressible boundary layer algorithm for use with SINDA '85

It is useful to interface a high-speed-flow solution and SINDA to analyze the thermal behavior of systems that include both conduction and high speed flows. When interfacing a high-speed-flow solution to SINDA, it may be necessary to include the viscous effects in the energy equations. Boundary layer effects of interest include heat transfer coefficients (including convection and viscous dissipation) and friction coefficients. To meet this need, a fast, uncoupled, compressible, two-dimensional, boundary layer algorithm was developed that can model flows with and without separation. This algorithm was used as a subroutine with SINDA. Given the core flow properties and the wall heat flux from SINDA, the boundary layer algorithm returns a wall temperature to SINDA and boundary layer algorithm are iterated until they predict the same wall temperature.

Sakowski, Barbara↗

Plume characteristics of an arcjet thruster

An experimental investigation is conducted of the plume of a 1.5 kW hydrogen-fueled arcjet thruster. The velocity and translational temperature are measured by laser induced fluorescence of excited-state atomic hydrogen at a variety of locations and operating conditions. The peak velocity and temperature at the center of the exit plane of 14.7 km/s and 5000 K lead to a Mach number of 1.9. Evidence is presented which demonstrates the extent to which viscous forces dominate the flow, such as the exit plane Reynolds number of 18. Data taken in the plume suggest a transition from supersonic to subsonic flow via viscous dissipation. In addition, the results suggest that high tank pressures (0.4 to 1.5 Torr) limit the ability to simulate space environments for different portions of the plume.

Liebeskind, John G.↗

Analysis of pultrusion processing for long fiber reinforced thermoplastic composite system

Pultrusion is one of the composite processing technology, commonly recognized as a simple and cost-effective means for the manufacturing of fiber-reinforced, resin matrix composite parts with different regular geometries. Previously, because the majority of the pultruded composite parts were made of thermosetting resin matrix, emphasis of the analysis on the process has been on the conservation of energy from various sources, such as heat conduction and the curing kinetics of the resin system. Analysis on the flow aspect of the process was almost absent in the literature for thermosetting process. With the increasing uses of thermoplastic materials, it is desirable to obtain the detailed velocity and pressure profiles inside the pultrusion die. Using a modified Darcy's law for flow through porous media, closed form analytical solutions for the velocity and pressure distributions inside the pultrusion die are obtained for the first time. This enables us to estimate the magnitude of viscous dissipation and it's effects on the pultruded parts. Pulling forces refined in the pultrusion processing are also analyzed. The analytical model derived in this study can be used to advance our knowledge and control of the pultrusion process for fiber reinforced thermoplastic composite parts.

Tso, W.↗

A new approach to the study of the turbulence layer

Turbulent thermal convection is of considerable importance in fluid dynamical transport phenomena occurring, for example, in the planetary boundary layer of the Earth, the interiors of stars, and accretion disks. In particular, during a significant portion of the evolutionary phase of many stars having convectively unstable cores or outer envelopes, a substantial fraction of energy is transported from the central layers to the outer layers by thermal convection. Moreover, as much of the interior of a star is in highly turbulent motion, a complete theory of stellar structure and evolution requires the explicit consideration of turbulence in order to have expressions for the turbulent quantities arising in the stellar structure equations, and particularly, the turbulent fluxes that appear in the total flux conservation equation, such as the convective flux, kinetic energy flux, etc. A reliable quantification of these fluxes continues to present a challenge in astrophysical fluid dynamics, primarily because astrophysical turbulence is almost always fully-developed and nearly inviscid, and therefore governed by strong nonlinear interactions that distribute the energy among a very wide spectrum of eddies with scales ranging from the characteristic dimension of the flow to those sufficiently small to be affected by viscous dissipation. Furthermore, astrophysical flows are invariably compressible, anisotropic, and inhomogeneous, which requires the consideration of the dynamics of longitudinal modes and their interaction with the transverse modes, as well as complicated boundary conditions. In order to reach a compromise between analytical and numerical tractability and the basic physics of turbulent convection, we have constructed a model of stationary turbulent convection that yields various turbulence statistics, including the convective flux, that are required in stellar evolution models.

Levenson, Ruth A.↗

Hydrodynamics of rotating stars and close binary interactions: Compressible ellipsoid models

We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids, and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest, such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.

Lai, Dong↗

Localization of toroidal motion and shear heating in 3-D high Rayleigh number convection with temperature-dependent viscosity

We have applied spectral-transform methods to study three-dimensional thermal convection with temperature-dependent viscosity. The viscosity varies exponentially with the form exp(-BT), where B controls the viscosity contrast and T is temperature. Solutions for high Rayleigh numbers, up to an effective Ra of 6.25 x 10(exp 6), have been obtained for an aspect-ratio of 5x5x1 and a viscosity contrast of 25. Solutions show the localization of toroidal velocity fields with increasing vigor of convection to a coherent network of shear-zones. Viscous dissipation increases with Rayleigh number and is particularly strong in regions of convergent flows and shear deformation. A time-varying depth-dependent mean-flow is generated because of the correlation between laterally varying viscosity and velocity gradients.

Balachandar, S.↗

Cryogenic Fuel Tank Draining Analysis Model

One of the technological challenges in designing advanced hypersonic aircraft and the next generation of spacecraft is developing reusable flight-weight cryogenic fuel tanks. As an aid in the design and analysis of these cryogenic tanks, a computational fluid dynamics (CFD) model has been developed specifically for the analysis of flow in a cryogenic fuel tank. This model employs the full set of Navier-Stokes equations, except that viscous dissipation is neglected in the energy equation. An explicit finite difference technique in two-dimensional generalized coordinates, approximated to second-order accuracy in both space and time is used. The stiffness resulting from the low Mach number is resolved by using artificial compressibility. The model simulates the transient, two-dimensional draining of a fuel tank cross section. To calculate the slosh wave dynamics the interface between the ullage gas and liquid fuel is modeled as a free surface. Then, experimental data for free convection inside a horizontal cylinder are compared with model results. Finally, cryogenic tank draining calculations are performed with three different wall heat fluxes to demonstrate the effect of wall heat flux on the internal tank flow field.

Greer, Donald↗

Gas-Particle Interactions in a Microgravity Flow Cell

We are developing a microgravity flow cell in which to study the interaction of a flowing gas with relatively massive particles that collide with each other and with the moving boundaries of the cell. The absence of gravity makes possible the independent control of the relative motion of the boundaries and the flow of the gas. The cell will permit gas-particle interactions to be studied over the entire range of flow conditions over which the mixture is not turbulent. Within this range, we shall characterize the viscous dissipation of the energy of the particle fluctuations, measure the influence of particle-phase viscosity on the pressure drop along the cell, and observe the development of localized inhomogeneities that are likely to be associated with the onset of clusters. These measurements and observations should contribute to an understanding of the essential physics of pneumatic transport.

Louge, Michel↗

Wave Driven Non-Linear Flow Oscillator for the 22-Year Solar Cycle

We propose that waves generate an oscillation in the Sun to account for the 22-year magnetic cycle. The mechanism we envision is analogous to that driving the Quasi Biennial Oscillation (QBO) observed in the terrestrial atmosphere, which is well understood in principal. Planetary waves and gravity waves deposit momentum in the background atmosphere and accelerate the flow under viscous dissipation. Analysis shows that such a momentum source represents a non-linearity of third or generally odd order, which generates also the fundamental frequency/period so that an oscillation is maintained without external time dependent forcing. For the Sun, we propose that the wave driven oscillation would occur just below the convection region, where the buoyancy frequency or convective stability becomes small to favor wave breaking and wave mean flow interaction. Using scale analysis to extrapolate from terrestrial to solar conditions, we present results from a simplified analytical model, applied to the equator, that incorporates Hines'Doppler Spread Parameterization for gravity waves (GW). Based on a parametric study, we conclude: (1) Depending on the adopted horizontal wavelengths of GW's, wave amplitudes < 10 m/s can be made to produce oscillating zonal winds of about 25 m/s that should be large enough to generate a corresponding oscillation in the main poloidal magnetic field; (2) The oscillation period can be made to be 22 years provided the buoyancy frequency (stability) is sufficiently small, which would place the oscillating wind field near the base of the convection region; (3) In this region, the turbulence associated with wave processes would be enhanced by low stability, and this also helps to produce the desired oscillation period and generate the dynamo currents that would produce the reversing magnetic field. We suggest that the above mechanism may also drive other long-period metronomes in planetary and stellar interiors.

Mayr, H. G.↗

Thermal Comptonization and Disk Thermal Reprocessing in NGC3516

We present an application of the thermal Comptonization/disk reprocessing model recently proposed by Zdziarski, Lubiniski, and Smith. We show that the absence of strong optical variations in the presence of strong concurrent X-ray variations, similar to those found by Hubble Space Telescope (HST)/Rossi X-Ray Timing Explorer (RXTE) monitoring observations of NGC 3516, can be explained by changing the geometry of the Comptonizing plasma rather than the accretion disk itself. The total X-ray luminosity of the Comptonizing plasma must decrease as its spatial extent increases. In contrast, the disk inner radius must be roughly fixed in order not to produce optical/ultraviolet color variations stronger than observed. By including emission due to internal viscous dissipation in the disk, we can roughly match the optical and X-ray flux levels and variability amplitudes seen from NGC 3516 during the HST/RXTE campaign.

Chiang, James↗

Using Multiwavelength Observations to Determine the Black Hole Mass and Accretion Rate in the Type 1 Seyfert Galaxy NGC 5548

We model the spectral energy distribution of the type 1 Seyfert galaxy NGC 5548, fitting data from simultaneous optical, UV, and X-ray monitoring observations. We assume a geometry consisting of a hot central Comptonizing region surrounded by a thin accretion disk. The properties of the disk and the hot central region are determined by the feedback occurring between the hot Comptonizing region and thermal reprocessing in the disk that, along with viscous dissipation, provides the seed photons for the Comptonization process. The constraints imposed upon this model by the multiwavelength data allow us to derive limits on the central black hole mass, Mu is approximately or less than 2x10(exp 7) solar mass, the accretion rate, Mu is approximately or less than 2.5x10(exp 5) sq solar mass per year/Mu, and the radius of the transition region between the thin outer disk and the geometrically thick, hot inner region, is approximately 2-5x10(exp 14) cm.

Chiang, James↗

Subgrid Scale Modeling in Solar Convection Simulations using the ASH Code

The turbulent solar convection zone has remained one of the most challenging and important subjects in physics. Understanding the complex dynamics in the solar con- vection zone is crucial for gaining insight into the solar dynamo problem. Many solar observatories have generated revealing data with great details of large scale motions in the solar convection zone. For example, a strong di erential rotation is observed: the angular rotation is observed to be faster at the equator than near the poles not only near the solar surface, but also deep in the convection zone. On the other hand, due to the wide range of dynamical scales of turbulence in the solar convection zone, both theory and simulation have limited success. Thus, cutting edge solar models and numerical simulations of the solar convection zone have focused more narrowly on a few key features of the solar convection zone, such as the time-averaged di erential rotation. For example, Brun & Toomre (2002) report computational finding of differential rotation in an anelastic model for solar convection. A critical shortcoming in this model is that the viscous dissipation is based on application of mixing length theory to stellar dynamics with some ad hoc parameter tuning. The goal of our work is to implement the subgrid scale model developed at CTR into the solar simulation code and examine how the differential rotation will be a affected as a result. Specifically, we implement a Smagorinsky-Lilly subgrid scale model into the ASH (anelastic spherical harmonic) code developed over the years by various authors. This paper is organized as follows. In x2 we briefly formulate the anelastic system that describes the solar convection. In x3 we formulate the Smagorinsky-Lilly subgrid scale model for unstably stratifed convection. We then present some preliminary results in x4, where we also provide some conclusions and future directions.

Young, Y.-N.↗

Planetary-Scale Inertio Gravity Waves in the Numerical Spectral Model

In the polar region of the upper mesosphere, horizontal wind oscillations have been observed with periods around 10 hours. Waves with such a period are generated in our Numerical Spectral Model (NSM), and they are identified as planetary-scale inertio gravity waves (IGW). These IGWs have periods between 9 and 11 hours and appear above 60 km in the zonal mean (m = 0), as well as in zonal wavenumbers m = 1 to 4. The waves can propagate eastward and westward and have vertical wavelengths around 25 km. The amplitudes in the wind field are typically between 10 and 20 m/s and can reach 30 m/s in the westward propagating component for m = 1 at the poles. In the temperature perturbations, the wave amplitudes above 100 km are typically 5 K and as large as 10 K for m = 0 at the poles. The IGWs are intermittent but reveal systematic seasonal variations, with the largest amplitudes occurring generally in late winter and spring. In the NSM, the IGW are generated like the planetary waves (PW). They are produced apparently by the instabilities that arise in the zonal mean circulation. Relative to the PWs, however, the IGWs propagate zonally with much larger velocities, such that they are not affected much by interactions with the background zonal winds. Since the IGWs can propagate through the mesosphere without much interaction, except for viscous dissipation, one should then expect that they reach the thermosphere with significant and measurable amplitudes.

Mayr, H. G.↗

Endothelial, cardiac muscle and skeletal muscle exhibit different viscous and elastic properties as determined by atomic force microscopy

This study evaluated the hypothesis that, due to functional and structural differences, the apparent elastic modulus and viscous behavior of cardiac and skeletal muscle and vascular endothelium would differ. To accurately determine the elastic modulus, the contribution of probe velocity, indentation depth, and the assumed shape of the probe were examined. Hysteresis was observed at high indentation velocities arising from viscous effects. Irreversible deformation was not observed for endothelial cells and hysteresis was negligible below 1 microm/s. For skeletal muscle and cardiac muscle cells, hysteresis was negligible below 0.25 microm/s. Viscous dissipation for endothelial and cardiac muscle cells was higher than for skeletal muscle cells. The calculated elastic modulus was most sensitive to the assumed probe geometry for the first 60 nm of indentation for the three cell types. Modeling the probe as a blunt cone-spherical cap resulted in variation in elastic modulus with indentation depth that was less than that calculated by treating the probe as a conical tip. Substrate contributions were negligible since the elastic modulus reached a steady value for indentations above 60 nm and the probe never indented more than 10% of the cell thickness. Cardiac cells were the stiffest (100.3+/-10.7 kPa), the skeletal muscle cells were intermediate (24.7+/-3.5 kPa), and the endothelial cells were the softest with a range of elastic moduli (1.4+/-0.1 to 6.8+/-0.4 kPa) depending on the location of the cell surface tested. Cardiac and skeletal muscle exhibited nonlinear elastic behavior. These passive mechanical properties are generally consistent with the function of these different cell types.

Non-NASA Center↗

Rotor-Stator Interaction Performance Effects

Decreased axial spacing between blade rows in an axial compressor stage is thought to increase stage performance because of an unsteady process that occurs in the downstream blade row and acts on the upstream blade row wakes. This process results in the "recovery" of part of the wake energy before all of this energy is irreversibly lost due to viscous diffusion. To study the wake-blade interaction mechanism, researchers at the NASA Lewis Research Center acquired two-component Laser Fringe Anemometer measurements of the rotor wake in the single-stage transonic compressor at two stage loading levels. The detailed measurements were acquired for one stator pitch in circumference at axial positions from the rotor trailing edge to 20 percent of the stator axial chord, at the exit of the stator passage, and downstream of the stator row including the stator wake. These data show that the changes in wake energy that occur inside the stator passage are not due to viscous dissipation alone, and thus the data provide evidence that "wake recovery" is occurring. A time-accurate, three-dimensional Navier Stokes simulation of the compressor stator was done at the corresponding stage loading levels. The measurements and simulations are being used in combination to show the effects of stator blade loading, quantify the effects of viscosity, and quantify the stage efficiency gain due to the wake recovery process. The accuracy of simple models of the wake recovery process is being evaluated in an effort to include the effects of wake recovery in the NASA-developed Average Passage code for multistage turbomachinery simulations.

VanZante, Dale E.↗

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗

V344 Lyrae: A Touchtone SU UMa Cataclysmic Variable in the Kepler Field

We report on the analysis of the Kepler short-cadence (SC) light curve of V344 Lyr obtained during 2009 June 20 through 2010 Mar 19 (Q2-Q4). The system is an SU UMa star showing dwarf nova outbursts and superoutbursts, and promises to be a touchstone for CV studies for the foreseeable future. The system displays both positive and negative super humps with periods of 2.20 and 2.06-hr, respectively, and we identify an orbital period of 2.11-hr. The positive super humps have a maximum amplitude of approximately 0.025-mag, the negative super humps a maximum amplitude of approximately 0.8 mag, and the orbital period at quiescence has an amplitude of approximately 0.025 mag. The quality of the Kepler data is such that we can test vigorously the models for accretion disk dynamics that have been emerging in the past several years. The SC data for V344 Lyr are consistent with the model that two physical sources yield positive superhumps: early in the superoutburst. the superhump signal is generated by viscous dissipation within the periodically flexing disk. but late in the superoutburst, the signal is generated as the accretion stream bright spot sweeps around the rim of the non-axisymmetric disk. The disk super humps are roughly anti-phased with the stream/late superhumps. The V344 Lyr data also reveal negative super humps arising from accretion onto a tilted disk precessing in the retrograde direction, and suggest that negative superhumps may appear during the decline of DN outbursts. The period of negative superhumps has a positive P(raised dot) in between outbursts.

Smale, Alan P.↗

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for high speed lunar return experiments conducted in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The nonequilibrium radiative transport equations are solved via NASA's in-house radiation code NEQAIR using flow field input from US3D flow solver. The post-shock flow properties for the 10 km/s Earth entry conditions are computed using the stagnation line of a blunt-body and a full facility CFD (Computational Fluid Dynamics) simulation of the EAST shock tube. The shocked gas in the blunt-body flow achieves a thermochemical equilibrium away from the shock front whereas EAST flow exhibits a nonequilibrium behavior due to strong viscous dissipation of the shock by boundary layer. The full-tube flow calculations capture the influence of the boundary layer on the shocked gas state and provide a realistic fluid dynamic input for the radiative predictions. The integrated radiance behind the shock is calculated in NEQAIR for wavelength regimes from Vacuum-UltraViolet (VUV) to InfraRed (IR), which are pertinent to the emission characteristics of high enthalpy shock waves in air. These radiance profiles are validated against corresponding EAST shots. The full-tube simulations successfully predict a sharp radiance peak at the shock front which gets smeared in the test data due to the spatial resolution in the measurements. The full facility based radiance behind the shock shows a slightly better match with the test data in the VUV and Red spectral regions, as compared to that from a blunt-body based predictions. The UV radiance is very similar for both geometries and under-predicts the test behavior. The IR test data matches better with the blunt-body based predictions where the full-tube simulations show a significant over-prediction.

Radiation↗