[aspect-devel] VectorTools::integrate_difference fails to return an error that decreases under grid refinement

[Elbridge Gerry Puckett]

An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161025/09b47c84/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: egp_hvl_annulus.cc
Type: text/x-c++src
Size: 19735 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161025/09b47c84/attachment-0001.cc>
-------------- next part --------------
set Additional shared libraries            = ./libegp_hvl_annulus.so

############### Global parameters

set Dimension                              = 2

set Start time                             = 0

set End time                               = 0

set Output directory                       = output_annulus_1

#set Pressure normalization                 = volume
set Nonlinear solver scheme                = Stokes only

set CFL number                             = 0.8

set Linear solver tolerance                    = 1e-14

############### Parameters describing the model

subsection Geometry model
  set Model name = spherical shell

  subsection Spherical shell
    set Inner radius = .4
    set Outer radius = 1
    set Opening angle = 360
  end
end

subsection Model settings
# Set to tangential velocity boundary conditions.
# Note that the function specified in "Boundary velocity model" has no influence unless the following line is unmommneted 
# and the parameter  "Tangential velocity boundary indicators" is reset.
  set Prescribed velocity boundary indicators = 0: function, 1: function 
#  set Prescribed velocity boundary indicators =
  set Tangential velocity boundary indicators = 
  set Zero velocity boundary indicators       =
end

subsection Boundary velocity model
  subsection Function
    set Function constants  = pi=3.1415926
    set Variable names      = x,y
    set Function expression = -y;x
  end
end

subsection Material model
  set Model name = EGPHVLMaterial

  subsection EGPHVL
    # Do not change these set of parameters! EGP && HL
    set Viscosity jump = 1
    set Reference density = 1
  end
end


subsection Gravity model
  set Model name = radial constant
  subsection Radial constant
    set Magnitude = 0
  end

  subsection Function
    set Variable names      = x,y
    set Function expression = (1/(sqrt(x*x+y*y)))*x;(1/(sqrt(x*x+y*y)))*y
  end
end


############### Parameters describing the temperature field

subsection Boundary temperature model
  set Model name = box
end


subsection Initial conditions
  set Model name = function
  subsection Function
    set Function expression = 0
  end
end

subsection Adiabatic conditions model
  set Model name = function

  subsection Function
    set Function constants  =
    set Function expression = 0; 0
    set Variable names      = depth
  end

end

subsection Compositional fields
#  set Number of fields = 0
#  set Names of fields = density_comp #, viscosity_comp
#  set Compositional field methods = particles #, particles
  # set Mapped particle properties = density_comp:function1 #, viscosity_comp:function2
  # This is what was recommended by Rene on Friday, August 11, 2016
#  set Mapped particle properties = density_comp:function1
end

subsection Compositional initial conditions
  set Model name = function
  subsection Function
    set Variable names      = x,z
    set Function expression = 0
  end
end

############### Parameters describing the discretization

subsection Discretization
  set Stokes velocity polynomial degree       = 2
  set Use locally conservative discretization = false
#  set Use discontinuous composition discretization = true
end


subsection Mesh refinement
  set Coarsening fraction                      = 0
  set Initial adaptive refinement              = 0
  set Initial global refinement = 7
end


############### Parameters describing what to do with the solution

subsection Postprocess
  set List of postprocessors = visualization, EGPHVLPostprocessor
  
  subsection Visualization
    set Output format                 = vtu
    set Number of grouped files       = 1
    set Time between graphical output = 0
    set List of output variables = density, viscosity, partition, gravity
    set Interpolate output            = true    
  end
end

subsection Termination criteria
  set Checkpoint on termination = true
  set Termination criteria      = user request

  subsection User request
    set File name = end
  end
end

subsection Checkpointing
  #set Steps between checkpoint = 5
end