10.08.2013 Views

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é ...

SHOW MORE
SHOW LESS

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>mi­maille 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 quasi­statique 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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!