Computing with ordinal expressions like $(\omega^{\omega^2+1}+\omega)^{\omega+2}$ is not very different from handling polynomials, at least if we restrict to ordinals below $\epsilon_0$ written in Cantor Normal Form.