I am trying to simulate a Bragg Fiber Grating in 2D (for now). I am using Wave Optics Module in comsol. My fiber for testing, has 20 periods. I am using ELECTROMAGNETIC WAVES, FREQUENCY DOMAIN (EWFD), adding a Scattering Boundary Condition where the wave will enter with a Wave given by E field X = 0, Y = 0 and Z = 1V/m. Other Scattering Boundary Condition was added where I want to read the wave that is transmitted. Also was added as Perfect Magnetic Conductor the Top and Bottom Exterior Boundaries.
In the Study, I add only the step 1: Frequency Domain with the range of frequency like range(0.5*f0,deltaf,1.5*f0), Where deltaf is 1THz and f0 is equal to c_const/1.45[um].
When I try to se the results using the formula: intop2(ewfd.nPoav)/intop1(0.5*E0^2/Z0_const*n0)
I got a weird graph.
What am I missing?