TY - JOUR
T1 - A high-order finite-difference method applied to large rayleigh number mantle convection
AU - Larsen, Tine B.
AU - Yuen, David A.
AU - Moser, Jiří
AU - Fornberg, Bengt
PY - 1997
Y1 - 1997
N2 - A variable-grid, high-order finite difference (FD) method is applied to the modeling of mantle convection in both two- and three-dimensional geometries. The algorithm combines extreme simplicity in programming with a very high degree of accuracy. Memory requirements are low and grow almost linearly with the total number of grid points in three dimensions, regardless of the increase in grid points in the vertical direction. Higher-order methods, such as eighth order, yield significantly better results than a second-order method for the same grid size, with only a modest increase in memory requirements. This is particularly important for high Rayleigh number convection, where the large number of grid points required to obtain an accurate enough solution with second-order schemes would make the computation extremely costly. The small-scale features in the hard-turbulent regime under high-Rayleigh number situations can greatly stress low-order methods, and in these situations a high-order method is definitely needed. We have numerically simulated three-dimensional time-dependent convection for constant properties up to Ra = 108, using an eighth-order FD scheme. Both purely base-heated and partially internally heated situations have been considered. The hot plumes are broader near the surface with internal heating. Detailed studies of the three-dimensional constant viscosity plumes indicate that no small-scale circulation takes place in the ascending plume heads regardless of the heating configuration in accordance with predictions from boundary layer theory.
AB - A variable-grid, high-order finite difference (FD) method is applied to the modeling of mantle convection in both two- and three-dimensional geometries. The algorithm combines extreme simplicity in programming with a very high degree of accuracy. Memory requirements are low and grow almost linearly with the total number of grid points in three dimensions, regardless of the increase in grid points in the vertical direction. Higher-order methods, such as eighth order, yield significantly better results than a second-order method for the same grid size, with only a modest increase in memory requirements. This is particularly important for high Rayleigh number convection, where the large number of grid points required to obtain an accurate enough solution with second-order schemes would make the computation extremely costly. The small-scale features in the hard-turbulent regime under high-Rayleigh number situations can greatly stress low-order methods, and in these situations a high-order method is definitely needed. We have numerically simulated three-dimensional time-dependent convection for constant properties up to Ra = 108, using an eighth-order FD scheme. Both purely base-heated and partially internally heated situations have been considered. The hot plumes are broader near the surface with internal heating. Detailed studies of the three-dimensional constant viscosity plumes indicate that no small-scale circulation takes place in the ascending plume heads regardless of the heating configuration in accordance with predictions from boundary layer theory.
KW - High Rayleigh number convection
KW - Higher-order finite-difference methods
KW - Three-dimensional convection
UR - http://www.scopus.com/inward/record.url?scp=0031437478&partnerID=8YFLogxK
U2 - 10.1080/03091929708208973
DO - 10.1080/03091929708208973
M3 - Article
AN - SCOPUS:0031437478
SN - 0309-1929
VL - 84
SP - 53
EP - 83
JO - Geophysical and Astrophysical Fluid Dynamics
JF - Geophysical and Astrophysical Fluid Dynamics
IS - 1-2
ER -