Final benchmark submitted results

benchmark/spe11c/r3_cp_50m-50m-8mish_convective.txt

@@ -1,6 +1,6 @@
 """Set the full path to the flow executable and flags"""
-mpirun -np 64 flow --relaxed-max-pv-fraction=0 --partition-method=2 --edge-weights-method=2 --imbalance-tol=1.1 --tolerance-mb=1e-7 --enable-tuning=true --enable-opm-rst-file=true --output-extra-convergence-info=steps,iterations --newton-min-iterations=1 --linear-solver=cpr_trueimpes --time-step-control=newtoniterationcount --newton-max-iterations=13 --solver-growth-factor=1.9 --time-step-control-growth-rate=1.5 --enable-drift-compensation=0
- 
+mpirun -np 64 flow --relaxed-max-pv-fraction=0 --partition-method=2 --edge-weights-method=2 --imbalance-tol=1.1 --tolerance-mb=1e-7 --enable-tuning=true --enable-opm-rst-file=true --output-extra-convergence-info=steps,iterations --newton-min-iterations=1 --linear-solver=cpr_trueimpes --time-step-control=newtoniterationcount --newton-max-iterations=13 --solver-growth-factor=1.9 --time-step-control-growth-rate=1.5
+
 """Set the model parameters"""
 spe11c master     #Name of the spe case (spe11a, spe11b, or spe11c) and OPM Flow version (master or release)
 convective gaswater #Name of the co2 model (immiscible, convective [convective requires a Flow version newer than 22-08-2024], or complete) and co2store implementation (gaswater or gasoil [oil properties are set to water internally in OPM flow])
@@ -15,13 +15,13 @@ corner-point      #Type of grid (cartesian, tensor, or corner-point)
 8.5e-1 2500       #Rock specific heat and density (for spe11b/c)
 0 5e4 1           #Added pore volume on top boundary (for spe11a [if 0, free flow bc]), pore volume on lateral boundaries, and width of buffer cell [m] (for spe11b/c)
 150 10            #Elevation of the parabola and back [m] (for spe11c)
- 
+
 """Set the saturation functions"""
 (max(0, (s_w - swi) / (1 - swi))) ** 1.5                                                        #Wetting rel perm saturation function [-]
 (max(0, (1 - s_w - sni) / (1 - sni))) ** 1.5                                                    #Non-wetting rel perm saturation function [-]
 penmax * math.erf(pen * ((s_w-swi) / (1.-swi)) ** (-(1.0 / 1.5)) * math.pi**0.5 / (penmax * 2)) #Capillary pressure saturation function [Pa]
 (np.exp(np.flip(np.linspace(0, 5.0, npoints))) - 1) / (np.exp(5.0) - 1)                         #Points to evaluate the saturation functions (s_w) [-]
- 
+
 """Properties sat functions"""
 """swi [-], sni [-], pen [Pa], penmax [Pa], npoints [-]"""
 SWI1 0.32 SNI1 0.1 PEN1 193531.39 PENMAX1 3e7 NPOINTS1 1000
@@ -31,7 +31,7 @@ SWI4 0.12 SNI4 0.1 PEN4   3870.63 PENMAX4 3e7 NPOINTS4 1000
 SWI5 0.12 SNI5 0.1 PEN5   3060.00 PENMAX5 3e7 NPOINTS5 1000
 SWI6 0.10 SNI6 0.1 PEN6   2560.18 PENMAX6 3e7 NPOINTS6 1000
 SWI7    0 SNI7   0 PEN7         0 PENMAX7 3e7 NPOINTS7    2
- 
+
 """Properties rock"""
 """K [mD], phi [-], disp [m], thconr [W m-1 K-1]"""
 PERM1 0.10132 PORO1 0.10 DISP1 10 THCONR1 1.90
@@ -41,12 +41,12 @@ PERM4 506.625 PORO4 0.20 DISP4 10 THCONR4 1.25
 PERM5 1013.25 PORO5 0.25 DISP5 10 THCONR5 0.92
 PERM6 2026.50 PORO6 0.35 DISP6 10 THCONR6 0.26
 PERM7    1e-5 PORO7 1e-6 DISP7  0 THCONR7 2.00
