Etude théorique de radars géologiques - Epublications - Université ...
Etude théorique de radars géologiques - Epublications - Université ...
Etude théorique de radars géologiques - Epublications - Université ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Samuel BESSE : Étu<strong>de</strong> Théorique <strong>de</strong> Radars Géologiques Analyses <strong>de</strong> sols, d'antennes et interprétation <strong>de</strong>s signaux 117<br />
L'intégration <strong>de</strong> l'équation précé<strong>de</strong>nte suivant ρ et la condition limite sur Ez à la surface du<br />
conducteur (Ez(ρ=a)=0) permet d'obtenir une expression <strong>de</strong> Ez fonction <strong>de</strong> ρ (Ez est indépendant <strong>de</strong><br />
θ et <strong>de</strong> z respectivement par symétrie <strong>de</strong> rotation et <strong>de</strong> translation selon l'axe z). On pose également<br />
L l'inductance linéique associée au fil.<br />
∂ I<br />
E z = L ∂ t v2 ∂Q<br />
∂ z <br />
L = <br />
ln 2 a <br />
(52.a)<br />
(52.b)<br />
Pour obtenir l'équivalence d'un champ moyen au sein <strong>de</strong> la cellule et éliminer la dépendance<br />
en 1/ρ, l'équation (52) doit être moyennée sur la surface d'une cellule ΔxΔy. On remarquera que<br />
cette équation (52) permet la prise en compte <strong>de</strong>s interactions entre le courant et le champ<br />
électrique.<br />
∂ I<br />
E z = L moy ∂ t v2 ∂Q<br />
∂ z <br />
L moy = <br />
2<br />
∫ ln x y , a<br />
<br />
a cos d d <br />
x y<br />
Par ailleurs, l'équation <strong>de</strong> conservation <strong>de</strong> la charge impose :<br />
(53.a)<br />
(53.b)<br />
∂ I ∂Q<br />
= 0 (54)<br />
∂ z ∂ t<br />
Finalement, les équations (53) et (54) constituent un système d'équations liées qui peut se<br />
résoudre via une métho<strong>de</strong> numérique telle que la FDTD. Les <strong>de</strong>nsités <strong>de</strong> charges Q et le courant I<br />
sont décalés spatialement d'une <strong>de</strong>mimaille et sont calculés à <strong>de</strong>s instants différents. Le fil seul se<br />
comporte comme un circuit série avec une self et une capacité mais le formalisme <strong>de</strong> Holland offre<br />
la possibilité d'ajouter <strong>de</strong>s circuits localisés tel que générateur <strong>de</strong> tension, résistances, selfs et<br />
capacités (en modifiant simplement les valeurs obtenues pour le fil seul). Il est même possible <strong>de</strong><br />
coupler la FDTD avec un co<strong>de</strong> <strong>de</strong> C.A.O. <strong>de</strong> type Spice pour introduire <strong>de</strong>s circuits complexes ou<br />
<strong>de</strong>s éléments non linéaires [76][77].<br />
Pour gérer l'extrémité d'un fil non connecté et la jonction entre un fil et un conducteur parfait<br />
on applique respectivement I=0 et Q=0.<br />
Il existe à ce jour <strong>de</strong>s formulations plus complexes qui permettent <strong>de</strong> prendre en compte<br />
plusieurs fils dans une même cellule, <strong>de</strong>s fils situés ailleurs que sur les arêtes, <strong>de</strong>s fils obliques dans<br />
un milieu homogène [73], <strong>de</strong>s jonctions entre plusieurs fils... On peut également citer une étu<strong>de</strong><br />
effectuée à l'INRIA [71] qui propose un modèle basé sur une approximation quasistatique du<br />
champ électrique au voisinage d'un fil. Cette étu<strong>de</strong> est complétée par la détermination <strong>de</strong><br />
l'inductance optimale qui est un paramètre clé <strong>de</strong> la métho<strong>de</strong>. Une formule explicite est donnée dans