3-D treatment of convective flow in the earth's mantle

Conference · Tue May 01 04:00:00 EDT 1984

OSTI ID:6237474

Baumgardner, J R

A three-dimensional finite-element method is used to investigate thermal convection in the earth's mantle. The equations of motion are solved implicitly by means of a fast multigrid technique. The computational mesh for the spherical problem is derived from the regular icosahedron. The calculation described use a mesh with 43,554 nodes and 81,920 elements and were run on a Cray X. The earth's mantle is modeled as a thick spherical shell with isothermal, free-slip boundaries. The infinite Prandtl number problem is formulated in terms of pressure, density, absolute temperature, and velocity and assumes an isotropic Newtonian rheology. Solutions are obtained for Rayleigh numbers up to approximately 10/sup 6/ for a variety of modes of heating. Cases initialized with a temperature distribution with warmer temperatures beneath speading ridges and cooler temperatures beneath present subduction zones yield whole-mantle convection solutions with surface velocities that correlate well with currently observed plate velocities. 8 references, 6 figures.

• Research Organization: Los Alamos National Lab., NM (USA)
• DOE Contract Number: W-7405-ENG-36
• OSTI ID: 6237474
• Report Number(s): LA-UR-84-3762; CONF-8405240-2; ON: DE85003734
• Country of Publication: United States
• Language: English

Related Subjects

58 GEOSCIENCES
580201* -- Geophysics-- Seismology & Tectonics-- (1980-1989)
CONVECTION
DIFFERENTIAL EQUATIONS
EARTH MANTLE
ENERGY TRANSFER
EQUATIONS
EQUATIONS OF MOTION
FINITE ELEMENT METHOD
FLUID FLOW
HEAT TRANSFER
MASS TRANSFER
NUMERICAL SOLUTION
PARTIAL DIFFERENTIAL EQUATIONS
THREE-DIMENSIONAL CALCULATIONS