- 
+
 """Wells radius and position"""
 """radius (0 to use the SOURCE keyword instead of well keywords), x, y, and z position [m] (final positions as well for spe11c)"""
 0 2700. 1000. 300. 2700. 4000. 300. #Well 1
 0 5100. 1000. 700. 5100. 4000. 700. #Well 2
- 
+
 """Define the injection values ([hours] for spe11a; [years] for spe11b/c)"""
 """injection time, time step size to write results, maximum solver time step, injected fluid (0 water, 1 co2) (well1), injection rate [kg/s] (well1), temperature [C] (well1), injected fluid (0 water, 1 co2) (well2), ..."""
 999.9 999.9    10 1  0 10 1  0 10

benchmark/spe11c/r4_cp_8m-8mish-8mish.txt

@@ -1,5 +1,5 @@
 """Set the full path to the flow executable and flags""" 
-mpirun -np 1024 flow --partition-method=2 --edge-weights-method=2 --imbalance-tol=1.1 --tolerance-mb=1e-7 --linear-solver=cprw --enable-tuning=true --enable-opm-rst-file=true --output-extra-convergence-info=steps,iterations --newton-min-iterations=1
+mpirun -np 1024 flow --relaxed-max-pv-fraction=0 --partition-method=2 --edge-weights-method=2 --imbalance-tol=1.1 --tolerance-mb=1e-7 --enable-tuning=true --enable-opm-rst-file=true --output-extra-convergence-info=steps,iterations --newton-min-iterations=1 --linear-solver=cpr_trueimpes --time-step-control=newtoniterationcount --newton-max-iterations=13 --solver-growth-factor=1.9 --time-step-control-growth-rate=1.5
 
 """Set the model parameters"""
 spe11c master     #Name of the spe case (spe11a, spe11b, or spe11c) and OPM Flow version (master or release)
@@ -49,8 +49,10 @@ PERM7    1e-5 PORO7 1e-6 DISP7  0 THCONR7 2.00
 
 """Define the injection values ([hours] for spe11a; [years] for spe11b/c)"""
 """injection time, time step size to write results, maximum solver time step, injected fluid (0 water, 1 co2) (well1), injection rate [kg/s] (well1), temperature [C] (well1), injected fluid (0 water, 1 co2) (well2), ..."""
- 25   5   5 1 50 10 1  0 10
- 25   5   5 1 50 10 1 50 10
- 50  25  25 1  0 10 1  0 10
-400  50  50 1  0 10 1  0 10
-500 100 100 1  0 10 1  0 10
+999.9 999.9    10 1  0 10 1  0 10
+  0.1   0.1   0.1 1  0 10 1  0 10
+   25     5     1 1 50 10 1  0 10
+   25     5   0.3 1 50 10 1 50 10
+   50    25   0.5 1  0 10 1  0 10
+  400    50   0.5 1  0 10 1  0 10
+  500   100     1 1  0 10 1  0 10

docs/_sources/benchmark.rst.txt

@@ -2,10 +2,6 @@
 Benchmark
 *********
 
-.. note::
-    These are preliminary results and will be updated up to the final submission deadline (September 20th 2024, see `here <https://connect.spe.org/discussion/update-on-spe11-submission-and-workshops#bm94ca322f-a2f3-421b-9f90-0191376f7b25>`_).
-    For example, new simulations for the spe11a r4 and spe11c r4 are still running to handle the mass issue and to include the 1000 years of initialization time, respectively. 
-
 All configuration files are located in the `benchmark <https://github.com/OPM/pyopmspe11/blob/main/benchmark>`_ folder.
 
 ======
@@ -137,4 +133,16 @@ Sparse data
 Spatial maps
 ------------
 
