Skip to content

Conversation

@mfournier01
Copy link
Collaborator

@mfournier01 mfournier01 commented Apr 15, 2025

The main modifications are:

  • Modification of the existing code to support 2D geometry. Quantities along the z-axis are only updated if there's more than 1 cell along that direction.
  • Development of a second method which uses face-centered velocity based on Riemann fluxes. This method only supports VL2 integrator for the moment (see Eq. 7 from https://ui.adsabs.harvard.edu/abs/2008ASPC..385..273L/abstract).
  • Splitting of the AdvectTracer routine into two separate methods. Method 0 corresponds to the already implemented velocity field interpolation method. Method 1 is the new flux-based method. The method can be specified in the input file within the <tracer> section, using the method parameter.

A comparison of the two methods is presented below, and can be compared to Fig. 9 from https://ui.adsabs.harvard.edu/abs/2013MNRAS.435.1426G/abstract.

hydro_density.mp4

A 2D test using the Kelvin-Helmholtz and method 1 is also presented below.

hydro_density.mp4

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

5 participants