  The present paper deals with the interior solid-fluid interaction problem in
harmonic regime with randomly perturbed boundaries. Analysis of the shape
derivative and shape Hessian of vector- and tensor-valued functions is
provided. Moments of the random solutions are approximated by those of the
shape derivative and shape Hessian, and the approximations are of third order
accuracy in terms of the size of the boundary perturbation. Our theoretical
results are supported by an analytical example on a square domain.
