1D Planetary Models

I switched to 1D - and set the ambient density to 1e4 what it was. However I still got odd behavior. After looking at the hypre matrices and vectors, I discovered that there was a typo in the code…

MinTemp*TempScale

should have been

MinTemp/TempScale

If TempScale is 1, this does not matter, but when it is 1e6 it creates problems :)

In 1D there is no 'self-gravity' but the simulations are able to evolve on a dynamical time scale.

2D results

Attachments (3)

Comments

No comments.