Étude de schémas numériques pour les écoulements diphasiques en milieu poreux déformable pour des maillages quelconques : application au stockage de déchets radioactifs, Study of numerical schemes for two-phase flow in porous media for any meshes : application to storage of nuclear waste

Ophélie Angelini

Description

Sous la direction de Robert Eymard
Thèse soutenue le 10 novembre 2010: Paris Est
Les écoulements diphasiques en milieu poreux sont des phénomènes complexes et qui concernent de nombreux problèmes industriels. EDF travaille sur la faisabilité et la sécurité d'un stockage en couche géologique profonde de déchets nucléaires. Dans ce domaine la simulation des écoulements diphasiques en milieu poreux est particulièrement importante dans au moins trois domaines : tout d'abord lors de la phase de ventilation des galeries du stockage qui pourrait désaturer la roche présente et ainsi en modifier ses propriétés de rétention, mais également lors de la phase de resaturation des matériaux et enfin lors de l'arrivée de l'eau sur les parties métalliques contenues dans le stockage qui entraînera alors des phénomènes de corrosion et un dégagement d'hydrogène. Dans ce contexte, EDF souhaite se doter de méthodes numériques performantes et robustes ne nécessitant pas de conditions restrictives sur la forme des mailles. Ce travail s'inscrivant dans cette problématique, est consacré dans un premier temps au développement du schéma volumes finis SUSHI (Scheme Using Stabilization and Hybrid Interfaces) dans le code de mécanique d'EDF, Code_Aster afin de modéliser les écoulements diphasique en milieu poreux. Ce schéma a été développé en 2D et en 3D. Parallèlement une nouvelle formulation qui permet de traiter de manière uniforme les écoulements en milieu saturé et insaturé pour des problèmes miscibles et immiscibles est proposée. Différentes études modélisant des difficultés liées aux problématiques du stockage de déchets radioactifs en couches géologiques profondes ont été traitées. On peut citer l'étude d'un bi-matériau qui met en avant le ré-équilibrage capillaire d'un matériau par un autre possédant des propriétés et des conditions initiales en saturation très hétérogènes. On citera également l'étude de l'injection d'hydrogène dans un milieu initialement saturé en eau pure qui est tirée du benchmark « Ecoulement diphasique » proposé par le GNR MOMAS. Cette étude avait pour objectif de mettre en évidence le bon traitement de l'apparition d'une phase dans un milieu saturé et donc la pertinence de notre nouvelle formulation à traiter d'une manière unifié un problème d'écoulement saturé et un problème d'écoulement insaturé
-Schémas numériques
-Volumes finis
-Milieu poreux
-Écoulements diphasiques
The two-phase flow in porous media is a complex phenomenon and which relate to many industrial problems. EDF works on the feasibility and the safety of a storage in deep geologic layer of nuclear waste. In this domain the simulation of the two-phase flow in porous media is particularly important in at least three domains : first of all during the phase of ventilation of the galleries of the storage which could desaturate the rock and so modify its properties, but also during the phase of resaturation of the materials and finally during the arrival of the water on the metal parts contained in the storage which will then involve phenomena of corrosion and a hydrogen release. In this context, EDF wishes to obtain robust numerical methods without restrictive condition on the mesh. This work is dedicated at first to the development of the finite volume scheme SUSHI (Scheme Using Stabilization and Hybrid Interfaces) in the code of mechanics of EDF, Code_Aster in order to simulate the two-phase flow in porous media. This scheme was developed in 2D and in 3D. At the same time a new formulation which allows to simulate in a uniform way the flows in saturated and unsaturated porous media for miscible and immiscible problems is proposed. Various studies simulating difficulties related to the problems of the storage of nuclear waste in deep geological layers were study. We can quote the study of a bi-material which advances the capillary rebalancing of a material by an other one possessing properties and initial very heterogeneous conditions in saturation. We will also quote the study of the injection of hydrogen in an porous media initially saturated in pure water which is proposed by the benchmark two-phase Flow proposed by the GNR MOMAS. This study had for objective to bring to light the good treatment of the appearance of a phase in a saturated porous media and thus the relevance of our new formulation to study with a way unified a problem of saturated flow and a problem of unsaturated flow
-Numerical scheme
-Finite volume
-Porous medium
-Two-phase flow
Source: http://www.theses.fr/2010PEST1004/document

