Modeling uncertainty in flow simulations via generalized polynomial chaos

We present a new algorithm to model the input uncertainty and its propagation in incompressible flow simulations. The stochastic input is represented spectrally by employing orthogonal polynomial functionals from the Askey scheme as trial basis to represent the random space. A standard Galerkin projection is applied in the random dimension to obtain the equations in the weak form. The resulting system of deterministic equations is then solved with standard methods to obtain the solution for each random mode. This approach can be considered as a generalization of the original polynomial chaos expansion, first introduced by Wiener [Am. J. Math. 60 (1938) 897]. The original method employs the Hermite polynomials (one of the 13 members of the Askey scheme) as the basis in random space. The algorithm is applied to micro-channel flows with random wall boundary conditions, and to external flows with random freestream. Efficiency and convergence are studied by comparing with exact solutions as well as numerical solutions obtained by Monte Carlo simulations. It is shown that the generalized polynomial chaos method promises a substantial speed-up compared with the Monte Carlo method. The utilization of different type orthogonal polynomials from the Askey scheme also provides a more efficient way to represent general non-Gaussian processes compared with the original Wiener–Hermite expansions.