version 3.9-dev
Dumux::FrictionLawManning< VolumeVariables > Class Template Reference

Implementation of the friction law after Manning. More...

#include <dumux/material/fluidmatrixinteractions/frictionlaws/manning.hh>

Inheritance diagram for Dumux::FrictionLawManning< VolumeVariables >:

Description

template<typename VolumeVariables>
class Dumux::FrictionLawManning< VolumeVariables >

The LET mobility model is used to limit the friction for small water depths if a roughness height > 0.0 is provided (default roughnessHeight = 0.0).

Public Member Functions

 FrictionLawManning (const Scalar gravity, const Scalar manningN, const Scalar roughnessHeight=0.0)
 Constructor. More...
 
Dune::FieldVector< Scalar, 2 > bottomShearStress (const VolumeVariables &volVars) const final
 Compute the bottom shear stress. More...
 
Scalar limitRoughH (const Scalar roughnessHeight, const Scalar waterDepth, const Scalar eps=1.0e-12) const
 Limit the friction for small water depth. More...
 

Constructor & Destructor Documentation

◆ FrictionLawManning()

template<typename VolumeVariables >
Dumux::FrictionLawManning< VolumeVariables >::FrictionLawManning ( const Scalar  gravity,
const Scalar  manningN,
const Scalar  roughnessHeight = 0.0 
)
inline
Parameters
gravityGravity constant (in m/s^2)
manningNManning friction coefficient (in s/m^(1/3)
roughnessHeightroughness height for limiting (in m) default = 0.0

Member Function Documentation

◆ bottomShearStress()

template<typename VolumeVariables >
Dune::FieldVector< Scalar, 2 > Dumux::FrictionLawManning< VolumeVariables >::bottomShearStress ( const VolumeVariables &  volVars) const
inlinefinalvirtual
Parameters
volVarsVolume 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.

Returns
shear stress in N/m^2. First entry is the x-component, the second the y-component.

Implements Dumux::FrictionLaw< VolumeVariables >.

◆ limitRoughH()

template<typename VolumeVariables >
Scalar Dumux::FrictionLaw< VolumeVariables >::limitRoughH ( const Scalar  roughnessHeight,
const Scalar  waterDepth,
const Scalar  eps = 1.0e-12 
) const
inlineinherited

Compute a small artificial water depth that is added to the actual water depth to avoid extreme friction values which can occur for small water depths.

The function is called with a roughness height, which can be seen as roughness height of the surface (e.g. grain size). For a zero roughnessHeight the artificial water depth will be zero.

A water depth minUpperH (minUpperH = 2 * roughnessHeight) is defined for the limiting. Limiting is applied between both 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 parameters 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 chosen. Therefore the calculation of the mobility is simplified significantly.

Parameters
roughnessHeightroughness height of the representative structure (e.g. largest grain size).
waterDepthwater depth.
epsIf the roughness height falls below this threshold, this function returns zero.

< maximal mobility


The documentation for this class was generated from the following file: