search

OpenSEMBA / dgtd

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

Maxwell DG Trace Integrator #13

Closed AlejandroMunozManterola closed 1 year ago

AlejandroMunozManterola commented 1 year ago

A customised DG Trace Integrator is required in order to create Discontinuous Galerkin flux matrices for communication of the elements between their interfaces.

The already available DG Trace Integrators existing in MFEM do not work, due to the nature of the expressions used.

A DG Trace Integrator based on the Jump Operator where the evaluated bilinear form is b * < n · u, w > is required.

AlejandroMunozManterola commented 1 year ago

A functional DG Trace Integrator has been implemented as an extension to the capabilities of MFEM.

The MaxwellDGTraceIntegrator takes into account up to two directions for multi-normal evaluated parameters, which exist in three dimensional problems, and a beta coefficient which can be used to scale the output matrix.

Like its MFEM counterparts, MaxwellDGTraceIntegrator can be used in the AddBdrFaceIntegrator method of a BilinearForm with a BoundaryMarker to only apply the Integrator in specific degrees of freedom or faces, this allows for an easy application of boundary conditions or specific modifications.