3. Modeling A#

Fluid behavior: THMC = LIQU_SATU

3.1. Characteristics of modeling A#

Plane modeling D_ PLAN_HM

_images/100010E20000159A000013A3E04FB40CCC36EFED.svg

1 mesh DPQ8de modeling D_ PLAN_HM: HM_ DPQ8

3.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/instant \((s)\)

Value

Reference \((\mathrm{Pa})\)

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

\(\mathrm{N1}/A\)

\(1(t=1s)\)

\(\mathrm{PRE1}\)

3,98.10-2

1.0

\(2(t=5s)\)

\(\mathrm{PRE1}\)

1,99.10-1

1.0

\(3(t=10s)\)

\(\mathrm{PRE1}\)

3,98.10-1

1.0

\(4(t=50s)\)

\(\mathrm{PRE1}\)

1.99

1.0

\(8(t={5.10}^{3}s)\)

\(\mathrm{PRE1}\)

1,95.10+2

1.0

\(16(t={10}^{10}s)\)

\(\mathrm{PRE1}\)

5.10+3

1.0

\(\mathrm{N3}/C\)

\(1(t=1s)\)

\(\mathrm{PRE1}\)

-3,98.10-2

1.0

\(2(t=5s)\)

\(\mathrm{PRE1}\)

-1,99.10-1

5.0

\(3(t=10s)\)

\(\mathrm{PRE1}\)

-3,98.10-1

2.0

\(4(t=50s)\)

\(\mathrm{PRE1}\)

-1.99

2.0

\(8(t={5.10}^{3}s)\)

\(\mathrm{PRE1}\)

-1,95.10+2

1.0

\(16(t={10}^{10}s)\)

\(\mathrm{PRE1}\)

-5.10+3

1.0