Polynomial bounds for the solutions of parametric transmission problems on smooth, bounded domains

We consider a family (P_ω)ω ∈ Ω of elliptic second order differential operators on a domain U₀ ⊂ ᵐ whose coefficients depend on the space variable x ∈ U₀ and on ω ∈ Ω, a parameter space. We allow the coefficients aᵢⱼ of P_ω to have jumps over a fixed interface Γ ⊂ U₀ (independent of ω ∈ Ω). We obtain estimates on the norm of the solution u_ω to the equation P_ω u_ω = f with transmission and mixed boundary conditions that are polynomial in the norms of the coefficients. In particular, we show that, if f and the coefficients aᵢⱼ are smooth enough and follow a log-normal-type distribution, then the map Ω ω → \|u_ω\|_Hᵏ⁺¹(U₀) is in Lᵖ(Ω), for all 1 ≤ p < ∞. The same is true for the norms of the inverses of the resulting operators. We also obtain similar integrability results for the parametric Finite Element approximations of the solution. We expect our estimates to be useful in Uncertainty Quantification.