Open Access
Int. J. Simul. Multidisci. Des. Optim.
Volume 8, 2017
Article Number A11
Number of page(s) 8
Published online 18 August 2017

© H. El Ouardi et al., published by EDP Sciences, 2017

Licence Creative Commons
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Let (N ≥ 1) be a bounded Lipshitz domain and 0 < T < ∞. It will be assumed throughout this paper that p(x) is continuous function defined in with logarithmic module of continuity: (1)

We set Q T  = Ω × (0, T) and Σ T  = ∂Ω × (0, T). Our aim is to prove the existence and uniqueness of solutions u = (u 1, u 2) to the nonlinear (p 1(x), p 2(x))-Laplacian system:(2)where is a function, (i = 1, 2).

The operator is called p(x)-Laplacian, which will be reduced to the p-Laplacian when p(x) = p a constant.

The (p 1(x), p 2(x))-Laplacian system (2) can be viewed as a generalization of (p, q)-Laplacian system(3)

For the case p i (x) = p i  > 2, and a i (x) = 1, (i = 1, 2), system (2) models as non-Newtonian fluids [2, 21] and nonlinear filtration [2], etc. In the non-Newtonian fluids theory, (p 1, p 2) is a characteristic quantity of the fluids, there have been many results about the existence, uniqueness of the solutions. We refer the readers to the bibliography given in [7, 9, 10, 11, 28, 31] and the references therein.

In recent years, the research of nonlinear problems with variable exponent growth conditions has been an interesting topic. p(·)-growth problems can be regarded as a kind of nonstandard growth problems and these problems possess very complicated nonlinearities, for instance, the p(x)-Laplacian operator is inhomogeneous. And these problems have many important applications in nonlinear elastic, electrorheological fluids and image restoration. The reader can find in [14, 22] several models in mathematical physics where this class of problem appears.

The case of a single equation of the type (2) has been studied in [46, 20] and the authors established the existence and uniqueness results, in [20], the authors use the difference scheme to transform the parabolic problem to a sequence of elliptic problems and then obtain the existence of solutions with less constraint to p i (x).

The more interesting question concerning parabolic systems of (p 1(x), p 2(x))-Laplacian type is to understand the asymptotic behavior of solutions when time goes to infinity. The study of the asymptotic behavior of the system is giving us relevant information about the structure of the phenomenon described in the model.

Concerning the elliptic systems with variable exponents, the results about existence and non-existence are proved in [26, 27, 29, 30].

Note that system (2) has a more complicated nonlinearity than the classical (p,q)-Laplacian system since it is nonhomogenous.

Recently, [24] study the equation the p(x)-Laplacian equationwhere λ > 0, m + λ − 2 > 0 and a(x) is a positive continuous function. They examine under which conditions on behavior of a(x) corresponding nonnegative solutions of the Cauchy problems possess the finite speed of propagations or interface blow-up phenomena.

In this paper, we consider the existence and uniqueness for the problem of the type (2) under some assumptions. The proof consists of two steps. First, we prove that the approximating problem admits a global solution; then we do some uniform estimates for these solutions. We mainly use skills of inequality estimation and the method of approximation solutions. By a standard limiting process, we obtain the existence to problem of the type (2).

The outline of this paper is the following: In Section 2, we introduce some basic Lebesgue and Sobolev spaces and state our main theorems. In Section 3, we give the existence and uniqueness of weak solutions. In Section 4, the blow-up results will be proved. The asymptotic behaviour of solution is established in Section 5.

2 Preliminaries

To consider problems with variable exponents, one needs the basic theory of spaces L p(x)(Ω) and W 1p(x)(Ω). For the convenience of readers, let us review them briefly here. The details and more properties of variable-exponent Lebesgue-Sobolev spaces can be found in [16, 17].

Let When p  > 1, one can introduce the variable-exponent Lebesgue spaceendowed with the Luxemburg norm.

The conjugate space is L q(x)(Ω), with

As in the case of a constant exponent, setendowed with the norm

Similarly we also denote by the closure of in and is the dual of with respect to the inner product in

In Propositions 2.1–2.3, we describe some results about the Luxembourg norm.

Proposition 2.1 [16, 17]

  1. The space is a separable, uniformly convex Banach space, and its conjugate space is where For any and we have the following Hölder-type inequality:

  2. If for any the imbedding is continuous, the norm of the imbedding does not exceed

Proposition 2.2 [16]

If we denotethen

  1. |w|(r(x)) < 1(= 1; > 1) ⇔ ρ(w) < 1 (= 1; > 1);

  2. |w|(r(x)) → 0 ⇔ ρ(w) → 0; |w|(r(x))ρ(w) → ∞.

