Existence of Solutions to a Nonlinear Parabolic Equation of Fourth-Order in Variable Exponent Spaces

Abstract

This paper is devoted to studying the existence and uniqueness of weak solutions for an initial boundary problem of a nonlinear fourth-order parabolic equation with variable exponent $v_t+{\mathrm{div}(|\nabla ▵v|}^{p(x)-2}{\nabla ▵v)-|▵v|}^{q(x)-2}▵v=g(x,v)$. By applying Leray-Schauder’s fixed point theorem, the existence of weak solutions of the elliptic problem is given. Furthermore, the semi-discrete method yields the existence of weak solutions of the corresponding parabolic problem by constructing two approximate solutions.

1. Introduction

We mainly study the following fourth-order parabolic equations with variable exponents:

$$\begin{array}{cccc}{v}_t+{\mathrm{div}(|\nabla ▵v|}^{p(x)-2}{\nabla ▵v)-|▵v|}^{q(x)-2}▵v=g(x,v),\hfill & & & \hfill (x,t)\in {\mathrm{\Omega }}_T,\end{array}$$

(1)

$$\begin{array}{cccc}v(x,t)=▵v(x,t)=0,\hfill & & & \hfill (x,t)\in {\mathrm{\Gamma }}_T,\end{array}$$

(2)

$$\begin{array}{cccc}v(x,0)=v_0(x),\hfill & & & \hfill x\in \mathrm{\Omega },\end{array}$$

(3)

where Ω is an open, bounded domain in ${\mathbb{R}}^N$, $\partial \mathrm{\Omega }\in C^1$. Define ${\mathrm{\Omega }}_T=\mathrm{\Omega }\times (0,T)$ and ${\mathrm{\Gamma }}_T=\partial \mathrm{\Omega }\times (0,T)$. If p is a constant (especially $p\equiv 2$ and $q\equiv 2$), the Equation (1) has the structure of the classical Cahn–Hilliard problem, which is often used to describe the evolution of a conserved concentration field during phase separation in physics. It is also related to the thin-film equation if ${|\nabla ▵v|}^{p(x)-2}$ becomes $v^p$, which can analyze the motion of a very thin layer of viscous incompressible fluids along an include plane.

There have been some results related to the existence, uniqueness and properties of solutions to the fourth-order degenerate parabolic equations (see [1,2]). The paper [3] has studied the existence of the Cahn–Hilliard equation and the reader may refer to [4] to obtain its physical background. For the constant exponent case of (1), the paper [5] has given the existence and uniqueness of weak solutions. For the problems in variable exponent spaces, the papers [6,7,8] have studied the existence of some fourth-order parabolic equations with a variable exponent, and [9] has given the Fujita type conditions for fast diffusion equation.

For the research of the existence and long-time behavior of the fourth-order partial differential equations, the entropy functional method is often applied in order to obtain the necessary estimates and to show the entropy dissipation. The large time behavior of solutions of the thin film equation $u_t+{(u^n{u}_{xxx})}_x=0$ was addressed in [10,11] by the entropy function method. For $0<n<3$, [12] proved the existence of (1) in the distributional sense and obtained the exponentially fast convergence in $L^{\infty }$-norm via the entropy method of a regularized problem. We apply the idea of the entropy method to deal with the corresponding problems with variable exponents.

In this paper, we apply the Leray-Schauder’s fixed point theorem to prove the existence of weak solutions of the corresponding elliptic problem of (1)–(3) in order to deal with the nonlinear source. Furthermore, the semi-discrete method yields the existence of weak solutions of the parabolic problem by constructing two approximate solutions. We will show the effect of the variable exponents and the second-order nonlinear diffusion to the degenerate parabolic Equation (1).

1.1. Preliminaries

We introduce some elementary concepts and lemmas related to the variable exponent spaces in this part.

Let $p(x)\ge 1$ be a continuous function in $\overline{\mathrm{\Omega }}$ and we define the variable exponent space as follows:

$${\mathcal{L}}^{p(x)}(\mathrm{\Omega })=\left\{v(x):v\text{ is measureable and }{\mathbb{T}}_{p(.)}(v)={\int }_{\mathrm{\Omega }}{\left|v\right|}^{p(x)}\mathrm{d}x<\infty \right\}$$

with the norm

$${\parallel v\parallel }_{p(x)}=inf\left\{\gamma >0:{\mathbb{T}}_{p(.)}\left(\frac{v}{\gamma }\right)\le 1\right\}.$$

It is easy to check that the variable exponent space ${\mathcal{L}}^{p(x)}(\mathrm{\Omega })$ becomes the classical Lebesgue space ${\mathcal{L}}^p(\mathrm{\Omega })$ when $p(x)$ is a positive constant.

For convenience, we list some definitions and notations of the generalized Lebesgue–Sobolev space ${\mathcal{W}}^{k,p(x)}(\mathrm{\Omega })$:

$$\begin{array}{c}{\mathcal{W}}^{k,p(x)}(\mathrm{\Omega })=:\left\{v(x)\in {\mathcal{L}}^{p(x)}(\mathrm{\Omega }):D^{\alpha }v\in {\mathcal{L}}^{p(x)}(\mathrm{\Omega }),\left|\alpha \right|\le k\right\},\hfill \\ {\parallel v\parallel }_{{\mathcal{W}}^{k,p(x)}}=:\sum _{\left|\alpha \right|\le k}{\left|D^{\alpha }v\right|}_{p(x)},\hfill \\ E_1=:\{v\in {\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })\cap {\mathcal{W}}^{2,p(x)}(\mathrm{\Omega })\cap {\mathcal{W}}^{2,q(x)}(\mathrm{\Omega })|▵v\in {\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })\},\hfill \\ E_2=:\{v\in {\mathcal{H}}_0^1(\mathrm{\Omega })\cap {\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })\cap {\mathcal{W}}^{2,p(x)}(\mathrm{\Omega })\cap {\mathcal{W}}^{2,q(x)}(\mathrm{\Omega })|▵v\in {\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })\},\hfill \\ {\mathcal{L}}^{p^{\prime }(x)}(\mathrm{\Omega })\text{ denotes the dual space with}\frac{1}{p^{\prime }(x)}+\frac{1}{p(x)}=1.\hfill \end{array}$$

Moreover, ${\mathcal{W}}_0^{k,p(x)}$ denotes the closure of $C_0^{\infty }(\mathrm{\Omega })$ in ${\mathcal{W}}^{k,p(x)}(\mathrm{\Omega })$-norm, ${\mathcal{W}}_0^{-1,p^{\prime }(x)}(\mathrm{\Omega })$ denotes the dual space of ${\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })$. For any positive continuous function $\theta (x)$, we define

$${\theta }^{-}=\underset{x\in \overline{\mathrm{\Omega }}}{inf}\theta (x),{\theta }^{+}=\underset{x\in \overline{\mathrm{\Omega }}}{sup}\theta (x).$$

Throughout the paper, C and $C_i(i=1,2,3,...)$ denote the general positive constants independent of solutions and may change from line to line.

In the following, we list some known results for the variable exponent spaces (see [13,14]).

Lemma 1.

Letting $f\in {\mathcal{L}}^{p(x)}(\mathrm{\Omega })$, one has

