pydune.physics.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.
