[aspect-devel] questions

Timo Heister heister at math.tamu.edu
Wed Apr 17 06:45:03 PDT 2013


> - during the integration process within an element (by means of a Gauss quadrature) how is the viscosity (and the density) computed at each Gauss point ? is it a constant within the element, or is it a different value for each quadrature point based
> on the local temp, press and composition fields ?

The viscosity is evaluated at each quadrature point. This is done in
the call  material_model->evaluate() in
Simulator<dim>::local_assemble_stokes_system.

> - turning now to the advection stabilisation parameters, and looking for instance at box.prm, I see three parameters: alpha, beta and cR.
> Now lookig at the GJI article about Aspect (kronbichler et al, 2012), section 3.2.6 of the article details the algorithm.
> As it turns out the alpha_max of the article seems to be beta in the code (a bit misleading ?)
> However, what are the equivalent of alpha and cR in the article ?

Yeah, the naming is somewhat misleading:
parameters.stabilization_beta = alpha_max in the paper
parameters.stabilization_c_R = alpha_E in the paper


--
Timo Heister
http://www.math.tamu.edu/~heister/


More information about the Aspect-devel mailing list