The PML is an absorbing material that surrounds the computational cell and acts to attenuate the waves that propagate through it. Thus you should see the amplitude of the wave decrease as it moves through the PML (although the phase should remain the same, see the following review of PML for more information: http://math.mit.edu/~stevenj/18.369/pml.pdf ).
However your ctl file specifies a line source along the Y direction that overaps the PML, I wonder if this is what you mean by the plane wave having a 'bulge at the center' that is removed when you remove the PML in the Y direction. Ardavan Ajayram S K S M wrote: > /Hello!! > > /*I am trying to simulate the propogation of a plane wave in a glass > of refractive index 1.47. > > But, the png image showed distortion of the plane wave with a bulge at > the center. > > When I remove the pml layer in Y direction, I am able to get the plane > wave in the glass. But, by removing the pml we increase the errors in > our simulation. any solution to this problem? Why is the phase of the > wave getting affected due to the pml boundary? * > / > > CTL CODE: > / > (set! geometry-lattice (make lattice (size .8 .8 no-size))) > > (set! geometry (list > (make block (center 0 0) (size .8 .8 .8) > (material(make dielectric (epsilon 2.1609)))))) > > (set! sources (list > (make source > (src (make continuous-src (frequency 2.1777))) > (direction X) > (component Hz) > (center -0.3 0) > (size 0 .8 0)))) > (set! pml-layers (list (make pml (thickness 0.1))) > (set! resolution 500) > (run-until 10 > (at-beginning output-epsilon) > > (run-until 10 > (at-every .5 output-efield-x) > (at-every .5 output-efield-y)) > ;(at-every .5 output-dpwr)) > > * > Thanks for your time > > regards > Ajay > * > -- > Ajay Ram S K S M > MSc(Physics) > Sri Sathya Sai University > Prasanthi Nilayam > India. > ------------------------------------------------------------------------ > > _______________________________________________ > meep-discuss mailing list > [email protected] > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
