Commit dc01d86f authored by Katharina Heck's avatar Katharina Heck
Browse files

update example 4 with deeper darcy domain

parent 6b904e8a
......@@ -44,7 +44,7 @@ AlphaBJ = 1.0
Swr = 0.057
Snr = 0.029
VgAlpha = 4.28e-5
VgN = 1.5
VgN = 1.32
[Problem]
Name = example1
......
[TimeLoop]
DtInitial = 0.1 # [s]
MaxTimeStepSize = 1800 # [s] (12 hours)
TEnd = 864000 # [s] (15 hours)
TEnd = 86400 # [s] (15 hours)
EpisodeLength = 1800
[Stokes.Grid]
......
......@@ -331,7 +331,8 @@ public:
FluidState fluidState;
fluidState.setPressure(0, 1e5);
fluidState.setTemperature(refTemperature());
fluidState.setMoleFraction(0, 0, 1);
fluidState.setMoleFraction(0, 0, 1-refMoleFrac());
fluidState.setMoleFraction(0,1,refMoleFrac());
Scalar density = FluidSystem::density(fluidState, 0);
PrimaryVariables values(0.0);
......
......@@ -3,5 +3,3 @@ add_input_file_links()
dune_add_test(NAME test_stokes1p2cdarcy2p2c_radiation_lysimeter
SOURCES main.cc
CMAKE_GUARD HAVE_UMFPACK)
set(CMAKE_BUILD_TYPE Debug)
......@@ -418,10 +418,10 @@ public:
values[pressureIdx] = pressure_ + 1. * this->spatialParams().gravity(globalPos)[1] * (globalPos[1] - this->gridGeometry().bBoxMax()[1]);
values[switchIdx] = initialSw_;
values[energyEqIdx] = temperature_; //20
// if (globalPos[1]>0.24)
// if (globalPos[1]>0.49)
// {
// values[energyEqIdx] = temperature_-4;
// values[switchIdx] = 0.6;
// values[switchIdx] = 0.7;
// }
return values;
}
......
......@@ -6,7 +6,7 @@ EpisodeLength = 3600
[Stokes.Grid]
Positions0 = -4 0.0 1.0 1.25
Positions1 = 0.25 0.5
Positions1 = 0.5 0.75
Grading0 = 1.0 1.0 1.0
Grading1 = 1.1
Cells0 = 40 25 10
......@@ -14,7 +14,7 @@ Cells1 = 40
[Darcy.Grid]
Positions0 = 0.0 1.0
Positions1 = 0.0 0.25
Positions1 = 0.0 0.5
Cells0 = 25
Cells1 = 40
Grading0 = 1.0
......@@ -29,11 +29,11 @@ InitializationTime = 0.01
Name = darcy
Pressure = 1.0e5
Saturation = 0.7567470 # initial Sw
Temperature = 291.55 # [K] 18.4 °C
Temperature = 293.15 # [K] 20 °C
InitPhasePresence = 3 # bothPhases
InitializationTime = 0.01
PlexiglassThermalConductivity = 2.8 # [W/(m*K)] 0.116
PlexiglassThickness = 0.005 # [m]
PlexiglassThermalConductivity = 50.2 # [W/(m*K)] steel
PlexiglassThickness = 0.01 # [m]
[SpatialParams]
Porosity = 0.376666666666668
......@@ -58,7 +58,7 @@ WriteFaceData = false
[Component]
SolidDensity = 1650
SolidThermalConductivity = 2.8
SolidThermalConductivity = 5.3
SolidHeatCapacity = 1381 #engineering toolbox
[RANS]
......
......@@ -84,6 +84,12 @@ public:
upperLayerParams_.setVgAlpha(7.54332e-4);
upperLayerParams_.setVgn(1.6);
lowerParams_.setSwr(0.14);
lowerParams_.setSnr(0.073);
// parameters for the vanGenuchten law
lowerParams_.setVgAlpha(6.6e-5);
lowerParams_.setVgn(1.59);
}
/*!
......@@ -94,8 +100,8 @@ public:
*/
PermeabilityType permeabilityAtPos(const GlobalPosition& globalPos) const
{
if (globalPos[1]>0.2)
return permeability_;
if (globalPos[1]<0.26)
return 7.7e-14;
else
return permeability_;
}
......@@ -105,7 +111,12 @@ public:
* \param globalPos The global position
*/
Scalar porosityAtPos(const GlobalPosition& globalPos) const
{ return porosity_; }
{
if (globalPos[1]<0.26)
return 0.398;
else
return porosity_;
}
/*! \brief Define the Beavers-Joseph coefficient in [-].
*
......@@ -122,8 +133,11 @@ public:
*/
const MaterialLawParams& materialLawParamsAtPos(const GlobalPosition& globalPos) const
{
if (globalPos[1]>0.22)
if (globalPos[1]>0.47)
return params_;
else if
(globalPos[1]<0.26)
return lowerParams_;
else
return params_;
}
......@@ -144,6 +158,7 @@ private:
Scalar porosity_;
Scalar alphaBJ_;
MaterialLawParams params_;
MaterialLawParams lowerParams_;
MaterialLawParams upperLayerParams_;
static constexpr Scalar eps_ = 1.0e-7;
......
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