Dust in AstroBEAR - Update 2020/03/16

Comment:

I wasn't able to make blog posts for the last couple of weeks so this is a summary of the progress made since then.

Objectives:

  • Finish active advection routines
  • Extend interpolation options to dust
  • Continue dust processing routines

Progress:

Interpolation

  • I've extended the interpolation scheme to include dust. So PCM, PLM, and PPM are now also available for the dust.

Dust Advection

  • I reduced the Cloud Crushing setup to just a clump moving through gas as I wasn't sure whether there's an issue with the boundaries, the wind, the fluxes, etc.
  • Found out that way that the x- and y-fluxes kept getting messed up → fixed in F_prim and included and added an HLL solver for the dust to both the HLL and HLLC gas solver.

Below are videos of the tests I ran to check the advection:

Simulations with the HLL Gas Solver (and HLL Dust Solver)

Simulations with the HLLC Gas Solver (and HLL Dust Solver)

This looks like it's working alright considering the dust is implemented as a pressure-less fluid. Once the dust and the gas are coupled via drag, the dust should start getting a velocity in y as well.

Next steps:

  • Now working on turning on the drag routines I have in the code (they still need to be properly integrated into the source routine)

Comments

No comments.