Entropy Formulation for Triply Nonlinear Degenerate Elliptic-Parabolic-Hyperbolic Equation with Zero-Flux Boundary Condition

Mohamed Karimou Gazibo

Laboratoire de Mathématiques Fondamentales et Applications (LMFA), Département de Mathématiques, Ecole Normale Supérieure, Université Abdou Moumouni de Niamey, Niamey, Niger.

DOI: 10.4236/jamp.2023.114063                

Abstract

In this note, we investigated existence and uniqueness of entropy solution for triply nonlinear degenerate parabolic problem with zero-flux boundary condition. Accordingly to the case of doubly nonlinear degenerate parabolic hyperbolic equation, we propose a generalization of entropy formulation and prove existence and uniqueness result without any structure condition.

Keywords

Degenerate Elliptic-Parabolic Hyerbolic Equation, Zero-Flux Boundary Condition, Structure Condition, Entropy Formulation, Multi-Step Approximation, Nonlinear Semigroup Theories, Integral and Mild Solution

1. Introduction

Let Ω be a bounded open set of ${ℝ}^{\mathcal{l}}$ with a Lipschitz boundary ∂Ω and $\eta$ the unit normal to ∂Ω outward to Ω. We consider the triply nonlinear degenerate parabolic-elliptic-hyperbolic problem with zero-flux boundary condition:

$\{\begin{array}{l}b{\left(u\right)}_t+\text{div}f\left(u\right)-\Delta \varphi \left(u\right)+\psi \left(u\right)=0\text{in}Q=\left(0,T\right)\times \Omega \\ b\left(u\right)\left(t=0,x\right)=b_0\left(x\right)\text{in}\Omega \\ \left(f\left(u\right)-\nabla \varphi \left(u\right)\right)\cdot \eta =0\text{on}\Sigma =\left(0,T\right)\times \partial \Omega \end{array}$ (P)

The particularity of this problem is its strong degeneracy. For practical reason and physical consideration, we suppose that [0, 1] will be the invariant domain of solution of (P) and that there exist two particular values of the unknown u. We denote by $u_s$ and $u_c$ (with $0<u_c\le u_s<1$ ) such that $\varphi \left(u\right)$ (resp $b\left(u\right)$ ) is strictly increasing only on $\left[u_c,1\right]$ (resp $\left[0,u_c\right]\cup \left[u_s,1\right]$ ) otherwise it has a flat region (see Figure 1). Then, problem (P) is of mixed elliptic parabolic hyperbolic type with absorption term, and thus combines the difficulties related to nonlinear conservation laws with those related to nonlinear degenerate diffusion equations. We refer to Kruzkov [1] for the case of hyperbolic problem ( $\varphi \equiv 0$ ) and Carrillo [2] for degenerate parabolic problem to inspiration. We need a notion of solution which is sufficient to deal with existence and uniqueness. One consequence is that the notion of weak solution generally leads to non-uniqueness, unless $\varphi$ is strictly increasing. It is necessary to adopt an entropy formulation. The notion of entropy solution we use is adapted from the founding paper of [3] in the case where $u_s=u_c$ and $\varphi =0$ . Several authors have studied the degenerate equation type we consider.

Some of these authors ( [3] [4] ) proved existence and uniqueness under the hypothesis that the convection flux f is a Lipschitz-continuous function and required that

$f\left(0\right)=0,f\left(1\right)=0.$ (H)

This hypothesis is necessary to obtain the solution in [0, 1] if the initial datum $u_0$ belongs in [0, 1] in the sense that $b\circ u_0\left(x\right)=b_0\left(x\right)$ and the hypothesis (H) is below. The main idea in the paper is to keep this hypothesis but we impose that initial datum $u_0$ belongs to [0, 1]. We suppose that the function b satisfies:

$\exists \left(\underset{\_}{\alpha },\bar{\alpha }\right)\in {ℝ}_{+}^2\text{suchthat}b\left(u\right)\in \left[\underset{\_}{\alpha },\bar{\alpha }\right].$ (1.1)

A simple choice is to take $b\left(0\right)=\underset{\_}{\alpha }$ and $b\left(1\right)=\bar{\alpha }$ . Our assumptions on b and $\varphi$ do not concern the case where the structure condition

$b\left(u\right)=b\left(v\right)\Rightarrow \varphi \left(u\right)=\varphi \left(v\right).$ (S.C)

holds. The presence of the absorption term in the equation requires us to assume it to be increasing. Further, we suppose:

$\psi \left(r\right)=sign\left(r\right)g\left(r\right)\text{where}g\text{isapositiffunction}\text{.}$ (1.2)

The case of triply nonlinear problems of the form (P) has been first addressed by Ouaro and Touré (see [5] ) and the references therein) and Ouaro [6] .

Figure 1. Convection and diffusion flow.

Well-posedness results are obtained in dimension $\mathcal{l}=1$ , under very general coercivity conditions; see also the works of Bénilan and Touré ( [7] and the references therein). Andreianov and Wittbold investigate in [8] about the continuous dependence of the solution of a degenerate elliptic-parabolic equation without structure condition related to b and f. They prove existence by passing to bi-monotonicity and penalization method as in [9] . Otherwise, in [10] , Andreianov et al. obtain a general continuous dependence result on data for our kind of triply nonlinear problem with help of structure condition. They showed similar result for the degenerate elliptic problem, which corresponds to the case of $b\equiv 0$ and general non-decreasing surjective $\psi$ . In our case the function $\psi$ is bounded continuous and strictly increasing. If (S.C) fails, the convergence of approximate solutions to (P) is known for a particular monotone approximation method developed by Ammar and Wittbold [9] . This approach leads to an existence result which bypasses (S.C). Notice that some essential arguments of uniqueness result in this works are specific to the case $\mathcal{l}=1$ . For Neumann boundary condition also called zero-flux boundary condition, it is easy to prove uniqueness of solutions such that the boundary condition is satisfied in the sense of strong boundary trace of the normal component of the flux $\left(f\left(u\right)-\nabla \varphi \left(u\right)\right)$ . Unfortunately, we are able to establish this additional solution regularity only for the stationary problem (S) associated to (P) and only in the case of one space dimension.

The paper is divided in three parts, in Section 2 we generalized the notion of entropy solution of paper [3] where $u_c<u_s$ and [4] in the pure hyperbolic case. In Section 3, we first prove existence and after uniqueness of entropy solution.

2. Formulation of Entropy Solution

2.1. Definition of Entropy Solution

We need the notion of weak solution for (P) with additional “entropy” conditions.

Definition 2.1. A measurable function u taking values on [0, 1] is called an entropy solution of the initial-boundary value problem (P) if satisfying the following conditions:

$\varphi \left(u\right)\in L^2\left(0,T;H^1\left(\Omega \right)\right)$ and $\forall \xi \in {\mathcal{C}}^{\infty }\left(\left[0,T\right)\times {ℝ}^{\mathcal{l}}\right)$ , with $\xi \ge 0$

${\displaystyle \underset{0}{\overset{T}{\int }}\langle b{\left(u\right)}_t,\xi \rangle \text{d}t}-{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\left(f\left(u\right)-\nabla \varphi \left(u\right)\right)\cdot \nabla \xi \text{d}x\text{d}t}}+{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\psi \left(u\right)\xi \text{d}x\text{d}t}}=0.$ (2.1)

$\varphi \left(u\right)\in L^2\left(0,T;H^1\left(\Omega \right)\right)$ and $\forall k\in \left[0,u_c\right]\cup \left[u_s,1\right]$ , $\forall \xi \in {\mathcal{C}}^{\infty }\left(\left[0,T\right)\times {ℝ}^{\mathcal{l}}\right)$ , with $\xi \ge 0$ ,

$\begin{array}{l}{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\left\{\left|b\left(u\right)-b\left(k\right)\right|{\xi }_t+sign\left(u-k\right)\left(f\left(u\right)-f\left(k\right)-\nabla \varphi \left(u\right)\right)\cdot \nabla \xi \right\}\text{d}x\text{d}t}}\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\partial \Omega }{\int }\left|f\left(k\right)\cdot \eta \left(x\right)\right|\xi \left(t,x\right)\text{d}{\mathcal{H}}^{\mathcal{l}-1}\left(x\right)\text{d}t}}-{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sign\left(u-k\right)\psi \left(u\right)\xi \text{d}x\text{d}t}}\\ +{\displaystyle \underset{\Omega }{\int }\left|b_0-b\left(k\right)\right|\xi \left(0,x\right)\text{d}x}\ge 0.\end{array}$ (2.2)

Here $\mathcal{H}$ represents the $\left(\mathcal{l}-1\right)$ -dimensional Hausdorff measure and $\langle .,.\rangle$ is the duality pairing between $L^2\left(0,T,{\left(H^1\left(\Omega \right)\right)}^{\star }\right)+L^1\left(Q\right)$ and $L^2\left(0,T,H^1\left(\Omega \right)\right)\cap L^{\infty }(\; Q\; )$

It is well known that the distributional derivative $b{\left(u\right)}_t$ of $b\left(u\right)$ can be identified with an element of the space $L^2\left(0,T,{\left(H^1\right)}^{\star }\left(\Omega \right)\right)+L^1\left(Q\right)$ . More exactly, we have

${\displaystyle \underset{0}{\overset{T}{\int }}\langle b{\left(u\right)}_t,\xi \rangle \text{d}t}={\displaystyle \underset{Q}{\int }b}\left(u\right){\xi }_t-{\displaystyle \underset{\Omega }{\int }{b}_0}\xi \left(0,x\right)\text{d}x$ (2.3)

