7. E modeling#

Fluid behavior: THMC = LIQU_SATU

7.1. Characteristics of modeling E#

  • Plane modeling

  • 1 mesh DPQ8 from the D_ PLAN_THM modeling: THM_ DPQ8

_images/100010E20000159A000013A3E04FB40CCC36EFED.svg

7.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 unknowns, fluid pressures 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/point

Order Number

Press

Reference

\((\mathrm{Pa})\) « 

Tolerance \((\text{\%})\)

\(\mathrm{N1}/A\)

1 (t=1 s)

\(\mathrm{PRE1}\)

3.98.10-2

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{N3}/B\)

1 (t=1 s)

\(\mathrm{PRE1}\)

-3,98.10-2

1.0

2 (t=5s)

\(\mathrm{PRE1}\)

-1,99.10-1

1.0

3 (t=10 s)

\(\mathrm{PRE1}\)

-3,98.10-1

2.0

4 (t=50s)

\(\mathrm{PRE1}\)

-1.99

2.0

8 (t=5.103s)

\(\mathrm{PRE1}\)

-1,95.10+2

1.0

16 (t=1010s)

\(\mathrm{PRE1}\)

-5.10+3

1.0