Sujets

• Méthode des volumes finis
• Porosité

Informations

• Publié par: Thesee
• Nombre de lectures: 54
• Langue: Français
• Poids de l'ouvrage: 10 Mo

Extrait

THÈSE
présentéepourobtenirlegradede
DOCTEUR DE L’UNIVERSITÉ PARIS-EST
Mention MATHÉMATIQUES APPLIQUÉES
par
OphélieANGELINI
intitulée
Étudedeschémasnumériques
pourlesécoulementsdiphasiques
enmilieuporeuxdéformable
pourdesmaillagesquelconques
Applicationaustockagededéchetsradioactifs
ÉcoleDoctorale: MSTIC
Encadrantsindustriels: ClémentCHAVANT
SylvieGRANET
Soutenuele10novembre2010devantlejurycomposéde
Rapporteurs: M.PascalOMNES CEA
M.RobertCHARLIER UniversitédeLiège
Examinateurs: M.RolandMASSON IFP
M.GuyBONNET UniversitéParis Est
M.ÉricCHÉNIER
Directeurdethèse: M.RobertEYMARD UniversitéParis Est
Invités: MmeSylvieGRANET EDFR&D
M.ClémentCHAVANT EDFR&D
tel-00587364, version 1 - 20 Apr 2011tel-00587364, version 1 - 20 Apr 2011Remerciements
Je tiens à remercier tout d’abord le département d’Analyse Mécanique et Accoustique d’EDF
R&D de m’avoir offert l’opportunité d’effectuer cette thèse dans des conditions idéales. Merci tout
particulièrementàFrancoisWaeckel,chefdudépartement,pourl’intérêtqu’ilasuporteràmontravail
et également à Pierre Badel et Erwann Galenne, chefs du groupe T64, pour la confiance qu’ils m’ont
accordée.
Mesremerciementss’adressentensuiteauxdifférentespersonnesquiontencadrémathèse.Mercià
RobertEymard,directeurdethèse,dem’avoirfaitpartagersesconnaissancesscientifiques.Jeremercie
ensuitemonco dirdethèse,EricChenier,pourlesnombreuxéchangestrèsconstructifsquenous
avonspuavoirtoutaulongdecestroisannées.Mesremerciementsvontensuiteàmesdeuxencadrants
industriels, Clément Chavant et Sylvie Granet, avec qui j’ai eu l’immense plaisir de partager mon
quotidien. Je tiens à leur exprimer toute ma reconnaissance pour tout ce que j’ai pu apprendre à leur
côtédurantcettepériode.Cefutunvéritableplaisirdetravaillerenleurcompagnie.
Je tiens maintenant à adresser mes remerciements aux personnes ayant fait partie de mon jury de
thèse; merci à Guy Bonnet d’avoir présider ce jury, merci à Pascal Omnes et Robert Charlier d’avoir
accepté de rapporter ma thèse. Je les remercie pour leurs remarques et questions très constructives.
MerciégalementàRolandMassond’avoiracceptédefairepartiedemonjurydethèse.
Durant ses trois années, j’ai eu le plaisir de travailler avec de nombreux collègues, thésards et
stagiairesetjelesremerciepourtousleséchangesquenousavonspuavoirainsiquelesbonsmoments
passésàEDFouendehors...
Je remercierai ensuite Renaud (et Boubou,... ) pour m’avoir encouragé, aidé et supporté durant
cetteétape.UngrandmerciàMarion,Fabien,Roro,Alex,Jojo,Sam,Jérôme,BricedeNice,Ons,Schérif,
MedhialiasHugoBossetbiend’autres...avecquij’aipassédetrèsbonsmomentsaucoursdecestrois
années et sans qui cette thèse n’aurait pas été aussi agréable. Merci aussi à Calimero pour son aide
précieuseetsonsoutien.MerciaussiàmoncollèguedebureauAlexaliasCoktailavecquij’aipasséde
trèsbonsmoments.
Pourfinirmerciàmafamille,àmesparents,àmasoeuretàmesgrands parentspourm’avoirépaulé,
encouragéetavoirtoujourscruenmoi.
tel-00587364, version 1 - 20 Apr 2011tel-00587364, version 1 - 20 Apr 2011Tabledesmatières
Introduction 17
Présentationdustockagedesdéchetsradioactifs . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
Difficultésliéesàlamodélisationdiphasiquedustockagededéchetsradioatifs . . . . . . . . . 21
Contexteindustriel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
Plandumanuscrit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
1 Modélisationdesécoulementsenmilieuporeux 25
1.1 GénéralitéssurlesmilieuxporeuxetlanotiondeV.E.R . . . . . . . . . . . . . . . . . . . . 26
1.2 Descriptiondesdifférentesphasesetdeleurspropriétés . . . . . . . . . . . . . . . . . . . 29
1.3 Propriétésetloisdecomportementdesdifférentscomposants . . . . . . . . . . . . . . . . 33
1.4 Conservationdelamasse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
1.5 LoideDarcy Muskat . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
1.6 Diffusionmoléculaire:laloideFick . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
1.7 Nouvelleformulationetchoixdesinconnues . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2 Schémasnumériques 45
2.1 Critèresdesélectiondelaméthodedediscrétisationspatiale . . . . . . . . . . . . . . . . . 46
2.2 Différentsschémasdedicrétisationspatiale . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
2.2.1 SchémasVolumesFinis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
2.2.2ÉlémentsFinisetleursextensions . . . . . . . . . . . . . . . . . . . . . . . 52
2.3 Définitionsetnotations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
2.3.1 Discrétisationentemps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
2.3.2enespace . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
2.3.3 Conventionsd’écriture . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
2.4 MéthodedeVolumesFinisSUSHI . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
2.4.1 Discrétisationdel’équationdediffusionlinéaireàl’aideduschémaSUSHI . . . . 56
2.4.2deséquationsd’écoulementsdiphasiquesenmilieuporeuxàl’aide
duschémaSUSHI . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
2.4.2.1 Choixdesinconnues . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
2.4.2.2 Ecritureduflux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
2.4.2.3 Décentrage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
3 Simulationnumérique:miseenœuvreetapplicationsnumériques 71
3.1 Miseenœuvre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
3.2 Applicationsnumériques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
3.2.1 Rééquilibragedelapressioncapillaired’unbi matériau:castestBO BG . . . . . . 75
3.2.1.1 Présentationduproblèmephysique . . . . . . . . . . . . . . . . . . . . . . 75
3.2.1.2 Discrétisations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
3.2.1.3 RésultatsobtenusavecleschémaVFDA . . . . . . . . . . . . . . . . . . . 76
3.2.1.4 Comparaisondesdifférentesformulations . . . . . . . . . . . . . . . . . 80
3.2.2 Injectiond’hydrogènedansunmilieuporeuxunidimensionnelinitialementsaturé
eneau . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
tel-00587364, version 1 - 20 Apr 20116 Tabledesmatières
3.2.2.1 Présentationduproblèmephysique . . . . . . . . . . . . . . . . . . . . . . 82
3.2.2.2 Discrétisations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
3.2.2.3 Résultats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
3.2.3 Injection d’hydrogène dans un milieu poreux bidimensionnel initialement saturé
eneau . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
3.2.3.1 Présentationduproblèmephysique . . . . . . . . . . . . . . . . . . . . . . 87
3.2.3.2 Discrétisations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
3.2.3.3 Résultats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
3.2.4 Drainaged’unecolonned’eau:castestdeLiakopoulos . . . . . . . . . . . . . . . . 93
3.2.4.1 Présentationduproblèmephysique . . . . . . . . . . . . . . . . . . . . . . 93
3.2.4.2 Discrétisations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
3.2.4.3 Résultats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
3.2.5 Injection d’hydrogène dans un milieu poreux bidimensionnel initialement saturé
eneauenprésenced’unezonedevide . . . . . . . . . . . . . . . . . . . . . . . . . 98
3.2.5.1 Présentationduproblèmephysique . . . . . . . . . . . . . . . . . . . . . . 98
3.2.5.2 Discrétisations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
3.2.5.3 Résultats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
3.2.6 Injectiond’hydrogènedansunmilieuporeuxtridimensionnelinitialementsaturé
eneau . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
3.2.6.1 Présentationduproblèmephysique . . . . . . . . . . . . . . . . . . . . . . 111
3.2.6.2 Discrétisations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
3.2.6.3 Résultats . . . . . . . . . .