\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2012 (2012), No. 116, pp. 1--9.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2012 Texas State University - San Marcos.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2012/116\hfil Solution to nonlocal problems]
{Solution to nonlocal problems of pseudohyperbolic equations}

\author[L. S. Pulkina \hfil EJDE-2012/116\hfilneg]
{Ludmila S. Pulkina}  % in alphabetical order

\address{Ludmila S. Pulkina \newline
 Samara State University, Samara, Russia}
\email{louise@samdiff.ru}

\thanks{Submitted January 25, 2014. Published April 21, 2014.}
\subjclass[2000]{35D05, 35L20, 35M99}
\keywords{Nonlocal; pseudohyperbolic equation; integral condition}

\begin{abstract}
 In this article considers a nonlocal problem with integral condition 
 for a fourth-order pseudohyperbolic equation. Existence and
 uniqueness of a generalized solution are proved.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{remark}[theorem]{Remark}
\newtheorem{definition}[theorem]{Definition}
\allowdisplaybreaks

\section{Introduction}

Currently, there is considerable interest in nonlocal problems for
evolution equations. One reason for this lies in the fact that
various phenomena of modern natural science can be described most conveniently
in terms of nonlocal problems. Problems with nonlocal integral conditions form
an important class of nonlocal problems. Recently, nonlocal boundary value
problems with integral conditions have been actively studied.
However, the majority of the works deals with second-order equations.
The initial works devoted to nonlocal problems for second-order partial
differential equations  with integral conditions go back to
Cannon \cite{cannon} and  Kamynin \cite{kamyn}.
 Note here some recent works:
\cite{A-A-G,A-A,bouzian,gia,K-P,kozh_S,Pu_3,Pu_04,Pu_08,Pu_2}.
 See also references therein.

Pseudohyperbolic equations form important and interesting subclass of Sobolev
type equations. Such equations may describe nonstationary waves in stratified
 and rotating liquid \cite{Korp}. The starting point in studying of Sobolev
type  equations is \cite{sobolev}. Now there are a lot of works devoted to
initial and boundary value problems for Sobolev type equations
(see \cite{sv} and references therein). One of recent works
dealing with some problems for pseudohyperbolic equations is \cite{Korp}.
In these work the author studies qualitative characteristics of solutions
to initial-boundary value problems. On the other hand, various physical
problems demand nonlocal conditions
\cite{can_est,cannon-lin,can_lin,kamV,Kozh_In,samar}.

Motivated by this, we consider a nonlocal problem with integral condition for
a pseudohyperbolic equation.

\section{Results}

Let $\Omega$ be a bounded domain in $R^n$ with smooth boundary
$\partial \Omega$, $Q_T=\Omega\times (0,T)$,
$S_T=\partial \Omega \times (0,T)$.  Consider an equation
\begin{equation} \label{eq1}
Lu\equiv\frac{\partial ^2}{\partial t^2}(u-\Delta u)
-(a_{ij}(x,t)u_{x_i})_{x_j}+c(x,t)u=f(x,t)
\end{equation}
and set a problem: Find a function $u(x,t)$  that is a solution of \eqref{eq1}
in $Q_T$,
satisfies initial data
\begin{equation}
 \label{eq2}
u(x,0)=0, \quad u_t(x,0)=0
\end{equation}
and nonlocal condition
\begin{equation}
 \label{eq3}
\Big(\frac{\partial ^2}{\partial t^2}\frac{\partial u}{\partial \nu}
+\frac{\partial u}{\partial N} +\int_{\Omega}K(x,y,t)u(y,t)dy\Big)\big|_{S_T}=0.
\end{equation}

We suppose throughout that repeated indices imply summation from $1$ to $n$, $\nu(x)=(\nu_1,\dots ,\nu_n)$
is outward unit normal vector to $\partial \Omega$ at the current point,
$\frac{\partial u}{\partial N}=a_{ij}u_{x_i}\cos (\nu,x_j)$,
$a_{ij}=a_{ji}$, $ K(x,y,t)$ is a given function.
We also use the following notation:
$$
u_x=(u_{x_1},\dots u_{x_n}), \quad
u_x^2=\sum_{i=1}^nu^2_{x_i}, \quad
u_xv_x=\sum_{i=1}^nu_{x_i}v_{x_i}.
$$
As the condition \eqref{eq3} does not look evident we give some explanations 
in an appendix at the end of the paper.

Let $W_2^1(Q_T)$ be the usual Sobolev space. We shall define
\begin{gather*}
W(Q_T)=\{u: u\in W_2^1(Q_T), \,  u_{xt}\in L_2(Q_T)\},\\
\hat{W}(Q_T)=\{v: v\in W(Q_T), \, v(x,T)=0\}.
\end{gather*}
First we give a definition of a generalized solution to the problem using 
the standard method \cite[p. 92]{lad}.
To this end we multiply \eqref{eq1} by $v\in \hat{W}(Q_T)$ and integrate over 
$Q_T$. It follows from \eqref{eq2}, \eqref{eq3} and an integration by parts that
\begin{equation} \label{eq4}
\begin{aligned}
&\int_0^T\int_{\Omega}(-u_tv_t-u_{xt}v_{xt}+a_{ij}u_{x_i}v_{x_j}+cuv)\,dx\,dt\\
&+\int_0^T\int_{\partial \Omega}v(0,t)\int_{\Omega}K(x,y,t)u(x,t)\,dy\,ds\,dt\\
&=\int_0^T\int_{\Omega}fv\,dx\,dt.
\end{aligned}
\end{equation}

\begin{definition} \label{def1} \rm
A function $u\in W(Q_T)$ is said to be a generalized solution
to the problem  \eqref{eq1}--\eqref{eq3} if $u(x,0)=0$ and for every 
$v\in \hat{W}(Q_T)$ the identity \eqref{eq4} holds.
\end{definition}

\begin{theorem} \label{thm1}
If the function $K(x,y,t)$ is continuous in $\bar \Omega \times \bar{Q}_T $, 
\begin{gather*}
f\in L_2(Q_T), \quad c\in C(\bar{Q}_T), \quad 
a_{ij}, \frac{\partial a_{ij}}{\partial t}\in C(\bar{Q}_T),\\
 \forall (x,t)\in \bar{Q}_T, \; \gamma \xi^2\leq a_{ij}(x,t)\xi_i\xi_j\leq \mu\xi^2,
 \; \gamma>0,
\end{gather*}
then there exists a unique generalized solution to the problem 
\eqref{eq1}--\eqref{eq3}.
\end{theorem}

\begin{proof}
First we prove the uniqueness. 
To this end we obtain a number of inequalities and
then use Gronwall's lemma. We prove the existence part in several steps. 
First, we construct approximations of the generalized solution by 
the Faedo-Galerkin method. Then we obtain
a priori estimates to guarantee convergence of approximations. 
Finally, we show that the limit of approximations is the required solution.

\noindent\textbf{Uniqueness.} 
Suppose that $u_1$ and  $u_2$ are two different solutions to \eqref{eq1}--\eqref{eq3}.
Then $u=u_1-u_2$ satisfies $u(x,0)=0$ and the identity
\begin{equation}\label{eq5}
\begin{aligned}
&\int_0^T\int_{\Omega}(-u_tv_t-u_{xt}v_{xt}+a_{ij}u_{x_i}v_{x_j}+cuv)\,dx\,dt\\
&+\int_0^T\int_{\partial \Omega}v(x,t)\int_{\Omega}K(x,y,t)u(y,t)\,dy\,ds\,dt=0
\end{aligned}
\end{equation}
holds for every $v\in \hat{W}(Q_T)$. For an arbitrary $\tau \in [0,T]$, take $v$ as
\begin{equation}  \label{eq6}
v(x,t)=\begin{cases}
\int_{\tau}^tu(x,\eta)d\eta, & 0\leq t\leq \tau,\\
0, & \tau \leq t\leq T.
\end{cases}
\end{equation}
It is easy to see that $v\in \hat{W}(Q_T)$, as well $v_t=u$, $v_{xt}=u_x$ in 
$Q_{\tau}=\Omega \times (0,\tau)$.

Substitute $v(x,t)$ from \eqref{eq6} in \eqref{eq5} and express $u$ in terms 
of $v$ and its derivatives.
As a result we obtain the equality
\begin{align*}
&\int_0^{\tau}\int_{\Omega}(-v_{tt}v_t+a_{ij}v_{x_it}v_{x_j}-v_{ttx}v_{xt}
 +cv_tv)\,dx\,dt\\
&+\int_0^{\tau}\int_{\partial \Omega}v(x,t)\int_{\Omega}K(x,y,t)v_t(y,t)\,dy\,ds\,dt=0.
\end{align*}
After integrating by parts first three terms, we obtain
\begin{equation}  \label{eq7}
\begin{aligned}
&\frac{1}{2}\int_{\Omega}[v_t^2(x,\tau)+a_{ij}(x,0)v_{x_i}(x,0)v_{x_j}(x,0)
 +v^2_{xt}(x,\tau)]dx \\
&=\int_0^{\tau}\int_{\Omega}c(x,t)v(x,t)v_t(x,t)\,dx\,dt
 -\frac{1}{2}\int_0^{\tau}\int_{\Omega}
 \frac{\partial a_{ij}}{\partial t}v_{x_i}v_{x_j}\,dx\,dt \\
&\quad +\int_0^{\tau}\int_{\partial \Omega}v(x,t)
\int_{\Omega}K(x,y,t)v_t(y,t)\,dy\,ds\,dt.
\end{aligned}
\end{equation}
Our next aim is to derive an estimate of a right-hand side of \eqref{eq7}.
Taking into account hypotheses of the theorem we can see that there exists 
positive number $c_0$ such that 
$$
\max_{\bar{Q}_T}\{\big|\frac{\partial a_{ij}}{\partial t}\big|, \ |c|\}\leq c_0.
$$
Let 
$$
k=\max_{\bar{Q}_T}\int_{\Omega}K^2(x,y,t)dy, \quad \omega=\int_{\partial \Omega}ds.
$$
Applying the Cauchy inequality we obtain
\begin{gather*}
\big|\int_0^{\tau}\int_{\Omega}cv_tv\,dx\,dt\big|
\leq \frac{c_0}{2}\int_0^{\tau}\int_{\Omega}(v_t^2+v^2)\,dx\,dt;
\\
\frac{1}{2}\big|\int_0^{\tau}\int_{\Omega}\frac{\partial a_{ij}}{\partial t}v_{x_i}v_{x_j}\,dx\,dt\big|
\leq c_0\int_0^{\tau}\int_{\Omega}v_x^2\,dx\,dt;
\\
\big|\int_0^{\tau}\int_{\partial \Omega}v\int_{\Omega}Kv_t\,dy\,ds\,dt\big|\leq
\frac{1}{2}\int_0^{\tau}\int_{\partial \Omega}v^2\,ds\,dt+\frac{\omega k}{2}
\int_0^{\tau}\int_{\Omega}v_t^2\,dy\,dt.
\end{gather*}
As by hypotheses $\partial \Omega$ is smooth then (see \cite[p. 77]{lad})
$$
\int_{\partial \Omega}v^2ds\leq c_1\int_{\Omega}(v_x^2+v^2)dx
$$
and we obtain the inequality
\begin{equation}  \label{eq8}
\begin{aligned}
&\int_{\Omega}[v_t^2(x,\tau)+a_{ij}(x,0)v_{x_i}(x,0)v_{x_j}(x,0)
 +v^2_{xt}(x,\tau)]dx \\
&\leq c_2\int_0^{\tau}\int_{\Omega}(v_t^2+v_x^2+v^2)\,dx\,dt,
\end{aligned}
\end{equation}
where $c_2$ depends only on $c_0$, $c_1$, $k$, and $\omega$.

Introduce now the functions $w_i(x,t)=\int_0^tu_{x_i}(x,\eta)d\eta$.
By  \eqref{eq6},
$$
v_{x_i}(x,t)=w_i(x,t)-w_i(x,\tau), \quad
 v_{x_i}(x,0)=-w_i(x,\tau).
$$
Furthermore, for a.e. $x\in \Omega$,
$$
\int_0^{\tau}v^2dt=\int_0^{\tau}\Big(\int_{\tau}^tu(x,\eta)d\eta\Big)^2dt
\leq \tau^2\int_0^{\tau}u^2dt.
$$
Thus,  from \eqref{eq8}, it follows that
\begin{align*}
&\int_{\Omega}[u^2(x,\tau)+a_{ij}(x,0)w_i(x,\tau)w_j(x,\tau)+u_x^2(x,\tau)]dx\\
&\leq2c_2\int_0^{\tau}\int_{\Omega}[(1+\tau^2/2)u^2
 +\sum_{i=1}^nw_i^2]\,dx\,dt+2c_2\tau \int_{\Omega}
\sum_{i=1}^nw_i^2(x,\tau)dx.
\end{align*}

Note that $a_{ij}(x,0)w_i(x,\tau)w_j(x,\tau)\geq \gamma w^2$.
As $\tau$ is arbitrary we choose it in such a way that an inequality
 $\gamma-2c_2\tau>0$ holds. 
Let  $\gamma-2c_2\tau\geq \gamma/2$. Then for every 
$\tau \in [0,\frac{\gamma}{4c_2}]$,
$$
\int_{\Omega}[u^2(x,\tau)+\sum_{i=1}^nw_i^2(x,\tau)+u_x^2]dx\leq
 c_3\int_0^{\tau}\int_{\Omega}(u^2(x,t)+\sum_{i=1}^nw_i^2(x,t))\,dx\,dt,
