Search NASA⌕ Search

Engineering topics

Yuen, D. A.

Publications and source records attributed to Yuen, D. A..

36 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Mantle convection with internal heating and pressure-dependent thermal expansivity

Recent laboratory work suggests that the thermal expansivity alpha of the mantle decreases strongly with pressure. alpha determines the buoyancy of plumes and the rate at which plumes lose their thermal signature as they rise, and so may be expected to have a strong influence on the temperature and velocity structure of the mantle. Numerical simulations were conducted on convection in an internally heated, compressible mantle including constant and pressure-dependent alpha and thermal conductivity kappa. They show that a pressure-dependent alpha allows the existence of thermal plumes rising from the core-mantle boundary, however the plumes are weak and dissipate higher in the mantle.

Leitch, A. M.↗
NASA NTRS1991

Post-glacial relaxation of a viscously stratified compressible mantle

The postglacial relaxation phenomenon of a viscously compressible nonself-gravitating spherical-shell model is investigated. Analytical solutions to an exponentially depth-dependent viscosity are developed with various types of compressible density models such as those with an exponential dependence of the radius and an algebraic-root dependence of the radius. The solutions of a viscously compressible multilayered model with constant thermodynamic properties and viscosity are developed by the propagator matrix method. The results show that inferences of deep mantle viscosity from postglacial rebound would be hampered by mantle compressibility for long-wavelength harmonics because of the smaller excitation of compressible eigenfunctions in the lower mantle. The relaxation times and velocity fields are more sensitive to higher viscosity contrast for the exponentially varying viscosity than for models with discrete jumps in the viscosity structure.

Wu, J.↗
NASA NTRS1989

Internal heating and thermal constraints on the mantle

