Because typical scalar wave equation has form,
$$ \tag 1 \ddot u = v^2 \nabla ^2 u ,$$
to which plane-wave eigenmodes (with a constant frequency) temporal part has form
$$ \tag 2 u(t) = C e^{-i\omega t} ,$$
where $C$ being some complex amplitude carrying wave phase information.
According to Euler's formula,
$$ \tag 3 e^{ix}=\cos x+i\sin x $$
real part of (2) can be rewritten as,
$$ \tag 4 Re[u(t)] = A \cos(\omega t) $$
Hence, the periodic function $\cos()$ in the scalar field temporal fluctuation. As an additional bonus, $\sin/\cos$ works very good in Fourier transform, i.e. signal can be split into different harmonics, where multiple sine/cosine waves interferes with each other producing final wave. Not to mention that for this analysis wave defining function $u(t)$ must be differentiable and continuous, which applies to periodic trigonometric functions, like sine.
Sure, general form of scalar wave solution does not have to be bound only by raw trigonometric functions, as it is in the form,
$$ \tag 5 u(x,t)=F(x-vt)+G(x+vt). $$
And so you can also have Sawtooth wave signal, like

but since $F(),G()$ are not differentiable here, you need to fallback to numerical methods like Discrete Fourier transform for component analysis, which may or may not be cumbersome. So once again- periodic differentiable functions like $\sin,\cos$ makes analytical wave analysis approachable.