# TEST for mesh
set Dimension = 3
set End time = 0
#set Use years in output instead of seconds = true
set Output directory = output-mesh-test-cac96
set Pressure normalization = surface
set Surface pressure = 0
set Nonlinear solver scheme = no Advection, no Stokes
# Here we specify the geometry of the domain, which is
# a spherical shell with inner radius of 3481km and
# outer radius of 6371km
subsection Geometry model
set Model name = spherical shell
subsection Spherical shell
set Inner radius = 59760000
set Outer radius = 63760000
set Cells along circumference = 96
end
end
# This section describes the gravity field, which is pointing
# towards the Earth's center with the same magnitude of 10 m/s^2
# everywhere
subsection Gravity model
set Model name = radial constant
subsection Radial constant
set Magnitude = 10
end
end
# This section prescribes the initial condition in the temperature
subsection Initial temperature model
set Model name = function
subsection Function
set Function expression = 0
end
end
# The material model is based on the simple material model, which assumes
# a constant density, and other parameters as stated below.
subsection Material model
set Model name = simple
subsection Simple model
set Reference density = 3300
set Viscosity = 1e21
set Thermal expansion coefficient = 3e-5
set Reference temperature = 1600
set Thermal conductivity = 4.125
set Reference specific heat = 1250
end
end
# For this calculation we only do 2 global refinement steps. This resolution
# is too low to fully resolve the mantle flow, however it does capture
# the main features.
subsection Mesh refinement
set Initial global refinement = 4
set Skip solvers on initial refinement = true
set Skip setup initial conditions on initial refinement = true
end
# We assume free slip at the inner and outer boundary
subsection Boundary velocity model
set Tangential velocity boundary indicators = top, bottom
end
subsection Postprocess
set List of postprocessors = visualization, basic statistics
subsection Visualization
set Output format = vtu
end
end