Hougen Piret formula