Updating docs (spe11b using sources and dispersion)

build_opm-flow_macOS.bash

@@ -26,14 +26,14 @@ for repo in common grid models
 do
     mkdir build/opm-$repo
     cd build/opm-$repo
-    cmake -DPYTHON_EXECUTABLE=$(which python) -DUSE_MPI=0 -DOPM_ENABLE_PYTHON=ON -DCMAKE_BUILD_TYPE=Release -DCMAKE_PREFIX_PATH="$CURRENT_DIRECTORY/dune-common/build-cmake;$CURRENT_DIRECTORY/dune-grid/build-cmake;$CURRENT_DIRECTORY/dune-geometry/build-cmake;$CURRENT_DIRECTORY/dune-istl/build-cmake;$CURRENT_DIRECTORY/build/opm-common;$CURRENT_DIRECTORY/build/opm-grid" $CURRENT_DIRECTORY/opm-$repo
+    cmake -DPYTHON_EXECUTABLE=$(which python) -DUSE_MPI=0 -DNDEBUG=1 -DOPM_ENABLE_PYTHON=ON -DCMAKE_BUILD_TYPE=Release -DCMAKE_PREFIX_PATH="$CURRENT_DIRECTORY/dune-common/build-cmake;$CURRENT_DIRECTORY/dune-grid/build-cmake;$CURRENT_DIRECTORY/dune-geometry/build-cmake;$CURRENT_DIRECTORY/dune-istl/build-cmake;$CURRENT_DIRECTORY/build/opm-common;$CURRENT_DIRECTORY/build/opm-grid" $CURRENT_DIRECTORY/opm-$repo
     make -j5
     cd ../..
 done    
 
 mkdir build/opm-simulators
 cd build/opm-simulators
-cmake -DUSE_MPI=0 -DCMAKE_BUILD_TYPE=Release -DCMAKE_PREFIX_PATH="$CURRENT_DIRECTORY/dune-common/build-cmake;$CURRENT_DIRECTORY/dune-grid/build-cmake;$CURRENT_DIRECTORY/dune-geometry/build-cmake;$CURRENT_DIRECTORY/dune-istl/build-cmake;$CURRENT_DIRECTORY/build/opm-common;$CURRENT_DIRECTORY/build/opm-grid;$CURRENT_DIRECTORY/build/opm-models" $CURRENT_DIRECTORY/opm-simulators
+cmake -DUSE_MPI=0 -DNDEBUG=1 -DCMAKE_BUILD_TYPE=Release -DCMAKE_PREFIX_PATH="$CURRENT_DIRECTORY/dune-common/build-cmake;$CURRENT_DIRECTORY/dune-grid/build-cmake;$CURRENT_DIRECTORY/dune-geometry/build-cmake;$CURRENT_DIRECTORY/dune-istl/build-cmake;$CURRENT_DIRECTORY/build/opm-common;$CURRENT_DIRECTORY/build/opm-grid;$CURRENT_DIRECTORY/build/opm-models" $CURRENT_DIRECTORY/opm-simulators
 make -j5 flow
 cd ../..
 

build_opm-flow_mpi.bash

@@ -12,13 +12,13 @@ for repo in common grid models
 do
     mkdir build/opm-$repo
     cd build/opm-$repo
-    cmake -DUSE_MPI=1 -DOPM_ENABLE_PYTHON=ON -DCMAKE_BUILD_TYPE=Release -DCMAKE_PREFIX_PATH="$CURRENT_DIRECTORY/build/opm-common;$CURRENT_DIRECTORY/build/opm-grid" $CURRENT_DIRECTORY/opm-$repo
+    cmake -DUSE_MPI=1 -DNDEBUG=1 -DOPM_ENABLE_PYTHON=ON -DCMAKE_BUILD_TYPE=Release -DCMAKE_PREFIX_PATH="$CURRENT_DIRECTORY/build/opm-common;$CURRENT_DIRECTORY/build/opm-grid" $CURRENT_DIRECTORY/opm-$repo
     make -j5
     cd ../..
 done    
 
 mkdir build/opm-simulators
 cd build/opm-simulators
-cmake -DUSE_MPI=1 -DCMAKE_BUILD_TYPE=Release -DCMAKE_PREFIX_PATH="$CURRENT_DIRECTORY/build/opm-common;$CURRENT_DIRECTORY/build/opm-grid;$CURRENT_DIRECTORY/build/opm-models" $CURRENT_DIRECTORY/opm-simulators
+cmake -DUSE_MPI=1 -DNDEBUG=1 -DCMAKE_BUILD_TYPE=Release -DCMAKE_PREFIX_PATH="$CURRENT_DIRECTORY/build/opm-common;$CURRENT_DIRECTORY/build/opm-grid;$CURRENT_DIRECTORY/build/opm-models" $CURRENT_DIRECTORY/opm-simulators
 make -j5 flow
 cd ../..

