pydune.physics.dune.bedinstability_1D.spatial_wavenumber#

spatial_wavenumber(w, A0, B0, mu, r)[source]#

Dune instability spatial wavenumber - real part of the complex wavenumber \(k_{+}\). Note that here, \(\mathcal{A} = \mathcal{A}_{0}\) where \(\mathcal{B} = \mathcal{B}_{0} - 1/(r^{2}\mu)\), taking into account slope effects.

Parameters:

  • w (scalar, numpy array) – Non dimensional pulsation \(\omega\).

  • A0 (scalar, numpy array) – Hydrodynamic coefficient \(\mathcal{A}_{0}\) (in-phase).

  • B0 (scalar, numpy array) – Hydrodynamic coefficient \(\mathcal{B}_{0}\) (in-quadrature).

  • mu (scalar, numpy array) – Friction coefficient \(\mu\).

  • r (scalar, numpy array) – Velocity ratio \(u_{*}/u_{\rm d} = \sqrt{\theta/\theta_{d}}\)

Returns:

Spatial dune wavenumber \(k\) calculated elementwise.

Return type:

scalar, numpy array

Notes

Note that all quantities are made non dimensional:

  • length scales by the saturation length \(L_{\rm sat}\).

  • time scales by \(L_{\rm sat}^{2}/Q\), where \(Q\) is the characteristic flux.

Examples

>>> import numpy as np
>>> w = np.linspace(0.001, 1, 1000)
>>> A0, B0, mu, r = 3.5, 2, 0.63, 2
>>> k = spatial_wavenumber(w, A0, B0, mu, r)