parthenon-hpc-lab / parthenon

Parthenon AMR infrastructure
https://parthenon-hpc-lab.github.io/parthenon/
Other
105 stars 33 forks source link

Add a vector wave equation example #1091

Closed lroberts36 closed 3 weeks ago

lroberts36 commented 1 month ago

PR Summary

This PR adds a vector wave equation to the fine_field example. That is, we solve the equations

\partial_t \vec{C} + \nabla \times \vec{D} = 0
\partial_t \vec{D} - \nabla \times \vec{C} = 0
\nabla \cdot \vec{C} = 0
\nabla \cdot \vec{D} = 0

using an upwind constrained transport scheme (currently just using piecewise constant reconstruction). The face centered fields are evolved and prolongation is done using the scheme of Toth & Roe.

TODO:

PR Checklist

Yurlungur commented 1 month ago

Include Toth & Roe prolongation so that divergence constraint is preserved in AMR

Is this needed if using edge fields?

lroberts36 commented 1 month ago

Include Toth & Roe prolongation so that divergence constraint is preserved in AMR

Is this needed if using edge fields?

@Yurlungur: No, it shouldn't be. I had started writing it using face fields though.

lroberts36 commented 3 weeks ago

@par-hermes format