I am working on absorption by graphene plasmonic with array of graphene discs. I have used surface boundary condition for modeling graphene in comsol wave optic. the problem is that when i attach port to disc and insulator underneath it, the absorption general form is like the result of article but with many sharp ripples. I assign port to an air-box above disc but the simulation time get much much more that I can't get the results.