Let me start with an example. If we base our calculations on the Newton's second law without any further mathematical treatment, then our equation describes equilibrium of forces, i.e. it is of the kind "force = force (e.g. as the product of mass and acceleration)". Similarly it could be momenta, angular momenta, electric charges, voltages, various quantities like "speed change per kilogram and tesla" etc. etc. etc…
How is that with the wave equation? It's crucial for interpretations of source terms. So, could you describe like that inhomogenous string equation:
$$
\frac{\partial^2y}{\partial x^2}-\frac{1}{c_0^2}\frac{\partial^2y}{\partial t^2}=q
$$
Wave equation for the acoustic pressure:
$$
\nabla^2 p-\frac{1}{c_0^2}\frac{\partial^2p}{\partial t^2}=q
$$
and for the acoustic velocity:
$$
\nabla^2 \vec{v}-\frac{1}{c_0^2}\frac{\partial^2 \vec{v}}{\partial t^2}=q
$$
where $q$ are source terms with required dimension and units. What is the physical meaning of those terms?
Note: If that would be considered helpful for the site, let's broaden the list of wave equations.
Edit to clarify: I am not looking for the steps in mathematical solution or the general principle of pressence a pulsating body in a fluid (in the case of acoustics). Question is how exactly represent $q$ that pulsation? What is the intuition and meaning of the term with dimension of $time^{-1} \times distance^{-1}$ (that's for acoustical velocity)? Etc.
Best Answer
To solve a non homogeneous wave equation you need to use Green's function. The Green function of the inhomogeneous wave equation is defined as $$ (\nabla^2-\frac{1}{c^2}\partial_t^2)G(x,x';t,t') =\delta(x-x')\delta(t-t') $$ It means that $G$ is the solution of the wave equation for a Dirac source localized at the position $r'$ and at a time $t'$. Why is this Green function important ? Because once you have it you can have any solution $\Phi(x,t)$ to the equation $$ (\nabla^2-\frac{1}{c^2}\partial_t^2)\Phi(x,t)=\rho(x,t) $$ by computing $$ \Phi=\int G(x,x';t,t')\rho(x',t')dx'dt' \space \space\space\space\space(1) $$
because
$$ (\nabla^2-\frac{1}{c^2}\partial_t^2)\Phi=\int (\nabla^2-\frac{1}{c^2}\partial_t^2)G(x,x';t,t')\rho(x',t')dx'dt' $$ using the definition of the Green function: $$ =\int \delta(x-x')\delta(t-t')\rho(x',t')dx'dt'= \rho(x,t) $$ In equation (1) lies the answer to your question: To obtain your solution you have to sum up the influence of every infinitesimal part of your source term $\rho$. This influence is given by the Green function, which is exactly the solution of the same problem for an infinitesimal pointlike source.
Its a little bit simpler in frequency domain, or even in static (where you have Poisson or Laplace equation instead of the wave equation) but the idea is always the same: obtain a Green function for a point-like (Dirac source) and then obtain the solution by integration.
A good example to physically illustrate this idea is the Huygens principle: every part of the phase front of a wave is a little source of a spherical wave. If you want the solution of a problem where your source ($\rho$) is a circular hole cut in an opaque screen you'll have to sum up little spherical wave (Green function) over all the hole surface. That's how diffraction problem are solved.