Compute the shallow water viscous momentum flux due to viscosity. More...
#include <dumux/flux/shallowwaterviscousflux.hh>
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).
Public Types | |
using | Cache = FluxVariablesCaching::EmptyDiffusionCache |
using | CacheFiller = FluxVariablesCaching::EmptyCacheFiller |
Static Public Member Functions | |
template<class Problem , class FVElementGeometry , class ElementVolumeVariables > | |
static NumEqVector | flux (const Problem &problem, const typename FVElementGeometry::GridGeometry::GridView::template Codim< 0 >::Entity &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const typename FVElementGeometry::SubControlVolumeFace &scvf) |
Compute the viscous momentum flux contribution from the interface shear stress. More... | |
using Dumux::ShallowWaterViscousFlux< NumEqVector, >::Cache = FluxVariablesCaching::EmptyDiffusionCache |
using Dumux::ShallowWaterViscousFlux< NumEqVector, >::CacheFiller = FluxVariablesCaching::EmptyCacheFiller |