for all $\xi \in L^2\left(0,T;H^1\left(\Omega \right)\right)$ with ${\xi }_t\in L^{\infty }\left(Q\right)$ and $\xi \left(T,x\right)=0$ .

We obtain notions of entropy sub-solution and entropy super-solution respectively if we replace (2.2) by one of the following inequalities

$\begin{array}{l}{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\left\{{\left(b\left(u\right)-b\left(k\right)\right)}^{+}{\xi }_t+sig{n}^{+}\left(u-k\right)\left(f\left(u\right)-f\left(k\right)-\nabla \varphi \left(u\right)\right)\cdot \nabla \xi \right\}\text{d}x\text{d}t}}\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }{\left(f\left(k\right)\cdot \eta \left(x\right)\right)}^{+}\xi \left(t,x\right)\text{d}{\mathcal{H}}^{\mathcal{l}-1}\left(x\right)\text{d}t}}-{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sig{n}^{+}\left(u-k\right)\psi \left(u\right)\xi \text{d}x\text{d}t}}\\ +{\displaystyle \underset{\Omega }{\int }{\left(b_0-b\left(k\right)\right)}^{+}\xi \left(0,x\right)\text{d}x}\ge 0\end{array}$ (2.4)

$\begin{array}{l}{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\left\{{\left(b\left(u\right)-b\left(k\right)\right)}^{-}{\xi }_t+sig{n}^{-}\left(u-k\right)\left(f\left(u\right)-f\left(k\right)-\nabla \varphi \left(u\right)\right)\cdot \nabla \xi \right\}\text{d}x\text{d}t}}\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }{\left(f\left(k\right)\cdot \eta \left(x\right)\right)}^{-}\xi \left(t,x\right)\text{d}{\mathcal{H}}^{\mathcal{l}-1}\left(x\right)\text{d}t}}-{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sig{n}^{-}\left(u-k\right)\psi \left(u\right)\xi \text{d}x\text{d}t}}\\ +{\displaystyle \underset{\Omega }{\int }{\left(b_0-b\left(k\right)\right)}^{-}\xi \left(0,x\right)\text{d}x}\ge 0.\end{array}$ (2.5)

Remark 2.2.

Obviously, a function u is an entropy solution if and only if u is entropy sub-solution and entropy super-solution simultaneously.

Our notion of entropy solution coincide with the Definition of [3] in the case $u_s=u_c$ , $\psi \equiv 0$ and assumption (S.C) is trivially satisfied.

Notice that if u satisfy (2.2), then use (1.1) and (H), we have also u verify (2.1).

Let us stress that, if (H) is satisfy, in particular, the zero-flux boundary condition $\left(f\left(u\right)-\nabla \varphi \left(u\right)\right)\cdot \eta =0$ is verified literally in the weak sense (see for exemple [3] and [11] ). A forthcoming work is to envisage envisage (P) if assumption (H) is dropped. We expect that the boundary condition should be relaxed.

2.2. Dissipative Property

We propose here an essential property of entropy solutions, based on the idea of J. Carrillo [2] .

Proposition 2.3. Let $\xi \in {\mathcal{C}}^{\infty }\left(\left[0,T\right[\times {ℝ}^{\mathcal{l}}\right)$ , $\xi \ge 0$ . Then for all $k\in \left[u_s,1\right]$ , for all $D\in {ℝ}^{\mathcal{l}}$ for all entropy solution u of (P), we have

$\begin{array}{l}{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\left\{\left|b\left(u\right)-b\left(k\right)\right|{\xi }_t+sign\left(u-k\right)\left(f\left(u\right)-f\left(k\right)-\nabla \varphi \left(u\right)-D\right)\cdot \nabla \xi \right\}}}\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\left|\left(f\left(k\right)-D\right)\cdot \eta \left(x\right)\right|\xi \left(t,x\right)\text{d}{\mathcal{H}}^{\mathcal{l}-1}\left(x\right)\text{d}t}}-{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sign\left(u-k\right)\psi \left(u\right)\xi \text{d}x\text{d}t}}\\ +{\displaystyle \underset{\Omega }{\int }\left|b_0-b\left(k\right)\right|\xi \left(0,x\right)\text{d}x}\ge \underset{\sigma \to 0}{\bar{\lim }}\frac{1}{\sigma }{\displaystyle \int {\displaystyle \underset{Q\cap A_{\varphi }^{\sigma }}{\int }\nabla \varphi \left(u\right)}}\cdot \left(\nabla \varphi \left(u\right)-D\right)\xi .\end{array}$ (2.6)

where $A_{\varphi }^{\sigma }=\left\{\left(t,x\right)\in Q\text{suchthat}-\sigma <\varphi \left(u\right)-\varphi \left(k\right)<\sigma \right\}$ .

Proof. The ingredient of the proof of Proposition 2.7 can from firstly, for all $u\in \left[0,u_c\right]\cup \left[u_s,1\right]$ and for all $k\in \left[u_s,1\right]$ , one has:

$sign\left(u-k\right)=sign\left(b\left(u\right)-b\left(k\right)\right)=sign\left(\varphi \left(u\right)-\varphi \left(k\right)\right)$ .

Secondly, $\nabla \varphi \left(u\right)=0$ ${\mathcal{L}}^N$ a.e. on the set $\left\{\left(t,x\right)\in Q\text{suchthat}u\left(t,x\right)\in \left[0,u_c\right]\right\}$ .

