Wiener chaos expansion and numerical solutions of the HJM interest rate model