Let $\Omega$ an open bounded in $\mathbb{R}^n$ and $I$ un interval in $\mathbb{R}$, let $t_0 \in I$, $f\in H^1_0(\Omega)$, $g\in L^2(\Omega)$. Let $(\lambda_n)_{n\in \mathbb{N}}$ the eigen values od the Laplacian operator $\Delta$, and $(e_n)_n$ the vector associates to the eigen values. We denote $$u=\sum_{n=1}^{+\infty} [(f,e_n) \cos[(t-t_0) \sqrt{\lambda_n}] + \dfrac{(g,e_n)}{\sqrt{\lambda_n}}\sin[(t-t_0)\sqrt{\lambda_n}]]e_n$$
How we can prouve that $u \in C^0(I,H^1_0(\Omega))$ and $\dfrac{\partial u}{\partial t} \in C^0(I,L^2(\Omega))$?
II
I undestand the answer of this question. I have an seconde question please.
How we prouve that $$\Box u = 0 \quad\mbox{in} \quad \mathcal{D}'(I^{0} \times \Omega)$$ where $\Box$ is wave operator, and $I^{0}$ is the interior od $u$.
Thank's for the help beaucause i lost.
I'll assume $\Omega$ has a smooth boundary (or something similar.) The Laplacian has a natural domain $\mathcal{D}(\Delta)=H^{2}(\Omega)\cap H^{1}_{0}(\Omega)$ because $$ -\Delta : \mathcal{D}(\Delta) \subset L^{2}(\Omega)\rightarrow L^{2}(\Omega) $$ is selfadjoint and non-negative on its domain. The domain of the positive square root of $-\Delta$ is $H_{0}^{1}(\Omega)$. Furthermore, $-\Delta$ has a complete orthonormal set of eigenfuntions $\{ e_{n}\}_{n=0}^{\infty}$ with corresponding eigenvalues $\{ \lambda_{n} \}_{n=0}^{\infty} \subset (0,\infty)$. By the spectral theorem, $$ f \in \mathcal{D}(-\Delta) \iff \sum_{n=0}^{\infty}|\lambda_{n}(f,e_{n})|^{2} < \infty,\;\;\; f \in \mathcal{D}(\sqrt{-\Delta}) \iff \sum_{n=0}^{\infty}|\sqrt{\lambda_{n}}(f,e_{n})|^{2} < \infty. $$ Using Friedrichs extensions, one finds that $\mathcal{D}(\sqrt{-\Delta})=H_{0}^{1}(\Omega)$. So the above characterization of $\mathcal{D}(\sqrt{-\Delta})$ becomes a characterization of $H_{0}^{1}$. Because $\sqrt{-\Delta}$ is selfadjoint, then $f(t)=e^{it\sqrt{-\Delta}}f_{0}$ is a continuous function from $[0,\infty)$ to $L^{2}(\Omega)$ which is given by Functional Calculus: $$ f(t) = e^{it\sqrt{-\Delta}}f_{0}=\sum_{n=0}^{\infty} e^{i\sqrt{\lambda_{n}}t}(f_{0},e_{n})e_{n}. $$ If $f_{0} \in \mathcal{D}(\sqrt{-\Delta})=H_{0}^{1}(\Omega)$, then $f$ is a continuous function from $[0,\infty)$ to $H_{0}^{1}(\Omega)$, and $f'(t)=\sqrt{-\Delta}f(t)=e^{it\sqrt{-\Delta}}(\sqrt{-\Delta}f_{0})$ is a continuous function from $[0,\infty)$ to $L^{2}(\Omega)$. Your situation is a mild variation of this.
Your function $u$: (Discussion added after your questions.) So you have defined $u=v+w$ where $$ v(t) = \sum_{n=0}^{\infty}(f,e_{n})\cos((t-t_{0})\sqrt{\lambda_{n}})e_{n}, \;\;\; f \in H_{0}^{1}(\Omega)\\ w(t) = \sum_{n=0}^{\infty}\frac{(g,e_{n})}{\sqrt{\lambda_{n}}}\sin((t-t_{0})\sqrt{\lambda_{n}})e_{n},\;\;\; g \in L^{2}(\Omega). $$ Both $v$ and $u$ are in $L^{2}(\Omega)$ for all $t\in\mathbb{R}$ because the Fourier coefficents with respect to the complete orthonormal set $\{ e_{n}\}_{n=0}^{\infty}$ in each series are square integrable: $$ \|v(t)\|^{2}_{L^{2}}=\sum_{n=0}^{\infty}|(f,e_{n})|^{2}\cos^{2}((t-t_{0})\sqrt{\lambda_{n}}) \le \sum_{n} |(f,e_{n})|^{2} = \|f\|^{2},\\ \|w(t)\|_{L^{2}}^{2}=\sum_{n=0}^{\infty}\frac{|(g,e_{n})|^{2}}{\lambda_{n}}\sin^{2}((t-t_{0})\sqrt{\lambda_{n}}) \le \sum_{n=0}^{\infty}\frac{1}{\lambda_{1}}|(g,e_{n})|^{2} = \frac{1}{\lambda_{1}}\|g\|^{2}. $$ Because $f \in H^{1}_{0}$, then the equivalent condition $\sum_{n=0}^{\infty}|\sqrt{\lambda_{n}}(f,e_{n})|^{2} < \infty$ holds. This forces $v(t) \in H^{1}_{0}$ for all $t$ because $$ \sum_{n=0}^{\infty}|\sqrt{\lambda_{n}}(v(t),e_{n})|^{2}= \sum_{n=0}^{\infty}|\sqrt{\lambda_{n}}(f,e_{n})|^{2}\sin^{2}((t-t_{0})\sqrt{\lambda_{n}})\le \sum_{n=0}^{\infty}|\sqrt{\lambda_{n}}(f,e_{n})|^{2} < \infty. $$ Similarly, $w(t) \in H^{1}_{0}$ for all $t$ because $g\in L^{2}$ and you're dividing the Fourier coefficients of $g$ by $\sqrt{\lambda_{n}}$, which puts the result in $H^{1}_{0}$. So everything is properly defined and $u(t) \in H^{1}_{0}$ for all $t$.
As I pointed out before, $h \in H^{1}_{0}$ iff $h\in\mathcal{D}(\sqrt{-\Delta})$, and, in that case, $$ \|h\|^{2}_{H^{1}_{0}}=\|h\|^{2}_{L^{2}}+\|\nabla h\|^{2}_{L^{2}} = \|h\|^{2}_{L^{2}}+\|\sqrt{-\Delta}h\|^{2}_{L^{2}}=\sum_{n=0}^{\infty}(1+\lambda_{n})|(f,e_{n})|^{2}. $$ In fact, for $h \in L^{2}$, the sum on the right is finite iff $h \in H^{1}_{0}$. To shorten the expressions, let $s_{n}(t)=\sin((t-t_{0})\sqrt{\lambda_{n}})$ and $c_{n}(t)=\cos((t-t_{0})\sqrt{\lambda_{n}})$. What you've got is $$ \begin{align} \|u(t)-u(t')\|_{H^{1}_{0}} & \le\left[\sum_{n=0}^{\infty}(1+\lambda_{n})|(f,e_{n})|^{2}|c_{n}(t)-c_{n}(t')|^{2}\right]^{1/2} \\ & +\left[\sum_{n=0}^{\infty}\frac{1+\lambda_{n}}{\lambda_{n}}|(g,e_{n})|^{2}|s_{n}(t)-s_{n}(t')|^{2}\right]^{1/2}. \end{align} $$ The sums are absolutely and uniformly convergent with respect to $t$, $t'$ because $f \in H^{1}_{0}$ and $g\in L^{2}$. It follows that the above is continuous in $t,t'$. So $u : \mathbb{R}\rightarrow H^{1}_{0}$ is a continuous vector function. And when you take a derivative of $u$, you end up with a continuous function $u' : \mathbb{R}\rightarrow L^{2}$. The reduced regularlity is because of the added $\sqrt{\lambda_{n}}$ which appears in the derived sums. I'll let you work out the details.
Your New Modifications: The new mention of $\Box$ requires that you be able to take two time derivatives. This can be done if $f\in \mathcal{D}(-\Delta)$ and $g \in \mathcal{D}(\sqrt{-\Delta})=H_{0}^{1}(\Omega)$. If $f$, $g$ are as stated, then taking two time derivatives gives the same result as applying $-\Delta$ (both multiply the n-th Fourier coefficient by $\lambda_{n}$.) The expression $\Box$ is now interpreted in terms of a derivative of a vector function in $t$ and in terms of an operator $-\Delta$, which is more subtle than classical derivatives. So $\frac{d^{2}}{dt^{2}}u(t) = \Delta u(t)$ in this operator language, where $u : [t_{0},\infty)\rightarrow L^{2}(\Omega)$ is a vector function whose value lies in the domain of $-\Delta$ for all $t > 0$ so that $-\Delta u(t)$ makes sense, provided $f \in \mathcal{D}(-\Delta)$ and $g\in H_{0}^{1}$.