$$\begin{array}{cc}\hfill (1)& {\parallel f\parallel }_{p(x)}<1(=1;>1)⟺{\mathbb{T}}_{p(.)}(f)<1(=1;>1);\hfill \\ \hfill (2)& {\parallel f\parallel }_{p(x)}<1⟹{\parallel f\parallel }_{p(x)}^{p^{+}}\le {\mathbb{T}}_{p(.)}(f)\le {\parallel f\parallel }_{p(x)}^{p^{-}};\hfill \\ & {\parallel f\parallel }_{p(x)}\ge 1⟹{\parallel f\parallel }_{p(x)}^{p^{-}}\le {\mathbb{T}}_{p(.)}(f)\le {\parallel f\parallel }_{p(x)}^{p^{+}};\hfill \\ \hfill (3)& {\parallel f\parallel }_{p(x)}\to 0⟺{\mathbb{T}}_{p(.)}(f)\to 0;{\parallel f\parallel }_{p(x)}\to \infty ⟺{\mathbb{T}}_{p(.)}(f)\to \infty .\hfill \end{array}$$

Lemma 2.

(Poincaré’s inequality) Letting $f\in {\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })$, there exists a positive constant C such that ${\parallel f\parallel }_{p(x)}\le C{\parallel \nabla f\parallel }_{p(x)}.$

Lemma 3.

(Hölder’s inequality) Letting $f\in {\mathcal{L}}^{p(x)}(\mathrm{\Omega })$ and $g\in {\mathcal{L}}^{p^{\prime }(x)}(\mathrm{\Omega })$, one has $\left|{\int }_{\mathrm{\Omega }}fg\mathrm{d}x\right|\le \left(\frac{1}{p^{-}}+\frac{1}{p^{\prime -}}\right){\parallel f\parallel }_{p(x)}{\parallel g\parallel }_{p^{\prime }(x)}\le {2\parallel f\parallel }_{p(x)}{\parallel g\parallel }_{p^{\prime }(x)}.$

1.2. Results

In (1), we require that $p(x)$ and $q(x)$ are two continuous functions in $\overline{\mathrm{\Omega }}$ and $p^{-},q^{-}>1$. Besides, the nonlinear source term $g(x,v)\in C^1(\mathrm{\Omega }\times R)$ satisfies the growth condition:

$$\begin{array}{ccc}\hfill \left|g(x,v)\right|\le {K\left|v\right|}^{l(x)}+s(x),& & \hfill v\in (-\infty ,+\infty ),x\in \mathrm{\Omega },\end{array}$$

(4)

where K is a positive constant, $l(x)$ is a continuous function in $\overline{\mathrm{\Omega }}$ and $s(x)\in {\mathcal{L}}^{p^{\prime }(x)}(\mathrm{\Omega })$. Furthermore, by letting $\pi (x)=:l(x)p^{\prime }(x)$, we require that

$$\frac{{\pi }^{+}}{p^{-}}<1,0\le l(x)\le \frac{Np(x)}{(N-p(x))p^{\prime }(x)}.$$

(5)

The corresponding steady-state problem of (1)–(3) has the form:

$$\begin{array}{c}{\mathrm{div}(|\nabla ▵v|}^{p(x)-2}{\nabla ▵v)-|▵v|}^{q(x)-2}▵v=g(x,v)\mathrm{in}\mathrm{\Omega },\hfill \end{array}$$

(6)

$$\begin{array}{c}v=▵v=0\mathrm{on}\partial \mathrm{\Omega }.\hfill \end{array}$$

(7)

The weak solution is defined in the following sense.

Definition 1.

A function $v\in E_1$ is said to be a weak solution of (6) and (7) provided that

$$\begin{array}{c}\hfill -{\int }_{\mathrm{\Omega }}{|\nabla ▵v|}^{p(x)-2}\nabla ▵v\nabla \varphi \mathrm{d}x-{\int }_{\mathrm{\Omega }}{|▵v|}^{q(x)-2}▵v\varphi \mathrm{d}x={\int }_{\mathrm{\Omega }}g(x,v)\varphi \mathrm{d}x\end{array}$$

(8)

for each $\varphi \in {\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })\cap {\mathcal{L}}^{q(x)}(\mathrm{\Omega })$ and $g(x,v(x))\in {\mathcal{L}}^{p^{\prime }(x)}(\mathrm{\Omega })$.

The following theorem gives the existence of solutions.

Theorem 1.

Let $v_0\in {\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })$. There exists at least a weak solution of (6) and (7) satisfying Definition 1.

For the evolution equation case, we define the weak solution of (1)–(3) as following.

Definition 2.

A function v is said to be a weak solution of (1)–(3) provided that

(i)

$v\in C([0,T];{\mathcal{L}}^2({\mathrm{\Omega }}_T))\cap {\mathcal{L}}^{p^{-}}(0,T;{\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega }))\cap {\mathcal{L}}^{q^{-}}(0,T;{\mathcal{W}}^{2,q(x)}(\mathrm{\Omega }))$, $▵v\in {\mathcal{L}}^{p^{\prime -}}(0,T;{\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega }))$, $\frac{\partial v}{\partial t}\in {\mathcal{L}}^{p^{\prime -}}(0,T;{\mathcal{W}}^{-1,p(x)}(\mathrm{\Omega }))$, $v(x,0)=v_0(x)$ a.e. in Ω;

(ii)

For any $\varphi \in {\mathcal{L}}^{p^{+}}(0,T;{\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega }))\cap {\mathcal{L}}^{q^{+}}(0,T;{\mathcal{W}}^{2,q(x)}(\mathrm{\Omega }))$, one has

$$\begin{array}{c}{\int }_0^T〈\frac{\partial v}{\partial t},\varphi 〉\mathrm{d}t={\int }_0^T{\int }_{\mathrm{\Omega }}{|\nabla ▵v|}^{p(x)-2}\nabla ▵v\nabla \varphi \mathrm{d}x\mathrm{d}t\hfill \\ +{\int }_0^T{\int }_{\mathrm{\Omega }}{|▵v|}^{q(x)-2}▵v\varphi \mathrm{d}x\mathrm{d}t+{\int }_0^T{\int }_{\mathrm{\Omega }}g(x,v)\varphi \mathrm{d}x\mathrm{d}t.\hfill \end{array}$$

The existence of solutions is the following theorem.

Theorem 2.

Let $p^{\prime -}>1$, $q^{\prime -}>1$, $v_0\in {\mathcal{H}}_0^1(\mathrm{\Omega })\cap {\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })$ and $▵v_0\in {\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })$. There exists at least a weak solution of (1)–(3).

Moreover, the solution of (1)–(3) is unique when $g(x,v)=\mu (v(x,t)-b(x))$ where μ is a constant and $b(x)\in {\mathcal{L}}^{p^{\prime }(x)}(\mathrm{\Omega })$.

This paper is organized as follows. In Section 2, we prove the existence and uniqueness of weak solution to the steady-state problem by using Leray-Schauder’s fixed point theorem. In Section 3, we prove the existence of the solution to an evolution equation by applying the semi-discrete method with necessary uniform estimates.

2. Steady-State Problem

In order to apply the fixed point theorem, we consider a steady-state problem with the source $g(x)$:

$$\begin{array}{c}{\mathrm{div}(|\nabla ▵v|}^{p(x)-2}{\nabla ▵v)-|▵v|}^{q(x)-2}▵v=g(x)\mathrm{in}\mathrm{\Omega },\hfill \end{array}$$

(9)

$$\begin{array}{c}v=▵v=0\mathrm{on}\partial \mathrm{\Omega }.\hfill \end{array}$$

(10)

By constructing an energy functional and obtaining its minimizer, we have the following existence of weak solutions.

Lemma 4.

Let $g\in {\mathcal{L}}^{p^{\prime }(x)}(\mathrm{\Omega })$. There exists a unique weak solution $v\in E_1$ of (9) and (10) satisfying

$$\begin{array}{cc}& -{\int }_{\mathrm{\Omega }}{|\nabla ▵v|}^{p(x)-2}\nabla ▵v\nabla ▵\varphi \mathrm{d}x-{\int }_{\mathrm{\Omega }}{|▵v|}^{q(x)-2}▵v▵\varphi \mathrm{d}x\hfill \\ \hfill =& {\int }_{\mathrm{\Omega }}g(x)▵\varphi \mathrm{d}x\hfill \end{array}$$

(11)

for any $\varphi \in E_1$.

Proof. 

Introduce a functional

$$\begin{array}{c}\hfill ϝ(v)={\int }_{\mathrm{\Omega }}\frac{1}{p(x)}{|\nabla ▵v|}^{p(x)}\mathrm{d}x+{\int }_{\mathrm{\Omega }}\frac{1}{q(x)}{|▵v|}^{q(x)}\mathrm{d}x-{\int }_{\mathrm{\Omega }}g▵v\mathrm{d}x.\end{array}$$

(12)

For the last term, Hölder’s inequality, the Young inequality, the Sobolev embedding theorem (see [15]) and the ${\mathcal{L}}^p$-theory of the second-order elliptic equation (see [16]) gives

$$\begin{array}{cc}\hfill |{\int }_{\mathrm{\Omega }}g▵v\mathrm{d}x|& \le {\parallel g\parallel }_{p^{\prime }(x)}{\parallel ▵v\parallel }_{p(x))}\hfill \\ & \le C\left(\frac{1}{\epsilon }\right)max\left\{{\parallel g\parallel }_{p^{\prime }(x)}^{p^{{}^{\prime }-}},{\parallel g\parallel }_{p^{\prime }(x)}^{p^{{}^{\prime }+}}\right\}+\epsilon {\int }_{\mathrm{\Omega }}{|▵v|}^{p(x)}\mathrm{d}x\hfill \\ & \le Cmax\left\{{\parallel g\parallel }_{p^{\prime }(x)}^{p^{{}^{\prime }-}},{\parallel g\parallel }_{p^{\prime }(x)}^{p^{{}^{\prime }+}}\right\}+\frac{1}{2}{\int }_{\mathrm{\Omega }}\frac{1}{p(x)}{|\nabla ▵v|}^{p(x)}\mathrm{d}x.\hfill \end{array}$$

(13)

On the other hand, (12) implies

$$\begin{array}{c}\hfill -C\le \underset{\nu \in E_1}{inf}ϝ(\nu )\le ϝ(0)=0.\end{array}$$

Hence there exists a sequence ${\left\{v_k\right\}}_{k=1}^{\infty }\in E_1$ such that

$$\begin{array}{c}\hfill ϝ(v_k)\to \underset{\nu \in E_1}{inf}ϝ(\nu )(k\to \infty ).\end{array}$$

(14)

Equations (13) and (14) give

$$\begin{array}{c}\hfill \parallel \nabla ▵v_k{\parallel }_{p(x)}\le C,{\parallel ▵v_k\parallel }_{q(x)}\le C,\end{array}$$

which implies that $ϝ(v_k)$ is bounded and thus Lemmas 1–3 yield

$$\begin{array}{c}{\parallel v_k\parallel }_{2,p(x)}\le \parallel ▵v_k{\parallel }_{p(x)}\le C{\parallel \nabla ▵v_k\parallel }_{p(x)}\hfill \\ \le C_{1}max\left\{{\left({\int }_{\mathrm{\Omega }}{|\nabla ▵v_k|}^{p(x)}\right)}^{\frac{1}{p^{-}}},{\left({\int }_{\mathrm{\Omega }}{|\nabla ▵v_k|}^{p(x)}\right)}^{\frac{1}{p^{+}}}\right\}\le C_2,\hfill \end{array}$$

and

$$\begin{array}{c}{\parallel v_k\parallel }_{2,q(x)}\le {\parallel ▵v_k\parallel }_{q(x)}\hfill \\ \le C_{3}max\left\{{\left({\int }_{\mathrm{\Omega }}{|▵v_k|}^{q(x)}\right)}^{\frac{1}{q^{-}}},{\left({\int }_{\mathrm{\Omega }}{|▵v_k|}^{q(x)}\right)}^{\frac{1}{q^{+}}}\right\}\le C_4.\hfill \end{array}$$

It shows that $v_k$ belongs to the space ${\mathcal{W}}_0^{1,p(x)}\bigcap {\mathcal{W}}^{2,p(x)}\bigcap {\mathcal{W}}^{2,q(x)}$ uniformly, and then there exists a function $v\in E_1$ such that

$$\begin{array}{cc}\hfill v_k& \rightharpoonup v\text{ weakly in }{\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })\cap {\mathcal{W}}^{2,p(x)}(\mathrm{\Omega })\cap {\mathcal{W}}^{2,q(x)}(\mathrm{\Omega }),\hfill \\ \hfill ▵v_k& \rightharpoonup ▵v\text{ weakly in }{\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })\cap {\mathcal{L}}^{q(x)}(\mathrm{\Omega }).\hfill \end{array}$$

Furthermore, since $ϝ(v)$ is weakly lower semi-continuous on $E_1$, we have

$$\begin{array}{c}\hfill \underset{\nu \in E_1}{inf}ϝ(\nu )\le ϝ(v)\le \underset{k\to \infty }{lim}infϝ(v_k)=\underset{\nu \in E_1}{inf}ϝ(\nu ),\end{array}$$

i.e., v is a minimizer of $ϝ(·)$ and $ϝ(v)={inf}_{\nu \in E_1}ϝ(\nu ).$ It guarantees that v is a weak solution of (9) and (10).

The uniqueness is obvious and we omit the details. ☐

Now, we consider the problem (6) and (7) with the nonlinear source $g(x,v)$.

Lemma 5.

Letting $v(x)\in E_1$ be a weak solution of (6) and (7), one has ${\parallel v\parallel }_{E_1}\le C$.

Proof. 

Multiplying (8) by v gives

$$\begin{array}{cc}& {\int }_{\mathrm{\Omega }}{|\nabla ▵v|}^{p(x)}\mathrm{d}x+{\int }_{\mathrm{\Omega }}{|▵v|}^{q(x)}\mathrm{d}x\hfill \\ \hfill =& -{\int }_{\mathrm{\Omega }}g(x,v)▵v\mathrm{d}x\le {2\parallel g\parallel }_{p^{\prime }(x)}{\parallel ▵v\parallel }_{p(x)}\hfill \\ \hfill \le & {C\parallel g\parallel }_{p^{\prime }(x)}{\parallel \nabla ▵v\parallel }_{p(x)}\hfill \\ \hfill \le & \frac{1}{4}{\int }_{\mathrm{\Omega }}{|\nabla ▵v|}^{p(x)}\mathrm{d}x+{Cmax\{\parallel g\parallel }_{p^{\prime }(x)}^{p^{\prime -}}{,\parallel g\parallel }_{p^{\prime }(x)}^{p^{\prime +}}\}.\hfill \end{array}$$

(15)

By Lemmas 2 and 3 and ${\mathcal{L}}^p$-estimate (see [16]), we conclude that

$$\begin{array}{cc}\hfill {\int }_{\mathrm{\Omega }}{\left|g(x,v)\right|}^{p^{\prime }(x)}\mathrm{d}x& \le C{\int }_{\mathrm{\Omega }}{(K|v|}^{l(x)}{+s(x))}^{p^{\prime }(x)}\mathrm{d}x\hfill \\ & \le C{M}^{p^{\prime +}}{\int }_{\mathrm{\Omega }}{\left|v\right|}^{l(x)p^{\prime }(x)}\mathrm{d}x+C{\int }_{\mathrm{\Omega }}{\left|s(x)\right|}^{p^{\prime (x)}}\mathrm{d}x+C\hfill \\ & \le {C\parallel v\parallel }_{{\mathcal{W}}^{1,p(x)}}^{l(x)p^{\prime }(x)}+C\le C{\parallel ▵v\parallel }_{p(x)}^{l(x)p^{\prime }(x)}+C\hfill \\ & \le {C\parallel \nabla ▵v\parallel }_{p(x)}^{l(x)p^{\prime }(x)}+C,\hfill \end{array}$$

(16)

and thus

$$\begin{array}{cc}\hfill max\{\parallel g{\parallel }_{p^{\prime }(x)}^{p^{\prime -}},\parallel g{\parallel }_{p^{\prime }(x)}^{p^{\prime +}}\}& \le Cmax\{\parallel \nabla ▵v{\parallel }_{p(x)}^{{\pi }^{+}},\parallel \nabla ▵v{\parallel }_{p(x)}^{{\pi }^{-}}\}+C\hfill \\ & \le \frac{1}{4}{\int }_{\mathrm{\Omega }}{|\nabla ▵v|}^{p(x)}\mathrm{d}x+C.\hfill \end{array}$$

(17)

Equations (15)–(17) yield

$$\begin{array}{c}\hfill {\int }_{\mathrm{\Omega }}{|\nabla ▵v|}^{p(x)}\mathrm{d}x+{\int }_{\mathrm{\Omega }}{|▵v|}^{q(x)}\mathrm{d}x\le C.\end{array}$$

(18)

It completes the proof of Lemma 5. ☐

Proof of Theorem 1.

Letting $\omega \in {\mathcal{L}}^{p^{*}(x)}(\mathrm{\Omega })$ and $\delta \in [0,1]$ where we choose $p^{*}(x)$ such that $E_1\hookrightarrow {\mathcal{L}}^{p^{*}(x)}(\mathrm{\Omega })$ is compact, we consider the auxiliary problem

$$\begin{array}{c}{\mathrm{div}(|\nabla ▵v|}^{p(x)-2}{\nabla ▵v)-|▵v|}^{q(x)-2}▵v=\delta g(x,\omega )\mathrm{in}\mathrm{\Omega },\hfill \\ v=▵v=0\mathrm{on}\partial \mathrm{\Omega }.\hfill \end{array}$$

Lemma 4 ensures its existence and so we can define the fixed point operator

$$\begin{array}{cc}\hfill T:[0,1]\times {\mathcal{L}}^{p^{*}(x)}(\mathrm{\Omega })& ⟶{\mathcal{L}}^{p^{*}(x)}(\mathrm{\Omega }),\hfill \\ \hfill (\delta ,\omega )& ⟼v\hfill \end{array}$$

and $T(\omega ,0)=0$.

If $\omega \in {\mathcal{L}}^{p^{*}(x)}(\mathrm{\Omega })$ satisfies $T(\omega ,\delta )=\omega$, we can check that ${\parallel \omega \parallel }_{E_1}\le C$ where $C>0$ is independence of ω and δ from the idea of Lemma 5. The compact embedding $E_1\hookrightarrow {\mathcal{L}}^{p^{*}(x)}(\mathrm{\Omega })$ can ensure that T is a continuous and compact operator. Leray-Schauder’s fixed point theorem yields the existence of solutions of (6) and (7). ☐

3. Evolution Equation

In this section, we study the existence solutions of (1)–(3). For this purpose, we establish a semi-discrete problem at first:

$$\begin{array}{cc}& \frac{1}{h}(v_k-v_{k-1})+\mathrm{div}(|\nabla ▵v_k{|}^{p(x)-2}\nabla ▵v_k)-|▵v_k{|}^{q(x)-2}▵v_k\hfill \\ \hfill =& g(x,v_{k-1}),x\in \mathrm{\Omega },\hfill \end{array}$$

(19)

$$\begin{array}{cc}& v_k(x,t)=▵v_k(x,t)=0,x\in \partial \mathrm{\Omega },\hfill \end{array}$$

(20)

where $v_k=v(x,kh)$, $h=\frac{T}{n},k=1,2,\cdots ,n$ and $n\in \mathbb{N}$.

Lemma 6.

Assume $v_0\in E_2$. (19) and (20) admits a unique weak solution $v_k\in E_2$ satisfying

$$\begin{array}{cc}& \sum _{k=1}^i{\int }_{\mathrm{\Omega }}|\nabla v_k{|}^2\mathrm{d}x+h\sum _{k=1}^i{\int }_{\mathrm{\Omega }}|\nabla ▵v_k{|}^{p(x)}\mathrm{d}x+h\sum _{k=1}^i{\int }_{\mathrm{\Omega }}{|▵v_k|}^{q(x)}\mathrm{d}x\hfill \\ \hfill \le & \sum _{k=1}^i{\int }_{\mathrm{\Omega }}|\nabla v_{k-1}{|}^2\mathrm{d}x+\frac{h}{2}\sum _{k=1}^i{\int }_{\mathrm{\Omega }}{|\nabla ▵v_{k-1}|}^{p(x)}\mathrm{d}x+CT,\hfill \end{array}$$

(21)

and

$$\begin{array}{cc}& {\int }_{\mathrm{\Omega }}{|\nabla v_i|}^2\mathrm{d}x+\frac{h}{2}\sum _{k=1}^i{\int }_{\mathrm{\Omega }}{|\nabla ▵v_i|}^{p(x)}\mathrm{d}x+\frac{h}{2}{\int }_{\mathrm{\Omega }}{|\nabla ▵v_i|}^{p(x)}\mathrm{d}x\hfill \\ & +h\sum _{k=1}^i{\int }_{\mathrm{\Omega }}{|▵v_k|}^{q(x)}\mathrm{d}x\hfill \\ & \le {\int }_{\mathrm{\Omega }}|\nabla v_0{|}^2\mathrm{d}x+\frac{h}{2}{\int }_{\mathrm{\Omega }}{|\nabla ▵v_0|}^{p(x)}\mathrm{d}x+CT.\hfill \end{array}$$

(22)

Proof. 

According to the argument of the Section 2, we conclude that the problem (19) and (20) has a unique weak solution $v_k\in E_2$ satisfying

$$\begin{array}{c}\frac{1}{h}{\int }_{\mathrm{\Omega }}\nabla (v_k-v_{k-1})·\nabla \varphi \mathrm{d}x+{\int }_{\mathrm{\Omega }}{|\nabla ▵v_k|}^{p(x)-2}|\nabla ▵v_k|·\nabla ▵\varphi \mathrm{d}x\hfill \\ +{\int }_{\mathrm{\Omega }}{|▵v_k|}^{q(x)-2}▵v_k▵\varphi \mathrm{d}x=-{\int }_{\mathrm{\Omega }}g(x,v_{k-1})▵\varphi \mathrm{d}x\hfill \end{array}$$

(23)

for any $\varphi \in C_0^{\infty }(\mathrm{\Omega })$. Letting $\varphi =v_k$ in (23), we have

