set Dimension = 2 set Use years in output instead of seconds = true subsection Termination criteria set Termination criteria = end step set End step = 4 end set Output directory = outputs/annulus_simple_Glob3Adapt2_newsolver # The following variables describe how the pressure should # be normalized. Here, we choose a zero average pressure # at the surface of the domain. The 'Surface pressure' and # 'Adiabatic surface temperature' are used to compute # the adiabatic reference profile. set Pressure normalization = surface set Surface pressure = 0 set Adiabatic surface temperature = 1573 set CFL number = 0.5 # Remove net rotation from solution subsection Nullspace removal set Remove nullspace = net rotation end ### SOLVER STUFF ### set Nonlinear solver scheme = single Advection, iterated defect correction Stokes subsection Solver parameters subsection Stokes solver parameters set Stokes solver type = block AMG set Number of cheap Stokes solver steps = 2000 set Use full A block as preconditioner = true set Linear solver tolerance = 1e-7 set Maximum number of expensive Stokes solver steps = 5000 end subsection Newton solver parameters set SPD safety factor = 0.9 set Maximum linear Stokes solver tolerance = 1e-2 set Use Eisenstat Walker method for Picard iterations = true end end ################### # 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 = 3481.0e3 set Outer radius = 6371.0e3 end end # This section specifies the temperature at the boundary of # the domain. Here we set the temperature to be constant, # but different from the reference temperature to approximate # boundary layers. subsection Boundary temperature model set Fixed temperature boundary indicators = top, bottom set List of model names = spherical constant subsection Spherical constant set Inner temperature = 2273 set Outer temperature = 273 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 subsection Initial temperature model set Model name = ascii data subsection Ascii data model set Data directory = /work2/08512/tg878111/stampede2/software/aspect_work/inputs/tomography/ set Data file name = annulus_T0.txt end end # The material model is based on the simple material model, which assumes # a constant viscosity, 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 = 1573 set Thermal conductivity = 4.125 set Reference specific heat = 1250 end end subsection Mesh refinement set Initial global refinement = 3 set Initial adaptive refinement = 2 set Refinement fraction = 0.2 set Strategy = minimum refinement function set Refinement criteria scaling factors = 1 set Time steps between mesh refinement = 0 subsection Minimum refinement function set Coordinate system = spherical set Function expression = ((r > (6371.e3-600.e3)) ? 5: \ (r > (6371.e3-1000.e3)) ? 4: \ 3) set Variable names = r,phi end end subsection Boundary velocity model set Tangential velocity boundary indicators = top, bottom # free slip end # We output the density, velocity, dynamic topography, geoid and heat flux density # for plotting. subsection Postprocess set List of postprocessors = velocity statistics, heat flux map, heat flux statistics, dynamic topography, visualization, basic statistics subsection Visualization set Output format = vtu set List of output variables = dynamic topography, density,viscosity set Time between graphical output = 1 # We only have dirichlet boundaries with tangential velocities, so we can # increase the output resolution as described in the documentation of the 'heat # flux map' postprocessor. subsection Heat flux map set Output point wise heat flux = true end end end