Taking as test function in (2.1) $sig{n}_{\sigma }\left(\varphi \left(u\right)-\varphi \left(k\right)\right)\xi$ with $\xi \in C^{\infty }\left(\left[0,T\right[\times {ℝ}^{\mathcal{l}}\right)$ and $sig{n}_{\sigma }$ is the approximation of $sign$ function, using Chain rule (see [12] ) and passing to the limit $\sigma \to 0$ . For more details see [3] . ▄

3. Existence and Uniqueness Result

3.1. Existence of Entropy Solution

The main result of this subsection is the following theorem:

Theorem 3.1. Assume that (1.1), (1.2), (H) holds. Then there exists an entropy solution u for (P).

3.1.1. Bi-Monotonicity Approach

Because of we are not in the case where (S.C) holds, we use the particular multi-step approximation approach of Ammar and Wittbold (see [9] ).

Theorem 3.2. Let ${\left(b_l,{\varphi }_l,{\psi }_{m,n},u_0^n\right)}_{\left(l,m,n\right)\in {\mathbb{N}}^{*3}}$ , a sequences converging to $\left(b,\varphi ,\psi ,u_0\right)$ in the following sense:

$b_l,{\varphi }_l,{\psi }_{m,n}\text{convergepointwiseto}b,\varphi ,\psi \text{as}l,m,n\to +\infty$ (3.1)

$u_0^n\in \left[0,1\right]\text{convergea}\text{.eto}{u}_0$ (3.2)

$b_l\left(u_0^n\right)\text{convergeto}b\left(u_0\right)\text{in}{L}^1\left(\Omega \right).$ (3.3)

There exists a weak solution $u_{m,n}^l$ of ${\left(\text{P}\right)}_{m,n}^l$ the analogue of (P) with corresponding data $\left(b_l,{\varphi }_l,{\psi }_{m,n},u_0^n\right)$ .

Proof. We consider the particular multi-step approximation approach of Ammar and Witt bold [9] . We replace b by $b_l:=b+\frac{1}{l}Id$ , $\varphi$ by ${\varphi }_l:=\varphi +\frac{1}{l}Id$ and $\psi$ , by ${\psi }_{m,n}:=\psi +\frac{1}{n}I{d}^{+}-\frac{1}{m}I{d}^{-}$ . Hence ${\varphi }_l,{\varphi }_l^{-1},b_l,b_l^{-1}$ are Lipschitz continuous on $ℝ$ .

We obtain the following equation:

$b_l{\left(u_{m,n}^l\right)}_t+\text{div}f\left(u_{m,n}^l\right)-\Delta {\varphi }_l\left(u_{m,n}^l\right)+{\psi }_{m,n}\left(u_{m,n}^l\right)=0\text{in}Q=\left(0,T\right)\times \Omega$

Take $v_{m,n}^l=b_l\left(u_{m,n}^l\right)$ , hence $b_l$ is invertible, one puts the problem into the doubly non-linear framework then we obtain the following problem

${\left(v_{m,n}^l\right)}_t+\text{div}\tilde{f}\left(v_{m,n}^l\right)-\Delta \widetilde{{\varphi }_l}\left(v_{m,n}^l\right)+{\tilde{\psi }}_{m,n}\left(v_{m,n}^l\right)=0\text{in}Q$

where $\tilde{f}=f°b_l^{-1}$ ; $\widetilde{{\varphi }_l}=\varphi °b_l^{-1}$ and ${\tilde{\psi }}_{m,n}={\psi }_{m,n}°b_l^{-1}$ . Using classical methods (cf. Andreianov and Gazibo [3] ), one shows that there exists a weak solution $u_{m,n}^l$ for the corresponding problem ${\left(\text{P}\right)}_{m,n}^l$ . ▄

Theorem 3.3. Let $u_{m,n}^l$ be the weak solution of ${\left(\text{P}\right)}_{m,n}^l$ the analogue of (P) with corresponding data $\left(b_l,{\varphi }_l,{\psi }_{m,n},u_0^n\right)$ . Then $u_{m,n}^l$ is also entropy solution of ${\left(\text{P}\right)}_{m,n}^l$ in the sense of Definition 2.1 and converge to u entropy solution of (P) in $L^{\infty }$ weakly star up to a subsequence. Furthermore:

${\varphi }_l\left(u_{m,n}^l\right)\to \varphi \left(u\right)\text{in}{L}^1\left(Q\right)$ (3.4)

$b_l\left(u_{m,n}^l\right)\to b\left(u\right)\text{in}{L}^1\left(Q\right)$ (3.5)

${\psi }_{m,n}\left(u_{m,n}^l\right)\to \psi \left(u\right)\text{in}{L}^1\left(Q\right).$ (3.6)

3.1.2. A priori Estimates

Lemma 3.4. Let ${\left(b_l,{\varphi }_l,{\psi }_{m,n},u_0^n\right)}_{\left(l,m,n\right)\in {\mathbb{N}}^{*3}}$ , be a sequence of data satisfying the as sumption of theorem 3.2. Assume that the corresponding data $\left(b,\varphi ,\psi ,u_0\right)$ verifies (1.1), (1.2), (H). Let $u_{m,n}^l$ be an entropy solution of ${\left(\text{P}\right)}_{m,n}^l$ then there exist $L>0$ such that:

$0\le u_{m,n}^l\le 1$ (3.7)

${\Vert {\varphi }_l\left(u_{m,n}^l\right)\Vert }_{L^2\left(0,T;H^1\left(\Omega \right)\right)}\le L$ (3.8)

${\Vert {\varphi }_l\left(u_{m,n}^l\right)\Vert }_{L^2\left(Q\right)}\le L$ (3.9)

${\Vert {\varphi }_l\left(u_{m,n}^l\right){\psi }_{m,n}\left(u_{m,n}^l\right)\Vert }_{L^1\left(Q\right)}\le L.$ (3.10)

Proof. Since $u_{m,n}^l$ is entropy solution, it is also entropy subsolution and entropy super solution of ${\left(\text{P}\right)}_{m,n}^l$ . Take $k=0$ in (2.5) and $\xi \left(t,x\right)=\xi \left(t\right)$ and (H) and (1.2) then

${\displaystyle {\int }_0^T\left({\displaystyle {\int }_{\Omega }\left({\left(b_l\left(u_{m,n}^l\right)-b_l\left(0\right)\right)}^{-}-{\left(b_l\left(u_n^0\right)-b_l\left(0\right)\right)}^{-}\right)\text{d}x}\right){\xi }_t\text{d}t}\ge 0$ (3.11)

Let us introduce the function

$H\left(t\right)={\displaystyle \underset{\Omega }{\int }\left({\left(b_l\left(u_{m,n}^l\right)-b_l\left(0\right)\right)}^{-}-{\left(b_l\left(u_n^0\right)-b_l\left(0\right)\right)}^{-}\right)\text{d}x}\text{if}t\in \left(0,T\right)$

Since (1.1), we have $u_{m,n}^l\ge 0$ . In the same way, $u_{m,n}^l$ satisfy (2.4), take now $k=1$ , we prove that $u_{m,n}^l\le 1$ i.e.

$0\le u_{m,n}^l\le 1$ (3.12)

We use the test function ${\varphi }_l\left(u_{m,n}^l\right)$ in the weak formulation of ${\left(\text{P}\right)}_{m,n}^l$ . The duality product between ${\varphi }_l\left(u_{m,n}^l\right)\in L^2\left(0,T;H^1\left(\Omega \right)\right)\cap L^{\infty }\left(Q\right)$ and ${\partial }_t{b}_l\left(u_{m,n}^l\right)\in L^2\left(0,T,{\left(H^1\right)}^{\star }\left(\Omega \right)\right)+L^1\left(Q\right)$ is treated via the standard chain rule argument

$\begin{array}{l}{\displaystyle \underset{0}{\overset{t}{\int }}\langle b_l{\left(u_{m,n}^l\right)}_t,{\varphi }_l\left(u_{m,n}^l\right)\rangle \text{d}t}-{\displaystyle \underset{0}{\overset{t}{\int }}{\displaystyle \underset{\Omega }{\int }\left(f\left(u_{m,n}^l\right)-\nabla {\varphi }_l\left(u_{m,n}^l\right)\right)\cdot \nabla {\varphi }_l\left(u_{m,n}^l\right)\text{d}x\text{d}t}}\\ +{\displaystyle \underset{0}{\overset{t}{\int }}{\displaystyle \underset{\Omega }{\int }{\psi }_{m,n}}}\left(u_{m,n}^l\right){\varphi }_l\left(u_{m,n}^l\right)\text{d}x\text{d}t=0.\end{array}$ (3.13)

Take $B_l\left(z\right)={\displaystyle {\int }_0^z{\varphi }_l}\left(s\right)\text{d}{b}_l\left(s\right)$ . Since $B_l\left(z\right)\le b_l\left(z\right){\varphi }_l\left(z\right)$ , we obtain the inequality

${\displaystyle \underset{\Omega }{\int }{B}_l}\left(u_{m,n}^l\left(t,.\right)\right)+{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }{\psi }_{m,n}}}\left(u_{m,n}^l\right){\varphi }_l\left(u_{m,n}^l\right)+C{\left|\nabla {\varphi }_l\left(u_{m,n}^l\right)\right|}^2\le {\displaystyle \underset{\Omega }{\int }{b}_l}\left(u_n^0\right){\varphi }_l\left(u_n^0\right).$

with some $C>0$ independent of n. Note that the functions $b_l,{\varphi }_l$ are locally uniformly bounded because they are monotone and converge pointwise to $b,\varphi$ respectively.

Therefore the right-hand side of the above inequality is bounded uniformly in $l$ , thanks to (3.12) and the uniform bounds on the data $u_n^0$ in $L^{\infty }\left(\Omega \right)$ . The uniform estimate of the left-hand side follows. We then estimate $\Vert {\varphi }_l\left(u_{m,n}^l\right)\Vert$ in $L^2\left(Q\right)$ by the Poincaré inequality and ${\Vert {\varphi }_l\left(u_{m,n}^l\right)\psi \left(u_{m,n}^l\right)\Vert }_{L^1}$ follows. ▄

Lemma 3.5. Let $u_{m,n}^l$ be the weak solution of ${\left(\text{P}\right)}_{m,n}^l$ the analogue of (P) with corresponding data $\left(b_l,{\varphi }_l,{\psi }_{m,n},u_0^n\right)$ . For $\tau >0$ , we have

${\displaystyle \underset{0}{\overset{T-\tau }{\int }}{\displaystyle \underset{\Omega }{\int }\left|{\varphi }_l\left(u_{m,n}^l\right)\left(t+\tau ,x\right)-{\varphi }_l\left(u_{m,n}^l\right)\left(t,x\right)\right|}}\le \omega \left(\tau \right)$ (3.14)

Proof. Let $\tau >0$ , Multiplying the first equation of ${\left(\text{P}\right)}_{m,n}^l$ by $\xi \in H^1\left(\Omega \right)\cap L^{\infty }\left(\Omega \right)$ and integer in $\left(t,t+\tau \right)\subset \left(0,T\right)$ we get:

$\begin{array}{l}{\displaystyle \underset{\Omega }{\int }\left(b_l\left(u_{m,n}^l\left(x,t+\tau \right)\right)-b_l\left(u_{m,n}^l\left(x,t\right)\right)\right)\xi }+{\displaystyle \underset{t}{\overset{t+\tau }{\int }}{\displaystyle \underset{\Omega }{\int }\left(f\left(u_{m,n}^l\right)-\nabla {\varphi }_l\left(u_{m,n}^l\right)\right)\cdot \nabla \xi }}\\ +{\displaystyle \underset{t}{\overset{t+\tau }{\int }}{\displaystyle \underset{\Omega }{\int }\psi \left(u_{m,n}^l\right)\xi }}=0.\end{array}$

Take now $\xi ={\varphi }_l\left(u_{m,n}^l\right)\left(t+\tau ,x\right)-{\varphi }_l\left(u_{m,n}^l\right)\left(t,x\right)$ and integrate in t. By Fubini theorem and estimates of Lemma 3.4, it appear a factor $\tau$ in the right hand side and we get

${\displaystyle \int {\displaystyle \underset{Q_{\tau }=\left(0,T-\tau \right)\times \Omega }{\int }\left|h_{\tau }\left(b_l\left(u_{m,n}^l\right)\left(t,x\right)\right)\right|\left|h_{\tau }\left({\varphi }_l\left(u_{m,n}^l\right)\left(t,x\right)\right)\right|}}\le C\tau$ (3.15)

with $h_{\tau }\left(g\right)\left(t,x\right)g\left(t+\tau ,x\right)-g\left(t,x\right)$ . Here C is a constant independent of $l,m,n$ . Take $\bar{\varphi }={\varphi }_l°b_l^{-1}$ which is continuous function, let $\bar{w}$ the modulus of continuity of $\bar{\varphi }$ on $\left[0,1\right]$ and $\bar{W}$ be its inverse and set $W\left(r\right)=r\bar{W}\left(r\right)$ and w be the inverse of W, notice that $w\left(0\right)=0$ . Denote by $s=b_l\left(u_{m,n}^l\right)\left(t+\tau ,x\right)$ and $r=b_l\left(u_{m,n}^l\right)\left(t,x\right)$ . Then

