pydune.physics.turbulent_flow.geometrical_model.basal_shear#

basal_shear(x, y, alpha, A0, B0, AR, theta)[source]#

Calculate the basal shear stress over a two dimensional sinusoidal topography for an arbitrary wind direction.

Parameters:

  • x (array, scalar) – Streamwise coordinate, non-dimensional (\(kx\)).

  • y (array, scalar) – Spanwise coordinate, non-dimensional (\(ky\)).

  • alpha (array, scalar) – Dune orientation with respect to the perpendicular to the flow direction (in degree).

  • A0 (array, scalar) – value of the hydrodynamic coefficient for \(\alpha = 0\), i.e. for a dune orientation perpendicular to the flow direction.

  • B0 (array, scalar) – value of the hydrodynamic coefficient for \(\alpha = 0\), i.e. for a dune orientation perpendicular to the flow direction.

  • AR (array, scalar) – dune aspect ratio, \(k\xi\).

  • theta (array, scalar) – wind direction, in degree, in the trigonometric convention.

Returns:

  • Taux (array, scalar) – streamwise component of the non-dimensional shear stress.

  • Tauy (array, scalar) – spanwise component of the non-dimensional shear stress