Skip to content
Snippets Groups Projects
Commit f7aa645e authored by Thomas Fetzer's avatar Thomas Fetzer
Browse files

[lineendings] Restore line ending type to Unix for all files

parent 3aa0edb7
No related branches found
No related tags found
2 merge requests!31Feature/colebrookwhiteboundarylayer,!12Feature/externalscript
Hide whitespace changes
Inline Side-by-side
Showing
with 895 additions and 895 deletions
test/decoupled/2p2c/test_adaptive2p2c2d.input
+ 29
− 29
View file @ f7aa645e
[Impet] [Impet]
CFLFactor = 0.8 CFLFactor = 0.8
ErrorTermLowerBound = 0.2 ErrorTermLowerBound = 0.2
ErrorTermUpperBound = 0.9 ErrorTermUpperBound = 0.9
ErrorTermFactor = 0.5 ErrorTermFactor = 0.5
EnableVolumeIntegral = 1 EnableVolumeIntegral = 1
RestrictFluxInTransport = 1 RestrictFluxInTransport = 1
[GridAdapt] [GridAdapt]
MinLevel = 0 MinLevel = 0
MaxLevel = 1 MaxLevel = 1
EnableMultiPointFluxApproximation = 1 EnableMultiPointFluxApproximation = 1
MaxInteractionVolumes = 2 MaxInteractionVolumes = 2
[Problem] [Problem]
SimulationName = test_adaptive2p2c2d # name of the output files SimulationName = test_adaptive2p2c2d # name of the output files
OutputInterval = 5 OutputInterval = 5
[TimeManager] [TimeManager]
TEnd= 3e3 TEnd= 3e3
DtInitial = 200 DtInitial = 200
[Grid] [Grid]
NumberOfCellsX = 10# [-] level 0 resolution in x-direction NumberOfCellsX = 10# [-] level 0 resolution in x-direction
NumberOfCellsY = 10# [-] level 0 resolution in y-direction NumberOfCellsY = 10# [-] level 0 resolution in y-direction
UpperRightX = 10# [m] dimension of the grid UpperRightX = 10# [m] dimension of the grid
UpperRightY = 10# [m] dimension of the grid # 2D UpperRightY = 10# [m] dimension of the grid # 2D
test/decoupled/2p2c/test_adaptive2p2c3d.input
+ 33
− 33
View file @ f7aa645e
[Impet] [Impet]
CFLFactor = 0.8 CFLFactor = 0.8
ErrorTermLowerBound = 0.02 ErrorTermLowerBound = 0.02
ErrorTermUpperBound = 0.9 ErrorTermUpperBound = 0.9
ErrorTermFactor = 0.5 ErrorTermFactor = 0.5
EnableVolumeIntegral = 1 EnableVolumeIntegral = 1
RestrictFluxInTransport = 1 RestrictFluxInTransport = 1
[GridAdapt] [GridAdapt]
EnableMultiPointFluxApproximation = 1 EnableMultiPointFluxApproximation = 1
MaxInteractionVolumes = 4 MaxInteractionVolumes = 4
MinLevel = 0 MinLevel = 0
MaxLevel = 1 MaxLevel = 1
[Grid] [Grid]
NumberOfCellsX = 10 # [-] resolution in x-direction NumberOfCellsX = 10 # [-] resolution in x-direction
NumberOfCellsY = 10 # [-] resolution in y-direction NumberOfCellsY = 10 # [-] resolution in y-direction
NumberOfCellsZ = 10 # [-] resolution in z-direction NumberOfCellsZ = 10 # [-] resolution in z-direction
UpperRightX = 10 # [m] length of the domain UpperRightX = 10 # [m] length of the domain
UpperRightY = 10 # [m] width of the domain UpperRightY = 10 # [m] width of the domain
UpperRightZ = 10 # [m] height of the domain UpperRightZ = 10 # [m] height of the domain
[BoundaryConditions] [BoundaryConditions]
Injectionrate = 1e-4 # kg / sec Injectionrate = 1e-4 # kg / sec
[Problem] [Problem]
Name = test_adaptive2p2c3d # name of the output files Name = test_adaptive2p2c3d # name of the output files
OutputInterval = 5 OutputInterval = 5
[TimeManager] [TimeManager]
TEnd= 3e3 #Simulated time TEnd= 3e3 #Simulated time
DtInitial = 500 #Initial time step DtInitial = 500 #Initial time step
test/implicit/2pnc/test_box2pnc.input
+ 56
− 56
View file @ f7aa645e
# Parameter file for test case 2pncmin. # Parameter file for test case 2pncmin.
# Everything behind a '#' is a comment. # Everything behind a '#' is a comment.
# Type "./test_2pncmin --help" for more information. # Type "./test_2pncmin --help" for more information.
[TimeManager] [TimeManager]
DtInitial = 5e-1 # initial time step size [s] DtInitial = 5e-1 # initial time step size [s]
MaxTimeStepSize = 1000 # maximum time step size MaxTimeStepSize = 1000 # maximum time step size
TEnd= 1e3 # duration of the simulation [s] TEnd= 1e3 # duration of the simulation [s]
[Grid] [Grid]
#Number of cells in X and Y directions #Number of cells in X and Y directions
NumberOfCellsX = 21 NumberOfCellsX = 21
NumberOfCellsY = 6 NumberOfCellsY = 6
#Extend of entire domain #Extend of entire domain
UpperRightX = 2.0e-3 # x-coordinate of the upper-right corner of the grid [m] UpperRightX = 2.0e-3 # x-coordinate of the upper-right corner of the grid [m]
UpperRightY = 5.5e-4 # y-coordinate of the upper-right corner of the grid [m] UpperRightY = 5.5e-4 # y-coordinate of the upper-right corner of the grid [m]
[Problem] [Problem]
Name = fuelcell Name = fuelcell
EnableGravity = 0 EnableGravity = 0
[FluidSystem] [FluidSystem]
NTemperature = 3 # [-] number of temperature table entries NTemperature = 3 # [-] number of temperature table entries
NPressure = 200 # [-] number of pressure table entries NPressure = 200 # [-] number of pressure table entries
PressureLow = 1e5 # [Pa] lower pressure limit for tabularization PressureLow = 1e5 # [Pa] lower pressure limit for tabularization
PressureHigh = 3e7 # [Pa] upper pressure limit for tabularization PressureHigh = 3e7 # [Pa] upper pressure limit for tabularization
TemperatureLow = 312.15 # [K] lower temperature limit for tabularization TemperatureLow = 312.15 # [K] lower temperature limit for tabularization
TemperatureHigh = 314.15 # [K] upper temperature limit for tabularization TemperatureHigh = 314.15 # [K] upper temperature limit for tabularization
InitialTemperature = 310 # [K] initial temperature for tabularization InitialTemperature = 310 # [K] initial temperature for tabularization
[ElectroChemistry] [ElectroChemistry]
SpecificResistance = 0.25 #[Ohm*cm^2] SpecificResistance = 0.25 #[Ohm*cm^2]
ReversibleVoltage = 1.191 #[V] ReversibleVoltage = 1.191 #[V]
CellVoltage = 0.5 #[V] CellVoltage = 0.5 #[V]
ThermoneutralVoltage = 1.4836 #[V] ThermoneutralVoltage = 1.4836 #[V]
RefCurrentDensity = 1.87e-8 #[A/cm] for calculating exchange current density at reference conditions RefCurrentDensity = 1.87e-8 #[A/cm] for calculating exchange current density at reference conditions
RefO2PartialPressure = 5.0e5 #[Pa] for calculating exchange current density at reference conditions RefO2PartialPressure = 5.0e5 #[Pa] for calculating exchange current density at reference conditions
RefTemperature = 353.25 #[K] for calculating exchange current density at reference conditions RefTemperature = 353.25 #[K] for calculating exchange current density at reference conditions
ActivationBarrier = 73.0e3 #[J/mol] ActivationBarrier = 73.0e3 #[J/mol]
TransferCoefficient = 0.5 #[-] (alpha) TransferCoefficient = 0.5 #[-] (alpha)
NumElectrons = 2 #[-] number of electrons transferred in reaction NumElectrons = 2 #[-] number of electrons transferred in reaction
SurfaceIncreasingFactor = 60 SurfaceIncreasingFactor = 60
TransportNumberH20 = 0.2327 #[-] account for osmotic H20 transport term from membrane, value Lena Walter TransportNumberH20 = 0.2327 #[-] account for osmotic H20 transport term from membrane, value Lena Walter
pO2Inlet = 0.3e5 #[Pa] partial pressure of oxygen at gas channel inlet pO2Inlet = 0.3e5 #[Pa] partial pressure of oxygen at gas channel inlet
MaxIterations = 300 #[-] Maximum number for iterations for solving electrochemical system MaxIterations = 300 #[-] Maximum number for iterations for solving electrochemical system
[Newton] [Newton]
RelTolerance = 1e-8 RelTolerance = 1e-8
TargetSteps = 10 TargetSteps = 10
MaxSteps = 15 MaxSteps = 15
WriteConvergence = false WriteConvergence = false
MaxTimeStepDivisions = 20 MaxTimeStepDivisions = 20
[LinearSolver] [LinearSolver]
ResidualReduction = 1e-8 ResidualReduction = 1e-8
Verbosity = 0 Verbosity = 0
MaxIterations = 100 MaxIterations = 100
\ No newline at end of file
test/implicit/mpnc/test_boxmpnckinetic.input
+ 75
− 75
View file @ f7aa645e
[TimeManager] [TimeManager]
DtInitial = 1.5 # [s] DtInitial = 1.5 # [s]
TEnd = 10 # [s] TEnd = 10 # [s]
MaxTimeStepSize = 2e20 # maximum allowed timestep size MaxTimeStepSize = 2e20 # maximum allowed timestep size
[Grid] [Grid]
NumberOfCellsX = 14 # NumberOfCellsX = 14 #
NumberOfCellsY = 30 # NumberOfCellsY = 30 #
LowerLeftX = 0.0 LowerLeftX = 0.0
UpperRightX = 1.0 UpperRightX = 1.0
LowerLeftY = 0.0 LowerLeftY = 0.0
UpperRightY = 0.5 UpperRightY = 0.5
InterfacePosY = 0.25 InterfacePosY = 0.25
GradingFactorY = 1.13 GradingFactorY = 1.13
RefineTop = false RefineTop = false
[InitialConditions] [InitialConditions]
SwFFInitial = 1e-5 # - 0.001# SwFFInitial = 1e-5 # - 0.001#
SwPMInitial = 0.999 # - SwPMInitial = 0.999 # -
pnInitial = 1e5 # 6.8e6 # 1e5 # 101475 # Pa pnInitial = 1e5 # 6.8e6 # 1e5 # 101475 # Pa
TInitial = 293 TInitial = 293
pnInjection = 100003 pnInjection = 100003
[BoundaryConditions] [BoundaryConditions]
percentOfEquil = .1 # percentage of the equilibrium concentration water in the nPhase influx percentOfEquil = .1 # percentage of the equilibrium concentration water in the nPhase influx
TInject = 293 TInject = 293
massFluxInjectedPhase = 0.75 massFluxInjectedPhase = 0.75
[SpatialParams.PorousMedium] [SpatialParams.PorousMedium]
permeability = 1e-11 # m^2 permeability = 1e-11 # m^2
porosity = 0.4 # porosity = 0.4 #
meanPoreSize = 1e-4 # characteristic length of the system meanPoreSize = 1e-4 # characteristic length of the system
factorEnergyTransfer = 1 # factorEnergyTransfer = 1 #
factorMassTransfer = 1 # factorMassTransfer = 1 #
[SpatialParams.FreeFlow] [SpatialParams.FreeFlow]
permeability = 1e-6 # m^2 permeability = 1e-6 # m^2
porosity = 0.99 # porosity = 0.99 #
meanPoreSize = 1e-2 # 4e-5 # characteristic length of the system meanPoreSize = 1e-2 # 4e-5 # characteristic length of the system
[SpatialParams.soil] [SpatialParams.soil]
density = 2600. # kg/m^3 http://www.agriinfo.in/default.aspx?page=topic&superid=4&topicid=271 density = 2600. # kg/m^3 http://www.agriinfo.in/default.aspx?page=topic&superid=4&topicid=271
thermalConductivity = 3 # W / (m K) thermalConductivity = 3 # W / (m K)
heatCapacity= 817 # J / (kg K) heatCapacity= 817 # J / (kg K)
# characteristic length, i.e. mean pore size # characteristic length, i.e. mean pore size
# 40 micrometer i.e. 1e-5 for the micromodel # 40 micrometer i.e. 1e-5 for the micromodel
VGAlpha = 3.512e-04 # VGAlpha = 3.512e-04 #
VGN = 4.716e+00 # VGN = 4.716e+00 #
Swr = 0 # Swr = 0 #
Snr = 0 # Snr = 0 #
specificSolidsurface = 4022.994# (1/m) specificSolidsurface = 4022.994# (1/m)
aWettingNonWettingA1 = -1.603e-01 # aWettingNonWettingA1 = -1.603e-01 #
aWettingNonWettingA2 = 1.429e-05 # aWettingNonWettingA2 = 1.429e-05 #
aWettingNonWettingA3 = 1.915e-01 # aWettingNonWettingA3 = 1.915e-01 #
BCPd = 2.290e+03 # BCPd = 2.290e+03 #
BClambda = 2.740e+00 # BClambda = 2.740e+00 #
aNonWettingSolidA1 = 1.369e+03 # aNonWettingSolidA1 = 1.369e+03 #
aNonWettingSolidA2 = -3.782e+00 # aNonWettingSolidA2 = -3.782e+00 #
aNonWettingSolidA3 = 1.063e-09 # aNonWettingSolidA3 = 1.063e-09 #
[SourceSink] [SourceSink]
heatIntoSolid = 0 # heatIntoSolid = 0 #
[Constants] [Constants]
outputName = evaporationatmosphere outputName = evaporationatmosphere
nRestart = 100 # after so many timesteps a restart file should be written nRestart = 100 # after so many timesteps a restart file should be written
[FluidSystem] [FluidSystem]
nTemperature = 100 # for the tabulation of water: so many interpolation points nTemperature = 100 # for the tabulation of water: so many interpolation points
nPressure = 100 # for the tabulation of water: so many interpolation points nPressure = 100 # for the tabulation of water: so many interpolation points
hammer = 1e4 hammer = 1e4
test/implicit/mpnc/test_boxmpncthermalnonequil.input
+ 55
− 55
View file @ f7aa645e
[TimeManager] [TimeManager]
DtInitial = 5e-1 # [s] DtInitial = 5e-1 # [s]
TEnd = 1e3 # [s] TEnd = 1e3 # [s]
MaxTimeStepSize = 2e20 # maximum allowed timestep size MaxTimeStepSize = 2e20 # maximum allowed timestep size
[Grid] [Grid]
File = ./grids/combustionOutflowGridLinNX100LogNx100.dgf File = ./grids/combustionOutflowGridLinNX100LogNx100.dgf
# File = ./grids/combustionGrid.dgf # File = ./grids/combustionGrid.dgf
lengthPM = 0.1 lengthPM = 0.1
[InitialConditions] [InitialConditions]
pnInitial = 1e5 # 6.8e6 # 1e5 # 101475 # Pa pnInitial = 1e5 # 6.8e6 # 1e5 # 101475 # Pa
TInitial = 372 # 373.15 # TInitial = 372 # 373.15 #
TRight = 400 # 373.15 # TRight = 400 # 373.15 #
[BoundaryConditions] [BoundaryConditions]
TBoundary =300 # TBoundary =300 #
SwBoundary = 1e-3 # SwBoundary = 1e-3 #
SwOneComponentSys = 1 # SwOneComponentSys = 1 #
massFluxInjectedPhase = 0.5 # massFluxInjectedPhase = 0.5 #
# heatFluxFromRight = 2e6 # fluid temperature on boundary: analytical approx 1080 K # heatFluxFromRight = 2e6 # fluid temperature on boundary: analytical approx 1080 K
heatFluxFromRight = 1.5e6 # fluid temperature on boundary: analytical approx 587.3 K heatFluxFromRight = 1.5e6 # fluid temperature on boundary: analytical approx 587.3 K
coldTime = 1 coldTime = 1
[SpatialParams.PorousMedium] [SpatialParams.PorousMedium]
porosity = 0.35 # porosity = 0.35 #
meanPoreSize = 5e-4 # characteristic length of the system meanPoreSize = 5e-4 # characteristic length of the system
factorEnergyTransfer = 1 # factorEnergyTransfer = 1 #
factorMassTransfer = 1 # factorMassTransfer = 1 #
[SpatialParams.OutFlow] [SpatialParams.OutFlow]
permeabilityOutFlow = 1e-6 # m^2 permeabilityOutFlow = 1e-6 # m^2
porosityOutFlow = 0.35 # porosityOutFlow = 0.35 #
soilThermalConductivityOutFlow = 0.01# soilThermalConductivityOutFlow = 0.01#
[SpatialParams.soil] [SpatialParams.soil]
density = 2600. # kg/m^3 density = 2600. # kg/m^3
thermalConductivity = 30 # W / (m K) soil:3, metal:30 thermalConductivity = 30 # W / (m K) soil:3, metal:30
heatCapacity = 466 # J / (kg K) steel:466 granite:800 heatCapacity = 466 # J / (kg K) steel:466 granite:800
Swr = 0.0 # 5e-3 Swr = 0.0 # 5e-3
Snr = 0 # Snr = 0 #
[Constants] [Constants]
outputName = combustion # outputName = combustion #
interfacialTension = 0.0589 # interfacial tension of water at 100 ° C interfacialTension = 0.0589 # interfacial tension of water at 100 ° C
nRestart = 10000 # after so many timesteps a restart file should be written nRestart = 10000 # after so many timesteps a restart file should be written
[Implicit] [Implicit]
MassUpwindWeight=1 MassUpwindWeight=1
MobilityUpwindWeight=1 MobilityUpwindWeight=1
test/multidomain/2cnistokes2p2cni/test_2cnistokes2p2cni.input
+ 82
− 82
View file @ f7aa645e
[TimeManager] [TimeManager]
DtInitial = 5e-1 # [s] DtInitial = 5e-1 # [s]
MaxTimeStepSize = 360 # [s] MaxTimeStepSize = 360 # [s]
InitTime = 864 # [s] Initialization time without coupling InitTime = 864 # [s] Initialization time without coupling
TEnd= 0.864e6 # [s] TEnd= 0.864e6 # [s]
EpisodeLength = 43200 # [s] EpisodeLength = 43200 # [s]
[Grid] [Grid]
# File = grids/test_2cnistokes2p2cni.dgf # unused, if interface mesh creator is used # File = grids/test_2cnistokes2p2cni.dgf # unused, if interface mesh creator is used
LowerLeftX = 0.0 # [m] LowerLeftX = 0.0 # [m]
UpperRightX = 0.25 # [m] UpperRightX = 0.25 # [m]
LowerLeftY = 0.0 # [m] LowerLeftY = 0.0 # [m]
UpperRightY = 0.5 # [m] UpperRightY = 0.5 # [m]
RunUpDistanceX = 0.0 # [m] Horizontal position without coupling to PM RunUpDistanceX = 0.0 # [m] Horizontal position without coupling to PM
NoDarcyX = 0.0 # [m] Horizontal position without PM below NoDarcyX = 0.0 # [m] Horizontal position without PM below
InterfacePosY = 0.25 # [m] Vertical position of coupling interface InterfacePosY = 0.25 # [m] Vertical position of coupling interface
# Number of elements in x-, y-, z-direction # Number of elements in x-, y-, z-direction
NumberOfCellsX = 30 NumberOfCellsX = 30
NumberOfCellsY = 62 NumberOfCellsY = 62
# Grading and refinement of the mesh in y direction # Grading and refinement of the mesh in y direction
GradingFactorY = 1.13 GradingFactorY = 1.13
RefineTopY = false RefineTopY = false
[Output] [Output]
NameFF = stokes2cni NameFF = stokes2cni
NamePM = 2p2cni NamePM = 2p2cni
#Frequency of restart file, flux and VTK output #Frequency of restart file, flux and VTK output
FreqRestart = 1000 # how often restart files are written out FreqRestart = 1000 # how often restart files are written out
FreqOutput = 10 # frequency of VTK output FreqOutput = 10 # frequency of VTK output
FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy) FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy)
FreqFluxOutput = 1000 # frequency of detailed flux output FreqFluxOutput = 1000 # frequency of detailed flux output
FreqVaporFluxOutput = 2 # frequency of summarized flux output FreqVaporFluxOutput = 2 # frequency of summarized flux output
LiveEvaporationRates = true # plot evaporation rates using gnuplot interface LiveEvaporationRates = true # plot evaporation rates using gnuplot interface
[Stokes] [Stokes]
StabilizationAlpha = -1.0 StabilizationAlpha = -1.0
[FreeFlow] [FreeFlow]
UseDirichletBC = true # dirichlet values are set at the inflow boundary UseDirichletBC = true # dirichlet values are set at the inflow boundary
RefVelocity = 3.5 # [m/s] RefVelocity = 3.5 # [m/s]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefMassfrac = 0.008 # [-] RefMassfrac = 0.008 # [-]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
SinusVelAmplitude = 0.0 # [m/s] SinusVelAmplitude = 0.0 # [m/s]
SinusVelPeriod = 3600 # [s] SinusVelPeriod = 3600 # [s]
SinusPressureAmplitude = 0.0 # [Pa] SinusPressureAmplitude = 0.0 # [Pa]
SinusPressurePeriod = 3600 # [s] SinusPressurePeriod = 3600 # [s]
SinusConcentrationAmplitude = 0.0 # [-] SinusConcentrationAmplitude = 0.0 # [-]
SinusConcentrationPeriod = 3600 # [s] SinusConcentrationPeriod = 3600 # [s]
SinusTemperatureAmplitude = 0.0 # [K] SinusTemperatureAmplitude = 0.0 # [K]
SinusTemperaturePeriod = 3600 # [s] SinusTemperaturePeriod = 3600 # [s]
[BoundaryLayer] [BoundaryLayer]
Model = 0 Model = 0
[MassTransfer] [MassTransfer]
Model = 0 Model = 0
[PorousMedium] [PorousMedium]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
InitialSw = 0.98 # [-] InitialSw = 0.98 # [-]
[SpatialParams] [SpatialParams]
AlphaBJ = 1.0 # [-] AlphaBJ = 1.0 # [-]
Permeability = 2.65e-10 # [m^2] Permeability = 2.65e-10 # [m^2]
Porosity = 0.41 # [-] Porosity = 0.41 # [-]
Swr = 0.005 # [-] Swr = 0.005 # [-]
Snr = 0.01 # [-] Snr = 0.01 # [-]
VgAlpha = 6.371e-4 # [1/Pa] VgAlpha = 6.371e-4 # [1/Pa]
VgN = 8.0 # [-] VgN = 8.0 # [-]
LambdaSolid = 5.26 # [W/(m*K)] LambdaSolid = 5.26 # [W/(m*K)]
[Newton] [Newton]
MaxRelativeShift = 1e-5 MaxRelativeShift = 1e-5
TargetSteps = 8 TargetSteps = 8
MaxSteps = 12 MaxSteps = 12
WriteConvergence = false WriteConvergence = false
MaxTimeStepDivisions = 20 MaxTimeStepDivisions = 20
[LinearSolver] [LinearSolver]
Verbosity = 0 Verbosity = 0
test/multidomain/2cnistokes2p2cni/test_2cnistokes2p2cni_boundarylayer.input
+ 85
− 85
View file @ f7aa645e
[TimeManager] [TimeManager]
DtInitial = 5e-1 # [s] DtInitial = 5e-1 # [s]
MaxTimeStepSize = 360 # [s] MaxTimeStepSize = 360 # [s]
InitTime = 0 # [s] Initialization time without coupling InitTime = 0 # [s] Initialization time without coupling
TEnd= 3600 # [s] TEnd= 3600 # [s]
EpisodeLength = 43200 # [s] EpisodeLength = 43200 # [s]
[Grid] [Grid]
# File = grids/test_2cnistokes2p2cni.dgf # unused, if interface mesh creator is used # File = grids/test_2cnistokes2p2cni.dgf # unused, if interface mesh creator is used
LowerLeftX = 0.0 # [m] LowerLeftX = 0.0 # [m]
UpperRightX = 0.25 # [m] UpperRightX = 0.25 # [m]
LowerLeftY = 0.0 # [m] LowerLeftY = 0.0 # [m]
UpperRightY = 0.5 # [m] UpperRightY = 0.5 # [m]
RunUpDistanceX = 0.0 # [m] Horizontal position without coupling to PM RunUpDistanceX = 0.0 # [m] Horizontal position without coupling to PM
NoDarcyX = 0.0 # [m] Horizontal position without PM below NoDarcyX = 0.0 # [m] Horizontal position without PM below
InterfacePosY = 0.25 # [m] Vertical position of coupling interface InterfacePosY = 0.25 # [m] Vertical position of coupling interface
# Number of elements in x-, y-, z-direction # Number of elements in x-, y-, z-direction
NumberOfCellsX = 15 NumberOfCellsX = 15
NumberOfCellsY = 20 NumberOfCellsY = 20
# Grading and refinement of the mesh in y direction # Grading and refinement of the mesh in y direction
GradingFactorY = 1.1 GradingFactorY = 1.1
RefineTopY = false RefineTopY = false
[Output] [Output]
NameFF = stokes2cni_boundarylayer NameFF = stokes2cni_boundarylayer
NamePM = 2p2cni_boundarylayer NamePM = 2p2cni_boundarylayer
#Frequency of restart file, flux and VTK output #Frequency of restart file, flux and VTK output
FreqRestart = 1000 # how often restart files are written out FreqRestart = 1000 # how often restart files are written out
FreqOutput = 10 # frequency of VTK output FreqOutput = 10 # frequency of VTK output
FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy) FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy)
FreqFluxOutput = 1000 # frequency of detailed flux output FreqFluxOutput = 1000 # frequency of detailed flux output
FreqVaporFluxOutput = 2 # frequency of summarized flux output FreqVaporFluxOutput = 2 # frequency of summarized flux output
LiveEvaporationRates = false # plot evaporation rates using gnuplot interface LiveEvaporationRates = false # plot evaporation rates using gnuplot interface
[Stokes] [Stokes]
StabilizationAlpha = -1.0 StabilizationAlpha = -1.0
[FreeFlow] [FreeFlow]
UseDirichletBC = true # dirichlet values are set at the inflow boundary UseDirichletBC = true # dirichlet values are set at the inflow boundary
RefVelocity = 3.5 # [m/s] RefVelocity = 3.5 # [m/s]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefMassfrac = 0.008 # [-] RefMassfrac = 0.008 # [-]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
SinusVelAmplitude = 0.0 # [m/s] SinusVelAmplitude = 0.0 # [m/s]
SinusVelPeriod = 3600 # [s] SinusVelPeriod = 3600 # [s]
SinusPressureAmplitude = 0.0 # [Pa] SinusPressureAmplitude = 0.0 # [Pa]
SinusPressurePeriod = 3600 # [s] SinusPressurePeriod = 3600 # [s]
SinusConcentrationAmplitude = 0.0 # [-] SinusConcentrationAmplitude = 0.0 # [-]
SinusConcentrationPeriod = 3600 # [s] SinusConcentrationPeriod = 3600 # [s]
SinusTemperatureAmplitude = 0.0 # [K] SinusTemperatureAmplitude = 0.0 # [K]
SinusTemperaturePeriod = 3600 # [s] SinusTemperaturePeriod = 3600 # [s]
[BoundaryLayer] [BoundaryLayer]
Model = 4 # Number/ID of the used model Model = 4 # Number/ID of the used model
Offset = 0.25 # Virtual run-up distance for BL models [m] Offset = 0.25 # Virtual run-up distance for BL models [m]
YPlus = 10 # Conversion value (model 4-6) [-] YPlus = 10 # Conversion value (model 4-6) [-]
[MassTransfer] [MassTransfer]
Model = 4 Model = 4
CharPoreRadius = 0.001 # Characteristic pore radius for Schluender model (model 2+4) [m] CharPoreRadius = 0.001 # Characteristic pore radius for Schluender model (model 2+4) [m]
[PorousMedium] [PorousMedium]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
InitialSw = 0.98 # [-] InitialSw = 0.98 # [-]
[SpatialParams] [SpatialParams]
AlphaBJ = 1.0 # [-] AlphaBJ = 1.0 # [-]
Permeability = 2.65e-10 # [m^2] Permeability = 2.65e-10 # [m^2]
Porosity = 0.41 # [-] Porosity = 0.41 # [-]
Swr = 0.005 # [-] Swr = 0.005 # [-]
Snr = 0.01 # [-] Snr = 0.01 # [-]
VgAlpha = 6.371e-4 # [1/Pa] VgAlpha = 6.371e-4 # [1/Pa]
VgN = 8.0 # [-] VgN = 8.0 # [-]
LambdaSolid = 5.26 # [W/(m*K)] LambdaSolid = 5.26 # [W/(m*K)]
[Newton] [Newton]
MaxRelativeShift = 1e-5 MaxRelativeShift = 1e-5
TargetSteps = 8 TargetSteps = 8
MaxSteps = 12 MaxSteps = 12
WriteConvergence = false WriteConvergence = false
MaxTimeStepDivisions = 20 MaxTimeStepDivisions = 20
[LinearSolver] [LinearSolver]
Verbosity = 0 Verbosity = 0
test/multidomain/2cnistokes2p2cni/test_2cnistokes2p2cni_reference.input
+ 82
− 82
View file @ f7aa645e
[TimeManager] [TimeManager]
DtInitial = 5e-1 # [s] DtInitial = 5e-1 # [s]
MaxTimeStepSize = 240 # [s] MaxTimeStepSize = 240 # [s]
InitTime = 0 # [s] Initialization time without coupling InitTime = 0 # [s] Initialization time without coupling
TEnd= 3600 # [s] TEnd= 3600 # [s]
EpisodeLength = 43200 # [s] EpisodeLength = 43200 # [s]
[Grid] [Grid]
# File = grids/test_2cnistokes2p2cni.dgf # unused, if interface mesh creator is used # File = grids/test_2cnistokes2p2cni.dgf # unused, if interface mesh creator is used
LowerLeftX = 0.0 # [m] LowerLeftX = 0.0 # [m]
UpperRightX = 0.25 # [m] UpperRightX = 0.25 # [m]
LowerLeftY = 0.0 # [m] LowerLeftY = 0.0 # [m]
UpperRightY = 0.5 # [m] UpperRightY = 0.5 # [m]
RunUpDistanceX = 0.0 # [m] Horizontal position without coupling to PM RunUpDistanceX = 0.0 # [m] Horizontal position without coupling to PM
NoDarcyX = 0.0 # [m] Horizontal position without PM below NoDarcyX = 0.0 # [m] Horizontal position without PM below
InterfacePosY = 0.25 # [m] Vertical position of coupling interface InterfacePosY = 0.25 # [m] Vertical position of coupling interface
# Number of elements in x-, y-, z-direction # Number of elements in x-, y-, z-direction
NumberOfCellsX = 15 NumberOfCellsX = 15
NumberOfCellsY = 20 NumberOfCellsY = 20
# Grading and refinement of the mesh in y direction # Grading and refinement of the mesh in y direction
GradingFactorY = 1.1 GradingFactorY = 1.1
RefineTopY = false RefineTopY = false
[Output] [Output]
NameFF = stokes2cni NameFF = stokes2cni
NamePM = 2p2cni NamePM = 2p2cni
#Frequency of restart file, flux and VTK output #Frequency of restart file, flux and VTK output
FreqRestart = 1000 # how often restart files are written out FreqRestart = 1000 # how often restart files are written out
FreqOutput = 10 # frequency of VTK output FreqOutput = 10 # frequency of VTK output
FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy) FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy)
FreqFluxOutput = 1000 # frequency of detailed flux output FreqFluxOutput = 1000 # frequency of detailed flux output
FreqVaporFluxOutput = 2 # frequency of summarized flux output FreqVaporFluxOutput = 2 # frequency of summarized flux output
LiveEvaporationRates = false # plot evaporation rates using gnuplot interface LiveEvaporationRates = false # plot evaporation rates using gnuplot interface
[Stokes] [Stokes]
StabilizationAlpha = -1.0 StabilizationAlpha = -1.0
[FreeFlow] [FreeFlow]
UseDirichletBC = true # dirichlet values are set at the inflow boundary UseDirichletBC = true # dirichlet values are set at the inflow boundary
RefVelocity = 3.5 # [m/s] RefVelocity = 3.5 # [m/s]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefMassfrac = 0.008 # [-] RefMassfrac = 0.008 # [-]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
SinusVelAmplitude = 0.0 # [m/s] SinusVelAmplitude = 0.0 # [m/s]
SinusVelPeriod = 3600 # [s] SinusVelPeriod = 3600 # [s]
SinusPressureAmplitude = 0.0 # [Pa] SinusPressureAmplitude = 0.0 # [Pa]
SinusPressurePeriod = 3600 # [s] SinusPressurePeriod = 3600 # [s]
SinusConcentrationAmplitude = 0.0 # [-] SinusConcentrationAmplitude = 0.0 # [-]
SinusConcentrationPeriod = 3600 # [s] SinusConcentrationPeriod = 3600 # [s]
SinusTemperatureAmplitude = 0.0 # [K] SinusTemperatureAmplitude = 0.0 # [K]
SinusTemperaturePeriod = 3600 # [s] SinusTemperaturePeriod = 3600 # [s]
[BoundaryLayer] [BoundaryLayer]
Model = 0 Model = 0
[MassTransfer] [MassTransfer]
Model = 0 Model = 0
[PorousMedium] [PorousMedium]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
InitialSw = 0.98 # [-] InitialSw = 0.98 # [-]
[SpatialParams] [SpatialParams]
AlphaBJ = 1.0 # [-] AlphaBJ = 1.0 # [-]
Permeability = 2.65e-10 # [m^2] Permeability = 2.65e-10 # [m^2]
Porosity = 0.41 # [-] Porosity = 0.41 # [-]
Swr = 0.005 # [-] Swr = 0.005 # [-]
Snr = 0.01 # [-] Snr = 0.01 # [-]
VgAlpha = 6.371e-4 # [1/Pa] VgAlpha = 6.371e-4 # [1/Pa]
VgN = 8.0 # [-] VgN = 8.0 # [-]
LambdaSolid = 5.26 # [W/(m*K)] LambdaSolid = 5.26 # [W/(m*K)]
[Newton] [Newton]
MaxRelativeShift = 1e-5 MaxRelativeShift = 1e-5
TargetSteps = 8 TargetSteps = 8
MaxSteps = 12 MaxSteps = 12
WriteConvergence = false WriteConvergence = false
MaxTimeStepDivisions = 20 MaxTimeStepDivisions = 20
[LinearSolver] [LinearSolver]
Verbosity = 0 Verbosity = 0
test/multidomain/2cnizeroeq2p2cni/test_2cnizeroeq2p2cni.input
+ 85
− 85
View file @ f7aa645e
[TimeManager] [TimeManager]
DtInitial = 5e-5 # [s] DtInitial = 5e-5 # [s]
MaxTimeStepSize = 360 # [s] MaxTimeStepSize = 360 # [s]
TEnd = 3600 # [s] TEnd = 3600 # [s]
EpisodeLength = 14400 # [s] # 14400s = 4h EpisodeLength = 14400 # [s] # 14400s = 4h
[Grid] [Grid]
LowerLeftX = 0.0 # [m] LowerLeftX = 0.0 # [m]
UpperRightX = 0.75 # [m] UpperRightX = 0.75 # [m]
LowerLeftY = 0.0 # [m] LowerLeftY = 0.0 # [m]
UpperRightY = 0.75 # [m] UpperRightY = 0.75 # [m]
# Number of elements in x-, y-direction # Number of elements in x-, y-direction
NumberOfCellsX = 15 NumberOfCellsX = 15
NumberOfCellsY = 30 NumberOfCellsY = 30
# Grading and refinement of the mesh # Grading and refinement of the mesh
GradingFactorY = 1.16 GradingFactorY = 1.16
RefineTopY = true RefineTopY = true
# Position information # Position information
NoDarcyX1 = 0.25 # [m] # Beginning of PM below NoDarcyX1 = 0.25 # [m] # Beginning of PM below
NoDarcyX2 = 0.5 # [m] # End of PM below NoDarcyX2 = 0.5 # [m] # End of PM below
RunUpDistanceX1 = 0.251 # [m] # Beginning of without Coupling to PM (x-coordinate) RunUpDistanceX1 = 0.251 # [m] # Beginning of without Coupling to PM (x-coordinate)
RunUpDistanceX2 = 0.499 # [m] # End of without Coupling to PM (x-coordinate) RunUpDistanceX2 = 0.499 # [m] # End of without Coupling to PM (x-coordinate)
InterfacePosY = 0.25 # [m] # Vertical position of coupling interface InterfacePosY = 0.25 # [m] # Vertical position of coupling interface
[Output] [Output]
NameFF = zeroeq2cni NameFF = zeroeq2cni
NamePM = 2p2cni NamePM = 2p2cni
# Frequency of restart file, flux and VTK output # Frequency of restart file, flux and VTK output
FreqRestart = 5 # how often restart files are written out FreqRestart = 5 # how often restart files are written out
FreqOutput = 5 # 10 # frequency of VTK output FreqOutput = 5 # 10 # frequency of VTK output
FreqMassOutput = 5 # 20 # frequency of mass and evaporation rate output (Darcy) FreqMassOutput = 5 # 20 # frequency of mass and evaporation rate output (Darcy)
[FreeFlow] [FreeFlow]
RefVelocity = 1.0 # [m/s] RefVelocity = 1.0 # [m/s]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefMassfrac = 0.008 # [-] RefMassfrac = 0.008 # [-]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
[BoundaryLayer] [BoundaryLayer]
Model = 0 # disable boundary layer models Model = 0 # disable boundary layer models
[MassTransfer] [MassTransfer]
Model = 0 # disable mass transfer models Model = 0 # disable mass transfer models
[PorousMedium] [PorousMedium]
RefPressurePM = 1e5 # [Pa] RefPressurePM = 1e5 # [Pa]
RefTemperaturePM = 298.15 # [K] RefTemperaturePM = 298.15 # [K]
RefSw = 0.98 # [-] RefSw = 0.98 # [-]
[SpatialParams] [SpatialParams]
AlphaBJ = 1.0 # [-] AlphaBJ = 1.0 # [-]
Permeability = 2.65e-10 # [m^2] Permeability = 2.65e-10 # [m^2]
Porosity = 0.41 # [-] Porosity = 0.41 # [-]
Swr = 0.005 # [-] Swr = 0.005 # [-]
Snr = 0.01 # [-] Snr = 0.01 # [-]
VgAlpha = 6.371e-4 # [1/Pa] VgAlpha = 6.371e-4 # [1/Pa]
VgN = 8.0 # [-] VgN = 8.0 # [-]
ThermalConductivitySolid = 5.26 # [W/(m*K)] ThermalConductivitySolid = 5.26 # [W/(m*K)]
[Newton] [Newton]
MaxRelativeShift = 1e-5 MaxRelativeShift = 1e-5
TargetSteps = 8 TargetSteps = 8
MaxSteps = 12 MaxSteps = 12
WriteConvergence = false WriteConvergence = false
[LinearSolver] [LinearSolver]
Verbosity = 0 Verbosity = 0
[ZeroEq] [ZeroEq]
# Eddy Viscosity Models # Eddy Viscosity Models
# 0 = none # 0 = none
# 1 = Prandtl # 1 = Prandtl
# 2 = modified Van Driest # 2 = modified Van Driest
# 3 = Baldwin Lomax # 3 = Baldwin Lomax
EddyViscosityModel = 2 EddyViscosityModel = 2
# Eddy Diffusivity and Eddy Conductivity Models # Eddy Diffusivity and Eddy Conductivity Models
# 0 = none # 0 = none
# 1 = Reynolds analogy # 1 = Reynolds analogy
# 2 = modified Van Driest # 2 = modified Van Driest
# 3 = Deissler # 3 = Deissler
# 4 = Meier and Rotta # 4 = Meier and Rotta
EddyDiffusivityModel = 3 EddyDiffusivityModel = 3
EddyConductivityModel = 3 EddyConductivityModel = 3
BBoxMinSandGrainRoughness = 0.0 # [m] BBoxMinSandGrainRoughness = 0.0 # [m]
BBoxMaxSandGrainRoughness = 0.0 # [m] BBoxMaxSandGrainRoughness = 0.0 # [m]
test/multidomain/2cstokes2p2c/test_2cstokes2p2c.input
+ 80
− 80
View file @ f7aa645e
[TimeManager] [TimeManager]
DtInitial = 5e-1 # [s] DtInitial = 5e-1 # [s]
MaxTimeStepSize = 360 # [s] MaxTimeStepSize = 360 # [s]
InitTime = 0 # [s] Initialization time without coupling InitTime = 0 # [s] Initialization time without coupling
TEnd = 7200 # [s] TEnd = 7200 # [s]
EpisodeLength = 3600 # [s] EpisodeLength = 3600 # [s]
[Grid] [Grid]
# File = grids/test_2cstokes2p2c.dgf # unused, if interface mesh creator is used # File = grids/test_2cstokes2p2c.dgf # unused, if interface mesh creator is used
LowerLeftX = 0.0 # [m] LowerLeftX = 0.0 # [m]
UpperRightX = 0.25 # [m] UpperRightX = 0.25 # [m]
LowerLeftY = 0.0 # [m] LowerLeftY = 0.0 # [m]
UpperRightY = 0.5 # [m] UpperRightY = 0.5 # [m]
RunUpDistanceX = 0.0 # [m] Horizontal position without coupling to PM RunUpDistanceX = 0.0 # [m] Horizontal position without coupling to PM
NoDarcyX = 0.0 # [m] Horizontal position without PM below NoDarcyX = 0.0 # [m] Horizontal position without PM below
InterfacePosY = 0.25 # [m] Vertical position of coupling interface InterfacePosY = 0.25 # [m] Vertical position of coupling interface
# Number of elements in x-, y-, z-direction # Number of elements in x-, y-, z-direction
NumberOfCellsX = 30 NumberOfCellsX = 30
NumberOfCellsY = 62 NumberOfCellsY = 62
# Grading and refinement of the mesh in y direction # Grading and refinement of the mesh in y direction
GradingFactorY = 1.13 GradingFactorY = 1.13
RefineTopY = false RefineTopY = false
[Output] [Output]
NameFF = stokes2c NameFF = stokes2c
NamePM = 2p2c NamePM = 2p2c
#Frequency of restart file, flux and VTK output #Frequency of restart file, flux and VTK output
FreqRestart = 1000 # how often restart files are written out FreqRestart = 1000 # how often restart files are written out
FreqOutput = 50 # frequency of VTK output FreqOutput = 50 # frequency of VTK output
FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy) FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy)
FreqFluxOutput = 1000 # frequency of detailed flux output FreqFluxOutput = 1000 # frequency of detailed flux output
FreqVaporFluxOutput = 2 # frequency of summarized flux output FreqVaporFluxOutput = 2 # frequency of summarized flux output
[Stokes] [Stokes]
StabilizationAlpha = -1.0 StabilizationAlpha = -1.0
[FreeFlow] [FreeFlow]
RefVelocity = 3.5 # [m/s] RefVelocity = 3.5 # [m/s]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefMassfrac = 0.008 # [-] RefMassfrac = 0.008 # [-]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
SinusVelAmplitude = 0.0 # [m/s] SinusVelAmplitude = 0.0 # [m/s]
SinusVelPeriod = 3600 # [s] SinusVelPeriod = 3600 # [s]
SinusPressureAmplitude = 0.0 # [Pa] SinusPressureAmplitude = 0.0 # [Pa]
SinusPressurePeriod = 3600 # [s] SinusPressurePeriod = 3600 # [s]
SinusConcentrationAmplitude = 0.0 # [-] SinusConcentrationAmplitude = 0.0 # [-]
SinusConcentrationPeriod = 3600 # [s] SinusConcentrationPeriod = 3600 # [s]
SinusTemperatureAmplitude = 0.0 # [K] SinusTemperatureAmplitude = 0.0 # [K]
SinusTemperaturePeriod = 3600 # [s] SinusTemperaturePeriod = 3600 # [s]
[BoundaryLayer] [BoundaryLayer]
Model = 0 Model = 0
[MassTransfer] [MassTransfer]
Model = 0 Model = 0
[PorousMedium] [PorousMedium]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
InitialSw = 0.98 # [-] InitialSw = 0.98 # [-]
[SpatialParams] [SpatialParams]
AlphaBJ = 1.0 # [-] AlphaBJ = 1.0 # [-]
Permeability = 2.65e-10 # [m^2] Permeability = 2.65e-10 # [m^2]
Porosity = 0.41 # [-] Porosity = 0.41 # [-]
Swr = 0.005 # [-] Swr = 0.005 # [-]
Snr = 0.01 # [-] Snr = 0.01 # [-]
VgAlpha = 6.371e-4 # [1/Pa] VgAlpha = 6.371e-4 # [1/Pa]
VgN = 6.9 # [-] VgN = 6.9 # [-]
LambdaSolid = 5.3 # [W/(m*K)] LambdaSolid = 5.3 # [W/(m*K)]
[Newton] [Newton]
MaxRelativeShift = 1e-5 MaxRelativeShift = 1e-5
TargetSteps = 8 TargetSteps = 8
MaxSteps = 12 MaxSteps = 12
WriteConvergence = false WriteConvergence = false
MaxTimeStepDivisions = 20 MaxTimeStepDivisions = 20
[LinearSolver] [LinearSolver]
Verbosity = 0 Verbosity = 0
test/multidomain/2cstokes2p2c/test_2cstokes2p2c_reference.input
+ 80
− 80
View file @ f7aa645e
[TimeManager] [TimeManager]
DtInitial = 5e-1 # [s] DtInitial = 5e-1 # [s]
MaxTimeStepSize = 240 # [s] MaxTimeStepSize = 240 # [s]
InitTime = 0 # [s] Initialization time without coupling InitTime = 0 # [s] Initialization time without coupling
TEnd = 7200 # [s] TEnd = 7200 # [s]
EpisodeLength = 3600 # [s] EpisodeLength = 3600 # [s]
[Grid] [Grid]
# File = grids/test_2cstokes2p2c.dgf # unused, if interface mesh creator is used # File = grids/test_2cstokes2p2c.dgf # unused, if interface mesh creator is used
LowerLeftX = 0.0 # [m] LowerLeftX = 0.0 # [m]
UpperRightX = 0.25 # [m] UpperRightX = 0.25 # [m]
LowerLeftY = 0.0 # [m] LowerLeftY = 0.0 # [m]
UpperRightY = 0.5 # [m] UpperRightY = 0.5 # [m]
RunUpDistanceX = 0.0 # [m] Horizontal position without coupling to PM RunUpDistanceX = 0.0 # [m] Horizontal position without coupling to PM
NoDarcyX = 0.0 # [m] Horizontal position without PM below NoDarcyX = 0.0 # [m] Horizontal position without PM below
InterfacePosY = 0.25 # [m] Vertical position of coupling interface InterfacePosY = 0.25 # [m] Vertical position of coupling interface
# Number of elements in x-, y-, z-direction # Number of elements in x-, y-, z-direction
NumberOfCellsX = 25 NumberOfCellsX = 25
NumberOfCellsY = 45 NumberOfCellsY = 45
# Grading and refinement of the mesh in y direction # Grading and refinement of the mesh in y direction
GradingFactorY = 1.1 GradingFactorY = 1.1
RefineTopY = false RefineTopY = false
[Output] [Output]
NameFF = stokes2c NameFF = stokes2c
NamePM = 2p2c NamePM = 2p2c
#Frequency of restart file, flux and VTK output #Frequency of restart file, flux and VTK output
FreqRestart = 1000 # how often restart files are written out FreqRestart = 1000 # how often restart files are written out
FreqOutput = 10 # frequency of VTK output FreqOutput = 10 # frequency of VTK output
FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy) FreqMassOutput = 2 # frequency of mass and evaporation rate output (Darcy)
FreqFluxOutput = 1000 # frequency of detailed flux output FreqFluxOutput = 1000 # frequency of detailed flux output
FreqVaporFluxOutput = 2 # frequency of summarized flux output FreqVaporFluxOutput = 2 # frequency of summarized flux output
[Stokes] [Stokes]
StabilizationAlpha = -1.0 StabilizationAlpha = -1.0
[FreeFlow] [FreeFlow]
RefVelocity = 3.5 # [m/s] RefVelocity = 3.5 # [m/s]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefMassfrac = 0.008 # [-] RefMassfrac = 0.008 # [-]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
SinusVelAmplitude = 0.0 # [m/s] SinusVelAmplitude = 0.0 # [m/s]
SinusVelPeriod = 3600 # [s] SinusVelPeriod = 3600 # [s]
SinusPressureAmplitude = 0.0 # [Pa] SinusPressureAmplitude = 0.0 # [Pa]
SinusPressurePeriod = 3600 # [s] SinusPressurePeriod = 3600 # [s]
SinusConcentrationAmplitude = 0.0 # [-] SinusConcentrationAmplitude = 0.0 # [-]
SinusConcentrationPeriod = 3600 # [s] SinusConcentrationPeriod = 3600 # [s]
SinusTemperatureAmplitude = 0.0 # [K] SinusTemperatureAmplitude = 0.0 # [K]
SinusTemperaturePeriod = 3600 # [s] SinusTemperaturePeriod = 3600 # [s]
[BoundaryLayer] [BoundaryLayer]
Model = 0 Model = 0
[MassTransfer] [MassTransfer]
Model = 0 Model = 0
[PorousMedium] [PorousMedium]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
InitialSw = 0.98 # [-] InitialSw = 0.98 # [-]
[SpatialParams] [SpatialParams]
AlphaBJ = 1.0 # [-] AlphaBJ = 1.0 # [-]
Permeability = 2.65e-10 # [m^2] Permeability = 2.65e-10 # [m^2]
Porosity = 0.41 # [-] Porosity = 0.41 # [-]
Swr = 0.005 # [-] Swr = 0.005 # [-]
Snr = 0.01 # [-] Snr = 0.01 # [-]
VgAlpha = 6.371e-4 # [1/Pa] VgAlpha = 6.371e-4 # [1/Pa]
VgN = 6.9 # [-] VgN = 6.9 # [-]
LambdaSolid = 5.26 # [W/(m*K)] LambdaSolid = 5.26 # [W/(m*K)]
[Newton] [Newton]
MaxRelativeShift = 1e-5 MaxRelativeShift = 1e-5
TargetSteps = 8 TargetSteps = 8
MaxSteps = 12 MaxSteps = 12
WriteConvergence = false WriteConvergence = false
MaxTimeStepDivisions = 20 MaxTimeStepDivisions = 20
[LinearSolver] [LinearSolver]
Verbosity = 0 Verbosity = 0
test/multidomain/2czeroeq2p2c/test_2czeroeq2p2c.input
+ 76
− 76
View file @ f7aa645e
[TimeManager] [TimeManager]
DtInitial = 5e-5 # [s] DtInitial = 5e-5 # [s]
MaxTimeStepSize = 900 # [s] MaxTimeStepSize = 900 # [s]
TEnd = 518400 # [s] # 518400 = 6 days TEnd = 518400 # [s] # 518400 = 6 days
EpisodeLength = 8e5 # [s] # 43200 EpisodeLength = 8e5 # [s] # 43200
[Grid] [Grid]
LowerLeftX = 0.0 # [m] LowerLeftX = 0.0 # [m]
UpperRightX = 0.5 # [m] UpperRightX = 0.5 # [m]
LowerLeftY = 0.0 # [m] LowerLeftY = 0.0 # [m]
UpperRightY = 0.75 # [m] UpperRightY = 0.75 # [m]
# Number of elements in x-, y-direction # Number of elements in x-, y-direction
NumberOfCellsX = 20 NumberOfCellsX = 20
NumberOfCellsY = 40 NumberOfCellsY = 40
# Grading and refinement of the mesh # Grading and refinement of the mesh
GradingFactorY = 1.16 GradingFactorY = 1.16
RefineTopY = false RefineTopY = false
# Position information # Position information
NoDarcyX = 0.25 # [m] # Horizontal position without PM below NoDarcyX = 0.25 # [m] # Horizontal position without PM below
RunUpDistanceX = 0.26 # [m] # Horizontal position without Coupling to PM RunUpDistanceX = 0.26 # [m] # Horizontal position without Coupling to PM
InterfacePosY = 0.25 # [m] # Vertical position of coupling interface InterfacePosY = 0.25 # [m] # Vertical position of coupling interface
[Output] [Output]
NameFF = zeroeq2c NameFF = zeroeq2c
NamePM = 2p2c NamePM = 2p2c
# Frequency of restart file, flux and VTK output # Frequency of restart file, flux and VTK output
FreqRestart = 50 # 500 # how often restart files are written out FreqRestart = 50 # 500 # how often restart files are written out
FreqOutput = 5 # 10 # frequency of VTK output FreqOutput = 5 # 10 # frequency of VTK output
FreqMassOutput = 5 # 20 # frequency of mass and evaporation rate output (Darcy) FreqMassOutput = 5 # 20 # frequency of mass and evaporation rate output (Darcy)
[FreeFlow] [FreeFlow]
RefVelocity = 3.5 # [m/s] RefVelocity = 3.5 # [m/s]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefMassfrac = 0.008 # [-] RefMassfrac = 0.008 # [-]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
[PorousMedium] [PorousMedium]
RefPressurePM = 1e5 # [Pa] RefPressurePM = 1e5 # [Pa]
RefTemperaturePM = 298.15 # [K] RefTemperaturePM = 298.15 # [K]
RefSw = 0.9 # [-] RefSw = 0.9 # [-]
[SpatialParams] [SpatialParams]
AlphaBJ = 1.0 # [-] AlphaBJ = 1.0 # [-]
Permeability = 2.65e-10 # [m^2] Permeability = 2.65e-10 # [m^2]
Porosity = 0.41 # [-] Porosity = 0.41 # [-]
Swr = 0.005 # [-] Swr = 0.005 # [-]
Snr = 0.01 # [-] Snr = 0.01 # [-]
VgAlpha = 6.371e-4 # [1/Pa] VgAlpha = 6.371e-4 # [1/Pa]
VgN = 6.9 # [-] VgN = 6.9 # [-]
# optional parameters # optional parameters
[Newton] [Newton]
MaxRelativeShift = 1e-5 MaxRelativeShift = 1e-5
TargetSteps = 10 TargetSteps = 10
MaxSteps = 15 MaxSteps = 15
WriteConvergence = false WriteConvergence = false
[LinearSolver] [LinearSolver]
Verbosity = 0 Verbosity = 0
[ZeroEq] [ZeroEq]
# Eddy Viscosity Models # Eddy Viscosity Models
# 0 = none # 0 = none
# 1 = Prandtl # 1 = Prandtl
# 2 = modified Van Driest # 2 = modified Van Driest
# 3 = Baldwin Lomax # 3 = Baldwin Lomax
EddyViscosityModel = 3 EddyViscosityModel = 3
# Eddy Diffusivity Models # Eddy Diffusivity Models
# 0 = none # 0 = none
# 1 = Reynolds analogy # 1 = Reynolds analogy
# 2 = modified Van Driest # 2 = modified Van Driest
# 3 = Deissler # 3 = Deissler
# 4 = Meier and Rotta # 4 = Meier and Rotta
EddyDiffusivityModel = 3 EddyDiffusivityModel = 3
BBoxMinSandGrainRoughness = 0.0 # [m] # 0.6e-3 BBoxMinSandGrainRoughness = 0.0 # [m] # 0.6e-3
BBoxMaxSandGrainRoughness = 0.0 # [m] # 0 BBoxMaxSandGrainRoughness = 0.0 # [m] # 0
test/multidomain/2czeroeq2p2c/test_2czeroeq2p2c_reference.input
+ 77
− 77
View file @ f7aa645e
[TimeManager] [TimeManager]
DtInitial = 5e-5 # [s] DtInitial = 5e-5 # [s]
MaxTimeStepSize = 900 # [s] MaxTimeStepSize = 900 # [s]
TEnd = 86400 # [s] TEnd = 86400 # [s]
EpisodeLength = 8e5 # [s] EpisodeLength = 8e5 # [s]
[Grid] [Grid]
LowerLeftX = 0.0 # [m] LowerLeftX = 0.0 # [m]
UpperRightX = 0.5 # [m] UpperRightX = 0.5 # [m]
LowerLeftY = 0.0 # [m] LowerLeftY = 0.0 # [m]
UpperRightY = 0.75 # [m] UpperRightY = 0.75 # [m]
# Number of elements in x-, y-direction # Number of elements in x-, y-direction
NumberOfCellsX = 10 NumberOfCellsX = 10
NumberOfCellsY = 25 NumberOfCellsY = 25
# Grading and refinement of the mesh # Grading and refinement of the mesh
GradingFactorY = 1.1 GradingFactorY = 1.1
RefineTopY = false RefineTopY = false
# Position information # Position information
NoDarcyX = 0.25 # [m] # Horizontal position without PM below NoDarcyX = 0.25 # [m] # Horizontal position without PM below
RunUpDistanceX = 0.26 # [m] # Horizontal position without Coupling to PM RunUpDistanceX = 0.26 # [m] # Horizontal position without Coupling to PM
InterfacePosY = 0.25 # [m] # Vertical position of coupling interface InterfacePosY = 0.25 # [m] # Vertical position of coupling interface
[Output] [Output]
NameFF = zeroeq2c NameFF = zeroeq2c
NamePM = 2p2c NamePM = 2p2c
# Frequency of restart file, flux and VTK output # Frequency of restart file, flux and VTK output
FreqRestart = 100 # 500 # how often restart files are written out FreqRestart = 100 # 500 # how often restart files are written out
FreqOutput = 5 # 10 # frequency of VTK output FreqOutput = 5 # 10 # frequency of VTK output
FreqMassOutput = 5 # 20 # frequency of mass and evaporation rate output (Darcy) FreqMassOutput = 5 # 20 # frequency of mass and evaporation rate output (Darcy)
[FreeFlow] [FreeFlow]
RefVelocity = 3.5 # [m/s] RefVelocity = 3.5 # [m/s]
RefPressure = 1e5 # [Pa] RefPressure = 1e5 # [Pa]
RefMassfrac = 0.008 # [-] RefMassfrac = 0.008 # [-]
RefTemperature = 298.15 # [K] RefTemperature = 298.15 # [K]
[PorousMedium] [PorousMedium]
RefPressurePM = 1e5 # [Pa] RefPressurePM = 1e5 # [Pa]
RefTemperaturePM = 298.15 # [K] RefTemperaturePM = 298.15 # [K]
RefSw = 0.9 # [-] RefSw = 0.9 # [-]
[SpatialParams] [SpatialParams]
AlphaBJ = 1.0 # [-] AlphaBJ = 1.0 # [-]
Permeability = 2.65e-10 # [m^2] Permeability = 2.65e-10 # [m^2]
Porosity = 0.41 # [-] Porosity = 0.41 # [-]
Swr = 0.005 # [-] Swr = 0.005 # [-]
Snr = 0.01 # [-] Snr = 0.01 # [-]
VgAlpha = 6.371e-4 # [1/Pa] VgAlpha = 6.371e-4 # [1/Pa]
VgN = 6.9 # [-] VgN = 6.9 # [-]
# optional parameters # optional parameters
[Newton] [Newton]
MaxRelativeShift = 1e-5 MaxRelativeShift = 1e-5
TargetSteps = 10 TargetSteps = 10
MaxSteps = 15 MaxSteps = 15
WriteConvergence = false WriteConvergence = false
MaxTimeStepDivisions = 20 MaxTimeStepDivisions = 20
[LinearSolver] [LinearSolver]
Verbosity = 0 Verbosity = 0
[ZeroEq] [ZeroEq]
# Eddy Viscosity Models # Eddy Viscosity Models
# 0 = none # 0 = none
# 1 = Prandtl # 1 = Prandtl
# 2 = modified Van Driest # 2 = modified Van Driest
# 3 = Baldwin Lomax # 3 = Baldwin Lomax
EddyViscosityModel = 3 EddyViscosityModel = 3
# Eddy Diffusivity Models # Eddy Diffusivity Models
# 0 = none # 0 = none
# 1 = Reynolds analogy # 1 = Reynolds analogy
# 2 = modified Van Driest # 2 = modified Van Driest
# 3 = Deissler # 3 = Deissler
# 4 = Meier and Rotta # 4 = Meier and Rotta
EddyDiffusivityModel = 3 EddyDiffusivityModel = 3
BBoxMinSandGrainRoughness = 0.0 # [m] # 0.6e-3 BBoxMinSandGrainRoughness = 0.0 # [m] # 0.6e-3
BBoxMaxSandGrainRoughness = 0.0 # [m] # 0 BBoxMaxSandGrainRoughness = 0.0 # [m] # 0
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