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