pyro2
pyro2 copied to clipboard
[WIP] implement one-sided stencils for 4th order
these are for symmetry boundaries.
this addresses #244
Here's the rt_multimode problem without the initial perturbation + reflecting BCs run with the 4th order solver, showing Mach number in a vertical slice:
blue is the current main and orange is this PR.
We see that this PR reduces the velocity noise near the boundary.