[multidomain] added reference solution for 2cnistokes2p2cni, call for

compare
can be found in CMakeLists.txt, paramater file is now called by
ParameterFile

reviewed by gruenich



git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@12859 2fb0f335-1f38-0410-981e-8018bf24f1b0

test/multidomain/2cnistokes2p2cni/CMakeLists.txt

+# Test for 2cnistokes2p2cni can be performed by the following commands:
+# make test_2cnistokes2p2cni && ./test_2cnistokes2p2cni -ParameterFile test_references/test_2cnistokes2p2cni_reference.input | tee log_test.out
+# ../../../bin/runTest.sh ../../../bin/fuzzycomparevtu.py test_references/stokes2cni-00043.vtu stokes2cni-00043.vtu /usr/bin/echo " "
+# ../../../bin/runTest.sh ../../../bin/fuzzycomparevtu.py test_references/2p2cni-00043.vtu 2p2cni-00043.vtu /usr/bin/echo " "
 add_dumux_test(test_2cnistokes2p2cni test_2cnistokes2p2cni test_2cnistokes2p2cni
   ${CMAKE_CURRENT_BINARY_DIR}/test_2cnistokes2p2cni
   -ParameterFile ${CMAKE_CURRENT_SOURCE_DIR}/test_2cnistokes2p2cni.input)

test/multidomain/2cnistokes2p2cni/test_2cnistokes2p2cni.cc

@@ -42,7 +42,7 @@
 void printUsage(const char *progName)
 {
     std::cout << "usage: " << progName
-            << " [--restart restartTime] -parameterFile test_2cnistokes2p2cni.input\n";
+            << " [--restart restartTime] -ParameterFile test_2cnistokes2p2cni.input\n";
     exit(1);
 }
 
@@ -66,7 +66,7 @@ int start_(int argc,
     typedef typename GET_PROP(TypeTag, ParameterTree) ParameterTree;
     Dune::ParameterTreeParser::readOptions(argc, argv, ParameterTree::tree());
 
-    if (ParameterTree::tree().hasKey("parameterFile") or argc==1)
+    if (ParameterTree::tree().hasKey("ParameterFile") or argc==1)
     {
         // read input file, but do not overwrite options specified
         // on the command line, since the latter have precedence.
@@ -81,7 +81,7 @@ int start_(int argc,
             inputFileName += ".input";
         }
         else
-            inputFileName = GET_RUNTIME_PARAM(TypeTag, std::string, parameterFile); // otherwise we read from the command line
+            inputFileName = GET_RUNTIME_PARAM(TypeTag, std::string, ParameterFile); // otherwise we read from the command line
 
         std::ifstream parameterFile;
 

test/multidomain/2cnistokes2p2cni/test_references/evaporation.out

+time; massChange; evaporationRate; energyChange; massH2O
+1.26042; 0.00194216; 0.000122908; 214.81; 1.27799
+2.41688; 6.70859e-05; 5.79143e-05; 46.8515; 1.27797
+4.17568; 5.63634e-05; 5.64577e-05; 45.2422; 1.27795
+7.03372; 5.61711e-05; 5.63023e-05; 45.0176; 1.2779
+11.678; 5.60652e-05; 5.61963e-05; 44.8374; 1.27784
+18.7413; 5.59159e-05; 5.60464e-05; 44.5865; 1.27774
+29.4833; 5.57186e-05; 5.58485e-05; 44.2602; 1.27759
+45.82; 5.54695e-05; 5.55987e-05; 43.8534; 1.27736
+70.6656; 5.51615e-05; 5.52897e-05; 43.3547; 1.27702
+108.451; 5.47848e-05; 5.49118e-05; 42.7482; 1.2765
+165.917; 5.43279e-05; 5.44534e-05; 42.0147; 1.27572
+253.314; 5.37777e-05; 5.39014e-05; 41.1325; 1.27454
+386.229; 5.312e-05; 5.32415e-05; 40.0772; 1.27277
+588.371; 5.23408e-05; 5.24597e-05; 38.8234; 1.27012
+895.795; 5.14272e-05; 5.1543e-05; 37.3462; 1.26616
+1255.79; 5.05468e-05; 5.06603e-05; 36.2589; 1.2616
+1615.79; 4.93673e-05; 4.94791e-05; 35.4807; 1.25715
+1975.79; 4.69622e-05; 4.70684e-05; 33.0892; 1.25291
+2106.7; 4.53043e-05; 4.54078e-05; 31.5652; 1.25143
+2185.79; 4.43415e-05; 4.44424e-05; 30.6802; 1.25055
+2289.6; 4.29944e-05; 4.30919e-05; 29.4589; 1.24943
+2357.72; 4.19451e-05; 4.20405e-05; 28.5185; 1.24871
+2461.33; 4.02924e-05; 4.03836e-05; 27.0544; 1.24767
+2586.52; 3.80778e-05; 3.81636e-05; 25.117; 1.24647
+2737.79; 3.46506e-05; 3.47285e-05; 22.1719; 1.24516
+2858.8; 3.03629e-05; 3.04309e-05; 18.4639; 1.24424
+3005.03; 2.64974e-05; 2.65543e-05; 15.4335; 1.24327
+3227.42; 2.431e-05; 2.43599e-05; 13.7831; 1.24191
+3565.63; 2.26387e-05; 2.26848e-05; 12.665; 1.23999
+3600; 2.22244e-05; 2.22815e-05; 12.384; 1.2398

test/multidomain/2cnistokes2p2cni/test_references/test_2cnistokes2p2cni_reference.input

+#############################################################
+#Configuration file for test_2cnistokes2p2cni
+#############################################################
+
+#############################################################
+[Problem]
+#############################################################
+#First part of the interface, where it is not coupled;
+#set to zero or negative if not desired
+RunUpDistanceX = 0.0
+NoDarcyX = 0.0
+
+#############################################################
+[TimeManager]
+#############################################################
+#Initial and maximum time-step size
+DtInitial = 5e-1 
+MaxTimeStepSize = 360
+
+#Initialization time without coupling
+#set to zero or negative if not desired
+InitTime = 0 
+
+#Simulation end
+TEnd = 3600
+
+#Define the length of an episode (for output)
+EpisodeLength = 43200
+
+#############################################################
+[Vtk]
+#############################################################
+#Names for VTK output
+NameFF = stokes2cni
+NamePM = 2p2cni
+AddVelocity = 1 
+
+#############################################################
+[Grid]
+#############################################################
+UseInterfaceMeshCreator = true
+
+#Name of the dgf file (grid)
+File = grids/interfacedomain.dgf
+Refinement = 0
+
+#Number of elements in x-, y-, z-direction
+CellsX = 15
+CellsY = 45
+CellsZ = 1
+#Grading and refinement of the mesh in y direction
+Grading = 1.13
+
+#Extend of the entire domain
+XMin = 0.0
+XMax = 0.25
+YMin = 0.0
+YMax = 0.5
+
+#Vertical position of coupling interface
+InterfacePos = 0.25
+
+#############################################################
+[Output]
+#############################################################
+#Frequency of restart file, flux and VTK output
+FreqRestart = 1000      # how often restart files are written out 
+FreqOutput = 1          # frequency of VTK output
+FreqMassOutput = 1      # frequency of mass and evaporation rate output (Darcy)
+FreqFluxOutput = 5      # frequency of detailed flux output
+FreqVaporFluxOutput = 1 # frequency of summarized flux output
+
+#############################################################
+[Stokes]
+#############################################################
+StabilizationAlpha = -1.0
+
+#############################################################
+[FreeFlow]
+#############################################################
+RefPressure = 1e5
+RefTemperature = 298.15
+RefMassfrac = 0.008
+VxMax = 3.5
+BeaversJosephSlipVel = 0.00134
+ExponentMTC = 0.0           # 1./6., Mass transfer coefficient for S^MTC
+UseBoundaryLayerModel = 0   # 0 for no boundary layer model, 1 for Blasius, 2 and 3 for turbulent BL, 9 for constant thickness
+BoundaryLayerOffset = 0.0   # for BL model like Blasius, determines a virtual run-up distance for the flow
+ConstThickness = 0.0016     # for a constant BL thickness, use BL model 9
+MassTransferModel = 0       # 0 for none, 1 for power law, 2 for Schluender model
+
+SinusVelVar = 0.0   # hourly velocity variation
+SinusPVar = 0.0     # diurnal pressure variation
+SinusTVar = 0.0     # diurnal temperature variation
+SinusXVar = 0.0     # diurnal variation of the vapor concentration (massfraction)
+
+#############################################################
+[PorousMedium]
+#############################################################
+RefPressurePM = 1e5
+RefTemperaturePM = 298.15
+InitialSw1 = 0.05
+InitialSw2 = 0.05
+CharPoreRadius = 1e-4       # for Schluender mass-transfer model
+PlotFluidProperties = false
+
+#############################################################
+[SpatialParams]
+#############################################################
+MaterialInterfaceX = -100.0
+AlphaBJ = 1.0
+
+# for homogeneous setups (0 for heterogeneous):
+SoilType = 2
+RegularizationThreshold = 1e-2 # linearization threshold for pc-Sw
+
+# GStat stuff, only required for SoilType=0
+GenerateNewPermeability = true
+GStatControlFileName = gstatControl_2D.txt
+GStatInputFileName = gstatInput.txt
+PermeabilityInputFileName = permeab.dat
+
+### SoilType = 1 ### (Zurich coarse, MUSIS sand)
+[SpatialParams.Coarse]
+Permeability1 = 7.0e-10 #Kozeny-Carman; 1.387e-11 measured?
+Porosity1 = 0.44
+Swr1 = 0.005
+Snr1 = 0.01
+VgAlpha1 = 8.74e-4 
+VgN1 = 3.35
+PlotMaterialLaw1 = false
+LambdaSolid1 = 5.26
+
+Pentry1 = 1012
+BCLambda1 = 3.277
+
+### SoilType = 2 ### (Zurich fine)
+[SpatialParams.Medium]
+Permeability2 = 2.65e-10
+Porosity2 = 0.41
+Swr2 = 0.005
+Snr2 = 0.01
+VgAlpha2 = 6.371e-4
+VgN2 = 8.0 # before: 6.9
+PlotMaterialLaw2 = false
+LambdaSolid2 = 5.26
+
+Pentry2 = 1357
+BCLambda2 = 6.960
+
+### SoilType = 4 ### (Zurich fine)
+[SpatialParams.LeverettJ]
+PermeabilityJ = 1.0e-09
+PorosityJ = 0.41
+PlotMaterialLawJ = false
+
+### SoilType = 3 ### (Colorado)
+[SpatialParams.Fine]
+Permeability3 = 1.06e-10
+Porosity3 = 0.334
+Swr3 = 0.028
+Snr3 = 0.01
+VgAlpha3 = 5.81e-4
+VgN3 = 17.8
+PlotMaterialLaw3 = false
+LambdaSolid3 = 5.26
+
+Pentry3 = 1012
+BCLambda3 = 3.277
+
+#############################################################
+[Newton]
+#############################################################
+RelTolerance = 1e-5
+TargetSteps = 8
+MaxSteps = 12
+WriteConvergence = false 
+MaxTimeStepDivisions = 20
+
+#############################################################
+[LinearSolver]
+#############################################################
+ResidualReduction = 1e-9
+Verbosity = 0
+MaxIterations = 200
+#GMResRestart = 100
+#PreconditionerIterations = 2