SEARCH · Search NASA

Results for “viscosity and deformation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Computations of Boiling in Microgravity

Tryggvason, Gretar↗

NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗

NASA NTRS1994

Control of cytoskeletal mechanics by extracellular matrix, cell shape, and mechanical tension

We have investigated how extracellular matrix (ECM) alters the mechanical properties of the cytoskeleton (CSK). Mechanical stresses were applied to integrin receptors on the apical surfaces of adherent endothelial cells using RGD-coated ferromagnetic microbeads (5.5-microns diameter) in conjunction with a magnetic twisting device. Increasing the number of basal cell-ECM contacts by raising the fibronectin (FN) coating density from 10 to 500 ng/cm2 promoted cell spreading by fivefold and increased CSK stiffness, apparent viscosity, and permanent deformation all by more than twofold, as measured in response to maximal stress (40 dyne/cm2). When the applied stress was increased from 7 to 40 dyne/cm2, the stiffness and apparent viscosity of the CSK increased in parallel, although cell shape, ECM contacts, nor permanent deformation was altered. Application of the same stresses over a lower number ECM contacts using smaller beads (1.4-microns diameter) resulted in decreased CSK stiffness and apparent viscosity, confirming that this technique probes into the depth of the CSK and not just the cortical membrane. When magnetic measurements were carried out using cells whose membranes were disrupted and ATP stores depleted using saponin, CSK stiffness and apparent viscosity were found to rise by approximately 20%, whereas permanent deformation decreased by more than half. Addition of ATP (250 microM) under conditions that promote CSK tension generation in membrane-permeabilized cells resulted in decreases in CSK stiffness and apparent viscosity that could be detected within 2 min after ATP addition, before any measurable change in cell size.(ABSTRACT TRUNCATED AT 250 WORDS).

Non-NASA Center↗

NASA NTRS1994

Modeling of periodic great earthquakes on the San Andreas fault: Effects of nonlinear crustal rheology

We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.

Reches, Ze'ev↗

NASA NTRS1994

Viscosity estimates for the crust and upper mantle from patterns of lacustrine shoreline deformation in the Eastern Great Basin

The deformed shorelines of Lake Bonneville constitute a classic source of information on lithospheric elastic thickness and upper mantle viscosity. We describe and apply a new model to a recently augmented data set. New data better constrain both the complex spatio-temporal pattern of the lake load and the crustal deformation response to that load. The history of lake level fluctuations has been significantly refined and somewhat modified. This is due to both more radiocarbon dates from within the Bonneville basin and to an improved calibration of the radiocarbon timescale itself. The data which constrain the crustal deformation pattern consist of ages and shoreline elevations from several hundred points which sample three major levels of Lake Bonneville and corresponding elevations from the high stands of three smaller lakes situated to the west of Lake Bonneville. The geometry of the Earth model incorporates an arbitrary number of layers overlying a half-space, and the rheology of each level can accommodate an arbitrary number of Maxwell viscoelastic elements in parallel. The inverse modeling comprises three complementary approaches: for the simplest configurations, we performed a direct search of the parameter space and delineated the irregular boundary of the subspace of acceptable models. For more complex configurations, we constrained the elastic parameters to their seismically determined values and then solved for viscosity versus depth profiles by either expressing the log(viscosity) versus log(depth) profile as a series of specially constructed orhtogonal polynomials, or by allowing each of 8-10 layers (plus the half-space) to have an independently determined viscosity. We found that the data do not strongly support (nor can they conclusively exclude) a more complex rheology than simple Maxwell viscoelasticity. The orthogonal polynomial solution exhibits an essentially monotonic decrease in viscosity with depth.

Bills, Bruce G.↗

NASA NTRS1980

Equilibrium geochemical modeling of a seasonal thermal energy storage aquifer field test

A geochemical mathematical modeling study designed to investigate the well plugging problems encountered at the Auburn University experimental field tests is summarized. The results, primarily of qualitative interest, include: (1) loss of injectivity was probably due to a combination of native particulate plugging and clay swelling and dispersion; (2) fluid-fluid incompatibilities, hydrothermal reactions, and oxidation reactions were of insignificant magnitude or too slow to have contributed markedly to the plugging; and (3) the potential for and contributions from temperature-induced dissolved gas solubility reductions, capillary boundary layer viscosity increases, and microstructural deformation cannot be deconvolved from the available data.

Stottlemyre, J. S.↗

NASA NTRS1980

Time dependent deformation and stress in the lithosphere