$$\begin{array}{cc}& \frac{1}{h}{\int }_{\mathrm{\Omega }}{|\nabla v_k|}^2\mathrm{d}x+{\int }_{\mathrm{\Omega }}{|\nabla ▵v_k|}^{p(x)}\mathrm{d}x+{\int }_{\mathrm{\Omega }}{|▵v_k|}^{q(x)}\mathrm{d}x\hfill \\ \hfill =& \frac{1}{h}{\int }_{\mathrm{\Omega }}\nabla v_k·\nabla v_{k-1}\mathrm{d}x-{\int }_{\mathrm{\Omega }}g(x,v_{k-1})▵v_k\mathrm{d}x\hfill \\ \hfill \le & \frac{1}{2h}{\int }_{\mathrm{\Omega }}{|\nabla v_{k-1}|}^2\mathrm{d}x+\frac{1}{2h}{\int }_{\mathrm{\Omega }}{|\nabla v_k|}^2\mathrm{d}x-{\int }_{\mathrm{\Omega }}g(x,v_{k-1})▵v_k\mathrm{d}x.\hfill \end{array}$$

(24)

Similar to the proof of (16) and (17), we get

$$\begin{array}{cc}& -{\int }_{\mathrm{\Omega }}g(x,v_{k-1})▵v_k\mathrm{d}x\hfill \\ \hfill \le & \frac{1}{2}{\int }_{\mathrm{\Omega }}|\nabla ▵v_k{|}^{p(x)}\mathrm{d}x+C{\int }_{\mathrm{\Omega }}{|g(x,v_{k-1})|}^{p^{\prime }(x)}\mathrm{d}x\hfill \\ \hfill \le & \frac{1}{2}{\int }_{\mathrm{\Omega }}|\nabla ▵v_k{|}^{p(x)}\mathrm{d}x+\frac{1}{4}{\int }_{\mathrm{\Omega }}{|\nabla ▵v_{k-1}|}^{p(x)}\mathrm{d}x+C.\hfill \end{array}$$

(25)

By (24) and (25), one has

$$\begin{array}{cc}& \frac{1}{2h}{\int }_{\mathrm{\Omega }}|\nabla v_k{|}^2\mathrm{d}x+\frac{1}{2}{\int }_{\mathrm{\Omega }}|\nabla ▵v_k{|}^{p(x)}\mathrm{d}x+{\int }_{\mathrm{\Omega }}{|▵v_k|}^{q(x)}\mathrm{d}x\hfill \\ \hfill \le & \frac{1}{2h}{\int }_{\mathrm{\Omega }}|\nabla v_{k-1}{|}^2\mathrm{d}x+\frac{1}{4}{\int }_{\mathrm{\Omega }}{|\nabla ▵v_{k-1}|}^{p(x)}\mathrm{d}x+C.\hfill \end{array}$$

(26)

Hence, for any $1\le i\le n$, we obtain

$$\begin{array}{cc}& \sum _{k=1}^i{\int }_{\mathrm{\Omega }}|\nabla v_k{|}^2\mathrm{d}x+h\sum _{k=1}^i{\int }_{\mathrm{\Omega }}|\nabla ▵v_k{|}^{p(x)}\mathrm{d}x+h\sum _{k=1}^i{\int }_{\mathrm{\Omega }}{|▵v_k|}^{q(x)}\mathrm{d}x\hfill \\ \hfill \le & \sum _{k=1}^i{\int }_{\mathrm{\Omega }}|\nabla v_{k-1}{|}^2\mathrm{d}x+\frac{h}{2}\sum _{k=1}^i{\int }_{\mathrm{\Omega }}{|\nabla ▵v_{k-1}|}^{p(x)}\mathrm{d}x+CT.\hfill \end{array}$$

(27)

It completes the proof of (21) and (22) obtained from (21). ☐

Now, we are in the position to define the first approximate solution of (1)–(3)

$$\begin{array}{cc}\hfill {\varpi }^{(n)}(x,t)=& \sum _{k=1}^n{\chi }_k(t)v_k(t)\hfill \end{array}$$

(28)

where ${\chi }_k(t)$ is the characteristic function over the interval $((k-1)h,kh]$ for $k=1,2,...,n$. For this approximate solution, we have the following uniform estimates.

Lemma 7.

One has

$$\begin{array}{c}\parallel {\varpi }^{(n)}{\parallel }_{{\mathcal{L}}^{\infty }(0,T;{\mathcal{H}}_0^1(\mathrm{\Omega }))}+{\parallel ▵{\varpi }^{(n)}\parallel }_{{\mathcal{L}}^{p^{-}}(0,T;{\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega }))}\hfill \\ +\parallel |\nabla ▵{\varpi }^{(n)}{|}^{p(x)-2}\nabla ▵{\varpi }^{(n)}{\parallel }_{{\mathcal{L}}^{p\prime -}({\mathrm{\Omega }}_T)}+\parallel |▵{\varpi }^{(n)}{{|}^{q(x)-2}▵{\varpi }^{(n)}\parallel }_{{\mathcal{L}}^{q^{\prime }}({\mathrm{\Omega }}_T)}\hfill \\ +\parallel g(·,{\varpi }^{(n)}){\parallel }_{{\mathcal{L}}^{p^{\prime -}}({\mathrm{\Omega }}_T)}\le C.\hfill \end{array}$$

(29)

Proof. 

By Lemma 6 and

$$\parallel \nabla {\varpi }^{(n)}{\parallel }_{{\mathcal{L}}^2(\mathrm{\Omega })}^2={\parallel \nabla v_k(x)\parallel }_{{\mathcal{L}}^2(\mathrm{\Omega })}^2\le C,$$

(30)

we have the estimate

$$\parallel \nabla {\varpi }^{(n)}{\parallel }_{{\mathcal{L}}^{\infty }(0,T;{\mathcal{L}}^2(\mathrm{\Omega }))}\le C.$$

(31)

On the other hand, we have

$$\begin{array}{c}\parallel ▵{\varpi }^{(n)}{\parallel }_{{\mathcal{L}}^{p^{-}}(0,T;{\mathcal{W}}_0^{1,p(x)})}={\left({\int }_0^T{\parallel ▵{\varpi }^{(n)}\parallel }_{1,p(x)}^{p^{-}}\mathrm{d}x\right)}^{\frac{1}{p^{-}}}\hfill \\ ={\left(h\sum _{k=1}^n{\parallel ▵v_k\parallel }_{1,p(x)}^{p^{-}}\mathrm{d}x\right)}^{\frac{1}{p^{-}}}\hfill \\ \le {\left(Ch\sum _{k=1}^n{\parallel \nabla ▵v_k\parallel }_{p(x)}^{p^{-}}\mathrm{d}x\right)}^{\frac{1}{p^{-}}}\le C{\left(h\sum _{k=1}^n({\int }_{\mathrm{\Omega }}|\nabla ▵v_k{|}^{p(x)}\mathrm{d}x+1)\right)}^{\frac{1}{p^{-}}}\hfill \\ \le C{(C+T)}^{\frac{1}{p^{-}}}.\hfill \end{array}$$

(32)

Letting $i=n$ in (22), we get

$$\begin{array}{cc}& {\int }_0^T{\int }_{\mathrm{\Omega }}|\nabla ▵{\varpi }^{(n)}{|}^{p(x)}\mathrm{d}x\mathrm{d}t+{\int }_0^T{\int }_{\mathrm{\Omega }}{|▵{\varpi }^{(n)}|}^{q(x)}\mathrm{d}x\mathrm{d}t\hfill \\ \hfill =& h\sum _{k=1}^n{\int }_{\mathrm{\Omega }}|\nabla ▵v_k{|}^{p(x)}\mathrm{d}x+h\sum _{k=1}^n{\int }_{\mathrm{\Omega }}{|▵v_k|}^{q(x)}\mathrm{d}x\le C.\hfill \end{array}$$

☐

Another approximate solution is defined as follows:

