Robust dynamic optimization in heterogeneous multiscale catalytic flow reactors using polynomial chaos expansion