Efficient computer programs incorporating frontal solution and time stepping procedure were developed for the modelling of geodynamic problems. This scheme allows for investigating the quasi static phenomena including the effects of the rheological structure of a tectonically active region. From three dimensional models of strike slip earthquakes, it was found that lateral variation of viscosity affects the characteristics of surface deformations. The vertical deformation is especially informative about the viscosity structure in a strike slip fault zone. A three dimensional viscoelastic model of a thrust earthquake indicated that the transient disturbance on plate velocity due to a great plate boundary earthquake is significant at intermediate distances, but becomes barely measurable 1000 km away from the source.

Yang, M.↗

NASA NTRS1980

Extensional flow of dilute polymer solution

Simultaneous effects of temperature and concentration on tensile viscosity are investigated using a semi-empirical scaling law established to relate the concentration, temperature and deformation rate by a reduced deformation rate variable. Results show the onset of high viscosity occurring in the vicinity of the onset of non-Newtonian viscosity observed in a simple shear flow experiment. High viscosity is also compared with steady and dynamic shear data. Finally, the asymptotic solution of Hassager's kinetic theory of a bead-rod model for macromolecular solutions is compared with the data, indicating that tensile viscosity depends strongly on the molecular weight, the higher the molecular weight, the more sensitive the tensile viscosity in terms of the deformation rate.

Peng, S. T. J.↗

NASA NTRS1975

Recent and late quaternary changes in water level

Water level changes of both the Great Lakes and the sea are described along with methods of analyzing water level data. The influence of elastic deformation of the earth and viscosity is discussed. Causes of water level changes reviewed include: earth movements, geoid changes, storm surges or meteorological phenomena, and melting ice in Antarctica, Greenland, and the mountain glaciers.

R. I. Walcott↗

NASA NTRS2000

Reduction of Effective Acceleration to Microgravity Levels

Acceleration due to earth's gravity causes buoyancy driven convection and sedimentation in solutions. In addition. pressure gradients occur as a function of the height within a liquid column. Hence gravity effects both equilbria conditions and phase transitions as a result of hydrostatic pressure gradients. The affect of gravity on the rate of heat and man transfer in solutal processes can be particularly important in polymer processing due to the high sensitivity of polymeric materials to processing conditions. The term microgravity has been coined to describe an environment in which the affects of gravitational acceleration am greatly reduced. It may seem odd to talk in term of reducing the effects of gravitational acceleration since gravitational attraction is a basic property of matter. However, die presence of gravity on in situ processing or measurements can be negated by achieving conditions in which the laboratory, or more specifically the container of the experimental materials, a subjected to the same acceleration as the materials themselves. With regard to the laboratory reference frame, there is virtually no force on the experimental solutions. This is difficult to achieve but can be done. A short review of Newtonian physics provides an explanation on both how processes we affected by gravity and how microgravity conditions are achieved. The fact that fluids deform when subject to a force bid solids do not indicates that solids have a structure able to exert an opposing force that negates an externally applied force. Liquids deform when a force is applied, indicating that a liquid structure cannot completely negate an applied force. Just how easily a liquid resists deformation is related to its viscosity. Spaceflight provides an environment in which the laboratory reference frame i.e. the spacecraft and all the equipment therein an experiencing virtually identical forces. There is no solid foundation underneath such a laboratory, so the laboratory accelerates according to the force of gravity as do the experimental fluids within the lab. Hence, the magnitude of the form excited by the laboratory on the experimental solutions within are greatly reduced. When compared with a laboratory on the ground and averaged over time, the fluids in a spaceflight laboratory experience approximately a 10 (sup -6)decrease in acceleration relative to their laboratory reference frame hence the term microgravity.

Downey, James P.↗

NASA NTRS2002

Laminar-Turbulent Transition in an Electromagnetically-Levitated Droplet

During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop can be inferred from motion of impurity particulates on the surface. In the experiments, subsequent to melting, Joule heating produces a continuous reduction of viscosity of the fluid resulting in an acceleration of the flow with time. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal drop that undergo flow instabilities. As the fluid temperature rises, the amplitude of the secondary flow increases; and beyond a point, the tracers exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. Assuming that the observed laminar-turbulent transition is shear-layer type, the internal structure of the toroidal loops is used to develop a semi- empirical correlation for the onset of turbulence. Our calculations indicate that the suggested correlation is in modest agreement with the experimental data, with the transition occurring at a Reynolds number of 600.

Hyers, R. W.↗

NASA NTRS2002

Laminar-Turbulent Transition in an Electromagnetically-Levitated Droplet

During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop can be inferred from motion of impurity particulates on the surface. In the experiments, subsequent to melting, Joule heating produces a continuous reduction of viscosity of the fluid resulting in an acceleration of the flow with time. These observations indicate formation of a pair of corotating toroidal flow structures inside the spheroidal drop that undergo flow instabilities. As the fluid temperature rises, the amplitude of the secondary flow increases; and beyond a point, the tracers exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. Assuming that the observed laminar-turbulent transition is shear-layer type, the internal structure of the toroidal loops is used to develop a semi-empirical correlation for the onset of turbulence. Our calculations indicate that the suggested correlation is in modest agreement with the experimental data, with the transition occurring at a Reynolds number of 600.

Hyers, Robert W.↗

NASA NTRS1984

Lubrication of Machine Elements

The understanding of hydrodynamic lubrication began with the classical experiments of Tower and Petrov. Reynolds used a reduced form of the Navier-Stokes equations and the continuity equation to generate a second order differential equation for the pressure in the narrow, converging gap of a bearing contact. Such a pressure enables a load to be transmitted between the surfaces with very low friction since the surfaces are completely separated by a film of fluid. In such a situation it is the physical properties of the lubricant, notably the dynamic viscosity, that dictate the behavior of the contact. The understanding of boundary lubrication is normally attributed to Hardy and Doubleday. In boundary lubrication it is the physical and chemical properties of thin films of molecular proportions and the surfaces to which they are attached that determine contact behavior. The lubricant viscosity is not an influential parameter. Research is devoted to a better understanding and more precise definition of other lubrication regimes between these extremes. One such regime, elastohydrodynamic lubrication, occurs in nonconformal contacts, where the pressures are high and the bearing surfaces deform elastically. In this situation the viscosity of the lubricant may raise considerably, and this further assists the formation of an effective fluid film. The science of these three lubrication regimes (hydrodynamic, elastohydrodynamic, and boundary) are described and the manner in which this science is used in the design of machine elements is examined.

Hamrock, Bernard J.↗

NASA NTRS1999

Correction to Transient Creep of a Composite Lower Crust 2: A Polymineralic Basis for Rapidly Evolving Postseismic Deformation Modes

Ivins used three relations for estimating the effective viscosity of quartzites and crustal rocks having strong quartz affinity.

movements seismic deformations↗

NASA NTRS2002

Spin-Up Instability of a Levitated Molten Drop in MHD-Flow Transition to Turbulence

When an alternating magnetic field interacts with induced eddy currents in a conducting body, there will be a repulsive force between the body and the driving coil system generating the field. This repulsive force is the basis of electromagnetic levitation, which allows containerless processing of different materials. The eddy currents in the conducting body also generate Joule heating. Axial rotation of electromagnetically levitated objects is a common observation in levitation systems and often an undesirable side effect of such experiments on 1-g and -g. There have been recent efforts to use magnetic damping and suppress this tendency of body rotation. The first report of rotation in EML drops was attributed to a slight asymmetry of the shape and location of the levitation coils could change the axis and speed of rotation. Other theories of sample rotation include a frequency difference in the traveling electromagnetic waves and a phase difference in two different applied fields of the same frequency. All of these different mechanisms share the following characteristics: the torque is small, constant for constant field strength, and very weakly dependent on the sample's temperature and phase (solid or liquid). During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop is inferred from motion of impurities on the surface. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal levitated drop that undergo secondary flow instabilities. As rise in the fluid temperature rises, the viscosity falls and the internal flow accelerates and becomes oscillatory; and beyond a point in the experiments, the surface impurities exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. In this work, a background of these set of observations will be given followed by a presentation of our results on the digital particle tracking analysis that has been performed on a number of available videos. The analysis indicates that the levitated drop attains a constant rotational speed during the melting phase and formation of the co-rotating axi-symmetric laminar toroidal structures. However, the rate of axial rotation increases dramatically during the deformation of the toroidal structures anti their breakup into chaotic entities. This new data suggests an interaction between the flow inside the levitated molten drop and the driving coils in the experiments. Possible mechanisms for this interaction will be reviewed. The data will also be used to make an assessment of existing theories on droplet rotation.

Abedian, B.↗

NASA NTRS1996

Multigrid Methods for EHL Problems

