Skip to content
Snippets Groups Projects
Commit 9d811602 authored by Timo Koch's avatar Timo Koch
Browse files

[test][generalproblem] Use vector for lens

parent f4edafde
No related branches found
No related tags found
2 merge requests!31Feature/colebrookwhiteboundarylayer,!19Fix/generallensproblem
...@@ -91,10 +91,8 @@ public: ...@@ -91,10 +91,8 @@ public:
GeneralLensSpatialParams(const GridView& gridView) GeneralLensSpatialParams(const GridView& gridView)
: ParentType(gridView), eps_(3e-6) : ParentType(gridView), eps_(3e-6)
{ {
lensLowerLeft_[0] = GET_RUNTIME_PARAM(TypeTag, Scalar, SpatialParams.LensLowerLeftX); lensLowerLeft_ = GET_RUNTIME_PARAM(TypeTag, GlobalPosition, SpatialParams.LensLowerLeft);
lensLowerLeft_[1] = GET_RUNTIME_PARAM(TypeTag, Scalar, SpatialParams.LensLowerLeftY); lensUpperRight_ = GET_RUNTIME_PARAM(TypeTag, GlobalPosition, SpatialParams.LensUpperRight);
lensUpperRight_[0] = GET_RUNTIME_PARAM(TypeTag, Scalar, SpatialParams.LensUpperRightX);
lensUpperRight_[1] = GET_RUNTIME_PARAM(TypeTag, Scalar, SpatialParams.LensUpperRightY);
// residual saturations // residual saturations
lensMaterialParams_.setSwr(0.18); lensMaterialParams_.setSwr(0.18);
...@@ -109,10 +107,12 @@ public: ...@@ -109,10 +107,12 @@ public:
outerMaterialParams_.setVgAlpha(0.0037); outerMaterialParams_.setVgAlpha(0.0037);
outerMaterialParams_.setVgn(4.7); outerMaterialParams_.setVgn(4.7);
// initialize with zero
lensK_ = 0.0; outerK_ = 0.0;
for (int i=0; i < dim; i++) for (int i=0; i < dim; i++)
{ {
lensK_[i][i] = 9.05e-12; lensK_[i][i] = 9.05e-12;
outerK_[i][i] = 4.6e-10; outerK_[i][i] = 4.6e-10;
} }
} }
......
...@@ -7,10 +7,8 @@ UpperRight = 6 4 # upper right corner coordinates [m] ...@@ -7,10 +7,8 @@ UpperRight = 6 4 # upper right corner coordinates [m]
Cells = 48 32 # number of cells in (x, y)-direction [-] Cells = 48 32 # number of cells in (x, y)-direction [-]
[SpatialParams] [SpatialParams]
LensLowerLeftX = 1.0 # [m] x-coordinate of the lower left lens corner LensLowerLeft = 1.0 2.0 # coordinates of the lower left lens corner [m]
LensLowerLeftY = 2.0 # [m] y-coordinate of the lower left lens corner LensUpperRight = 4.0 3.0 # coordinates of the upper right lens corner [m]
LensUpperRightX = 4.0 # [m] x-coordinate of the upper right lens corner
LensUpperRightY = 3.0 # [m] y-coordinate of the upper right lens corner
[Vtk] [Vtk]
OutputLevel = 1 OutputLevel = 1
......
...@@ -7,10 +7,8 @@ UpperRight = 6 4 # upper right corner coordinates [m] ...@@ -7,10 +7,8 @@ UpperRight = 6 4 # upper right corner coordinates [m]
Cells = 48 32 # number of cells in (x, y)-direction [-] Cells = 48 32 # number of cells in (x, y)-direction [-]
[SpatialParams] [SpatialParams]
LensLowerLeftX = 1.0 # [m] x-coordinate of the lower left lens corner LensLowerLeft = 1.0 2.0 # coordinates of the lower left lens corner [m]
LensLowerLeftY = 2.0 # [m] y-coordinate of the lower left lens corner LensUpperRight = 4.0 3.0 # coordinates of the upper right lens corner [m]
LensUpperRightX = 4.0 # [m] x-coordinate of the upper right lens corner
LensUpperRightY = 3.0 # [m] y-coordinate of the upper right lens corner
[Vtk] [Vtk]
OutputLevel = 1 OutputLevel = 1
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment