The overall efficiency and accuracy of standard finite element methods may be severely reduced if the solution of the boundary value problem entails singularities. In the particular case of time-harmonic Maxwell’s equations in nonconvex polygonal domains Ω, H1-conforming nodal finite element methods may even fail to converge to the physical solution. In this paper, we present a new nodal finite element adaptation for solving time-harmonic Maxwell’s equations with perfectly conducting electric boundary condition in general polygonal domains. The originality of the present algorithm lies in the use of explicit extraction formulas for the coefficients of the singularities to define an iterative procedure for the improvement of the finite element solutions. A priori error estimates in the energy norm and in the L2 norm show that the new algorithm exhibits the same convergence properties as it is known for problems with regular solutions in the Sobolev space H2Ω2 in convex and nonconvex domains without the use of graded mesh refinements or any other modification of the bilinear form or the finite element spaces. Numerical experiments that validate the theoretical results are presented.
Analysis and Application of Single Level, Multi-Level Monte Carlo and Quasi-Monte Carlo Finite Element Methods for Time-Dependent Maxwell's Equations with Random Inputs