$\begin{array}{c}{\displaystyle \int {\displaystyle \underset{Q_{\tau }}{\int }\left|\bar{\varphi }\left(s\right)-\bar{\varphi }\left(r\right)\right|}}={\displaystyle \int {\displaystyle \underset{Q_{\tau }}{\int }w\left(W\left(\left|\bar{\varphi }\left(s\right)-\bar{\varphi }\left(r\right)\right|\right)\right)}}\\ \le \left|Q_{\tau }\right|w\left(\frac{1}{\left|Q_{\tau }\right|}{\displaystyle \int {\displaystyle \underset{Q_{\tau }}{\int }W}}\left(\left|\bar{\varphi }\left(s\right)-\bar{\varphi }\left(r\right)\right|\right)\right)\end{array}$

Since

$\left|\bar{\varphi }\left(s\right)-\bar{\varphi }\left(r\right)\right|\le \bar{w}\left(\left|s-r\right|\right)$

we have:

$\bar{W}\left(\left|\bar{\varphi }\left(s\right)-\bar{\varphi }\left(r\right)\right|\right)\le \left|s-r\right|$

$\begin{array}{c}W\left(\left|\bar{\varphi }\left(s\right)-\bar{\varphi }\left(r\right)\right|\right)=\bar{W}\left(\left|\bar{\varphi }\left(s\right)-\bar{\varphi }\left(r\right)\right|\right)\left|\bar{\varphi }\left(s\right)-\bar{\varphi }\left(r\right)\right|\\ \le \left|s-r\right|\left|\bar{\varphi }\left(s\right)-\bar{\varphi }\left(r\right)\right|\\ =\left|h_{\tau }\left(b_l\left(u_{m,n}^l\right)\right)\left(t,x\right)\right|\left|h_{\tau }\left({\varphi }_l\left(u_{m,n}^l\right)\left(t,x\right)\right)\right|.\end{array}$

Therefore (3.15) implies

$\begin{array}{c}{\displaystyle \int {\displaystyle \underset{Q_{\tau }}{\int }\left|h_{\tau }\left({\varphi }_l\left(u_{m,n}^l\right)\right)\left(t,x\right)\right|}}\le \left|Q_{\tau }\right|w\left(\frac{1}{\left|Q_{\tau }\right|}{\displaystyle \int {\displaystyle \underset{Q_{\tau }}{\int }\left|h_{\tau }\left(u_{m,n}^l\right)\left(t,x\right)\right|\left|h_{\tau }\left({\varphi }_l\left(u_{m,n}^l\right)\left(t,x\right)\right)\right|}}\right)\\ =\left|Q_{\tau }\right|w\left(\frac{1}{\left|Q_{\tau }\right|}\tau \right)\le Cw\left(C\tau \right)\end{array}$ (3.16)

the left-hand side of (3.16) tends to zero as $\tau \to 0$ , we deduce (3.14). ▄

The proof of Theorem 3.1 is a direct consequence of Theorem 3.2 and Theorem 3.3.

Since, we have establish the proof of theorem 3.2, let us demonstrate Theorem 3.3.

Proof of Theorem 3.3. There exists a function $u_{m,n}^k$ constructed by means of the nonlinear semigroup theory (see, e.g., [3] [13] ), such that $v_k^{m,n}=b^k\left(u_{m,n}^k\right)$ is the unique integral solution to the abstract evolution problem associated with ${\left(P\right)}_{m,n}^k$ (here and below, we refer to Andreianov and Gazibo [3] , Ammar and Wittbold [9] , Ammar and Redwane [14] for details). One then shows that $u_{m,n}^k$ coincides with the unique entropy solution of ${\left(P\right)}_{m,n}^k$ , the existence of this entropy solution being already shown. Further, the whole set ${\left(u_{m,n}^k\right)}_{k,m,n}$ verifies the uniform a priori estimates of Lemma 3.4 and 3.5. We then pass to the limit in $u_{m,n}^k$ in the following order: first $k\to +\infty$ , then $n\to +\infty ,m\to +\infty$ .

While letting $k\to +\infty$ , we use the fact that ${\psi }_{m,n}^{-1}$ is Lipschitz continuous. The fundamental estimates for the semigroup solutions permit to show that ${\psi }_{m,n}\left(u_{m,n}\right)$ are uniformly continuous on $\left(0,T\right)$ with values in $L^1\left(\Omega \right)$ ; thus we get the strong precompactness of ${\left(u_{m,n}^k\right)}_k\in L^1\left(Q\right)$ . Thus, up to a subsequence, $u_{m,n}^k$ converge to $u_{m,n}$ which is an entropy solution of problem ${\left(\text{P}\right)}_{m,n}^l$ corresponding to the data $\left(b,{\psi }_{m,n},\varphi ,u_0\right)$ . Finally, we use the inequalities $u_{m+1,n}\le u_{m,n}\le u_{m,n+1}$ which follow readily form the comparison principle (this come from uniqueness of ${\left(P\right)}_{m,n}^k$ ). The monotonicity argument yields the strong convergence of $u_{m,n}$ . Passing to the limit in $u_{m,n}$ we conclude that the limit u is an entropy solution of the original problem (P) (one can use Lemma 3.6 and 3.7 of [3] ). ▄

3.2. Uniqueness of Entropy Solution in One Space Dimension

3.2.1. Stationary Problem

Let us stress that to our knowledge the problem of uniqueness is still open in multiple space dimensions. The definition of strong traces of the solution with respect to the lateral boundary of the domain Ω is possible if for example the diffusion term $\varphi \left(u\right)$ is such that $f\left(u\right)-\nabla \varphi \left(u\right)$ is continuous up to the boundary $\partial \Omega$ . If there existed “sufficiently many” solutions (in the sense of [3] , [4] , see Definition 3.9 below) having this regularity, uniqueness would follow. Unfortunate, we could obtain this regularity for this moment only for the stationary problem associate to (P) and in one space dimension.

Now, we consider the stationary problem associated to (P)

$\{\begin{array}{l}b{\left(u\right)}_t+\text{div}f\left(u\right)-\Delta \varphi \left(u\right)+\psi \left(u\right)=s\text{in}\Omega \\ \left(f\left(u\right)-\nabla \varphi \left(u\right)\right)\cdot \eta =0\text{on}\partial \Omega \end{array}$ (S)

where $s\in L^{\infty }\left(\Omega \right)$ .

Remark 3.6.

1) If, we suppose that $\left(b+\psi \right)\left(u\right)$ is bijective, then performing a change of the unknown one puts the problem into the doubly nonlinear framework in the form $u+\text{div}\left(f\left(u\right)-\nabla \varphi \left(u\right)\right)=s$ . Existence and uniqueness follows (see [3] ).

2) If $b\left(u\right)$ independent of t and u is solution of (S) it is also solution of (P) with source term $s-b\left(u\right)$ . Then, we can deduce from Definition 2.1 and Proposition 2.3 their equivalent form for the stationary problem.

Definition 3.7. A measurable function u taking values in [0, 1] is an entropy solution of (S) if u is a weak solution of (S) and $\varphi \left(u\right)\in H^1\left(\Omega \right)$ and for all $\xi \in {\mathcal{C}}^{\infty }{\left({ℝ}^{\mathcal{l}}\right)}^{+}$ , $\forall k\in \left[0,u_c\right]\cup \left[u_s,1\right]$ ,

$\begin{array}{l}{\displaystyle \underset{\Omega }{\int }sign\left(u-k\right)\left(f\left(u\right)-f\left(k\right)-\nabla \varphi \left(u\right)\right)\cdot \nabla \xi \text{d}x}+{\displaystyle \underset{\partial \Omega }{\int }\left|f\left(k\right)\cdot \eta \left(y\right)\right|\xi \text{d}{\mathcal{H}}^{\mathcal{l}-1}\left(x\right)}\\ -{\displaystyle \underset{\Omega }{\int }sign\left(u-k\right)\left(b\left(u\right)+\psi \left(u\right)-s\right)\xi \text{d}x}\ge 0.\end{array}$ (3.17)

Proposition 3.7. Let $\xi \in {\mathcal{C}}^{\infty }\left({ℝ}^{\mathcal{l}}\right)$ ; then for all $k\in \left[u_s,1\right]$ , for all $D\in {ℝ}^{\mathcal{l}}$ , for all entropy solution u of (S), we have:

$\begin{array}{l}{\displaystyle \underset{\Omega }{\int }sign}\left(u-k\right)\left(\text{}f\left(u\right)-f\left(k\right)-\nabla \varphi \left(u\right)+D\right)\cdot \nabla \xi \text{d}x+{\displaystyle \underset{\partial \Omega }{\int }\left|\left(f\left(k\right)-D\right)\eta \right|\xi \text{d}{\mathcal{H}}^{\mathcal{l}-1}\left(x\right)}\\ -{\displaystyle \underset{\Omega }{\int }sign}\left(u-k\right)\left(b\left(u\right)+\psi \left(u\right)-s\right)\xi \text{d}x\ge \underset{\sigma \to 0}{\bar{\lim }}\frac{1}{\sigma }{\displaystyle \underset{\Omega \cap A_{\varphi }^{\sigma }}{\int }\nabla }\varphi \left(u\right)\left(\nabla \varphi \left(u\right)-D\right)\xi .\end{array}$ (3.18)

In the next subsection, we give a Definition of so called “trace regular entropy solution”

3.2.2. Trace Regular Entropy Solution

Definition 3.8. An entropy solution is called trace regular solution of (P) if the normal component of the total flux $\left(f\left(u\right)-\nabla \varphi \left(u\right)\right)\cdot \eta$ , has $L^1$ strong trace $\gamma \left(f\left(u\right)-\nabla \varphi \left(u\right)\right)$ at boundary of Lipschitz domain i.e.: for $\xi \in L^{\infty }(\; \partial \; \Omega \; )$

$\underset{\sigma \to 0}{\lim }\frac{1}{\sigma }{\displaystyle \underset{0}{\overset{1}{\int }}{\displaystyle \underset{\partial \Omega }{\int }\xi }}\left|\left(f\left(u\right)-\nabla \varphi \left(u\right)\right)\cdot \eta \left(x\right)-\gamma \left(f\left(u\right)-\nabla \varphi \left(u\right)\right)\cdot \eta \left(x,\sigma \right)\right|\text{d}x\text{d}\tau =0.$ (3.19)

The difficulty is that the regularized zero-flux boundary condition does not permit control over the tangential derivatives (with respect to $\partial \Omega$ ) of the solution. Thus, boundary flux traces of solution seem hard to obtain and we need the concepts of domains with Lipschitz deformable boundaries and traces (see [15] , [16] for more details).

Remark 3.9. Notice that if the normal component of the flux $\left(f\left(u\right)-\nabla \varphi \left(u\right)\right)\cdot \eta$ is continuous function then it satisfy (3.19).

From now on, we will suppose that $\Omega =\left(\text{}a,b\right)$ is a bounded interval of $ℝ$ . We have this property

Proposition 3.10. For all $s\in \left[0,1\right]$ , the problem (S) admits a solution u such that $\left(f\left(u\right)-\varphi {\left(u\right)}_y\right)$ is continuous up the boundary, i.e., $\left(f\left(u\right)-\varphi {\left(u\right)}_y\right)\in \mathcal{C}\left(\left[a,b\right]\right)$ .

Moreover, $f\left(u\right)-\varphi {\left(u\right)}_y$ is zero at $y=a$ and $y=b$ .

Since $\psi$ is bijective, the proof of Proposition 3.11 is identical to the proof of Proposition 4.8 of [3] .

The main result of this section is the following theorem:

Theorem 3.11. Suppose that $\Omega =\left(a,b\right)$ is a bounded interval of $ℝ$ , then (P) admits a unique $b\left(u\right)$ such that u is entropy solution of (P).

3.2.3. Abstract Evolution Problem

We present now the problem (P) under the abstract form of an evolution equation governed by an accretive operator, in order to apply classical results of the nonlinear semigroup theory (see, e.g., [17] ).

Let us define the (possibly multivalued) operator $A_{f,\varphi ,\psi }^b$ by it resolvent

$\left(u,z\right)\in A_{f,\varphi ,\psi }^b=\left\{\begin{array}{l}u\backslash u\text{entropysolutionof}\left(\text{S}\right),\text{with}z=s-b\left(u\right)-\psi \left(u\right)\text{}\text{}\text{}\hfill \end{array}\right\}.$

Consider the abstract equation:

$b{\left(u\right)}_t+A_{f,\varphi ,\psi }^b\left(u\right)\ni 0,b\left(u\right)\left(t=0\right)=b_0$ (E)

For an operator $A:L^1\left(\Omega \right)\to L^1\left(\Omega \right)$ , denote by $R\left(A\right)$ its range, by $D\left(A\right)$ its domain and by $\overline{R\left(A\right)}$ , $\overline{D\left(A\right)}$ their closures in $L^1\left(\Omega \right)$ respectively

Let us stress that for $u\in D\left(A_{f,\varphi ,\psi }^b\right)$ , $f\left(u\right)-\varphi {\left(u\right)}_y\in {\mathcal{C}}_0\left(\left[a,b\right]\right)$ due to Proposition 3.11.

Recall (cf. [17] ) that an operator A is accretive if ${\left[\beta -\hat{\beta },\alpha -\hat{\alpha }\right]}_{L^1\left(\Omega \right)}\ge 0$ for all $\left(\beta ,\alpha \right),\left(\hat{\beta },\hat{\alpha }\right)\in A$ , where for $\beta ,\alpha \in L^1\left(\Omega \right)$ the bracket ${\left[.,.\right]}_{L^1\left(\Omega \right)}$ is defined by ${\left[\beta ,\alpha \right]}_{L^1\left(\Omega \right)}={\displaystyle \underset{\Omega }{\int }sign}\left(\beta \right)\alpha +{\displaystyle \underset{\left[\beta =0\right]}{\int }\left|\alpha \right|}$ .

If A is accretive and $R\left(I+\lambda A\right)=L^1\left(\Omega \right)$ for some $\lambda >0$ , then A is m-accretive.

Proposition 3.13. Let $\left(u,z\right)\in A_{f,\varphi ,\psi }^b$ , $\left(\hat{u},\hat{z}\right)\in A_{f,\varphi ,\psi }^b$ . Then for $\xi \in {\mathcal{C}}^{\infty }{\left(\bar{\Omega }\right)}^{+}$

$\begin{array}{l}{\displaystyle \underset{\Omega }{\int }\left|b\left(u\right)-b\left(\hat{u}\right)\right|\xi \text{d}y}+{\displaystyle \underset{\Omega }{\int }\left|\psi \left(u\right)-\psi \left(\hat{u}\right)\right|\xi \text{d}y}\\ +{\displaystyle \underset{\Omega }{\int }sign}\left(u-\hat{u}\right)\left(f\left(u\right)-f\left(\hat{u}\right)-\varphi {\left(u\right)}_y+\varphi {\left(\hat{u}\right)}_y\right)\cdot {\xi }_y\text{d}y\\ \le {\displaystyle \underset{\Omega }{\int }sign}\left(u-\hat{u}\right)\left(s-\hat{s}\right)\xi \text{d}y+{\displaystyle \underset{\left[u=\hat{u}\right]}{\int }\text{}}\left|s-\hat{s}\right|\xi \text{d}y={[}_{\text{}}.\end{array}$ (3.20)

Proof. (Sketched) The proof of Proposition 3.13 is actually contained in the proof of Theorem 3.17 below, due to Remark 2.2. Actually a simpler argument applies, because both $f\left(\hat{u}\right)-\varphi {\left(\hat{u}\right)}_y$ and $f\left(u\right)-\varphi {\left(u\right)}_y$ have strong trace in the context of the stationary problem (S). ▄

Somewhat abusively, we will write $L^1\left(\Omega ;\left[0,1\right]\right)$ for the set of all measurable functions from [a, b] to [0, 1].

Proposition 3.14. The following properties hold true.

1) $A_{f,\varphi ,\psi }^b$ is accretive in $L^1\left(\Omega \right)$ .

2) For all $\lambda$ sufficiently small, $R\left(I+\lambda A_{f,\varphi ,\psi }^b\right)$ contains $L^1\left(\Omega ;\left[0,1\right]\right)$ .

3) $\overline{D\left(A_{f,\varphi ,\psi }^b\right)}=L^1\left(\Omega ;\left[0,1\right]\right)$ .

Proof.

1) Let $\left(u,z\right)\in A_{f,\varphi ,\psi }^b$ , $\left(\hat{u},\hat{z}\right)\in A_{f,\varphi ,\psi }^b$ . Applying Proposition 3.13 with $\xi =1$ in (3.20) and the standard properties of the bracket (see [17] ), we get

$\begin{array}{l}{\Vert b\left(u\right)-b\left(\hat{u}\right)\Vert }_{L^1\left(\Omega \right)}+{\Vert \psi \left(u\right)-\psi \left(\hat{u}\right)\Vert }_{L^1\left(\Omega \right)}\le {\left[u-\hat{u},s-\hat{s}\right]}_{L^1\left(\Omega \right)}\\ \le {\left[u-\hat{u},b\left(u\right)-b\left(\hat{u}\right)+\psi \left(u\right)-\psi \left(\hat{u}\right)+z-\hat{z}\right]}_{L^1\left(\Omega \right)}\\ \le {\Vert b\left(u\right)-b\left(\hat{u}\right)\Vert }_{L^1\left(\Omega \right)}+{\Vert \psi \left(u\right)-\psi \left(\hat{u}\right)\Vert }_{L^1\left(\Omega \right)}+{\left[u-\hat{u},z-\hat{z}\right]}_{L^1\left(\Omega \right)}.\end{array}$

We deduce that ${\left[u-\hat{u},z-\hat{z}\right]}_{L^1\left(\Omega \right)}\ge 0$ , so that $A_{f,\varphi ,\psi }^b$ is accretive.

2) For $\lambda >0$ , consider the problem

