[aspect-devel] Some interesting (pathological?!) behaviour of particles ... Do we need CVI here?

[Elbridge Gerry Puckett]

Hi Rene and Cedric,

Check out the SolKz benchmark as in the ASPECT manual, but with 
particles running around in the flow.

Do we need CVI here?

- G
-------------- next part --------------
# A description of the SolKZ benchmark for which a known solution
# is available. See the manual for more information.

set Additional shared libraries            = /home/egpuckett/Computations/2016_Summer_REU/SolKz/Original/SolKz_Library/libsolkz.so

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

set Dimension                              = 2

set Start time                             = 0
set End time                               = 250

set Output directory                       = output

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


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

subsection Geometry model
  set Model name = box

  subsection Box
    set X extent = 1
    set Y extent = 1
  end
end


subsection Model settings
  set Prescribed velocity boundary indicators =
  set Tangential velocity boundary indicators = left, right, bottom, top
  set Zero velocity boundary indicators       =
end


subsection Material model
  set Model name = SolKzMaterial
end


subsection Gravity model
  set Model name = vertical
end


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

subsection Boundary temperature model
  set Model name = box
end


subsection Initial conditions
  set Model name = perturbed box
end


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

subsection Discretization
  set Stokes velocity polynomial degree       = 2
  set Use locally conservative discretization = false
end


subsection Mesh refinement
  set Initial adaptive refinement              = 0
  set Initial global refinement                = 5
end


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

subsection Postprocess
  set List of postprocessors = SolKzPostprocessor, visualization, tracers

  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
  end

  subsection Tracers
    set Number of tracers         = 1024
    set Time between data output  = 0
    set Data output format        = vtu

    set Particle generator name = uniform box

    subsection Generator
      subsection Uniform box
        set Minimum x = 0.00
        set Maximum x = 1.0
        set Minimum y = 0.0
        set Maximum y = 1.0
      end
    end
  end
end

subsection Termination criteria
  set End step = 1
  set Termination criteria = end time
end