I.4.2.4 Diffraction sur un profil NACA.

Le domaine de calcul est maillé en trois parties.

• Domaine 1. Autour du profil, ressemble à une couche limite, épaisse d'environ $3\%$ de la longueur de l'aile. Le paramètre de taille du maillage représente $1\%$ de la longueur de l'aile.
• Domaine 2. Ce domaine entoure la couche limite (domaine 1) dans un domaine limité, à peu près de même surface que le profil. Le paramètre de taille du maillage représente $5\%$ de la longueur de l'aile. Les domaines de calcul 1 et 2 sont les domaines dans lesquels nous représentons la carte des iso-valeurs des champs, figures I.4.10 à I.4.13.
• Domaine 3. C'est le domaine principal en ce qui concerne la surface. Limité par le domaine précédent et le cercle de diamètre trois fois la longueur du profil. Le paramètre de taille du maillage représente $20 \%$ de la longueur de l'aile. Soulignons que ce maillage est extrêmement grossier puisque à 1500 MHz il correspond à $h \approx 2\lambda$, qui est très loin de la loi empirique (I.0.3) entre le pas de discrétisation et la fréquence.

Nous obtenons 2976 éléments et 1615 nuds. Le nombre d'arêtes libres extérieures (frontière absorbante artificielle) est de 48, le nombre d'arêtes libres intérieures est de 202. Pour l'aile NACA, nous donnons une représentation du champ u autour du profil. Ceci est réalisé à l'aide d'une grille plus fine que celle du maillage de discrétisation. Nous calculons u_h à l'aide de la formule (I.2.27) section I.2.1.5. Le maillage utilisé pour la représentation du champ est constitué des domaines 1 et 2 décrits ci-dessus, mais avec un paramètre de taille du maillage uniforme qui vaut $1\%$ de la longueur de l'aile. Nous obtenons 3572 éléments et 2010 nuds. Le nombre d'arêtes libres extérieures (frontière graphique artificielle) est de 246, le nombre d'arêtes libres intérieures est de 202. Les valeurs des parties réelles ou imaginaires du champ diffracté à $1 \, 500$ MHz sont données par les courbes de niveau figures I.4.10, I.4.11, I.4.12 et I.4.13. L'onde incidente aborde le profil en incidence rasante, soit par le bord de fuite (figures I.4.10 et I.4.11), soit par le bord d'attaque (figures I.4.12 et I.4.13).

Figure I.4.10: Bord de fuite, Dirichlet TM, $\Im {(u_D)}$.

Figure I.4.11: Bord de fuite, Neumann TE, $\Im {(u_D)}$.

Figure I.4.12: Bord d'attaque, Dirichlet TM, $\Re {(u_D)}$.

Figure I.4.13: Bord d'attaque, Neumann TE, $\Re {(u_D)}$.

Nous obtenons, figures I.4.14 et I.4.15, les calculs de SER comparés à ceux du code SHF2D en balayage bistatique. La figure I.4.15 présente la comparaison sous la forme d'un diagramme d'antenne, c'est-à-dire en coordonnées polaires : ce genre de représentation permet de voir facilement l'énergie diffusée selon la direction d'observation. Le lecteur notera que, à 1500 MHz, nos calculs donnent de légères oscillations autour des valeurs calculées par SHF2D dans le mode TM: ceci provient d'une contribution de la CLA, contribution parasite faible, qui est plus apparente dans le mode TM.

Figure I.4.14: NACA 0012 SER, UWVF-IE comparaison à 1500 MHz.

Bord de fuite,	bord d'attaque.

Figure I.4.15: Diagramme d'antenne NACA 0012 SER, UWVF-IE comparaison à 150 MHz.

Bord de fuite,	bord d'attaque.