Critical heat flux prediction in rod bundles under upward low mass flux densities