temp

test/multidomain/boundary/ransdarcy/1p2c_2p2c/problem_darcy.hh

@@ -161,7 +161,7 @@ public:
         initialWaterContent_ = evaluateWaterMassStorageTerm(curSol, gridVariables);
         lastWaterMass_ = initialWaterContent_;
 
-        storageFileName_ = "storage_" + ransName + ".csv";
+        storageFileName_ = "storage_" + problemName_ + "_" + ransName + ".csv";
         storageFile_.open(storageFileName_, std::ios::app);
         storageFile_ << "#Time[s]" << ";"
                      << "WaterMass[kg]" << ";"

test/multidomain/boundary/ransdarcy/1p2c_2p2c/scripts/evapplot.gp

-reset
-set term pngcairo size 800,600 solid
-set output "./EvaporationRate.png"
-set datafile separator ';'
-DATA='./'
-set ylabel "Evaporation Rate [mm/d]"
-set xlabel "time [days]"
-set yrange [0:7]
-set xrange [0:8]
-set key right top samplen 1
-
-plot DATA.'storage_KEpsilon.csv'  u ($1/86400):4  w l lc rgb "red"    t 'KEpsilon',\
- DATA.'storage_LowReKEpsilon.csv' u ($1/86400):4  w l lc rgb "orange" t 'Low Re KEpsilon',\
- DATA.'storage_KOmega.csv'        u ($1/86400):4  w l lc rgb "green"  t 'KOmega',\
- DATA.'storage_OneEq.csv'         u ($1/86400):4  w l lc rgb "cyan"   t 'OneEq',\
- DATA.'storage_ZeroEq.csv'        u ($1/86400):4  w l lc rgb "purple" t 'ZeroEq'

test/multidomain/boundary/ransdarcy/1p2c_2p2c/scripts/params_evaporation.input

 [TimeLoop]
 DtInitial = 1e-4 # [s]
-MaxTimeStepSize = 3600 # [s] 1 hr
+MaxTimeStepSize = 350 # [s]
 TEnd = 691200 # [s] 8 days
 
 [RANS.Grid]
 Verbosity = true
-Positions0 = -1.0 0.0 1.0 1.5
+Positions0 = -1.0 0.0 1.0 2.0
 Positions1 = 0.5 1.0
 Grading0 = -1.1 1.0 1.1
 Grading1 = 1.3
-Cells0 = 6 20 4
+Cells0 = 15 30 15
 Cells1 = 20
 
 [Darcy.Grid]
 Verbosity = true
 Positions0 = 0.0 1.0
 Positions1 = 0.0 0.5
-Cells0 = 20
+Cells0 = 30
 Cells1 = 20
 Grading0 = 1.0
 Grading1 = -1.3
@@ -33,7 +33,7 @@ Name = darcy
 Pressure = 1e5 # [Pa]
 Temperature = 298.15 # [K]
 InitPhasePresence = 3 # bothPhases
-InitialSaturation = 0.85 # [-]
+InitialSaturation = 0.90 # [-]
 
 [Darcy.SpatialParams]
 Porosity = 0.41
@@ -51,6 +51,8 @@ InterfaceDiffusionCoefficientAvg = FreeFlowOnly # Harmonic
 [RANS.RANS]
 UseStoredEddyViscosity = false
 EddyViscosityModel = "baldwinLomax"
+TurbulentSchmidtNumber = 0.7
+TurbulentPrandtlNumber = 0.85
 
 [RANS.KEpsilon]
 YPlusThreshold = 50. # should be small (10-30) for coarse grids large (30-60) for fine grids
@@ -58,7 +60,7 @@ YPlusThreshold = 50. # should be small (10-30) for coarse grids large (30-60) fo
 [Vtk]
 AddVelocity = true
 WriteFaceData = false
-OutputName = rectangle_interface
+OutputName = flat_interface
 
 [Component]
 SolidDensity = 2700
@@ -66,9 +68,9 @@ SolidThermalConductivity = 2.8
 SolidHeatCapacity = 790
 
 [Newton]
-TargetSteps = 5
-MaxSteps = 10
-MaxRelativeShift = 1e-6
+TargetSteps = 6
+MaxSteps = 8
+MaxRelativeShift = 1e-7
 
 [Assembly]
 NumericDifferenceMethod = 0

test/multidomain/boundary/ransdarcy/1p2c_2p2c/scripts/runevaporation.sh

-runSim () {
-./$1 $INPUT -Vtk.OutputName $2 | tee -a logfile.out
-}
-
-### Evaporation Tests (Drying)
-INPUT=scripts/params_evaporation.input
-
-runSim test_md_boundary_darcy2p2cni_rans1p2cnikepsilon       evaporation_kepsilon
-runSim test_md_boundary_darcy2p2cni_rans1p2cnikomega         evaporation_komega
-runSim test_md_boundary_darcy2p2cni_rans1p2cnilowrekepsilon  evaporation_lowrekepsilon
-runSim test_md_boundary_darcy2p2cni_rans1p2cnioneeq          evaporation_oneeq
-runSim test_md_boundary_darcy2p2cni_rans1p2cnizeroeq         evaporation_zeroeq
-
-gnuplot scripts/evapplot.gp

test/multidomain/boundary/ransdarcy/1p2c_2p2c/scripts/runevaporationturbulencemodels.sh

 runSim () {
-./$1 $INPUT -Vtk.OutputName $2 | tee -a logfile.out
+./$1 $INPUT -Vtk.OutputName $2 -RANS.Problem.InletVelocity $3| tee -a logfile.out
 }
 
 ### Evaporation Tests (Drying)
 INPUT=scripts/params_evaporation.input
 
-runSim test_md_boundary_darcy2p2cni_rans1p2cnikepsilon       evaporation_kepsilon
-runSim test_md_boundary_darcy2p2cni_rans1p2cnikomega         evaporation_komega
-runSim test_md_boundary_darcy2p2cni_rans1p2cnilowrekepsilon  evaporation_lowrekepsilon
-runSim test_md_boundary_darcy2p2cni_rans1p2cnioneeq          evaporation_oneeq
-runSim test_md_boundary_darcy2p2cni_rans1p2cnizeroeq         evaporation_zeroeq
+runSim test_md_boundary_darcy2p2cni_rans1p2cnikepsilon       evaporation_kepsilon_3       3.0
+runSim test_md_boundary_darcy2p2cni_rans1p2cnikomega         evaporation_komega_3         3.0
+runSim test_md_boundary_darcy2p2cni_rans1p2cnilowrekepsilon  evaporation_lowrekepsilon_3  3.0
+runSim test_md_boundary_darcy2p2cni_rans1p2cnioneeq          evaporation_oneeq_3          3.0
+runSim test_md_boundary_darcy2p2cni_rans1p2cnizeroeq         evaporation_zeroeq_3         3.0
+runSim test_md_boundary_darcy2p2cni_rans1p2cnikomega         evaporation_komega_0p3       0.3
+runSim test_md_boundary_darcy2p2cni_rans1p2cnikomega         evaporation_komega_15        15.0
+runSim test_md_boundary_darcy2p2cni_rans1p2cnikomega         evaporation_komega_30        30.0
+
 
 gnuplot scripts/plotevaporationrates.gp