In many bearings and contacts, forces are transmitted through thin continuous fluid films which separate two contacting elements. Objects in contact are normally subjected to friction and wear which can be reduced effectively by using lubricants. If the lubricant film is sufficiently thin to prevent the opposing solids from coming into contact and carries the entire load, then we have hydrodynamic lubrication, where the lubricant film is determined by the motion and geometry of the solids. However, for loaded contacts of low geometrical conformity, such as gears, rolling contact bearings and cams, this is not the case due to high pressures and this is referred to as Elasto-Hydrodynamic Lubrication (EHL) In EHL, elastic deformation of the contacting elements and the increase in fluid viscosity with pressure are very significant and cannot be ignored. Since the deformation results in changing the geometry of the lubricating film, which in turn determines the pressure distribution, an EHL mathematical model must simultaneously satisfy the complex elasticity (integral) and the Reynolds lubrication (differential) equations. The nonlinear and coupled nature of the two equations makes numerical calculations computationally intensive. This is especially true for highly loaded problems found in practice. One novel feature of these problems is that the solution may exhibit sharp pressure spikes in the outlet region. To this date both finite element and finite difference methods have been used to solve EHL problems with perhaps greater emphasis on the use of the finite difference approach. In both cases, a major computational difficulty is ensuring convergence of the nonlinear equations solver to a steady state solution. Two successful methods for achieving this are direct iteration and multigrid methods. Direct iteration methods (e.g Gauss Seidel) have long been used in conjunction with finite difference discretizations on regular meshes. Perhaps one of the best examples of the application of such methods is the recent Effective Influence Method of Dowson and Wang. Multigrid methods have also been used with great success by Venner and Venner and Lubrecht with a good summary being given by Venner. As both these finite difference discretization based approaches appear to provide an efficient way of solving EHL problems, it is important to understand their relative merits. This paper is a first attempt at providing such an understanding in the context of EHL point contact problem, (contact of two spheres), in which the contact zone is a point and an ellipse or circle for unloaded and loaded dry contacts respectively. Since the film thickness and the contact width are generally small compared to the local radius of curvature of the two surfaces, the reduced geometry of the surfaces in the contact area can be accurately approximated to the contact between a paraboloid and a flat surface. The layout of the remainder of this paper is as follows. In section 2 we introduce the form of the equations to be solved. The Effective Influence Newton Method is described in Section 3 while Section 4 describes the Multigrid method to be used. Sections 5 and 6 describe the test problems to be used in the comparison between the two methods and compare the performance of the two methods. Section 7 concludes the paper with an argument of the two methods and suggests some future research directions.

Nurgat, Elyas↗

NASA NTRS1995

On Lateral Viscosity Contrast in the Mantle and the Rheology of Low-Frequency Geodynamics

Mantle-wide heterogeneity is largely controlled by deeply penetrating thermal convective currents. These thermal currents are likely to produce significant lateral variation in rheology, and this can profoundly influence overall material behaviour. How thermally related lateral viscosity variations impact models of glacio-isostatic and tidal deformation is largely unknown. An important step towards model improvement is to quantify, or bound, the actual viscosity variations that characterize the mantle. Simple scaling of viscosity to shear-wave velocity fluctuations yields map-views of long- wavelength viscosity variation. These give a general quantitative description and aid in estimating the depth dependence of rheological heterogeneity throughout the mantle. The upper mantle is probably characterized by two to four orders of magnitude variation (peak-to-peak). Discrepant time-scales for rebounding Holocene shorelines of Hudson Bay and southern Iceland are consistent with this characterization. Results are given in terms of a local average viscosity ratio, (Delta)eta(bar)(sub i), of volumetric concentration, phi(sub i). For the upper mantle deeper than 340 km the following reasonable limits are estimated for (delta)eta(bar) approx. equal 10(exp -2): 0.01 less than or equal to phi less than or equal to 0.15. A spectrum of ratios (Delta)eta(bar)(sub i) less than 0.1 at concentration level eta(sub i) approx. equal 10(exp -6) - 10(exp -1) in the lower mantle implies a spectrum of shorter time-scale deformational response modes for second-degree spherical harmonic deformations of the Earth. Although highly uncertain, this spectrum of spatial variation allows a purely Maxwellian viscoelastic rheology simultaneously to explain all solid tidal dispersion phenomena and long-term rebound-related mantle viscosity. Composite theory of multiphase viscoelastic media is used to demonstrate this effect.

Ivins, Erik R.↗

NASA NTRS1992

Effects on post-glacial rebound from the hard rheology in the transition zone

The effect of an increase of viscosity in the earth-mantle's transition zone between 420 and 670 km on the geophysical signatures induced by postglacial rebound was analyzed using a self-gravitating model consisting of an elastic lithosphere, a three-layer viscoelastic mantle, and an inviscid core. Results show that the horizontal displacements and velocities as well as the stress pattern are extremely sensitive to the viscosity increase and to the chemical stratification of the transition zone. The findings indicate that geodetic data on intraplate deformation can be used to put constraints on the viscosity increase in the transition zone and on the amount of chemical stratification in the mantle.

Spada, Giorgio↗