Computer Science
Permanent URI for this collection
Browse
Browsing Computer Science by Subject "2D Thermochemical Mantle Convection Model"
Now showing 1 - 1 of 1
Results Per Page
Sort Options
Item Manufacturing an Exact Solution for 2D Thermochemical Mantle Convection Models(Wiley Open Access [Commercial Publisher], American Geophysical Union [Society Publisher], Geochemical Society [Society Publisher], 2023) Trim, Sean; Butler, Samuel; McAdam, Shawn S.C.; Spiteri, RaymondIn this study, we manufacture an exact solution for a set of 2D thermochemical mantle convection problems. The derivation begins with the specification of a stream function corresponding to a non-stationary velocity field. The method of characteristics is then applied to determine an expression for composition consistent with the velocity field. The stream function formulation of the Stokes equation is then applied to solve for temperature. The derivation concludes with the application of the advection-diffusion equation for temperature to solve for the internal heating rate consistent with the velocity, composition, and temperature solutions. Due to the large number of terms, the internal heating rate is computed using Maple™, and code is also made available in Fortran and Python. Using the method of characteristics allows the compositional transport equation to be solved without the addition of diffusion or source terms. As a result, compositional interfaces remain sharp throughout time and space in the exact solution. The exact solution presented allows for precision testing of thermochemical convection codes for correctness and accuracy.