Let $u(x,y,t)$ be the solution of wave equation $u_{tt}=u_{xx}+u_{yy}, 0\lt x\lt 1, 0\lt y\lt 1, t\ge 0,$ $u(x,y,0)=(xx^2)(yy^2), u_t(x,y,0)=0$ and $ u(x,y,t)=0$ on the boundary of the square. Then $$u(1/2,1/2,t)=f(t)=\frac{64}{\pi^6}\sum_{n=0}^\infty\sum_{m=0}^\infty \frac{(1)^{m+n}}{(2m+1)^3(2n+1)^3}\cos\left(\pi t \sqrt{(2m+1)^2+(2n+1)^2}\right)$$ It is unlikely that the function has a closed form expression. My question is: Is there an analytical way to find $\sup_{t\ge 0} f(t)$ and $\limsup_{t\to\infty} f(t)$?

An easy to find upper bound for the expression without the constant term is $$\sum_{n=0}^{\infty}\sum_{m=0}^{\infty}\frac{1}{(2m+1)^3(2n+1)^3}\approx 1.1063.$$ For each pair of a pair of odd numbers $(n_1,m_1)$ and $(n_2,m_2)$ with $n_1^2+m_1^2=n_2^2+m_2^2$ and $n_1+m_1 \not \equiv n_2+m_2\ (\mod 4)$ (e.g. (7,1) and (5,5)) we can improve this bound a little bit by $$2 \min\left(\frac{1}{n_1^3m_1^3},\frac{1}{n_2^3m_2^3}\right)$$ as $$(1)^{\frac{n_1+m_12}{2}}\cos(\pi t \sqrt{n_1^2+m_1^2})=(1)^{\frac{n_2+m_22}{2}} \cos(\pi t \sqrt{n_2^2+m_2^2}).$$ Of course this is far from being a solution. 

