Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

electromagntic propagation through three layers

Please login with a confirmed email address before reporting spam

hello
I am student on electronic and I am trying to simulate with comsol on rf module , in summarize my work is to simulate the sun light propagation through three layers bonded to each other and I had used the floquet periodic condition and this witch include this equations :

Parameter

lambda0= 0.5e-6[m]
k0_emwh = 2*pi/lambda0
alpha1 =30[deg]
n1 1.97235
k1x= k0_emwh*real(n1)*cos(alpha1)
k1y =k0_emwh*real(n1)*sin(alpha1)
freq c/lambda0
c= 3*10^8[m/s]
n2 0.125
alpha2= asin(sin(alpha1)*real(n1/n2))
k2x= k0_emwh*real(n2)*cos(alpha2)
K2y= k0_emwh*real(n2)*sin(alpha2)
n3 1.97235
alpha3 =asin(sin(alpha3)*n2/real(n3))
E0z= exp(-i*k1x*x)
Ez1= E0z[V/m]
k3x= k0_emwh*real(n3)*cos(alpha3)
k3y= k0_emwh*real(n3)*sin(alpha3)

and I put the boundary condition above and below the layers
Periodic condition 1
K vector for folquet periodicity :
k1x
k1y
Port 1
In the front of the layer
Put the wave excitation in this port (on)
Pin=1
Téta(in) = 0
E0= Ez1
Beta =k1y
Port 2
Put the wave excitation in this port (off)
E0= Ez1
Beta =k1y

in this simulation i would like to study the propagation of the electromagnetic wave through three micrometer thickness of material, and each material have property , and my problem when reduce the thickness the electromagnetic wave propagate with difficulty and when i increase the thickness the wave propagate easily, i don't know if my simulation is true.

Thank you for your help



15 Replies Last Post 7 mar 2013, 03:32 GMT-5
Sergei Yushanov Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 18 gen 2013, 11:22 GMT-5
Belbey,

Several notes about your model:

1) Expressions such as E0z= exp(-i*k1x*x) should be defined under Model>Definitions>Variables but not under Global Definitions>Parameters as in your model, since E0x depends on spatial coordinates. That is why your parameter is highlighted as red.

2) Propagation constant in the port conditions should be normal component of the wavector, not parallel. That is, beta=k1x and not beta=k1y as in your model.

3) I assume you want to analyze TE polarization. In that case, incident wave should be defined as E={0,0, exp(-i*k1y*y) }, not E={0,0,exp(-i*k1x*x)} as in your model.

4) Consider creating mesh with several elements through the thickness of the thin middle layer using Mapped mesh.

5) If you are not sure about your results, compare Comsol results with analytical solution for multilayer structures available in textbooks. Numerical and analytical solutions should be in perfect agreement for all incidence angles for 0 to 90 degrees.

Regards,
Sergei
Belbey, Several notes about your model: 1) Expressions such as E0z= exp(-i*k1x*x) should be defined under Model>Definitions>Variables but not under Global Definitions>Parameters as in your model, since E0x depends on spatial coordinates. That is why your parameter is highlighted as red. 2) Propagation constant in the port conditions should be normal component of the wavector, not parallel. That is, beta=k1x and not beta=k1y as in your model. 3) I assume you want to analyze TE polarization. In that case, incident wave should be defined as E={0,0, exp(-i*k1y*y) }, not E={0,0,exp(-i*k1x*x)} as in your model. 4) Consider creating mesh with several elements through the thickness of the thin middle layer using Mapped mesh. 5) If you are not sure about your results, compare Comsol results with analytical solution for multilayer structures available in textbooks. Numerical and analytical solutions should be in perfect agreement for all incidence angles for 0 to 90 degrees. Regards, Sergei

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 19 gen 2013, 02:50 GMT-5
thank you sergai

I'll try i'll try to change what are you had tell , just I would like tell you i use the floquet periodic condition, i suppose it's the same , i ll give you the results that i found
thank you sergai I'll try i'll try to change what are you had tell , just I would like tell you i use the floquet periodic condition, i suppose it's the same , i ll give you the results that i found

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 19 gen 2013, 06:29 GMT-5
hello sergai

I redid my simulation with your modification but it's the same results , and in this simulation i want to study the electromagnetic wave propagation through x axes , can you help me please .

thank you
hello sergai I redid my simulation with your modification but it's the same results , and in this simulation i want to study the electromagnetic wave propagation through x axes , can you help me please . thank you

Sergei Yushanov Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 21 gen 2013, 12:56 GMT-5
Belbey,

I don’t have 4.1, so I am attaching report of your model created with later version. Comparison with analytical solution for the case of TM polarization is also included.

