Commit eecfcdde authored by Timo Koch's avatar Timo Koch
Browse files

[fuelcell] Reduce test time

At 10s we almost reached equilibrium already.
Turn off line search--not needed--gained a few seconds.
Change input file to unix style line endings.
parent 12053fc6
......@@ -4,8 +4,8 @@ add_dumux_test(fuelcell fuelcell fuelcell.cc
python ${dumux_INCLUDE_DIRS}/bin/runtest.py
--script fuzzy
--files ${CMAKE_SOURCE_DIR}/lecture/references/fuelcell-reference.vtu
${CMAKE_CURRENT_BINARY_DIR}/fuelcell-00041.vtu
--command "${CMAKE_CURRENT_BINARY_DIR}/fuelcell")
${CMAKE_CURRENT_BINARY_DIR}/fuelcell-00017.vtu
--command "${CMAKE_CURRENT_BINARY_DIR}/fuelcell -TimeManager.TEnd 10")
# headers for installation and headercheck
install(FILES
......
[TimeManager]
DtInitial = 1e-2 # small because of Acosta pc law
MaxTimeStepSize = 20
TEnd = 20 #equilibrium time
[Grid]
NumberOfCellsX = 35
NumberOfCellsY = 140
LowerLeftX = 0.0
LowerLeftY = 0.0
UpperRightX = 5.0e-4
UpperRightY = 2.0e-3
[Problem]
Name = fuelcell
NTemperature = 3 # [-] number of temperature table entries
NPressure = 200 # [-] number of pressure table entries
PressureLow = 1e5 # [Pa] lower pressure limit for tabularization
PressureHigh = 3e7 # [Pa] upper pressure limit for tabularization
TemperatureLow = 300.15 # 330.15 # [Pa] lower temperature limit for tabularization
TemperatureHigh = 400.15 # 360.15 # [Pa] upper temperature limit for tabularization
InitialTemperature = 343.15 #[K]
ReactionLayerWidth = 0.00005 #[m]
[Vtk]
AddVelocity = true
[ElectroChemistry]
# parameters from Acosta paper
SpecificResistance = 0.25 #[Ohm*cm^2]
ReversibleVoltage = 1.191 #[V]
CellVoltage = 0.5 #[V]
ThermoneutralVoltage = 1.4836 #[V]
#Parameters to calculate exchange current density for reference conditions
RefCurrentDensity = 1.87e-8 #[A/cm^2]
RefO2PartialPressure = 5.0e5 #[Pa]
RefTemperature = 353.15 #[K]
ActivationBarrier = 73.0e3 #[J/mol]
TransferCoefficient = 0.5 #[] (alpha)
NumElectrons = 2 #[] number of electrons tranferred in reaction
TransportNumberH20 = 0.2327 #[] tH2O electro-osmotic drag oefficient
pO2Inlet = 35730.75 #[Pa] partial pressure of oxygen at gas channel inlet
SurfaceIncreasingFactor = 60 # fv
MaxIterations = 300
[Newton]
MaxRelativeShift = 1e-8
TargetSteps = 10
MaxSteps = 15
WriteConvergence = false
UseLineSearch = 1
[LinearSolver]
ResidualReduction = 1e-8
Verbosity = 0
MaxIterations = 100
[OperationalConditions]
GasPressureInlet1 = 2.013e5 # pressure at the upper inlet
GasPressureInlet2 = 2.013e5 # pressure at the lower inlet #2.013e5
LiquidWaterSaturationInlet = 0.055
#LiquidWaterSaturationInlet = 0.41 #drainage
[Output]
FreqOutput = 1 # frequency of VTK output
[TimeManager]
DtInitial = 1e-2 # small because of Acosta pc law
MaxTimeStepSize = 20
TEnd = 20 #equilibrium time
[Grid]
NumberOfCellsX = 35
NumberOfCellsY = 140
LowerLeftX = 0.0
LowerLeftY = 0.0
UpperRightX = 5.0e-4
UpperRightY = 2.0e-3
[Problem]
Name = fuelcell
NTemperature = 3 # [-] number of temperature table entries
NPressure = 200 # [-] number of pressure table entries
PressureLow = 1e5 # [Pa] lower pressure limit for tabularization
PressureHigh = 3e7 # [Pa] upper pressure limit for tabularization
TemperatureLow = 300.15 # 330.15 # [Pa] lower temperature limit for tabularization
TemperatureHigh = 400.15 # 360.15 # [Pa] upper temperature limit for tabularization
InitialTemperature = 343.15 #[K]
ReactionLayerWidth = 0.00005 #[m]
[Vtk]
AddVelocity = true
[ElectroChemistry]
# parameters from Acosta paper
SpecificResistance = 0.25 #[Ohm*cm^2]
ReversibleVoltage = 1.191 #[V]
CellVoltage = 0.5 #[V]
ThermoneutralVoltage = 1.4836 #[V]
#Parameters to calculate exchange current density for reference conditions
RefCurrentDensity = 1.87e-8 #[A/cm^2]
RefO2PartialPressure = 5.0e5 #[Pa]
RefTemperature = 353.15 #[K]
ActivationBarrier = 73.0e3 #[J/mol]
TransferCoefficient = 0.5 #[] (alpha)
NumElectrons = 2 #[] number of electrons tranferred in reaction
TransportNumberH20 = 0.2327 #[] tH2O electro-osmotic drag oefficient
pO2Inlet = 35730.75 #[Pa] partial pressure of oxygen at gas channel inlet
SurfaceIncreasingFactor = 60 # fv
MaxIterations = 300
[Newton]
MaxRelativeShift = 1e-5
[OperationalConditions]
GasPressureInlet1 = 2.013e5 # pressure at the upper inlet
GasPressureInlet2 = 2.013e5 # pressure at the lower inlet #2.013e5
LiquidWaterSaturationInlet = 0.055
#LiquidWaterSaturationInlet = 0.41 #drainage
[Output]
FreqOutput = 1 # frequency of VTK output
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