awslabs / palace

3D finite element solver for computational electromagnetics
https://awslabs.github.io/palace/dev
Apache License 2.0
224 stars 50 forks source link

Upgrade error estimation for time-dependent problems #233

Closed sebastiangrimberg closed 1 month ago

sebastiangrimberg commented 2 months ago

Consider both electric flux density and magnetic field recovery error for transient problems. Note: The improved recovery is incorrect for "infinitely thin" interior boundaries applying a PEC or impedance BC. In this case these BCs break the continuity of the recovered H and D fields. This will be resolved in a future PR and is a separate improvement to the changes here.

TODO

Note: Based on https://github.com/awslabs/palace/pull/237

sebastiangrimberg commented 2 months ago

Visualizations of the error indicator for the CPW example with lumped ports, at 60 GHz (to show the effect of multiple wavelengths):

Screenshot 2024-04-19 at 2 04 32 PM

Screenshot 2024-04-19 at 2 02 16 PM

Screenshot 2024-04-19 at 2 02 04 PM

We see that the correction for both energy terms (kinetic + potential, or magnetic + electric fields) makes the indicator more uniform over the length of the transmission line as we might expect, instead of preferentially indicating high errors at nodes of the electric field.