[aspect-devel] questions

Thieulot, C. (Cedric) c.thieulot at uu.nl
Wed Apr 17 01:25:51 PDT 2013


Hi all, 
I have two unrelated questions:

- 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 ?

- 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 ?

cheers, 

Cedric.




More information about the Aspect-devel mailing list