-.. figure:: figs/massfractc.png
+.. figure:: figs/massfractc.png
+
+.. note::
+    To show the high resolution results, all spatial maps (spe11a, spe11b, and spe11c) corresponds to the actual simulation grid (not the benchmark reporting grid), and can be generated by typing in the terminal:
+
+    .. code-block:: bash
+
+        plopm -v xco2l -u opm -i 'r1_Cart_1cm/flow/R1_CART_1CM r2_Cart_1cm_capmax2500Pa/flow/R2_CART_1CM_CAPMAX2500PA r3_cp_1cmish_capmax2500Pa/flow/R3_CP_1CMISH_CAPMAX2500PA r4_Cart_1mm_capmax2500Pa/flow/R4_CART_1MM_CAPMAX2500PA' -dpi 2000 -c cet_diverging_protanopic_deuteranopic_bwy_60_95_c32 -cnum 3 -xlnum 8 -clabel 'SPE11A: CO$_2$ mass fraction (liquid phase) after 1 day' -d 16,6.5 -t "r1 Cart 1cm  r2 Cart 1cm capmax 2500 Pa  r3 cp 1cmish capmax 2500 Pa  r4 Cart 1mm capmax 2500 Pa" -yunits cm -xunits cm -yformat .0f -xformat .0f -r 29 -save massfracta -cformat .2e -mask satnum -maskthr 7e-5 -suptitle 0 -subfigs 2,2 -cbsfax 0.35,0.97,0.3,0.02
+        plopm -v xco2l -u opm -i "r1_Cart_10m/flow/R1_CART_10M r2_cp_10mish/flow/R2_CP_10MISH r3_cp_10mish_convective/flow/R3_CP_10MISH_CONVECTIVE r4_Cart_1m/flow/R4_CART_1M" -dpi 2000 -c cet_diverging_protanopic_deuteranopic_bwy_60_95_c32 -cnum 3 -xlnum 8 -clabel 'SPE11B: CO$_2$ mass fraction (liquid phase) after 500 years' -d 16,3 -t "r1 Cart 10m  r2 cp 10mish  r3 cp 10mish convective  r4 Cart 1m" -yunits km -xunits km -yformat .1f -xformat .1f -r 98 -save massfractb -cformat .2e -mask satnum -maskthr 5e-3 -suptitle 0 -subfigs 2,2 -cbsfax 0.35,0.97,0.3,0.02
+        plopm -v xco2l -u opm -i "r1_Cart_50m-50m-10m/flow/R1_CART_50M-50M-10M r2_cp_50m-50m-8mish/flow/R2_CP_50M-50M-8MISH r3_cp_50m-50m-8mish_convective/flow/R3_CP_50M-50M-8MISH_CONVECTIVE r4_cp_8m-8mish-8mish/flow/R4_CP_8M-8MISH-8MISH" -dpi 2000 -c cet_diverging_protanopic_deuteranopic_bwy_60_95_c32 -cnum 3 -xlnum 8 -clabel 'SPE11C: CO$_2$ mass fraction (liquid phase) after 1000 years (y=2.5 [km])' -d 16,3 -t "r1 Cart [50m,50m,10m]  r2 cp [50m,50m,8mish]  r3 cp [50m,50m,8mish] convective  r4 cp [8m,8mish,8mish]" -yunits km -xunits km -yformat .1f -xformat .1f -r 27 -save massfractc -cformat .2e -mask satnum -maskthr 1e-4 -suptitle 0 -subfigs 2,2 -cbsfax 0.30,0.97,0.4,0.02 -s ',51, ,51, ,51, ,304,'
+
+.. tip::
+    You can install plopm by executing in the terminal: pip install git+https://github.com/cssr-tools/plopm.git.

docs/benchmark.html

@@ -104,11 +104,6 @@
 
   <section id="benchmark">
 <h1>Benchmark<a class="headerlink" href="#benchmark" title="Link to this heading"></a></h1>
-<div class="admonition note">
-<p class="admonition-title">Note</p>
-<p>These are preliminary results and will be updated up to the final submission deadline (September 20th 2024, see <a class="reference external" href="https://connect.spe.org/discussion/update-on-spe11-submission-and-workshops#bm94ca322f-a2f3-421b-9f90-0191376f7b25">here</a>).
-For example, new simulations for the spe11a r4 and spe11c r4 are still running to handle the mass issue and to include the 1000 years of initialization time, respectively.</p>
-</div>
 <p>All configuration files are located in the <a class="reference external" href="https://github.com/OPM/pyopmspe11/blob/main/benchmark">benchmark</a> folder.</p>
 <section id="spe11a">
 <h2>SPE11A<a class="headerlink" href="#spe11a" title="Link to this heading"></a></h2>
