ENSM.SE , 3MI, 158, cours Fauriel, 42023 Saint Etienne cedex 2
Olivier Gipouloux
Institut Camille Jordan UMR 5208, Batiment ISTIL, 15 Boulevard Latarjet, 69622 Villeurbanne Cedex
Olivier.Gipouloux@cdcsp.univ-lyon1.fr
Le but de ce projet est de
fournir une mesure quantitative des incertitudes affectant le comportement d’un
stockage de déchets nucléaires qui est imparfaitement connu. Plus précisément
on veut calculer la probabilité que la(les) concentration(s) en radionucléides
accessibles à l’homme dépassent un certain seuil.
Pour faire ce calcul on dispose
d’un modèle simplifié d’écoulement et de transport en milieu poreux, ainsi que
d’une géométrie 2D simplifiée. Le code de simulation numérique est CASTEM
développé au CEA.
Le calcul de la probabilité qui nous intéresse pourrait se faire en utilisant une méthode de Monte Carlo en réalisant un grand nombre de simulations. Mais devant le temps de calcul important pour chaque simulation, on a choisi une méthode alternative, i.e. les plans d’expériences numériques.
Les résultats présentés ici concernent un nombre de facteurs (paramètres d’entrée du code numérique) égal à six (différentes perméabilités, porosités et coefficients de diffusion) et une seule réponse (sortie du code numérique) représentée par la concentration du radionucléide à un endroit fixé, appelé exutoire.
Le plan utilisé est un plan de type hyper cube latin à 30 niveaux. La surface de réponse calculée à partir de ces 30 expériences est une surface obtenue par krigeage, avec une fonction de corrélation de type exponentielle.
Une fois la surface de réponse calculée, on calcule la probabilité d’intérêt par Monte Carlo.