Numerical simulations have been carried out to study the average temperature and thermal structure of an internally heated mantle. It is found that in an incompressible mantle the hottest parts are localized in regions of slow flow in the middle and upper mantle, whereas in a compressible mantle they occur over laterally extensive regions near the core-mantle boundary (CMB). For chondritic concentrations of heat-producing elements (hpe's), the average thermal conductivity in the mantle must be high to avoid heating the core, and the temperature of the CMB is low to avoid large-scale melting. The mantle may have been extensively molten in the Archaean.

Leitch, A. M.↗
NASA NTRS1986

Excitation of the earth's rotational axis by recent glacial discharges

The effects of present-day glacial discharges and the growth of the Antarctic ice sheet on exciting the earth's rotational axis are studied. Glacial forcing could cause a maximum change in J2 of about one-third of the observed amount, for the Maxwell rheology and for Burgers' body models with a long-term, lower-mantle viscosity greater than about 10 to the 23rd P. For transient rheologies the amount of excitation due to glacial melting decreases. Polar wander is not much excited by recent glacial melting for the various types of rheologies examined.

Gasperini, P.↗
NASA NTRS1986

Adiabaticity and viscosity in deep mantle convection

A study has been conducted of steady convection with adiabatic and viscous heating for variable viscosity in the Boussinesq limit using the mean-field theory. A strong nonlinear coupling is found between the thermodynamic constants governing adiabatic heating and the rheological parameters. The range of rheological values for which adiabaticity would occur throughout the mantle has been established. Too large an activation volume, greater than 6 cu cm/mol for the cases examined, would produce unreasonably high temperature at the bottom of the mantle (greater than 6000 K) and superadiabatic gradients, especially in the lower mantle. Radiogenic heating plays a profound role in controlling dynamically mantle temperatures. Present values for the averaged mantle heat production would yield objectionably high temperatures in the lower mantle.

Quareni, F.↗
NASA NTRS1985

The role of a pressure-dependent rheology in the dynamics of mantle circulation

A thermomechanical model for upper mantle convection was constructed such that the thickness and the structure of the lithosphere are determined self consistently by the heat transported by convection. In this study of the interaction between the lithosphere and upper mantle, strongly temperature and pressure dependent rheologies for both Newtonian and non-Newtonian creep mechanisms are employed. For a strictly temperature dependent rheology an insignificant amount of heat, less than 12.5 mW/sq m, can be transported convectively for an interior viscosity, 0(10 sup 21 Pas), compatible with post glacial rebound. On the other hand, for similar values of the interior viscosity, steady heat fluxes between 20 and 40 mW/sq m are produced by introducing pressure dependence into the rheology. For the temperature and pressure dependent flow law the horizontally averaged interior temperature displays very little variation with the amount of heat evacuated, once all of the rheological parameters are fixed. This finding may have important ramifications for parameterized convection.

Yuen, D. A.↗
NASA NTRS1985

Surface deformation and geoid anomalies over single and double-layered convective systems

Using a primitive variable formulation of the finite-element method, the differences in the surface observables, such as topography and geoid, produced by single- and double-layered thermal convection, were compared. Both constant and depth-dependent viscosities have been considered. For the same Rayleigh number, larger surface perturbations are produced by single-cell convection. For the same Nusselt number, the magnitudes of the surface observables are greater for double-layered convection. For the same surface heat-flux, surface topographies have similar magnitudes, but the relative amplitudes of geoid anomalies depend greatly on the style of viscosity stratification. This difference in the geoid between the two systems increases with greater surface heat-flow, regardless of viscosity structure.

Koch, M.↗
NASA NTRS1985

Source phase shift - A new phenomenon in wave propagation due to anelasticity

The free oscillations of an anelastic earth model due to earthquakes were calculated directly by means of the correspondence principle from wave propagation theory. The formulation made it possible to find the source phase which is not predictable using first order perturbation theory. The predicted source phase was largest for toroidal modes with source components proportional to the radial strain scalar instead of the radial displacement scalar. The source phase increased in relation to the overtone number. In addition, large relative differences were found in the excitation modulus and the phase when the elastic excitation was small. The effect was sufficient to bias estimates of source properties and elastic structure.

Buland, R.↗
NASA NTRS1985

The effects of transient rheology on the interpretation of lower mantle viscosity

The role played by transient rheology in the interpretation of mantle viscosity is reexamined. The investigation has been carried out by comparing the amplitude responses with the data of secular variation of J(2), the relative sea-level histories at sites well within the ice margins and at the ice margin like the city of Boston. A linear Burgers body rheology has been assumed in ther lower mantle. The data near the edge of the ice load proves most sensitive to the transient viscosity structure. The non-monotonic behavior of sea-level data near Boston can be explained both by a steady-state lower mantle viscosity of 10 to the 22nd P with a thick lithosphere and by a transient lower mantle rheology but with a thin lithosphere. The long-term viscosity of the lower mantle in this second model has a steady-state value of around 5 x 10 to the 23rd P.

Sabadini, R.↗
NASA NTRS1985

Constraints on short-term mantle rheology from the J2 observation and the dispersion of the 18.6 y tidal Love number

Information derived from data recently acquired from the LAGEOS satellite is used to place some constraints on the rheological parameters of short-term mantle rheology. The validity of Lambeck and Nakiboglu's (1983) rheological model is assessed by formally developing an expression for the transformed shear modulus using a truncated retardation spectrum. This analytical formula is used to show that the parameters of the above mentioned model are not consistent at all with the amount of anelastic dispersion expected in the Chandler wobble and with the attenuation of seismic normal modes. The feasibility of a standard linear solid (SLS) rheology operating over intermediate timescales between 1 and 100 yr is investigated to determine whether the tidal dispersion at 18.6 yr can be explained by this model. An attempt is made to place some constraints on the parameters of the SLS model and the nature of short-term mantle rheology for timescales of less than 100 yr is discussed.

Sabadini, R.↗
NASA NTRS1985

Transient polar motions and the nature of the asthenosphere for short time scales

A uniformly valid mathematical formalism is developed to study the secular motions of the rotational axis of a layered viscoelastic earth due to seismic excitation. The changes required for implementing the formulation within the framework of the faulting problem. The rationale of adopting the chosen nrheological model, which contains a low-viscosity zone beneath the lithosphere and is based on linear Maxell constitutive relationship, is discussed. The impact of this low-viscosity channel on thhe two families of relaxation time, governing both isostatic readjustment and rotational processes, is considered. It is found that the polar motions depend sensitively on the viscosity structure of the asthenosphere and not at all on the underlying mantle. A gloal low-velocity zone with short-term asthenospheric viscosities less than about 5 x 10 to the 18th Pa-s and widths greater than 50 km is ruled out.

Boschi, E.↗
NASA NTRS1984

Secular rotational motions and the mechanical structure of a dynamical viscoelastic earth

A survey is presented of analytical methods for computing the linear responses of the rotational axis of a layered viscoelastic earth to surface loading. Theoretical research in this area is first summarized, and the differences between the mechanical boundary conditions to be applied at the interface separating the upper and lower mantles for an adiabatically and chemically stratified mantle are discussed. Some examples of polar wander and secular variation of the spin rate from glacial excitation are presented for various types of chemical and viscosity stratifications. The effects of an artificial density jump at the base of the lithosphere in models are examined, and certain issues concerning the fluid tidal Love number for different types of density stratification are addressed. The meaning of effective plate thickness over geological time scales for rotational dynamics is discussed.

Yuen, D. A.↗
NASA NTRS1978

Shear flow beneath oceanic plates - Local nonsimilarity boundary layers for olivine rheology

The principle of local similarity, which has been used to model the two-dimensional boundary layers in the oceanic upper mantle, permits calculation of the temperature, velocity, and stress fields with essentially analytic techniques. Finite difference numerical methods are hard pressed to resolve the detail required by the large variation of viscosity between the lithosphere and the asthenosphere. In this paper the local similarity approximation has been justified by quantitatively evaluating the effect of nonsimilarity due to viscous heating, nonlinear temperature- and pressure-dependent rheology, buoyancy, adiabatic cooling, etc. Nonsimilar effects produce only small modifications of the locally similar boundary layers; important geophysical observables such as surface heat flux and ocean floor topography are given to better than 10 percent by the locally similar solution. A posteriori evaluations of the terms neglected in the boundary layer simplification of the complete equations have been conducted on the locally similar temperature and velocity profiles close to the spreading ridge. The boundary layer models are valid to depths of 100 km at 3 m.y. and 10 km at 0.3 m.y.

Yuen, D. A.↗
NASA NTRS1976

Oceanic lithosphere and asthenosphere - Thermal and mechanical structure

A coupled thermomechanical subsolidus model of the oceanic lithosphere and asthenosphere is developed which includes vertical heat conduction, a temperature-dependent thermal conductivity, heat advection by a horizontal and vertical mass flow that depends on depth and age, contributions of viscous dissipation or shear heating, a linear or nonlinear deformation law relating shear stress and strain rate, as well as a temperature- and pressure-dependent viscosity. The model requires a constant horizontal velocity and temperature at the surface, but zero horizontal velocity and constant temperature at great depths. The depth- and age-dependent temperature, horizontal and vertical velocities, and viscosity structure of the lithosphere and asthenosphere are determined along with the age-dependent shear stress in those two zones. The ocean-floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of ocean-floor age; seismic velocity profiles which exhibit a marked low-velocity zone are constructed from the age-dependent geotherms and assumed values of the elastic parameters. It is found that simple boundary-layer cooling determines the thermal structure at young ages, while effects of viscous dissipation become more important at older ages.

Schubert, G.↗
NASA NTRS1976

Mantle plumes - A boundary layer approach for Newtonian and non-Newtonian temperature-dependent rheologies

Stress is placed on the temperature dependence of both a linear Newtonian rheology and a nonlinear olivine rheology in accounting for narrow mantle flow structures. The boundary-layer theory developed incorporates an arbitrary temperature-dependent power-law rheology for the medium, in order to facilitate the study of mantle plume dynamics under real conditions. Thermal, kinematic, and dynamic structures of mantle plumes are modelled by a two-dimensional natural-convection boundary layer rising in a fluid with a temperature-dependent power-law relationship between shear stress and strain rate. An analytic similarity solution is arrived at for upwelling adjacent to a vertical isothermal stress-free plane. Newtonian creep as a deformation mechanism, thermal anomalies resulting from chemical heterogeneity, the behavior of plumes in non-Newtonian (olivine) mantles, and differences in the dynamics of wet and dry olivine are discussed.

Yuen, D. A.↗
NASA NTRS1976

Oceanic lithosphere and asthenosphere: The thermal and mechanical structure

A coupled thermal and mechanical solid state model of the oceanic lithosphere and asthenosphere is presented. The model includes vertical conduction of heat with a temperature dependent thermal conductivity, horizontal and vertical advection of heat, viscous dissipation or shear heating, and linear or nonlinear deformation mechanisms with temperature and pressure dependent constitutive relations between shear stress and strain rate. A constant horizontal velocity u sub 0 and temperature t sub 0 at the surface and zero horizontal velocity and constant temperature t sub infinity at great depth are required. In addition to numerical values of the thermal and mechanical properties of the medium, only the values of u sub 0, t sub 0 and t sub infinity are specified. The model determines the depth and age dependent temperature horizontal and vertical velocity, and viscosity structures of the lithosphere and asthenosphere. In particular, ocean floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of the age of the ocean floor.

Schubert, G.↗
NASA NTRS1976

Thermal and mechanical structure of the upper mantle: A comparison between continental and oceanic models

Temperature, velocity, and viscosity profiles for coupled thermal and mechanical models of the upper mantle beneath continental shields and old ocean basins show that under the continents, both tectonic plates and the asthenosphere, are thicker than they are beneath the oceans. The minimum value of viscosity in the continental asthenosphere is about an order of magnitude larger than in the shear zone beneath oceans. The shear stress or drag underneath continental plates is also approximately an order of magnitude larger than the drag on oceanic plates. Effects of shear heating may account for flattening of ocean floor topography and heat flux in old ocean basins.

Froidevaux, C.↗