$\{\begin{array}{l}b\left(u_{\lambda }\right)+\lambda {\left(f\left(u_{\lambda }\right)-{\left(\varphi \left(u_{\lambda }\right)\right)}_y\right)}_y+\psi \left(u_{\lambda }\right)=s\text{in}\Omega ,\\ \lambda \left(f\left(u_{\lambda }\right)-\varphi {\left(u_{\lambda }\right)}_y\right)\cdot \eta \left(y\right)=0\text{on}\partial \Omega \end{array}(\; S\; \lambda \; )$

Notice that the notion of solution for $\left(\text{S}{}_{\lambda }\right)$ is like the Definition 3.7. Let $s\in L^1\left(\Omega ;\left[0,1\right]\right)$ and $\lambda >0$ then, there exists $u_{\lambda }$ entropy solution of $\left(\text{S}{}_{\lambda }\right)$ (see Proposition 3.11) such that $\left(u_{\lambda },\frac{s-b\left(u_{\lambda }\right)-\psi \left(u_{\lambda }\right)}{\lambda }\right)\in A_{f,\varphi ,\psi }^b$ .

Hence $s\in R\left(I+\lambda A_{f,\varphi ,\psi }^b\right)$ and therefore $R\left(I+\lambda A_{f,\varphi ,\psi }^b\right)\supset L^1\left(\Omega ;\left[0,1\right]\right)$ which was to be shown.

3) Let $PC\left(\left[a,b\right];\left[0,1\right]\right)$ be the set of piecewise constant functions from [a, b] to [0, 1]. Then $PC\left(\left[a,b\right];\left[0,1\right]\right)$ is dense in $L^1\left(\left[a,b\right];\left[0,1\right]\right)$ . Take $s\in PC\left(\left[a,b\right];\left[0,1\right]\right)$ , $s={\displaystyle \underset{i}{\sum }{s}_i}{1}_{\left(a_i,b_i\right)}$ where the $\left(a_i,b_i\right)$ are disjoint intervals. There exists $u_n\in L^{\infty }\left(a,b\right)$ entropy solution of $\left(S_{1/n}\right)$ , i.e., we have $\left(u_n,n\left(s-b\left(u_n\right)-\psi \left(u_n\right)\right)\right)\in A_{f,\varphi ,\psi }^b$ . For $k\in \left[0,u_c\right]\cup \left[u_s,1\right]$ , for all $\xi \in {\mathcal{C}}_0^{\infty }\left(ℝ\right)$ . We get

$\begin{array}{l}\frac{1}{n}{\displaystyle \underset{a}{\overset{b}{\int }}sign}\left(u_n-k\right)\left(f\left(u_n\right)-f\left(k\right)-{\partial }_y\varphi \left(u_n\right)\right)\cdot {\partial }_y\xi \text{d}y\\ +{\displaystyle \underset{a}{\overset{b}{\int }}sign}\left(u_n-k\right)\left(s-b\left(u_n\right)-\psi \left(u_n\right)\right)\xi \text{d}y+\frac{1}{n}{\displaystyle \underset{\left[a,b\right]}{\int }\left|f\left(k\right)\cdot \eta \left(y\right)\right|\xi \text{d}\sigma }\ge 0\end{array}$ (3.21)

For every i, one can construct ${\xi }_i^n$ such that ${\xi }_i^n\to 1_{\left(a_i,b_i\right)}$ , as $n\to \infty$ , supp ${\xi }_i^n\subset \left(a_i,b_i\right)$ , ${\Vert {\partial }_y{\xi }_i^n\Vert }_{L^{\infty }}\le 2\sqrt[3]{n}$ and ${\xi }_i^n\equiv 1$ in $\left(a_i+{\delta }_n^i,b_i-{\delta }_n^i\right)$ with ${\delta }_n^i=\frac{b_i-a_i}{2\sqrt[3]{n}}$ . Take $k=c_i$ and ${\xi }_i^n$ in (3.21).

$\begin{array}{l}{\displaystyle {\int }_{a_i+{\delta }_n^i}^{b_i-{\delta }_n^i}\left|b\left(u_n\right)+\psi \left(u_n\right)-c_i\right|\text{d}y}\\ \le \frac{1}{n}{\displaystyle {\int }_{a_i+{\delta }_n^i}^{b_i-{\delta }_n^i}sign}\left(u_n-c_i\right)\left(f\left(u_n\right)-f\left(c_i\right)-{\partial }_y\varphi \left(u_n\right)\right){\partial }_y{\xi }_i^n\text{d}y\\ \le \frac{2}{n}\left|b-a\right|{\Vert f\Vert }_{L^{\infty }}\cdot {\Vert {\partial }_y{\xi }_i^n\Vert }_{L^1}+|{\Vert \frac{1}{\sqrt{n}}{\partial }_y\varphi \left(u_n\right)\Vert }_{L^2}\cdot {\Vert {\partial }_y{\xi }_i^n\Vert }_{L^2}.\end{array}$

Then, for all $\delta >{\delta }_n^i$ , $u_n\to g$ a.e. on ${\displaystyle {\cup }_i\left(a_i+\delta ,b_i-\delta \right)}$ . We conclude by the Lebesgue theorem that $u_n\to g$ in $L^1\left(\left[a,b\right]\right)$ . In conclusion, $D\left(A_{f,\varphi ,\psi }^b\right)$ is dense in $PC\left(\left[a,b\right];\left[0,1\right]\right)$ and therefore, it is also dense in $L^1\left(\Omega ;\left[0,1\right]\right)$ . ▄

3.2.4. Integral Solution and Uniqueness

Now, we can exploit the notion of integral solution (see, e.g., [7] [17] ).

Definition 3.15. Suppose that $u_0\in L^1\left(\Omega \right)$ function $b\left(v\right)\in \mathcal{C}\left(\left[0,T\right];L^1\left(\left[a,b\right];\left[\underset{\_}{\alpha },\bar{\alpha }\right]\right)\right)$ is an integral solution of (E) if

$b\left(v\right)\left(0,.\right)=b_0(.)\text{andforall}\left(u,z\right)\in A_{f,\varphi ,\psi }^b$

${\Vert b\left(v\right)\left(t\right)-b\left(\hat{v}\right)\left(t\right)\Vert }_{L^1}+{\Vert \psi \left(v\right)\left(t\right)-\psi \left(u\right)\Vert }_{L^1\left(\Omega \right)}\le {\Vert b\left(u_0\right)-b\left({\hat{u}}_0\right)\Vert }_{L^1}.$

In particular, the integral solution is unique.

Theorem 3.17. Let $\Omega =\left[a,b\right]$ . Let v be an entropy solution of (P) and u be an entropy solution of (S). Then

$\begin{array}{l}\frac{\text{d}}{\text{d}t}{\Vert b\left(v\right)\left(t\right)-b\left(u\right)\Vert }_{\text{}{L}^1\left(\Omega \right)}+{\Vert \psi \left(v\right)-\psi \left(u\right)\Vert }_{\text{}{L}^1\left(\Omega \right)}\\ \le {\displaystyle \underset{\Omega }{\int }sign}\left(v-u\right)\left(b\left(u\right)-s\right)\text{d}x\text{in}{\mathcal{D}}^{\prime }\left(0,T\right).\end{array}$ (3.22)

In particular, $b\left(v\right)$ is an integral solution of (E).

Proof of Theorem 3.17 We adopt the doubling of variables of Kruzkhov [1] in the sense of [3] [18] . We compare regular solution and entropy solution. Keep in mind that by the result of [19] an entropy solution v of (S) is automatically time-continuous with values in $L^1\left(\Omega ;\left[0,1\right]\right)$ . We consider $v=v\left(t,x\right)$ an entropy solution of (P) and $u=u\left(y\right)$ an entropy solution of (S). Consider nonnegative function $\xi =\xi \left(t,x,y\right)$ having the property that $\xi \left(.,.,y\right)\in {\mathcal{C}}^{\infty }\left(\left[0,T\right)\times \bar{\Omega }\right)$ for each $y\in \bar{\Omega }$ , $\xi \left(t,x,.\right)\in {\mathcal{C}}_0^{\infty }\left(\bar{\Omega }\right)$ for each $\left(t,x\right)\in \left[0,T\right)\times \bar{\Omega }$ . As in [3] , we denote ${\Omega }_x=\left\{x\in \Omega ;v\left(t,x\right)\in \left[0,u_c\right]\right\}$ , ${\Omega }_y=\left\{y\in \Omega ;u\left(y\right)\in \left[0,u_c\right]\right\}$ and ${\Omega }_x^c$ , ${\Omega }_y^c$ their complementaries in $\Omega$ . To simplify the notations, take $\mathcal{F}\left[v\right]=f\left(v\right)-\varphi {\left(v\right)}_x$ , $\mathcal{F}\left[u\right]=f\left(u\right)-\varphi {\left(u\right)}_y$ and $\Phi \left(w\left(r\right),\hat{w}\left(s\right)\right)=\varphi {\left(w\right)}_r\left(\varphi {\left(w\right)}_r-\varphi {\left(\hat{w}\right)}_s\right)$ .

In (2.6), take $\xi =\xi \left(t,x,y\right)$ , $k=u\left(y\right)$ , $D=\varphi {\left(u\right)}_y$ and integrate over ${\Omega }_y^c$ . We get

