where u is the displacement vector, k and l are Lame’s constants
for the material, q is the density, and r2 is the three-dimensional
Laplace operator.
By introducing Helmholtz decomposition, the vector u is
decomposed in terms of a scalar potential function and a vector
potential function as