OpenSEMBA / dgtd

Maxwell's curl equations solver using discontinuous Galerkin methods.
BSD 3-Clause "New" or "Revised" License
15 stars 2 forks source link

Upwind Flux #14

Closed AlejandroMunozManterola closed 1 year ago

AlejandroMunozManterola commented 1 year ago

While Centered Flux is enough to define simple problems, a modification of said flux, known as Upwind Flux, allows for the definition of more complex problems while increasing the general stability of the model.

An Upwind Flux, based on the Strong Form described by Tim Warburton and Jan S. Hesthaven in their book Nodal Discontinuous Galerkin Methods, is the chosen one to be implemented in our model.

This Upwind Flux implementation, based on jump operators, requires from a customised Discontinuous Galerkin Trace Integrator to work.

Required for #10 to work.

AlejandroMunozManterola commented 1 year ago

The addition of specific FiniteElementOperator methods focused on the implementation of the Upwind Flux parameters, along with a modification of the Mult method in our MaxwellEvolution : TimeDependentOperator custom classes, now allows for the use of Upwind Flux for the problem definition.