Pytesimal Metal INTrusion model
A Python-powered 3D thermal evolution model of the pallasite formation region in the palalsite parent body. Designed to be used with parent-body background mantle temperatures from Pytesimal [1].
PytesiMINT uses the semi-implicit fraction-step method to apply Crank-Nicolson to 3D [2, 3]. Spatial diffusivity is incorporated and the heat capacity method tracks latent heat during crystallisation of the metal.
While PytesiMINT is designed to model the formation of pallasite meteorites, specifically the cooling following injection of molten metal into the solid olivine mantle of the parent body, it can be applied to any conductively-cooling intrusion problem, especially those where spatially varying diffusivity between the intrusive material and country rock is relevant.
The following web resources were useful in the initial development stages of this project:
- The Crank-Nicolson method implemented from scratch in Python, Georg Walther (2013). Archived version, 2023-03-27.
- Python implementation of Crank-Nicolson scheme, Claudio Bellei (2019). Archived version, 2023-03-27.
- Practical Numerical Methods with Python, L.A. Barba, C.D. Cooper, G.F. Forsyth (2014). Archived version, 2023-03-27.
-
Murphy Quinlan, M., Walker, A. M., Davies, C. J., Mound, J. E., Müller, T., & Harvey, J. (2021). The conductive cooling of planetesimals with temperature-dependent properties. Journal of Geophysical Research: Planets, 126, e2020JE006726. https://doi.org/10.1029/2020JE006726
-
W. Cen, R. Hoppe, N. Gu. Fast and accurate determination of 3D temperature distribution using fraction-step semi-implicit method AIP Adv., 6 (2016), Article 095305, https://doi.org/10.1063/1.4962665
-
S. Sahijpal, Thermal evolution of non-spherical asteroids in the early solar system, Icarus, Volume 362, (2021). https://doi.org/10.1016/j.icarus.2021.114439.