$$
with $c_3=c_2\max\{1+\tau^2, 2\}/\min\{1,\gamma/2\}$, and in particular,
$$
\int_{\Omega}[u^2(x,\tau)+\sum_{i=1}^nw_i^2(x,\tau)]dx\leq
 c_3\int_0^{\tau}\int_{\Omega}(u^2(x,t)+\sum_{i=1}^nw_i^2(x,t))\,dx\,dt.
$$
Now by Gronwall's lemma we conclude that, for $\tau \in [0,\frac{\gamma}{4c_2}]$,
$$
\int_{\Omega}(u^2(x,\tau)+\sum_{i=1}^nw_i^2(x,\tau))dx\leq 0.
$$
It follows immediately  that $u(x,\tau)=0$ for $\tau \in [0,\frac{\gamma}{4c_2}]$.

Following \cite{lad} we repeat these arguments for
 $\tau \in [\frac{\gamma}{4c_2}, \frac{\gamma}{2c_2}]$
and then continue this procedure. It follows that  $u(x,\tau)=0$ for all
 $\tau \in[0,T]$.
It implies that there exists at most one solution to  \eqref{eq1}--\eqref{eq3}.
\smallskip

\noindent\textbf{Existence.}
Let $w_k(x)\in C^2(\bar \Omega)$ be a basis in $W_2^1(\Omega)$.
We define the approximations
\begin{equation}
 \label{eq9}
u^m(x,t)=\sum_{k=1}^mc_k(t)w_k(x),
\end{equation}
where $c_k(t)$ are solutions to the Cauchy problem
\begin{gather}  \label{eq10}
\begin{aligned}
&\int_{\Omega}(u^m_{tt}w_p+a_{ij}u^m_{x_i}w_{px_j}+u^m_{x_itt}w_{px_j}+cu^mw_p)dx\\
&+\int_{\partial \Omega}w_p(x)\int_{\Omega}K(x,y,t)u^m\,dy\,ds\\
&=\int_{\Omega}fw_p\,dx,
\end{aligned}
\\  \label{eq11}
c_k(0)=0, \quad c'_k(0)=0.
\end{gather}
We write the Cauchy problem \eqref{eq10}--\eqref{eq11}
such that:
\begin{equation} \label{eq12}
\begin{gathered}
\sum_{k=1}^mc_k''(t)A_{kp}+\sum_{k=1}^mc_k(t)B_{kp}(t)=f_p(t),\\
A_{kp}=(w_k, w_p)_{W_2^1(\Omega)},
\end{gathered}
\end{equation}
where
\begin{gather*}
\begin{aligned}
B_{kp}(t)&=\int_{\Omega}[a_{ij}(x,t)w_{kx_i}w_{px_j}+c(x,t)w_kw_p]dx\\
&\quad +\int_{\partial \Omega}w_p(x)\int_{\Omega}K(x,y,t)w_k(y)\,dy\,ds,
\end{aligned}
\\
f_p(t)=\int_{\Omega}f(x,t)w_p(x)dx.
\end{gather*}

Note that the matrix $\|(w_k, w_j)_{W_2^1(\Omega)}\|$ is Gramian matrix as the
functions $w_k$ are linearly independent, hence the system \eqref{eq12} is normal.
Under the hypothesis of the theorem  coefficients $A_{kp}$, $B_{kp}$ are bounded 
and $f_j\in L_1(0,T)$.
Thus the Cauchy problem has a unique solution $c_k\in W_2^2(0,T)$ for every $m$ and
all approximations \eqref{eq9} are defined.

Next, we need a priori estimates to pass to the limit as $m\to \infty$.

Multiplying \eqref{eq10} by $c'_p(t)$, summing from $p=1$ to $p=m$
and integrating with respect to $t$ from $0$ to $\tau$, we obtain
\begin{equation}\label{eq13}
\begin{aligned}
&\int_0^{\tau}\int_{\Omega}(u^m_{tt}u^m_t+a_{ij}u^m_{x_i}u^m_{x_jt}
 +u^m_{x_itt}u^m_{x_jt} +cu^mu_t^m)\,dx\,dt\\
&+ \int_0^{\tau}\int_{\Omega}\frac{\partial a_{ij}}{\partial x_i}
 u^m_{x_i}u^m_t\,dx\,dt
 +\int_0^{\tau}\int_{\partial \Omega}u_t^m\int_{\Omega}K(x,y,t)u^m\,dy\,ds\,dt\\
&=\int_0^{\tau}\int_{\Omega}f(x,t)u_t^m(x,t)\,dx\,dt.
\end{aligned}
\end{equation}
Integrating by parts on the first term of the left-hand side of \eqref{eq13},
 we obtain
\begin{equation} \label{eq14}
\begin{aligned}
&\int_{\Omega}[(u^m_t)^2+a_{ij}u^m_{x_i}u^m_{x_j}+(u^m_{xt})^2]\big|_{t=\tau}dx\\
&=2\int_0^{\tau}\int_{\Omega}fu^m_t\,dx\,dt
-2\int_0^{\tau}\int_{\Omega}cu^mu^m_t\,dx\,dt
+\int_0^{\tau}\int_{\Omega}\frac{\partial a_{ij}}{\partial t}u^m_{x_i}u^m_{x_j}
 \,dx\,dt \\
&\quad -2\int_0^{\tau}\int_{\partial \Omega}u^m_t(x,t)
\int_{\Omega}K(x,y,t)u^m(x,t)\,dy\,ds\,dt.
\end{aligned}
\end{equation}
Consider the right-hand side of \eqref{eq14} and focus our attention 
on the term generated by nonlocal conditions.
By applying Cauchy and Cauchy-Bunyakovskii inequalities, we obtain
\begin{align*}
&\big|2\int_0^{\tau}\int_{\partial \Omega}u^m_t(x,t)
\int_{\Omega}K(x,y,t)u^m(y,t)\,dy\,ds\,dt\big|\\
&\leq\int_0^{\tau}\int_{\partial \Omega}(u^m_t(x,t))^2\,ds\,dt+
k\omega\int_0^{\tau}\int_{\Omega}(u^m(x,t))^2\,dx\,dt,
\end{align*}
where $k=\max_{[0,T]}\int_{\Omega}K^2(x,y,t)dy$,
$\omega=\int_{\partial \Omega}ds$.

As the boundary $\partial \Omega$ is smooth \cite{lad}, we have
$$
\int_{\partial \Omega}(u^m_t)^2ds\leq c_1\int_{\Omega}[(u^m_{xt})^2+(u_t^m)^2]dx.
$$
Hence
\begin{equation} \label{eq15}
\begin{aligned}
&2\big|\int_0^{\tau}\int_{\Omega}(u^m_t(x,t))\int_{\Omega}K(x,y,t)u^m(y,t)
 \,dy\,ds\,dt\big| \\
&\leq c_1\int_0^{\tau}\int_{\Omega}[(u_t^m)^2+(u^m_{xt})^2]\,dx\,dt+
k\omega\int_0^{\tau}\int_{\Omega}(u^m)^2\,dx\,dt.
\end{aligned}
\end{equation}
Continue our estimates of right-hand side of \eqref{eq14}. 
As mentioned above there exists $c_0>0$ such that $ |a_{ijt}|, |c|\leq c_0$
and $a_{ij}\xi_i\xi_j\geq \gamma \xi^2$ with $\gamma >0$.
Now we apply Cauchy inequality to estimate the second and the third terms in 
the right-hand side of \eqref{eq14} and obtain
\begin{gather*}
2\big|\int_0^{\tau}\int_{\Omega}cu^mu_t^m\,dx\,dt\big|\leq
 c_0\int_0^{\tau}\int_{\Omega}((u^m)^2+(u_t^m)^2)\,dx\,dt,
\\
2\big|\int_0^{\tau}\int_{\Omega}fu_t^m\,dx\,dt\big|
\leq \int_0^{\tau}\int_{\Omega}f^2\,dx\,dt+
\int_0^{\tau}\int_{\Omega}(u_t^m)^2\,dx\,dt.
\end{gather*}
With this result, from \eqref{eq14} and \eqref{eq15}, we can now obtain
\begin{equation} \label{eq16}
\begin{aligned}
&\int_{\Omega}[(u^m_t)^2+\gamma(u^m_x)^2+(u^m_{xt})^2]\big|_{t=\tau}dx\\
&\leq c_2\int_0^{\tau}\int_{\Omega}[(u^m)^2+(u^m_t)^2+(u^m_x)^2+(u^m_{xt})^2]\,dx\,dt
+\int_0^{\tau}\int_{\Omega}f^2(x,t)\,dx\,dt
\end{aligned}
\end{equation}
It easy to see that the relation
$$
u^m(x,\tau)=\int_0^{\tau}u^m_t(x,t)dt+u^m(x,0)
$$
implies (as $u^m(x,0)=0$) the  inequality
$$
\int_0^l(u^m(x,\tau))^2\,dx\,dt
\leq \tau\int_0^{\tau}\int_{\Omega}(u_t^m(x,t))^2\,dx\,dt.
$$
Adding this inequality to  \eqref{eq16}, we obtain
\begin{equation} \label{eq17}
\begin{aligned}
&m_0\int_{\Omega}[(u^m)^2+(u^m_t)^2+(u^m_x)^2+(u^m_{xt})^2]\big|_{t=\tau}dx\\
&\leq M\int_0^{\tau}\int_{\Omega}[(u^m)^2+(u^m_t)^2+(u^m_x)^2+(u^m_{xt})^2]\,dx\,dt\\
&\quad +N\int_0^l[(u^m(x,0))^2+(u^m_t(x,0))^2+(u^m_x(x,0))^2+(u^m_{xt}(x,0))^2]dx\\
&\quad +\int_0^{\tau}\int_{\Omega}f^2(x,t)\,dx\,dt,
\end{aligned}
\end{equation}
where $M>0$, $N>0$ depend only on $c_0$, $c_1$, $\omega$, $\gamma$, $T$.
By Gronwall's lemma, we conclude that for all $m\geq 1$,
\begin{equation} \label{eq19}
\|u^m\|_{W(Q_T)}\leq P,
\end{equation}
where $P>0$ and does not depend on $m$.

Note that $W(Q_T)$ is Hilbert space. Therefore, because of  \eqref{eq19},
 we can extract from $\{u^m\}$ a subsequence that convergence weakly in $W(Q_T)$
and uniformly with respect to $t\in [0,T]$ in the norm of 
$L_2(\Omega)$ to $u\in W(Q_T)$.
We need only to show that this limit function is a required generalized solution.

Initial condition  $u(x,0)=0$ is fulfilled as $u^m(x,t)\to u(x,t)$
 in $L_2(\Omega)$ uniformly for every $t\in[0,T]$ and
 $u^m(x,0)\to 0$ in $L_2(\Omega)$.
To show that \eqref{eq4} is valid we multiply \eqref{eq10}
 by $d_p\in C^1[0,T]$, $d_p(T)=0$, take sum from $p=1$ to $p=m$ and integrate
 with respect to $t$ from $0$ to $T$.
This leads us to the equality
\begin{align*}
&\int_0^T\int_{\Omega}\left(u^m_{tt}\eta+a_{ij}u^m_{x_i}\eta_{x_j}
+u^m_{x_itt}\eta_{x_j}+cu^m\eta \right)\,dx\,dt
+\int_0^T\int_{\partial \Omega}\eta\int_{\Omega}Ku^m\,dy\,ds\,dt\\
&= \int_0^T\int_{\Omega}f\eta \,dx\,dt.
\end{align*}
Denote $\eta(x,t)=\sum_{p=1}^m d_p(t)w_p(x)$. After integrating by parts 
the terms containing $u^m_{tt}$ and $u^m_{xtt}$,  we obtain
\begin{equation}  \label{eq20}
\begin{aligned}
&\int_0^T\int_{\Omega}\left(-u^m_t\eta_t-u^m_{xt}\eta_{xt}+a_{ij}u^m_{x_i}\eta_{x_j}
+cu^m\eta\right)\,dx\,dt \\
&+\int_0^T\int_{\partial \Omega}\eta(x,t)\int_{\Omega}K(x,y,t)u^m(y,t)\,dy\,ds\,dt\\
&=\int_0^T\int_{\Omega}f\eta \,dx\,dt.
\end{aligned}
\end{equation}

Taking into account the convergence proved above one can pass to the
limit in  \eqref{eq20} as $m\to \infty$ for any fixed $\eta$.
Denote the set of functions $\eta=\sum_{p=1}^m d_p(t)w_p(x)$ by $\mathcal{N}_m$.
As $\cup_{m=1}^{\infty}\mathcal{N}_m$ is dense in $\hat{W}$ \cite{lad},
 it follows that the limit relation is fulfilled
for every function $v\in \hat{W}(Q_T)$, hence,
$u$ is the solution of \eqref{eq1}--\eqref{eq3}.
\end{proof}

\begin{remark} \label{rmk2} \rm
We use homogeneous initial data \eqref{eq2} and a nonlocal condition \eqref{eq3} 
for technical reasons. This involves no loss of generality but simplifies 
computational work. Nonhomogeneous conditions with usual properties 
can also be considered .
In fact, suppose
\begin{gather*}
u(x,0)=\varphi(x), \quad u_t(x,0)=\psi(x), \\
\frac{\partial ^2}{\partial t^2}\frac{\partial u}{\partial \nu}
+\frac{\partial u}{\partial N}
+\int_{\Omega}K(x,y,t)u(y,t)dy=g(x,t), \quad (x,t)\in S_T.
\end{gather*}
The identity \eqref{eq4} becomes
\begin{align*}
&\int_0^T\int_{\Omega}(-u_tv_t-u_{xt}v_{xt}+a_{ij}u_{x_i}v_{x_j}+cuv)\,dx\,dt
+\int_0^T\int_{\partial \Omega}v(0,t)\int_{\Omega}Ku\,dy\,ds\,dt\\
&=\int_0^T\int_{\Omega}fv\,dx\,dt+\int_{\Omega}(\psi v+\psi_xv_x)dx
+\int_0^T\int_{\partial \Omega}gv\,ds\,dt.
\end{align*}
If $\varphi, \psi \in W_2^2(Q_T), g\in L_2(\partial Q_T)$, we are able 
to obtain necessary a priori estimates as above.
\end{remark}

\section{Appendix}

 Here we give some reasons of arising nonlocal condition \eqref{eq3}.
Consider very simple particular case of \eqref{eq1},
\begin{equation}  \label{a1}
u_{tt}-u_{ttxx}-u_{xx}+c(x,t)u=0
\end{equation}
and set a following problem: Find a solution to \eqref{a1} in the domain 
$Q_T=(0,l)\times (0,T)$ such that
\begin{gather}  \label{a2}
u(x,0)=\varphi(x), \quad u_t(x,0)=\psi(x), \quad u_x(0,t)=0, \\
 \label{a3}
\int_0^lu(x,t)dx=0.
\end{gather}
Note that \eqref{a3} is a nonlocal condition of the first kind.
On physical grounds this one or more general integral conditions of the 
form $\int_0^lK(x,t)u(x,t)\,dx$
are very natural (see \cite{samar}--\cite{kamV})  but give rise some
difficulties when we try to prove a solvability of a nonlocal problem 
(see \cite{Pu_1,Pu_2} and references therein).
One method has been advanced for overcoming these difficulties
 in \cite{Pu_1}, \cite{Pu_2} for hyperbolic
equations. The main idea of the procedure is as follows. 
We  reduce the nonlocal condition of the first kind
to a certain nonlocal condition of the second kind.
This method may be applied in a similar way to the problem \eqref{a1}--\eqref{a3}.
We show it in brief.

Let $u(x,t)$ be a solution to  \eqref{a1}--\eqref{a3}. Integrating  
\eqref{a1} with respect to $x$ from $0$ to $l$ we obtain
\begin{equation}  \label{a4}
u_{ttx}(l,t)-u_{x}(l,t)-\int_0^lc(x,t)u(x,t)dx=0.
\end{equation}
It is easy to see that \eqref{a4} is a nonlocal condition of the 
second kind as this relation involves terms outside the integral.

If we assume now that $u(x,t)$ satisfies \eqref{a1}, \eqref{a2},  \eqref{a4} and
the compatibility conditions
$$
\int_0^l\varphi(x)dx=0, \quad \int_0^l\psi(x)dx=0,
$$
then after integrating  \eqref{a1} with respect to $x$
from $0$ to $l$ we obtain
$$
\frac{d^2}{dt^2}\int_0^lu(x,t)dx=0.
$$
The compatibility conditions give us zero initial data for an unknown 
function $\int_0^lu(x,t)dx$, hence $\int_0^lu(x,t)dx=0$.

Now it may be concluded that  problems \eqref{a1}--\eqref{a3} and 
\eqref{a1}, \eqref{a2},  \eqref{a4} are equivalent.

In addition, we can now consider the nonlocal condition \eqref{eq3} 
as a generalization of \eqref{a4}.


\begin{thebibliography}{99}

\bibitem{A-A-G} Avalishvili, G.; Avalishvili, M.; Gordeziani, D.;
\emph{On Integral Nonlocal Boundary Value Problems for some Partial Differential 
Equations}. Bulletin of the Georgian National Academy
of Sciences, 2011, vol. 5, no. 1, p. 31-37.

\bibitem{A-A} Ashyralyev, A.; Aggez, N.;
\emph{On the Solution of Multipoint NBVP for Hyperbolic Equation with
Integral Condition}. Malaysian Journal of Mathematical Sciences 2012, 6(S), p. 111-121.