@@ -266,6 +261,19 @@ <h3>Spatial maps<a class="headerlink" href="#id7" title="Link to this heading">
 <figure class="align-default">
 <img alt="_images/massfractc.png" src="_images/massfractc.png" />
 </figure>
+<div class="admonition note">
+<p class="admonition-title">Note</p>
+<p>To show the high resolution results, all spatial maps (spe11a, spe11b, and spe11c) corresponds to the actual simulation grid (not the benchmark reporting grid), and can be generated by typing in the terminal:</p>
+<div class="highlight-bash notranslate"><div class="highlight"><pre><span></span>plopm<span class="w"> </span>-v<span class="w"> </span>xco2l<span class="w"> </span>-u<span class="w"> </span>opm<span class="w"> </span>-i<span class="w"> </span><span class="s1">&#39;r1_Cart_1cm/flow/R1_CART_1CM r2_Cart_1cm_capmax2500Pa/flow/R2_CART_1CM_CAPMAX2500PA r3_cp_1cmish_capmax2500Pa/flow/R3_CP_1CMISH_CAPMAX2500PA r4_Cart_1mm_capmax2500Pa/flow/R4_CART_1MM_CAPMAX2500PA&#39;</span><span class="w"> </span>-dpi<span class="w"> </span><span class="m">2000</span><span class="w"> </span>-c<span class="w"> </span>cet_diverging_protanopic_deuteranopic_bwy_60_95_c32<span class="w"> </span>-cnum<span class="w"> </span><span class="m">3</span><span class="w"> </span>-xlnum<span class="w"> </span><span class="m">8</span><span class="w"> </span>-clabel<span class="w"> </span><span class="s1">&#39;SPE11A: CO$_2$ mass fraction (liquid phase) after 1 day&#39;</span><span class="w"> </span>-d<span class="w"> </span><span class="m">16</span>,6.5<span class="w"> </span>-t<span class="w"> </span><span class="s2">&quot;r1 Cart 1cm  r2 Cart 1cm capmax 2500 Pa  r3 cp 1cmish capmax 2500 Pa  r4 Cart 1mm capmax 2500 Pa&quot;</span><span class="w"> </span>-yunits<span class="w"> </span>cm<span class="w"> </span>-xunits<span class="w"> </span>cm<span class="w"> </span>-yformat<span class="w"> </span>.0f<span class="w"> </span>-xformat<span class="w"> </span>.0f<span class="w"> </span>-r<span class="w"> </span><span class="m">29</span><span class="w"> </span>-save<span class="w"> </span>massfracta<span class="w"> </span>-cformat<span class="w"> </span>.2e<span class="w"> </span>-mask<span class="w"> </span>satnum<span class="w"> </span>-maskthr<span class="w"> </span>7e-5<span class="w"> </span>-suptitle<span class="w"> </span><span class="m">0</span><span class="w"> </span>-subfigs<span class="w"> </span><span class="m">2</span>,2<span class="w"> </span>-cbsfax<span class="w"> </span><span class="m">0</span>.35,0.97,0.3,0.02
+plopm<span class="w"> </span>-v<span class="w"> </span>xco2l<span class="w"> </span>-u<span class="w"> </span>opm<span class="w"> </span>-i<span class="w"> </span><span class="s2">&quot;r1_Cart_10m/flow/R1_CART_10M r2_cp_10mish/flow/R2_CP_10MISH r3_cp_10mish_convective/flow/R3_CP_10MISH_CONVECTIVE r4_Cart_1m/flow/R4_CART_1M&quot;</span><span class="w"> </span>-dpi<span class="w"> </span><span class="m">2000</span><span class="w"> </span>-c<span class="w"> </span>cet_diverging_protanopic_deuteranopic_bwy_60_95_c32<span class="w"> </span>-cnum<span class="w"> </span><span class="m">3</span><span class="w"> </span>-xlnum<span class="w"> </span><span class="m">8</span><span class="w"> </span>-clabel<span class="w"> </span><span class="s1">&#39;SPE11B: CO$_2$ mass fraction (liquid phase) after 500 years&#39;</span><span class="w"> </span>-d<span class="w"> </span><span class="m">16</span>,3<span class="w"> </span>-t<span class="w"> </span><span class="s2">&quot;r1 Cart 10m  r2 cp 10mish  r3 cp 10mish convective  r4 Cart 1m&quot;</span><span class="w"> </span>-yunits<span class="w"> </span>km<span class="w"> </span>-xunits<span class="w"> </span>km<span class="w"> </span>-yformat<span class="w"> </span>.1f<span class="w"> </span>-xformat<span class="w"> </span>.1f<span class="w"> </span>-r<span class="w"> </span><span class="m">98</span><span class="w"> </span>-save<span class="w"> </span>massfractb<span class="w"> </span>-cformat<span class="w"> </span>.2e<span class="w"> </span>-mask<span class="w"> </span>satnum<span class="w"> </span>-maskthr<span class="w"> </span>5e-3<span class="w"> </span>-suptitle<span class="w"> </span><span class="m">0</span><span class="w"> </span>-subfigs<span class="w"> </span><span class="m">2</span>,2<span class="w"> </span>-cbsfax<span class="w"> </span><span class="m">0</span>.35,0.97,0.3,0.02
+plopm<span class="w"> </span>-v<span class="w"> </span>xco2l<span class="w"> </span>-u<span class="w"> </span>opm<span class="w"> </span>-i<span class="w"> </span><span class="s2">&quot;r1_Cart_50m-50m-10m/flow/R1_CART_50M-50M-10M r2_cp_50m-50m-8mish/flow/R2_CP_50M-50M-8MISH r3_cp_50m-50m-8mish_convective/flow/R3_CP_50M-50M-8MISH_CONVECTIVE r4_cp_8m-8mish-8mish/flow/R4_CP_8M-8MISH-8MISH&quot;</span><span class="w"> </span>-dpi<span class="w"> </span><span class="m">2000</span><span class="w"> </span>-c<span class="w"> </span>cet_diverging_protanopic_deuteranopic_bwy_60_95_c32<span class="w"> </span>-cnum<span class="w"> </span><span class="m">3</span><span class="w"> </span>-xlnum<span class="w"> </span><span class="m">8</span><span class="w"> </span>-clabel<span class="w"> </span><span class="s1">&#39;SPE11C: CO$_2$ mass fraction (liquid phase) after 1000 years (y=2.5 [km])&#39;</span><span class="w"> </span>-d<span class="w"> </span><span class="m">16</span>,3<span class="w"> </span>-t<span class="w"> </span><span class="s2">&quot;r1 Cart [50m,50m,10m]  r2 cp [50m,50m,8mish]  r3 cp [50m,50m,8mish] convective  r4 cp [8m,8mish,8mish]&quot;</span><span class="w"> </span>-yunits<span class="w"> </span>km<span class="w"> </span>-xunits<span class="w"> </span>km<span class="w"> </span>-yformat<span class="w"> </span>.1f<span class="w"> </span>-xformat<span class="w"> </span>.1f<span class="w"> </span>-r<span class="w"> </span><span class="m">27</span><span class="w"> </span>-save<span class="w"> </span>massfractc<span class="w"> </span>-cformat<span class="w"> </span>.2e<span class="w"> </span>-mask<span class="w"> </span>satnum<span class="w"> </span>-maskthr<span class="w"> </span>1e-4<span class="w"> </span>-suptitle<span class="w"> </span><span class="m">0</span><span class="w"> </span>-subfigs<span class="w"> </span><span class="m">2</span>,2<span class="w"> </span>-cbsfax<span class="w"> </span><span class="m">0</span>.30,0.97,0.4,0.02<span class="w"> </span>-s<span class="w"> </span><span class="s1">&#39;,51, ,51, ,51, ,304,&#39;</span>
+</pre></div>
+</div>
+</div>
+<div class="admonition tip">
+<p class="admonition-title">Tip</p>
+<p>You can install plopm by executing in the terminal: pip install git+https://github.com/cssr-tools/plopm.git.</p>
+</div>
 </section>
 </section>
 </section>

