pydune.physics.dune.bedinstability_1D.temporal_pulsation#

temporal_pulsation(k, A0, B0, mu, r)[source]#

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

Parameters:

  • k (scalar, numpy array) – non dimensional wavenumber \(k\).

  • 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:

omega_r – temporal dune pulsation \(\omega_{\rm r}\).

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_{*\).

Examples

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