The Andra Couplex 1 Test Case: Comparisons Between Finite-Element, Mixed Hybrid Finite Element and Finite Volume Element Discretizations
Issue Title: Simulation of Transport Around a Nuclear Waste Disposal Site: The COUPLEX Test Cases We present here results for the Andra Couplex 1 test case, obtained with the code Cast3m. This code is developped at the CEA (Commissariat a l'energie atomique) and is used mainly to solve problems...
Saved in:
Published in | Computational geosciences Vol. 8; no. 2; pp. 187 - 201 |
---|---|
Main Authors | , , , , |
Format | Journal Article |
Language | English |
Published |
Dordrecht
Springer Nature B.V
2004
Springer Verlag |
Subjects | |
Online Access | Get full text |
Cover
Loading…
Summary: | Issue Title: Simulation of Transport Around a Nuclear Waste Disposal Site: The COUPLEX Test Cases We present here results for the Andra Couplex 1 test case, obtained with the code Cast3m. This code is developped at the CEA (Commissariat a l'energie atomique) and is used mainly to solve problems of solid mechanics, fluid mechanics and heat transfers. Different types of discretization are available, among them finite element, finite volume and mixed hybrid finite element method. Cast3m is also a componant of the platteform Alliances (co-developped by Andra, CEA), which will be used by Andra for the safety calculation of an underground waste disposal in year 2004. We solve the Darcy equation for the water flow and a convection-diffusion transport equation for the Iodine 129 which escapes from a repository cave into the water. The water flow is calculated with a MHFE discretization. It is shown that this method provides sharp results even on relatively coarse grids. The convection-diffusion transport equation is discretized with FE (Finite Element), MHFE (Mixed Hybrid Finite Element) and FV (Finite Volume) methods. In our comparison, we point out the differences of these methods in term of accuracy, respect of the maximum principle and calculations cost. Neither the finite element nor the mixed hybrid finite element approach respects the maximum principle. This results in the presence of negative concentrations near the repository cave, whereas FV calculations respect the monotonicity. We show that mass lumping techniques suppress this problem but with strong restrictions on the grid. FE and MHFE approaches are more accurate than FV for the diffusion equation, but the overall results are equivalent since the advective terms are dominant in the far field and are discretized with centered schemes. We conclude by studying the influence of the grid: a very fine grid near the repository solves almost all the problems of monotonicity, without employing mass lumping techniques. We also observed a very important increase of the accuracy on a structured grid made up of rectangles.[PUBLICATION ABSTRACT] |
---|---|
ISSN: | 1420-0597 1573-1499 |
DOI: | 10.1023/B:COMG.0000035079.68284.49 |