docs/text/benchmark.rst

@@ -2,10 +2,6 @@
 Benchmark
 *********
 
-.. note::
-    These are preliminary results and will be updated up to the final submission deadline (September 20th 2024, see `here <https://connect.spe.org/discussion/update-on-spe11-submission-and-workshops#bm94ca322f-a2f3-421b-9f90-0191376f7b25>`_).
-    For example, new simulations for the spe11a r4 and spe11c r4 are still running to handle the mass issue and to include the 1000 years of initialization time, respectively. 
-
 All configuration files are located in the `benchmark <https://github.com/OPM/pyopmspe11/blob/main/benchmark>`_ folder.
 
 ======
@@ -137,4 +133,16 @@ Sparse data
 Spatial maps
 ------------
 
-.. figure:: figs/massfractc.png
+.. figure:: figs/massfractc.png
+
+.. note::
+    To show the high resolution results, all spatial maps (spe11a, spe11b, and spe11c) corresponds to the actual simulation grid (not the benchmark reporting grid), and can be generated by typing in the terminal:
+
+    .. code-block:: bash
+
+        plopm -v xco2l -u opm -i 'r1_Cart_1cm/flow/R1_CART_1CM r2_Cart_1cm_capmax2500Pa/flow/R2_CART_1CM_CAPMAX2500PA r3_cp_1cmish_capmax2500Pa/flow/R3_CP_1CMISH_CAPMAX2500PA r4_Cart_1mm_capmax2500Pa/flow/R4_CART_1MM_CAPMAX2500PA' -dpi 2000 -c cet_diverging_protanopic_deuteranopic_bwy_60_95_c32 -cnum 3 -xlnum 8 -clabel 'SPE11A: CO$_2$ mass fraction (liquid phase) after 1 day' -d 16,6.5 -t "r1 Cart 1cm  r2 Cart 1cm capmax 2500 Pa  r3 cp 1cmish capmax 2500 Pa  r4 Cart 1mm capmax 2500 Pa" -yunits cm -xunits cm -yformat .0f -xformat .0f -r 29 -save massfracta -cformat .2e -mask satnum -maskthr 7e-5 -suptitle 0 -subfigs 2,2 -cbsfax 0.35,0.97,0.3,0.02
+        plopm -v xco2l -u opm -i "r1_Cart_10m/flow/R1_CART_10M r2_cp_10mish/flow/R2_CP_10MISH r3_cp_10mish_convective/flow/R3_CP_10MISH_CONVECTIVE r4_Cart_1m/flow/R4_CART_1M" -dpi 2000 -c cet_diverging_protanopic_deuteranopic_bwy_60_95_c32 -cnum 3 -xlnum 8 -clabel 'SPE11B: CO$_2$ mass fraction (liquid phase) after 500 years' -d 16,3 -t "r1 Cart 10m  r2 cp 10mish  r3 cp 10mish convective  r4 Cart 1m" -yunits km -xunits km -yformat .1f -xformat .1f -r 98 -save massfractb -cformat .2e -mask satnum -maskthr 5e-3 -suptitle 0 -subfigs 2,2 -cbsfax 0.35,0.97,0.3,0.02
+        plopm -v xco2l -u opm -i "r1_Cart_50m-50m-10m/flow/R1_CART_50M-50M-10M r2_cp_50m-50m-8mish/flow/R2_CP_50M-50M-8MISH r3_cp_50m-50m-8mish_convective/flow/R3_CP_50M-50M-8MISH_CONVECTIVE r4_cp_8m-8mish-8mish/flow/R4_CP_8M-8MISH-8MISH" -dpi 2000 -c cet_diverging_protanopic_deuteranopic_bwy_60_95_c32 -cnum 3 -xlnum 8 -clabel 'SPE11C: CO$_2$ mass fraction (liquid phase) after 1000 years (y=2.5 [km])' -d 16,3 -t "r1 Cart [50m,50m,10m]  r2 cp [50m,50m,8mish]  r3 cp [50m,50m,8mish] convective  r4 cp [8m,8mish,8mish]" -yunits km -xunits km -yformat .1f -xformat .1f -r 27 -save massfractc -cformat .2e -mask satnum -maskthr 1e-4 -suptitle 0 -subfigs 2,2 -cbsfax 0.30,0.97,0.4,0.02 -s ',51, ,51, ,51, ,304,'
+
+.. tip::
+    You can install plopm by executing in the terminal: pip install git+https://github.com/cssr-tools/plopm.git.