\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2014 (2014), No. 74, pp. 1--11.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2014 Texas State University - San Marcos.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2014/74\hfil Pullback attractor]
{Pullback attractor for  non-autonomous
$p$-Laplacian equations with dynamic flux boundary conditions}

\author[B. You, F. Li \hfil EJDE-2014/74\hfilneg]
{Bo You, Fang Li}

\address{Bo You \newline
School of Mathematics and Statistics, Xi'an Jiaotong University\\
Xi'an, 710049,  China}
\email{youb03@126.com}

\address{Fang Li \newline
Department of Mathematics, Nanjing University\\
Nanjing, 210093,  China}
\email{lifang101216@126.com}

\thanks{Submitted June 26, 2013. Published March 18, 2014.}
\subjclass[2000]{35B40, 37B55}
\keywords{Pullback attractor; Sobolev compactness embedding;  
$p$-Laplacian; 
\hfill\break\indent norm-to-weak continuous process;
 asymptotic a priori estimate; non-autonomous;
\hfill\break\indent nonlinear flux boundary conditions}

\begin{abstract}
 This article studies the long-time asymptotic behavior of solutions for
 the non-autonomous $p$-Laplacian equation
 \[
 u_t-\Delta_pu+ |u|^{p-2}u+f(u)=g(x,t)
 \]
 with dynamic flux boundary conditions
 \[
 u_t+|\nabla u|^{p-2}\frac{\partial u}{\partial\nu}+f(u)=0
 \]
 in a $n$-dimensional bounded smooth  domain $\Omega$ under some
 suitable assumptions. We prove the existence of a pullback
 attractor in $\big(W^{1,p}(\Omega)\cap L^q(\Omega)\big)\times L^q(\Gamma)$
 by asymptotic a priori estimate.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{remark}[theorem]{Remark}
\allowdisplaybreaks

\section{Introduction}

We are concerned with the existence of a pullback attractor in 
$\big(W^{1,p}(\Omega)\cap L^q(\Omega)\big)\times L^q(\Gamma)$ 
for the process $\{U(t,\tau)\}_{t\geq\tau}$ associated
with solutions of the following non-autonomous $p$-Laplacian equation
\begin{equation}\label{1}
u_t-\Delta_pu+ |u|^{p-2}u+f(u)=g(x,t),\quad (x,t)\in\Omega\times
\mathbb{R}_{\tau}.
\end{equation}
This equation is subject to the dynamic flux boundary condition
\begin{equation}\label{2}
u_t+|\nabla u|^{p-2}\frac{\partial u}{\partial\nu}+f(u)=0,\quad
\quad (x,t)\in\Gamma\times\mathbb{R}_{\tau}
\end{equation}
and the initial conditions
\begin{gather}
\label{3} u(x,\tau)=u_\tau(x),\quad x\in\Omega,\\
\label{4} u(x,\tau)=\theta_\tau(x),\quad x\in\Gamma,
\end{gather}
where  $\Omega \subset \mathbb{R}^n$ $(n\geq 3)$ is a bounded domain
with smooth boundary $\Gamma$, $\nu$ denotes the outer unit
normal on $\Gamma$, $p\geq 2$, $\mathbb{R}_\tau=[\tau,+\infty)$, 
the nonlinearity $f$ and the external force $g$ satisfy
some  conditions, specified later.

 To study problem \eqref{1}-\eqref{4}, we assume the following conditions:
\begin{itemize}
\item[(H1)] the function $f\in C^1(\mathbb{R},\mathbb{R})$
 and satisfies
\begin{equation}\label{5}
f'(u)\geq -l
\end{equation}
for some $l\geq 0$, and
\begin{equation}\label{6}
 c_1|u|^q-k\leq f(u)u\leq c_2|u|^q+k,
\end{equation}
where $c_i> 0$ ($i=1,2$), $q\geq 2$, $ k> 0$.

\item[(H2)] The external force $g:\Omega\times \mathbb{R}\to \mathbb{R}$ is locally
Lipschitz continuous,  $ g$ belongs to $H_{\rm loc}^1(\mathbb{R}, L^2(\Omega))$,
  and satisfies 
  \begin{equation}\label{7}
\int_{-\infty}^te^{c_1 s}\|g(s)\|_{L^2(\Omega)}^2\,ds
+\int_{-\infty}^te^{c_1 s}\|g_t(s)\|_{L^2(\Omega)}^2\,ds<\infty
\end{equation}
 for all $t\in\mathbb{R}$.
\end{itemize}

Dynamic boundary conditions are very natural in many mathematical
 models such as heat transfer in a solid in contact with a moving fluid,
thermoelasticity, diffusion phenomena, heat transfer in two
 medium, problems in fluid dynamics (see  
\cite{ ajm, ajm1, bar, ca, ca1, fzh, pj, pl, yl, yl1}). 
The understanding of the asymptotic behavior of dynamical systems is one of 
the most important problems of modern mathematical physics. One way to
treat this problem for a dissipative system is to analyze the existence and 
structure of its attractor. Generally speaking, the attractor has a very 
complicated geometry which reflects the complexity of the long-time behavior of 
the system.  There are many authors who have considered the long-time behavior 
of solutions for the problems of dynamic boundary conditions. For example, 
the authors considered the existence of global attractors, respectively, 
in $L^2(\bar{\Omega},d\mu)$, $L^q(\bar{\Omega},d\mu)$ and 
$\big(H^1(\Omega)\cap L^q(\Omega)\big)\times L^q(\Gamma)$ for the 
reaction-diffusion equation with dynamic flux boundary conditions in \cite{fzh}.  
The existence of uniform attractors in $L^2(\bar{\Omega},d\mu)$ and 
$\big(H^1(\Omega)\cap L^q(\Omega)\big)\times L^q(\Gamma)$ 
for the reaction-diffusion equation with dynamic flux boundary conditions was 
proved in \cite{yl}.
 In \cite{yb}, the authors proved the existence of global attractors for 
the autonomous $p$-Laplacian equation with dynamic flux boundary conditions 
in $L^2(\bar{\Omega},d\mu)$, $L^q(\bar{\Omega},d\mu)$ by the Sobolev compactness 
embedding theorem and the existence of a global attractor for the autonomous 
$p$-Laplacian equation with dynamic flux boundary conditions in 
$\big(W^{1,p}(\Omega)\cap L^q(\Omega)\big)\times L^q(\Gamma)$ by asymptotical 
a priori estimate.  Recently, the existence of uniform attractors in 
$L^2(\bar{\Omega},d\mu)$ and 
$\big(W^{1,p}(\Omega)\cap L^q(\Omega)\big)\times L^q(\Gamma)$ 
for the non-autonomous $p$-Laplacian equation with dynamic flux boundary 
conditions was obtained in \cite{lk}.

Non-autonomous equations appear in many applications in  natural
sciences, so they are of great importance and interest. 
The long-time behavior of solutions for the non-autonomous equations 
has been studied
extensively in recent years (see 
\cite{ dn,hc, hc1,ct, pe, pe1, lss, bs, yl}). 
For instance, the existence of a pullback attractor in $L^2(\Omega)$ 
was studied in \cite{ct1}. The authors obtained the existence of 
a pullback attractor in $H_0^1(\Omega)$ in \cite{sht1}. 
The existence of a pullback attractor in $H_0^1(\Omega)$ was considered 
in \cite{ly}. The authors proved the existence of a pullback attractor 
in $L^p(\Omega)$ for a reaction-diffusion equation in \cite{ly1} under 
the assumption
\begin{equation*}
\|g(s)\|_2^2\leq Me^{\alpha |s|}
\end{equation*}
for all $s\in\mathbb{R}$ and $0\leq\alpha<\lambda_1$, where $\lambda_1$ 
is the first eigenvalue of $-\Delta$ with Dirichlet boundary condition. 
In \cite{yl1}, the authors used a new type of uniform Gronwall inequality 
and proved the existence of a pullback attractor in 
$L^{r_1}(\Omega)\times L^{r_2}(\Gamma)$ for the  equation
  \begin{gather*}
u_t-\Delta_pu+ |u|^{p-2}u+f(u)=h(t),\quad (x,t)\in\Omega\times\mathbb{R}_\tau,\\
u_t+|\nabla u|^{p-2}\frac{\partial u}{\partial\nu}+g(u)=0,\quad
(x,t)\in\Gamma\times\mathbb{R}_\tau,\\
u(x,\tau)=u_0(x),\quad x\in\bar{\Omega}
\end{gather*}
under the assumptions that $f$, $g$ satisfy the polynomial growth condition 
with orders $r_1$, $r_2$ and $\|h(t)\|_{L^2(\Omega)}$ satisfies some weak 
assumption
\begin{equation*}
\int_{-\infty}^te^{\theta s}\|h(s)\|_{L^2(\Omega)}^2\,ds<\infty
\end{equation*}
for all $t\in\mathbb{R}$, where $\theta$ is some positive constant. 
By using their main result, we can get the following result.

  \begin{corollary}\label{51}
  Let $\Omega$ be a bounded domain in $\mathbb{R}^n$ with smooth boundary 
$\Gamma$, let $f$ and $g$ satisfy {\rm (H1)--(H2)}.
Then the process $\{U(t,\tau)\}_{t\geq \tau}$ corresponding to
\eqref{1}-\eqref{4} has a pullback $\mathcal{D}$-attractor $\mathcal{A}_q$ 
in $L^q(\bar{\Omega},d\mu)$, which is
 pullback $\mathcal{D}$-attracting in the topology of 
$L^q(\bar{\Omega},d\mu)$-norm.
 \end{corollary}

The study of non-autonomous dynamical systems is an important subject, 
it is necessary to study the existence of pullback attractors for the 
non-autonomous $p$-Laplacian equation with dynamic flux boundary conditions. 
Nevertheless, there are few results about the existence of a pullback 
attractor in $(W^{1,p}(\Omega)\cap L^q(\Omega))\times L^q(\Gamma)$ 
for the non-autonomous $p$-Laplacian equation with dynamic flux boundary 
conditions. The main difficulty is that in our case of the equation with 
$p$-Laplacian operator for $p>2$, we cannot use $-\Delta u_2$ as the test 
function to verify pullback $\mathcal{D}$-condition, which increases the 
difficulty in getting an appropriate form of compactness. 
To overcome this difficulty, we combine the idea of norm-to-weak process 
with asymptotic a priori estimates to prove the existence of a pullback 
attractor for the non-autonomous $p$-Laplacian equation with dynamic flux 
boundary conditions in $\big(W^{1,p}(\Omega)\cap L^q(\Omega)\big)\times L^q(\Gamma)$.



   The main purpose of this paper is to study the existence of a pullback 
attractor in $(W^{1,p}(\Omega)\cap L^q(\Omega))\times L^q(\Gamma)$ 
for the non-autonomous $p$-Laplacian evolutionary equation \eqref{1}-\eqref{4} 
under quite general assumptions \eqref{5}-\eqref{7}. 
Here, we state our main result as follows.

\begin{theorem}\label{main}
Assume that {\rm (H1)--(H2)} hold. Then the process $\{U(t,\tau)\}_{t\geq \tau}$
corresponding to problem \eqref{1}-\eqref{4}
has a pullback $\mathcal{D}$-attractor $\mathcal {A}$ in 
$(W^{1,p}(\Omega)\cap L^q(\Omega))\times L^q(\Gamma)$.
\end{theorem}


This article is organized as follows: In the next section, we give some 
notation and lemmas used in the sequel.
Section 3 is devoted to proving the existence of a pullback absorbing 
set in $\big(L^2(\Omega)\cap W^{1,p}(\Omega)\cap L^q(\Omega)\big)
\times\big(L^2(\Gamma)\cap L^q(\Gamma)\big)$ and the existence of a 
pullback attractor in $\big(L^q(\Omega)\cap W^{1,p}(\Omega)\big)\times L^q(\Gamma)$.

Throughout this paper, let $C$ be a positive constant, which may be different 
from line to line (and even in the same line),
we denote the trace of $u$ by $v$.

\section{Preliminaries}

To study  \eqref{1}-\eqref{4}, we recall the Sobolev space
 $W^{1,p}(\Omega)$ defined as the closure of 
$C^{\infty}(\Omega)\cap W^{1,p}(\Omega)$ in the norm
\begin{equation*}
\|u\|_{1,p}=\Big(\int_{\Omega}|\nabla u|^p+|u|^p\,dx\Big)^{1/p}
\end{equation*}
and denote by $X^*$ the dual space of $X$. We also define the Lebesgue 
spaces as follows
\begin{equation*}
L^r(\Gamma)=\{v:\|v\|_{L^r(\Gamma)}<\infty\},
\end{equation*}
where
\begin{equation*}
\|v\|_{L^r(\Gamma)}=\Big(\int_{\Gamma}|v|^r\,dS\Big)^{1/r}
\end{equation*}
for $r\in [1,\infty)$. Moreover, we have
\begin{gather*}
L^s(\Omega)\oplus L^s(\Gamma)=L^s(\bar{\Omega},d\mu),\quad s\in [1,\infty),\\
\|U\|_{L^s(\bar{\Omega},d\mu)}=\Big(\int_{\Omega}|u|^s\,dx\Big)^{1/s}
+\Big(\int_{\Gamma}|v|^s\,dS\Big)^{1/s}
\end{gather*}
for any $U=\begin{pmatrix}u\\ v\end{pmatrix}\in L^s(\bar{\Omega},d\mu)$,
 where the measure 
$d\mu=dx|_{\Omega}\oplus dS|_{\Gamma}$ on $\bar{\Omega}$ is defined for 
any measurable set $A\subset \bar{\Omega}$ by 
$\mu(A)=|A\cap\Omega|+S(A\cap \Gamma)$. In general, any vector 
$\theta\in L^s(\bar{\Omega},d\mu)$ will be of the form $\begin{pmatrix} \theta_1 \\
\theta_2\end{pmatrix}$ 
with $\theta_1\in L^s(\Omega,dx)$ and $\theta_2\in L^s(\Gamma,dS)$, 
and there need not be any connection between $\theta_1$ and $\theta_2$.

\begin{remark}[\cite{gcg}] \label{80} \rm
$C(\bar{\Omega})$ is a dense subspace of $L^2(\bar{\Omega},d\mu)$ and a 
closed subspace of $L^{\infty}(\bar{\Omega},d\mu)$.
\end{remark}

Next, we recall briefly some lemmas used to prove the existence of pullback 
absorbing sets for \eqref{1}-\eqref{4} under some suitable assumptions.

\begin{lemma}[\cite{bt}] \label{11}
 Let $x,y\in\mathbb{R}^n$ and let $\langle\cdot,\cdot\rangle$ be the standard 
scalar product in $\mathbb{R}^n$. Then for any $p\geq 2$, there exist two 
positive constants $C_1$, $C_2$ which depend on $p$ such that
\begin{gather*}
\langle|x|^{p-2}x-|y|^{p-2}y,x-y\rangle\geq C_1|x-y|^p,\\
\big|\,|x|^{p-2}x-|y|^{p-2}y\big|\leq C_2(|x|+|y|)^{p-2}|x-y|.
\end{gather*}
\end{lemma}

\section{Existence of pullback attractors}

In this section, we prove the existence of pullback attractors of solutions 
for problem \eqref{1}-\eqref{4}.

\subsection{Well-posedness of solutions for problem \eqref{1}-\eqref{4}}

In this subsection, we give the well-posedness of solutions for problem 
\eqref{1}-\eqref{4} which can be obtained by the Faedo-Galerkin method 
(see \cite{tr}). Here, we only state it as follows.

\begin{theorem}\label{14}
Under the assumptions {\rm (H1)--(H2)}, for any initial data 
$(u_\tau,\theta_\tau)\in L^2(\bar{\Omega},d\mu)$, there exists a unique weak
 solution $u(x,t)\in C(\mathbb{R}_\tau;L^2(\bar{\Omega},d\mu))$ of
 problem \eqref{1}-\eqref{4} and the mapping
\begin{equation*}
(u_\tau,\theta_\tau)\to (u(t),v(t))
\end{equation*}
is continuous on $L^2(\bar{\Omega},d\mu)$.
\end{theorem}


 By Theorem \ref{14}, we can define a family of continuous processes
$\{U(t,\tau):-\infty<\tau\leq t<\infty\}$ in $L^2(\bar{\Omega},d\mu)$ as follows: 
for all $t\geq\tau$,
\begin{equation*}
U(t,\tau)(u_\tau,\theta_\tau)=(u(t),v(t))
:=(u(t;\tau,(u_\tau,\theta_\tau)),v(t;\tau,(u_\tau,\theta_\tau))),
\end{equation*}
where $u(t)$ is the solution of problem \eqref{1}-\eqref{4} with initial data 
$(u(\tau),v(\tau))=(u_\tau, \theta_\tau) \in L^2(\bar{\Omega},d\mu)$.
 That is, a family of mappings 
$U(t,\tau): L^2(\bar{\Omega},d\mu)\to L^2(\bar{\Omega},d\mu)$ satisfies
\begin{gather*}
U(\tau,\tau)=id\quad \text{(identity)},\\
U(t,\tau)=U(t,r)U(r,\tau)\quad\text{for all }\tau\leq r\leq t.
\end{gather*}

\subsection{Existence of a pullback absorbing set}

In this subsection, we recall some basic definitions and abstract results 
about pullback attractors.

\begin{definition}[\cite{ly, yl}] \rm
 Let $X$ be a Banach space. A process $\{U(t,\tau)\}_{t\geq\tau}$ 
is said to be norm-to-weak
continuous on $X$, if for any $t$, $\tau\in\mathbb{R}$ with $t\geq\tau$ 
and for every sequence $x_n\in X$, from the condition $x_n\to x$ strongly 
in $X$, it follows that   $U(t,\tau)x_n\to U(t,\tau)x$ weakly in $X$.
\end{definition}

\begin{lemma}[\cite{ly, yl}] \label{33}
Let $X$ and $Y$ be two Banach spaces, and let $X^{*}$ and $Y^{*}$ be the dual 
spaces of $X$ and $Y$, respectively. If $X$ is dense in $Y$, the
injection $i: X\to Y$ is continuous and its adjoint $i^{*}: Y^{*}\to X^{*}$ is dense.
In addition, assume that $\{U(t,\tau)\}_{t\geq \tau}$ is a continuous or 
weak continuous process on $Y$. Then $\{U(t,\tau)\}_{t\geq \tau}$ is a norm-to-weak
continuous process on $X$ if and only if $\{U(t,\tau)\}_{t\geq \tau}$ maps 
compact sets of $X$ into bounded sets of $X$ for any $t$, $\tau\in\mathbb{R}$, 
$t\geq\tau$.
\end{lemma}

Let $\mathcal{D}$ be a nonempty class of families 
$\hat{D}=\{D(t):t\in\mathbb{R}\}$ of nonempty subsets of $X$.

\begin{definition}[\cite{ct}] \rm
The process $\{U(t,\tau)\}_{t\geq\tau}$ is said to be pullback 
$\mathcal{D}$-asymp\-totically compact, if for any $t\in\mathbb{R}$ and any 
$\hat{D}\in\mathcal{D}$, any sequence $\tau_n\to -\infty$ and any
sequence $x_n\in D(\tau_n)$, the sequence $\{U(t,\tau_n)x_n\}_{n=1}^{\infty}$ 
is relatively compact in $X$.
\end{definition}

\begin{definition}[\cite{yl}] \rm
A family $\hat{\mathcal{A}}=\{A(t):t\in\mathbb{R}\}$ of nonempty subsets of 
$X$ is said to be a pullback $\mathcal{D}$-attractor for the process 
$\{U(t,\tau)\}_{t\geq\tau}$ in $X$, if
\begin{itemize}
\item [(i)] $A(t)$ is  compact in $X$ for any $t\in\mathbb{R}$,
\item [(ii)] $\hat{\mathcal{A}}$ is invariant, i.e., $U(t,\tau)A(\tau)=A(t)$ 
 for any $\tau\leq t$,
\item [(iii)] $\hat{\mathcal{A}}$ is pullback $\mathcal{D}$-attracting, i.e.,
\begin{equation*}
\lim_{\tau\to -\infty}\operatorname{dist}(U(t,\tau)D(\tau),A(t))=0
\end{equation*}
for any $t\in\mathbb{R}$ and any $\hat{D}=\{D(t):t\in\mathbb{R}\}\in\mathcal{D}$.
\end{itemize}
Such a family $\hat{\mathcal{A}}$ is called minimal if $A(t)\subset C(t)$ 
for any family $\hat{C}=\{C(t):t\in\mathbb{R}\}$ of closed subsets of $X$ 
such that $\lim_{\tau\to -\infty}\operatorname{dist}(U(t,\tau)B(\tau),C(t))=0$ 
for any $\hat{B}=\{B(t):t\in\mathbb{R}\}\in\mathcal{D}$.
\end{definition}

\begin{definition}[\cite{ct, yl}] \rm
It is said that $\hat{B}\in\mathcal{D}$ is pullback $\mathcal{D}$-absorbing 
for the process $\{U(t,\tau)\}_{t\geq\tau}$, if for any $\hat{D}\in\mathcal{D}$ 
and any $t\in\mathbb{R}$, there exists a $\tau_0(t,\hat{D})\leq t$ such that
 $U(t,\tau)D(\tau)\subset B(t)$ for any $\tau\leq\tau_0(t,\hat{D})$.
\end{definition}

\begin{lemma}[\cite{ct, ly, yl}] \label{34}
Let $\{U(t,\tau)\}_{t\geq\tau}$ be a process in $X$ satisfying the following 
conditions:
\begin{itemize}
\item [(1)] $\{U(t,\tau)\}_{t\geq\tau}$ be norm-to-weak continuous in $X$.
\item [(2)] There exists a family $\hat{B}$ of pullback $\mathcal{D}$-absorbing 
 sets $\{B(t):t\in\mathbb{R}\}$ in $X$.
\item [(3)] $\{U(t,\tau)\}_{t\geq\tau}$ is pullback $\mathcal{D}$-asymptotically
  compact.
\end{itemize}
Then there exists a minimal pullback $\mathcal{D}$-attractor 
$\hat{\mathcal{A}}=\{A(t):t\in\mathbb{R}\}$ in $X$ given by
\begin{equation*}
A(t)= \cap_{s\leq t}\overline{\cup_{\tau\leq s}U(t,\tau)B(\tau)}.
\end{equation*}
\end{lemma}

\begin{lemma}[\cite{yl}]\label{36}
Suppose that
\begin{equation*}
y'(s)+\delta y(s)\leq b(s)
\end{equation*}
for some $\delta>0$, $t_0\in\mathbb{R}$ and for any $s\geq t_0$, where the 
functions $y$, $y'$, $b$ are assumed to be locally integrable and $y$, $b$ 
are nonnegative on the interval $t<s<t+r$ for some $t\geq t_0$. Then
 \begin{equation*}
y(t+r)\leq e^{-\frac{\delta r}{2}}\frac{2}{r}\int_t^{t+\frac{r}{2}}y(s)\,ds
+e^{-\delta(t+r)}\int_t^{t+r} e^{\delta s}b(s)\,ds
 \end{equation*}
 for all $t\geq t_0$.
\end{lemma}


In the following, let $\mathcal{D}$ be the class of all families 
$\{D(t):t\in\mathbb{R}\}$ of nonempty subsets of $L^2(\bar{\Omega},d\mu)$ such that
\begin{equation*}
\lim_{t\to-\infty} e^{c_1 t}[D(t)]=0,
\end{equation*}
where $[D(t)]=\sup\{\|(u,v)\|_{L^2(\bar{\Omega},d\mu)}: (u,v)\in D(t)\}$. 
We prove the existence of a pullback absorbing
set for the process $\{U(t,\tau)\}_{t\geq\tau}$  corresponding to problem 
\eqref{1}-\eqref{4}.


\begin{theorem}\label{37}
Under assumptions {\rm (H1)--(H2)}. Let $\{U(t,\tau)\}_{t\geq\tau}$ 
be a process associated with problem \eqref{1}-\eqref{4}. 
Then there exists a pullback $\mathcal{D}$-absorbing set in
 $\big(L^2(\Omega)\cap W^{1,p}(\Omega)\cap L^q(\Omega)\big)
\times\big(L^2(\Gamma)\cap L^q(\Gamma)\big)$.
\end{theorem}

\begin{proof}
 Taking the inner product of \eqref{1} with $u$, we
deduce that
\begin{equation}
\begin{aligned}
 &\frac{1}{2}\frac{d}{dt}\left(\|u\|_{L^2
(\Omega)}^2+\|v\|_{L^2
(\Gamma)}^2\right)+\|u\|_{W^{1,p}}^p+\int_{\Omega}f(u)u\,dx
 + \int_{\Gamma}f(v)v\,dS\\
 &=\int_{\Omega}g(t)u\,dx.
\end{aligned}
\end{equation}
By \eqref{6}, H\"{o}lder inequality and Young inequality, we obtain
\begin{align*}
&\frac{1}{2}\frac{d}{dt}\left(\|u\|_{L^2
 (\Omega)}^2+\|v\|_{L^2
 (\Gamma)}^2\right)+\|u\|_{W^{1,p}(\Omega)}^p+c_1\|u\|_{L^q(\Omega)}^q
 +c_1\|v\|_{L^q(\Gamma)}^q\\
&\leq\frac{1}{2}\|
 g(t)\|^2 _{L^2(\Omega)}+\frac{1}{2}\|
 u\|^2_{L^2(\Omega)}+k|\Omega|+k|\Gamma|\\
&\leq \frac{1}{2}\| g(t)\|^2
_{L^2(\Omega)}+\frac{1}{2}\|
u\|^2_{L^2(\Omega)}+\frac{1}{2}\|
v\|^2_{L^2(\Gamma)}+k|\Omega|+k|\Gamma|.
\end{align*}
Therefore,
\begin{equation}\label{40}
\begin{aligned}
 &\frac{d}{dt}\left(\|u\|_{L^2
(\Omega)}^2+\|v\|_{L^2
(\Gamma)}^2\right)+2\|u\|_{W^{1,p}(\Omega)}^p+2c_1\|u\|_{L^q(\Omega)}^q
 +2c_1\|v\|_{L^q(\Gamma)}^q\\
&\leq\| g(t)\|^2_{L^2(\Omega)}+\|u\|^2_{L^2(\Omega)}+\|
v\|^2_{L^2(\Gamma)}+2k|\Omega|+2k|\Gamma|.
\end{aligned}
\end{equation}
It follows from \eqref{40} that
\begin{equation}\label{41}
\begin{aligned}
&\frac{d}{dt}\left(\|u\|_{L^2
(\Omega)}^2+\|v\|_{L^2 (\Gamma)}^2\right)+c_1\left(\|u\|_{L^2
(\Omega)}^2+\|v\|_{L^2 (\Gamma)}^2\right)\\
&\quad +2\|u\|_{W^{1,p}(\Omega)}^p+c_1\|u\|_{L^q(\Omega)}^q
 +c_1\|v\|_{L^q(\Gamma)}^q\\
&\leq \| g(t)\|^2_{L^2(\Omega)}+C.
\end{aligned}
\end{equation}
From the classical Gronwall inequality, we find that
\begin{equation}\label{42}
\begin{aligned}
 &\|u(t)\|_{L^2(\Omega)}^2+\|v(t)\|_{L^2(\Gamma)}^2\\
&\leq \left(\| u_\tau\|_{L^2(\Omega)}^2+\|\theta_\tau\|_{L^2(\Gamma)}^2\right)
 e^{c_1(\tau-t)}+e^{-c_1t}\int_{-\infty}^te^{c_1s}\|g(s)\|^2_{L^2(\Omega)}\,ds+C,
\end{aligned}
\end{equation}
which implies
\begin{equation}\label{43}
\| u(t)\|_{L^2(\Omega)}^2+\|v(t)\|_{L^2(\Gamma)}^2
\leq \mathcal{C}_0\Big(e^{-c_1t}\int_{-\infty}^te^{c_1s}
\|g(s)\|^2_{L^2(\Omega)}\,ds+1\Big)
\end{equation}
 uniformly with respect to all initial conditions $(u_\tau,v_\tau)\in D(\tau)$ 
for $\tau\leq\tau_0(t,\hat{D})$, where $\mathcal{C}_0$ is a positive constant.

 Let $F(s)=\int_{0}^{s}f(\theta)\,d\theta$, we deduce from \eqref{6} that 
there exist three positive constants $\alpha_1$, $\alpha_2$, $\beta$ such that
\begin{gather}
\alpha_1|u|^q-\beta\leq F(u)\leq \alpha_2|u|^q+\beta, \nonumber \\
\label{44}\alpha_1|u|_{L^q(\Omega)}^q-\beta|\Omega|\leq\int_{\Omega}F(u)\,dx
\leq \alpha_2|u|_{L^q(\Omega)}^q+\beta|\Omega|,\\
\label{45}\alpha_1|v|_{L^q(\Gamma)}^q-\beta|\Gamma|\leq\int_{\Gamma} F(v)\,dS
\leq \alpha_2|v|_{L^q(\Gamma)}^q+\beta|\Gamma|.
\end{gather}
Integrating \eqref{41} from $t$ to $t+1$ and combining \eqref{42} with 
\eqref{44}-\eqref{45}, we obtain
 %\label{46}
\begin{align*}
&2\int_t^{t+1}\|u(s)\|_{W^{1,p}(\Omega)}^p\,ds
 +\frac{c_1}{\alpha_2}\int_t^{t+1}\int_{\Omega}F(u(s))\,dx\,ds
 +\frac{c_1}{\alpha_2}\int_t^{t+1} \int_{\Gamma}F(v(s))\,dS\,ds\\
&\leq \mathcal{C}_0\Big(e^{-c_1t}\int_{-\infty}^te^{c_1s}\|g(s)\|^2_{L^2(\Omega)}\,ds
+1\Big)+\int_t^{t+1}\| g(s)\|^2_{L^2(\Omega)}\,ds+C\\
&\leq \mathcal{C}_1\Big(e^{-c_1t}\int_{-\infty}^te^{c_1s}\|g(s)\|^2_{L^2(\Omega)}\,ds
+1\Big)
\end{align*}
uniformly with respect to all initial conditions $(u_\tau,v_\tau)\in D(\tau)$ 
for $\tau\leq\tau_0(t,\hat{D})$, where $\mathcal{C}_1$ is a positive constant.

Taking the inner product of \eqref{1} with $u_{t}$, we obtain
\begin{align*}
&\| u_t\|_{L^2 (\Omega)}^2+\|v_t\|_{L^2(\Gamma)}^2+\frac{d}{dt}
\Big(\frac{1}{p}\|u\|_{W^{1,p}(\Omega)}^p+\int_{\Omega}F(u)\,dx
+\int_{\Gamma}F(v)\,dS\Big)\\
&= \int_{\Omega}g(x,t)u_t\,dx\\
&\leq \frac{1}{2}\|g(t)\|^2 _{L^2(\Omega)}+\frac{1}{2}\|u_t\|^2_{L^2(\Omega)},
\end{align*}
which implies
\begin{equation}\label{47}
\begin{aligned}
&\| u_t\|_{L^2 (\Omega)}^2+\|v_t\|_{L^2(\Gamma)}^2
 +\frac{d}{dt}\Big(\frac{2}{p}\|u\|_{W^{1,p}(\Omega)}^p
 +2\int_{\Omega}F(u)\,dx+2\int_{\Gamma}F(v)\,dS\Big)\\
&\leq \|g(t)\|^2 _{L^2(\Omega)}.
\end{aligned}
\end{equation}
It follows from the uniform Gronwall inequality that
\begin{equation}\label{48}
\begin{aligned}
 &\|u(t+1)\|_{W^{1,p}(\Omega)}^p+\int_{\Omega}F(u(t+1))\,dx
+\int_{\Gamma}F(v(t+1))\,dS\\
&\leq \mathcal{C}_2\Big(e^{-c_1t}\int_{-\infty}^te^{c_1s}\|g(s)\|^2_{L^2(\Omega)}\,ds
+1\Big)
\end{aligned}
\end{equation}
 uniformly with respect to all initial conditions $(u_\tau,v_\tau)\in D(\tau)$ 
for $\tau\leq\tau_0(t,\hat{D})$, where $\mathcal{C}_2$ is a positive constant.


We infer from \eqref{44}-\eqref{45} and \eqref{48} that
\begin{equation}\label{49}
\begin{aligned}
&\|u(t+1)\|_{W^{1,p}(\Omega)}^p+\|u(t+1)\|_{L^q(\Omega)}^q
 +\|v(t+1)\|_{L^q(\Gamma)}^q\\
&\leq \mathcal{C}_3\Big(e^{-c_1t}\int_{-\infty}^te^{c_1s}
\|g(s)\|^2_{L^2(\Omega)}\,ds +1\Big)
\end{aligned}
\end{equation}
 uniformly with respect to all initial conditions $(u_\tau,v_\tau)\in D(\tau)$ 
for $\tau\leq\tau_0(t,\hat{D})$, where $\mathcal{C}_3$ is a positive constant.
\end{proof}

Since $W^{1,p}(\Omega)\hookrightarrow L^2(\Omega)$ and 
$W^{1,p}(\Omega)\hookrightarrow L^2(\Gamma)$ is compact, we obtain the 
following result.

\begin{theorem}\label{50}
Under the assumptions {\rm (H1)--(H2)}, the process 
$\{U(t,\tau)\}_{t\geq \tau}$ corresponding to problem \eqref{1}-\eqref{4}
 has a pullback $\mathcal{D}$-attractor $\mathcal {A}_2$ in 
$L^2(\bar{\Omega},d\mu)$, which is compact, connected and invariant.
\end{theorem}

\subsection{Existence of a pullback attractor in
 $(W^{1,p}(\Omega)\cap L^q(\Omega))\times L^q(\Gamma)$}
 From Lemma \ref{33} and Theorem \ref{37}, we know that
the process $\{U(t, \tau)\}_{t\geq \tau}$ corresponding to problem 
\eqref{1}-\eqref{4} is norm-to-weak continuous in 
$(W^{1,p}(\Omega)\cap L^q(\Omega))\times
L^q(\Gamma)$. In this subsection, we prove the existence of a pullback 
$\mathcal{D}$-attractor in $(W^{1,p}(\Omega)\cap L^q(\Omega))\times L^q(\Gamma)$ 
by verifying asymptotic a priori estimates.

Next, we give an auxiliary theorem to prove the pullback 
$\mathcal{D}$-asymptotical compactness of the process $\{U(t, \tau)\}_{t\geq \tau}$ 
in $(W^{1,p}(\Omega)\cap L^q(\Omega))\times L^q(\Gamma)$.

\begin{theorem}\label{60}
Under assumptions {\rm (H1)--(H2)},
for any $\hat{D}\in\mathcal{D}$ and $t\in\mathbb{R}$, there exists a family 
of positive constants $\{\rho(t):t\in\mathbb{R}\}$ and $\tau_1(t,\hat{D})\leq t$ 
such that
\begin{equation*}
\|u_t(t)\|_{L^2(\Omega)}^2+\|v_t(t)\|_{L^2(\Gamma)}^2\leq \rho(t)
\end{equation*}
for any $(u_\tau,\theta_\tau)\in D(t)$ and $\tau\leq \tau_1(t,\hat{D})$, where
$$
(u_t(s),v_t(s))=\frac{d}{dt}\left(U(t,\tau)(u_\tau,\theta_\tau)\right)\big| _{t=s}
$$ 
and $\rho(t)$ is a positive constant which is independent of the initial data.
\end{theorem}

\begin{proof}
Differentiating \eqref{1} and \eqref{2} with respect to $t$,
and denoting by $\zeta=u_t$, $\eta=v_t$, we obtain
\begin{gather} \label{61}
\begin{aligned}
& \zeta_t-\operatorname{div}(|\nabla u|^{p-2}\nabla \zeta)-(p-2)
\operatorname{div}\left(|\nabla u|^{p-4}(\nabla
u\cdot\nabla \zeta)\nabla u\right)\\
&+(p-1)|u|^{p-2}\zeta+f'(u)\zeta=\frac{dg}{dt},
\end{aligned} \\
\label{62}
\eta_t+(p-2)|\nabla v|^{p-4}(\nabla v\cdot\nabla \eta)\frac{\partial
v}{\partial\nu}+|\nabla v|^{p-2}\frac{\partial
\eta}{\partial\nu}+f'(v)\eta=0,
\end{gather}
where ``$\cdot$'' denotes the dot product in $\mathbb{R}^n$.


Multiplying \eqref{61} by $\zeta$ and integrating over $\Omega$, and
combining \eqref{5} with  \eqref{62}, we obtain
\begin{align*}
&\frac{1}{2}\frac{d}{dt}\left(\|\zeta\|_{L^2(\Omega)}^2+\|
\eta\|_{L^2(\Gamma)}^2\right)+\int_{\Omega}|\nabla
u|^{p-2}|\nabla \zeta|^2\,dx\\
&+(p-2)\int_{\Omega}|\nabla u|^{p-4}(\nabla
u\cdot \nabla \zeta)^2\,dx+(p-1)\int_{\Omega}|u|^{p-2}|\zeta|^2\,dx\\
&\leq l\left(\|\zeta\|_{L^2(\Omega)}^2+\|\eta\|_{L^2(\Gamma)}^2\right)
+\|\frac{dg}{dt}(t)\|_{L^2(\Omega)}\|\zeta\|_{L^2(\Omega)}.
\end{align*}
Integrating \eqref{47} from $t$ to $t+1$ and
using \eqref{48}, we find that
\begin{align*}
&\int_t^{t+1}\|\zeta(s)\|_{L^2(\Omega)}^2\,ds
 +\int_t^{t+1}\|\eta(s)\|_{L^2(\Gamma)}^2\,ds\\
&\leq \mathcal{C}_4(e^{-c_1t}\int_{-\infty}^{t+1}e^{c_1s}
 \| g(s)\|^2_{L^2(\Omega)}\,ds+1)
\end{align*}
uniformly with respect to all initial conditions $(u_\tau,v_\tau)\in D(\tau)$ 
for $\tau\leq\tau_0(t,\hat{D})$, where $\mathcal{C}_4$ is a positive constant.

Therefore, we deduce from the uniform Gronwall inequality that
\begin{align*}
&\|u_t(t+2)\|_{L^2(\Omega)}^2+\|v_t(t+2)\|_{L^2(\Gamma)}^2\\
&\leq \mathcal{C}_5\Big(e^{-c_1t}\int_{-\infty}^{t+1}e^{c_1s}\| g(s)\|^2
_{L^2(\Omega)}\,ds+1+\int_{t-1}^t\|\frac{dg}{dt}(t)\|_{L^2(\Omega)}^2\,ds\Big),
\end{align*}
uniformly with respect to all initial conditions $(u_\tau,v_\tau)\in D(\tau)$ 
for $\tau\leq\tau_0(t,\hat{D})$, where $\mathcal{C}_5$ is a positive constant.
\end{proof}

Next, we prove the process $\{U(t, \tau)\}_{t\geq \tau}$ is pullback 
$\mathcal{D}$-asymptotically compact in 
$(W^{1,p}(\Omega)\cap L^q(\Omega))\times L^q(\Gamma)$.

\begin{theorem}\label{thm62}
 Assume that $f$ and $g$ satisfy conditions {\rm (H1)--(H2)}.
 Then the process $\{U(t, \tau)\}_{t\geq \tau}$
corresponding to problem \eqref{1}-\eqref{4}
is pullback $\mathcal{D}$-asymptotically compact in 
$(W^{1,p}(\Omega)\cap L^q(\Omega))\times L^q(\Gamma)$.
\end{theorem}

\begin{proof}
 Let $B_0=\{B(t):t\in\mathbb{R}\}$ be a pullback $\mathcal{D}$-absorbing 
set in $(W^{1,p}(\Omega)\cap L^q(\Omega))\times L^q(\Gamma)$ obtained in Theorem
\ref{37}, then we need only to show that
 for any $t\in\mathbb{R}$, any $\tau_n\to -\infty$ and 
$(u_{\tau_n},v_{\tau_n})\in B(\tau_n)$, 
$\{(u_n(\tau_n), v_n(\tau_n))\}_{n=0}^{\infty}$ is pre-compact in 
$(W^{1,p}(\Omega)\cap L^q(\Omega))\times L^q(\Gamma)$, where 
$$
(u_n(\tau_n), v_n(\tau_n))=(u(t; \tau_n,(u_{\tau_n},v_{\tau_n})),
v(t;\tau_n, (u_{\tau_n},v_{\tau_n})))=U(t,\tau_n)(u_{\tau_n},v_{\tau_n}).
$$
Note that for Corollary \ref{51}, it remains to prove that
for any $(u_{\tau_n},v_{\tau_n})\in B(\tau_n)$ and $\tau_n\to
 -\infty$, $\{u_n(\tau_n)\}_{n=0}^{\infty}$ is pre-compact in
 $W^{1,p}(\Omega)$.

From Theorem \ref{50} and Corollary \ref{51}, we know that
$\{(u_n(\tau_n), v_n(\tau_n))\}_{n=0}^{\infty}$  is
pre-compact in $L^2(\bar{\Omega},d\mu)$ and
$L^q(\bar{\Omega},d\mu)$. Without loss of
generality, we assume that
$\{(u_n(\tau_n), v_n(\tau_n))\}_{n=0}^{\infty}$ is a Cauchy
sequence in $L^2(\bar{\Omega},d\mu)$ and $L^q(\bar{\Omega},d\mu)$.

 In the following, we prove that $\{u_n(\tau_n)\}_{n=0}^{\infty}$ 
is a Cauchy sequence in $W^{1,p}(\Omega)$.
Then, by simply calculations, we deduce from Lemma \ref{11} that
\begin{align*}
 &\|u_{n_k}(\tau_{n_k})-u_{n_j}(\tau_{n_j})\|_{W^{1,p}(\Omega)}^p\\
 &\leq (-\frac{d}{dt}u_{n_k}(\tau_{n_k})-f(u_{n_k}(\tau_{n_k}))
+\frac{d}{dt}u_{n_j}(\tau_{n_j})
+f(u_{n_j}(\tau_{n_j})),u_{n_k}(\tau_{n_k})-u_{n_j}(\tau_{n_j}))\\
 &\quad +(-\frac{d}{dt}v_{n_k}(\tau_{n_k})-f(v_{n_k}(\tau_{n_k}))
+\frac{d}{dt}v_{n_j}(\tau_{n_j})
+f(v_{n_j}(\tau_{n_j})),v_{n_k}(\tau_{n_k})-v_{n_j}(\tau_{n_j}))\\
 &= I_1+I_2.
\end{align*}
We now estimate separately the two terms $I_{1}$ and
$I_{2}$. By simply calculations and H\"{o}lder's inequality, we
deduce that
\begin{align*}
  I_1
&\leq \| \frac{d}{dt}u_{n_k}(\tau_{n_k})
 -\frac{d}{dt}u_{n_j}(\tau_{n_j})\|_{L^2(\Omega)}
\| u_{n_k}(\tau_{n_k})-u_{n_j}(\tau_{n_j})\|_{L^2(\Omega)}\\
&\quad +C(1+\| u_{n_k}(\tau_{n_k})\|_{L^q(\Omega)}^{q-1}+
\| u_{n_j}(\tau_{n_j})\|_{L^q(\Omega)}^{q-1})\|
u_{n_k}(\tau_{n_k})-u_{n_j}(\tau_{n_j})\|_{L^q(\Omega)}
\end{align*}
and
\begin{align*}
  I_2
&\leq \| \frac{d}{dt}u_{n_k}(\tau_{n_k})-\frac{d}{dt}u_{n_j}
(\tau_{n_j})\|_{L^2(\Gamma)}
\|u_{n_k}(\tau_{n_k})-u_{n_j}(\tau_{n_j})\|_{L^2(\Gamma)}\\
&+C(1+\|u_{n_k}(\tau_{n_k})\|_{L^q(\Gamma)}^{q-1}+
\|u_{n_j}(\tau_{n_j})\|_{L^q(\Gamma)}^{q-1})\|
u_{n_k}(\tau_{n_k})-u_{n_j}(\tau_{n_j})\|_{L^q(\Gamma)}.
\end{align*}
Combining Theorem \ref{50}, Corollary \ref{51} with Theorem \ref{60}, 
yields Theorem \ref{thm62} immediately.
\end{proof}


From Lemma \ref{34} and Theorems \ref{37}, \ref{thm62}, we immediately obtain
Theorem \ref{main}.


\subsection*{Acknowledgments}
The authors would like to thank the anonymous referees for their many helpful 
comments and suggestions. This work was partially supported by the China 
Postdoctoral Science Foundation funded project (No. 2013M532026).


\begin{thebibliography}{00}

\bibitem{ajm} J. M. Arrieta, A. N. Carvalho, A. R. Bernal;
\emph{Parabolic Problems with Nonlinear Boundary Conditions
and Critical Noninearities}, Journal of Differential Equations. 156 (1999)
376-406.

\bibitem{ajm1} J. M. Arrieta, A. N. Carvalho, An\'{i}bal Rodr\'{i}guez-Bernal;
\emph{Attractors of parabolic problems with
nonlinear boundary conditions uniform bounds}, Partial Differential
Equations. 25 (2000) 1-37.

\bibitem{bar} A. R. Bernal;
\emph{Attractors for Parabolic Equations with
Nonlinear Boundary Conditions}, Critical Exponents and Singular
Initial Data, Journal of Differential Equations. 181 (2002) 165-196.

\bibitem{bav} A. V. Babin, M. I. Vishik;
\emph{Attractors of Evolution Equations},
North-Holland, Amsterdam, 1992.

\bibitem{bt} T. Bartsch, Z. Liu;
\emph{On a superlinear elliptic $p$-Laplacian equation}, 
Journal of Differential Equations. 198 (2004) 149-175.

\bibitem{ca} A. Constantin, J. Escher;
\emph{Global existence for fully parabolic
boundary value problems}, Nonlinear Differential Equations and Applications. 13
(2006) 91-118.

\bibitem{ca1} A. Constantin, J. Escher, Z. Yin;
\emph{Global solutions for quasilinear parabolic systems},
Journal of Differential Equations. 197 (2004) 73-84.

\bibitem{dn} D. N. Cheban, P. E. Kloeden, B. Schmalfu{\ss};
\emph{The relationship between pullback, forwards and global attractors of
nonautonomous dynamical systems}, Nonlinear Dynamics and Systems Theory. 2
(2002) 9-28.

\bibitem{hc} H. Crauel, F. Flandoli;
\emph{Attractors for random dynamical systems}, Probability Theory 
and Related Fields. 100 (1994) 365-393.

\bibitem{hc1} H. Crauel, A. Debussche, F. Flandoli;
\emph{Random attractors}, Journal of Dynamics and Differential Equations.
 9 (1997) 307-341.

\bibitem{ct} T. Caraballo, G. {\L}ukasiewicz, J. Real;
\emph{Pullback attractors for asymptotically compact non-autonomous dynamical 
systems}, Nonlinear Analysis. 64 (2006) 484-498.

\bibitem{ct1} T. Caraballo, J. A. Langa, J. Valero;
\emph{The dimension of attractors of non-autonomous partial differential equations}, 
ANZIAM Journal. 45 (2003) 207-222.

\bibitem{cvv} V. V. Chepyzhov, M. I. Vishik;
\emph{Attractors for Equations of Mathematical Physics}, 
American Mathematical Society,  Rhode Isand, 2002.

\bibitem{fzh} Z. H. Fan, C. K. Zhong;
\emph{Attractors for parabolic equations
with dynamic boundary conditions}, Nonlinear Analysis. 68 (2008) 1723-1732.

 \bibitem{gcg} C. G. Gal;
\emph{On a class of degenerate parabolic equations with dynamic
boundary conditions}, Journal of Differential Equations. 253 (2012) 126-166.

\bibitem{pe} P. E. Kloeden, B. Schmalfu{\ss};
\emph{Non-autonomous systems, cocycle attractors and variable time-step
discretization}, Numerical Algorithms. 14 (1997) 141-152.

\bibitem{pe1} P. E. Kloeden, D. J. Stonier;
\emph{Cocycle attractors in nonautonomously perturbed differential equations}, 
Dynamics of Continuous, Discrete and Impulsive Systems. 4 (1998) 211-226.

\bibitem{lk} K. Li, B. You;
\emph{Uniform attractors for the non-autonomous $p$-Laplacian equations 
with dynamic flux boundary conditions}, Boundary Value Problems. 2013, 2013:128.

\bibitem{lss} S. S. Lu, H. Q. Wu, C. K. Zhong;
\emph{Attractors for nonautonomous 2D Navier-Stokes equations with normal 
external force}, Discrete and Continuous Dynamical Systems A. 23 (2005) 701-719.

\bibitem{ly} Y. Li, C. K. Zhong;
\emph{Pullback attractor for the norm-to-weak  continuous process and application 
to the non-autonomous reaction-diffusion  equations},
Applied Mathematics and Computations. 190 (2007) 1020-1029.

\bibitem{ly1} Y. Li, S. Wang, H. Wu;
\emph{Pullback attractors for non-autonomous reaction-diffusion equations 
in $L^p$}, Applied Mathematics and Computation. 207 (2009) 373-379.

\bibitem{pj} J. Petersson;
\emph{A note on quenching for parabolic equations with dynamic boundary conditions},
 Nonlinear Analysis. 58 (2004) 417-423.

\bibitem{pl} L. Popescu, A. R. Bernal;
\emph{On a singularly perturbed wave equation with dynamical boundary
conditions}, Proceedings of the Royal Society of Edinburgh, Section. 
134 (2004) 389-413.

\bibitem{bs} B. Schmalfu{\ss};
\emph{Attractors for non-autonomous dynamical systems}, in: B. Fiedler,
K. Gr\"{o}er, J. Sprekels (Eds.), Proc. Equadiff 99, Berlin, World
Scientific, Singapore, 2000,  684-689.

\bibitem{sht1} H. T. Song, H. Q. Wu;
\emph{Pullback attractors of non-autonomous reaction-diffusion equations}, 
Journal of Mathematical Analysis and Applications. 325 (2007) 1200-1215.

\bibitem{tr}  R. Temam;
\emph{Infinite-dimensional systems in mechanics and physics}, 
Springer-Verlag, New York, 1997.

\bibitem{yb} B. You, C. K. Zhong;
\emph{Global attractors for $p$-Laplacian equations with dynamic
flux boundary conditions}, Advanced Nonlinear Studies. 13 (2013) 391-410.

\bibitem{yl} L. Yang;
\emph{Uniform attractors for the closed process and
applications to the reaction-diffusion with dynamical boundary
condition}, Nonlinear Analysis. 71 (2009) 4012-4025.

\bibitem{yl1} L. Yang, M. H. Yang, P. E. Kloeden;
\emph{Pullback attractors for non-autonomous quasi-linear parabolic 
equations with a dynamical boundary condition}, 
Discrete and Continuous Dynamical Systems B. 17 (2012) 2635-2651.

\bibitem{yw} Y. Wang, C. Zhong;
\emph{Pullback D-attractors for nonautonomous sine-Gordon
equations}, Nonlinear Analysis. 67 (2007) 2137-2148.

\end{thebibliography}

\end{document}