$\begin{array}{l}{\displaystyle \underset{{\Omega }_y^c}{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\left\{\left|b\left(v\right)-b\left(u\right)\right|{\xi }_t+sign\left(v-u\right)\left(\mathcal{F}\left[v\right]-\mathcal{F}\left[u\right]\right)\cdot {\xi }_x\right\}\text{d}x\text{d}t\text{d}y}}}\\ +{\displaystyle \underset{{\Omega }_y^c}{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{x\in \partial \Omega }{\int }\left|\mathcal{F}\left[u\right]\cdot \eta \left(x\right)\right|\xi \text{d}\sigma \text{d}t\text{d}y}}}+{\displaystyle \underset{{\Omega }_y^c}{\int }{\displaystyle \underset{\Omega }{\int }\left|b_0-b\left(u\right)\right|\xi \left(0,x,y\right)}}\\ -{\displaystyle \underset{{\Omega }_y^c}{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\psi \left(v\right)\xi \text{d}x\text{d}t\text{d}y\ge \underset{\sigma \to 0}{\bar{\lim }}\frac{1}{\sigma }{\displaystyle \underset{{\Omega }_y^c}{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{x\in \Omega \cap A_{\varphi }^{\sigma }}{\int }\Phi }\left(v,u\right)\xi }}\end{array}$ (3.23)

In the same way, in (2.2) take $\xi =\xi \left(t,x,y\right)$ , $k=u\left(y\right)$ , integrate over ${\Omega }_y$ , and use the fact that $\varphi {\left(u\right)}_y=0$ in ${\Omega }_y$ . We get

$\begin{array}{l}{\displaystyle \underset{{\Omega }_y}{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\left\{\left|b\left(v\right)-b\left(u\right)\right|{\xi }_t+sign\left(v-u\right)\left(\mathcal{F}\left[v\right]-\mathcal{F}\left[u\right]\right)\cdot {\xi }_x\right\}}}}\\ +{\displaystyle \underset{{\Omega }_y}{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{x\in \partial \Omega }{\int }\left|\mathcal{F}\left[u\right]\cdot \eta \left(x\right)\right|\xi }}}+{\displaystyle \underset{{\Omega }_y}{\int }\text{}}{\displaystyle \underset{\Omega }{\int }\left|b_0-b\left(u\right)\right|\xi \left(0,x,y\right)}\\ -{\displaystyle \underset{{\Omega }_y}{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\psi \left(v\right)\xi \text{d}x\text{d}t\text{d}y\ge 0.\end{array}$ (3.24)

Since $\Omega ={\Omega }_x\cup {\Omega }_x^c$ , by adding (3.23) to (3.24) we obtain

$\begin{array}{l}{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\left\{\left|b\left(v\right)-b\left(u\right)\right|{\xi }_t+sign\left(v-u\right)\left(\mathcal{F}\left[v\right]-\mathcal{F}\left[u\right]\right)\cdot {\xi }_x\right\}\text{d}x\text{d}t\text{d}y}}}\\ +{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{x\in \partial \Omega }{\int }\left|\mathcal{F}\left[u\right]\cdot \eta \left(x\right)\right|\xi \text{d}\sigma \text{d}t\text{d}y}}}+{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{\Omega }{\int }\left|b_0-b\left(u\right)\right|\xi \left(0,x,y\right)\text{d}x\text{d}y}}\\ -{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sign}\left(v-u\right)\psi \left(v\right)\xi }}\ge \underset{\sigma \to 0}{\bar{\lim }}\frac{1}{\sigma }{\displaystyle \underset{{\Omega }_y^c}{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{x\in \Omega \cap A_{\varphi }^{\sigma }}{\int }\Phi }\left(v,u\right)\xi }}\text{}\end{array}$ (3.25)

In (3.18), take $\xi =\xi \left(t,x,y\right)$ , $k=v\left(t,x\right)$ , $D=\varphi {\left(v\right)}_x$ and integrate over $\left(0,T\right)\times {\Omega }_x^c$

$\begin{array}{l}{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{{\Omega }_x^c}{\int }{\displaystyle \underset{\Omega }{\int }s}ign}}\left(v-u\right)\left(\mathcal{F}\left[v\right]-\mathcal{F}\left[u\right]\right)\cdot {\xi }_y\text{d}y\text{d}x\text{d}t\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{{\Omega }_x^c}{\int }{\displaystyle \underset{y\in \partial \Omega }{\int }\left|\mathcal{F}\left[v\right]\cdot \eta \left(y\right)\right|\xi \text{d}\sigma \text{d}x\text{d}t}}}+{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\psi \left(u\right)\xi \text{d}x\text{d}t\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{{\Omega }_x^c}{\int }{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\left(b\left(u\right)-s\left(y\right)\right)\xi \ge \underset{\sigma \to 0}{\bar{\lim }}\frac{1}{\sigma }{\displaystyle \underset{{\Omega }_x^c}{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{y\in \Omega \cap A_{\varphi }^{\sigma }}{\int }\Phi }}}\left(u,v\right)\xi \end{array}$ (3.26)

Since $u\left(y\right)$ is entropy solution, then take in (3.17) $\xi =\xi \left(t,x,y\right)$ , integrate over and use the fact that in $\varphi {\left(v\right)}_x=0$ in $\left(0,T\right)\times {\Omega }_x$ .

$\begin{array}{l}{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{{\Omega }_x}{\int }{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\left(\mathcal{F}\left[v\right]-\mathcal{F}\left[u\right]\right)\cdot {\xi }_y\text{d}y\text{d}x\text{d}t\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{{\Omega }_x}{\int }{\displaystyle \underset{y\in \partial \Omega }{\int }\left|\mathcal{F}\left[v\right]\cdot \eta \left(y\right)\right|\xi \text{d}\sigma \text{d}x\text{d}t}}}+{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\psi \left(u\right)\xi \text{d}y\text{d}x\text{d}t\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{{\Omega }_x}{\int }{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\left(b\left(u\right)-s\left(y\right)\right)\xi \text{d}x\text{d}t\text{d}y\ge 0.\end{array}$ (3.27)

By adding (3.26) to (3.27), we obtain

$\begin{array}{l}{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\text{}}}{\displaystyle \underset{\Omega }{\int }s}ign\left(v-u\right)\left(\mathcal{F}\left[v\right]-\mathcal{F}\left[u\right]\right)\cdot {\xi }_y\text{d}y\text{d}x\text{d}t\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\text{}}}{\displaystyle \underset{y\in \partial \Omega }{\int }\left|\mathcal{F}\left[v\right]\cdot \eta \left(y\right)\right|\xi \text{d}\sigma \text{d}x\text{d}t}+{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\psi \left(u\right)\xi \text{d}y\text{d}x\text{d}t\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\left(b\left(u\right)-s\right)\xi \ge \underset{\sigma \to 0}{\bar{\lim }}\frac{1}{\sigma }{\displaystyle \underset{{\Omega }_x^c}{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{y\in \Omega \cap A_{\varphi }^{\sigma }}{\int }\Phi }}}\left(u,v\right)\xi .\end{array}$ (3.28)

Now, sum (3.25) and (3.28) to obtain

$\begin{array}{l}{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{\Omega }{\int }\left|b\left(v\right)-b\left(u\right)\right|{\xi }_t\text{d}y\text{d}x\text{d}t}}}+{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{\Omega }{\int }\left|b_0-b\left(u\right)\right|\xi \left(0,x,y\right)\text{d}x\text{d}y}}\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\left(\mathcal{F}\left[v\right]-\mathcal{F}\left[u\right]\right)\cdot \left({\xi }_x+{\xi }_y\right)\text{d}y\text{d}x\text{d}t\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{x\in \partial \Omega }{\int }\text{}}}{\displaystyle \underset{\Omega }{\int }\left|\mathcal{F}\left[u\right]\cdot \eta \left(x\right)\right|\xi \text{d}\sigma \text{d}t\text{d}y}+{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\text{}}{\displaystyle \underset{y\in \partial \Omega }{\int }\left|\mathcal{F}\left[v\right]\cdot \eta \left(y\right)\right|\xi \text{d}y\text{d}\sigma \text{d}t}}\\ +{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\left(b\left(u\right)-s\left(y\right)\right)\xi -{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }\left|\psi \left(v\right)-\psi \left(u\right)\right|\xi \text{d}x\text{d}t}}}\\ \ge \underset{\sigma \to 0}{\bar{\lim }}\frac{1}{\sigma }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \int {\displaystyle \underset{{\Omega }_x^c\times {\Omega }_y^c\cap A_{\varphi }^{\sigma }}{\int }{\left|\varphi {\left(v\right)}_x-\varphi {\left(u\right)}_y\right|}^2\xi \text{d}y\text{d}x\text{d}t}}}\ge 0.\end{array}$ (3.29)

Next, following the idea of [3] we consider the test function $\xi \left(t,x,y\right)=\theta \left(t\right){\rho }_n\left(x,y\right)$ , where $\theta \in {\mathcal{C}}_0^{\infty }\left(0,T\right)$ , $\theta \ge 0$ , ${\rho }_n\left(x,y\right)={\delta }_n\left(\Delta \right)$ and $\Delta =\left(1-\frac{1}{n\left(b-a\right)}\right)x-y+\frac{a+b}{2n\left(b-a\right)}$ . Then, ${\rho }_n\in \mathcal{D}\left(\bar{\Omega }\times \bar{\Omega }\right)$ and ${\rho }_{n_{{|}_{\Omega \times \partial \Omega }}}\left(x,y\right)=0$ .

Due to this choice ${\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{x\in \Omega }{\int }{\displaystyle \underset{y\in \partial \Omega }{\int }\left|\mathcal{F}\left[v\right]\cdot \eta \left(y\right)\right|{\rho }_n\theta \text{d}y\text{d}\sigma \text{d}t}}}=0$ .

By Proposition 3.11, $\left(f\left(u\right)-\varphi {\left(u\right)}_y\right)\in {\mathcal{C}}_0\left(\left[a,b\right]\right)$ . Therefore we have $\left|\mathcal{F}\left[u\right]\cdot \eta \left(x\right)\right|\to 0$ when $x\to y$ , i.e., as $n\to \infty$ . We conclude that $\underset{n\to \infty }{\lim }{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{x\in \partial \Omega }{\int }{\displaystyle \underset{y\in \Omega }{\int }\left|\mathcal{F}\left[u\right]\cdot \eta \left(x\right)\right|{\rho }_n\theta \text{d}y\text{d}\sigma \text{d}t}}}=0.$