$$\begin{array}{cc}\hfill v^{(n)}(x,t)& =\sum _{k=1}^n{\chi }_k(t)[{\vartheta }_k(t)v_k(x)+(1-{\vartheta }_k(t))v_{k-1}(x)],\hfill \end{array}$$

(33)

where

$$\begin{array}{c}\hfill {\vartheta }_k(t)=\left\{\begin{array}{cc}\frac{t}{h}-(k-1),\hfill & \mathrm{if} t\in ((k-1)h,kh],\hfill \\ 0,\hfill & \mathrm{otherwise}.\hfill \end{array}\right.\end{array}$$

We also obtain some uniform estimates for this approximate solution.

Lemma 8.

One has

$$\begin{array}{c}\hfill {\parallel \frac{\partial v^{(n)}}{\partial t}\parallel }_{{\mathcal{L}}^{p^{\prime -}}(0,T;{\mathcal{W}}^{-1,p^{\prime }(x)}(\mathrm{\Omega }))}+\parallel v^{(n)}{\parallel }_{{\mathcal{L}}^{\infty }(0,T;{\mathcal{H}}_0^1(\mathrm{\Omega }))}\le C.\end{array}$$

(34)

Proof. 

By (33), we get

$$\begin{array}{c}\hfill \frac{\partial v^{(n)}}{\partial t}=\frac{1}{h}\sum _{k=1}^n{\chi }_k(v_k-v_{k-1})\end{array}$$

(35)

and then

$$\begin{array}{cc}\hfill {\int }_0^T\parallel \frac{\partial v^{(n)}}{\partial t}{\parallel }_{{\mathcal{W}}^{-1,p^{\prime }(x)}(\mathrm{\Omega })}^{p^{\prime -}}\mathrm{d}t\le & {\int }_0^T|〈\frac{\partial v^{(n)}}{\partial t},\psi 〉{|}^{p^{\prime -}}\mathrm{d}t\hfill \\ \hfill \le & C{\int }_0^T|\sum _{k=1}^n{\chi }_k(t){\int }_{\mathrm{\Omega }}{|\nabla ▵v_k|}^{p(x)-2}\nabla ▵v_k\nabla \psi \mathrm{d}x{|}^{p^{\prime -}}\mathrm{d}t\hfill \\ & +C{\int }_0^T|\sum _{k=1}^n{\chi }_k(t){\int }_{\mathrm{\Omega }}g(x,v_{k-1})\psi \mathrm{d}x{|}^{p^{\prime -}}\mathrm{d}t\hfill \\ & +C{\int }_0^T|\sum _{k=1}^n{\chi }_k(t){\int }_{\mathrm{\Omega }}{|▵v_k|}^{q(x)-2}▵v_k\psi \mathrm{d}x{|}^{p^{\prime -}}\mathrm{d}t\hfill \\ \hfill \le & C\parallel \nabla ▵{\varpi }^{(n)}{\parallel }_{{\mathcal{L}}^{p^{\prime -}}({\mathrm{\Omega }}_T)}^{p^{\prime -}}+C{\parallel g(x,v_{k-1})\parallel }_{{\mathcal{L}}^{p^{\prime -}}({\mathrm{\Omega }}_T)}^{p^{\prime -}}\hfill \\ \hfill \le & C,\hfill \end{array}$$

for any $\psi \in {\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })$ with ${\parallel \psi \parallel }_{{\mathcal{W}}_0^{1,p(x)}(\mathrm{\Omega })}\le 1$. It follows from (22) and (33) that

$$\begin{array}{cc}& \parallel v^{(n)}{\parallel }_{{\mathcal{L}}^m(0,T;{\mathcal{H}}_0^1(\mathrm{\Omega }))}^m\hfill \\ \hfill \le & C^m{\int }_0^T{\left({\int }_{\mathrm{\Omega }}{|\nabla v^{(n)}|}^2\mathrm{d}x\right)}^{\frac{m}{2}}\mathrm{d}t\hfill \\ \hfill =& C^m{\int }_0^T{\left({\int }_{\mathrm{\Omega }}|\sum _{k=1}^n{\chi }_k(t)[{\vartheta }_k(t)\nabla v_k(x)+(1-{\vartheta }_k(t))\nabla v_{k-1}(x)]{|}^2\mathrm{d}x\right)}^{\frac{m}{2}}\mathrm{d}t\hfill \\ \hfill =& C^m\sum _{k=1}^n{\int }_{(k-1)h}^{(kh)}{\left({\int }_{\mathrm{\Omega }}|[{\vartheta }_k(t)\nabla v_k(x)+(1-{\vartheta }_k(t))\nabla v_{k-1}(x)]{|}^2\mathrm{d}x\right)}^{\frac{m}{2}}\mathrm{d}t\hfill \\ \hfill \le & C^m\sum _{k=1}^{n}h{\left({\int }_{\mathrm{\Omega }}(|\nabla v_k{(x)|}^2+|\nabla v_{k-1}(x){|}^2)\mathrm{d}x\right)}^{\frac{m}{2}}\hfill \\ \hfill \le & C^{{}^m+\frac{m}{2}}T\hfill \end{array}$$

with $C>0$ independent of m. Perform the limit $m\to \infty$ to get

$$\begin{array}{c}\hfill {\parallel v^{(n)}\parallel }_{{\mathcal{L}}^{\infty }(0,T;{\mathcal{H}}_0^1(\mathrm{\Omega }))}=\underset{m\to \infty }{lim}\parallel v^{(n)}{\parallel }_{{\mathcal{L}}^m(0,T;{\mathcal{H}}_0^1(\mathrm{\Omega }))}\le C.\end{array}$$

☐

Proof of Theorem 2.

By (29), we can seek a subsequence of ${\varpi }^{(n)}$ (still denoted by itself) and two functions $\nu \in {\mathcal{L}}^{p^{\prime }(x)}({\mathrm{\Omega }}_T),{\nu }^{\prime }\in {\mathcal{L}}^{q^{\prime }(x)}({\mathrm{\Omega }}_T)$ such that

$$\begin{array}{c}{\varpi }^{(n)}\rightharpoonup v\mathrm{weakly}*\mathrm{in}{\mathcal{L}}^{\infty }(0,T;{\mathcal{H}}_0^1(\mathrm{\Omega })),\hfill \\ ▵{\varpi }^{(n)}\rightharpoonup ▵v\mathrm{weakly}\mathrm{in}{\mathcal{L}}^{p^{-}}(0,T;{\mathcal{W}}_0^{1,p}(\mathrm{\Omega })),\hfill \\ |\nabla ▵{\varpi }^{(n)}{|}^{p(x)-2}\nabla ▵{\varpi }^{(n)}\rightharpoonup \nu \mathrm{weakly}\mathrm{in}{\mathcal{L}}^{p^{\prime }(x)}({\mathrm{\Omega }}_T),\hfill \\ |▵{\varpi }^{(n)}{|}^{q(x)-2}▵{\varpi }^{(n)}\rightharpoonup {\nu }^{\prime }\mathrm{weakly}\mathrm{in}{\mathcal{L}}^{q^{\prime }(x)}({\mathrm{\Omega }}_T),\hfill \end{array}$$

as $n\to \infty$.

It is easy to check that there exists a positive integer r such that ${\mathcal{W}}^{-1,p^{\prime }}(\mathrm{\Omega })\hookrightarrow {\mathcal{H}}^{-r}(\mathrm{\Omega })$ and thus the embedding ${\mathcal{H}}_0^1(\mathrm{\Omega })\stackrel{\mathrm{compact}}{\hookrightarrow }{\mathcal{L}}^2(\mathrm{\Omega })\hookrightarrow {\mathcal{H}}^{-r}(\mathrm{\Omega })$, the uniform estimate (34) and the Aubin lemma [17] yield the existence of a subsequence of $v^{(n)}$ and a function ϱ such that, as $n\to \infty$,

