|
Read Full Article Cited by
JOURNAL OF GEOPHYSICAL RESEARCH,
VOL. 99, NO. B8,
PAGES 15,919–15,938,
1994
Sensitivity of convection with an endothermic phase change to the form of governing equations, initial conditions, boundary
conditions, and equation of state
Joel Ita
Seismological Laboratory, California Institute of Technology, Pasadena
Scott D. King
Department of Earth and Atmospheric Sciences, Purdue University, West Lafayette, Indiana
Abstract
Recent convection calculations have demonstrated that an endothermic phase transition can greatly decrease the vertical flow
through the transition in a convecting system, in some cases leading to a layered flow. Using reasonable estimates of both
the Rayleigh number and Clapeyron slope of the spinel to perovskite plus magnesiowüstite phase change, these results suggest
that the 670-km phase change has a strong effect on mantle convection. This so-called “dynamic layering” phenomenon is further
investigated with a compressible finite element code using a two-dimensional, Cartesian geometry. We find a weak sensitivity
of the pattern of flow to the form of the equations, considering Boussinesq, extended Boussinesq, and anelastic compressible
forms of the governing equations, assuming that the thermodynamic properties (thermal expansivity, heat capacity, and latent
heat) remain constant. The pattern of flow, however, depends strongly on the initial conditions, boundary conditions and equation
of state. We compare the simple equation-of-state formulations used in previous work with a self-consistent equation of state
based on Debye and Birch-Murnaghan finite strain theory under a Mie-Grüneisen formulation. A thermal expansion coefficient
that decreases monotonically with depth and is unaffected by changes in phase or temperature greatly enhances dynamic layering.
This trend is reversed when the temperature, pressure, and phase dependence of thermodynamic properties such as thermal expansivity,
entropy, and heat capacity is introduced. At moderate Rayleigh numbers, the pattern of the flow is strongly influenced by
the pattern of the initial condition (i.e., the location of upwellings and downwellings); however, it is not sensitive to
the thickness of the initial thermal boundary layers. The sensitivity of the flow to the pattern of the initial condition
can potentially bias mass fluxes, especially for moderate Rayleigh number calculations.
Received 27
September
1993;
accepted 31
March
1994.
Read Full Article Cited by
Citation: Ita, J., and S. D. King
(1994),
Sensitivity of convection with an endothermic phase change to the form of governing equations, initial conditions, boundary
conditions, and equation of state,
J. Geophys. Res.,
99(B8),
15,919–15,938.
Copyright 1994 by the American Geophysical Union.
|