pydune.physics.sedtransport.transport_laws.cubic_transport_law#

cubic_transport_law(theta, theta_d, omega)[source]#

Cubic transport law \(q_{\rm sat}/Q = \Omega \sqrt{\theta}(\theta - \theta_{\rm th})\), from Duràn et al. 2011.

Parameters:
  • theta (scalar, numpy array) – Shield number.

  • theta_d (scalar, numpy array) – Threshold Shield number.

  • omega (scalar, numpy array) – Prefactor of the transport law.

Returns:

Sediment fluxes calculated elementwise using the cubic transport law.

Return type:

scalar, numpy array

Examples

>>> import numpy as np
>>> theta = np.random.random((2000, ))
>>> theta_d, omega = 0.0053, 7.8
>>> qsat = cubic_transport_law(theta, theta_d, omega)

References

[1] Durán, O., Claudin, P., & Andreotti, B. (2011). On aeolian transport: Grain-scale interactions,

dynamical mechanisms and scaling laws. Aeolian Research, 3(3), 243-270.