5. C modeling#
Fluid behavior: THMC = LIQU_SATU
5.1. Characteristics of C modeling#
Volume modeling: 3D_HM
1 mesh HEXA20 from the 3D_HM modeling: HM_HEXA20
5.2. Tested sizes and results#
Discretization in time: Several time steps (16) to study the evolution of pressure during the transition phase until it stabilizes. The time pattern is implicit \((\vartheta =1)\).
List of calculation times in seconds:
\(\mathrm{1,}\mathrm{5,}\mathrm{10,}\mathrm{50,}\mathrm{100,}\mathrm{500,}{10}^{3},{5.10}^{3},{10}^{4},5.{10}^{4},{10}^{5},5.{10}^{5},{10}^{6},5.{10}^{6},{10}^{7},{10}^{10}\)
The nodal fluid pressure unknowns evaluated in Code_Aster are variations from the initial reference pressures defined under the keyword THM_INIT, which is why this table shows pressure variations in our comparison between the Code_Aster calculation and the reference solution.
Node |
Order Number |
Press |
Reference |
||
\((\mathrm{Pa})\) « |
Tolerance \((\text{\%})\) |
||||
\(\mathrm{NO20}\) |
1 (t=1 s) |
\(\mathrm{PRE1}\) |
|
1.0 |
|
2 (t=5s) |
\(\mathrm{PRE1}\) |
1,99.10-1 |
1.0 |
||
3 (t=10 s) |
\(\mathrm{PRE1}\) |
3,98.10-1 |
1.0 |
||
4 (t=50 s) |
\(\mathrm{PRE1}\) |
1.99 |
1.0 |
||
8 (t=5.103s) |
\(\mathrm{PRE1}\) |
1,95.10+2 |
1.0 |
||
16 (t=1010s) |
\(\mathrm{PRE1}\) |
5.10+3 |
1.0 |
||
\(\mathrm{NO1}\) |
1 (t=1 s) |
\(\mathrm{PRE1}\) |
|
1.0 |
|
2 (t=5s) |
\(\mathrm{PRE1}\) |
|
1.0 |
||
3 (t=10 s) |
\(\mathrm{PRE1}\) |
|
2.0 |
||
4 (t=50s) |
\(\mathrm{PRE1}\) |
|
2.0 |
||
8 (t=5.103s) |
\(\mathrm{PRE1}\) |
|
1.0 |
||
16 (t=1010s) |
\(\mathrm{PRE1}\) |
|
1.0 |