LE GPR MAX

GprMax a été initialement développé en 1996, lorsque la modélisation numérique utilisant la

méthode FDTD et, en général, la modélisation numérique du GPR étaient à leurs
balbutiements. Depuis lors, un certain nombre de solveurs commerciaux, et d'autres solutions
librement disponibles FDTD sont devenus disponibles, mais gprMax est resté l'un des outils
de simulation les plus utilisés dans la communauté GPR. Il a été utilisé avec succès pour un
large éventail d'applications dans les universités et l'industrie, et a été cité plus de 200 fois
depuis 2005.
Comprendre comment les ondes électromagnétiques se propagent dans des environnements
hétérogènes naturels ou créés par l'homme est un problème difficile. Par conséquent,
l'interprétation des données acquises à l'aide de GPR est souvent difficile en raison des
interactions complexes entre le système GPR, la ou les cibles d'intérêt et l'environnement.
Cela est particulièrement évident lorsque vous essayez d'interpréter des informations
quantitatives à partir de données GPR. Une interprétation réussie des données GPR repose
généralement sur une expérience considérable acquise grâce à une expérimentation
approfondie, et même alors se limite souvent à identifier les domaines d'intérêt ou les
anomalies dans les données. Pour faire progresser notre compréhension du GPR et fournir un
moyen de tester de nouvelles techniques de traitement des données et algorithmes
d'interprétation, il est important de disposer d'un logiciel de simulation précis et robuste.
GprMax est un logiciel open source qui simule la propagation des ondes électromagnétiques.
Il résout les équations de Maxwell en 3D en utilisant la méthode FDTD (Finite-Difference
Time-Domain). GprMax a été conçu pour modéliser le radar de pénétration du sol (GPR) mais
peut également être utilisé pour modéliser la propagation des ondes électromagnétiques pour
de nombreuses autres applications.
GprMax possède de nombreuses fonctionnalités puissantes et avancées telles que:
 Modélisation des matériaux anisotropes
 Modélisation de matériaux dispersifs (avec des formulations multipolaires Debye,
Lorenz ou Drude)
 Modélisation de sols aux propriétés diélectriques et géométriques réalistes
 Construction d'objets hétérogènes
 Construction d'objets à surfaces rugueuses
 Bibliothèques intégrées de modèles d'antennes
Fondements des méthodes temporelles
Les fondements des méthodes temporelles sont exposés dans cette section par une description
générale et par une présentation des méthodes de résolution temporelles recensées lors de la
revue de la littérature. Les difficultés de mise en œuvre, les avantages et inconvénients par
rapport aux méthodes fréquentielles sont ensuite indiqués, et un inventaire de logiciels de
calcul basés sur des méthodes temporelles est dressé.
Un système dynamique est décrit par des équations gouvernantes découlant des relations
d'équilibre. Ces équations gouvernantes résultent généralement en des équations aux dérivées
partielles dans lesquelles certains termes peuvent varier en fonction du temps.
Les méthodes dites temporelles sont des méthodes qui permettent d'analyser directement dans
le temps un phénomène dynamique, à partir de ses équations gouvernantes, de ses conditions
aux frontières et de ses conditions initiales. La résolution du problème se fait en fonction de la
variable temporelle, pas à pas, suivant un schéma temporel donné, c'est-à-dire que l'état du
système à un instant donné est évalué à partir de l'état du système déjà évalué aux instants
précédents.
Pour chaque pas de temps, la solution est calculée simultanément pour l’ensemble des points
du domaine considéré. Pour les problématiques de bruit et de vibrations, les méthodes
temporelles sont utilisées pour résoudre l'équation de propagation d'ondes acoustiques et
l'équation l astrodynamique, dans le temps. Si ces méthodes s'affranchissent positivement des
hypothèses simplificatrices propres aux méthodes fréquentielles, deux critères sont toutefois
incontournables pour assurer leur utilisation bénéfique et avantageuse :
Disposer d'un modèle mathématique exact pour décrire le problème et disposer de la capacité
de résolution du modèle, analytiquement, ou plutôt numériquement dans la plupart des cas.
Les méthodes de résolution temporelles

 Les méthodes Finite Difference Time Domain (FDTD)

L’approche FDTD de la solution numérique des équations de Maxwell consiste à
discrétiser à la fois l'espace et le temps continuent. Ainsi la discrétisation spatiale ∆x,
∆y et ∆z et les pas temporels ∆t jouent un rôle très important - car plus ils sont petits,
plus le modèle FDTD est proche d’une réelle représentation du problème. Cependant,
les valeurs des étapes de discrétisation doivent toujours être finies, car les ordinateurs
ont une quantité limitée de stockage et une vitesse de traitement finie.
Par conséquent, le modèle FDTD représente une version discrétisée du problème réel
et de taille limitée. Le batiment bloc de cette grille FDTD discrétisée est la cellule Yee
nommée d'après Kane Yee qui a été le pionnier de la Méthode FDTD. Ceci est illustré
pour le cas 3D de la figure.
 Figure : La cellule 3D FDTD Yee
Dans un milieu continu, isotrope et linéaire, les composantes des champs électrique et
magnétique, des équations de Maxwell à 3D suivant (x, y, z), s'écrivent alors :

Naturellement, plus nous utilisons un grand nombre de valeurs de F dans les cellules autour
de la position étudiée, plus l'évaluation de la dérivée sera précise. La précision et la stabilité
de l'algorithme dépendent bien sûr de l'ordre des dérivées : un algorithme d'ordre est plus
stable et plus précis mais il nécessite de connaître à chaque instant plusieurs valeurs
précédentes de la fonction F qu'il faut donc stocker dans la mémoire de l'ordinateur (d'où des
problèmes lors de l'utilisation de modèles de grandes dimensions) et il utilise un plus grand
nombre d'opérations

Figure : schéma numérique 3D de la maille élémentaire selon Yee.

Il y’a aussi d’autres méthodes comme :

 Les méthodes Time Domain Finite Element Method (TD-FEM)

 Les méthodes Time Domain Boundary Element Method (TD-BEM)
 Les méthodes TD-BEM/FEM
 Les méthodes Local Boundary Integral Equations (LBIE)