[aspect-devel] S40RTS and thermal boundary layers

Marie Kajan marie.kajan at gmail.com
Thu Jul 21 09:28:32 PDT 2016

Hey all,

Quick question. I'm not sure if this is specific to how ASPECT works or
just my chosen parameters, but I'm using the S40RTS perturbation with the
simple compressible material model and can already see boundary layers in
my temperature profile at the first time step. Can anyone tell me where
these come from? Are they imposed by ASPECT due to my temperature boundary
conditions, and how?

I'll paste my parameter file minus output settings below. Thanks!



set Dimension                              = 3
set Start time                             = 0
set End time                               = 0
set Linear solver tolerance                = 1e-5
set Max nonlinear iterations               = 100
set Nonlinear solver scheme                = iterated IMPES
set Nonlinear solver tolerance             = 1e-5
set Number of cheap Stokes solver steps    = 0
set Output directory                       =
set Use years in output instead of seconds = true

subsection Boundary temperature model
  set Model name = spherical constant

  subsection Spherical constant
    set Inner temperature = 4000
    set Outer temperature = 300

subsection Discretization
  set Use locally conservative discretization = false

  subsection Stabilization parameters
    set alpha = 2
    set beta  = 0.117
    set cR    = 0.5

subsection Geometry model
   set Model name = spherical shell

  subsection Spherical shell
    set Inner radius = 3481000
    set Outer radius = 6371000

subsection Gravity model
  set Model name = radial earth-like

subsection Heating model
  set List of model names = adiabatic heating

  subsection Adiabatic heating
    set Use simplified adiabatic heating = true

subsection Initial conditions
  set Model name = S40RTS perturbation

  subsection S40RTS perturbation
    set Data directory                                               =
    set Initial condition file name                                  =
    set Remove degree 0 from perturbation                            = true
    set Spline knots depth file name                                 =
    set Thermal expansion coefficient in initial temperature scaling = 3e-5
    set Vs to density scaling                                        = 0.2

subsection Material model
  set Model name = simple compressible

subsection Mesh refinement
  set Initial adaptive refinement = 0
  set Initial global refinement   = 4

subsection Model settings
  set Fixed temperature boundary indicators   = inner, outer
  set Tangential velocity boundary indicators = inner, outer
  set Remove nullspace                        = net rotation
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20160721/200a66b3/attachment-0001.html>

More information about the Aspect-devel mailing list