Build the Upwind stabilized stiffness matrix for an advection problem in cylindrical coordinates with axisymmetric configuration.
The equation taken into account is:
1/r * (r * beta u)’ = f
where beta is an element-wise constant.
Instead of passing the vector field beta directly one can pass a piecewise linear conforming scalar function phi as the last input. In such case beta = grad phi is assumed.
If phi is a single scalar value beta is assumed to be 0 in the whole domain.
See also: bim1a_axisymmetric_advection_diffusion, bim1a_axisymmetric_rhs, bim1a_axisymmetric_reaction, bim1a_axisymmetric_laplacian.
Package: bim