Commit 6b904e8a authored by Katharina Heck's avatar Katharina Heck
Browse files

update example 4

parent 861b8bca
......@@ -205,7 +205,9 @@ public:
const auto radiation = Radiation::radationEquilibrium(element, elemVolVars[scvf.insideScvIdx()], scvf, couplingManager().couplingData().stokesTemperature(element, scvf), time);
evaporation += flux[0] * scvf.area() * elemVolVars[scvf.insideScvIdx()].extrusionFactor();
const auto& globalPos = scvf.center();
if (globalPos[0] > 0.0)
evaporation += flux[0] * scvf.area() * elemVolVars[scvf.insideScvIdx()].extrusionFactor();
energyFlux += fluxEnergy* scvf.area() * elemVolVars[scvf.insideScvIdx()].extrusionFactor();
averageTemperature += elemVolVars[scvf.insideScvIdx()].temperature();
......
......@@ -29,10 +29,10 @@ InitializationTime = 0.01
Name = darcy
Pressure = 1.0e5
Saturation = 0.7567470 # initial Sw
Temperature = 291.55 # [K] 17 °C
Temperature = 291.55 # [K] 18.4 °C
InitPhasePresence = 3 # bothPhases
InitializationTime = 0.01
PlexiglassThermalConductivity = 0.184 # [W/(m*K)] 0.116
PlexiglassThermalConductivity = 2.8 # [W/(m*K)] 0.116
PlexiglassThickness = 0.005 # [m]
[SpatialParams]
......@@ -58,7 +58,7 @@ WriteFaceData = false
[Component]
SolidDensity = 1650
SolidThermalConductivity = 1.8
SolidThermalConductivity = 2.8
SolidHeatCapacity = 1381 #engineering toolbox
[RANS]
......
......@@ -79,7 +79,7 @@ public:
params_.setVgn(getParam<Scalar>("SpatialParams.VgN"));
upperLayerParams_.setSwr(getParam<Scalar>("SpatialParams.Swr"));
upperLayerParams_.setSnr(0.2);
upperLayerParams_.setSnr(getParam<Scalar>("SpatialParams.Snr"));
// parameters for the vanGenuchten law
upperLayerParams_.setVgAlpha(7.54332e-4);
upperLayerParams_.setVgn(1.6);
......@@ -95,7 +95,7 @@ public:
PermeabilityType permeabilityAtPos(const GlobalPosition& globalPos) const
{
if (globalPos[1]>0.2)
return 1.1e-12;
return permeability_;
else
return permeability_;
}
......@@ -122,7 +122,7 @@ public:
*/
const MaterialLawParams& materialLawParamsAtPos(const GlobalPosition& globalPos) const
{
if (globalPos[1]>0.2)
if (globalPos[1]>0.22)
return params_;
else
return params_;
......
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