Surface Impedance Boundary Conditions

Does anyone have tried Surface Impedance Boundary Conditions (SIBC) on NGSolve?
I want to do the induction heating analysis on NGsolve, but the frequency is high and I need it. I might be able to do it also with NGXFEM, but SIBC is easier.

you simply add a term like

bfa += sqrt(1j*omega*sigma/mu) * u.Trace()*v.Trace()*ds("SIBCboundary") 

to your complex bilinear-from.

Joachim

1 Like