A numerical model for steady-state temperature distributions in solid-medium high-pressure cell assemblies

We present a numerical model for calculating the temperature distribution inside resistance-heated high-pressure solid-medium axi-symmetric cell assemblies that incorporates both composition- and temperature-dependent thermal conductivity. The code was validated using both analytic solutions of simplified thermal diffusion problems and comparisons to actual laboratory experiments and was found to be reliable in matching the temperature characteristics of multi-anvil experiments. Calculations for various cell assembly designs resulted in temperature fields that are consistent with experimental measurements of thermal gradients. These calculations also illustrated the influence of temperature-dependence of thermal conductivity, all important and often-overlooked property, oil the thermal profiles. This model may be used to fine-tune the design of cell assemblies, either to minimize thermal gradients or to produce a desired temperature distribution. The four "typical" multi-anvil cells that we used to demonstrate this technique have temperature profiles across the sample that range from 25 to 75 degrees C when the thermocouple temperature is 1200 degrees C. The thermocouple in all four is in a region where the temperature gradient is oil the order of 100 degrees C per millimeter, which could lead to experimental temperature uncertainties that are cot-related with the thermocouple location.