This example is based on the The TASK D_THM1 of the DECOVALEXTHMC project. The TASK D_THM1 of the DECOVALEXTHMC project studies the coupled thermal hydraulic and mechanical (THM) processes in the FEBEX type nuclear waster repository. In the FEBEX type repository, nuclear waster canisters are stored in the drifts excavated in the deep rock mass of granite, which is fully water saturated, and they are sealed with bentonite. Initially, the bentonite, the sealing material, is partially saturated. TASK D_THM1 defines a 2D model about the FEBEX type repository \cite BirEtAl:2008. In this example, TASK D_THM1 is simplified in order to test the staggered scheme for TH process in OGS. The simplifications are
resizing the domain to an area that can represent the near field of an installed nuclear water canister,
assuming the bentonite is fully saturated from the beginning,
ignoring the mechanical process.
With such simplifications, the geometry of the present example is illustrated in the following figure:
In the above figure, the domain in the annulus sector represents the sealing material, bentonite. A heat power, which is generated by the nuclear waste with one million year variation, is applied onto the inner arc of the annulus sector. On the top boundary, the boundary conditions are p = 4.3 ⋅ 10^{6} Pa, T = 294 K. While on the bottom boundary, the boundary conditions are set as p = 4.7 ⋅ 10^{6} Pa, T = 310 K. The initial conditions are given as p = 4.7 ⋅ 10^{6} Pa, T = 298 K.
The material properties are shown in the following table:
Property  Value  Unit 

Bentonite  
Density  kg/m^{3}  1600 
Porosity    0.01 
Thermal conductivity  W/(mK)  3 
Specific heat capacity  J/(kgK)  3 
Saturated permeability  m^{2}  2.0 ⋅ 10^{−21} 
Granite  
Density  kg/m^{3}  2700 
Porosity    0.41 
Thermal conductivity  W/(mK)  3 
Specific heat capacity  J/(kgK)  900 
Saturated permeability  m^{2}  10^{−17} 
As the reference results, the temperature and pressure distributions in the domain at the time of 18 years are shown in the following figure, in which the thermal convection effective can be seen clearly.
This article was written by Wenqing Wang. If you are missing something or you find an error please let us know.
Generated with Hugo 0.101.0
in CI job 282302

Last revision: November 22, 2022
Commit: [PL/THM] Implement freezing for temperature eq. 68ebbec
 Edit this page on