Implementation of the abstract base class for friction laws. More...
#include <dumux/material/fluidmatrixinteractions/frictionlaws/frictionlaw.hh>
Implementation of the abstract base class for friction laws.
The LET mobility model is used to limit the friction for small water depths.
Public Member Functions | |
virtual Dune::FieldVector< Scalar, 2 > | bottomShearStress (const VolumeVariables &volVars) const =0 |
Compute the bottom shear stress. More... | |
virtual Dune::FieldVector< Scalar, 2 > | shearStress (const VolumeVariables &volVars) const |
Scalar | limitRoughH (const Scalar roughnessHeight, const Scalar waterDepth) const |
Limit the friction for small water depth. More... | |
virtual | ~FrictionLaw ()=default |
|
virtualdefault |
|
pure virtual |
Compute the bottom shear stress.
volVars | Volume variables |
Compute the bottom shear stress due to bottom friction. The bottom shear stress is a projection of the shear stress tensor onto the river bed. It can therefore be represented by a (tangent) vector with two entries.
Implemented in Dumux::FrictionLawManning< VolumeVariables >, Dumux::FrictionLawNikuradse< VolumeVariables >, and Dumux::FrictionLawNoFriction< VolumeVariables >.
|
inline |
Limit the friction for small water depth.
We define a water depth minUpperH. If the water depth is smaller, we start to limit the friction. So the friction term get's not extreme large for small water depths.
----------------------— minUpperH --------—
---------------------—roughnessHeight ------------— /\ /\ roughness /grain\ ----------------------------—bottom ---------------— /////////////////////////////////////////////////
For the limiting the LET model is used, which is usually applied in the porous media flow to limit the permeability due to the saturation. It employs the three empirical paramaters L, E and T, which describe the limiting curve (mobility).
auto mobility = (mobility_max * pow(sw,L))/(pow(sw,L) + E * pow(1.0-sw,T));
For the limitation of the roughness height L = 0.0, T = 2.0 and E = 1.0 are choosen. Therefore the calculation of the mobility is simplified significantly.
roughnessHeight | roughness height of the representative structure (e.g. largest grain size). |
waterDepth | water depth. |
< maximal mobility
|
inlinevirtual |