Second order directional shape derivatives of integrals on submanifolds

We compute first and second order shape sensitivities of integrals on smooth submanifolds using a variant of shape differentiation. The result is a quadratic form in terms of one perturbation vector field that yields a second order quadratic model of the perturbed functional. We discuss the structure of this derivative, derive domain expressions and Hadamard forms in a general geometric framework, and give a detailed geometric interpretation of the arising terms.