This talk is concerned with the numerical solution of partial differential equations on random domains. On the one hand, we employ multilevel quadrature techniques for the stochastic variable in order to speed-up the computations. On the other hand, we discuss the use of Euler and Lagrange coordinates for the physical variable. This choice is reflected by the numerical solver for particular realizations of the underlying boundary value problem. In particular, we consider boundary element methods, finite element methods, and the coupling of finite element and boundary element methods.