  In this paper, we study the inverse source problem for the biharmonic wave
equation. Mathematically, we characterize the radiating sources and
non-radiating sources at a fixed wavenumber. We show that a general source can
be decomposed into a radiating source and a non-radiating source. The radiating
source can be uniquely determined by Dirichlet boundary measurements at a fixed
wavenumber. Moreover, we derive a Lipschitz stability estimate for determining
the radiating source. On the other hand, the non-radiating source does not
produce any scattered fields outside the support of the source function.
Numerically, we propose a novel source reconstruction method based on Fourier
series expansion by multi-wavenumber boundary measurements. Numerical
experiments are presented to verify the accuracy and efficiency of the proposed
method.
