Skip to content
Snippets Groups Projects
Select Git revision
  • feature/install-runtest
  • releases/3.10 protected
  • fix/test-donea-box-stabilized
  • fix/free-flow-failed-tests-debug
  • master default protected
  • doc/update-changelog
  • test-floating-point-compiler-flags
  • fix/test-reduce-timestep-variation
  • doc/newton-documentation
  • feature/hybrid-pq1-new-problem
  • feature/hybrid-pq1
  • feature/cvfe-localvariables
  • feature/cvfe-localdofs
  • cleanup/py-compatibility-code
  • releases/3.9 protected
  • fix/intersections-3d3d
  • feature/constant-extrusion
  • fix/checkpoints_failed_for_small_epsilon
  • feature/air-henry-implementation
  • feature/test-el1p-mandel
  • 3.10.0-rc1
  • 3.9.0
  • 3.9.0-rc1
  • 3.8.0
  • 3.8.0-rc1
  • 3.7.0
  • 3.7.0-rc1
  • 3.6.0
  • 3.6.0-rc1
  • 3.5.0
  • 3.5.0-rc1
  • 3.4.0
  • 3.4.0-rc1
  • 3.3.0
  • 3.3.0-rc1
  • 3.2.0
  • 3.2.0-rc2
  • 3.2.0-rc1
  • 3.1.0
  • 3.1.0-rc1
40 results
Blame Permalink
  • Martin Beck's avatar
    91d66f10
    [effectivepermeabilityrutqvisttsang] Outsourced permeability calculation · 91d66f10
    Martin Beck authored 
    Reviewed by heck

Instead of calculating it in the local residual, the effective permeability
is now calculated within a new file permeabilityrutqvisttsang.hh in the
folder dumux/material/fluidmatrixinteractions.
This is just one possible effective permeability model. Several others can be
added and set via the property EffectivePermeabilityModel.

Up to now, the permeabilities were averaged twice, first by calculating
the mean of the mean of the intrinsic permeabilities in the flux variables
and later by averaging the effective permeabilities in the local residual.
This is redundant. Now, only the mean of the effective permeability is
calculated in the local residual.
This change required to change the reference files for the sequential and
the parallel test.
    91d66f10
    History
    [effectivepermeabilityrutqvisttsang] Outsourced permeability calculation
    Martin Beck authored 
    Reviewed by heck

Instead of calculating it in the local residual, the effective permeability
is now calculated within a new file permeabilityrutqvisttsang.hh in the
folder dumux/material/fluidmatrixinteractions.
This is just one possible effective permeability model. Several others can be
added and set via the property EffectivePermeabilityModel.

