Dust in AstroBEAR - Update 2021/07/16

Progress:

For a while now I've been working on a new setup based on papers investigating dust destruction in so-called "windblown bubbles" around SNe. Most studies use several codes for this (one to simulate the immediate explosion aftermath to generate the shock structure, one to simulate the expansion of remnant for a few decades, one to simulate the expansion for several thousand years after & interaction with ISM). I'm trying to do the whole within one simulation. The most crucial part seems to be reproducing the correct shock speeds and structures. I think I'm pretty close now but will probably increase the initial velocity a bit to better match observations of Cas A. There are also some cosmetical issues (instabilities forming in the CSM where AMR levels change) plus issues at the boundaries. Will re-introduce dust once that's done.

Here's a gif of the current simulation going to almost 2k years: ​gif

And the individual outputs: ​1, ​2, ​3, ​4, ​5, ​6, ​7, ​8, ​9, ​10, ​11, ​12, ​13, ​14, ​15, ​16, ​17, ​18, ​19, ​20.

Two questions about this:

  • are the instabilities developing along the y-axis caused by numerics?
  • I'm getting string reflections from the centre which I'm not sure are physical, any way to suppress that?

Comments

No comments.