Skip to content
GitLab
Commit b74229a4 authored by Bernd Flemisch's avatar Bernd Flemisch
Browse files

Merge branch 'fix/1p-model' into 'master' 

Fix/1p model

* Remove permeability output (make 1p work with tensor permeabilities)
* Only compute velocity when needed for output

See merge request !136
parents e7d5507b a395cd04
Expand all Hide whitespace changes
Inline Side-by-side
dumux/porousmediumflow/1p/implicit/model.hh
View file @ b74229a4
......@@ -85,7 +85,6 @@ public:
// create the required scalar fields
unsigned numDofs = this->numDofs();
ScalarField *p = writer.allocateManagedBuffer(numDofs);
ScalarField *K = writer.allocateManagedBuffer(numDofs);
VectorField *velocity = writer.template allocateManagedBuffer<double, dimWorld>(numDofs);
ImplicitVelocityOutput<TypeTag> velocityOutput(this->problem_());
......@@ -120,22 +119,18 @@ public:
for (int scvIdx = 0; scvIdx < fvGeometry.numScv; ++scvIdx)
{
int dofIdxGlobal = this->dofMapper().subIndex(element, scvIdx, dofCodim);
const SpatialParams &spatialParams = this->problem_().spatialParams();
(*p)[dofIdxGlobal] = elemVolVars[scvIdx].pressure();
(*K)[dofIdxGlobal] = spatialParams.intrinsicPermeability(element,
fvGeometry,
scvIdx);
}
// velocity output
velocityOutput.calculateVelocity(*velocity, elemVolVars, fvGeometry, element, /*phaseIdx=*/0);
if (velocityOutput.enableOutput())
{
velocityOutput.calculateVelocity(*velocity, elemVolVars, fvGeometry, element, /*phaseIdx=*/0);
}
}
}
writer.attachDofData(*p, "p", isBox);
writer.attachDofData(*K, "K", isBox);
if (velocityOutput.enableOutput())
{
writer.attachDofData(*velocity, "velocity", isBox, dim);
......
test/references/1pniboxconduction-reference.vtu
View file @ b74229a4
......@@ -39,42 +39,6 @@
100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
100000 100000 100000 100000 100000 100000
</DataArray>
<DataArray type="Float32" Name="K" NumberOfComponents="1" format="ascii">
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
</DataArray>
<DataArray type="Float32" Name="velocity" NumberOfComponents="3" format="ascii">
-1.17342e-09 8.52243e-19 0 -1.16807e-09 1.68217e-18 0 -1.17342e-09 8.52243e-19 0 -1.16807e-09 1.68217e-18 0
-1.15247e-09 3.1127e-18 0 -1.15247e-09 3.1127e-18 0 -1.1275e-09 1.84323e-18 0 -1.1275e-09 1.84323e-18 0
......
test/references/1pniboxconvection-reference.vtu
View file @ b74229a4
......@@ -19,22 +19,6 @@
102168 102168 101897 101897 101626 101626 101355 101355 101084 101084 100813 100813
100542 100542 100271 100271 100000 100000
</DataArray>
<DataArray type="Float32" Name="K" NumberOfComponents="1" format="ascii">
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
</DataArray>
<DataArray type="Float32" Name="velocity" NumberOfComponents="3" format="ascii">
0.0001 -2.06574e-18 0 0.0001 0 0 0.0001 -2.06574e-18 0 0.0001 0 0
0.0001 1.37713e-18 0 0.0001 1.37713e-18 0 0.0001 8.60723e-19 0 0.0001 8.60723e-19 0
......
test/references/1pniccconduction-reference.vtu
View file @ b74229a4
......@@ -22,25 +22,6 @@
100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
100000 100000 100000 100000 100000 100000 100000 100000
</DataArray>
<DataArray type="Float32" Name="K" NumberOfComponents="1" format="ascii">
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
</DataArray>
<DataArray type="Float32" Name="velocity" NumberOfComponents="3" format="ascii">
-1.17311e-09 0 0 -1.1625e-09 0 0 -1.14204e-09 0 0 -1.11268e-09 0 0
-1.07545e-09 0 0 -1.03153e-09 0 0 -9.82132e-10 0 0 -9.28491e-10 0 0
......
test/references/1pniccconvection-reference.vtu
View file @ b74229a4
......@@ -12,15 +12,6 @@
105284 105013 104742 104471 104200 103929 103658 103387 103116 102845 102574 102303
102032 101761 101490 101219 100948 100677 100406 100135
</DataArray>
<DataArray type="Float32" Name="K" NumberOfComponents="1" format="ascii">
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
</DataArray>
<DataArray type="Float32" Name="velocity" NumberOfComponents="3" format="ascii">
0.0001 0 0 0.0001 0 0 0.0001 0 0 0.0001 0 0
0.0001 0 0 0.0001 0 0 0.0001 0 0 0.0001 0 0
......
test/references/1ptestbox-reference.vtu
View file @ b74229a4
......@@ -16,19 +16,6 @@
108419 108793 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
100000
</DataArray>
<DataArray type="Float32" Name="K" NumberOfComponents="1" format="ascii">
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-12 1e-12 1e-12 1e-12 1e-12 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-12
1e-12 1e-12 1e-12 1e-12 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-12 1e-12
1e-12 1e-12 1e-12 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-12 1e-12 1e-12
1e-12 1e-12 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-12 1e-12 1e-12 1e-12
1e-12 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10
</DataArray>
</PointData>
<CellData Scalars="process rank">
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
......
test/references/1ptestbox1d3d-reference.vtp
View file @ b74229a4
......@@ -6,9 +6,6 @@
<DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
200000 182416 150000 117584 100000
</DataArray>
<DataArray type="Float32" Name="K" NumberOfComponents="1" format="ascii">
1e-10 1e-12 1e-12 1e-12 1e-10
</DataArray>
</PointData>
<CellData Scalars="process rank">
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
......
test/references/1ptestbox2d3d-reference.vtu
View file @ b74229a4
......@@ -6,9 +6,6 @@
<DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
200000 200000 150000 150000 200000 150000 100000 100000 100000
</DataArray>
<DataArray type="Float32" Name="K" NumberOfComponents="1" format="ascii">
1e-10 1e-10 1e-12 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
</DataArray>
</PointData>
<CellData Scalars="process rank">
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
......
test/references/1ptestcc-reference.vtu
View file @ b74229a4
......@@ -14,17 +14,6 @@
103161 103161 104019 106256 111315 112828 104090 103532 102255 101487 101162 101162
101487 102255 103532 104090
</DataArray>
<DataArray type="Float32" Name="K" NumberOfComponents="1" format="ascii">
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-12 1e-12
1e-12 1e-12 1e-12 1e-12 1e-10 1e-10 1e-10 1e-10 1e-12 1e-12 1e-12 1e-12
1e-12 1e-12 1e-10 1e-10 1e-10 1e-10 1e-12 1e-12 1e-12 1e-12 1e-12 1e-12
1e-10 1e-10 1e-10 1e-10 1e-12 1e-12 1e-12 1e-12 1e-12 1e-12 1e-10 1e-10
1e-10 1e-10 1e-12 1e-12 1e-12 1e-12 1e-12 1e-12 1e-10 1e-10 1e-10 1e-10
1e-12 1e-12 1e-12 1e-12 1e-12 1e-12 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
1e-10 1e-10 1e-10 1e-10
</DataArray>
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
......
test/references/1ptestcc1d3d-reference.vtp
View file @ b74229a4
......@@ -6,9 +6,6 @@
<DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
198550 173550 126449 101449
</DataArray>
<DataArray type="Float32" Name="K" NumberOfComponents="1" format="ascii">
1e-10 1e-12 1e-12 1e-10
</DataArray>
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
0 0 0 0
</DataArray>
......
test/references/1ptestcc2d3d-reference.vtu
View file @ b74229a4
......@@ -6,9 +6,6 @@
<DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
180435 167391 180435 167391 132609 119565 132609 119565
</DataArray>
<DataArray type="Float32" Name="K" NumberOfComponents="1" format="ascii">
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
</DataArray>
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
0 0 0 0 0 0 0 0
</DataArray>
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment