new models and small improvements

Universite-Gustave-Eiffel · Nov 24, 2023 · d9d4ed6 · d9d4ed6
1 parent faf0e14
commit d9d4ed6
Show file tree

Hide file tree

Showing 118 changed files with 117,248 additions and 465 deletions.
diff --git a/BilConfig.h b/BilConfig.h
diff --git a/BilExtraLibs.h b/BilExtraLibs.h
diff --git a/OPTIONS b/OPTIONS
@@ -1,4 +1,4 @@
 # This is a comment. Put a # in the first column to comment the line.
 ENABLE_OPENMP   = OFF
 ENABLE_PTHREAD  = OFF
-ENABLE_MPI      = OFF
+ENABLE_MPI      = ON
diff --git a/base/SulfacoESA3d/Ma30 b/base/SulfacoESA3d/Ma30
@@ -0,0 +1,122 @@
+Units
+Length = decimeter
+Time   = second
+Mass   = hectogram
+
+
+
+Geometry
+Dimension = 1  plan
+
+
+
+Mesh
+2
+0. 1
+1.
+1
+1
+
+
+
+Material
+Model = Sulfaco
+porosity = 0.2
+N_CH     = 1.53
+N_CSH    = 1.39
+N_AH3    = 1.e-6
+N_CSH2   = 0
+N_AFm    = 0
+N_AFt    = 0
+N_C3AH6  = 0.114
+R_CSH2   = 1.e-12
+R_AFm    = 1.e-6
+R_C3AH6  = 1.e-6
+Curves_log = Sat  r = Range{r1 = 1.e-8 , r2 = 1.e-3 , n = 1000} S_r = Expressions(1){r0=155.82e-8; m = 0.2516 ; S_r = (1 + (r0/r)**(1/(1-m)))**(-m)}
+A_i = 2.3e-8
+A_p = 1.5e-9
+AlphaCoef = 0
+BetaCoef = 0.87
+K_bulk = 20.e9
+Biot = 1
+Strain0 = 6.e-4
+Strainf = 3.9e-3
+
+
+
+Fields
+1
+Value = 1       Gradient = 0 0 0      Point = 0 0 0
+
+
+
+Initialization
+6
+Region = 1 Unknown = logc_so4    Field = 1 Function = 1
+Region = 1 Unknown = psi         Field = 0 Function = 0
+Region = 1 Unknown = z_ca        Field = 1 Function = 4
+Region = 1 Unknown = z_al        Field = 1 Function = 3
+Region = 1 Unknown = logc_k      Field = 1 Function = 2
+Region = 1 Unknown = logc_oh     Field = 1 Function = 5
+
+
+
+Functions
+5
+N = 2 F(0) = -6        F(86400) = -0.678
+N = 2 F(0) = -5.7      F(86400) = -0.377
+N = 2 F(0) = -7.32     F(86400) = -6.32
+N = 2 F(0) = 1         F(86400) = 1
+N = 2 F(0) = -2        F(86400) = -2    # pH(10) = -4; pH(11) = -3; pH(12) = -2
+#N = 2 F(0) = -8        F(86400) = -0.456
+
+
+
+Boundary Conditions
+4
+Region = 1 Unknown = logc_so4   Field = 1 Function = 1
+Region = 1 Unknown = psi        Field = 0 Function = 0
+#Region = 1 Unknown = z_ca       Field = 1 Function = 4
+#Region = 1 Unknown = z_al       Field = 1 Function = 3
+Region = 1 Unknown = logc_k     Field = 1 Function = 2
+Region = 1 Unknown = logc_oh    Field = 1 Function = 5
+
+
+
+Loads
+0
+
+
+
+Points
+1
+0.5
+
+
+
+Dates
+2
+0.  3.03e6 
+
+
+
+Objective Variations
+logc_so4     = 1.e-3
+z_ca         = 1.e-3
+psi          = 1.e-1
+logc_k       = 1.e-1
+z_al         = 1.e-1
+logc_oh      = 1.e-1
+
+
+
+Iterative Process
+Iter = 20 
+Tol = 1.e-4
+Repetition = 0
+
+
+
+Time Steps
+Dtini = 1
+Dtmax = 1.e4
diff --git a/base/SulfacoESA3d/Ma30.p1 b/base/SulfacoESA3d/Ma30.p1
diff --git a/base/SulfacoESA3d/Ma30.t0 b/base/SulfacoESA3d/Ma30.t0
@@ -0,0 +1,8 @@
+# Version 2.8.2, Fri Dec  4 08:41:50 2020
+# Time = 0.000000e+00
+# Model = Sulfaco
+# Number of views = 59
+# Numbers of components per view = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+# Coordinates(1) ph(4) c_ca(5) c_caoh(6) c_h2sio4(7) c_h3sio4(8) c_h4sio4(9) c_cah2sio4(10) c_cah3sio4(11) c_h2so4(12) c_hso4(13) c_so4(14) c_caso4aq(15) c_cahso4(16) c_al(17) c_alo4h4(18) c_k(19) c_koh(20) C_Ca(21) C_Si(22) C_S(23) C_Al(24) C_K(25) s_ch(26) s_csh2(27) s_ah3(28) s_afm(29) s_aft(30) s_c3ah6(31) n_ch(32) n_csh2(33) n_csh(34) n_ah3(35) n_afm(36) n_aft(37) n_c3ah6(38) n_so4^ads(39) porosite(40) potentiel_electrique(41) charge(42) V_CSH(43) C/S(44) W_Si(45) W_Ca(46) W_S(47) W_Al(48) W_q(49) P_CSH2(50) Damage(51) N_Ca(52) N_Si(53) N_S(54) N_Al(55) N_K(56) N_Cl(57) Saturation degree of crystal(58) Interface equilibrium saturation index of AFt(59) Pore wall equilibrium saturation index of AFt(60) Crystallization pressure(61) Strain(62)
+ 0.000000e+00   0.000000e+00   0.000000e+00   1.259289e+01   1.004081e-02   6.526398e-03   5.530350e-08   4.445893e-07   1.079357e-09   2.210657e-05   7.075020e-08   5.512522e-36   2.158963e-17   1.000000e-06   1.539393e-06   2.766075e-18   1.163105e-33   7.687651e-07   1.995262e-06   1.839898e-08   1.659803e-02   2.267829e-05   2.539393e-06   7.687651e-07   2.013661e-06   1.000000e+00   3.817409e-04   4.786301e-08   1.934294e-03   3.148982e-03   8.536541e-02   1.530000e+00   0.000000e+00   1.390000e+00   0.000000e+00   0.000000e+00   1.255772e-05   1.140000e-01   0.000000e+00   1.999911e-01   0.000000e+00   2.428985e-06   7.992496e-02   1.793498e+00  -0.000000e+00  -0.000000e+00  -0.000000e+00  -0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00   4.368357e+00   1.390005e+00   3.818104e-05   2.280253e-01   4.027323e-07   2.000000e-100  4.460000e-05   1.000583e+00   1.000000e+00   0.000000e+00   0.000000e+00  
+ 1.000000e+00   0.000000e+00   0.000000e+00   1.259289e+01   1.004081e-02   6.526398e-03   5.530350e-08   4.445893e-07   1.079357e-09   2.210657e-05   7.075020e-08   5.512522e-36   2.158963e-17   1.000000e-06   1.539393e-06   2.766075e-18   1.163105e-33   7.687651e-07   1.995262e-06   1.839898e-08   1.659803e-02   2.267829e-05   2.539393e-06   7.687651e-07   2.013661e-06   1.000000e+00   3.817409e-04   4.786301e-08   1.934294e-03   3.148982e-03   8.536541e-02   1.530000e+00   0.000000e+00   1.390000e+00   0.000000e+00   0.000000e+00   1.255772e-05   1.140000e-01   0.000000e+00   1.999911e-01   0.000000e+00   2.428985e-06   7.992496e-02   1.793498e+00  -0.000000e+00  -0.000000e+00  -0.000000e+00  -0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00   4.368357e+00   1.390005e+00   3.818104e-05   2.280253e-01   4.027323e-07   2.000000e-100  4.460000e-05   1.000583e+00   1.000000e+00   0.000000e+00   0.000000e+00  
diff --git a/base/SulfacoESA3d/Ma30.t1 b/base/SulfacoESA3d/Ma30.t1
@@ -0,0 +1,8 @@
+# Version 2.8.2, Fri Dec  4 08:41:53 2020
+# Time = 3.030000e+06
+# Model = Sulfaco
+# Number of views = 59
+# Numbers of components per view = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+# Coordinates(1) ph(4) c_ca(5) c_caoh(6) c_h2sio4(7) c_h3sio4(8) c_h4sio4(9) c_cah2sio4(10) c_cah3sio4(11) c_h2so4(12) c_hso4(13) c_so4(14) c_caso4aq(15) c_cahso4(16) c_al(17) c_alo4h4(18) c_k(19) c_koh(20) C_Ca(21) C_Si(22) C_S(23) C_Al(24) C_K(25) s_ch(26) s_csh2(27) s_ah3(28) s_afm(29) s_aft(30) s_c3ah6(31) n_ch(32) n_csh2(33) n_csh(34) n_ah3(35) n_afm(36) n_aft(37) n_c3ah6(38) n_so4^ads(39) porosite(40) potentiel_electrique(41) charge(42) V_CSH(43) C/S(44) W_Si(45) W_Ca(46) W_S(47) W_Al(48) W_q(49) P_CSH2(50) Damage(51) N_Ca(52) N_Si(53) N_S(54) N_Al(55) N_K(56) N_Cl(57) Saturation degree of crystal(58) Interface equilibrium saturation index of AFt(59) Pore wall equilibrium saturation index of AFt(60) Crystallization pressure(61) Strain(62)
+ 0.000000e+00   0.000000e+00   0.000000e+00   1.216810e+01   7.101598e-02   1.735670e-02   7.819253e-09   1.671727e-07   1.079357e-09   2.210657e-05   1.881574e-07   8.183474e-30   1.205144e-11   2.098940e-01   2.285267e+00   1.092056e-11   4.069414e-34   5.376898e-09   4.197590e-01   1.455459e-03   2.373669e+00   2.247080e-05   2.495161e+00   5.376898e-09   4.212144e-01   1.000000e+00   5.667037e+02   1.656013e-11   9.935129e-01   3.564477e+12   2.953559e-05   9.175283e-01   1.669154e-03   1.390000e+00   0.000000e+00   9.436284e-02   1.965050e-02   0.000000e+00   0.000000e+00   1.941468e-01   0.000000e+00   2.899770e-02   7.992496e-02   1.793498e+00  -0.000000e+00  -0.000000e+00  -0.000000e+00  -0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00   4.368362e+00   1.390004e+00   6.394189e-01   2.280267e-01   8.177879e-02   1.941500e-100  6.973043e-02   1.178536e+00   3.294340e+12   9.885118e+07   3.446468e-04  
+ 1.000000e+00   0.000000e+00   0.000000e+00   1.216810e+01   7.101598e-02   1.735670e-02   7.819253e-09   1.671727e-07   1.079357e-09   2.210657e-05   1.881574e-07   8.183474e-30   1.205144e-11   2.098940e-01   2.285267e+00   1.092056e-11   4.069414e-34   5.376898e-09   4.197590e-01   1.455459e-03   2.373669e+00   2.247080e-05   2.495161e+00   5.376898e-09   4.212144e-01   1.000000e+00   5.667037e+02   1.656013e-11   9.935129e-01   3.564477e+12   2.953559e-05   9.175283e-01   1.669154e-03   1.390000e+00   0.000000e+00   9.436284e-02   1.965050e-02   0.000000e+00   0.000000e+00   1.941468e-01   0.000000e+00   2.899770e-02   7.992496e-02   1.793498e+00  -0.000000e+00  -0.000000e+00  -0.000000e+00  -0.000000e+00   0.000000e+00   0.000000e+00   0.000000e+00   4.368362e+00   1.390004e+00   6.394189e-01   2.280267e-01   8.177879e-02   1.941500e-100  6.973043e-02   1.178536e+00   3.294340e+12   9.885118e+07   3.446468e-04