Numerical studies of thermal convection with temperatureand pressure-dependent viscosity at extreme viscosity contrasts

Tania S. Khaleque, A. C. Fowler, P. D. Howell, M. Vynnycky

Research output: Contribution to journalArticlepeer-review

Abstract

Motivated by convection of planetary mantles, we consider a mathematical model for Rayleigh-Bénard convection in a basally heated layer of a fluid whose viscosity depends strongly on temperature and pressure, defined in an Arrhenius form. The model is solved numerically for extremely large viscosity variations across a unit aspect ratio cell, and steady solutions for temperature, isotherms, and streamlines are obtained. To improve the efficiency of numerical computation, we introduce a modified viscosity law with a low temperature cutoff. We demonstrate that this simplification results in markedly improved numerical convergence without compromising accuracy. Continued numerical experiments suggest that narrow cells are preferred at extreme viscosity contrasts, and this conclusion is supported by a linear stability analysis.

Original languageEnglish
Article number076603
Pages (from-to)-
JournalPhysics of Fluids
Volume27
Issue number7
DOIs
Publication statusPublished - Jul 2015

Fingerprint

Dive into the research topics of 'Numerical studies of thermal convection with temperatureand pressure-dependent viscosity at extreme viscosity contrasts'. Together they form a unique fingerprint.

Cite this