docs/_sources/configuration_file.rst.txt

@@ -8,7 +8,7 @@ The first input parameter in the configuration file is:
     :linenos:
 
     """Set the full path to the flow executable and flags"""
-    flow --linear-solver=cprw --newton-min-iterations=1 --enable-tuning=true --enable-opm-rst-file=true --output-extra-convergence-info=steps,iterations  --enable-well-operability-check=false --min-time-step-before-shutting-problematic-wells-in-days=1e-99 
+    flow --linear-solver=cprw --newton-min-iterations=1 --enable-tuning=true --enable-opm-rst-file=true --output-extra-convergence-info=steps,iterations 
 
 If **flow** is not in your path, then write the full path to the executable
 (e.g., /Users/dmar/opm/build/opm-simulators/bin/flow). We also add in the same 
@@ -38,7 +38,7 @@ The following input lines in the configuration file are:
     12                #If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (18 entries)
     70 36.25          #Temperature bottom and top rig [C]            
     19620000 0.1      #Pressure on the top [Pa] and multiplier for the permeability in the z direction [-]
-    1e-9 2e-8 0       #Diffusion (in liquid and gas) [m^2/s] and dispersion [m] (dispersion only available in Flow master)
+    1e-9 2e-8 10      #Diffusion (in liquid and gas) [m^2/s] and dispersion [m] (dispersion requires a Flow version newer than 17-11-2023)
     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) 
@@ -61,11 +61,7 @@ resulting in 18 levels. Then, the z-refinment in each of these levels is set. Se
 
     Examples of a cartesian (top), tensor (middle), and corner-point (bottom) grids. The cartesian is generated using 170, 100, and 120
     elements, the tensor grid by setting the x-array to 50,100,20; the y-array to 5,10,70,10,5; and the z-array to 5,10,20,20,20,20,20,5; and 
-    the corner-point grid using the same xy-arrays as in the tensor grid and for the z-refinment 4,4,3,3,5,3,7,5,13,5,4,3,5,3,9,21,21,2. 
-
-.. warning::
-    Dispersion seems to work fine for the spe11a case where there are no thermal effects nor water evaporation. However, for spe11b/c there are
-    issues with the current implementation (we are working in this, then for now disperison should be set to 0 for spe11b/c).  
+    the corner-point grid using the same xy-arrays as in the tensor grid and for the z-refinment 4,4,3,3,5,3,7,5,13,5,4,3,5,3,9,21,21,2.  
 
 ***********************
 Soil-related parameters
@@ -117,16 +113,16 @@ The last part of the configuration file sets the wells radius, location, and the
     :lineno-start: 45
 
     """Wells radius and position"""
