Modélisation du spectro-photomètre lyxs EN SPECTRAL avec Starlyx
Simuler la diffusion de Mie
Important: avant de traiter ce TP voir le sous-TP "TP1".
Ici nous simulons une détection spectrale de lyxs. Les capteurs (les sphères intégrantes) ont plusieurs canaux de détection pour les longueurs d'onde:
415 445 480 515 555 590 630 680 910 nm
Dans le fichier scene.xml, l´attribut "FILTERS" des capteurs a changé:
Par exemple pour la sphère intégrante BS:
<sensor
FILTERS="415 445 480 515 555 590 630 680 910 "
VOLUME="World"
ALBEDO="0.0"
NAME="sphereBS"
FILE="sphere_2_trous.obj"
ANGLE="180.0"
>
<translation X="-67" Y="0" Z="0" ></translation>
</sensor>
La source émet la même puissance pour chacune de ces longueurs d'onde:
<spectrum
DATA="415 1 445 1 480 1 515 1 555 1 590 1 630 1 680 1 910 1 "
NAME="spectrumSource"
>
</spectrum>
Les paramètres de transfert radiatif sont déterminés dans le noeud xml "mie1".
<Mie
NAME="mie1"
D_UM="1"
NR="1.5"
NI="0"
PHI="0.001"
>
</Mie>
D_UM : diamètre des particules sphériques en um.
NR: indice de réfraction, partie réelle
NI: indice de réfraction, partie imaginaire
PHI: fraction volumique
Starlyx détermine les coefficients ks,k et g pour les longueurs d'onde du système à partir de la théorie de Mie.
Exercice
Configurer l'environnement:
cd /home/lafrier/environnement_tp_starlyx/tp_starlyx
source profile
Se mettre dans le dossier de ce tp:
cd tp_starlyx/TP2
Exécuter starlyx:
starlyx2 -x scene.xml
On peut voir le résultat brut avec la commande:
more sensors.csv
On peut utiliser gnuplot pour afficher les spectres:
gnuplot plot.gnuplot
Modifier le diamètre des particules D_UM par exemple à 2 um. Lancer le calcul et afficher les nouveaux spectres.