$$\begin{array}{c}\frac{\partial v^{(n)}}{\partial t}\rightharpoonup \frac{\partial \varrho }{\partial t}\mathrm{weakly}\mathrm{in}{\mathcal{L}}^{p^{\prime -}}(0,T;{\mathcal{W}}^{-1,p(x)}(\mathrm{\Omega })),\hfill \\ v^{(n)}\rightharpoonup \varrho \mathrm{weakly}*\mathrm{in}{\mathcal{L}}^{\infty }(0,T;{\mathcal{H}}_0^1(\mathrm{\Omega })),\hfill \\ v^{(n)}\to \varrho \mathrm{strongly}\mathrm{in}C([0,T];{\mathcal{L}}^2(\mathrm{\Omega })),\hfill \\ v^{(n)}\to \varrho \mathrm{a}.\mathrm{e}.\mathrm{in}{\mathrm{\Omega }}_T.\hfill \end{array}$$

Moreover, (23) gives, for any $\varphi \in C_0^{\infty }({\mathrm{\Omega }}_T)$,

$$\begin{array}{cc}& \left|{\int }_0^T{\int }_{\mathrm{\Omega }}({\varpi }^{(n)}-v^{(n)})\varphi \mathrm{d}x\mathrm{d}t\right|\hfill \\ \hfill =& \left|{\int }_0^T{\int }_{\mathrm{\Omega }}\sum _{k=1}^n{\chi }_k(t)(1-{\vartheta }_k(t))(v_k-v_{k-1})\varphi \mathrm{d}x\mathrm{d}t\right|\hfill \\ \hfill =& |{\int }_0^{T}h\sum _{k=1}^n{\chi }_k(1-{\vartheta }_k){\int }_{\mathrm{\Omega }}(|\nabla ▵v_k{|}^{p(x)-2}\nabla ▵v_k\nabla \varphi \hfill \\ & +|▵v_k{|}^{q(x)-2}▵v_k\varphi +g(x,v_{k-1})\varphi )\mathrm{d}x\mathrm{d}t|\hfill \\ \hfill \le & h{\int }_0^T(|{\int }_{\mathrm{\Omega }}|\nabla ▵{\varpi }^{(n)}{|}^{p(x)-2}\nabla ▵{\varpi }^{(n)}\nabla \varphi \mathrm{d}x|+|{\int }_{\mathrm{\Omega }}{|▵{\varpi }^{(n)}|}^{q(x)-2}▵{\varpi }^{(n)}\varphi \mathrm{d}x|\hfill \\ & +\left|{\int }_{\mathrm{\Omega }}g(x,{\varpi }^{(n)})\varphi \mathrm{d}x\right|+\left|{\int }_{\mathrm{\Omega }}g(x,v_0)\varphi \mathrm{d}x\right|)\mathrm{d}t\hfill \\ \hfill \le & Ch\to 0,\hfill \end{array}$$

as $n\to \infty$.

By the continuity of g and $\varrho =v$ a.e. in ${\mathrm{\Omega }}_T$, we have $g(x,{\varpi }^{(n)})\to g(x,v)$ a.e. in ${\mathrm{\Omega }}_T$. Furthermore, the estimate (see Lemma 7) $\parallel g(x,{\varpi }^{(n)}){\parallel }_{{\mathcal{L}}^{p^{\prime }(x)}({\mathrm{\Omega }}_T)}<\infty$ gives

$$\begin{array}{c}\hfill g(x,{\varpi }^{(n)})\rightharpoonup g(x,v)\mathrm{weakly}\mathrm{in}{\mathcal{L}}^{p^{\prime }(x)}({\mathrm{\Omega }}_T).\end{array}$$

Applying (23) and (35), we obtain, for any test function ϕ,

$$\begin{array}{cc}& {\int }_0^T〈\frac{\partial v^{(n)}}{\partial t},\varphi 〉\mathrm{d}t\hfill \\ \hfill =& {\int }_0^T{\int }_{\mathrm{\Omega }}{|\nabla ▵{\varpi }^{(n)}|}^{p(x)-2}\nabla ▵{\varpi }^{(n)}\nabla \varphi \mathrm{d}x\mathrm{d}t\hfill \\ & +{\int }_0^T{\int }_{\mathrm{\Omega }}{|▵{\varpi }^{(n)}|}^{q(x)-2}▵{\varpi }^{(n)}\varphi \mathrm{d}x\mathrm{d}t\hfill \\ & +{\int }_0^T{\int }_{\mathrm{\Omega }}g(x,{\varpi }^{(n)})\varphi \mathrm{d}x\mathrm{d}t+h{\int }_{\mathrm{\Omega }}g(x,v_0)\varphi (x,0)\mathrm{d}x\hfill \\ & -{\int }_{(n-1)h}^T{\int }_{\mathrm{\Omega }}g(x,v_n)\varphi \mathrm{d}x\mathrm{d}t.\hfill \end{array}$$

By taking $n\to \infty$, we have

$$\begin{array}{c}\hfill {\int }_0^T〈\frac{\partial v}{\partial t},\varphi 〉\mathrm{d}t={\int }_0^T{\int }_{\mathrm{\Omega }}\nu \nabla \varphi \mathrm{d}x\mathrm{d}t+{\int }_0^T{\int }_{\mathrm{\Omega }}{\nu }^{\prime }\varphi \mathrm{d}x\mathrm{d}t+{\int }_0^T{\int }_{\mathrm{\Omega }}g(x,v)\varphi \mathrm{d}x\mathrm{d}t.\end{array}$$

It remains to prove $\nu ={|\nabla ▵v|}^{p(x)-2}\nabla ▵v$ and ${\nu }^{\prime }={|▵v|}^{q(x)-2}▵v$. By using v as a test function in (23) and integrating by part, we get

$$\begin{array}{cc}& \frac{1}{2}{\int }_{\mathrm{\Omega }}{|\nabla v(x,T)|}^2\mathrm{d}x-\frac{1}{2}{\int }_{\mathrm{\Omega }}{|\nabla v_0|}^2\mathrm{d}x+{\int }_0^T{\int }_{\mathrm{\Omega }}\nu \nabla ▵v\mathrm{d}x\mathrm{d}t+{\int }_0^T{\int }_{\mathrm{\Omega }}{\nu }^{\prime }▵v\mathrm{d}x\mathrm{d}t\hfill \\ \hfill =& -{\int }_0^T{\int }_{\mathrm{\Omega }}g(x,v)▵v\mathrm{d}x\mathrm{d}t.\hfill \end{array}$$

(36)

On the other hand, (23) implies

$$\begin{array}{cc}& \frac{1}{2}{\int }_{\mathrm{\Omega }}|\nabla {\varpi }^{(n)}{(x,T)|}^2\mathrm{d}x-\frac{1}{2}{\int }_{\mathrm{\Omega }}{|\nabla v_0|}^2\mathrm{d}x\hfill \\ & +{\int }_0^T{\int }_{\mathrm{\Omega }}|\nabla ▵{\varpi }^{(n)}{|}^{p(x)}\mathrm{d}x\mathrm{d}t+{\int }_0^T{\int }_{\mathrm{\Omega }}{|▵{\varpi }^{(n)}|}^{q(x)}\mathrm{d}x\mathrm{d}t\hfill \\ \hfill \le & -{\int }_0^T{\int }_{\mathrm{\Omega }}g(x,{\varpi }^{(n)})▵{\varpi }^{(n)}\mathrm{d}x\mathrm{d}t-h{\int }_{\mathrm{\Omega }}g(x,v_0)v_0\mathrm{d}x\hfill \\ & +{\int }_{(n-1)h}^T{\int }_{\mathrm{\Omega }}g(x,v_n)v_n\mathrm{d}x\mathrm{d}t.\hfill \end{array}$$

