internal_heating_statistics.prm
############### Global parameters
# This is a test for the radiogenic heating and the internal
# heating statistics postprocessor. Since density and radiogenic
# heating are constant there is an analytical solution, which
# is reproduced in this testcase.
set Dimension = 2
set Start time = 0
set End time = 1000
set Use years in output instead of seconds = false
set Output directory = radiogenic_heating
subsection Heating model
set List of model names = constant heating
subsection Constant heating
set Radiogenic heating rate = 1
end
end
subsection Geometry model
set Model name = box
subsection Box
set X extent = 1000000
set Y extent = 1000000
end
end
subsection Gravity model
set Model name = vertical
subsection Vertical
set Magnitude = 10.0
end
end
# The parameters below this comment were created by the update script
# as replacement for the old 'Model settings' subsection. They can be
# safely merged with any existing subsections with the same name.
subsection Boundary temperature model
# We only fix the temperature at the upper boundary, the other boundaries
# are isolating. To guarantuee a steady downward flow, we fix the velocity
# at the top and bottom, and set it to free slip on the sides.
set Fixed temperature boundary indicators = 3
end
subsection Boundary velocity model
set Prescribed velocity boundary indicators = 2:function, 3:function
end
subsection Boundary velocity model
set Tangential velocity boundary indicators = 0, 1
end
############### Boundary conditions
# We set the top temperature to T1=1000K.
subsection Boundary temperature model
set List of model names = box
subsection Box
set Top temperature = 1000
end
end
# We prescribe a constant downward flow.
subsection Boundary velocity model
subsection Function
set Function expression = 0;-1
set Variable names = x,y
end
end
subsection Initial temperature model
set Model name = function
subsection Function
set Function expression = 1000.0
set Variable names = x,y
end
end
subsection Material model
set Model name = simple
subsection Simple model
set Reference density = 3400
set Reference specific heat = 1000
set Reference temperature = 1000
set Thermal conductivity = 0.0
set Thermal expansion coefficient = 0.0
set Thermal viscosity exponent = 0.0
set Viscosity = 1e21
end
end
subsection Mesh refinement
set Initial adaptive refinement = 0
set Initial global refinement = 5
set Time steps between mesh refinement = 0
end
subsection Postprocess
set List of postprocessors = temperature statistics, heating statistics
subsection Visualization
set Interpolate output = false
set List of output variables = heating
set Time between graphical output = 0
end
end