Proposition 2.3 [16]

For there exists a constant such that:

This implies that and are equivalent norms of

System (2) does not admit classical solutions in general. So, we introduce weak solutions in the following sence.

Definition 2.4 A function u = (u 1 , u 2) is said to be a weak solution of equation (2) , if the following conditions are satisfied:

  1. such that:

  2. For any

In the study of the global existence of solutions, we need the following hypotheses (H):

  • (H1)

  • (H2)

  • (H3)

3 Main results

Remark 3.1 In this paper, we shall denote by c i , C i differents constants, depending on p i (x), T, Ω, but not on n, which may vary from line to line. Sometimes we shall refer to a constant depending on specific parameters C i (T), etc.

Our main existence result is the following:

Theorem 3.2 Let (H1)–(H3) hold. Then system (2) admits a unique solution u = (u 1 ,u 2  (C((0, T); L 2(Ω)))2 . Moreover, the mapping (φ 1, φ 2) → (u 1(t), u 2(t)) is continuous in L 2(Ω) × L 2(Ω).

Proof of the main results.

3.1 Existence

We will semi-discrete (2) in time and solve the corresponding elliptic problem. Based on the semi-discrete problem, we construct the corresponding approximate solutions. The key procedure is to establish necessary a priori estimates for finding the limit of the approximate solutions via a compactness argument.

We first consider the discrete scheme (4)(4)where  = T and T is a fixed positive real, and 1 ≤ n ≤ N.

Lemma 3.3 For any fixed n, if Problem (4) admits a weak solution

Proof. On the space we consider the functionalwhere is a known function. Using Young’s inequality and Proposition 2.1, there exist constants C 1, C 2 > 0, such that:hence Φ(v) → ∞, as Since the norm is lower semi-continuous and is continuous functional, Φ(v) is weakly lower semi-continuous on and satisfy the coercive condition. From [14] we conclude that there exists such that:and v* is the weak solutions of the Euler equation corresponding to Φ(v),

Choosing we obtain a weak solution of (4).(5)

Since , we may prove by induction that (4) has a solution in We put and for any integer k > 0, we may take as a test function in (5) to get

By using the Hölder inequality and , we have

We deduce

Letting k, we get Consider , we get easily that , i.e. and if we choose such that , we obtain

This completes the proof of lemma 3.3.

Now, we define the approximate solutions as set by: for all n ∈ {1, …, N}.are well defined and satisfied in addition(6)

We first establish some energy estimates of .

We need several lemmas to complete the proof of Theorem 3.2.

Lemma 3.4 There exists a positive constant C(T, u 0) such that, for all n = 1, …, N (7) (8) (9)and(10)

Proof. (a) By lemma 3.3, for any n ∈ N, is bounded; whence (7)

(b) Multiplying (4) by , summing from n = 1 to N and integrating over Ω, we obtain(11)

By using the Young inequality, for small, there exists such that(12)

With the aid of the identity 2α(α − β) = α 2 − β 2 + (α − β)2, we get

With the above estimates, we get (8).

(c) Multiplying the equation (4) by and summing from n = 1 to N, we get

By using the Young inequality, we get(13)

From the convexity of the expression we get the following inequality:(14)

which imply with (12) and (13) that

By lemma 3.4, there exists M i  > 0 independent of such that:(15)

Therefore, taking and up to subsequence, we get that there exists such that , and as ,(16) (17)

From (14), it follows that u i  = v i . From (15), we get that(18)

By Aubin-Simon’s compactness results [22], we have(19)

Now, multiplying (4) by and using (14) and (15), we get by straightforward calculations:where as .

Thus, we get that(20)and from (16) we have thus,and consequently by the same as that in [24]

Therefore, u i satisfies (3).

3.2 Uniqueness

Let (H1)–(H3) be satisfied. Then system (2) has a unique solution u = (u 1, u 2) in Q T.

Proof. Let u = (u 1, u 2) and v = (v 1, v 2) be solutions of (2), we have:

Since f i (x,.,.) is locally Lipschitz uniformly in Ω, the difference w i  = u i  − v i satisfieswe observe that w → −Δ p(x) w is monotone from to (21)

We finally deduce from Gronwall’s lemma,

Thus, we deduce that u i  = v i .

Thus the solution is unique. The continuity of the the mapping can be obtained similarly.

Remark 3.5 From Theorem 3.2, the solution of system (2) generates a semigroup in L 2(Ω) × L 2(Ω).

Remark 3.6 If we assume that f i (x, s 1, s 2) = g i (x, s 1, s 2) − h i (x, s i ) where is Carathéodory mapping and that there exist positive constants L j , c j , m j , C j such that:

  1. for any x ∈ Ω and

  2. for any x ∈ Ω and where with and

where with and

By the same argument as that in [12, 32], One can show in the same way that the semigroup associated with system (2) admits an absorbing set in there is a bounded set such that, for any bounded set B in L 2(Ω) × L 2(Ω), there exists a T 0 > 0 such that for any t ≥ T 0. Where T 0 depends only on B.

4 Blow-up of solutions

In this section, we shall investigate the blow-up properties of solutions to system (2) using energy methods. To this end, we consider the following hypotheses on the data.

(H4) such that:

(H5) and H is such that:

Throughout this section, we define for t ≥ 0

Theorem 4.1 Let (H1)–(H5) be satisfied, then the solutions of system (2) blow up in finite time, namely, there exists a T* < ∞ such that as t → T*.

Proof. We define .

Multiplying the first equation of (2) by , the second by , integrating by parts, we have(22)which implies that E(t) ≤ E(0).

Next define

Multiplying the first equation of (2) by u 1, the second by u 2, integrating by parts, we have(23)

By using Hölder’s inequality, we have(24)where

By the formula we have by combining (23), (24)where

A direct integration of the above inequality over (0, t) then yieldswhich implies that g(t) bows up at a finite time and so does u i .

5 Asymptotic behaviour

This section is devoted to the asymptotic behaviour of solutions. In order to prove the asymptotic behaviour, we assume


Theorem 5.1 The weak solution u = (u 1(t), u 2(t)) obtained in Theorem 3.2, satifies: where C i  > 0 (i = 1, 2, 3), or or .

Proof. Let u i be solution of (2).

Multiplying the first equation in (2) by u 1 and integrating over Q T ,(25)

Multiplying the second equation in (2) by u 2 and integrating over Q T ,(26)

Summing up (25) and (26), we have from hypothesis (H5) that(27)

By u i using Poincaré’s inequality, we obtain(28)

If and by Proposition 2.2,(29)

According to the assumption that p 1(x) ≤ p 2(x), Then

Hence, we get from (28) that(30)

By the formula we have(31)this implies that(32)where C 3 = min(C 1, C 2).


Then, we obtain from (32) and (H2) that(33)

If and by Proposition 2.2,

Then we get from (28) that(34)

That is(35)

Again we have

Similarly, if and or and we can also obtain the similar results


The proof is complete.


  1. Astrita G, Marrucci G. 1974. Principles of non-Newtonian fluid mechanics. McGraw-Hill: New York.
  2. Acerbi E, Mingine G. 2002. Regularity results for stationary electrorheological fluids. Arch. Ration. Mech. Anal., 164, 213–259. [CrossRef] [MathSciNet] (In the text)
  3. Acerbi E, Mingione G, Seregin GA. 2004. Regularity results for parabolic systems related to a class of non Newtonian fluids. Ann. Ins. H. Poincaré Anal. Non Linéaire, 21, 25–60.
  4. Antontsev SN, Shmarev SI. 2009. Anisotropic parabolic equations with variable nonlinearity. Publ. Mat., 53(2), 355–399. [CrossRef] [MathSciNet] (In the text)
  5. Antontsev SN, Shmarev SI. 2005. A model porous medium equation with variable exponent of nonlinearity: existence, uniqueness and localization properties of solutions. Nonlinear Anal., 60, 515–545. [CrossRef] [MathSciNet]
  6. Antontsev SN, Shmarev SI. 2009. Localization of solutions of anisotropic parabolic equations. Nonlinear Anal., 71, 725–737. [CrossRef] (In the text)
  7. El Ouardi H, de Thelin F. 1989. Supersolutions and stabilization of the solutions of a nonlinear parabolic system. Publicacions Mathematiques, 33, 369–381. [CrossRef] [MathSciNet] (In the text)
  8. El Ouardi H, El Hachimi A. 2002. Existence and regularity of a global attractor for doubly nonlinear parabolic equations. Electron. J. Diff. Eqns., 2002(45), 1–15.
  9. El Ouardi H, El Hachimi A. 2001. Existence and attractors of solutions for nonlinear parabolic systems. Electron. J. Qual. Theory Differ. Equ., 5, 1–16. [CrossRef] [EDP Sciences] [MathSciNet] (In the text)
  10. El Ouardi H, El Hachimi A. 2006. Attractors for a class of doubly nonlinear parabolic systems. Electron. J. Qual. Theory Differ. Equ., 1, 1–15. [CrossRef] [EDP Sciences] [MathSciNet] (In the text)
  11. El Ouardi H. 2011. Global attractor for quasilinear parabolic systems involving weighted p-Laplacian operators. J. Pure App. Mat. Advances App., 5(2), 79–97. [MathSciNet] (In the text)
  12. El Ouardi H. 2013. Long-time behavior for a nonlinear parabolic systems with variable exponent of nonlinearity. Int. J. Adv. Math. Math. Sci., 2(1), 27–36. (In the text)
  13. Esteban JR, Vazquez JL. 1982. On the equation of turbulent filtration in one-dimensional porous media. Nonlinear Anal., 10, 1303–1325. [CrossRef]
  14. Kristaly A, Radulescu V, Varga C. 2010. Variational principles in mathematical physics, geometry, and economics. Qualitative analysis of nonlinear equations and unilateral problems in Encyclopedia of mathematics and its applications, vol. 136. Cambridge University Press: Cambridge. (In the text)
  15. Chang KC. 1986. Critical point theory and its applications. Shanghai Scientific and Technology Press: Shanghai.
  16. Fan X. 2007. Global C1,α regularity for variable exponent elliptic equations in divergence form. J. Diff. Equ., 235, 397–417. [CrossRef] (In the text)
  17. Fan X, Zhao D. 2001. On the space Lp(x)(Ω) and Wm,p(x) (Ω). J. Math. Anal. Appl., 263, 749–760. [CrossRef] (In the text)
  18. Fan X, Zhang QH. 2003. Existence of solutions for p(x)-Laplacian Dirichlet problem. Nonlinear Anal., 52, 1843–1852. [CrossRef] [MathSciNet]
  19. Lions JL. 1969. Quelques méthodes de résolution des problèmes aux limites non linéaires. Dunod: Paris.
  20. Lian S, Gao W, Yuan H, Cao C. 2012. Existence of solutions to an initial problem of evolutional p(x)-Laplace equations. Ann. I. H. Poincaré – AN, 29(3), 377–399, DOI: 10.1016/j.anihpc.2012.01.001. [CrossRef] [MathSciNet] (In the text)
  21. Martinson LK, Pavlov KB. 1971. Unsteady shear flows of a conducting fluid with a rheological power law. Magnitnaya Gidrodinamika, 2, 50–58. (In the text)
  22. Ruzicka M. 2000. Electrorheological fluids: Modeling and mathematical theory. Lecture Notes in Mathematics. vol. 1748. Berlin: Springer-Verlag. [CrossRef] (In the text)
  23. Simon J. 1987. Compact sets in Lp(0, T; B). Ann. Mat. Pura Appl., 146, 65–96. [CrossRef] [MathSciNet]
  24. Tedeev AF. 2007. The interface blow-up phenomenon and local estimates for doubly degenerate parabolic equation. Appl. Anal., 86(6), 755. [CrossRef] [MathSciNet] (In the text)
  25. Zhao JN. 1993. Existence and nonexistence of solutions for ut = div(|∇u|p − 2u) + f(∇u, u, x, t). J. Math. Anal. Appl., 172(1), 130–146. [CrossRef] [MathSciNet]
  26. Fu Y, Zhang X. 2009. Multiple solutions for a class of p(x)-Laplacian systems. J. Inequal. Appl. Art., 2009(191649), 1–12. (In the text)
  27. El Hamidi A. 2004. Existence results to elliptic systems with non standard growth condition. J. Math. Anal. Appl., 300, 30–42. [CrossRef] [MathSciNet] (In the text)
  28. Marion M. 1987. Attractors for reaction-diffusion equation: existence and estimate of their dimension. Appl. Anal., 25, 101–147. [CrossRef] [MathSciNet] (In the text)
  29. Ogras S, Mashiyev A, Avci M, Yucedag Z. 2008. Existence of solution for a class of elliptic systems in RN involving the (p(x), q(x))-Laplacian. J. Inequal. Appl., 612938, 16. [MathSciNet] (In the text)
  30. Xu X, An Y. 2008. Existence and multiplicity of solutions for elliptic systems with nonstandard growth conditions in RN. Nonlinear Anal., 68, 956–968. [CrossRef] [MathSciNet] (In the text)
  31. Wei Y, Gao W. 2007. Existence and uniqueness of local solutions to a class of quasilinear degenerate parabolic systems. Appl. Math. Comput., 190, 1250–1257. [CrossRef] [MathSciNet] (In the text)
  32. Niu W. 2012. Long-time behavior for a nonlinear parabolic problem with variable exponents. J. Math. Anal. Appl., 393, 56–65. [CrossRef] [MathSciNet] (In the text)

Cite this article as: El Ouardi H & Ghabbar Y: Study of solutions to a class of certain parabolic systems with variable exponents. Int. J. Simul. Multisci. Des. Optim., 2017, 8, A11.

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.