Dumux::FrictionLaw< VolumeVariables > Class Template Reference

Implementation of the abstract base class for friction laws. More...

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

Inheritance diagram for Dumux::FrictionLaw< VolumeVariables >:

Description

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

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 >	shearStress (const VolumeVariables &volVar) const =0
	Compute the shear stress. More...
Scalar	limitRoughH (const Scalar roughnessHeight, const Scalar waterDepth) const
	Limit the friction for small water depth. More...
virtual	~FrictionLaw ()=default

Constructor & Destructor Documentation

~FrictionLaw()

template<typename VolumeVariables >

virtual Dumux::FrictionLaw< VolumeVariables >::~FrictionLaw ( )

virtualdefault

Member Function Documentation

limitRoughH()

template<typename VolumeVariables >

Scalar Dumux::FrictionLaw< VolumeVariables >::limitRoughH ( const Scalar  roughnessHeight, const Scalar  waterDepth  ) const

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 friciton. So the friciton 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.

Parameters

roughnessHeight	roughness height of the representive structure (e.g. largest grain size).
waterDepth	water depth.

< maximal mobility

shearStress()

template<typename VolumeVariables >

virtual Dune::FieldVector< Scalar, 2 > Dumux::FrictionLaw< VolumeVariables >::shearStress ( const VolumeVariables &  volVar ) const

pure virtual

Compute the shear stress.

Parameters

volVar

Volume Variables.

Compute the shear stress due to friction. The shear stress is not a tensor as know from contiuums mechanics, but a force projected on an area. Therefore it is a vector with two entries.

Returns

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

Implemented in Dumux::FrictionLawManning< VolumeVariables >, Dumux::FrictionLawNikuradse< VolumeVariables >, and Dumux::FrictionLawNoFriction< VolumeVariables >.

---

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

• frictionlaw.hh