-    """radius, x, y, and z position [m] (final positions as well for spe11c)"""
+    """radius (0 to use the SOURCE keyword instead of well keywords, this requires a Flow version newer than 23-01-2024), x, y, and z position [m] (final positions as well for spe11c)"""
     0.15 2700. 1000. 300. 2700. 4000. 300. #Well 1 
     0.15 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), ..."""
-    1000 1000   1 0  0 10 0  0 10
-      25    5 0.1 1 50 10 0  0 10
+    1000 1000 0.1 1  0 10 1  0 10
+      25    5 0.1 1 50 10 1  0 10
       25    5 0.1 1 50 10 1 50 10
-     950    5 0.1 0  0 10 0  0 10
+     950    5 0.1 1  0 10 1  0 10
 
 .. warning::
     Keep the linebreak between the sections (in the current implementation this is used for the reading of the parameters).

docs/configuration_file.html

@@ -89,7 +89,7 @@
 <h1>Configuration file<a class="headerlink" href="#configuration-file" title="Link to this heading"></a></h1>
 <p>The first input parameter in the configuration file is:</p>
 <div class="highlight-python notranslate"><div class="highlight"><pre><span></span><span class="linenos">1</span><span class="sd">&quot;&quot;&quot;Set the full path to the flow executable and flags&quot;&quot;&quot;</span>
-<span class="linenos">2</span><span class="n">flow</span> <span class="o">--</span><span class="n">linear</span><span class="o">-</span><span class="n">solver</span><span class="o">=</span><span class="n">cprw</span> <span class="o">--</span><span class="n">newton</span><span class="o">-</span><span class="nb">min</span><span class="o">-</span><span class="n">iterations</span><span class="o">=</span><span class="mi">1</span> <span class="o">--</span><span class="n">enable</span><span class="o">-</span><span class="n">tuning</span><span class="o">=</span><span class="n">true</span> <span class="o">--</span><span class="n">enable</span><span class="o">-</span><span class="n">opm</span><span class="o">-</span><span class="n">rst</span><span class="o">-</span><span class="n">file</span><span class="o">=</span><span class="n">true</span> <span class="o">--</span><span class="n">output</span><span class="o">-</span><span class="n">extra</span><span class="o">-</span><span class="n">convergence</span><span class="o">-</span><span class="n">info</span><span class="o">=</span><span class="n">steps</span><span class="p">,</span><span class="n">iterations</span>  <span class="o">--</span><span class="n">enable</span><span class="o">-</span><span class="n">well</span><span class="o">-</span><span class="n">operability</span><span class="o">-</span><span class="n">check</span><span class="o">=</span><span class="n">false</span> <span class="o">--</span><span class="nb">min</span><span class="o">-</span><span class="n">time</span><span class="o">-</span><span class="n">step</span><span class="o">-</span><span class="n">before</span><span class="o">-</span><span class="n">shutting</span><span class="o">-</span><span class="n">problematic</span><span class="o">-</span><span class="n">wells</span><span class="o">-</span><span class="ow">in</span><span class="o">-</span><span class="n">days</span><span class="o">=</span><span class="mf">1e-99</span>
+<span class="linenos">2</span><span class="n">flow</span> <span class="o">--</span><span class="n">linear</span><span class="o">-</span><span class="n">solver</span><span class="o">=</span><span class="n">cprw</span> <span class="o">--</span><span class="n">newton</span><span class="o">-</span><span class="nb">min</span><span class="o">-</span><span class="n">iterations</span><span class="o">=</span><span class="mi">1</span> <span class="o">--</span><span class="n">enable</span><span class="o">-</span><span class="n">tuning</span><span class="o">=</span><span class="n">true</span> <span class="o">--</span><span class="n">enable</span><span class="o">-</span><span class="n">opm</span><span class="o">-</span><span class="n">rst</span><span class="o">-</span><span class="n">file</span><span class="o">=</span><span class="n">true</span> <span class="o">--</span><span class="n">output</span><span class="o">-</span><span class="n">extra</span><span class="o">-</span><span class="n">convergence</span><span class="o">-</span><span class="n">info</span><span class="o">=</span><span class="n">steps</span><span class="p">,</span><span class="n">iterations</span>
 </pre></div>
 </div>
 <p>If <strong>flow</strong> is not in your path, then write the full path to the executable
@@ -113,7 +113,7 @@ <h2>Reservoir-related parameters<a class="headerlink" href="#reservoir-related-p
 <span class="linenos">11</span><span class="mi">12</span>                <span class="c1">#If cartesian, number of z cells [-]; if tensor, variable array of z-refinment; if corner-point, fix array of z-refinment (18 entries)</span>
 <span class="linenos">12</span><span class="mi">70</span> <span class="mf">36.25</span>          <span class="c1">#Temperature bottom and top rig [C]</span>
 <span class="linenos">13</span><span class="mi">19620000</span> <span class="mf">0.1</span>      <span class="c1">#Pressure on the top [Pa] and multiplier for the permeability in the z direction [-]</span>
-<span class="linenos">14</span><span class="mf">1e-9</span> <span class="mf">2e-8</span> <span class="mi">0</span>       <span class="c1">#Diffusion (in liquid and gas) [m^2/s] and dispersion [m] (dispersion only available in Flow master)</span>
+<span class="linenos">14</span><span class="mf">1e-9</span> <span class="mf">2e-8</span> <span class="mi">10</span>      <span class="c1">#Diffusion (in liquid and gas) [m^2/s] and dispersion [m] (dispersion requires a Flow version newer than 17-11-2023)</span>
 <span class="linenos">15</span><span class="mf">8.5e-1</span> <span class="mi">2500</span>       <span class="c1">#Rock specific heat and density (for spe11b/c)</span>
 <span class="linenos">16</span><span class="mi">0</span> <span class="mf">5e4</span> <span class="mi">1</span>           <span class="c1">#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)</span>
 <span class="linenos">17</span><span class="mi">150</span> <span class="mi">10</span>            <span class="c1">#Elevation of the parabola and back [m] (for spe11c)</span>
@@ -144,11 +144,6 @@ <h2>Reservoir-related parameters<a class="headerlink" href="#reservoir-related-p
 the corner-point grid using the same xy-arrays as in the tensor grid and for the z-refinment 4,4,3,3,5,3,7,5,13,5,4,3,5,3,9,21,21,2.</span><a class="headerlink" href="#id1" title="Link to this image"></a></p>
 </figcaption>
 </figure>
-<div class="admonition warning">
-<p class="admonition-title">Warning</p>
-<p>Dispersion seems to work fine for the spe11a case where there are no thermal effects nor water evaporation. However, for spe11b/c there are
-issues with the current implementation (we are working in this, then for now disperison should be set to 0 for spe11b/c).</p>
-</div>
 </section>
 <section id="soil-related-parameters">
 <h2>Soil-related parameters<a class="headerlink" href="#soil-related-parameters" title="Link to this heading"></a></h2>
@@ -194,16 +189,16 @@ <h2>Soil-related parameters<a class="headerlink" href="#soil-related-parameters"
 <h2>Well-related parameters<a class="headerlink" href="#well-related-parameters" title="Link to this heading"></a></h2>
 <p>The last part of the configuration file sets the wells radius, location, and the injection schedule.</p>
 <div class="highlight-python notranslate"><div class="highlight"><pre><span></span><span class="linenos">45</span><span class="sd">&quot;&quot;&quot;Wells radius and position&quot;&quot;&quot;</span>
-<span class="linenos">46</span><span class="sd">&quot;&quot;&quot;radius, x, y, and z position [m] (final positions as well for spe11c)&quot;&quot;&quot;</span>
+<span class="linenos">46</span><span class="sd">&quot;&quot;&quot;radius (0 to use the SOURCE keyword instead of well keywords, this requires a Flow version newer than 23-01-2024), x, y, and z position [m] (final positions as well for spe11c)&quot;&quot;&quot;</span>
 <span class="linenos">47</span><span class="mf">0.15</span> <span class="mf">2700.</span> <span class="mf">1000.</span> <span class="mf">300.</span> <span class="mf">2700.</span> <span class="mf">4000.</span> <span class="mf">300.</span> <span class="c1">#Well 1</span>
 <span class="linenos">48</span><span class="mf">0.15</span> <span class="mf">5100.</span> <span class="mf">1000.</span> <span class="mf">700.</span> <span class="mf">5100.</span> <span class="mf">4000.</span> <span class="mf">700.</span> <span class="c1">#Well 2</span>
 <span class="linenos">49</span>
 <span class="linenos">50</span><span class="sd">&quot;&quot;&quot;Define the injection values ([hours] for spe11a; [years] for spe11b/c)&quot;&quot;&quot;</span>
 <span class="linenos">51</span><span class="sd">&quot;&quot;&quot;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), ...&quot;&quot;&quot;</span>
-<span class="linenos">52</span><span class="mi">1000</span> <span class="mi">1000</span>   <span class="mi">1</span> <span class="mi">0</span>  <span class="mi">0</span> <span class="mi">10</span> <span class="mi">0</span>  <span class="mi">0</span> <span class="mi">10</span>
-<span class="linenos">53</span>  <span class="mi">25</span>    <span class="mi">5</span> <span class="mf">0.1</span> <span class="mi">1</span> <span class="mi">50</span> <span class="mi">10</span> <span class="mi">0</span>  <span class="mi">0</span> <span class="mi">10</span>
+<span class="linenos">52</span><span class="mi">1000</span> <span class="mi">1000</span> <span class="mf">0.1</span> <span class="mi">1</span>  <span class="mi">0</span> <span class="mi">10</span> <span class="mi">1</span>  <span class="mi">0</span> <span class="mi">10</span>
+<span class="linenos">53</span>  <span class="mi">25</span>    <span class="mi">5</span> <span class="mf">0.1</span> <span class="mi">1</span> <span class="mi">50</span> <span class="mi">10</span> <span class="mi">1</span>  <span class="mi">0</span> <span class="mi">10</span>
 <span class="linenos">54</span>  <span class="mi">25</span>    <span class="mi">5</span> <span class="mf">0.1</span> <span class="mi">1</span> <span class="mi">50</span> <span class="mi">10</span> <span class="mi">1</span> <span class="mi">50</span> <span class="mi">10</span>
-<span class="linenos">55</span> <span class="mi">950</span>    <span class="mi">5</span> <span class="mf">0.1</span> <span class="mi">0</span>  <span class="mi">0</span> <span class="mi">10</span> <span class="mi">0</span>  <span class="mi">0</span> <span class="mi">10</span>
+<span class="linenos">55</span> <span class="mi">950</span>    <span class="mi">5</span> <span class="mf">0.1</span> <span class="mi">1</span>  <span class="mi">0</span> <span class="mi">10</span> <span class="mi">1</span>  <span class="mi">0</span> <span class="mi">10</span>
 </pre></div>
 </div>
 <div class="admonition warning">