Robust optimization of a multiscale heterogeneous catalytic reactor system with spatially-varying uncertainty descriptions using polynomial chaos expansions