mpfa/computetransmissibility.hh

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
//
// SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
// SPDX-License-Identifier: GPL-3.0-or-later
//
#ifndef DUMUX_DISCRETIZATION_CC_MPFA_COMPUTE_TRANSMISSIBILITY_HH
#define DUMUX_DISCRETIZATION_CC_MPFA_COMPUTE_TRANSMISSIBILITY_HH
 
#include <dune/common/typetraits.hh>
#include <dune/common/fvector.hh>
 
#include <dumux/common/math.hh>
#include <dumux/discretization/extrusion.hh>
 
namespace Dumux {
 
template<class EG, class IVSubControlVolume, class Tensor>
Dune::FieldVector<typename Tensor::field_type, IVSubControlVolume::myDimension>
computeMpfaTransmissibility(const EG& fvGeometry,
                            const IVSubControlVolume& scv,
                            const typename EG::SubControlVolumeFace& scvf,
                            const Tensor& t,
                            typename IVSubControlVolume::ctype extrusionFactor)
{
    Dune::FieldVector<typename Tensor::field_type, IVSubControlVolume::myDimension> wijk;
    for (unsigned int dir = 0; dir < IVSubControlVolume::myDimension; ++dir)
        wijk[dir] = vtmv(scvf.unitOuterNormal(), t, scv.nu(dir));
 
    using Extrusion = Extrusion_t<typename EG::GridGeometry>;
    wijk *= Extrusion::area(fvGeometry, scvf)*extrusionFactor;
    wijk /= scv.detX();
 
    return wijk;
}
 
template< class EG, class IVSubControlVolume, class Tensor, std::enable_if_t< Dune::IsNumber<Tensor>::value, int > = 1 >
Dune::FieldVector<Tensor, IVSubControlVolume::myDimension>
computeMpfaTransmissibility(const EG& fvGeometry,
                            const IVSubControlVolume& scv,
                            const typename EG::SubControlVolumeFace& scvf,
                            const Tensor& t,
                            typename IVSubControlVolume::ctype extrusionFactor)
{
    Dune::FieldVector<Tensor, IVSubControlVolume::myDimension> wijk;
    for (unsigned int dir = 0; dir < IVSubControlVolume::myDimension; ++dir)
        wijk[dir] = vtmv(scvf.unitOuterNormal(), t, scv.nu(dir));
 
    using Extrusion = Extrusion_t<typename EG::GridGeometry>;
    wijk *= Extrusion::area(fvGeometry, scvf)*extrusionFactor;
    wijk /= scv.detX();
 
    return wijk;
}
 
} // end namespace Dumux
 
#endif