move computePoreThroatRadius() to 2PFlow