Débuts avec FreeFem++
Aide Mémoire pour les fenêtres de graphes :
f : lignes de niveau ou remplissage
l : lumière
v : affiche/masque les valeurs des lignes de niveau
m : affiche le mesh
w : sauvegarde le graphe
Exemple trouvé sur le site de Grudinin : http://grudinin.com/?p=508 cité en lien sur le site de FreeFem++
Pour lancer FreeFem++ : aller dans le dossier contenant le .edp puis taper FreeFem++ Truc.edp
Le fichier Mustang/FreeFem/Sphere.edp pour une sphère de rayon 36 µm, indice 1.46 et un nombre azimutal m=213 donne
le mode 0 et le mode 2 suivant : (comparé à ce qu'a calculé Yannick)
Longueur d'onde de la résonance : 1.4786 µm pour le mode 0 et 1.47387 µm pour le mode 2.
Yannick a obtenu 1.4738 µm, pas mal !
Par contre, il y a un mode 1 :
Le mode 3 est le suivant, à λ = 1.47184 µm
Remarque : Si on rétrécit la fenêtre sur la gauche (xl=radius-5 au lieu de xl=radius-10), le mode 0 à 1.4786 µm ne sort pas, on a
directement mode 0 à 1.47387 µm, le mode 1 est l'ancien mode 3, le mode 2 ressemble à l'ancien mode 1 (pourri) et on a un mode 3 à 1.46721 µm comme ceci :
En fait, le problème vient du nombre de points servant à la fabrication de la grille. En transformant la ligne 40 par
(cf Sphere2.edp) :
Th = buildmesh(W1r(Nm/2)+W1l(Nm)+W4(Nm)+W3l(Nm)+W3r(Nm/2)+W2(Nm)+S(Nm)) au lieu de
Th = buildmesh(W1r(Nm/4)+W1l(Nm/4)+W4(Nm)+W3l(Nm/4)+W3r(Nm/4)+W2(Nm)+S(Nm))
On retrouve les mêmes modes qu'avec la fenêtre plus petite.
Le programme est Mustang/FreeFem/16fev15/Disque.edp
Caractéristiques du disque (celles de l'article de Barclay) : r = 4 µm, e= 250 nm.
indice : 3.0564502060633147 celui du GaP calculé avec les calculs rapides à 1.522 µm
La fenêtre de calcul va de -3 à 0.5 en horizontal et -0.7 à 0.7 en vertical
Il y a 60 points sur la grille en horizontal et 40 en vertical.
En fixant m=23, on obtient ceci, pour le mode 2 : n=4, λ = 1.522 µm et mode 5 : n= 1, λ = 1.483 :
Quand on augmente la hauteur de la fenêtre de calcul (en augmentant aussi le nombre de points de la grille verticale) on se retrouve avec des
modes pour le moins exotiques :
la grille ressemble à ça avec 40 points en vertical (les calculs affichés juste au-dessus ont été faits avec 50):
