search

FormingWorlds / PROTEUS

Coupled atmosphere-interior framework to simulate the temporal evolution of rocky planets.
https://fwl-proteus.readthedocs.io
Apache License 2.0
11 stars 1 forks source link

Rock vapour calculation #215

Open nichollsh opened 4 hours ago

nichollsh commented 4 hours ago

Currently, PROTEUS only supports volatile outgassing, but at higher temperatures we should expect significant amounts of rock vapour to be produced. Integrating these into the radiative transfer would be complex, but they could quite easily be incorporated into calculations of MMW and bulk density.

We can calculate vapour partial (surface) pressures at a given temperature and fO2 using codes such as Vaporock and Lavatmos. These codes are somewhat annoying to setup, since they depend on ThermoEngine, but I have been successful in doing this on my Mac laptop.

We would have to think carefully about how to include these gases alongside the currently supported volatiles.

timlichtenberg commented 4 hours ago

Agree, we had similar issues in Björn Koops' BSc thesis earlier this year when trying out Vaporock. At the meeting on Wednesday, we should discuss whether pre-computed look-up tables may provide an entry.

Why do you think that integrating them into the radiative transfer would be complex? If we have the opacities for the major gases, SOCRATES should be able to compute k-tables?

nichollsh commented 4 hours ago

Pre-computing this is a nice idea, particularly if we are assuming BSE composition since that would only leave 2 variables (Tsurf, fO2) for the table's axes. We would probably have to store the data in units of surface density (kg/m^2), rather than pressure, so that it's agnostic to the surface gravity.

I think it might be complex because SOCRATES doesn't currently support many refractory species. It uses hardcoded gas IDs defined in gas_list_pcf.f90, so we would need to manually add all of the required gases to this list. This is easy to do, but it would mean diverging from the Met Office SOCRATES distribution somewhat.

nichollsh commented 4 hours ago

The currently supported refractory species are: TiO, VO, Na, K, FeH, CrH, Li, Rb, Cs.

timlichtenberg commented 4 hours ago

Mh, yeah, these are clearly not sufficient. But yes, perhaps assuming BSE and condensing the axes would make sense. I wonder if the axes should be (Tsurf, Psurf, fO2), but that's a thing to discuss on Wednesday with Yamila's group further.