Up to now, the permeabilities were averaged twice, first by calculating
the mean of the mean of the intrinsic permeabilities in the flux variables
and later by averaging the effective permeabilities in the local residual.
This is redundant. Now, only the mean of the effective permeability is
calculated in the local residual.
This change required to change the reference files for the sequential and
the parallel test.
el2p-parallel-reference.vtu 15.44 KiB 
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
  <UnstructuredGrid>
    <Piece NumberOfCells="16" NumberOfPoints="50">
      <PointData Scalars="T" Vectors="u">
        <DataArray type="Float32" Name="T" NumberOfComponents="1" format="ascii">
          343.15 343.15 343.15 343.15 335.65 335.65 335.65 335.65 343.15 343.15 335.65 335.65
          343.15 343.15 335.65 335.65 343.15 343.15 335.65 335.65 328.15 328.15 328.15 328.15
          328.15 328.15 328.15 328.15 328.15 328.15 320.65 320.65 320.65 320.65 320.65 320.65
          320.65 320.65 320.65 320.65 313.15 313.15 313.15 313.15 313.15 313.15 313.15 313.15
          313.15 313.15
        </DataArray>
        <DataArray type="Float32" Name="pW" NumberOfComponents="1" format="ascii">
          2.08789e+07 2.08789e+07 2.11188e+07 2.12054e+07 1.82817e+07 1.82817e+07 1.86687e+07 1.90162e+07 2.12054e+07 2.13659e+07 1.90162e+07 1.91503e+07
          2.11188e+07 2.12054e+07 1.86687e+07 1.90162e+07 2.08789e+07 2.08789e+07 1.82817e+07 1.82817e+07 1.56845e+07 1.56845e+07 1.63659e+07 1.65145e+07
          1.65145e+07 2.44988e+07 1.63659e+07 1.65145e+07 1.56845e+07 1.56845e+07 1.30873e+07 1.30873e+07 1.34366e+07 1.37793e+07 1.37793e+07 1.3974e+07
          1.34366e+07 1.37793e+07 1.30873e+07 1.30873e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07
          1.04901e+07 1.04901e+07
        </DataArray>
        <DataArray type="Float32" Name="pN" NumberOfComponents="1" format="ascii">
          2.08964e+07 2.08964e+07 2.11364e+07 2.1223e+07 1.82993e+07 1.82993e+07 1.86863e+07 1.90338e+07 2.1223e+07 2.13835e+07 1.90338e+07 1.91678e+07
          2.11364e+07 2.1223e+07 1.86863e+07 1.90338e+07 2.08964e+07 2.08964e+07 1.82993e+07 1.82993e+07 1.57021e+07 1.57021e+07 1.63835e+07 1.65321e+07
          1.65321e+07 2.45192e+07 1.63835e+07 1.65321e+07 1.57021e+07 1.57021e+07 1.31049e+07 1.31049e+07 1.34541e+07 1.37969e+07 1.37969e+07 1.39916e+07
          1.34541e+07 1.37969e+07 1.31049e+07 1.31049e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07
          1.05077e+07 1.05077e+07
        </DataArray>
        <DataArray type="Float32" Name="pC" NumberOfComponents="1" format="ascii">
          17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.8
          17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
          17569.9 20468.6 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17570.1
          17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
          17569.9 17569.9
        </DataArray>
        <DataArray type="Float32" Name="SW" NumberOfComponents="1" format="ascii">
          1 1 1 1 1 1 1 1 1 1 1 1
          1 1 1 1 1 1 1 1 1 1 1 0.999999
          0.999999 0.938081 1 0.999999 1 1 1 1 1 1 1 0.999994
          1 1 1 1 1 1 1 1 1 1 1 1
          1 1
        </DataArray>
        <DataArray type="Float32" Name="SN" NumberOfComponents="1" format="ascii">
          0 0 0 0 0 0 -2.79107e-16 -3.03227e-11 0 0 -3.03227e-11 -1.9951e-06
          0 0 -2.79107e-16 -3.03227e-11 0 0 0 0 0 0 1.26192e-11 1.16498e-06
          1.16498e-06 0.0619192 1.26192e-11 1.16498e-06 0 0 0 0 1.73792e-15 1.02144e-10 1.02144e-10 5.95419e-06
          1.73792e-15 1.02144e-10 0 0 0 0 0 0 0 0 0 0
          0 0
        </DataArray>
        <DataArray type="Float32" Name="rhoW" NumberOfComponents="1" format="ascii">
          1056.55 1056.55 1056.63 1056.67 1059.66 1059.66 1059.81 1059.94 1056.67 1056.72 1059.94 1059.99
          1056.63 1056.67 1059.81 1059.94 1056.55 1056.55 1059.66 1059.66 1062.49 1062.49 1062.75 1062.8
          1062.8 1065.82 1062.75 1062.8 1062.49 1062.49 1064.99 1064.99 1065.12 1065.26 1065.26 1065.33
          1065.12 1065.26 1064.99 1064.99 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12
          1067.12 1067.12
        </DataArray>
        <DataArray type="Float32" Name="rhoN" NumberOfComponents="1" format="ascii">
          695.741 695.741 699.765 701.181 675.345 675.345 683.412 690.278 701.181 703.787 690.278 692.844
          699.765 701.181 683.412 690.278 695.741 695.741 675.345 675.345 673.462 673.462 690.273 693.671
          693.671 806.172 690.273 693.671 673.462 673.462 669.574 669.574 681.12 691.343 691.343 696.773
          681.12 691.343 669.574 669.574 660.265 660.265 660.265 660.265 660.265 660.265 660.265 660.265
          660.265 660.265
        </DataArray>
        <DataArray type="Float32" Name="u" NumberOfComponents="3" format="ascii">
          -0.00277252 0 0 9.38795e-15 0 0 -0.0051998 -0.0051998 0 1.34598e-14 -0.00758703 0
          -0.00358451 0 0.0138498 9.54299e-15 0 0.0124885 -0.00763547 -0.00763547 0.0113756 1.64572e-14 -0.0123225 0.00747015
          -0.00758703 -1.63127e-14 0 1.52919e-14 -2.07936e-14 0 -0.0123225 -1.62525e-14 0.00747015 2.01394e-14 -2.36144e-14 -0.0187165
          -0.0051998 0.0051998 0 1.37546e-14 0.00758703 0 -0.00763547 0.00763547 0.0113756 1.96227e-14 0.0123225 0.00747015
          -0.00277252 0 0 1.0583e-14 0 0 -0.00358451 0 0.0138498 1.20176e-14 0 0.0124885
          -0.00482824 0 0.0322354 5.01065e-15 0 0.0330174 -0.0114511 -0.0114511 0.0332345 7.97628e-15 -0.0385885 0.0340769
          -0.0385885 -1.72168e-14 0.0340769 1.30699e-14 -2.25566e-14 0.0347518 -0.0114511 0.0114511 0.0332345 2.02051e-14 0.0385885 0.0340769
          -0.00482824 0 0.0322354 1.04974e-14 0 0.0330174 -0.00262126 0 0.0482306 -1.75082e-15 0 0.0521098
          -0.00589184 -0.00589184 0.0548623 -6.46787e-15 -0.0103042 0.0622795 -0.0103042 -2.6559e-14 0.0622795 -7.66454e-15 -2.84865e-14 0.0930209
          -0.00589184 0.00589184 0.0548623 2.06746e-15 0.0103042 0.0622795 -0.00262126 0 0.0482306 2.90755e-15 0 0.0521098
          -0.00333445 0 0.0545541 -3.54304e-15 0 0.0583753 -0.0050437 -0.0050437 0.062298 -5.40032e-15 -0.00717425 0.0708312
          -0.00717425 -2.87929e-14 0.0708312 -4.69979e-15 -2.90889e-14 0.0801422 -0.0050437 0.0050437 0.062298 1.34065e-15 0.00717425 0.0708312
          -0.00333445 0 0.0545541 3.05686e-15 0 0.0583753
        </DataArray>
      </PointData>
      <CellData Scalars="principal stress 1" Vectors="effective stress changes X">
        <DataArray type="Float32" Name="effective stress changes X" NumberOfComponents="3" format="ascii">
          -148681 33881.3 21285.8 -285539 35371.4 93158.1 -285539 -35371.4 93158.1 -148681 -33881.3 21285.8
          -212300 106246 21753.9 -797547 159122 146641 -797547 -159122 146641 -212300 -106246 21753.9
          -196420 103608 -51717.6 -779134 157749 -183799 -779134 -157749 -183799 -196420 -103608 -51717.6
          -112879 28806.9 -59464.3 -256217 32714.7 -149952 -256217 -32714.7 -149952 -112879 -28806.9 -59464.3
        </DataArray>
        <DataArray type="Float32" Name="effective stress changes Y" NumberOfComponents="3" format="ascii">
          33881.3 111005 36659.3 35371.4 -285539 93158.1 -35371.4 -285539 -93158.1 -33881.3 111005 -36659.3
          106246 275185 88789 159122 -797547 146641 -159122 -797547 -146641 -106246 275185 -88789
          103608 258720 -131759 157749 -779134 -183799 -157749 -779134 183799 -103608 258720 131759
          28806.9 113647 -102448 32714.7 -256217 -149952 -32714.7 -256217 149952 -28806.9 113647 102448
        </DataArray>
        <DataArray type="Float32" Name="effective stress changes Z" NumberOfComponents="3" format="ascii">
          21285.8 36659.3 -278639 93158.1 93158.1 -159812 93158.1 -93158.1 -159812 21285.8 -36659.3 -278639
          21753.9 88789 -511703 146641 146641 -1.09026e+06 146641 -146641 -1.09026e+06 21753.9 -88789 -511703
          -51717.6 -131759 -497048 -183799 -183799 -1.12947e+06 -183799 183799 -1.12947e+06 -51717.6 131759 -497048
          -59464.3 -102448 -171305 -149952 -149952 -172449 -149952 149952 -172449 -59464.3 102448 -171305
        </DataArray>
        <DataArray type="Float32" Name="principal stress 1" NumberOfComponents="1" format="ascii">
          4.32126e+07 4.32438e+07 4.32438e+07 4.32126e+07 3.77761e+07 3.81095e+07 3.81095e+07 3.77761e+07 3.20309e+07 3.23961e+07 3.23961e+07 3.20309e+07
          2.59781e+07 2.58894e+07 2.58894e+07 2.59781e+07
        </DataArray>
        <DataArray type="Float32" Name="principal stress 2" NumberOfComponents="1" format="ascii">
          4.29529e+07 4.31181e+07 4.31181e+07 4.29529e+07 3.72886e+07 3.81095e+07 3.81095e+07 3.72886e+07 3.15757e+07 3.23961e+07 3.23961e+07 3.15757e+07
          2.57516e+07 2.58056e+07 2.58056e+07 2.57516e+07
        </DataArray>
        <DataArray type="Float32" Name="principal stress 3" NumberOfComponents="1" format="ascii">
          4.28229e+07 4.31181e+07 4.31181e+07 4.28229e+07 3.69892e+07 3.78168e+07 3.78168e+07 3.69892e+07 3.12751e+07 3.20457e+07 3.20457e+07 3.12751e+07
          2.56932e+07 2.58056e+07 2.58056e+07 2.56932e+07
        </DataArray>
        <DataArray type="Float32" Name="total stresses X" NumberOfComponents="3" format="ascii">
          4.29529e+07 33881.3 21285.8 4.31181e+07 35371.4 93158.1 4.31181e+07 -35371.4 93158.1 4.29529e+07 -33881.3 21285.8
          3.72886e+07 106246 21753.9 3.81095e+07 159122 146641 3.81095e+07 -159122 146641 3.72886e+07 -106246 21753.9
          3.15757e+07 103608 -51717.6 3.23961e+07 157749 -183799 3.23961e+07 -157749 -183799 3.15757e+07 -103608 -51717.6
          2.57516e+07 28806.9 -59464.3 2.58056e+07 32714.7 -149952 2.58056e+07 -32714.7 -149952 2.57516e+07 -28806.9 -59464.3
        </DataArray>
        <DataArray type="Float32" Name="total stresses Y" NumberOfComponents="3" format="ascii">
          33881.3 4.32126e+07 36659.3 35371.4 4.31181e+07 93158.1 -35371.4 4.31181e+07 -93158.1 -33881.3 4.32126e+07 -36659.3
          106246 3.77761e+07 88789 159122 3.81095e+07 146641 -159122 3.81095e+07 -146641 -106246 3.77761e+07 -88789
          103608 3.20309e+07 -131759 157749 3.23961e+07 -183799 -157749 3.23961e+07 183799 -103608 3.20309e+07 131759
          28806.9 2.59781e+07 -102448 32714.7 2.58056e+07 -149952 -32714.7 2.58056e+07 149952 -28806.9 2.59781e+07 102448
        </DataArray>
        <DataArray type="Float32" Name="total stresses Z" NumberOfComponents="3" format="ascii">
          21285.8 36659.3 4.28229e+07 93158.1 93158.1 4.32438e+07 93158.1 -93158.1 4.32438e+07 21285.8 -36659.3 4.28229e+07
          21753.9 88789 3.69892e+07 146641 146641 3.78168e+07 146641 -146641 3.78168e+07 21753.9 -88789 3.69892e+07
          -51717.6 -131759 3.12751e+07 -183799 -183799 3.20457e+07 -183799 183799 3.20457e+07 -51717.6 131759 3.12751e+07
          -59464.3 -102448 2.56932e+07 -149952 -149952 2.58894e+07 -149952 149952 2.58894e+07 -59464.3 102448 2.56932e+07
        </DataArray>
        <DataArray type="Float32" Name="initial stresses X" NumberOfComponents="3" format="ascii">
          4.28905e+07 0 0 4.28905e+07 0 0 4.28905e+07 0 0 4.28905e+07 0 0
          3.71718e+07 0 0 3.71718e+07 0 0 3.71718e+07 0 0 3.71718e+07 0 0
          3.14531e+07 0 0 3.14531e+07 0 0 3.14531e+07 0 0 3.14531e+07 0 0
          2.57343e+07 0 0 2.57343e+07 0 0 2.57343e+07 0 0 2.57343e+07 0 0
        </DataArray>
        <DataArray type="Float32" Name="initial stresses Y" NumberOfComponents="3" format="ascii">
          0 4.28905e+07 0 0 4.28905e+07 0 0 4.28905e+07 0 0 4.28905e+07 0
          0 3.71718e+07 0 0 3.71718e+07 0 0 3.71718e+07 0 0 3.71718e+07 0
          0 3.14531e+07 0 0 3.14531e+07 0 0 3.14531e+07 0 0 3.14531e+07 0
          0 2.57343e+07 0 0 2.57343e+07 0 0 2.57343e+07 0 0 2.57343e+07 0
        </DataArray>
        <DataArray type="Float32" Name="initial stresses Z" NumberOfComponents="3" format="ascii">
          0 0 4.28905e+07 0 0 4.28905e+07 0 0 4.28905e+07 0 0 4.28905e+07
          0 0 3.71718e+07 0 0 3.71718e+07 0 0 3.71718e+07 0 0 3.71718e+07
          0 0 3.14531e+07 0 0 3.14531e+07 0 0 3.14531e+07 0 0 3.14531e+07
          0 0 2.57343e+07 0 0 2.57343e+07 0 0 2.57343e+07 0 0 2.57343e+07
        </DataArray>
        <DataArray type="Float32" Name="delta pEff" NumberOfComponents="1" format="ascii">
          211012 513102 513102 211012 329132 1.73522e+06 1.73522e+06 329132 319095 1.72213e+06 1.72213e+06 319095
          130161 327499 327499 130161
        </DataArray>
        <DataArray type="Float32" Name="effectivePressure" NumberOfComponents="1" format="ascii">
          1.97913e+07 2.00934e+07 2.00934e+07 1.97913e+07 1.73122e+07 1.87183e+07 1.87183e+07 1.73122e+07 1.4705e+07 1.6108e+07 1.6108e+07 1.4705e+07
          1.19188e+07 1.21162e+07 1.21162e+07 1.19188e+07
        </DataArray>
        <DataArray type="Float32" Name="Pcr_tensile" NumberOfComponents="1" format="ascii">
          -2.30316e+07 -2.30247e+07 -2.30247e+07 -2.30316e+07 -1.9677e+07 -1.90985e+07 -1.90985e+07 -1.9677e+07 -1.65701e+07 -1.59377e+07 -1.59377e+07 -1.65701e+07
          -1.37743e+07 -1.36894e+07 -1.36894e+07 -1.37743e+07
        </DataArray>
        <DataArray type="Float32" Name="Pcr_shear" NumberOfComponents="1" format="ascii">
          -2.28368e+07 -2.29619e+07 -2.29619e+07 -2.28368e+07 -1.92836e+07 -1.89521e+07 -1.89521e+07 -1.92836e+07 -1.61922e+07 -1.57625e+07 -1.57625e+07 -1.61922e+07
          -1.36319e+07 -1.36475e+07 -1.36475e+07 -1.36319e+07
        </DataArray>
        <DataArray type="Float32" Name="effective Kxx" NumberOfComponents="1" format="ascii">
          1.00281e-14 1.00651e-14 1.00651e-14 1.00281e-14 1.004e-14 1.02414e-14 1.02414e-14 1.004e-14 1.00387e-14 1.02416e-14 1.02416e-14 1.00387e-14
          1.00152e-14 1.0061e-14 1.0061e-14 1.00152e-14
        </DataArray>
        <DataArray type="Float32" Name="effective Porosity" NumberOfComponents="1" format="ascii">
          0.200025 0.200058 0.200058 0.200025 0.200036 0.200215 0.200215 0.200036 0.200035 0.200215 0.200215 0.200035
          0.200014 0.200055 0.200055 0.200014
        </DataArray>
      </CellData>
      <Points>
        <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
          250 0 0 500 0 0 250 250 0 500 250 0
          250 0 250 500 0 250 250 250 250 500 250 250
          250 500 0 500 500 0 250 500 250 500 500 250
          250 750 0 500 750 0 250 750 250 500 750 250
          250 1000 0 500 1000 0 250 1000 250 500 1000 250
          250 0 500 500 0 500 250 250 500 500 250 500
          250 500 500 500 500 500 250 750 500 500 750 500
          250 1000 500 500 1000 500 250 0 750 500 0 750
          250 250 750 500 250 750 250 500 750 500 500 750
          250 750 750 500 750 750 250 1000 750 500 1000 750
          250 0 1000 500 0 1000 250 250 1000 500 250 1000
          250 500 1000 500 500 1000 250 750 1000 500 750 1000
          250 1000 1000 500 1000 1000
        </DataArray>
      </Points>
      <Cells>
        <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
          0 1 3 2 4 5 7 6 2 3 9 8
          6 7 11 10 8 9 13 12 10 11 15 14
          12 13 17 16 14 15 19 18 4 5 7 6
          20 21 23 22 6 7 11 10 22 23 25 24
          10 11 15 14 24 25 27 26 14 15 19 18
          26 27 29 28 20 21 23 22 30 31 33 32
          22 23 25 24 32 33 35 34 24 25 27 26
          34 35 37 36 26 27 29 28 36 37 39 38
          30 31 33 32 40 41 43 42 32 33 35 34
          42 43 45 44 34 35 37 36 44 45 47 46
          36 37 39 38 46 47 49 48
        </DataArray>
        <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
          8 16 24 32 40 48 56 64 72 80 88 96
          104 112 120 128
        </DataArray>
        <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
          12 12 12 12 12 12 12 12 12 12 12 12
          12 12 12 12
        </DataArray>
      </Cells>
    </Piece>
  </UnstructuredGrid>
</VTKFile>