Solving random diffusion models with nonlinear perturbations by the Wiener–Hermite expansion method

This paper deals with the construction of approximate series solutions of random nonlinear diffusion equations where nonlinearity is considered by means of a frank small parameter and uncertainty is introduced through white noise in the forcing term. For the simpler but important case in which the diffusion coefficient is time independent, we provide a Gaussian approximation of the solution stochastic process by taking advantage of the Wiener–Hermite expansion together with the perturbation method. In addition, approximations of the main statistical functions associated with a solution, such as the mean and variance, are computed. Numerical values of these functions are compared with respect to those obtained by applying the Runge–Kutta second-order stochastic scheme as an illustrative example.