3.3.0
DUNE for Multi-{Phase, Component, Scale, Physics, ...} flow and transport in porous media
manning.hh
Go to the documentation of this file.
1// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2// vi: set et ts=4 sw=4 sts=4:
3/*****************************************************************************
4 * See the file COPYING for full copying permissions. *
5 * *
6 * This program is free software: you can redistribute it and/or modify *
7 * it under the terms of the GNU General Public License as published by *
8 * the Free Software Foundation, either version 3 of the License, or *
9 * (at your option) any later version. *
10 * *
11 * This program is distributed in the hope that it will be useful, *
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of *
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
14 * GNU General Public License for more details. *
15 * *
16 * You should have received a copy of the GNU General Public License *
17 * along with this program. If not, see <http://www.gnu.org/licenses/>. *
18 *****************************************************************************/
24#ifndef DUMUX_MATERIAL_FLUIDMATRIX_FRICTIONLAW_MANNING_HH
25#define DUMUX_MATERIAL_FLUIDMATRIX_FRICTIONLAW_MANNING_HH
26
27#include <algorithm>
28#include <cmath>
29#include <dune/common/math.hh>
30#include "frictionlaw.hh"
31
32namespace Dumux {
40template <typename VolumeVariables>
41class FrictionLawManning : public FrictionLaw<VolumeVariables>
42{
43 using Scalar = typename VolumeVariables::PrimaryVariables::value_type;
44public:
51 FrictionLawManning(const Scalar gravity, const Scalar manningN)
52 : gravity_(gravity), manningN_(manningN) {}
53
65 Dune::FieldVector<Scalar, 2> shearStress(const VolumeVariables& volVars) const final
66 {
67 using std::pow;
68 using Dune::power;
69 using std::hypot;
70
71 Dune::FieldVector<Scalar, 2> shearStress(0.0);
72
73 Scalar roughnessHeight = power(25.68/(1.0/manningN_),6);
74 roughnessHeight = this->limitRoughH(roughnessHeight, volVars.waterDepth());
75 const Scalar c = pow((volVars.waterDepth() + roughnessHeight),1.0/6.0) * 1.0/(manningN_);
76 const Scalar uv = hypot(volVars.velocity(0),volVars.velocity(1));
77
78 shearStress[0] = -gravity_/(c*c) * volVars.velocity(0) * uv;
79 shearStress[1] = -gravity_/(c*c) * volVars.velocity(1) * uv;
80
81 return shearStress;
82 }
83private:
84 Scalar gravity_;
85 Scalar manningN_;
86};
87
88} // end namespace Dumux
89
90#endif
Definition: adapt.hh:29
Implementation of the abstract base class for friction laws.
Definition: frictionlaw.hh:42
Scalar limitRoughH(const Scalar roughnessHeight, const Scalar waterDepth) const
Limit the friction for small water depth.
Definition: frictionlaw.hh:88
Implementation of the friction law after Manning.
Definition: manning.hh:42
Dune::FieldVector< Scalar, 2 > shearStress(const VolumeVariables &volVars) const final
Compute the shear stress.
Definition: manning.hh:65
FrictionLawManning(const Scalar gravity, const Scalar manningN)
Constructor.
Definition: manning.hh:51