37                                                        const typename FVElementGeometry::SubControlVolumeFace& scvf,
 
   38                                                        const typename FVElementGeometry::SubControlVolume& scv,
 
   40                                                        typename FVElementGeometry::SubControlVolume::Traits::Scalar extrusionFactor)
 
   42    using GlobalPosition = 
typename FVElementGeometry::SubControlVolumeFace::Traits::GlobalPosition;
 
   43    GlobalPosition Knormal;
 
   44    T.mv(scvf.unitOuterNormal(), Knormal);
 
   46    auto distanceVector = scvf.ipGlobal();
 
   47    distanceVector -= scv.center();
 
   48    distanceVector /= distanceVector.two_norm2();
 
   50    return (Knormal*distanceVector) * extrusionFactor;
 
 
   70                                   const typename FVElementGeometry::SubControlVolumeFace& scvf,
 
   71                                   const typename FVElementGeometry::SubControlVolume& scv,
 
   73                                   typename FVElementGeometry::SubControlVolumeFace::Traits::Scalar extrusionFactor)
 
   75    auto distanceVector = scvf.ipGlobal();
 
   76    distanceVector -= scv.center();
 
   77    distanceVector /= distanceVector.two_norm2();
 
   79    return t * extrusionFactor * (distanceVector * scvf.unitOuterNormal());
 
 
Tensor::field_type computeTpfaTransmissibility(const FVElementGeometry &fvGeometry, const typename FVElementGeometry::SubControlVolumeFace &scvf, const typename FVElementGeometry::SubControlVolume &scv, const Tensor &T, typename FVElementGeometry::SubControlVolume::Traits::Scalar extrusionFactor)
Free function to evaluate the Tpfa transmissibility associated with the flux (in the form of flux = T...
Definition tpfa/computetransmissibility.hh:36