Shape optimization problems with random coefficients via the penalty method

For shape optimization problems, governed by elliptic equations with Dirichlet boundary condition and random coefficients, we utilize a penalization technique to get the approximate problem. We consider that uncertainties exists in the diffusion coefficients and minimize objective functions in mean value form. Finite element method, Monte Carlo method and accelerated version of the gradient descent method are applied to solve the corresponding discretized problem. The convergence analysis and numerical results are included.