\bibitem{bouzian} Bouziani, A.;
\emph{Solution Forte d'un Problem Mixte avec Condition  Non Locales pour une
 Classe d'equations Hyperboliques}.
//Bull. de la Classe des Sciences, Academie Royale de Belgique, 1997, no.8, p. 53-70.

 \bibitem{cannon} Cannon, J. R.;
\emph{The solution of the heat equation subject to the specification of energy}. 
Quart. Appl. Math. 1963, vol.21, no.2, p. 155--160.

\bibitem{can_est}  Cannon, John R.;  Perez, Esteva Salvador; van der Hoek;
\emph{A Galerkin Procedure for the diffusion equation subject to the specification 
of mass}. Siam J. Numer. Anal. 1987, vol.24, no.3, p. 449-515.

\bibitem{cannon-lin} Cannon, J. R.; Lin, Y.;
\emph{An Inverse Problem of Finding a Parameter in a Semi-linear Heat
 Equation}.
Journal of Mathematical Analysis and Applications, 1990, no. 145, p. 470-484.

\bibitem{can_lin} Cannon, J. R., Lin, Y.; Wang, Shingmin;
\emph{Determination of a control parameter
in a parabolic partial differential equation}. 
J. Austral. Math. Soc. Ser. B. 1991, no. 33, p. 149-163.

\bibitem{gia}  Gordeziani, D. G.; Avalishvili, G. A.;
\emph{Solutions of Nonlocal Problems for One-dimensional
 Oscillations of the Medium}.
Mat. Modelir., 2000, vol.12, no.1, p.94-103.

\bibitem{kamyn} Kamynin, L. I.;
\emph{On a boundary problem in the theory of heat conduction
with a nonclassical boundary conditions}. 
Zh. Vychisl. Math. Math. Fiz. 1964, vol. 4, no. 6, p. 1006--1024.

\bibitem{kamV} Kamynin, V. L.;
\emph{Unique solvability of the inverse problem of determination of the leading
coefficient in a parabolic equation}. Differential equations,
 2011, vol.47, no. 1, p. 91--101.

\bibitem{Kozh_In} Kozhanov, A. I.;
\emph{Parabolic equations with an unknown time-dependent coefficient}.
 Comput. Math. and Math. Phys. 2005, no.12, p. 2085--2101.

 \bibitem{K-P} Kozhanov, A. I.; Pulkina, L. S.;
\emph{On the Solvability of Boundary Value Problems
with a Nonlocal Boundary Condition of Integral Form for Multidimentional 
Hyperbolic Equations}.
Differential Equations, 2006, vol.42, no.9, p. 1233--1246.

\bibitem{kozh_S} Kozhanov, A. I.;
\emph{On Solvability of Certain Spatially Nonlocal Boundary Problems for Linear
 Parabolic Equations}.
Vestnik of Samara State University, 2008, no.3, p. 165-174.

\bibitem{Korp} Korpusov, O. M.;
\emph{Blow-up in nonclassical wave equations}. Moscow, URSS. 2010.

\bibitem{lad} Ladyzhenskaya, O. A.;
\emph{Boundary-value problems of mathematical physics}. Moscow: Nauka, 1973.

\bibitem{Pu_3} Pulkina, L. S.;
\emph{A mixed problem with integral condition for the hyperbolic equation} 
Mathematical Notes, 2003, vol. 74, no.3, p. 411--421.

\bibitem{Pu_04} Pulkina, L. S.;
\emph{A Nonlocal Problem with Integral Conditions for a Hyperbolic Equation}.
 Differential equations, 2004, vol.40, no. 7, p. 947-953.


\bibitem{Pu_08} Pulkina, L. S.;
\emph{Initial-Boundary Value Problem with a Nonlocal Boundary Condition
 for a Multidimensional Hyperbolic equation}.
Differential equations, 2008, vol.44, no. 8, p. 1119-1125.


\bibitem{Pu_4} Pulkina, L.S. A nonlocal problem with integral conditions for the hyperbolic equation //
Nanosystems: physics, chemistry, mathematics. 2011, no.2(4), p. 61--70.

\bibitem{Pu_1} Pulkina, L. S.;
\emph{Boundary value problems for a hyperbolic equation with nonlocal conditions 
of the I and II kind},  Russian Mathematics (Iz. VUZ) 2012, vol. 56, no. 4, p. 62--69.

\bibitem{Pu_2} Pulkina, L.S.;
\emph{A nonlocal problem for a hyperbolic equation with integral conditions 
of the Ist  kind with time-dependent kernels}. 
Russian Mathematics (Iz. VUZ) 2012, vol. 56, no.10. P. 26--37.

\bibitem{samar} Samarskii, A. A.;
\emph{On Certain Problems of the Modern Theory of Differential equations}.
Differ. Uravn., 1980, vol. 16, no. 11, p. 1221-1228.

\bibitem{sobolev} Sobolev, S. L.;
\emph{On certain new problem of mathematical physics}.
 Izvestia Acad. Nauk.  Mathem. 1954, no. 18, p. 3--50.

\bibitem{sv} Sveshnikov, A. G.; Alshin, A. B.; Korpusov, O. M.; Pletner, Yu. D.;
\emph{Linear and nonlinear Sobolev type equations}. Moscow., Phizmatlit. 2007.

\end{thebibliography}

\end{document}
