Simulation numérique d'écoulements en milieu poreux avec l'équation de Richards

La théorie des écoulements saturés et non saturés en milieu poreux repose sur l'équation de Richards, une équation différentielle partielle non linéaire. L'équation de Richards peut être résolue efficacement avec un schéma itératif de type Picard adapté de manière à être conservatif. Une approche numérique par éléments finis permet alors de résoudre l'équation dans des domaines où la répartition des propriétés géologiques est complexe. Le milieu naturel contient souvent des conditions à la limite de type infiltration ou suintement qui ne sont pas des réflexions directes des conditions à la limite mathématiques pour une équation du type Richards. Une méthode itérative permettant de simuler ce type de conditions est présentée puis validée avec des mesures sur un modèle physique.