#include <cmath>
#include <algorithm>
#include <utility>
#include <type_traits>
#include <array>
#include <dune/common/std/type_traits.hh>
#include <dune/common/exceptions.hh>
#include <dumux/common/parameters.hh>
#include <dumux/flux/fluxvariablescaching.hh>
#include <dumux/flux/shallowwater/fluxlimiterlet.hh>
Go to the source code of this file.
Compute the shallow water viscous momentum flux due to viscosity.
The viscous momentum flux
\[ \int \int_{V} \mathbf{\nabla} \cdot \nu_t h \mathbf{\nabla} \mathbf{u} dV \]
is re-written using Gauss' divergence theorem to:
\[ \int_{S_f} \nu_t h \mathbf{\nabla} \mathbf{u} \cdot \mathbf{n_f} dS \]
The effective kinematic viscosity \( \nu_t \) can be calculated by adding a vertical (Elder-like) and a horizontal (Smagorinsky-like) part. This enabled by setting "ShallowWater.UseMixingLengthTurbulenceModel" to "true".
For now the calculation of the shallow water viscous momentum flux is implemented strictly for 2D depth-averaged models (i.e. 3 equations).
Classes | |
class | Dumux::ShallowWaterViscousFlux< NumEqVector, > |
Compute the shallow water viscous momentum flux due to viscosity. More... | |
Namespaces | |
namespace | Dumux |