Re: [Meep-discuss] question on effect of pml layer on the phase of the EM wave

[Ardavan F. Oskooi]

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
> meep-discuss@ab-initio.mit.edu
> http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss


___
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

[Meep-discuss] question on effect of pml layer on the phase of the EM wave

[Ajayram S K S M]

*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
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss