-
Filter by Topic
Most Popular
All Topics
- List all discussions
Facebook
Twitter
Delicious
Digg
Stumbleupon
More services...How to set epsilon/mu parameters in PML(4.0 RF)
|
Thread index | Previous thread | Next thread | Start a new discussion |
RSS feed |
Turn on email notifications |
10 Replies
Last post: September 23, 2010 3:38pm UTC
July 29, 2010 8:26pm UTC
How to set epsilon/mu parameters in PML(4.0 RF)
Hi,
I tried to simulate a TM plane wave by setting background electric field Ey=exp(-j*emw.k0*x).
PML is surrounded by the rectangle simulatd area.
When I set parameters epsilon/mu/conductivity as air (epsilon=1, mu=1,conductivity=0)all over the geometry(PML+simulated area), I can get pefect plane wave. However, when I set parameters as other material(epsilon=2, mu=1,conductivity=0) all over the geometry,I can not get pefect plane wave.
BTW, the thickness of PML is one wavelength.
The attached images are the results of these two different settings.
My question is:
in order to get perfect plane wave, how to set PML parameters if the simulated area is not air?
Reply | Reply with Quote | Send private message | Report Abuse
August 1, 2010 11:49am UTC in response to Jie Xiang
Re: How to set epsilon/mu parameters in PML(4.0 RF)
Hi
for me it looks like your lateral sides are not "Perfectly" matched and you have reflexion from mthe lateral sides interacting, but this is just a guess.
In which version are you ? use the new forum tags, it helps us others udnerstand better too ;)
--
Good luck
Ivar
Reply | Reply with Quote | Send private message | Report Abuse
September 8, 2010 3:25pm UTC in response to Ivar Kjelberg
Re: How to set epsilon/mu parameters in PML(4.0 RF)
I am experiencing the same issue in a scattering model using Comsol 4.0.0.937. The particle is located at the origin, then we have some space and the outer subdomain is a PML boundary (See attached file pml_air.png).
If I use air in the surrounding and the PML subdomains, the solution looks nice. If I switch both permittivities to some value unequal 1, I get severe errors in the solution. (See pml_glass.png). I also tried setting the permittivity of the PML region equal to 1 and the surrounding of the particle to some value unequal 1. Same wrong results.
Seems like these are UPML (UnPerfectly Matched Layers ;o)). Unfortunately, the documentation is not of much help here. Has anyone experienced the same issue or does anyone have an idea?
Best regards!
Reply | Reply with Quote | Send private message | Report Abuse
September 10, 2010 5:39am UTC in response to Christian Matyssek
Re: How to set epsilon/mu parameters in PML(4.0 RF)
For the expression of incident wave in 'Scalar Variables', the wavevector k0 should be corrected as sqrt(epsilon)*k0, in which the epsilon is the permittivity of surrounding. My COMSOL version is 3.5a, and I did not know if this solution would work in COMSOL4.0. Good luck!
Reply | Reply with Quote | Send private message | Report Abuse
September 10, 2010 9:20am UTC in response to Kaiting He
Re: How to set epsilon/mu parameters in PML(4.0 RF)
Hi
an check that you have latest patch, fo me today this means 4.0.0.993, see tha main COMSOL web site
--
Good luck
Ivar
Reply | Reply with Quote | Send private message | Report Abuse
September 15, 2010 9:35am UTC in response to Ivar Kjelberg
Re: How to set epsilon/mu parameters in PML(4.0 RF)
Hi!
Thanks for hinting me at the wave vector definition. sqrt(epsilon)*k0 indeed did the trick. This is somehow inconvenient, because when changing the material of the surrounding domain, you'll also have to change the incident field definition. The only solution I see so far is to adding a global variable with the refractive index.
By the way: Is there a complete variable reference for emw.*? Haven't been able to find one so far.
Thanks and best regards!
Christian
Reply | Reply with Quote | Send private message | Report Abuse
September 15, 2010 4:49pm UTC in response to Christian Matyssek
Re: How to set epsilon/mu parameters in PML(4.0 RF)
I also have not seen variable of this kind. :)
He Kaiting
Reply | Reply with Quote | Send private message | Report Abuse
September 16, 2010 6:44am UTC in response to Christian Matyssek
Re: How to set epsilon/mu parameters in PML(4.0 RF)
Hi
but epsilon and "n" the refractive index are linked to the spped of light no ?
--
Good luck
Ivar
Reply | Reply with Quote | Send private message | Report Abuse
September 16, 2010 7:15am UTC in response to Ivar Kjelberg
Re: How to set epsilon/mu parameters in PML(4.0 RF)
Hi!
In a medium with refractive index n the speed of light is c=c0/n, where c0 is the _vacuum_ speed of light. In terms of the wave vector, in a medium with refractive index n, the wave vector is k=k0*n. The refractive index and the permittivity epsilon are related by sqrt(epsilon)=n.
So in the Comsol model, we use the free space wave vector emw.k0 to specify the wave. This is somewhat confusing, because I would have expected Comsol to scale the wave vector automatically in a medium with n!=1, but it seems like it does not.
Reply | Reply with Quote | Send private message | Report Abuse
September 16, 2010 9:21am UTC in response to Christian Matyssek
Re: How to set epsilon/mu parameters in PML(4.0 RF)
Hi
ideally we should have this confirmed by COMSOL support, no ?
(PS I'n not by COMSOL ;)
--
Good luck
Ivar
Reply | Reply with Quote | Send private message | Report Abuse
September 23, 2010 3:38pm UTC in response to Ivar Kjelberg
Re: How to set epsilon/mu parameters in PML(4.0 RF)
Yes, you are right. I had this issue in 3.5a already and the support told me to rescale the wave vector. I'll contact the support.
Reply | Reply with Quote | Send private message | Report Abuse
Rules and guidelines
