Commit e4b41878 authored by Thomas Fetzer's avatar Thomas Fetzer
Browse files

[multidomain]

added cmake and ctest compatibility for multidomain problems

reviewed by gruenich


git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@13371 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent fa054ba9
......@@ -4,5 +4,6 @@ add_subdirectory("freeflow")
add_subdirectory("geomechanics")
add_subdirectory("implicit")
add_subdirectory("material")
add_subdirectory("multidomain")
add_directory_test_target(_test_target)
# 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
# test: test_2cnistokes2p2cni_1ff
add_dumux_test(test_2cnistokes2p2cni_1ff test_2cnistokes2p2cni test_2cnistokes2p2cni.cc
${CMAKE_SOURCE_DIR}/bin/runTest.sh
${CMAKE_SOURCE_DIR}/bin/fuzzycomparevtu.py
${CMAKE_SOURCE_DIR}/test/references/2cnistokes2p2cni-ff-reference.vtu
${CMAKE_CURRENT_BINARY_DIR}/stokes2cni-00004.vtu
${CMAKE_CURRENT_BINARY_DIR}/test_2cnistokes2p2cni
-ParameterFile ${CMAKE_CURRENT_SOURCE_DIR}/test_2cnistokes2p2cni.input)
-ParameterFile ${CMAKE_CURRENT_SOURCE_DIR}/test_2cnistokes2p2cni_reference.input
)
add_dune_superlu_flags(test_2cnistokes2p2cni)
add_dune_ug_flags(test_2cnistokes2p2cni)
# test: test_2cnistokes2p2cni_2pm
add_dumux_test(test_2cnistokes2p2cni_2pm test_2cnistokes2p2cni test_2cnistokes2p2cni.cc
${CMAKE_SOURCE_DIR}/bin/runTest.sh
${CMAKE_SOURCE_DIR}/bin/fuzzycomparevtu.py
${CMAKE_SOURCE_DIR}/test/references/2cnistokes2p2cni-pm-reference.vtu
${CMAKE_CURRENT_BINARY_DIR}/2p2cni-00004.vtu
echo "test_2cnistokes2p2cni_pm"
)
\ No newline at end of file
#############################################################
#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
UseDirichletFF = true
#############################################################
[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 = 20
CellsZ = 1
#Grading and refinement of the mesh in y direction
Grading = 1.1
#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 = 10 # frequency of VTK output
FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy)
FreqFluxOutput = 1000 # frequency of detailed flux output
FreqVaporFluxOutput = 2 # frequency of summarized flux output
#############################################################
[Stokes]
#############################################################
StabilizationAlpha = -1.0
#############################################################
[FreeFlow]
#############################################################
RefVelocity = 3.5
RefPressure = 1e5
RefMassfrac = 0.008 # 0.01
RefTemperature = 298.15
SinusVelAmplitude = 0.0 # 0.1
SinusVelPeriod = 3600 # in seconds
SinusPressureAmplitude = 0.0 # 0.5
SinusPressurePeriod = 3600 # in seconds
SinusConcentrationAmplitude = 0.0 # 0.5
SinusConcentrationPeriod = 3600 # in seconds
SinusTemperatureAmplitude = 0.0 # 0.5
SinusTemperaturePeriod = 3600 # in seconds
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
#############################################################
[PorousMedium]
#############################################################
RefPressurePM = 1e5
RefTemperaturePM = 298.15
InitialSw1 = 0.98
InitialSw2 = 0.98
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
# Test for 2cstokes2p2c can be performed by the following commands:
# make test_2cstokes2p2c && ./test_2cstokes2p2c -ParameterFile test_references/test_2cstokes2p2c_reference.input | tee log_test.out
# ../../../bin/runTest.sh ../../../bin/fuzzycomparevtu.py test_references/stokes2c-00043.vtu stokes2c-00043.vtu /usr/bin/echo " "
# ../../../bin/runTest.sh ../../../bin/fuzzycomparevtu.py test_references/2p2c-00043.vtu 2p2c-00043.vtu /usr/bin/echo " "
add_dumux_test(test_2cstokes2p2c test_2cstokes2p2c test_2cstokes2p2c
# test: test_2cstokes2p2c_1ff
add_dumux_test(test_2cstokes2p2c_1ff test_2cstokes2p2c test_2cstokes2p2c.cc
${CMAKE_SOURCE_DIR}/bin/runTest.sh
${CMAKE_SOURCE_DIR}/bin/fuzzycomparevtu.py
${CMAKE_SOURCE_DIR}/test/references/2cstokes2p2c-ff-reference.vtu
${CMAKE_CURRENT_BINARY_DIR}/stokes2c-00007.vtu
${CMAKE_CURRENT_BINARY_DIR}/test_2cstokes2p2c
-ParameterFile ${CMAKE_CURRENT_SOURCE_DIR}/test_2cstokes2p2c.input)
-ParameterFile ${CMAKE_CURRENT_SOURCE_DIR}/test_2cstokes2p2c_reference.input
)
add_dune_superlu_flags(test_2cstokes2p2c)
add_dune_ug_flags(test_2cstokes2p2c)
# test: test_2cstokes2p2c_2pm
add_dumux_test(test_2cstokes2p2c_2pm test_2cstokes2p2c test_2cstokes2p2c.cc
${CMAKE_SOURCE_DIR}/bin/runTest.sh
${CMAKE_SOURCE_DIR}/bin/fuzzycomparevtu.py
${CMAKE_SOURCE_DIR}/test/references/2cstokes2p2c-pm-reference.vtu
${CMAKE_CURRENT_BINARY_DIR}/2p2c-00007.vtu
echo "test_2cstokes2p2c_pm"
)
\ No newline at end of file
#############################################################
#Configuration file for test_2cstokes2p2c
#############################################################
#############################################################
[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= 7200 # 2 hours
#Define the length of an episode (for output)
EpisodeLength = 3600
#############################################################
[Vtk]
#############################################################
#Names for VTK output
NameFF = stokes2c
NamePM = 2p2c
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 = 25
CellsY = 45
CellsZ = 1
#Grading and refinement of the mesh in y direction
Grading = 1.1
#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 = 10 # frequency of VTK output
FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy)
FreqFluxOutput = 1000 # frequency of detailed flux output
FreqVaporFluxOutput = 2 # frequency of summarized flux output
#############################################################
[Stokes]
#############################################################
StabilizationAlpha = -1.0
#############################################################
[FreeFlow]
#############################################################
RefVelocity = 3.5
RefPressure = 1e5
RefMassfrac = 0.008 # 0.01
RefTemperature = 298.15
SinusVelAmplitude = 0.0 # 0.1
SinusVelPeriod = 3600 # in seconds
SinusPressureAmplitude = 0.0 # 0.5
SinusPressurePeriod = 3600 # in seconds
SinusConcentrationAmplitude = 0.0 # 0.5
SinusConcentrationPeriod = 3600 # in seconds
SinusTemperatureAmplitude = 0.0 # 0.5
SinusTemperaturePeriod = 3600 # in seconds
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
#############################################################
[PorousMedium]
#############################################################
RefPressurePM = 1e5
RefTemperaturePM = 298.15
InitialSw1 = 0.98
InitialSw2 = 0.98
CharPoreDiameter = 1e-4 # for Schluender mass-transfer model
#############################################################
[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 ###
[SpatialParams.Coarse]
Permeability1 = 7.0e-10
Porosity1 = 0.44
Swr1 = 0.005
Snr1 = 0.01
VgAlpha1 = 8.74e-4
VgN1 = 3.35
PlotMaterialLaw1 = false
LambdaSolid1 = 5.3
### SoilType = 2 ###
[SpatialParams.Medium]
Permeability2 = 2.65e-10
Porosity2 = 0.41
Swr2 = 0.005
Snr2 = 0.01
VgAlpha2 = 6.371e-4
VgN2 = 6.9
PlotMaterialLaw2 = false
LambdaSolid2 = 5.3
### SoilType = 3 ###
[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.3
#############################################################
[Newton]
#############################################################
RelTolerance = 1e-5
TargetSteps = 8
MaxSteps = 12
WriteConvergence = false
MaxTimeStepDivisions = 20
#############################################################
[LinearSolver]
#############################################################
ResidualReduction = 1e-10
Verbosity = 0
MaxIterations = 200
add_subdirectory(2cstokes2p2c)
add_subdirectory(2cnistokes2p2cni)
add_subdirectory("2cstokes2p2c")
add_subdirectory("2cnistokes2p2cni")
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
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