(37)

For any test functions $\varphi ,{\varphi }_1$ and constant $\epsilon >0$, we have

$$\begin{array}{c}\hfill {\int }_0^T{\int }_{\mathrm{\Omega }}(|{\rho }^{(n)}{|}^{p(x)-2}{\rho }^{(n)}-|{\eta }_{\epsilon }{|}^{p(x)-2}{\eta }_{\epsilon })({\rho }^{(n)}-{\eta }_{\epsilon })\mathrm{d}x\mathrm{d}t\ge 0,\end{array}$$

(38)

and

$$\begin{array}{c}\hfill {\int }_0^T{\int }_{\mathrm{\Omega }}(|{\rho }_1^{(n)}{|}^{q(x)-2}{\rho }_1^{(n)}-|{\eta }_{\epsilon }{|}^{q(x)-2}{\eta }_{1\epsilon })({\rho }_1^{(n)}-{\eta }_{1\epsilon })\mathrm{d}x\mathrm{d}t\ge 0,\end{array}$$

(39)

where ${\rho }^{(n)}=\nabla ▵{\varpi }^{(n)}$, ${\eta }_{\epsilon }=\nabla ▵(v-\epsilon \varphi ))$ and ${\rho }_1^{(n)}=▵{\varpi }^{(n)}$, ${\eta }_{1\epsilon }=▵(v-\epsilon {\varphi }_1)$.

By (37)–(39), we arrive at

$$\begin{array}{c}\frac{1}{2}{\int }_{\mathrm{\Omega }}{|\nabla v(x,T)|}^2\mathrm{d}x-\frac{1}{2}{\int }_{\mathrm{\Omega }}|\nabla v_0{|}^2\mathrm{d}x-{\int }_0^T{\int }_{\mathrm{\Omega }}{|\nabla ▵(v-\epsilon \varphi )|}^{p(x)}\mathrm{d}x\mathrm{d}t\hfill \\ +{\int }_0^T{\int }_{\mathrm{\Omega }}{|\nabla ▵(v-\epsilon \varphi )|}^{p(x)-2}\nabla ▵(v-\epsilon \varphi )\nabla ▵v\mathrm{d}x\mathrm{d}t\hfill \\ +{\int }_0^T{\int }_{\mathrm{\Omega }}\nu \nabla ▵(v-\epsilon \varphi )\mathrm{d}x\mathrm{d}t-{\int }_0^T{\int }_{\mathrm{\Omega }}{|▵(v-\epsilon {\varphi }_1)|}^{q(x)}\mathrm{d}x\mathrm{d}t\hfill \\ +{\int }_0^T{\int }_{\mathrm{\Omega }}{|▵(v-\epsilon {\varphi }_1)|}^{q(x)-2}▵(v-\epsilon {\varphi }_1)▵v\mathrm{d}x\mathrm{d}t\hfill \\ +{\int }_0^T{\int }_{\mathrm{\Omega }}{\nu }^{\prime }▵(v-\epsilon {\varphi }_1)\mathrm{d}x\mathrm{d}t\le -{\int }_0^T{\int }_{\mathrm{\Omega }}g(x,v)▵v\mathrm{d}x\mathrm{d}t.\hfill \end{array}$$

Moreover, (36) gives

$$\begin{array}{c}{\int }_0^T{\int }_{\mathrm{\Omega }}{[|\nabla ▵(v-\epsilon \varphi )|}^{p(x)-2}]\nabla ▵(v-\epsilon \varphi )-\nu ]\nabla ▵\varphi \mathrm{d}x\mathrm{d}t\hfill \\ +{\int }_0^T{\int }_{\mathrm{\Omega }}[|▵(v-\epsilon {\varphi }_1){|}^{q(x)-2}]▵(v-\epsilon \varphi )-{\nu }^{\prime }]▵\varphi \mathrm{d}x\mathrm{d}t\le 0.\hfill \end{array}$$

By letting $\epsilon \to 0$, we obtain

$$\begin{array}{c}{\int }_0^T{\int }_{\mathrm{\Omega }}{\left[\right|\nabla ▵v|}^{p(x)-2}]\nabla ▵v-\nu ]\nabla ▵\varphi \mathrm{d}x\mathrm{d}t+{\int }_0^T{\int }_{\mathrm{\Omega }}{\left[\right|▵v|}^{q(x)-2}▵v-{\nu }^{\prime }]▵{\varphi }_1\mathrm{d}x\mathrm{d}t\hfill \\ \le 0.\hfill \end{array}$$

By letting ${\varphi }_1=0$, we have

$$\begin{array}{c}{\int }_0^T{\int }_{\mathrm{\Omega }}{\left[\right|\nabla ▵v|}^{p(x)-2}]\nabla ▵v-\nu ]\nabla ▵\varphi \mathrm{d}x\mathrm{d}t\le 0.\hfill \end{array}$$

The arbitrariness of ϕ yields ${|\nabla ▵v|}^{p(x)-2}\nabla ▵v=\nu$, a.e. in ${\mathrm{\Omega }}_T$. Similarly, we can obtain ${|▵v|}^{q(x)-2}▵v={\nu }^{\prime }$, a.e. in ${\mathrm{\Omega }}_T$. ☐

Proof of Uniqueness.

Let $v_1(x)$ and $v_2(x)$ be two weak solutions to (1)–(3) and $\varphi =v_1-v_2$. By taking $▵\varphi$ as the test function, we get

$$\begin{array}{c}\frac{1}{2}{\int }_0^t{\int }_{\mathrm{\Omega }}\frac{d}{dt}{|\nabla \varphi |}^2\mathrm{d}x\mathrm{d}t\hfill \\ +{\int }_0^t{\int }_{\mathrm{\Omega }}[|\nabla ▵v_1{|}^{p(x)-2}]\nabla ▵v_1-|\nabla ▵v_2{|}^{p(x)-2}]\nabla ▵v_2]\nabla ▵\varphi \mathrm{d}x\mathrm{d}t\hfill \\ +{\int }_0^t{\int }_{\mathrm{\Omega }}[|▵v_1{|}^{q(x)-2}▵v_1-|▵v_2{|}^{q(x)-2}▵v_2]▵\varphi \mathrm{d}x\mathrm{d}t=\mu {\int }_0^t{\int }_{\mathrm{\Omega }}{|\nabla \varphi |}^2\mathrm{d}x\mathrm{d}t.\hfill \end{array}$$

It implies

$$\begin{array}{c}\hfill {\int }_{\mathrm{\Omega }}{|\nabla \varphi |}^2\mathrm{d}x\le 2\mu {\int }_0^t{\int }_{\mathrm{\Omega }}{|\nabla \varphi |}^2\mathrm{d}x\mathrm{d}t,\end{array}$$

where we have used the fact $(|x{|}^{m-2}x-|y{|}^{m-2}y)(x-y)\ge 0$ for $m>1$ and $x,y\in R$ (or $R^N$). By Gronwall’s inequality, we obtain $v_1(x,t)=v_2(x,t)$ a.e. in ${\mathrm{\Omega }}_T$. ☐