6. D modeling

6.1. Characteristics of modeling D

The results presented come from 3D modeling_ THH2MS. This is the same modeling as modeling A but with a non-linear exchange term in HR and not in pressure.

To complete, the calculations are also carried out in 3D_ HH2MS, in 3D_ THH2S and in 3D_ HH2S.

Since movements and temperatures are blocked, the results should be the same in all cases. The aim here is to verify that this trade term has been taken into account.

The mesh is composed of 640 HEXA20 elements.

6.2. Tested values and results

We check that the water flow (“FH11X” + “FH12X” recovered by SIEF_NOEU) obtained on the taxation edge corresponds to its analytical expression, namely \(q_{w}^{ext}{.}n= \alpha.\rho_{vp}.\left( exp\left[ \frac{Pc.M_{H2O}}{\rho_{l}.R.T}\right] -{HR}^{ext} \right)\). A tolerance of 10% is allowed due to the size of the mesh and the error generated by the interpolation at the node.

Points \((x,y)\)	Time ( \(s\) )	“FH11X” + “FH12X”	Analytical value
\((\mathrm{0,2};0)\)	1800	7.3754E-09	8.09e-09

For reasons of saving time, models other than 3D_ THH2MS will only be tested for 60 s. It is verified that the pressures at 60 s are all the same regardless of the models.

Points \((x,y)\)	Time ( \(s\) )	PRE1
\((\mathrm{0,2};0)\)	1800	10017