[aspect-devel] VectorTools::integrate_difference fails to return an error that decreases under grid refinement
Elbridge Gerry Puckett
egpuckett at ucdavis.edu
Tue Oct 25 13:17:59 PDT 2016
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
More information about the Aspect-devel
mailing list