13#ifndef DUMUX_DISCRETIZATION_BOX_FOURIERS_LAW_HH 
   14#define DUMUX_DISCRETIZATION_BOX_FOURIERS_LAW_HH 
   25template<
class TypeTag, 
class DiscretizationMethod>
 
   32template <
class TypeTag>
 
   38    using FVElementGeometry = 
typename GridGeometry::LocalView;
 
   39    using SubControlVolumeFace = 
typename GridGeometry::SubControlVolumeFace;
 
   44    using Element = 
typename GridView::template Codim<0>::Entity;
 
   54    static Scalar flux(
const Problem& problem,
 
   55                       const Element& element,
 
   56                       const FVElementGeometry& fvGeometry,
 
   57                       const ElementVolumeVariables& elemVolVars,
 
   58                       const SubControlVolumeFace& scvf,
 
   59                       const ElementFluxVariablesCache& elemFluxVarsCache)
 
   62        const auto& insideScv = fvGeometry.scv(scvf.insideScvIdx());
 
   63        const auto& outsideScv = fvGeometry.scv(scvf.outsideScvIdx());
 
   64        const auto& insideVolVars = elemVolVars[insideScv];
 
   65        const auto& outsideVolVars = elemVolVars[outsideScv];
 
   68        auto insideLambda = insideVolVars.effectiveThermalConductivity();
 
   69        auto outsideLambda = outsideVolVars.effectiveThermalConductivity();
 
   72        insideLambda *= insideVolVars.extrusionFactor();
 
   73        outsideLambda *= outsideVolVars.extrusionFactor();
 
   76        const auto lambda = 
faceTensorAverage(insideLambda, outsideLambda, scvf.unitOuterNormal());
 
   79        const auto& fluxVarsCache = elemFluxVarsCache[scvf];
 
   82        Dune::FieldVector<Scalar, GridView::dimensionworld> gradTemp(0.0);
 
   83        for (
auto&& scv : scvs(fvGeometry))
 
   84            gradTemp.axpy(elemVolVars[scv].temperature(), fluxVarsCache.gradN(scv.indexInElement()));
 
   87        return -1.0*
vtmv(scvf.unitOuterNormal(), lambda, gradTemp)*Extrusion::area(fvGeometry, scvf);
 
Traits extracting the public Extrusion type from T Defaults to NoExtrusion if no such type is found.
Definition extrusion.hh:155
forward declaration of the method-specific implementation
Definition flux/box/fourierslaw.hh:26
Defines all properties used in Dumux.
Helper classes to compute the integration elements.
A free function to average a Tensor at an interface.
Dune::DenseMatrix< MAT >::value_type vtmv(const Dune::DenseVector< V1 > &v1, const Dune::DenseMatrix< MAT > &M, const Dune::DenseVector< V2 > &v2)
Evaluates the scalar product of a vector v2, projected by a matrix M, with a vector v1.
Definition math.hh:880
typename GetProp< TypeTag, Property >::type GetPropType
get the type alias defined in the property
Definition propertysystem.hh:296
Define some often used mathematical functions.
The available discretization methods in Dumux.
Scalar faceTensorAverage(const Scalar T1, const Scalar T2, const Dune::FieldVector< Scalar, dim > &normal)
Average of a discontinuous scalar field at discontinuity interface (for compatibility reasons with th...
Definition facetensoraverage.hh:29
typename Extrusion< T >::type Extrusion_t
Convenience alias for obtaining the extrusion type.
Definition extrusion.hh:166