I set the conductivity of the first ZnO layer to zero because it does not make much sense to specify evanescent wave at the input port. If you want to study wave attenuation in the first ZnO layer, then you might want to model structure with four layers: first layer is lossless dielectric (say, air) where you specify incident wave, second layer is ZnO with specified thickness, etc.

Regards,
Sergei
Belbey, I don’t have 4.1, so I am attaching report of your model created with later version. Comparison with analytical solution for the case of TM polarization is also included. I set the conductivity of the first ZnO layer to zero because it does not make much sense to specify evanescent wave at the input port. If you want to study wave attenuation in the first ZnO layer, then you might want to model structure with four layers: first layer is lossless dielectric (say, air) where you specify incident wave, second layer is ZnO with specified thickness, etc. Regards, Sergei


Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 27 gen 2013, 06:21 GMT-5
thank you sergai for your answer , i will trying this and i give the news
thank you sergai for your answer , i will trying this and i give the news

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 27 gen 2013, 13:58 GMT-5
Could you please attach the mph file; thank you!
Could you please attach the mph file; thank you!

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 27 gen 2013, 14:31 GMT-5
of course
of course

Sergei Yushanov Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 28 gen 2013, 12:56 GMT-5
SR,

Attached is model file. Please note that this is 4.3a.

Regards,
Sergei
SR, Attached is model file. Please note that this is 4.3a. Regards, Sergei


Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 29 gen 2013, 04:09 GMT-5
thank you sergai

I don't have 4.3 A version , but I am trying to work under 4.1 version with you last exemple , when I finish I post the result
thank you sergai I don't have 4.3 A version , but I am trying to work under 4.1 version with you last exemple , when I finish I post the result

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 29 gen 2013, 08:42 GMT-5
hello ,


I have a small problem with parameter equation when i try to enter this equation in parameter case

sqrt(k0^2*em - kay^2 - j*eps^3)

x-component of wavevector of middle layer (evanescent square root)

the comsol don't accept this , and I would like to know what the
eps

thank you
hello , I have a small problem with parameter equation when i try to enter this equation in parameter case sqrt(k0^2*em - kay^2 - j*eps^3) x-component of wavevector of middle layer (evanescent square root) the comsol don't accept this , and I would like to know what the eps thank you

Sergei Yushanov Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 29 gen 2013, 09:41 GMT-5
Belbey,

eps is small number. Under Global Parameters, define a=eps and you will see that it's numerical value is 2.22e-16.
If you define eps1=1e-12, it will be OK.

Sergei
Belbey, eps is small number. Under Global Parameters, define a=eps and you will see that it's numerical value is 2.22e-16. If you define eps1=1e-12, it will be OK. Sergei

Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 29 gen 2013, 14:03 GMT-5
Hi

indeed "eps" is a predefined number in COMSOL (as in other numerical tools) it's the smalles real number one can distinguish from "1.0" in the binary representation of real numbers. It's rather FPU dependent, so you can forget to put it to any positive power it will wanish.
Even, as most physics equations are 2nd order PDE's then sqrt(eps) represents, the absolute best, relative ratio, or number of significant digits you could expect from your results

--
Good luck
Ivar
Hi indeed "eps" is a predefined number in COMSOL (as in other numerical tools) it's the smalles real number one can distinguish from "1.0" in the binary representation of real numbers. It's rather FPU dependent, so you can forget to put it to any positive power it will wanish. Even, as most physics equations are 2nd order PDE's then sqrt(eps) represents, the absolute best, relative ratio, or number of significant digits you could expect from your results -- Good luck Ivar

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 30 gen 2013, 05:31 GMT-5
Hi sergai

this is what I've done but it's still wrong , but the notation in certain parameter indicate with red color , please if some one can correct me ,i am always trying to find a solution , in parallel if i find something i poste

thank you
Hi sergai this is what I've done but it's still wrong , but the notation in certain parameter indicate with red color , please if some one can correct me ,i am always trying to find a solution , in parallel if i find something i poste thank you

Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 30 gen 2013, 13:04 GMT-5
Hi

red means normally typos, or undefined variable name:s so those should be tracked and corrected,
orange is unit inconsistencies, thise might be ignored, with the risk of scaling errors

--
Good luck
Ivar
Hi red means normally typos, or undefined variable name:s so those should be tracked and corrected, orange is unit inconsistencies, thise might be ignored, with the risk of scaling errors -- Good luck Ivar

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 7 mar 2013, 03:32 GMT-5
Hi ,

just I finished my simulation through the last exemple , but in last part from the example i don't understand how to obtain the diagram refraction in 1 D plot groupe can help you me to have the diagrme
thank you
Hi , just I finished my simulation through the last exemple , but in last part from the example i don't understand how to obtain the diagram refraction in 1 D plot groupe can help you me to have the diagrme thank you

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.