COMSOL-nanoparticles-PML-scattering formalism in RF module

[pusk.lt]

MOEMS membrane as a pressure sensor

Dear team, minds :)

The thing is that I have micrcro membrane with optical gratings and I need to have it with optical signal registration: i.e. some pressure to the membrane is applied and using laser light interference of diffraction we need to model it To callibrate diffraction minimums and maximums to pressure units.

P.S. this should be done in COMSOL

Sincerely

PHD . Karolis

 

[anupam g]

hi guys - i am simulating scattering by gold nano-sphere in comsol 4.2. I am encountering following problem:-

Problem - In the far-field analysis : When I tried to calculate the scattering cross section using the near field the results are incorrect but for the same program the absorption cross section calculated using the near field match perfectly.

I am trying to simulate the scattering in the same way as it is explained in mie scattering 3.5 example.
some info :-
radius of : scatterer -100nm,auxiliary boundary=200nm,inner pml=250nm,outer pml=300nm. Its the aux boundary where Efar is being defined.

any help will be appreciated...i wonder if someone has encountered same problem because in 3.5 it works well but the same formulation is giving incorrect result as explained above but absorption cross section is fine.

thanks.

 

[anupam g]

hi guys - i am simulating scattering by gold nano-sphere in comsol 4.2. I am encountering following problem:-

Problem - In the far-field analysis : When I tried to calculate the scattering cross section using the near field the results are incorrect but for the same program the absorption cross section calculated using the near field match perfectly.

I am trying to simulate the scattering in the same way as it is explained in mie scattering 3.5 example.
some info :-
radius of : scatterer -100nm,auxiliary boundary=200nm,inner pml=250nm,outer pml=300nm. Its the aux boundary where Efar is being defined.

any help will be appreciated...i wonder if someone has encountered same problem because in 3.5 it works well but the same formulation is giving incorrect result as explained above but absorption cross section is fine.

thanks.

 

[zhoujingj]

Hello Srbasket,

Thank you so much for sharing such useful information about how to solve a scattering problem by Comsol! The website http://srdjancomsol.weebly.com is also very helpful.

Actually, I followed the 2-step procedure as you provided to solve my case which is light scattering from a waveguide. The waveguide made of polymer (refractive index 1.58) has a square cross-section (1 micrometer times 1 micrometer) and it is supported by a layer of the same polymer with the thickness of 0.3 micrometer. Beneath the polymer layer, there is SiO2 substrate (refractive index 1.45). In this case, The effect of the 2-layer substrate on the background excitation field should be taken into account. However, it is not straightforward to get the analytical expression of the (Fresnel formulation) background field in the presence of the polymer layer and the SiO2 substrate. On the other hand, it is either incorrect to simply put a plane wave onto the whole structure including the scatter and the supporting layers since the supporting layers are then limited in lateral size and it also causes scattering.

Fortunately, I found your message online and tried to follow your directions. First, I calculated the background field with the supporting layers (by setting periodic boundary conditions to ensure the infinite lateral size) but without the scatter. Second, I stored the solution (Solve Manager/Store Solution) of the first step and then used it as an initial value (Solve Manager/Initial Values/Stored Solution) for the second-step calculation. However, I found some problem because the initial value is not regarded as the excitation field by Comsol so that this 2-step calculation just gives the same result as simply putting a plane excitation wave, exp(-j*k0_rfweh*y), onto the scatter and the supporting layers which also serve as scatters due to its limited lateral size.

I also tried total field calculation, but the limited lateral size of the substrate is always the problem. Vice verse, if I eliminate the problem of the limited lateral size of the substrate by setting periodic boundary conditions, the solution is according to an array of scatters instead of a single one.

Would you please give me some more suggestions?

Thank you!
Jing Zhou
jizhou@umich.edu

One clarification:
In 3) second time you are solving press RESTART, not the SOLVE

Hi, all,
folks keep contacting me about different things regarding this subject. I have advanced a bit, but still some issues remaining.
1) In excitation field definition you can forget about using complicated formulae with different expressions for z>0 and z<0. It looks like that if u use k_rfw, COMSOL will automatically consider refractive index of material into k vector in different media.
2) THERE is a big difference if u are solving for SCATTERD or TOTAL field. You should read Jiaming Jin's book to find out more. In short, if u are solving for SCATTERD field your excitation should be defined in the most precize analytical way, using Fresnel coefficients, and in that case you need complicated formulae with z>0 and z<0, since in SCATTERING MODULE you have to define field present everywhere except PML. In total field calculation, you have to insert incident wave on one boundary (port or scattering boundary condition with excitation), and propagation should be calculated... There are some advantages/disadvantages in both cases, depending what you whant to calculatet, how much memory you have, and some FEM issues that you can read in Jim's book...
3) About 2 step procedure: When u create new model name your dependents variables to scEx1,etc, and Application mode name to rfw1. Then you can make your geometry. After, you go to Multyphysics/Module navigator and ADD new scEx2,...,rfw2, and now you have two models on the same geometry. In both models you define separately Subdomain setting, Boundary settings, etc. Meshing is the same, since you have one geometry, and 2 models defioned on the same geometry. So, first you solve situation when your particle is made of air let's say, and later use that solution as excitation in your second model. Thus, in Physics/Scalar Variables you set Eoix... with Fresnel coefficients (or not if u don't want), and for E0ix2 you set Ex, E0iy2 set Ey, etc... First step is to SOLVER MANAGER/Solve for you highlight rfw1, and for the output rfw1 (not necessary), and when it is done, you go to SOLVER MANAGER/Initial value, there you first Store solution, and check STORED SOLUTION option as initial values, SOLVE FOR you highlight rfw2, and the same for OUTPUT... That should work. Theoretical benefits of this procedure are not so clear to me. IF SOMEONE COMPARES IT; PLEASE PUT YOUR RESULTS AND CONCLUSION HERE...
4) Someone asked me "How did you separate the scattered field from the substrate and from the nanoparticle?" Maybe this 2 step procedure is the answer to that. Maybe you can normalized what you get for particle air and partice made of metal case, I don't know exactely.
5)Scattering cross-section: Read RF module reference guide. There is a Stratton-Chu formula implemented ion COMSOL for calculating far-fileld components. Onece you have them, you need spherical surface (it doesn't have to be the same as one for Stratton-Chu, can be smaller to save you integratin time) and one should do boundary integration of normEfar*normEfar, and that is proportional to scattering cross section. If you want to be more precise calculate this expresion: (normEfar*normEfar)/(Surface of particle*Ein*Ein*R*R) where R is the radius of your integrating spherical surface.
6) IF ANYONE HAVE GOOD RECIPE FOR PML THAT WORKS FINE, please POST IT!

Cheers