Navigation

Operators and Keywords

Function List:

C++ API

Function File: [A] = bim2a_axisymmetric_advection_upwind (mesh, beta)

Build the Upwind stabilized stiffness matrix for an advection problem in cylindrical coordinates with axisymmetric configuration.

The equation taken into account is:

1/r * d/dr (r * beta_r u) + d/dz (beta_z u) = f

where beta is an element-wise constant vector function.

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: bim2a_axisymmetric_rhs, bim2a_axisymmetric_reaction, bim2a_axisymmetric_advection_diffusion, bim2c_mesh_properties.

Package: bim