Quasi-static Maxwell's equations with a dissipative non-linear boundary condition: Full discretization

We study a time dependent eddy current equation for the magnetic field H accompanied with a nonlinear boundary condition, which is a generalization of the classical Silver–Müller condition for a non-perfect conductor. More exactly, the relation between the normal components of electric (E) and magnetic (H  ) fields obeys the following power law (linearized for small and large values) ν×E=ν×(|H×ν|^α−1H×ν)$\mathit{\nu }\times \mathit{E}=\mathit{\nu }\times ({|\mathit{H}\times \mathit{\nu }|}^{\alpha -1}\mathit{H}\times \mathit{\nu })$ for some α∈(0,1]$\alpha \in (0,1]$. We design a linear fully discrete approximation scheme to solve this nonlinear problem. The convergence of the approximations to a weak solution is proved, error estimates describing the dependence of the error on discretization parameters are derived as well. The efficiency of the proposed method is supported by numerical experiments.