5. C modeling#
5.1. Characteristics of C modeling#
Same modeling as A but with a vapour density exchange term (via relative humidity RH) in order to test the “ECHANGE_THM_HR” function.
Modeling in 2D plane deformations D_ PLAN_THH2S
Hydraulic law: “LIQU_AD_GAZ_VAPE”
50 \(Q8\) elements refined at the edge of application of the trade term in order to verify the value of the flow at the edge in relation to the analytical solution.
The simulation is carried out over one month. The external humidity is fixed at HR=0.7, the trade term at \({10}^{-3} m/s\).
3 calculations are carried out in succession: the first (U1) with a saturating vapor pressure pvsat=2460 Pa, the second (U2) with the same value but defined as a constant function (validation of AFFE_CHAR_MECA_F) and the third (U3) with a saturating vapor pressure depending on the temperature: \(Pvsat ={10}^{(-2.78+ (T-273.)/(31.559+0.1354*(T-273.)))}\)
5.2. Result of the C modeling#
The capillary pressure values at node 3 (x=0; y=0) are checked as a non-regression value and then it is verified that the water flow calculated by the test case is indeed equal to:
Time |
Calculation |
PRE1 (Pa) |
1 hour |
U1 |
950517 |
1 hour |
U2 |
950517 |
1 hour |
U3 |
1366503 |
1 month |
U1 |
941978 |
1 month |
U2 |
941978 |
It is verified that the flows are indeed equal to the analytical solution.
Points \((x,y)\) |
Time ( \(m\) ) |
“FH11X” + “FH12X” |
Analytical value |
\((\mathrm{0,0};0)\) |
1 |
4.681E-09 |
4.708e-09 |