It remains to study the limit, as $n\to \infty$

$I_n={\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }{\displaystyle \underset{\Omega }{\int }\theta sign}}}\left(v-u\right)\left(\mathcal{F}\left[v\right]-\mathcal{F}\left[u\right]\right)\cdot \left({\left({\rho }_n\right)}_x+{\left({\rho }_n\right)}_y\right)\text{d}y\text{d}x\text{d}t$

We use the change of variable $\left(x,y\right)\mapsto \left(x,z\right)$ with $z=n\left(x-y\right)-\frac{1}{b-a}x+\frac{a+b}{b-a}$ ,

$\begin{array}{c}{I}_n=\frac{2}{b-a}{\displaystyle \underset{-1}{\overset{1}{\int }}{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{\Omega }{\int }sign}}}\left(v-u\right)\left(\mathcal{F}\left[v\right]-\mathcal{F}\left[u\right]\right)\cdot {{\delta }^{\prime }}_n\left(z\right)\theta \\ =\frac{2}{b-a}{\displaystyle \underset{-1}{\overset{1}{\int }}{\displaystyle \underset{0}{\overset{T}{\int }}{\displaystyle \underset{a}{\overset{b}{\int }}sign}}}\left(v\left(t,x\right)-u_n\left(x,z\right)\right)\left(p\left(t,x\right)-q_n\left(x,z\right)\right){{\delta }^{\prime }}_n\left(z\right)\theta \left(t\right),\end{array}$

where $u_n\left(x,z\right):=u\left(y\right)$ , $p\left(t,x\right):=\mathcal{F}\left[v\right]$ and $q_n:=\mathcal{F}\left[u\right]$ .

For z given, $u_n\left(.,z\right)$ converges to $u(.)$ in $L^1$ and $q_n\left(.,z\right)$ converges to $q(.):=f\left(u\right)-\varphi {\left(u\right)}_x$ in $L^1$ . From Lemma 4.14 of [3] , we deduce that for all $z\in \left[-1,1\right]$

$K_n\left(z\right):={\displaystyle \underset{Q}{\int }sign}\left(v_n\left(t,x,z\right)\right)h_n\left(t,x,z\right)\text{d}x\text{d}t{\to }_{n\to \infty }{\displaystyle \underset{Q}{\int }sign}\left(v\right)h\text{d}x\text{d}t=:K,$

where $v_n:=v-u_n$ , $h_n:=p-q_n$ and $h:=p-q$ .

Then $K_n(.)$ converges to $K=\text{const}$ independently of z. Moreover, from the definition of $K_n$ one finds easily the uniform $L^{\infty }$ bound $\left|K_n\right|\le 2\left({\Vert p\Vert }_{L^1\left(Q\right)}+T{\Vert q\Vert }_{L^1\left(\Omega \right)}\right)$ , for n large enough. Hence by the Lebesgue theorem,

$\underset{n\to \infty }{\lim }{\displaystyle {\int }_{-1}^1{K}_n}\left(z\right){\delta }^{\prime }\left(z\right)=K{\displaystyle {\int }_{-1}^1{\delta }^{\prime }}\left(z\right)=0.$

We have shown that the limit of $I_n$ equals zero. The passage to the limit in other terms in (3.29) is straightforward. Finally (3.29) gives for $n\to \infty$

$\begin{array}{l}{\displaystyle {\int }_0^T\text{}}{\displaystyle {\int }_{\Omega }\left|b\left(v\left(t,x\right)\right)-b\left(u\left(y\right)\right)\right|{\theta }^{\prime }\left(t\right)\text{d}x\text{d}t}-{\displaystyle {\int }_0^T\text{}}{\displaystyle {\int }_{\Omega }\left|\psi \left(v\right)-\psi \left(u\right)\right|\theta }\\ +{\displaystyle {\int }_0^T\text{}}{\displaystyle {\int }_{\Omega }s}ign\left(v-u\right)\left(b\left(u\right)-s\right)\theta \ge 0.\end{array}$

Hence

$\begin{array}{l}\frac{\text{d}}{\text{d}t}{\Vert b\left(v\right)\left(t\right)-b\left(u\right)\Vert }_{L^1\left(\Omega \right)}\\ \le {\displaystyle \underset{\Omega }{\int }\left[sign\left(v-u\right)\left(b\left(u\right)-s\right)-\left|\psi \left(v\right)-\psi \left(u\right)\right|\right]\text{d}x}\text{in}{\mathcal{D}}^{\prime }\left(0,T\right).\end{array}$

Thus, $b\left(v\right)$ is an integral solution of (E).

Now, the claim of Theorem 3.12 is a direct consequence of the fact that if u is the entropy solution then $b\left(u\right)$ is an integral solution, and of Corollary 3.16. ▄

Acknowledgements

The author is grateful to Boris Andreianov for her attention to this work and for very stimulating discussions. He also thanks Stanislas Ouaro for discussions and remarks.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1]	Kruzhkov, S.N. (1970) First Order Quasilinear Equations with Several Independent Variables. Matematicheskii Sbornik, 81, 228-255.
[2]	Carrillo, J. (1999) Entropy Solutions for Nonlinear Degenerate Problems. Archive for Rational Mechanics and Analysis, 147, 269-361. https://doi.org/10.1007/s002050050152
[3]	Andreianov, B. and Karimou Gazibo, M. (2013) Entropy Formulation of Degenerate Parabolic Equation with Zero-Flux Boundary Condition. Zeitschrift für Angewandte Mathematik und Physik, 64, 1471-1491. https://doi.org/10.1007/s00033-012-0297-6
[4]	Bürger, R., Frid, H., and Karlsen, K.H. (2007) On the Well-Posedness of Entropy Solution to Conservation Laws with a Zero Fux Boundary Condition. Journal of Mathematical Analysis and Applications, 326, 108-120. https://doi.org/10.1016/j.jmaa.2006.02.072
[5]	Ouaro, S. and Touré, H. (2007) Uniqueness of Entropy Solution to Nonlinear Elliptic-Parabolic Problems. Electronic Journal of Differential Equations, No. 82, 1-15.
[6]	Ouaro, S. (2007) Entropy Solutions of Nonlinear Elliptic-Parabolic-Hyperbolic Degenerate Problems in One Dimension. International Journal of Evolution Equations, 3, 1-17.
[7]	Bénilan, Ph. and Touré, H. (1995) Sur l’équation générale ut = a(., u, φ(., u)x)x + v dans L1. I. Etude du problème stationnaire. In: Evolution Equations (Baton Rouge, LA, 1992), Dekker, New York, 35-62.
[8]	Andreianov, B. and Wittbold, P. (2012) Convergence of Approximate Solutions to an Elliptic-Parabolic Equation without Structure Condition. Nonlinear Differential Equations and Applications, 19, 695-717. https://doi.org/10.1007/s00030-011-0148-8
[9]	Ammar, K. and Wittbold, P. (2003) Existence of Renormalized Solutions of Degenerate Elliptic Parabolic Problems. Proceedings of the Royal Society of Edinburgh Section A, 133, 477-496. https://doi.org/10.1017/S0308210500002493
[10]	Andreianov, B., Bendahmane, M., Karlsen, K.H. and Ouaro, S. (2009) Well-Posedness Results for Triply Nonlinear Degenerate Parabolic Equations. Journal of Differential Equations, 247, 277-302. https://doi.org/10.1016/j.jde.2009.03.001
[11]	Andreianov, B. and Gazibo, M.K. (2016) Explicit Formulation for the Dirichlet Problem for Parabolic-Hyperbolic Conservation Laws. Networks and Heterogeneous Media, 11, 203-222. https://doi.org/10.3934/nhm.2016.11.203
[12]	Alt, H.W. and Luckhaus, S. (1981) Quasilinear Elliptic-Parabolic Differential Equations. Mathematische Zeitschrift, 183, 311-341. https://doi.org/10.1007/BF01176474
[13]	Ammar, K., Carrillo, J. and Willbold, P. (2006) Scalar Conservation Laws with General Boundary Condition and Continuous Flux Function. Journal of Differential Equations, 228, 111-139. https://doi.org/10.1016/j.jde.2006.05.002
[14]	Ammar, K. and Redwane, H. (2008) Degenerate Stationary Problems with Homogeneous Boundary Conditions. Electronic Journal of Differential Equations, 30, 1-18.
[15]	Vasseur, A. (2001) Strong Traces for Solutions of Multidimensional Scalar Conservation Laws. Archive for Rational Mechanics and Analysis, 160, 181-193. https://doi.org/10.1007/s002050100157
[16]	Gazibo Karimou, M. (2023) Degenerate Parabolic Equation with Zero Flux Boundary Condition and Its Approximations. Preprint.
[17]	Bénilan, Ph., Crandall, M.G. and Pazy, A. (1994) Nonlinear Evolution Equations in Banach Spaces. Preprint Book.
[18]	Andreianov, B. and Bouhsiss, F. (2004) Uniqueness for an Elliptic-Parabolic Problem with Neumann Boundary Condition. Journal of Evolution Equations, 4, 273-295. https://doi.org/10.1007/978-3-0348-7924-8_37
[19]	Cances, C. and Gallouet, Th. (2011) On the Time Continuity of Entropy Solutions. Journal of Evolution Equations, 11, 43-45. https://doi.org/10.1007/s00028-010-0080-0