This　article　mainly　concerns　modeling　the　stochastic　input　and　its　propagation　in　incompressible　Navier-Stokes(N-S)　flow　simulations.　The　stochastic　input　is　represented　spectrally　by　employing　orthogonal　polynomial　functionals　front　the　Askey　scheme　as　trial　basis　to　represent　the　random　space.　A　standard　Galerkin　projection　is　applied　in　the　random　dimension　to　derive　the　equations　in　the　weak　form.　The　resulting　set　of　deterministic　equations　is　then　solved　with　standard　methods　to　obtain　the　mean　solution.　Ill　this　article,　the　main　method　employs　the　Hermite　polynomial　as　the　basis　in　random　space.　Numerical　examples　are　given　and　the　error　analysis　is　demonstrated　for　a　model　problem.　(C)　2008　Wiley　Periodicals,　Inc.　Numer　Methods　Partial　Differential　Eq　26:　14-23,　2010