Migrating HSE from 2D to 3D

This morning I ran a simulation of the HydroStatic Star module in 3D. The simulation was successful and stable.
A remarkable fact is the change in the pressure profile that leads to a more natural temperature profile. Here are the two models compared, they have an identical density profile but differ in pressure profiles and therefore temperature profiles. This is most likely due to the switch from plummersoftening2D to its 3D version.
Edit (Jonathan explained the reason for the change): in 2D the point gravity object is actually a line charge, so the force falls with 1/r rate. in 3D it is 1/r2

I ran the 3D model with 4 AMR levels, and this led to less mass loss and increased stability overall. Computational time was not too long (~2 hours on 32 cores). ​movie

As I did with our 2D simulation, I compared the values of this new 3D model with the refence model, I was able to load all the values in the same plot this time for a better understanding. Legend:

  • continuos lines: reference model
  • dots: our model
  • blue/purple: pressure
  • red/orange: temperature
  • gray/green: density


Overall, this looks like a good result to me. The next step is to put our accretion disk back in the picture.

Comments

No comments.