Interpolation inequalities for weak solutions of nonlinear parabolic systems

Abstract

The authors investigate differentiability of the solutions of nonlinear parabolic systems of order 2 m in divergence form of the following type

$$\sum _{\left|\alpha \right|\le m}{\left(-1\right)}^{\left|\alpha \right|}{D}^{\alpha }{a}^{\alpha }\left(X,Du\right)+\frac{\partial u}{\partial t}=0.$$

The achieved results are inspired by the paper of Marino and Maugeri 2008, and the methods there applied.

This note can be viewed as a continuation of the study of regularity properties for solutions of systems started in Ragusa 2002, continued in Ragusa 2003 and Floridia and Ragusa 2012 and also as a generalization of the paper by Capanato and Cannarsa 1981, where regularity properties of the solutions of nonlinear elliptic systems with quadratic growth are reached.

Mathematics Subject Classification (2000)

Primary 35K41, 35K55. Secondary 35B65, 35B45, 35D10

1 Introduction

The study of regularity for solutions of partial differential equations and systems has received considerable attention over the last thirty years. On the other hand, little is known concerning parabolic systems in divergence form of order 2m with quadratic growth and the corresponding analytic properties of solutions. To such classes of systems, our attention is devoted.

This note is a natural continuation of the study, carried out in the last decade and a half, of embedding results of Gagliardo-Nirenberg type from which we deduce local differentiability theorems, making use of interpolation theory in Besov spaces (see e.g. [1–6] and [7]).

In this respect, we mention at first the note [8] where the author proves that, let Ω ⊂ ℝⁿ an open set, 0 < T < ∞ and Q = Ω × (-T, 0), $x^0=\left(x_1^0,x_2^0,\dots ,x_n^0\right)\in \Omega$ ρ > 0 and $B\left(\rho \right)=B\left(x^0,\rho \right)=\left\{x=\left(x_1,x_2,\dots ,x_n\right):|x_i-x_i^0|<\rho ,i=1,2,\dots ,n\right\}$, if

$$u\in L^2(-T,0,H^1(\Omega ,{ℝ}^N))\cap C^{0,\lambda }(Q,{ℝ}^N)),\forall 0<\lambda <1$$

(1.1)

is a solution of a second order nonlinear parabolic system of variational type and under the assumptions that the coefficients a^α(x, Du) have quadratic growth is obtained that

$$u\in L^2\left(-a,0,H^{1+\theta }\left(B\left(\sigma \right),{ℝ}^N\right)\right),$$

for every $a\in \left(0,\frac{T}{2}\right)$, ∀θ ∈ (0, 1) and for each cube B(2σ) ⊂⊂ Ω.

In the same paper, Fattorusso stressed that it is not possible to improve this result in such a way to achieve, for each solution u to the above system, the differentiability

$$u\in L^2\left(-a,0,H^2\left(B\left(\sigma \right),{ℝ}^N\right)\right),$$

(1.2)

if, preliminarily, is not ensured the regularity

$$D_{i}u\in L^4\left(-a,0,L^4\left(B\left(\sigma \right),{ℝ}^N\right)\right),i=1,\dots ,n$$

(1.3)

for every a ∈ (0, T), and for every B(2σ) ⊂⊂ Ω,

The technique used in [8] allows the author to achieve, instead of (1.3), the condition

$$D_{i}u\in L^{2\left(1+\theta \right)}\left(-a,0,L^4\left(B\left(\sigma \right),{ℝ}^N\right)\right),i=1,\dots ,n,$$

for every a ∈ (0, T), ∀B(σ) ⊂⊂ Ω and every $\theta \in \left(\frac{n}{n+4\lambda },1\right)$, which is not enough to ensure that is true (1.2)

In [9], under the same assumptions of the previous result [8], the differentiability result (1.2) is proved, for u satisfying (1.1).

Key of this note is the use of interpolation theorems of Gagliardo-Nirenberg type.

The use of interpolation theory, made in [9] and in [1] with montonicity assumption and quadratic growth, as illustrated in [10], has recently allowed Fattorusso and Marino to obtain differentiability also for weak solutions of nonlinear parabolic systems of second order having nonlinearity 1 < q < 2 (see for details [11]).

Inspired by the note mentioned above by Marino and Maugeri, in the present note, the authors extend differentiability properties to the case of parabolic systems of order 2m. More precisely, let Ω be an open subset of ℝⁿ, n > 2, and 0 < T < ∞, aim of this note is to study, in the cylinder Q = Ω × (-T, 0), the problem of interior local differentiability for solutions

$$u\in L^2\left(-T,0,H^m\left(\Omega ,{ℝ}^N\right)\right)\cap C^{m-1,\lambda }\left(Q,{ℝ}^N\right),0<\lambda <1$$

of the nonlinear parabolic systems of order 2m of variational type

$$\sum _{\left|\alpha \right|\le m}{\left(-1\right)}^{\left|\alpha \right|}{D}^{\alpha }{a}^{\alpha }\left(X,Du\right)+\frac{\partial u}{\partial t}=0.$$

Using the above explained idea is proved the following local differentiability with respect to the spatial derivatives

$$u\in L^2\left(-a,0,H^{m+1}\left(B\left(\sigma \right),{ℝ}^N\right)\right)\cap H^1\left(-a,0,L^2\left(B\left(\sigma \right),{ℝ}^N\right)\right),\forall a\in \left(0,T\right),\forall B\left(\sigma \right)\subset \subset \Omega .$$

Let us also mention the considerable note by [1] where the authors prove that a solution u of nonlinear parabolic systems of order 2 with natural growth and coefficients uniformly monotone in D u belongs to

$$L^2\left(-a,0,H^2\left(B\left(\sigma \right),{ℝ}^N\right)\right)\cap H^1\left(-a,0,L^2\left(B\left(\sigma \right),{ℝ}^N\right)\right)$$

Results similar to those obtained by Marino and Maugeri in [9], with stronger assumptions, are obtained by Naumann in [12] and by Naumann and Wolf in [13]. Let us also bear in mind the study made by Campanato in [14] on parabolic systems in divergence form.

We want to finish this historical overview, concerning interior differentiability of weak solutions, recalling the recent note [4] where similar results are achieved for elliptic systems of order 2m.

2 Useful assumptions and results

Let Ω be an bounded open set in ℝⁿ, n > 2, x = (x₁, x₂, ..., x _n ) denotes a generic point therein, 0 < T < ∞ and Q the cylinder Ω × (-T, 0), let N be a positive integer. In Q, we consider the following parabolic metric

$$d\left(X,Y\right)=max\left\{\left|\right|x-y|{|}_n,|t-\tau {|}^{\frac{1}{2}}\right\},X=\left(x,t\right),Y=\left(y,\tau \right).$$

Let us set k a positive integer greater than 1, (·|·) _k and ||·|| _k respectively the scalar product and the norm in ℝ^k. If there is no ambiguity, we omit the index k.

Let k be a nonnegative integer and λ ∈]0, 1]. We denote by $C^{k,\lambda }\left(\bar{Q},{ℝ}^N\right)$ the subspace of $C^k\left(\bar{Q},{ℝ}^N\right)$ of functions $u:\bar{Q}\to {ℝ}^N$ that satisfy a Hölder condition of exponent λ, together with all their derivatives D^αu, |α| ≤ k. If $u\in C^{k,\lambda }\left(\bar{Q},{ℝ}^N\right)$, then we set

$$\left|\right|u|{|}_{C^{k,\lambda }\left(\bar{Q},{ℝ}^N\right)}=\sum _{\left|\alpha \right|\le k}\underset{\bar{Q}}{sup}|\left|D^{\alpha }u\right|{|}_N+\sum _{\left|\alpha \right|=k}{\left[D^{\alpha }u\right]}_{\lambda ,\bar{Q}}$$

where

$${\left[D^{\alpha }u\right]}_{\lambda ,\bar{Q}}=\underset{\begin{array}{c}\hfill X,Y\in \bar{Q}\hfill \\ \hfill X\ne Y\hfill \end{array}}{sup}\frac{\left|\right|D^{\alpha }u\left(X\right)-D^{\alpha }u\left(Y\right)|{|}_N}{d^{\lambda }\left(X,Y\right)}<+\infty ,\forall \alpha :\left|\alpha \right|=k.$$

The space $C^{k,\lambda }\left(\bar{Q},{ℝ}^N\right)$ is a Banach space, provided with the norm

$$\left|\right|u|{|}_{C^{k,\lambda }\left(\bar{Q},{ℝ}^N\right)}=|\left|u\right|{|}_{C^k\left(\bar{Q},{ℝ}^N\right)}+\sum _{\left|\alpha \right|=k}{\left[D^{\alpha }u\right]}_{\lambda ,\bar{Q}}.$$

Definition 2.1 (see e.g. [15, 16]). Let Ω be an bounded open set in ℝⁿ, let k and j be two positive integers, k ≥ j. If p ∈ [1, +∞ [and $u\in C^{\infty }\left(\bar{\Omega },{ℝ}^N\right)$, so we set

$$|u{|}_{j,p,\Omega }={\left(\underset{\Omega }{\int }\sum _{\left|\alpha \right|=j}\left|\right|D^{\alpha }u|{|}_N^p\mathsf{\text{d}}x\right)}^{\frac{1}{p}},|\left|u\right|{|}_{k,p,\Omega }={\left(\sum _{j=0}^k|u{|}_{j,p,\Omega }^p\right)}^{\frac{1}{p}}$$

(2.1)

and denote respectively by H^{k, p}(Ω, ℝ^N) and $H_0^{k,p}\left(\Omega ,{ℝ}^N\right)$ the spaces obtained as closure of $C^{\infty }\left(\bar{\Omega },{ℝ}^N\right)$ and $C_0^{\infty }\left(\Omega ,{ℝ}^N\right)$ regarding the norm ||u||_{k, p,Ω}. The spaces H^{k, p}(Ω, ℝ^N) and $H_0^{k,p}\left(\Omega ,{ℝ}^N\right)$ are known in literature as Sobolev Spaces.

We remark that H^0,p(Ω, ℝ^N) = L^p(Ω, ℝ^N), 1 ≤ p < +∞.

If p = 2, then we shall simply write H^k(Ω, ℝ^N), $H_0^k\left(\Omega ,{ℝ}^N\right)$, |u|_j,Ω, ||u||_k,Ω.

Let Ω be an bounded open set in ℝⁿ, let us set ϑ ∈ (0, 1), p ∈ [1, +∞ [.

Definition 2.2. We say that a function u defined in Ω having values in ℝ^N belongs to H^{ϑ, p}(Ω, ℝ^N) if u ∈ L^p(Ω, ℝ^N) and is finite

$$|u{|}_{\vartheta ,p,\Omega }^p=\underset{\Omega }{\int }\mathsf{\text{d}}x\underset{\Omega }{\int }\frac{\left|\right|u\left(x\right)-u\left(y\right)|{|}_N^p}{\left|\right|x-y|{|}_n^{n+\vartheta p}}\mathsf{\text{d}}y.$$

Definition 2.3. If k is a nonnegative integer, we mean for H^{k+ϑ, p}(Ω, ℝ^N) the subspace of H^{k, p}(Ω, ℝ^N) of functions u ∈ H^k,p(Ω, ℝ^N) such that

$$D^{\alpha }u\in H^{\vartheta ,p}\left(\Omega ,{ℝ}^N\right),\forall \alpha :\left|\alpha \right|=k.$$

We stress that H^{k+ϑ, p}(Ω, ℝ^N) is a Banach space equipped with the following norm

$$\left|\right|u|{|}_{k+\vartheta ,p,\Omega }={\left(\left|\right|u|{|}_{k,p,\Omega }^p+\sum _{\left|\alpha \right|=k}|D^{\alpha }u{|}_{\vartheta ,p,\Omega }^p\right)}^{\frac{1}{p}}.$$

If p = 2, then we shall simply write H^k+ϑ(Ω, ℝ^N) and ||u||_k+ϑ,Ω.

Let k a positive integer, p ∈ [1, +∞ [, ϑ ∈ (0, 1), in the following, we will consider the spaces

$$\begin{array}{c}{L}^p\left(-T,0,H^{k,p}\left(\Omega ,{ℝ}^N\right)\right)\\ =\left\{u\left(x,t\right)|u\left(\cdot ,t\right)\in H^{k,p}\left(\Omega ,{ℝ}^N\right)\mathsf{\text{ for a}}\mathsf{\text{.e}}\mathsf{\text{. }}t\in \left(-T,0\right)\mathsf{\text{ and }}\underset{-T}{\overset{0}{\int }}{∥u\left(\cdot ,t\right)∥}_{k,p,\Omega }^p\mathsf{\text{d}}t<\infty \right\}\end{array}$$

and

$$\begin{array}{c}{L}^p\left(-T,0,H^{k+\theta ,p}\left(\Omega ,{ℝ}^N\right)\right)\\ =\left\{u\left(x,t\right)|u\left(\cdot ,t\right)\in H^{k+\theta ,p}\left(\Omega ,{ℝ}^N\right)\mathsf{\text{for}}\mathsf{\text{a}}\mathsf{\text{.e}}\mathsf{\text{.}}t\in \left(-T,0\right)\mathsf{\text{ and }}\underset{-T}{\overset{0}{\int }}|\left|u\left(\cdot ,t\right)\right|{|}_{k+\theta ,p,\Omega }^p\mathsf{\text{d}}t<\infty \right\}.\end{array}$$

We say a function u ∈ L²(-T, 0, H^m(Ω, ℝ^N) ∩ C^m-1,λ(Q, ℝ^N), N positive integer and 0 < λ < 1, weak solution in Q to the nonlinear parabolic system of order 2m

$$\sum _{\left|\alpha \right|\le m}{\left(-1\right)}^{\left|\alpha \right|}{D}^{\alpha }{a}^{\alpha }\left(X,Du\right)+\frac{\partial u}{\partial t}=0$$

if

$$\underset{Q}{\int }\left\{\sum _{\left|\alpha \right|\le m}\left(a^{\alpha }\left(X,Du\right)|D^{\alpha }\phi \right)-\left(u|\frac{\partial \phi }{\partial t}\right)\right\}\mathsf{\text{d}}X=0,\forall \phi \in C_0^{\infty }\left(Q,{ℝ}^N\right).$$

(2.2)

Let us now state some properties useful in the sequel.

Let τ ∈]0, 1[, ρ and a two positive numbers and h ∈ ℝ\{0}, where |h| < (1 - τ)ρ.

If u is a function from B(ρ) × (-a, 0) in ℝ^N and X = (x, t) ∈ B(τρ) × (-a, 0), we set

$${\tau }_{i,h}u\left(X\right)=u\left(x+h{e}^i,t\right)-u\left(X\right),i=1,2,\dots ,n,$$

(2.3)

where {eⁱ}_i=1,2,...,nis the canonic basis of ℝⁿ.

Let us now state the following results, proved in [17, 18] and [19], useful to achieve the main result of the note.

Theorem 2.1. If u ∈ L^p(-a, 0, L^p(B(2ρ), ℝ^N)), a, ρ > 0, 1 < p < +∞, N is a positive integer and exists M > 0 such that

$$\underset{-a}{\overset{0}{\int }}\mathsf{\text{d}}t\underset{B\left(\rho \right)}{\int }\left|\right|{\tau }_{i,h}u|{|}^p\mathsf{\text{d}}x\le |h{|}^{p}M,\forall \left|h\right|<\left(1-\tau \right)\rho ,\forall i=1,\dots ,n,$$

then u ∈ L^p(-a, 0, H^1,p(B(ρ), ℝ^N)) and

$$\underset{-a}{\overset{0}{\int }}\mathsf{\text{d}}t\underset{B\left(\rho \right)}{\int }\left|\right|D_{i}u|{|}^p\mathsf{\text{d}}x\le M,\forall i=1,\dots ,n.$$

Theorem 2.2. Let u ∈ H^1,p(B(ρ), ℝ^N) for a, ρ > 0, 1 ≤ p < +∞ and N be a positive integer. Then, for every τ ∈ (0, 1) and every h ∈ ℝ, |h| < (1 - τ)ρ, we have

$$\left|\right|{\tau }_{i,h}u|{|}_{0,p,B\left(\tau \rho \right)}\le |h\left|\right|\left|D_{i}u\right|{|}_{0,p,B\left(\rho \right)},i=1,2,\dots ,n.$$

Theorem 2.3 (see [18, 20]). Let N be a positive integer and Ω a cube of ℝⁿ. If

$$u\in W^{m,r}\left(\Omega ,{ℝ}^N\right)\cap C^{s,\lambda }\left(\Omega ,{ℝ}^N\right),$$

with m ≥ 2, m integer, 1 < r < ∞, s ≥ 0, s integer, 0 < λ < 1, s < m - 1, then, for each integer j with s +λ < j < m, there exists two constants c₁ and c₂ (depending on Ω, m, r, s, λ, j) such that

$$\underset{\left|\alpha \right|=j}{max}|D^{\alpha }u{|}_{0,p,\Omega }\le c_1{\left(\underset{\left|\alpha \right|=m}{max}|D^{\alpha }u{|}_{0,r,\Omega }\right)}^{\delta }\cdot {\left(\underset{\left|\alpha \right|=s}{max}{\left[D^{\alpha }u\right]}_{\lambda ,\Omega }\right)}^{1-\delta }+c_2\underset{\left|\alpha \right|=s}{max}{\left[D^{\alpha }u\right]}_{\lambda ,\Omega }$$

where $\frac{1}{p}=\frac{j}{n}+\delta \left(\frac{1}{r}-\frac{m}{n}\right)-\left(1-\delta \right)\frac{s+\lambda }{n}$, $\forall \delta \in \left[\frac{j-s-\lambda }{m-s-\lambda },1\right[$.

Theorem 2.4 (see [9]). Let N be a positive integer and Ω a cube of ℝⁿ. If

$$u\in W^{m+\theta ,r}\left(\Omega ,{ℝ}^N\right)\cap C^{s,\lambda }\left(\Omega ,{ℝ}^N\right),$$

with m ≥ 1, m integer, 0 < θ < 1, 1 < r < ∞, s ≥ 0, s integer, 0 < λ < 1, s < m, then, for each integer j with $max\left(s+\lambda ,m+\theta -\frac{n}{r}\right)<j<m+\theta$, it results

$$u\in W^{j,p}\left(\Omega ,{ℝ}^N\right)$$

and there exists a constant c (depending on Ω, m, θ, r, s, λ, j, n, δ) such that

$$\left|\right|u|{|}_{j,p,\Omega }\le c|\left|u\right|{|}_{m+\theta ,r,\Omega }^{\delta },\left|\right|u|{|}_{C^{s,\lambda }\left(\Omega ,{ℝ}^N\right)}^{1-\delta },$$

where $\frac{1}{p}=\frac{j}{n}+\delta \left(\frac{1}{r}-\frac{1m+\theta }{n}\right)-\left(1-\delta \right)\frac{s+\lambda }{n}$, $\forall \delta \in \left[\frac{j-s-\lambda }{m+\theta -s-\lambda }\right.$, 1[with (1 - δ)(s + λ) + δ(m + θ) noninteger.

3 Interior differentiability of the solutions

Let us set m, N positive integers, α = (α₁, ..., α _n ) a multi-index and |α| = α₁ + · · · + α _n the order of α. We denote by $ℛ$ the Cartesian product

$$ℛ=\prod _{\left|\alpha \right|\le m}{ℝ}_{\alpha }^N$$

and $p={\left\{p^{\alpha }\right\}}_{\left|\alpha \right|\le m}$, p^α ∈ ℝ^N, the generic point of $ℛ$. If $p\in ℛ$, we set p = (p', p") where $p^{\prime }={\left\{p^{\alpha }\right\}}_{\left|\alpha \right|<m}\in {ℛ}^{\prime }={\prod }_{\left|\alpha \right|<m}{ℝ}_{\alpha }^N$, $p^{″}={\left\{p^{\alpha }\right\}}_{\left|\alpha \right|=m}\in {ℛ}^{″}={\prod }_{\left|\alpha \right|=m}{ℝ}_{\alpha }^N$, and

$$\left|\right|p|{|}^2=\sum _{\left|\alpha \right|\le m}|\left|p^{\alpha }\right|{|}_N^2,\left|\right|p^{\prime }|{|}^2=\sum _{\left|\alpha \right|<m}|\left|p^{\alpha }\right|{|}_N^2,\left|\right|p^{″}|{|}^2=\sum _{\left|\alpha \right|=m}|\left|p^{\alpha }\right|{|}_N^2.$$

We consider, as usual,

$$\begin{array}{lll}\hfill D_i& =\frac{\partial }{\partial x_i},i=1,\dots ,n;D^{\alpha }=D_1^{{\alpha }_1}{D}_2^{{\alpha }_2}\dots D_n^{{\alpha }_n};& \hfill \\ \hfill Du& ={\left\{D^{\alpha }u\right\}}_{\left|\alpha \right|\le m},D^{\prime }u={\left\{D^{\alpha }u\right\}}_{\left|\alpha \right|<m},D^{″}u={\left\{D^{\alpha }u\right\}}_{\left|\alpha \right|=m}.\end{array}$$

Let us consider the following differential nonlinear variational parabolic system of order 2m :

$$\sum _{\left|\alpha \right|\le m}{\left(-1\right)}^{\left|\alpha \right|}{D}^{\alpha }{a}^{\alpha }\left(X,Du\right)+\frac{\partial u}{\partial t}=0$$

(3.1)

where a^α(X, p) = a^α(X, p', p") are functions of $\Lambda =Q\times ℛ$ in ℝ^N, satisfying the following conditions:

(3.2) for every α : |α| < m and every $p\in ℛ$, the function X → a^α(X, p), defined in Q having values in ℝ^N, is measurable in X;

(3.3) for every α : |α| < m and every X ∈ Q, the function p → a^α(X, p), defined in $ℛ$ having values in ℝ^N, is continuous in p;

(3.4) for every α : |α| < m and every (X, p) ∈ Λ, such that ||p'|| ≤ K, we have

$$\left|\right|a^{\alpha }\left(X,p\right)\left|\right|\le M\left(K\right)\left(\left|f^{\alpha }\left(X\right)\right|+\left|\right|p^{″}|{|}^2\right),$$

where f^α ∈ L²(Q);

(3.5) for every α : |α| = m, the function a^α(X, p', p"), defined in $Q\times ℛ$ having values in ℝ^N, are of class C¹ in $Q\times ℛ$ and, for every $\left(X,p^{\prime },p^{″}\right)\in Q\times ℛ$ with ||p'|| ≤ K, we have

$$\begin{array}{lll}\hfill \left|\right|a^{\alpha }\left|\right|+& \sum _{r=1}^n∥\frac{\partial a^{\alpha }}{\partial x_r}∥+\sum _{k=1}^N\sum _{\left|\beta \right|<m}∥\frac{\partial a^{\alpha }}{\partial p_k^{\beta }}∥\le M\left(K\right)\left(1+\left|\right|p^{″}\left|\right|\right),& \hfill \\ \sum _{k=1}^N\sum _{\left|\beta \right|=m}∥\frac{\partial a^{\alpha }}{\partial p_k^{\beta }}∥\le M\left(K\right);& \hfill \end{array}$$

(3.6) ∃ ν = ν(K) > 0 such that:

$$\sum _{h,k=1}^N\sum _{\left|\alpha \right|=\left|\beta \right|=m}\frac{\partial a_h^{\alpha }\left(X,p\right)}{\partial p_k^{\beta }}{\xi }_h^{\alpha }{\xi }_k^{\beta }\ge \nu \left(K\right)\sum _{\left|\beta \right|=m}\left|\right|{\xi }^{\beta }|{|}_N^2=\nu |\left|\xi \right|{|}^2,$$

for every ξ = (ξ^α) ∈ R" and for every $\left(X,p\right)\in Q\times ℛ$, with ||p'|| ≤ K. If the coefficients a^α satisfy condition (3.6) we say that the system (3.1) is strictly elliptic in Ω.

Theorem 3.1. If u ∈ L²(-T, 0, H^m(Ω, ℝ^N)) ∩ C^m-1,λ(Q, ℝ^N), 0 < λ < 1, is a weak solution of the system (3.1) and if the assumptions (3.2)-(3.6) hold, then ∀B(3σ) = B(x⁰, 3σ) ⊂⊂ Ω, ∀a, b ∈ (0, T), a < b, it results

$$u\in L^4\left(-a,0,H^{m,4}\left(B\left(\sigma \right),{ℝ}^N\right)\right)$$

(3.7)

and the following estimate holds

$$\underset{-a}{\overset{0}{\int }}\left|\right|u|{|}_{m,4,B\left(\sigma \right)}^{4}dt\le c\left(\nu ,K,U,\lambda ,\sigma ,a,b,m,n\right)\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^{\frac{1+\vartheta }{2}}+\underset{-b}{\overset{0}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\}$$

(3.8)

where K = sup _Q ||D' u|| and $U=\left|\right|u|{|}_{C^{m-1,\lambda }\left(\bar{Q},{ℝ}^N\right)}$.

Proof Let us observe that, using Theorem 2.III in [21], for every 0 < ϑ < 1 and $b^{*}=\frac{a+b}{2}$, we have

$$u\in L^2\left(-b^{*},0,H^{m+\vartheta }\left(B\left(\frac{5}{2}\sigma \right),{ℝ}^N\right)\right),$$

and

$$\begin{array}{c}\underset{-b^{*}}{\overset{0}{\int }}|D^{″}u{|}_{\vartheta ,B\left(\frac{5}{2}\sigma \right)}^2\mathsf{\text{d}}t\\ \le c\left(\nu ,K,U,\vartheta ,\lambda ,\sigma ,a,b,m,n\right)\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^{\frac{1+\vartheta }{2}}+\underset{-b^{*}}{\overset{0}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\}.\end{array}$$

(3.9)

Hence, we remark that $u\in C^{m-1,\lambda }\left(\bar{\Omega },{ℝ}^N\right)$, then, it results, for a. e. t ∈ (-b*, 0),

$$u\left(x,t\right)\in H^{m+\vartheta }\left(B\left(\frac{5}{2}\sigma \right),{ℝ}^N\right)\cap C^{m-1,\lambda }\left(\overline{B\left(\frac{5}{2}\sigma \right)},{ℝ}^N\right),\forall 0<\vartheta <1,\forall B\left(3\sigma \right)\subset \subset \Omega .$$

Then, from Theorem 2.4 with $\Omega =B\left(\frac{5}{2}\sigma \right)$, 1 - λ < θ < 1, for $\delta =\frac{1}{2}$, and for a.e. t ∈ (-b*, 0):

$$u\left(x,t\right)\in H^{m,p}\left(B\left(\frac{5}{2}\sigma \right),{ℝ}^N\right),$$

and there exists a constant c = c(θ, λ, σ, m, n) such that

$$\left|\right|u|{|}_{m,p,B\left(\frac{5}{2}\sigma \right)}\le c|\left|u\right|{|}_{m+\theta ,B\left(\frac{5}{2}\sigma \right)}^{\frac{1}{2}}\left|\right|u|{|}_{C^{m-1,\lambda }\left(B\left(\frac{5}{2}\sigma \right),{ℝ}^N\right)}^{\frac{1}{2}},$$

where $p=4+\frac{8\left(\theta +\lambda -1\right)}{n-2\left(\theta +\lambda -1\right)}>4$.

The choice $\theta =1-\frac{\lambda }{2}\left(>1-\lambda \right)$ ensures that for a. e. t ∈ (-b*, 0) we have

$$u\left(x,t\right)\in H^{m,p}\left(B\left(\frac{5}{2}\sigma \right),{ℝ}^N\right),\mathsf{\text{with}}p=4+\frac{4\lambda }{n-\lambda },\forall B\left(3\sigma \right)\subset \subset \Omega .$$

(3.10)

and

$$\left|\right|u|{|}_{m,p,B\left(\frac{5}{2}\sigma \right)}\le c\left(\theta ,\lambda ,\sigma ,m,n\right)|\left|u\right|{|}_{m+1-\frac{\lambda }{2},B\left(\frac{5}{2}\sigma \right)}^{\frac{1}{2}}\left|\right|u|{|}_{C^{m-1,\lambda }\left(B\left(\frac{5}{2}\sigma \right),{ℝ}^N\right)}^{\frac{1}{2}},$$

(3.11)

where $p=4+\frac{4\lambda }{n-\lambda }>4$.

Then we have, for a. e. t ∈ (-b*, 0), the following inclusion between Sobolev spaces

$$u\left(x,t\right)\in H^{m,p}\left(B\left(\frac{5}{2}\sigma \right),{ℝ}^N\right)\subset \subset H^{m,4}\left(B\left(\frac{5}{2}\sigma \right),{ℝ}^N\right)$$

(3.12)

then, using (3.9), written with $\theta =1-\frac{\lambda }{2}$, and (3.10)-(3.12), we have

$$\begin{array}{c}\underset{-b^{*}}{\overset{0}{\int }}\left|\right|u|{|}_{m,4,B\left(\frac{5}{2}\sigma \right)}^4\mathsf{\text{d}}t\le c\left(\sigma \right)\underset{-b^{*}}{\overset{0}{\int }}|\left|u\right|{|}_{m,4+\frac{4\lambda }{n-\lambda },B\left(\frac{5}{2}\sigma \right)}^4\mathsf{\text{d}}t\\ \le c\left(\theta ,\lambda ,\sigma ,m,n\right)\underset{-b^{*}}{\overset{0}{\int }}\left|\right|D^{″}u\left|{|}_{1-\frac{\lambda }{2},B\left(\frac{5}{2}\sigma \right)}^2\right|\left|u\right|{|}_{C^{m-1,\lambda }\left(B\left(\frac{5}{2}\sigma \right),{ℝ}^N\right)}^2\mathsf{\text{d}}t\\ \le c\left(\nu ,K,U,\lambda ,\sigma ,m,n\right)\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^{\frac{1+\vartheta }{2}}+\underset{-b}{\overset{0}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\},\end{array}$$

(3.13)

then it follows the requested inequality (3.8).   ■

Theorem 3.2 (main result). If u ∈ L²(-T, 0, H^m(Ω, ℝ^N)) ∩ C^m-1,λ(Q, ℝ^N), 0< λ < 1, is a weak solution of the system (3.1) and if the assumptions (3.2)-(3.6) hold, then ∀B(3σ) = B(x⁰, 3σ) ⊂⊂ Ω, ∀a, b ∈ (0, T), a < b it results

$$u\in L^2\left(-a,0,H^{m+1}\left(B\left(\sigma \right),{ℝ}^N\right)\right)\cap H^1\left(-a,0,L^2\left(B\left(\sigma \right),{ℝ}^N\right)\right)$$

(3.14)

and the following estimate holds

$$\begin{array}{c}\underset{-a}{\overset{0}{\int }}\left\{|u{|}_{m+1,B\left(\sigma \right)}^2+{\left|\frac{\partial u}{\partial t}\right|}_{0,B\left(\sigma \right)}^2\right\}\mathsf{\text{d}}t\\ \le c\left(\nu ,K,U,\lambda ,\sigma ,a,b,m,n\right)\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^2+\underset{-b}{\overset{0}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\}\end{array}$$

(3.15)

where K = sup _Q ||D' u|| and $U=\left|\right|u|{|}_{C^{m-1,\lambda }\left(\bar{Q},{ℝ}^N\right)}$.

Proof Let us fix B(3σ) = B(x⁰, 3σ) ⊂⊂ Ω, a, b ∈ (0, T) with a < b and h ∈ ℝ such that $\left|h\right|<\frac{\sigma }{2}$, set $b^{*}=\frac{a+b}{2}$, and let $\psi \left(x\right)\in C_0^{\infty }\left({ℝ}^n\right)$ a real function satisfying the following properties 0 ≤ ψ ≤ 1 in ℝⁿ, ψ = 1 in B(σ), ψ = 0 in ℝⁿ\B(2σ), $\left|\right|D\psi \left|\right|\le \frac{c}{\sigma }$ in ℝⁿ.

Let us also define the function ρ _μ (t), for $\mu >\frac{2}{a}$, μ integer, the following real function

$${\rho }_{\mu }\left(t\right)=\left\{\begin{array}{cc}\hfill 1\hfill & \hfill \mathsf{\text{if }}-a\le t\le -\frac{2}{\mu }\hfill \\ \hfill 0\hfill & \hfill \mathsf{\text{if }}t\le -b\mathsf{\text{and }}t\ge -\frac{1}{\mu }\hfill \\ \hfill \frac{t+b}{b-a}\hfill & \hfill \mathsf{\text{if }}-b<t<-a\hfill \\ \hfill -\left(\mu t+1\right)\hfill & \hfill \mathsf{\text{if }}-\frac{2}{\mu }<t<-\frac{1}{\mu }.\hfill \end{array}\right.$$

(3.16)

Moreover set {g _s (t)} the sequence of symmetric regularizing functions such that

$$\begin{array}{c}{g}_s\left(t\right)\in C_0^{\infty }\left(ℝ\right),g_s\left(t\right)\ge 0,g_s\left(t\right)=g_s\left(-t\right),\\ \mathsf{\text{supp }}{g}_s\subset \left[-\frac{1}{s},\frac{1}{s}\right],{\int }_{ℝ}{g}_s\left(t\right)\mathsf{\text{d}}t=1.\end{array}$$

Let i be a positive integer, i ≤ n, and h a real number such that $\left|h\right|<\frac{\sigma }{2}$. For every $\mu >\frac{2}{a}$ and for every $s>max\left\{\mu ,\frac{1}{T-b}\right\}$, let us define the following "test function"

$$\phi \left(X\right)={\tau }_{i,-h}\left\{{\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]\right\},\forall X=\left(x,t\right)\in Q.$$

(3.17)

Substituting in (2.2) the above defined function φ, we have

$$\begin{array}{c}\underset{Q}{\int }\sum _{\left|\alpha \right|=m}\left({\tau }_{i,h}{a}^{\alpha }\left(X,Du\right)|D^{\alpha }\left({\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]\right)\right)\mathsf{\text{d}}X\\ =\underset{Q}{\int }\left({\tau }_{i,h}u|{\psi }^{2m}{\left\{{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]\right\}}^{\prime }\right)\mathsf{\text{d}}X-\\ -\sum _{\left|\alpha \right|<m}\underset{Q}{\int }\left(a^{\alpha }\left(X,Du\right)|{\tau }_{i,-h}{D}^{\alpha }\left\{{\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]\right\}\right)\mathsf{\text{d}}X.\end{array}$$

(3.18)

For every α : |α| = m and a. e. X = (x, t) ∈ Q, we have

$$\begin{array}{c}{\tau }_{i,h}{a}^{\alpha }\left(X,Du\left(X\right)\right)=a^{\alpha }\left(x+h{e}^i,t,Du\left(x+h{e}^i,t\right)\right)-a^{\alpha }\left(X,Du\left(X\right)\right)\\ =\underset{0}{\overset{1}{\int }}\frac{d}{d\eta }{a}^{\alpha }\left(x+\eta h{e}^i,t,Du\left(X\right)+\eta {\tau }_{i,h}Du\left(X\right)\right)d\eta \\ =h\underset{0}{\overset{1}{\int }}\frac{\partial }{\partial x_i}{a}^{\alpha }\left(x+\eta h{e}^i,t,Du\left(X\right)+\eta {\tau }_{i,h}Du\left(X\right)\right)d\eta \\ +\sum _{\left|\beta \right|\le m}\sum _{k=1}^N\left({\tau }_{i,h}{D}^{\beta }{u}_k\left(X\right)\right)\underset{0}{\overset{1}{\int }}\frac{\partial }{\partial p_k^{\beta }}{a}^{\alpha }\left(x+\eta h{e}^i,t,Du\left(X\right)+\eta {\tau }_{i,h}Du\left(X\right)\right)d\eta \\ =h\frac{\widetilde{\partial a^{\alpha }}}{\partial x_i}+\sum _{\left|\beta \right|\le m}\sum _{k=1}^N\left({\tau }_{i,h}{D}^{\beta }{u}_k\left(X\right)\right)\frac{\widetilde{\partial a^{\alpha }}}{\partial p_k^{\beta }},\end{array}$$

where, if b = b(X, p), for simplicity of notation, we set

$$\tilde{b}\left(X\right)=\underset{0}{\overset{1}{\int }}b\left(x+\eta h{e}^i,t,Du\left(X\right)+\eta {\tau }_{i,h}Du\left(X\right)\right)d\eta .$$

Then, equality (3.18) becomes

$$\begin{array}{c}\underset{Q}{\int }\sum _{\left|\alpha \right|=m}\left(h\frac{\widetilde{\partial a^{\alpha }}}{\partial x_i}+\sum _{\left|\beta \right|\le m}\sum _{k=1}^N\left({\tau }_{i,h}{D}^{\beta }{u}_k\left(X\right)\right)\left.\frac{\widetilde{\partial a^{\alpha }}}{\partial p_k^{\beta }}\right|D^{\alpha }\left({\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]\right)\right)\mathsf{\text{d}}X\\ =\underset{Q}{\int }{\psi }^{2m}{\rho }^{\prime }\mu \left({\tau }_{i,h}u|\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right)\mathsf{\text{d}}X+\underset{Q}{\int }\left({\tau }_{i,h}u|{\psi }^{2m}{\rho }_{\mu }{\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]}^{\prime }\right)\mathsf{\text{d}}X-\\ \sum _{\left|\alpha \right|<m}\underset{Q}{\int }\left(a^{\alpha }\left(X,Du\right)|{\tau }_{i,-h}{D}^{\alpha }\left\{{\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]\right\}\right)\mathsf{\text{d}}X.\end{array}$$

Taking into account, for α : |α| = m, that

$$D^{\alpha }\left({\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]\right)={\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}{D}^{\alpha }u\right)*g_s\right]+{\psi }^m{\rho }_{\mu }\sum _{\gamma <\alpha }{c}_{\alpha \gamma }\left(\psi \right)\left[\left({\rho }_{\mu }{\tau }_{i,h}{D}^{\gamma }u\right)*g_s\right]$$

where

$$\left|c_{\alpha \gamma }\left(\psi \right)\right|\le \frac{c\left(m,n\right)}{{\sigma }^{m-\left|\gamma \right|}},$$

we obtain

$$\begin{array}{c}\underset{Q}{\int }{\psi }^{2m}{\rho }_{\mu }\sum _{\left|\alpha \right|=\left|\beta \right|=m}\sum _{k=1}^N\left(\left({\tau }_{i,h}{D}^{\beta }{u}_k\left(X\right)\right)\left.\frac{\widetilde{\partial a^{\alpha }}}{\partial p_k^{\beta }}\right|\left({\rho }_{\mu }{\tau }_{i,h}{D}^{\alpha }u\right)*g_s\right)\mathsf{\text{d}}X\\ =-\sum _{\left|\alpha \right|=\left|\beta \right|=m}\sum _{\gamma <\alpha }\sum _{k=1}^N\underset{Q}{\int }\left(\left({\tau }_{i,h}{D}^{\beta }{u}_k\left(X\right)\right)\left.\frac{\widetilde{\partial a^{\alpha }}}{\partial p_k^{\beta }}\right|{\psi }^m{\rho }_{\mu }{c}_{\alpha \gamma }\left(\psi \right)\left[\left({\rho }_{\mu }{\tau }_{i,h}{D}^{\gamma }u\right)*g_s\right]\right)\mathsf{\text{d}}X\\ -\sum _{\left|\alpha \right|=m}\sum _{\left|\beta \right|<m}\sum _{k=1}^N\underset{Q}{\int }\left(\left({\tau }_{i,h}{D}^{\beta }{u}_k\left(X\right)\right)\left.\frac{\widetilde{\partial a^{\alpha }}}{\partial p_k^{\beta }}\right|D^{\alpha }\left({\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]\right)\right)\mathsf{\text{d}}X\\ -h\sum _{\left|\alpha \right|=m}\underset{Q}{\int }\left(\left.\frac{\widetilde{\partial a^{\alpha }}}{\partial x_i}\right|D^{\alpha }\left({\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]\right)\right)\mathsf{\text{d}}X+\underset{Q}{\int }{\psi }^{2m}{\rho }^{\prime }\mu \left({\tau }_{i,h}u|\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right)\mathsf{\text{d}}X\\ +\underset{Q}{\int }\left({\tau }_{i,h}u|{\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g^{\prime }s\right]\right)\mathsf{\text{d}}X-\sum _{\left|\alpha \right|<m}\underset{Q}{\int }\left(a^{\alpha }\left(X,Du\right)|{\tau }_{i,-h}{D}^{\alpha }\left\{{\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g_s\right]\right\}\right)\mathsf{\text{d}}X.\end{array}$$

For s → +∞, using ellipticity condition (3.6), symmetry hypothesis, convolution property of g _s and that

$$\underset{s\to +\infty }{lim}\underset{Q}{\int }\left({\tau }_{i,h}u|{\psi }^{2m}{\rho }_{\mu }\left[\left({\rho }_{\mu }{\tau }_{i,h}u\right)*g^{\prime }s\right]\right)\mathsf{\text{d}}X=0,$$

we have

$$\begin{array}{c}\nu \underset{-b}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x=\nu \underset{-b}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\sum _{\left|\alpha \right|=m}|\left|{\tau }_{i,h}{D}^{\alpha }u\right|{|}^2\mathsf{\text{d}}x\\ \le \underset{Q}{\int }{\psi }^{2m}{\rho }_{\mu }\sum _{\left|\alpha \right|=\left|\beta \right|=m}\sum _{k=1}^N\left(\left({\tau }_{i,h}{D}^{\beta }{u}_k\left(X\right)\right)\left.\frac{\widetilde{\partial a^{\alpha }}}{\partial p_k^{\beta }}\right|{\rho }_{\mu }{\tau }_{i,h}{D}^{\alpha }u\right)\mathsf{\text{d}}X\\ \le A+B+C+D+E,\end{array}$$

(3.19)

where

$$A=-\sum _{\left|\alpha \right|=\left|\beta \right|=m}\sum _{\gamma <\alpha }\sum _{k=1}^N\underset{Q}{\int }{c}_{\alpha \gamma }\left(\psi \right){\psi }^m{\rho }_{\mu }^2\left(\left({\tau }_{i,h}{D}^{\beta }{u}_k\left(X\right)\right)\left.\frac{\widetilde{\partial a^{\alpha }}}{\partial p_k^{\beta }}\right|\left({\tau }_{i,h}{D}^{\gamma }u\right)\left(X\right)\right)\mathsf{\text{d}}X,$$

(3.20)

$$B=-\sum _{\left|\alpha \right|=m}\sum _{\left|\beta \right|<m}\sum _{k=1}^N\underset{Q}{\int }\left(\left({\tau }_{i,h}{D}^{\beta }{u}_k\left(X\right)\right)\left.\frac{\widetilde{\partial a^{\alpha }}}{\partial p_k^{\beta }}\right|D^{\alpha }\left({\psi }^{2m}{\rho }_{\mu }^2{\tau }_{i,h}u\left(X\right)\right)\right)\mathsf{\text{d}}X,$$

(3.21)

$$C=-h{\displaystyle \sum _{\left|\alpha \right|=m}{\displaystyle \underset{Q}{\int }\left(\frac{\widetilde{\partial a^{\alpha }}}{\partial x_i}|D^{\alpha }\left({\psi }^{2m}{\rho }_{\mu }^2{\tau }_{i,h}u)(X)\right)\right)}}\text{d}X,$$

(3.22)

$$D=\underset{Q}{\int }{\psi }^{2m}{\rho }^{\prime }\mu {\rho }_{\mu }\left|\right|{\tau }_{i,h}u|{|}^2\mathsf{\text{d}}X,$$

(3.23)

$$E=-\sum _{\left|\alpha \right|<m}\underset{Q}{\int }\left(a^{\alpha }\left(X,Du\right)|{\tau }_{i,-h}{D}^{\alpha }\left({\psi }^{2m}{\rho }_{\mu }^2{\tau }_{i,h}u\right)\right)\mathsf{\text{d}}X.$$

(3.24)

We observe that, for every ε > 0, we have

$$\left|A\right|\le \epsilon {\displaystyle \underset{-b}{\overset{-\frac{1}{\mu }}{\int }}\text{d}t}{\displaystyle \underset{B(2\sigma )}{\int }{\psi }^{2m}}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\text{d}x+c(K,\sigma ,m,n,\epsilon )h^2{\displaystyle \underset{-b}{\overset{-\frac{1}{\mu }}{\int }}\text{d}t}{\displaystyle \underset{B(3\sigma )}{\int }(1+|\left|D^{″}u\right|{|}^2\text{)}}\text{d}x.$$

(3.25)

The term B can be estimated, for every ε > 0, as follows

$$\begin{array}{c}\left|B\right|\le \left\{\epsilon +c\left(K,U,m,n\right)\left(|h{|}^{\lambda }+|h{|}^{2\lambda }\right)\right\}\underset{-b}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x\\ +c\left(K,\sigma ,m,n,\epsilon \right)h^2\underset{-b}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(3\sigma \right)}{\int }\left|\right|D^{″}u|{|}^2\mathsf{\text{d}}x\\ +c\left(K,m,n,\epsilon \right)\underset{-b}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{\prime }u\left|{|}^2\right|\left|D^{″}u\right|{|}^2\mathsf{\text{d}}x.\end{array}$$

(3.26)

Let us consider the term C, for every ε > 0, we have

$$\begin{array}{c}\left|C\right|\le \left\{\epsilon +c\left(K,m,n\right)\left(h^2+\left|h\right|\right)\right\}\underset{-b}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x\\ +c\left(K,\sigma ,m,n,\epsilon \right)h^2\underset{-b}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(3\sigma \right)}{\int }\left(1+\left|\right|D^{″}u|{|}^2\right)\mathsf{\text{d}}x.\end{array}$$

To estimate the term D, we firstly observe that

$$\left({\rho }^{\prime }\mu {\rho }_{\mu }\right)\left(t\right)\left\{\begin{array}{cc}\hfill =0\hfill & \hfill \mathsf{\text{if }}t\le -b,-a\le t\le -\frac{2}{\mu },t\ge -\frac{1}{\mu }\hfill \\ \hfill \le \frac{1}{b-a}\hfill & \hfill \mathsf{\text{if }}-b\le t\le -a\hfill \\ \hfill \le 0\hfill & \hfill \mathsf{\text{if }}-\frac{2}{\mu }\le t\le -\frac{1}{\mu }\hfill \end{array}\right.$$

(3.27)

then, using Theorem 2.2, we obtain

$$\begin{array}{c}D=\underset{Q}{\int }{\psi }^{2m}{\rho }^{\prime }\mu {\rho }_{\mu }\left|\right|{\tau }_{i,h}u|{|}^2\mathsf{\text{d}}X\\ =\underset{-b}{\overset{-a}{\int }}\mathsf{\text{d}}t\underset{B_{\left(2\sigma \right)}}{\int }{\psi }^{2m}{\rho }^{\prime }\mu {\rho }_{\mu }\left|\right|{\tau }_{i,h}u|{|}^2\mathsf{\text{d}}x+\underset{-\frac{2}{\mu }}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B_{\left(2\sigma \right)}}{\int }{\psi }^{2m}{\rho }^{\prime }\mu {\rho }_{\mu }|\left|{\tau }_{i,h}u\right|{|}^2\mathsf{\text{d}}x\\ \le \frac{1}{b-,a}\underset{-b}{\overset{-a}{\int }}\mathsf{\text{d}}t\underset{B_{\left(2\sigma \right)}}{\int }\left|\right|{\tau }_{i,h}u|{|}^2\mathsf{\text{d}}x\le \frac{h^2}{b-a}\underset{-b}{\overset{-a}{\int }}\mathsf{\text{d}}t\underset{B_{\left(3\sigma \right)}}{\int }|\left|D_{i}u\right|{|}^2\mathsf{\text{d}}x.\end{array}$$

(3.28)

Finally, using (3.4) condition, the term E can be expressed as follows

$$\left|E\right|\le c\left(K,m,n\right)\sum _{\left|\alpha \right|<m}\underset{-b}{\overset{-\frac{1}{\mu }}{\int }}{\rho }_{\mu }^2\mathsf{\text{d}}t\underset{B\left(\frac{5}{2}\sigma \right)}{\int }\left(\left|f^{\alpha }\right|+\left|\right|D^{″}u|{|}^2\right)\left|\right|{\tau }_{i,-h}{D}^{\alpha }\left({\psi }^{2m}{\tau }_{i,h}u\right)\left|\right|\mathsf{\text{d}}x.$$

(3.29)

Then, from (3.19) estimating the terms A, B, C, D, and E, for every ε > 0, we have

$$\begin{array}{c}\nu \underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x\\ \le \left\{3\epsilon +c\left(K,U,m,n\right)\left(\left|h\right|+h^2+|h{|}^{\lambda }+|h{|}^{2\lambda }\right)\right\}\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x\\ +c\left(K,\sigma ,a,b,m,n,\epsilon \right)h^2\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(3\sigma \right)}{\int }\left(1+\left|\right|D^{″}u|{|}^2\right)\mathsf{\text{d}}x+c\left(\sigma ,a,b,n,\right)K{h}^2\\ +c\left(K,\sigma ,m,n,\epsilon \right)\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{\prime }u\left|{|}^2\right|\left|D^{″}u\right|{|}^2\mathsf{\text{d}}x\\ +c\left(K,m,n\right)\sum _{\left|\alpha \right|<m}\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}{\rho }_{\mu }^2\mathsf{\text{d}}t\underset{B\left(\frac{5}{2}\sigma \right)}{\int }\left(\left|f^{\alpha }\right|+\left|\right|D^{″}u|{|}^2\right)\left|\right|{\tau }_{i,-h}{D}^{\alpha }\left({\psi }^{2m}{\tau }_{i,h}u\right)\left|\right|\mathsf{\text{d}}x.\end{array}$$

(3.30)

We observe that the function

$$h\to c\left(K,U,\sigma ,m,n\right)\left(\left|h\right|+h^2+|h{|}^{\lambda }+|h{|}^{2\lambda }\right)$$

is continuous in the origin, then ∃h₀(ν, K, U, λ, σ, m, n), $0<h_0<min\left\{1,\frac{\sigma }{2}\right\}$, such that for every |h| < h₀, we have

$$c\left(K,U,\sigma ,m,n\right)\left(\left|h\right|+h^2+|h{|}^{\lambda }+|h{|}^{2\lambda }\right)<\frac{\nu }{4}.$$

For each integer i = 1, ..., n, for $\epsilon =\frac{\nu }{12}$ and every h such that |h| < h₀(< 1), it follows

$$\begin{array}{c}\frac{\nu }{2}\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x\\ \le c\left(\nu ,K,\sigma ,a,b,m,n\right)|h{|}^2\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(3\sigma \right)}{\int }\left(1+\left|\right|D^{″}u|{|}^2\right)\mathsf{\text{d}}x\\ +c\left(\nu ,K,\sigma ,m,n\right)\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{\prime }u\left|{|}^2\right|\left|D^{″}u\right|{|}^2\mathsf{\text{d}}x\\ +c\left(K,m,n\right)\sum _{\left|\alpha \right|<m}\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}{\rho }_{\mu }^2\mathsf{\text{d}}t\underset{B\left(\frac{5}{2}\sigma \right)}{\int }\left(\left|f^{\alpha }\right|+\left|\right|D^{″}u|{|}^2\right)\left|\right|{\tau }_{i,-h}{D}^{\alpha }\left({\psi }^{2m}{\tau }_{i,h}u\right)\left|\right|\mathsf{\text{d}}x.\end{array}$$

(3.31)

Let us focus our attention on the last term, taking into account that from (3.12), for a. e. t ∈ (-b*, 0), we have

$$u\left(\cdot ,t\right)\in H^{m,4}\left(B\left(\frac{5}{2}\sigma \right),{ℝ}^N\right)$$

then using Hölder and Young inequalities, for every α such that |α| < m, for every ε > 0, it follows

$$\begin{array}{c}\underset{B\left(\frac{5}{2}\sigma \right)}{\int }\left(\left|f^{\alpha }\right|+\left|\right|D^{″}u|{|}^2\right)\left|\right|{\tau }_{i,-h}{D}^{\alpha }\left({\psi }^{2m}{\tau }_{i,h}u\right)\left|\right|\mathsf{\text{d}}x\\ \le {\left(\underset{B\left(3\sigma \right)}{\int }\left|h{|}^{-2}\right|\left|{\tau }_{i,-h}{D}^{\alpha }\left({\psi }^{2m}{\tau }_{i,h}u\right)\right|{|}^2\mathsf{\text{d}}x\right)}^{\frac{1}{2}}{\left(\underset{B\left(\frac{5}{2}\sigma \right)}{\int }{h}^2{\left(\left|f^{\alpha }\right|+\left|\right|D^{″}u|{|}^2\right)}^2\mathsf{\text{d}}x\right)}^{\frac{1}{2}}\\ \le \frac{\epsilon }{2}\left|h{|}^{-2}\underset{B\left(3\sigma \right)}{\int }\right|\left|{\tau }_{i,-h}{D}^{\alpha }\left({\psi }^{2m}{\tau }_{i,h}u\right)\right|{|}^2\mathsf{\text{d}}x+c\left(\epsilon \right)h^2\underset{B\left(\frac{5}{2}\sigma \right)}{\int }\left(|f^{\alpha }{|}^2+|\left|D^{″}u\right|{|}^4\right)\mathsf{\text{d}}x.\end{array}$$

Furthermore, for every α such that |α| < m, from Theorem 2.2 for every h ∈ ℝ with |h| < h₀ and for every ε > 0, we have

$$\begin{array}{c}\frac{\epsilon }{2}{\left|h\right|}^{-2}\underset{B\left(3\sigma \right)}{\int }\left|\right|{\tau }_{i,-h}{D}^{\alpha }\left({\psi }^{2m}{\tau }_{i,h}u\right)|{|}^2\mathsf{\text{d}}x\le \frac{\epsilon }{2}\underset{B\left(\frac{7}{2}\sigma \right)}{\int }|\left|D^{″}\left({\psi }^{2m}{\tau }_{i,h}u\right)\right|{|}^2\mathsf{\text{d}}x\\ \le \epsilon \underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x+c\left(\sigma ,\epsilon \right)\underset{B\left(2\sigma \right)}{\int }|\left|{\tau }_{i,h}{D}^{\prime }u\right|{|}^2\mathsf{\text{d}}x\\ \le \epsilon \underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x+c\left(\sigma ,\epsilon \right)h^2\underset{B\left(3\sigma \right)}{\int }|\left|D^{″}u\right|{|}^2\mathsf{\text{d}}x\end{array}$$

the last inequality follows, as before, applying Theorem 2.2 for p = 2.

Let us now choose $\epsilon =\frac{\nu }{4c\left(K,m,n\right)}$, it ensures

$$\begin{array}{c}\underset{B\left(\frac{5}{2}\sigma \right)}{\int }\left(\left|f^{\alpha }\right|+\left|\right|D^{″}u|{|}^2\right)\left|\right|{\tau }_{i,-h}{D}^{\alpha }\left({\psi }^{2m}{\tau }_{i,h}u\right)\left|\right|\mathsf{\text{d}}x\\ \le \frac{\nu }{4c\left(K,m,n\right)}\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x+c\left(\nu ,K,\sigma ,m,n\right)h^2\left\{\underset{B\left(3\sigma \right)}{\int }|f^{\alpha }{|}^2\mathsf{\text{d}}x+|u{|}_{m,B\left(3\sigma \right)}^2+|u{|}_{m,4,B\left(\frac{5}{2}\sigma \right)}^4\right\}.\end{array}$$

Multiplying each term for ${\rho }_{\mu }^2$ and integrating respect to $\left(-b^{*},-\frac{1}{\mu }\right)$ and applying (3.13), we achieve

$$\begin{array}{c}\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}{\rho }_{\mu }^2\mathsf{\text{d}}t\underset{B\left(\frac{5}{2}\sigma \right)}{\int }\left(\left|f^{\alpha }\right|+\left|\right|D^{″}u|{|}^2\right)\left|\right|{\tau }_{i,-h}{D}^{\alpha }\left({\psi }^{2m}{\tau }_{i,h}u\right)\left|\right|\mathsf{\text{d}}x\\ \le \frac{\nu }{4c\left(K,m,n\right)}\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}{\rho }_{\mu }^2\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x\\ +c\left(\nu ,K,U,\lambda ,\sigma ,a,b,m,n\right)h^2\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^{\frac{1+\vartheta }{2}}+\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\}.\end{array}$$

Taking into consideration the last inequality and the properties of the function ψ, from (3.31) we deduce

$$\begin{array}{c}\underset{-a}{\overset{-\frac{2}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(\sigma \right)}{\int }\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x\\ \le c\left(\nu ,K,U,\lambda ,\sigma ,a,b,m,n\right)h^2\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^{\frac{1+\vartheta }{2}}+\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\}\\ +c\left(\nu ,K,\sigma ,m,n\right)\underset{-b^{*}}{\overset{-\frac{1}{\mu }}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{\prime }u\left|{|}^2\right|\left|D^{″}u\right|{|}^2\mathsf{\text{d}}x.\end{array}$$

From which, passing the limit μ → ∞, we get

$$\begin{array}{c}\underset{-a}{\overset{0}{\int }}\mathsf{\text{d}}t\underset{B\left(\sigma \right)}{\int }\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x\\ \le c\left(\nu ,K,\lambda ,\sigma ,m,n\right)h^2\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^{\frac{1+\vartheta }{2}}+\underset{-b^{*}}{\overset{0}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\}\\ +c\left(\nu ,K,\sigma ,m,n\right)\underset{-b^{*}}{\overset{0}{\int }}\mathsf{\text{d}}t\underset{B\left(2\sigma \right)}{\int }{\psi }^{2m}{\rho }_{\mu }^2\left|\right|{\tau }_{i,h}{D}^{\prime }u\left|{|}^2\right|\left|D^{″}u\right|{|}^2\mathsf{\text{d}}x.\end{array}$$

(3.32)

Let us now estimate the last term in (3.32). Using Hölder inequality, applying Theorem 2.2 (for p = 4, $B\left(\frac{5}{2}\sigma \right)$ instead of B(σ) and $t=\frac{4}{5}$) and formula (3.13), for every |h| < h₀, it follows

$$\begin{array}{c}\underset{B\left(2\sigma \right)}{\int }\left|\right|{\tau }_{i,h}{D}^{\prime }u\left|{|}^2\right|\left|D^{″}u\right|{|}^2\mathsf{\text{d}}x\le {\left(\underset{B\left(2\sigma \right)}{\int }\left|\right|{\tau }_{i,h}{D}^{\prime }u|{|}^4\mathsf{\text{d}}x\right)}^{\frac{1}{2}}{\left(\underset{B\left(2\sigma \right)}{\int }\left|\right|D^{″}u|{|}^4\mathsf{\text{d}}x\right)}^{\frac{1}{2}}\\ \le \left|h{|}^2\right|\left|D^{″}u\right|{|}_{0,4,B\left(\frac{5}{2}\sigma \right)}^2\left|\right|D^{″}u|{|}_{0,4,B\left(2\sigma \right)}^2\le |h{|}^2{\left|u\right|}_{m,4,B\left(\frac{5}{2}\sigma \right)}^4.\end{array}$$

Integrating in (-b*, 0), from (3.32), it follows

$$\begin{array}{c}\underset{-a}{\overset{0}{\int }}\mathsf{\text{d}}t\underset{B\left(\sigma \right)}{\int }\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x\\ \le c\left(\nu ,K,U,\lambda ,\sigma ,a,b,m,n\right)|h{|}^2\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^{\frac{1+\vartheta }{2}}+\underset{-b}{\overset{0}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\}.\end{array}$$

(3.33)

If $h_0\le \left|h\right|<\frac{\sigma }{2}$, for every i = 1, 2, ..., n we easily obtain

$$\begin{array}{l}{\displaystyle \underset{-a}{\overset{0}{\int }}\text{d}t{\int }_{B\sigma )}\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\text{d}x}\le 4{\displaystyle \underset{-a}{\overset{0}{\int }}\text{d}}t{\displaystyle \underset{B(3\sigma )}{\int }|}\left|D^{″}u\right|{|}^2\text{d}x\le 4\frac{h^2}{h_0^2}{\displaystyle \underset{-a}{\overset{0}{\int }}\text{d}t{\displaystyle \underset{B(3\sigma )}{\int }|}\left|D^{″}u\right|{|}^2\text{d}x}\\ \le c(\nu ,K,U,\lambda ,\sigma ,a,b,m,n)h^2{\displaystyle \underset{-b}{\overset{0}{\int }}|u{|}_{m,B(3\sigma )}^2\text{d}t}\\ \le c(\nu ,K,U,\lambda ,\sigma ,a,b,m,n)|h{|}^2\left\{1+{\displaystyle \sum _{\left|\alpha \right|<m}{\left({\displaystyle \underset{-b}{\overset{0}{\int }}|}\left|f^{\alpha }\right|{|}_{0,B(3\sigma )}\text{d}t\right)}^{\frac{1+\vartheta }{2}}}+{\displaystyle \underset{-b}{\overset{0}{\int }}|}u{|}_{m,B(3\sigma )}^2\text{d}t\right\}.\end{array}$$

It is then proved, for every $\left|h\right|<\frac{\sigma }{2}$ and every i ∈ {1, 2, ..., n}, that

$$\begin{array}{c}\underset{-a}{\overset{0}{\int }}\mathsf{\text{d}}t\underset{B\left(\sigma \right)}{\int }\left|\right|{\tau }_{i,h}{D}^{″}u|{|}^2\mathsf{\text{d}}x\\ \le c\left(\nu ,K,U,\lambda ,\sigma ,a,b,m,n\right)|h{|}^2\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^{\frac{1+\vartheta }{2}}+\underset{-b}{\overset{0}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\},\end{array}$$

applying Theorem 2.1, it follows

$$u\in L^2\left(-a,0,H^{m+1}\left(B\left(\sigma \right),{ℝ}^N\right)\right)$$

and

$$\begin{array}{c}\underset{-a}{\overset{0}{\int }}|u{|}_{m+1,B\left(\sigma \right)}^2\mathsf{\text{d}}t\\ \le c\left(\nu ,K,U,\lambda ,\sigma ,a,b,m,n\right)\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^{\frac{1+\vartheta }{2}}+\underset{-b}{\overset{0}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\}.\end{array}$$

(3.34)

Finally we have to prove that u ∈ H¹(-a, 0, L²(B(σ), ℝ^N)) and inequality (3.15).

From inequality (3.8), we have

$$\begin{array}{c}\underset{-a}{\overset{0}{\int }}\mathsf{\text{d}}t\underset{B\left(\sigma \right)}{\int }\left|\right|D^{″}u|{|}^4\mathsf{\text{d}}x\\ \le c\left(\nu ,K,U,\lambda ,\sigma ,a,b,m,n\right)\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^{\frac{1+\vartheta }{2}}+\underset{-b}{\overset{0}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\}\end{array}$$

(3.35)

then we have

$$D^{″}u\in L^4\left(B\left(\sigma \right)\times \left(-a,0\right),{ℛ}^{″}\right).$$

(3.36)

Moreover, bearing in mind that, for |α| < m, a^α(X, p) satisfies (3.4), and for |α| = m, a^α(X, p) satisfies (3.5), we have

$$\begin{array}{c}\hfill D^{\alpha }{a}^{\alpha }\left(X,p\right)\in L^2\left(B\left(\sigma \right)\times \left(-a,0\right),R^N\right)\forall \alpha :\left|\alpha \right|\le m\hfill \end{array}$$

(3.37)

Recalling the definition of weak solution, for every $\varphi \in C_0^{\infty }\left(Q,{ℝ}^N\right)$, proceeding as in [22], we have

$$\underset{-a}{\overset{0}{\int }}\mathsf{\text{d}}t\underset{B\left(\sigma \right)}{\int }\left(u\left|\frac{\partial \phi }{\partial t}\right.\right)\mathsf{\text{d}}x=\sum _{\left|\alpha \right|\le m}\underset{-a}{\overset{0}{\int }}\mathsf{\text{d}}t\underset{B\left(\sigma \right)}{\int }\left(D^{\alpha }{a}^{\alpha }\left(X,Du\right)|\phi \right)\mathsf{\text{d}}x,$$

(3.38)

and, bearing in mind (3.37), we obtain that

$$\exists \frac{\partial u}{\partial t}\in L^2\left(B\left(\sigma \right)\times \left(-a,0\right),{ℝ}^N\right).$$

(3.39)

From (3.4), (3.5) and (3.38), we get

$$\begin{array}{c}\underset{-a}{\overset{0}{\int }}\mathsf{\text{d}}t\underset{B\left(\sigma \right)}{\int }{∥\frac{\partial u}{\partial t}∥}^2\mathsf{\text{d}}x\\ \le c\left(\nu ,K,U,\lambda ,\sigma ,a,b,m,n\right)\left\{1+\sum _{\left|\alpha \right|<m}{\left(\underset{-b}{\overset{0}{\int }}\left|\right|f^{\alpha }|{|}_{0,B\left(3\sigma \right)}\mathsf{\text{d}}t\right)}^2+\underset{-b}{\overset{0}{\int }}|u{|}_{m,B\left(3\sigma \right)}^2\mathsf{\text{d}}t\right\}.\end{array}$$

The last inequality and (3.34) allows us to conclude the proof.   ■