1 / 9

COUPLEX 1 Test Case Nuclear Waste Disposal Far Field Simulation D. VEYRET – M. GEORGIADOU

COUPLEX 1 Test Case Nuclear Waste Disposal Far Field Simulation D. VEYRET – M. GEORGIADOU I.U.S.T.I. – C.N.R.S. UMR 6595 Université de Provence 5 rue Enrico Fermi 13453 Marseille cedex 13 France. Méthode :

kamuzu
Télécharger la présentation

COUPLEX 1 Test Case Nuclear Waste Disposal Far Field Simulation D. VEYRET – M. GEORGIADOU

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. COUPLEX 1 Test Case Nuclear Waste Disposal Far Field Simulation D. VEYRET – M. GEORGIADOU I.U.S.T.I. – C.N.R.S. UMR 6595 Université de Provence 5 rue Enrico Fermi 13453 Marseille cedex 13 France

  2. Méthode : • - Le code utilise la méthode des éléments finis. Il est basé sur un algorithme de type projection pour la résolution des équations de la mécanique des fluides. Le système matriciel qui en résulte est résolu à l’aide de méthodes directes ou bien itératives (type BICGSTAB). • Schéma d’intégration en temps : un schéma de deuxième ordre (de type « trapezoidal rule ») est employé. Il se décompose en deux étapes : • a) b) où • La matrice masse M est transformée de manière à se présenter sous forme diagonale. • - Aucun traitement particulier n’est apporté aux termes convectifs, c’est à dire qu’aucune diffusivité artificielle n’est ajoutée. • - La résolution des systèmes d’équations s’effectue à l’aide d’une méthode de type directe (décomposition LU), préférée à une méthode de type itérative. Elle nécessite plus de ressources mémoire mais délivre un résultat dans un temps CPU inférieur.

  3. Méthode (suite) : • - Pas de temps : le problème a été résolu avec des pas de temps différents suivant la présence ou non du terme source. A titre d’exemple, dans le cas de l’iodine : • un pas de 1 an pour la période allant de 1000 à 3000 ans • un pas de 10 ans pour la période allant de 3000 ans à 10000 ans • un pas de 10000 ans pour la période allant de 10000 ans à 10 000000 ans • Le matériel informatique : les calculs ont été menés sur un ordinateur de type PC (Pentium III - 750 MHz). Le temps de résolution pour un pas de temps est de 12 secondes CPU. • - Le maillage : il est composé exclusivement d’éléments bilinéaires quadrilatéraux à 4 nœuds. Il comprend 48069 nœuds soit 441 nœuds sur X * 109 nœuds sur Z. Le maillage est raffiné près de la zone de stockage où se trouve le terme source et à la frontière entre les différentes couches soit :

  4. Résultats

  5. Résultats

  6. Résultats Vue à l’échelle : Facteur d’échelle sur z : Figure 3 : Champ de Pression

  7. Résultats : Iodine t = 50110 t = 10110 t = 107 t = 106 Figure 4 : Iso-concentrations (facteur d’échelle sur z)

  8. Résultats : Iodine Figure 5 : Flux cumulés

  9. Résultats : Plutonium • Les concentration / flux obtenus sont très inférieurs : • Rw est 106 supérieur => les concentrations sont bien plus faibles • Figure 6 est le seul cas où les concentrations sont supérieures à 10-12 • - Les flux cumulés n’ont pas pu être quantifiés Figure 6 : Iso-concentrations (t = 106)

More Related