Let $\Omega$ be domain of $\mathbb{R}^n$ and $\Phi : \Omega \to \Phi(\Omega)$ a deformation. Consider the Stokes equations written in the deformed configuration
\begin{align} - 2\mu \operatorname{div}(D(u)) + \nabla p &= f, \quad \text{in } \Phi(\Omega) \\ \operatorname{div}(u) &= 0, \quad \text{in } \Phi(\Omega) \\ \end{align}
where $u$ is the velocity of the fluid, $p$ the pressure, $\mu>0$ is the constant viscosity, $f$ is an external force and $D$ is the operator defined by
$$ D(u) = \frac{1}{2} (\nabla u + \nabla u^T). $$
How can the Stokes equations be written in the domain $\Omega$ using a change of variable ?