Build the standard finite element stiffness matrix for a diffusion problem in cylindrical coordinates with axisymmetric configuration. Rotational symmetry is assumed with respect to be the vertical axis r=0. Only grids that DO NOT contain r=0 are admissible.
| |-------| OK |--|-----| NO! r=0 r=0
The equation taken into account is:
- 1/r * (r * epsilon * kappa ( u’ ))’ = f
where epsilon is an element-wise constant scalar function, while kappa is a piecewise linear conforming scalar function.
See also: bim1a_axisymmetric_rhs, bim1a_axisymmetric_reaction, bim1a_axisymmetric_advection_diffusion, bim2a_laplacian, bim3a_laplacian.
Package: bim