12#ifndef DUMUX_ONEEQ_STAGGERED_FLUXVARIABLES_HH
13#define DUMUX_ONEEQ_STAGGERED_FLUXVARIABLES_HH
32template<
class TypeTag,
class BaseFluxVariables,
class DiscretizationMethod>
33class OneEqFluxVariablesImpl;
35template<
class TypeTag,
class BaseFluxVariables>
37:
public BaseFluxVariables
39 using ParentType = BaseFluxVariables;
43 using GridVolumeVariables =
typename GridVariables::GridVolumeVariables;
44 using ElementVolumeVariables =
typename GridVolumeVariables::LocalView;
45 using VolumeVariables =
typename GridVolumeVariables::VolumeVariables;
47 using GridFluxVariablesCache =
typename GridVariables::GridFluxVariablesCache;
48 using FluxVariablesCache =
typename GridFluxVariablesCache::FluxVariablesCache;
50 using GridFaceVariables =
typename GridVariables::GridFaceVariables;
51 using ElementFaceVariables =
typename GridFaceVariables::LocalView;
52 using FaceVariables =
typename GridFaceVariables::FaceVariables;
57 using FVElementGeometry =
typename GridGeometry::LocalView;
58 using SubControlVolumeFace =
typename GridGeometry::SubControlVolumeFace;
60 using GridView =
typename GridGeometry::GridView;
62 using Element =
typename GridView::template Codim<0>::Entity;
66 static constexpr int viscosityTildeEqIdx = Indices::viscosityTildeEqIdx - ModelTraits::dim();
74 const Element &element,
75 const FVElementGeometry& fvGeometry,
76 const ElementVolumeVariables& elemVolVars,
77 const ElementFaceVariables& elemFaceVars,
78 const SubControlVolumeFace &scvf,
79 const FluxVariablesCache& fluxVarsCache)
81 CellCenterPrimaryVariables flux = ParentType::computeMassFlux(problem, element, fvGeometry,
82 elemVolVars, elemFaceVars, scvf, fluxVarsCache);
85 auto upwindTermK = [](
const auto& volVars)
87 return volVars.viscosityTilde() * volVars.density();
90 flux[viscosityTildeEqIdx]
91 = ParentType::advectiveFluxForCellCenter(problem, fvGeometry, elemVolVars, elemFaceVars, scvf, upwindTermK);
94 const auto& insideScv = fvGeometry.scv(scvf.insideScvIdx());
95 const auto& outsideScv = fvGeometry.scv(scvf.outsideScvIdx());
96 const auto& insideVolVars = elemVolVars[scvf.insideScvIdx()];
97 const auto& outsideVolVars = elemVolVars[scvf.outsideScvIdx()];
100 Scalar insideCoeff = (insideVolVars.viscosity() + insideVolVars.viscosityTilde() * insideVolVars.density())
101 / insideVolVars.sigma();
102 Scalar outsideCoeff = (outsideVolVars.viscosity() + outsideVolVars.viscosityTilde() * outsideVolVars.density())
103 / outsideVolVars.sigma();
106 insideCoeff *= insideVolVars.extrusionFactor();
107 outsideCoeff *= outsideVolVars.extrusionFactor();
115 distance = (insideScv.dofPosition() - scvf.ipGlobal()).two_norm();
121 (outsideScv.dofPosition() - scvf.ipGlobal()).two_norm(),
122 (insideScv.dofPosition() - scvf.ipGlobal()).two_norm());
123 distance = (outsideScv.dofPosition() - insideScv.dofPosition()).two_norm();
126 const auto bcTypes = problem.boundaryTypes(element, scvf);
127 if (!(scvf.boundary() && (bcTypes.isOutflow(Indices::viscosityTildeEqIdx)
128 || bcTypes.isSymmetry())))
130 flux[viscosityTildeEqIdx]
132 * (insideVolVars.viscosityTilde() - outsideVolVars.viscosityTilde())
133 * Extrusion::area(fvGeometry, scvf);
CellCenterPrimaryVariables computeMassFlux(const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFaceVariables &elemFaceVars, const SubControlVolumeFace &scvf, const FluxVariablesCache &fluxVarsCache)
Computes the flux for the cell center residual.
Definition: freeflow/rans/oneeq/staggered/fluxvariables.hh:73
The flux variables class for the one-equation model by Spalart-Allmaras using the staggered grid disc...
Definition: freeflow/rans/oneeq/fluxvariables.hh:23
Defines all properties used in Dumux.
Helper classes to compute the integration elements.
Base class for the flux variables living on a sub control volume face.
constexpr Scalar arithmeticMean(Scalar x, Scalar y, Scalar wx=1.0, Scalar wy=1.0) noexcept
Calculate the (weighted) arithmetic mean of two scalar values.
Definition: math.hh:38
static ctype distance(const Dune::FieldVector< ctype, dimWorld > &a, const Dune::FieldVector< ctype, dimWorld > &b)
Compute the shortest distance between two points.
Definition: distance.hh:282
typename GetProp< TypeTag, Property >::type GetPropType
get the type alias defined in the property
Definition: propertysystem.hh:296
The available discretization methods in Dumux.
typename Extrusion< T >::type Extrusion_t
Convenience alias for obtaining the extrusion type.
Definition: extrusion.hh:166