viscosity_refinement.prm
# Test the viscosity refinement criterion
set Adiabatic surface temperature = 1623 # default: 0
set Dimension = 2
set End time = 0
subsection Boundary temperature model
set List of model names = initial temperature
subsection Initial temperature
set Maximal temperature = 3773
set Minimal temperature = 273
end
end
subsection Geometry model
set Model name = box # default:
subsection Box
set X extent = 500000
set Y extent = 500000
end
end
subsection Gravity model
set Model name = vertical # default:
subsection Vertical
set Magnitude = 10.0
end
end
subsection Initial temperature model
set Model name = function # default:
subsection Function
set Function constants = l=250000
set Function expression = if(x < l, 1873, 1623)
set Variable names = x,y
end
end
subsection Material model
set Model name = simple # default:
subsection Simple model
set Thermal viscosity exponent = 5.0
set Reference temperature = 1750
end
end
subsection Mesh refinement
set Coarsening fraction = 0.05
set Refinement fraction = 0.3
set Initial adaptive refinement = 3 # default: 2
set Initial global refinement = 2 # default: 2
set Strategy = viscosity
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
set Fixed temperature boundary indicators = 0,3 # default:
end
subsection Boundary velocity model
set Tangential velocity boundary indicators = 0,1,2,3
end
subsection Boundary velocity model
set Zero velocity boundary indicators = # default:
end
subsection Postprocess
set List of postprocessors =
end