Adding random noise to the solution of 2D stokes problem

The easiest solution is probably to use an Expression or a UserExpression to directly interpolate a random variable into an appropriate FunctionSpace. See this thread for examples: