Degenerated and Competing Dirichlet Problems with Weights and Convection

Abstract

This paper focuses on two Dirichlet boundary value problems whose differential operators in the principal part exhibit a lack of ellipticity and contain a convection term (depending on the solution and its gradient). They are driven by a degenerated $(p,q)$-Laplacian with weights and a competing $(p,q)$-Laplacian with weights, respectively. The notion of competing $(p,q)$-Laplacians with weights is considered for the first time. We present existence and approximation results that hold under the same set of hypotheses on the convection term for both problems. The proofs are based on weighted Sobolev spaces, Nemytskij operators, a fixed point argument and finite dimensional approximation. A detailed example illustrates the effective applicability of our results.

1. Introduction

Consider a bounded domain $\Omega$ in ${\mathbb{R}}^N$ ($N\ge 1$) with a Lipschitz boundary $\partial \Omega$, numbers $1<q<p<\infty$, functions $a,b\in L^1(\Omega )$ with $a(x),b(x)>0$ for a.e. $x\in \Omega$ and a Carathéodory function $f:\Omega \times \mathbb{R}\times {\mathbb{R}}^N\to \mathbb{R}$ (i.e., $f(·,t,\xi )$ is measurable on $\Omega$ for each $(t,\xi )\in \mathbb{R}\times {\mathbb{R}}^N$ and $f(x,·,·)$ is continuous on $\mathbb{R}\times {\mathbb{R}}^N$ for a.e. $x\in \Omega$). The aim of this paper is to investigate the quasilinear Dirichlet problems

$$\left\{\begin{array}{cc}-\mathrm{div}(a(x)|\nabla u{|}^{p-2}\nabla u+b(x)|\nabla u{|}^{q-2}\nabla u)=f(x,u,\nabla u)\hfill & \mathrm{in}\Omega \hfill \\ u=0\hfill & \mathrm{on}\partial \Omega \hfill \end{array}\right.$$

(1)

and

$$\left\{\begin{array}{cc}-\mathrm{div}(a(x)|\nabla u{|}^{p-2}\nabla u-b(x)|\nabla u{|}^{q-2}\nabla u)=f(x,u,\nabla u)\hfill & \mathrm{in}\Omega \hfill \\ u=0\hfill & \mathrm{on}\partial \Omega .\hfill \end{array}\right.$$

(2)

Notice that problem (1) is driven by a sum of weighted p-Laplacians, whereas problem (2) by a difference of weighted p-Laplacians. The weights $a\in L^1(\Omega )$ and $b\in L^1(\Omega )$ are strongly related to the ellipticity property, but act in a fundamentally different way in these problems. The celebrated p-Laplacian and q-Laplacian are used instead of more general operators in the above formulations just to highlight the main ideas.

The differential operator in the principal part of Equation (1) is the sum

$$u\mapsto {\mathrm{div}(a(x)|\nabla u|}^{p-2}{\nabla u)+\mathrm{div}(b(x)|\nabla u|}^{p-2}\nabla u)$$

of the degenerated p-Laplacian with weight $a\in L^1(\Omega )$ and the degenerated q-Laplacian with weight $b\in L^1(\Omega )$ that should be consistent. This operator was introduced in [1] where it was called the degenerated $(p,q)$-Laplacian with weights $a,b\in L^1(\Omega )$. Its construction is reviewed in Section 2. The characteristic property of this operator is the degeneracy, meaning that one cannot guarantee the existence of a constant $k>0$ to have

$${\langle -\mathrm{div}(a(x)|\nabla u|}^{p-2}+{b(x)|\nabla u|}^{q-2})\nabla u),u\rangle \ge k{\int }_{\Omega }{(|\nabla u(x)|}^p+{|\nabla u(x)|}^q)dx.$$

Due to this, one cannot apply the classical elliptic theory.

The differential operator in the principal part of Equation (2) is the difference

$$u\mapsto {\mathrm{div}(a(x)|\nabla u|}^{p-2}{\nabla u)-\mathrm{div}(b(x)|\nabla u|}^{p-2}\nabla u)$$

of the degenerated p-Laplacian with weight $a\in L^1(\Omega )$ and of the degenerated q-Laplacian with weight $b\in L^1(\Omega )$. Such a nonlinear operator with weights is considered for the first time. We call it the competing $(p,q)$-Laplacian with weights $a,b\in L^1(\Omega )$. In this case, we go beyond the degeneracy, actually completely dropping the ellipticity because the quantity

$${\langle -\mathrm{div}(a(x)|\nabla u|}^{p-2}\nabla u-{b(x)|\nabla u|}^{q-2}\nabla u),u\rangle ={\int }_{\Omega }{(a(x)|\nabla u(x)|}^p-{b(x)|\nabla u(x)|}^q)dx$$

can have an arbitrary sign (note that $a(x)$ and $b(x)$ are positive). For problem (2), any method of monotone type, including the use of pseudomonotone operators, fails to apply.

The right-hand side $f(x,u,\nabla u)$ of the equations in (1) and (2) is a convection term; that is, it depends on the solution u and on its gradient $\nabla u$. The dependence on the gradient $\nabla u$ generally prevents having a variational structure for problems (1) and (2), so the variational methods are not applicable. In order to find the needed estimates, an essential part of our development is devoted to the Nemytskij operator associated with the convection term $f(x,u,\nabla u)$ under an appropriate growth condition for the function $f(x,t,\xi )$ on $\Omega \times \mathbb{R}\times {\mathbb{R}}^N$. Different results regarding unweighted problems involving $(p,q)$-Laplacian and convection terms can be found in [2].

The problems (1) and (2) have only recently been regarded in their generality. To the best of our knowledge, there is solely the existence theorem for problem (1), obtained in [1] through the theory of pseudomonotone operators. For the particular case of (1) where the equation is governed by a degenerated p-Laplacian (i.e., $b=0$ in (1)), existence results based on minimization and degree theoretic methods can be found in [3] and a method to create a sub-supersolution was developed in [4]. Concerning problem (2) driven by competing operators, there is no available result except for the most particular situation where $a(x)=b(x)\equiv 1$ in $\Omega$ (i.e., the problem without weights), whose study was initiated in [5] and continued in [6,7].

In the present paper, we overcome the lack of ellipticity, monotonicity and variational structure in problems (1) and (2) by means of a passing to limit process involving approximate solutions generated through fixed point arguments on finite dimensional spaces. This approach was implemented in [6,7] for unweighted problems (i.e., $a(x)=b(x)\equiv 1$ in $\Omega$). Here, the development is substantially modified due to the completely different functional setting under the weights $a\in L^1(\Omega )$ and $b\in L^1(\Omega )$.

For problem (1), we are able to establish the existence of a solution in a weak sense, whereas for problem (2), we prove the existence of a solution in a generalized sense. It is worth noting that in the case of problem (1) any generalized solution is a weak solution. Moreover, our results can be viewed as providing approximations in the sense of strong convergence for solutions to problems (1) and (2) by finite dimensional approximate solutions.

Inspired by [3], a major step in our treatment is a reduction within the framework of classical Sobolev spaces. We impose a suitable growth condition for the convection term $f(x,u,\nabla u)$ to match this reduction. The growth condition is expressed using a positive quantity ($p_s$ in the text) described by the weights $a\in L^1(\Omega )$ and $b\in L^1(\Omega )$, which provide the best integrability rate.

We plan to use the present work for studying evolutionary counterparts for problems (1) and (2).

The rest of the paper is organized as follows. Section 2 presents the degenerated and competing $(p,q)$-Laplacians with weights. Section 3 sets forth the associated Nemytskij operator. Section 4 contains our main result on the solvability and approximation for problem (1). Section 5 focuses on the solvability of problem (2). Section 6 illustrates by an example the effective applicability of our theorems.

2. Degenerated and Competing $(\mathit{P},\mathit{Q})$-Laplacians with Weights

Throughout the text, we denote by → the strong convergence and by ⇀ the weak convergence in any normed space X under consideration. The norm on X is denoted by ${\parallel ·\parallel }_X$, while the notation ${\langle ·,·\rangle }_X$ stands for the duality pairing between X and its dual $X^{*}$. For the rest of the paper, by a bounded map we understand a map between normed spaces that maps bounded sets to bounded sets.

We fix the framework for the underlying weighted Sobolev spaces related to problems (1) and (2). For a systematic study of weighted Sobolev spaces, we refer to [3,8]. The completeness property for such spaces is discussed in [9]. This functional setting was also discussed in [1].

Given a real number $p\in (1,+\infty )$ and a positive function $a\in L^1(\Omega )$, the weighted space

$$W^{1,p}(a,\Omega ):=\{u\in L^p(\Omega ):{\int }_{\Omega }{a(x)|\nabla u(x)|}^{p}dx<\infty \},$$

is endowed with the norm

$${\parallel u\parallel }_{W^{1,p}(a,\Omega )}:={\left({\parallel u\parallel }_{L^p(\Omega )}^p+{\int }_{\Omega }a(x){|\nabla u(x)|}^{p}dx\right)}^{\frac{1}{p}},\forall u\in W^{1,p}(a,\Omega ).$$

We note that $C_0^{\infty }(\Omega )\subset W^{1,p}(a,\Omega )$. The closure of $C_0^{\infty }(\Omega )$ in $W^{1,p}(a,\Omega )$ with respect to the norm ${\parallel ·\parallel }_{W^{1,p}(a,\Omega )}$ is the space $W_0^{1,p}(a,\Omega )$. The dual spaces of $W^{1,p}(a,\Omega )$ and $W_0^{1,p}(a,\Omega )$ are denoted by $W^{1,p}{(a,\Omega )}^{*}$ and $W_0^{1,p}{(a,\Omega )}^{*}$, respectively.

A reduction in the setting of classical Sobolev spaces is based on the following condition from [3] (p. 26):

(H1).

$a^{-s}\in L^1(\Omega )$ for some $s\in \left(max\left\{\frac{N}{p},\frac{1}{p-1}\right\},+\infty \right)$.

Proposition 1.

Under condition (H1), there are the continuous embeddings

$$\begin{array}{c}\hfill W^{1,p}(a,\Omega )\hookrightarrow W^{1,p_s}(\Omega )\hookrightarrow L^p(\Omega ),\end{array}$$

(3)

where

$$\begin{array}{c}\hfill p_s={\displaystyle \frac{ps}{s+1}}.\end{array}$$

(4)

In addition, the embedding $W^{1,p_s}(\Omega )\hookrightarrow L^p(\Omega )$ is compact. Furthermore,

$${\parallel u\parallel }_{W_0^{1,p}(a,\Omega )}:={\left({\int }_{\Omega }a(x){|\nabla u(x)|}^{p}dx\right)}^{\frac{1}{p}},\forall u\in W_0^{1,p}(a,\Omega ),$$

is an equivalent norm on $W_0^{1,p}(a,\Omega )$ for which $W_0^{1,p}(a,\Omega )$ becomes a uniformly convex Banach space.

Proof. 

The proof is essentially completed in [3]. For the sake of clarity, we highlight aspects relevant for problems (1) and (2).

It can be seen from (4) that $p_s>1$ if and only if $s>1/(p-1)$, which by assumption (H1) is true. In order to prove the first inclusion in (3), let $u\in W^{1,p}(a,\Omega )$. Using Hölder’s inequality, hypothesis (H1) and (4) (note $p_s<p$), we infer that

$$\begin{array}{ccc}& & {\int }_{\Omega }{|\nabla u(x)|}^{p_s}dx={\int }_{\Omega }{(a(x)}^{\frac{p_s}{p}}{|\nabla u(x)|}^{p_s}{)a(x)}^{-\frac{p_s}{p}}dx\hfill \\ & & \le {\left({\int }_{\Omega }a(x){|\nabla u(x)|}^{p}dx\right)}^{\frac{p_s}{p}}{\left({\int }_{\Omega }a{(x)}^{-\frac{p_s}{p-p_s}}dx\right)}^{\frac{p-p_s}{p}}\hfill \\ & & \le \parallel a^{-s}{\parallel }_{L^1(\Omega )}^{\frac{1}{s+1}}{\parallel u\parallel }_{W^{1,p}(a,\Omega )}^{p_s},\forall u\in W_0^{1,p}(a,\Omega ).\hfill \end{array}$$

The continuous inclusion $W^{1,p}(a,\Omega )\hookrightarrow W^{1,p_s}(\Omega )$ is proven.

The Rellich–Kondrachov embedding theorem ensures the compact embedding $W^{1,p_s}$$(\Omega )\hookrightarrow L^r(\Omega )$, with $1\le r<p_s^{*}$, where $p_s^{*}$ is the critical exponent corresponding to $p_s$, that is,

$$\begin{array}{c}\hfill p_s^{*}:=\left\{\begin{array}{cc}\frac{N{p}_s}{N-p_s}\hfill & \mathrm{if}N>p_s(\iff ps<N(s+1))\hfill \\ \\ +\infty \hfill & \mathrm{if}N\le p_s(\iff ps\ge N(s+1)).\hfill \end{array}\right.\end{array}$$

We have that $p_s^{*}>p$ if and only if $s>N/p$. Since the latter holds by assumption (H1), the compactness of the second inclusion in (3) follows.

The desired equivalence of norms is a consequence of (3) and the Poincaré inequality on $W_0^{1,p_s}(\Omega )$ because with a positive constant C,

$$\begin{array}{c}\hfill {\parallel u\parallel }_{L^p(\Omega )}\le C{\parallel u\parallel }_{W^{1,p}(a,\Omega )},\forall u\in W_0^{1,p}(a,\Omega ).\end{array}$$

It remains to show that $W_0^{1,p}(a,\Omega )$ is a uniformly convex Banach space. It suffices to have $a^{-\frac{1}{p-1}}\in L^1(\Omega )$ (see [3] Theorem 1.3). From hypothesis (H1), it is known that $a^{-s}\in L^1(\Omega )$ with $s>1/(p-1)$, which results in

$$\begin{array}{ccc}& & {\int }_{\Omega }a{(x)}^{-\frac{1}{p-1}}dx={\int }_{\{a(x)<1\}}a{(x)}^{-\frac{1}{p-1}}dx+{\int }_{\{a(x)\ge 1\}}a{(x)}^{-\frac{1}{p-1}}dx\hfill \\ & & \le {\int }_{\Omega }a{(x)}^{-s}dx+\mathrm{meas}(\Omega )<\infty ,\hfill \end{array}$$

thus completing the proof. □

The degenerated p-Laplacian with the weight $a\in L^1(\Omega )$ is defined as the map ${\Delta }_p^a:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ given by ${\Delta }_p^a(u)={\mathrm{div}(a(x)|\nabla u|}^{p-2}\nabla u)$ for all $u\in W_0^{1,p}(a,\Omega )$, i.e.,

$$\begin{array}{c}\hfill {\langle -{\Delta }_p^a(u),v\rangle }_{W_0^{1,p}(a,\Omega )}={\int }_{\Omega }a(x){|\nabla u(x)|}^{p-2}\nabla u(x)\nabla v(x)dx,\forall u,v\in W_0^{1,p}(a,\Omega ).\end{array}$$

The definition makes sense as can be seen through Hölder’s inequality

$$\begin{array}{ccc}& & \left|{\int }_{\Omega }a(x){|\nabla u(x)|}^{p-2}\nabla u(x)\nabla v(x)dx\right|\hfill \\ & & \le {\int }_{\Omega }{(a(x)}^{\frac{p-1}{p}}{|\nabla u(x)|}^{p-1}{)(a(x)|}^{\frac{1}{p}}|\nabla v(x)|)dx\hfill \\ & & \le {\left({\int }_{\Omega }a(x){|\nabla u(x)|}^{p}dx\right)}^{\frac{p-1}{p}}{\left({\int }_{\Omega }a(x){|\nabla v(x)|}^{p}dx\right)}^{\frac{1}{p}}<\infty ,\forall u,v\in W_0^{1,p}(a,\Omega ).\hfill \end{array}$$

The ordinary p-Laplacian is recovered when $a(x)\equiv 1$ in $\Omega$.

The degenerated p-Laplacian ${\Delta }_p^a:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ is continuous and bounded. We denote by ${\lambda }_1$ the first eigenvalue of $-{\Delta }_p^a:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ corresponding to the weight $a\in L^1(\Omega )$ with $a^{-\frac{1}{p-1}}\in L^1(\Omega )$. Specifically, ${\lambda }_1$ is the least $\lambda >0$ for which the problem

$$\left\{\begin{array}{cc}-{\mathrm{div}(a(x)|\nabla u|}^{p-2}{\nabla u)=\lambda |u|}^{p-2}u\hfill & \mathrm{in}\Omega \hfill \\ u=0\hfill & \mathrm{on}\partial \Omega \hfill \end{array}\right.$$

possesses a nontrivial solution. It can be variationally characterized as

$$\begin{array}{c}\hfill {\lambda }_1=\underset{u\in W_0^{1,p}(a,\Omega )\setminus \left\{0\right\}}{inf}\frac{{\int }_{\Omega }a(x){|\nabla u(x)|}^{p}dx}{{\parallel u\parallel }_{L^p(\Omega )}^p}.\end{array}$$

(5)

More details on the degenerated p-Laplacian with weight can be seen in [3].

For the positive weights $a\in L^1(\Omega )$ and $b\in L^1(\Omega )$ entering problems (1) and (2), we have the degenerated p-Laplacian ${\Delta }_p^a:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ with weight $a\in L^1(\Omega )$ and the degenerated q-Laplacian ${\Delta }_q^b:W_0^{1,2}(b,\Omega )\to W_0^{1,2}{(b,\Omega )}^{*}$ with weight $b\in L^1(\Omega )$. The two operators need to be consistent, which is achieved under the following compatibility condition for the weights:

(H2).

$1<q<p<+\infty$ and $a^{-\frac{q}{p-q}}{b}^{\frac{p}{p-q}}\in L^1(\Omega )$.

Proposition 2.

Assume that condition (H2) holds. Then, one has the continuous embedding $W_0^{1,p}(a,\Omega )\hookrightarrow W_0^{1,q}(b,\Omega )$.

Proof. 

By hypothesis (H2) and Hölder’s inequality, we infer that

$$\begin{array}{ccc}& & {\int }_{\Omega }{b(x)|\nabla u(x)|}^{q}dx={\int }_{\Omega }(a{(x)}^{-\frac{q}{p}}b(x)){(a(x)}^{\frac{q}{p}}{|\nabla u(x)|}^q)dx\hfill \\ & & \le {\left({\int }_{\Omega }a{(x)}^{-\frac{q}{p-q}}b{(x)}^{\frac{p}{p-q}}dx\right)}^{\frac{p-q}{p}}{\left({\int }_{\Omega }a(x){|\nabla u(x)|}^{p}dx\right)}^{\frac{q}{p}}\hfill \\ & & \le \parallel a^{-\frac{q}{p-q}}{b}^{\frac{p}{p-q}}{\parallel }_{L^1(\Omega )}^{\frac{p-q}{p}}{\parallel u\parallel }_{W^{1,p}(a,\Omega )}^q,\forall u,v\in W_0^{1,p}(a,\Omega ),\hfill \end{array}$$

which proves the result. □

Under condition (H2), on the basis of Proposition 2, the map ${\Delta }_p^a+{\Delta }_q^b:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ called the degenerated $(p,q)$-Laplacian with weights $a,b\in L^1(\Omega )$ is well-defined. It is given by

$$\begin{array}{ccc}\hfill & & {\langle -({\Delta }_p^a+{\Delta }_q^b)u,v\rangle }_{W_0^{1,p}(a,\Omega )}\hfill \\ & & ={\int }_{\Omega }{(a(x)|\nabla u(x)|}^{p-2}\nabla u(x)+{b(x)|\nabla u(x)|}^{q-2}\nabla u(x))\nabla v(x)dx,\forall u,v\in W_0^{1,p}(a,\Omega ).\hfill \end{array}$$

(6)

The degenerated $(p,q)$-Laplacian with weights $a,b\in L^1(\Omega )$ was introduced in [1].

Again on the basis of Proposition 2, the map ${\Delta }_p^a-{\Delta }_q^b:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ given by

$$\begin{array}{ccc}\hfill & & {\langle -({\Delta }_p^a-{\Delta }_q^b)u,v\rangle }_{W_0^{1,p}(a,\Omega )}\hfill \\ & & ={\int }_{\Omega }{(a(x)|\nabla u(x)|}^{p-2}\nabla u(x)+{b(x)|\nabla u(x)|}^{q-2}\nabla u(x))\nabla v(x)dx,\forall u,v\in W_0^{1,p}(a,\Omega ),\hfill \end{array}$$

(7)

is well-defined provided condition (H2) is satisfied. We call it the competing $(p,q)$-Laplacian with weights $a,b\in L^1(\Omega )$ and is introduced here for the first time.

Proposition 3.

Under assumption (H2), the maps ${\Delta }_p^a+{\Delta }_q^b:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ and ${\Delta }_p^a-{\Delta }_q^b:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ are continuous and bounded. In addition, under (H1) and (H2), the ${(S)}_{+}$ property holds for the map $-({\Delta }_p^a+{\Delta }_q^b):W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$; that is, any sequence $\left\{u_n\right\}\subset W_0^{1,p}(a,\Omega )$ satisfying $u_n\rightharpoonup u$ in $W_0^{1,p}(a,\Omega )$ and

$$\begin{array}{c}\hfill \underset{n\to \infty }{lim\; sup}{\langle -({\Delta }_p^a+{\Delta }_q^b)u_n,u_n-u\rangle }_{W_0^{1,p}(a,\Omega )}\le 0\end{array}$$

(8)

is strongly convergent. Thus, $u_n\to u$ in $W_0^{1,p}(a,\Omega )$.

Proof. 

Due to the continuous embedding $W^{1,p}(a,\Omega )\hookrightarrow W^{1,q}(b,\Omega )$ in Proposition 2, ${\Delta }_p^a+{\Delta }_q^b$ and ${\Delta }_p^a-{\Delta }_q^b$ inherit the continuity and boundedness from ${\Delta }_p^a$ and ${\Delta }_q^b$.

For the second part of the statement, let a sequence $\left\{u_n\right\}\subset W_0^{1,p}(a,\Omega )$ with the required properties. By (6), the monotonicity of $-{\Delta }_q^b$ and Hölder’s inequality, we obtain

$$\begin{array}{ccc}& & {\langle -({\Delta }_p^a+{\Delta }_q^b)(u_n)+({\Delta }_p^a+{\Delta }_q^b)(u),u_n-u\rangle }_{W_0^{1,p}(a,\Omega )}\hfill \\ & & \ge {\langle -{\Delta }_p^a(u_n)+{\Delta }_p^a(u),u_n-u\rangle }_{W_0^{1,p}(a,\Omega )}\hfill \\ & & \ge (\parallel u_n{\parallel }_{W_0^{1,p}(a,\Omega )}-{\parallel u\parallel }_{W_0^{1,p}(a,\Omega )})(\parallel u_n{\parallel }_{W_0^{1,p}(a,\Omega )}^{p-1}-{\parallel u\parallel }_{W_0^{1,p}(a,\Omega )}^{p-1})\ge 0.\hfill \end{array}$$

It follows from the above estimate, (8) and $u_n\rightharpoonup u$ in $W_0^{1,p}(a,\Omega )$ that there holds $\underset{n\to +\infty }{lim}\parallel u_n{\parallel }_{W_0^{1,p}(a,\Omega )}={\parallel u\parallel }_{W_0^{1,p}(a,\Omega )}$. From Proposition 1, we know that the space $W_0^{1,p}(a,\Omega )$ is uniformly convex. Therefore, we can conclude that $u_n\to u$ in $W_0^{1,p}(a,\Omega )$. □

3. An Associated Nemytskij Operator

In this section we focus on the right-hand side of the equations in (1) and (2), i.e., the convection term $f(x,u,\nabla u)$. Our goal is to identify the growth condition for the function $f(x,t,\xi )$ to match the reduction in Proposition 1 to the unweighted Sobolev space $W_0^{1,p_s}(\Omega )$. The appropriate growth for $f(x,t,\xi )$ is the one used in [1].

In order to simplify the presentation, for any real number $r>1$, we denote $r^{\prime }:=r/(r-1)$ (the Hölder conjugate of r). This convention will be preserved for the rest of the paper.

Lemma 1.

Assume (H1) and (H2) and in addition that the Carathéodory function $f:\Omega \times \mathbb{R}\times {\mathbb{R}}^N\to \mathbb{R}$ satisfies the growth condition:

(H3). 

$$\begin{array}{c}\hfill \left|f(x,t,\xi )\right|\le \sigma (x)+c_1{\left|t\right|}^{\alpha }+c_2{\left|\xi \right|}^{\beta }\mathit{for}a.ex\in \Omega ,\forall (t,\xi )\in \mathbb{R}\times {\mathbb{R}}^N,\end{array}$$

(9)

with $\sigma \in L^{{\gamma }^{\prime }}(\Omega )$ for $\gamma \in (1,p_s^{*})$ and constants $c_1>0$, $c_2>0$, $\alpha \in [0,p_s^{*}-1)$, $\beta \in [0,\frac{p_s}{{(p_s^{*})}^{\prime }})$.

Set

$$\begin{array}{c}\hfill \theta :=min\left\{{\gamma }^{\prime },\frac{p_s^{*}}{\alpha },\frac{p_s}{\beta }\right\}.\end{array}$$

(10)

Then, the Nemytskij operator ${\mathcal{N}}_f:L^{p_s^{*}}(\Omega )\times {(L^{p_s}(\Omega ))}^N\to L^{\theta }(\Omega )$ associated with the function f which is given by

$$\begin{array}{c}\hfill {\mathcal{N}}_f(v,z)=f(·,v(·),z(·)),\forall v,z\in L^{p_s^{*}}(\Omega )\times {(L^{p_s}(\Omega ))}^N\end{array}$$

(11)

is well-defined, continuous and bounded.

Proof. 

The requirements in (H3) postulate ${\gamma }^{\prime }>1$ (note $\gamma >1$),

$$\begin{array}{c}\hfill \frac{p_s^{*}}{\alpha }>\frac{p_s^{*}}{p_s^{*}-1}={(p_s^{*})}^{\prime }>1,\frac{p_s}{\beta }>{(p_s^{*})}^{\prime }>1\end{array}$$

(note that $p_s^{*}>p>1$). Then, a consequence of (10) is that $\theta >1$.

We observe that (9) yields

$$\begin{array}{c}\hfill \left|f(x,t,\xi )\right|\le \tilde{\sigma }(x)+c_1{\left|t\right|}^{\frac{p_s^{*}}{\theta }}+c_2{\left|\xi \right|}^{\frac{p_s}{\theta }}\mathrm{for}\mathrm{a}.\mathrm{e}x\in \Omega ,\forall (t,\xi )\in \mathbb{R}\times {\mathbb{R}}^N,\end{array}$$

(12)

with $\tilde{\sigma }\in L^{\theta }(\Omega )$. Indeed, (10) gives

$$\begin{array}{c}\hfill \alpha \le \frac{p_s^{*}}{\theta }and\beta \le \frac{p_s}{\theta }.\end{array}$$

Hence (12) is derived from (9) with $\tilde{\sigma }(x)=\sigma (x)+c_1+c_2$ for a.e. $x\in \Omega$ obtaining $\tilde{\sigma }\in L^{{\gamma }^{\prime }}(\Omega )\subset L^{\theta }(\Omega )$.

Using Krasnoselskij’s theorem, we infer from (12) that ${\mathcal{N}}_f$ introduced in (11) has the required properties, thus proving the result. □

Let $N_f:W^{1,p}(a,\Omega )\to W^{1,p}{(a,\Omega )}^{*}$ be defined by

$$\begin{array}{c}\hfill {\langle N_f(u),v\rangle }_{W^{1,p}(a,\Omega )}={\int }_{\Omega }f(x,u(x),\nabla u(x))v(x)dx,\forall u,v\in W_0^{1,p}(a,\Omega ).\end{array}$$

Due to the first inclusion in (3), it holds that $(u,\nabla u)\in L^{p_s^{*}}(\Omega )\times {(L^{p_s}(\Omega ))}^N$ whenever $u\in W_0^{1,p}(a,\Omega )$. It turns out

$$\begin{array}{c}\hfill {\langle N_f(u),v\rangle }_{W^{1,p}(a,\Omega )}={\langle {\mathcal{N}}_f(u,\nabla u),v\rangle }_{L^{\theta }(\Omega )},\forall u,v\in W_0^{1,p}(a,\Omega ).\end{array}$$

(13)

The assertion below provides a key tool for investigating problems (1) and (2).

Proposition 4.

Assume (H1)–(H3). If $u_n\rightharpoonup u$ in $W_0^{1,p}(a,\Omega )$, it holds that

$$\begin{array}{c}\hfill \underset{n\to \infty }{lim}{\langle N_f(u_n),u_n-u\rangle }_{W_0^{1,p}(a,\Omega )}=0.\end{array}$$

(14)

Proof. 

Recalling the convention made in the beginning of this section, (10) entails

$$\begin{array}{c}\hfill {\theta }^{\prime }=max\left\{\gamma ,{\left(\frac{p_s^{*}}{\alpha }\right)}^{\prime },{\left(\frac{p_s}{\beta }\right)}^{\prime }\right\}.\end{array}$$

(15)

By (H3) we have $\gamma <p_s^{*}$,

$$\begin{array}{c}\hfill {\left(\frac{p_s^{*}}{\alpha }\right)}^{\prime }=\frac{\frac{p_s^{*}}{\alpha }}{\frac{p_s^{*}}{\alpha }-1}=\frac{p_s^{*}}{p_s^{*}-\alpha }<p_s^{*},\end{array}$$

$$\begin{array}{c}\hfill {\left(\frac{p_s}{\beta }\right)}^{\prime }=\frac{\frac{p_s}{\beta }}{\frac{p_s}{\beta }-1}=\frac{p_s}{p_s-\beta }<\frac{p_s}{p_s-\frac{p_s}{{(p_s^{*})}^{\prime }}}=\frac{{(p_s^{*})}^{\prime }}{{(p_s^{*})}^{\prime }-1}=p_s^{*}.\end{array}$$

Hence, (15) yields $1<{\theta }^{\prime }<p_s^{*}$ and we can apply the Rellich–Kondrachov compact embedding theorem to deduce that the embedding $W^{1,p_s}(\Omega )\hookrightarrow L^{{\theta }^{\prime }}(\Omega )$ is compact, which results in $u_n\to u$ in $L^{{\theta }^{\prime }}(\Omega )$.

Since (13) implies that

$$\begin{array}{c}\hfill |{\langle N_f(u_n),u_n-u\rangle }_{W^{1,p}(a,\Omega )}|\le \parallel {\mathcal{N}}_f(u_n,\nabla u_n){\parallel }_{L^{\theta }(\Omega )}{\parallel u_n-u\parallel }_{L^{{\theta }^{\prime }}(\Omega )}.\end{array}$$

(16)

and Lemma 1 ensures that $\left\{{\mathcal{N}}_f(u_n,\nabla u_n)\right\}$ is bounded in $L^{\theta }(\Omega )$, from (16) we arrive at (14), as desired. □

4. Solvability and Approximation for the Degenerate Elliptic Problem (1)

The object of this section is to develop an approach based on finite dimensional approximations for problem (1).

Since the Banach space $W_0^{1,p}(a,\Omega )$ is separable (see Section 2), there exists a Galerkin basis for it. This amounts to saying that there is a sequence $\left\{X_n\right\}$ of vector subspaces of $W_0^{1,p}(a,\Omega )$ such that

(i)

$\dim (X_n)<\infty ,\forall n;$

(ii)

$X_n\subset X_{n+1},\forall n;$

(iii)

$$\overline{\bigcup _n{X}_n}=W_0^{1,p}(a,\Omega ).$$

We fix such a sequence of subspaces $\left\{X_n\right\}$. Each approximate problem on $X_n$ will be resolved by means of a consequence of Brouwer’s fixed point theorem.

Proposition 5.

Assume the conditions (H1)–(H3) and in addition

(H4).

there exists $\rho \in L^1(\Omega )$ and constants $d_1>0$ and $d_2>0$ provided ${\lambda }_1^{-1}{d}_1+d_2<1$, where ${\lambda }_1$ denotes the first eigenvalue of $-{\Delta }_p^a$ on $W_0^{1,p}(a,\Omega )$, such that

$$\begin{array}{c}\hfill f(x,t,\xi )t\le \rho (x)+d_1{\left|t\right|}^p+d_{2}a(x){\left|\xi \right|}^p\end{array}$$

(17)

for a.e $x\in \Omega$ and all $(t,\xi )\in \mathbb{R}\times {\mathbb{R}}^N$.

Then for each n there exists $u_n\in X_n$ such that

$${\langle -({\Delta }_p^a+{\Delta }_q^b)(u_n),v\rangle }_{W_0^{1,p}(a,\Omega )}={\int }_{\Omega }f(x,u_n(x),\nabla u_n(x))v(x)dx,\forall v\in X_n.$$

(18)

Proof. 

For each n, consider the continuous map $A_n:X_n\to X_n^{*}$ defined by

$${\langle A_n(u),v\rangle }_{X_n}={\langle -({\Delta }_p^a+{\Delta }_q^b)(u),v\rangle }_{W_0^{1,p}(a,\Omega )}-{\int }_{\Omega }f(x,u(x),\nabla u(x))v(x)dx,\forall v\in X_n.$$

The definition of the operator $A_n$, (17) and (5) lead to

$$\begin{array}{ccc}& & {\langle A_n(v),v\rangle }_{X_n}={\int }_{\Omega }{(a(x)|\nabla v|}^p+{b(x)|\nabla v|}^q-f(x,v,\nabla v)v)dx\hfill \\ & & \ge {\parallel v\parallel }_{W_0^{1,p}(a,\Omega )}^p-{\parallel \rho \parallel }_{L^1(\Omega )}-d_1{\parallel v\parallel }_{L^p(\Omega )}^p-d_2{\parallel v\parallel }_{W_0^{1,p}(a,\Omega )}^p\hfill \\ & & \ge (1-d_1{\lambda }_1^{-1}-d_2){\parallel v\parallel }_{W_0^{1,p}(a,\Omega )}^p-{\parallel \rho \parallel }_{L^1(\Omega )},\forall v\in X_n.\hfill \end{array}$$

Thanks to the assumption $1-d_1{\lambda }_1^{-1}-d_2>0$ in (H4), it follows that

$$\begin{array}{c}\hfill {\langle A_n(v),v\rangle }_{X_n}\ge 0\mathrm{whenever}v\in X_n\mathrm{with}{\parallel v\parallel }_{W_0^{1,p}(a,\Omega )}=R\end{array}$$

provided $R=R(n)>0$ is sufficiently large. In view of the fact that $X_n$ is a finite dimensional space, by a well-known consequence of Brouwer’s fixed point theorem (see, e.g., [10] (p. 37)) there exists $u_n\in X_n$ solving the equation $A_n(u_n)=0$. This means exactly that $u_n\in X_n$ is a solution for problem (18), which completes the proof. □

We are in a position to state our main result on problem (1).

Theorem 1.

Assume that the conditions (H1)–(H4) are fulfilled. Then, the sequence $\left\{u_n\right\}$, with $u_n\in X_n$ constructed in Proposition 5, contains a subsequence which is strongly convergent in $W_0^{1,p}(a,\Omega )$ to a weak solution of problem (1) meaning that

$$\begin{array}{c}\hfill {\int }_{\Omega }{(a(x)|\nabla u|}^{p-2}\nabla u+{b(x)|\nabla u|}^{q-2}\nabla u)\nabla vdx={\int }_{\Omega }f(x,u,\nabla u)vdx\end{array}$$

(19)

for all $v\in W_0^{1,p}(a,\Omega )$.

Proof. 

We claim that the sequence $\left\{u_n\right\}$ built in Proposition 5 is bounded in $W_0^{1,p}(a,\Omega )$. Acting with $v=u_n$ in (18) gives

$$\begin{array}{c}\hfill \parallel u_n{\parallel }_{W_0^{1,p}(a,\Omega )}^p+{\parallel u_n\parallel }_{W_0^{1,q}(a,\Omega )}^q={\int }_{\Omega }f(x,u_n,\nabla u_n)u_{n}dx.\end{array}$$

Then, through (17) and (5) we obtain

$$\begin{array}{ccc}& & \parallel u_n{\parallel }_{W_0^{1,p}(a,\Omega )}^p\le {\parallel \rho \parallel }_{L^1(\Omega )}+d_1{\parallel u\parallel }_{L^p(\Omega )}^p+d_2{\parallel u\parallel }_{W_0^{1,p}(a,\Omega )}^p\hfill \\ & & \le {\parallel \rho \parallel }_{L^1(\Omega )}+(d_1{\lambda }_1^{-1}+d_2){\parallel u\parallel }_{W_0^{1,p}(a,\Omega )}^p.\hfill \end{array}$$

Thanks to ${\lambda }_1^{-1}{d}_1+d_2<1$, as known from hypothesis (H4), the claim is verified.

Recall from Proposition 1 that $W_0^{1,p}(a,\Omega )$ is a uniformly convex Banach space, so it is reflexive. Hence, the bounded sequence $\left\{u_n\right\}$ possesses a subsequence denoted again $\left\{u_n\right\}$ such that for some $u\in W_0^{1,p}(a,\Omega )$ it holds $u_n\rightharpoonup u$ in $W_0^{1,p}(a,\Omega )$.

Proposition 3 and Lemma 1 ensure that the operators ${\Delta }_p^a+{\Delta }_q^b:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ and $N_f:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ are bounded. Then, in view of the reflexivity of $W_0^{1,p}(a,\Omega )$ along a relabeled subsequence, one has

$$\begin{array}{c}\hfill -({\Delta }_p^a+{\Delta }_q^b)(u_n)-N_f(u_n)\rightharpoonup \eta \mathrm{in}{W}_0^{1,p}{(a,\Omega )}^{*}\end{array}$$

(20)

for some $\eta \in W_0^{1,p}{(a,\Omega )}^{*}$.

Let us prove that $\eta =0$. For $v\in {\displaystyle \bigcup _n}{X}_n$ choose m with $v\in X_m$. According to Proposition 5 and property (ii) in the definition of Galerkin basis, we may apply (18) for all $n\ge m$, which reads as

$$\begin{array}{c}\hfill {\langle -({\Delta }_p^a+{\Delta }_q^b)(u_n)-N_f(u_n),v\rangle }_{W_0^{1,p}(a,\Omega )}=0\mathrm{for}\mathrm{all}n\ge m.\end{array}$$

Letting $n\to \infty$ enables us to derive from (20) that

$$\begin{array}{c}\hfill {\langle \eta ,v\rangle }_{W_0^{1,p}(a,\Omega )}=0,\forall v\in \bigcup _n{X}_n.\end{array}$$

The property (iii) in the definition of Galerkin basis $\left\{X_n\right\}$ highlights the density of the set ${\displaystyle \bigcup _n}{X}_n$ in $W_0^{1,p}(a,\Omega )$. As $\eta$ vanishes on ${\displaystyle \bigcup _n}{X}_n$, it follows that $\eta =0$.

Therefore, (20) becomes

$$\begin{array}{c}\hfill -({\Delta }_p^a+{\Delta }_q^b)(u_n)-N_f(u_n)\rightharpoonup 0\mathrm{in}{W}_0^{1,p}{(a,\Omega )}^{*}.\end{array}$$

(21)

In particular, we have

$$\underset{n\to \infty }{lim}{\langle -({\Delta }_p^a+{\Delta }_q^b)(u_n)-N_f(u_n),u\rangle }_{W_0^{1,p}(a,\Omega )}=0.$$

(22)

Now, we return to (18) and insert $v=u_n$, obtaining

$${\langle -({\Delta }_p^a+{\Delta }_q^b)(u_n)-N_f(u_n),u_n\rangle }_{W_0^{1,p}(a,\Omega )}=0,\forall n,$$

which in conjunction with (22) yields

$$\underset{n\to \infty }{lim}{\langle -({\Delta }_p^a+{\Delta }_q^b)(u_n)-N_f(u_n),u_n-u\rangle }_{W_0^{1,p}(a,\Omega )}=0.$$

Taking into account Proposition 4, this amounts to saying that

$$\underset{n\to \infty }{lim}{\langle -({\Delta }_p^a+{\Delta }_q^b)(u_n),u_n-u\rangle }_{W_0^{1,p}(a,\Omega )}=0.$$

Consequently, the sequence $\left\{u_n\right\}$ satisfies (8). We are thus allowed to apply Proposition 3 which provides the strong convergence $u_n\to u$ in $W_0^{1,p}(a,\Omega )$.

Using the continuity of the nonlinear operators ${\Delta }_p^a+{\Delta }_q^b:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ and $N_f:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ as known by Proposition 2 and Lemma 1, we infer from the strong convergence $u_n\to u$ in $W_0^{1,p}(a,\Omega )$ that

$$\begin{array}{c}\hfill -({\Delta }_p^a+{\Delta }_q^b)(u_n)-N_f(u_n)\to -({\Delta }_p^a+{\Delta }_q^b)(u)-N_f(u)\mathrm{in}{W}_0^{1,p}{(a,\Omega )}^{*}.\end{array}$$

A simple comparison with (21) confirms that

$$\begin{array}{c}\hfill -({\Delta }_p^a+{\Delta }_q^b)(u)-N_f(u)=0,\end{array}$$

which is just (19). The proof is complete. □

5. Resolving the Non-Elliptic Problem (2)

Due to the total lack of ellipticity of the competing $(p,q)$-Laplacian ${\Delta }_p^a-{\Delta }_q^b$ with weights $a\in L^1(\Omega )$ and $b\in L^1(\Omega )$ as introduced in (7), i.e., the differential operator $\mathrm{div}(a(x)|\nabla u{|}^{p-2}\nabla u-b(x)|\nabla u{|}^{q-2}\nabla u)$, when the weights $a(x)$ and $b(x)$ are positive, we are not able to prove the existence of a weak solution for problem (2) in the weak sense. For this reason, we seek a solution in the following generalized sense.

Definition 1.

An element $u\in W_0^{1,p}(a,\Omega )$ is called a generalized solution to problem (2) if there exists a sequence $\left\{u_n\right\}\subset W_0^{1,p}(a,\Omega )$ such that

(j)

$u_n\rightharpoonup u$ in $W_0^{1,p}(\Omega )$;

(jj)

for every $v\in W_0^{1,p}(a,\Omega )$, it holds that

$$\begin{array}{c}\hfill \underset{n\to \infty }{lim}{\int }_{\Omega }((a(x)|\nabla u_n{|}^{p-2}\nabla u_n-b(x)|\nabla u_n{|}^{q-2}\nabla u_n)\nabla v-f(x,u_n,\nabla u_n)v)dx=0;\end{array}$$

(jjj)

$$\begin{array}{c}\hfill \underset{n\to \infty }{lim}{\int }_{\Omega }((a(x)|\nabla u_n{|}^{p-2}\nabla u_n-b(x)|\nabla u_n{|}^{q-2}\nabla u_n)\nabla (u_n-u)dx=0.\end{array}$$

Our result for the non-elliptic problem (2) is as follows.

Theorem 2.

Assume for $1<q<p<+\infty$ that the positive weights $a\in L^1(\Omega )$ and $b\in L^1(\Omega )$ and the Carathéodory function $f:\Omega \times \mathbb{R}\times {\mathbb{R}}^N\to \mathbb{R}$ that the conditions (H1)–(H4) hold. Then, there exists at least a generalized solution $u\in W_0^{1,p}(a,\Omega )$ of problem (2) in the sense of Definition 1.

Proof. 

The proof is carried over along the pattern of Theorem 1. Fix a Galerkin basis $\left\{X_n\right\}$ of $W_0^{1,p}(a,\Omega )$, i.e., a sequence of finite dimensional vector subspaces of $W_0^{1,p}(a,\Omega )$ such that the properties (i)–(iii) in Section 4 hold.

We claim that for each n there exists $u_n\in X_n$ such that

$${\int }_{\Omega }((a(x)|\nabla u_n{|}^{p-2}\nabla u_n-b(x)|\nabla u_n{|}^{q-2}\nabla u_n)\nabla v-f(x,u_n,\nabla u_n)v)dx=0,\forall v\in X_n.$$

(23)

To this end, define the continuous map $B_n:X_n\to X_n^{*}$ by

$${\langle B_n(u),v\rangle }_{X_n}={\int }_{\Omega }{((a(x)|\nabla u|}^{p-2}\nabla u-{b(x)|\nabla u|}^{q-2}\nabla u)\nabla v-f(x,u,\nabla u)v)dx$$

for all $u,v\in X_n$. The continuous embedding $W^{1,p}(a,\Omega )\hookrightarrow W^{1,q}(b,\Omega )$ (see Proposition 2), (17) and (5) imply that

$$\begin{array}{c}\hfill {\langle B_n(v),v\rangle }_{X_n}\ge (1-d_1{\lambda }_1^{-1}-d_2){\parallel v\parallel }_{W_0^{1,p}(a,\Omega )}^p-{C\parallel v\parallel }_{W_0^{1,p}(a,\Omega )}^q-{\parallel \rho \parallel }_{L^1(\Omega )},\forall v\in X_n.\end{array}$$

with a constant $C>0$. By the assumption $1-d_1{\lambda }_1^{-1}-d_2>0$ in (H4) and the fact that $p>q$, it turns out

$$\begin{array}{c}\hfill {\langle B_n(v),v\rangle }_{X_n}\ge 0\mathrm{for}\mathrm{all}v\in X_n\mathrm{with}{\parallel v\parallel }_{W_0^{1,p}(a,\Omega )}=R\end{array}$$

if $R=R(n)>0$ is sufficiently large. According to condition (i) in the definition of Galerkin basis, the space $X_n$ is finite dimensional. This enables us to apply a well-known consequence of Brouwer’s fixed point theorem (see, e.g., [10] (p. 37)) obtaining a $u_n\in X_n$ with $B_n(u_n)=0$. Therefore, we obtain (23), thus proving the claim.

Next, we show that the sequence $\left\{u_n\right\}$ is bounded in $W_0^{1,p}(a,\Omega )$. Since $u_n\in X_n$, we can take $v=u_n$ as a test function in (23), where

$$\begin{array}{c}\hfill \parallel u_n{\parallel }_{W_0^{1,p}(a,\Omega )}^p={\parallel u_n\parallel }_{W_0^{1,q}(a,\Omega )}^q+{\int }_{\Omega }f(x,u_n,\nabla u_n)u_{n}dx.\end{array}$$

(24)

The continuous embedding $W^{1,p}(a,\Omega )\hookrightarrow W^{1,q}(b,\Omega )$ in Proposition 2, in conjunction with (17) and (5), ensures the estimate

$$\begin{array}{ccc}& & \parallel u_n{\parallel }_{W_0^{1,p}(a,\Omega )}^p\le C\parallel u_n{\parallel }_{W_0^{1,p}(a,\Omega )}^q+{\parallel \rho \parallel }_{L^1(\Omega )}+d_1{\parallel u\parallel }_{L^p(\Omega )}^p+d_2{\parallel u\parallel }_{W_0^{1,p}(a,\Omega )}^p\hfill \\ & & \le C\parallel u_n{\parallel }_{W_0^{1,p}(a,\Omega )}^q+{\parallel \rho \parallel }_{L^1(\Omega )}+(d_1{\lambda }_1^{-1}+d_2){\parallel u\parallel }_{W_0^{1,p}(a,\Omega )}^p,\hfill \end{array}$$

with a constant $C>0$. On account of $p>q$ and assumption ${\lambda }_1^{-1}{d}_1+d_2<1$ in (H4), we conclude that the sequence $\left\{u_n\right\}$ is bounded in $W_0^{1,p}(a,\Omega )$.

Proposition 1 guarantees the reflexivity of the space $W_0^{1,p}(a,\Omega )$. We are thus allowed to extract a subsequence still denoted as $\left\{u_n\right\}$ such that $u_n\rightharpoonup u$ in $W_0^{1,p}(\Omega )$ for some $u\in W_0^{1,p}(a,\Omega )$. The requirement (j) in Definition 1 is fulfilled.

Equality (23) expresses that

$$\begin{array}{c}\hfill {\langle -{\Delta }_p^a(u_n)+{\Delta }_q^b(u_n)-N_f(u_n),v\rangle }_{W_0^{1,p}(a,\Omega )}=0,\forall v\in X_n.\end{array}$$

(25)

Inserting $v=u_n$ in (25) leads to

$$\begin{array}{c}\hfill {\langle -{\Delta }_p^a(u_n)+{\Delta }_q^b(u_n)-N_f(u_n),u_n\rangle }_{W_0^{1,p}(a,\Omega )}=0,\forall n.\end{array}$$

(26)

The sequence $\left\{(-{\Delta }_p^a+{\Delta }_q^b-N_f)(u_n)\right\}$ is bounded in $W_0^{1,p}{(a,\Omega )}^{*}$ because the nonlinear operators ${\Delta }_p^a,{\Delta }_q^b,N_f:W_0^{1,p}(a,\Omega )\to W_0^{1,p}{(a,\Omega )}^{*}$ are bounded. Due to the reflexivity of the space $W_0^{1,p}{(\Omega )}^{*}$, we can pass to a relabeled subsequence such that for a $\zeta \in W_0^{1,p}{(\Omega )}^{*}$ it holds that

$$\begin{array}{c}\hfill (-{\Delta }_p^a+{\Delta }_q^b-N_f)(u_n)\rightharpoonup \zeta \mathrm{in}{W}_0^{1,p}{(\Omega )}^{*}.\end{array}$$

(27)

Let $v\in X_m$ for some m. Assertion (ii) in the definition of Galerkin basis renders $v\in X_n$ for every $n\ge m$. Then, (25) and (27) imply

$$\begin{array}{c}\hfill {\langle \zeta ,v\rangle }_{W_0^{1,p}(a,\Omega )}=0.\end{array}$$

By (iii) in the definition of Galerkin basis $\left\{X_n\right\}$, the set ${\displaystyle \bigcup _n}{X}_n$ is dense in $W_0^{1,p}(a,\Omega )$. Therefore, $\zeta =0$, so that (27) becomes

$$\begin{array}{c}\hfill (-{\Delta }_p^a+{\Delta }_q^b-N_f)(u_n)\rightharpoonup 0\mathrm{in}{W}_0^{1,p}{(\Omega )}^{*},\end{array}$$

which establishes property (jj) in Definition 1.

Setting $v=u$ in (jj) provides

$$\begin{array}{c}\hfill \underset{n\to \infty }{lim}{\langle -{\Delta }_p^a(u_n)+{\Delta }_q^b(u_n)-N_f(u_n),u\rangle }_{W_0^{1,p}(a,\Omega )}=0,\end{array}$$

which, with (24), produces

$$\begin{array}{c}\hfill \underset{n\to \infty }{lim}{\langle -{\Delta }_p^a(u_n)+{\Delta }_q^b(u_n)-N_f(u_n),u_n-u\rangle }_{W_0^{1,p}(a,\Omega )}=0.\end{array}$$

(28)

Proposition 4 and (28) ensure that

$$\begin{array}{c}\hfill \underset{n\to \infty }{lim}{\langle -{\Delta }_p^a(u_n)+{\Delta }_q^b(u_n),u_n-u\rangle }_{W_0^{1,p}(a,\Omega )}=0\end{array}$$

which shows the validity of part (jjj) in Definition 1. Summarizing, $u\in W_0^{1,p}(a,\Omega )$ is a generalized solution to problem (2) in the sense of Definition 1. □

Remark 1.

The notion of a generalized solution can be introduced for problem (2), too. Precisely, $u\in W_0^{1,p}(a,\Omega )$ is called a generalized solution to problem (1) if there exists a sequence $\left\{u_n\right\}\subset W_0^{1,p}(a,\Omega )$ such that (j) in Definition 1 holds with

(jj)′

for every $v\in W_0^{1,p}(a,\Omega )$ one has

$$\begin{array}{c}\hfill \underset{n\to \infty }{lim}{\int }_{\Omega }((a(x)|\nabla u_n{|}^{p-2}\nabla u_n+b(x)|\nabla u_n{|}^{q-2}\nabla u_n)\nabla v-f(x,u_n,\nabla u_n)v)dx=0;\end{array}$$

(jjj)″

with $u\in W_0^{1,p}(a,\Omega )$ in (j),

$$\begin{array}{c}\hfill \underset{n\to \infty }{lim}{\int }_{\Omega }((a(x)|\nabla u_n{|}^{p-2}\nabla u_n+b(x)|\nabla u_n{|}^{q-2}\nabla u_n)\nabla (u_n-u)dx=0.\end{array}$$

In the case of problem (1), $u\in W_0^{1,p}(a,\Omega )$ is a generalized solution if and only if it is a weak solution in the sense of (19). Indeed, if $u\in W_0^{1,p}(a,\Omega )$ is a weak solution to problem (1), then the constant sequence $\{u_n=u\}\subset W_0^{1,p}(a,\Omega )$ verifies (j), (jj)′, (jjj)″; thus, u is a generalized solution. Conversely, let $u\in W_0^{1,p}(a,\Omega )$ be a generalized solution for (1) with the sequence $\left\{u_n\right\}\subset W_0^{1,p}(a,\Omega )$ satisfying (j), (jj)′, (jjj)″. Condition (jjj)″ reads as

$$\begin{array}{c}\hfill \underset{n\to \infty }{lim}{\langle -({\Delta }_p^a+{\Delta }_q^b)(u_n),u_n-u\rangle }_{W_0^{1,p}(a,\Omega )}=0.\end{array}$$

By (j) and Proposition 3 we deduce that $u_n\to u$ in $W_0^{1,p}(\Omega )$. Then from (jj)′ we get (19), where u is a weak solution of problem (1).

6. An Application

The goal of this section is to illustrate the effective applicability of our results. For the sake of simplicity, we focus on problems of types (1) and (2) on the unit open ball $B=\{x\in {\mathbb{R}}^3:|x|<1\}$ in ${\mathbb{R}}^3$ and for degenerated and competing $(3,2)$-Laplacians with weights.

Consider the Dirichlet problems

$$\left\{\begin{array}{cc}-{\mathrm{div}(|x|}^r{|\nabla u|\nabla u+(1-|x|}^2{)}^h\nabla u)=g(x,u)+k_0\frac{{\left|x\right|}^{r}u}{1+u^2}{|\nabla u|}^{\mu }-k_{1}u{|\nabla u|}^{\nu }\hfill & \mathrm{in}B\hfill \\ \\ u=0\hfill & \mathrm{on}\partial B\hfill \end{array}\right.$$

(29)

and

$$\left\{\begin{array}{cc}-{\mathrm{div}(|x|}^r{|\nabla u|\nabla u-(1-|x|}^2{)}^h\nabla u)=g(x,u)+k_0\frac{{\left|x\right|}^{r}u}{1+u^2}{|\nabla u|}^{\mu }-k_{1}u{|\nabla u|}^{\nu }\hfill & \mathrm{in}B\hfill \\ \\ u=0\hfill & \mathrm{on}\partial B\hfill \end{array}\right.$$

(30)

on B, with constants $r\in (0,\frac{3}{2})$, $h\ge 0$, $k_0\in |0,1)$, $k_1\ge 0$, $\mu \in [0,\frac{5}{3})$, $\nu \in [0,\frac{5}{6})$, and a Carathéodory function $g:\Omega \times \mathbb{R}\to \mathbb{R}$ satisfying

$$\begin{array}{c}\hfill \left|g(x,t)\right|\le a_0{t}^2+b_0\mathrm{for}\mathrm{a}.\mathrm{e}x\in \Omega ,\forall t\in \mathbb{R},\end{array}$$

with constants $a_0\ge 0$ and $b_0\ge 0$ provided $(a_0+b_0){\lambda }_1^{-1}+k_0<1$, where ${\lambda }_1$ represents the first eigenvalue of $-{\Delta }_3^a$ on $W_0^{1,3}(a,B)$ with $a(x)={\left|x\right|}^r$. Notice that (29) and (30) are particular cases of problems (1) and (2), respectively, with $N=3$, $p=3$, $q=2$, $\Omega =B$, $a(x)={\left|x\right|}^r$, $b(x)={(1-|x|}^2{)}^h$ and

$$\begin{array}{c}\hfill f(x,t,\xi )=g(x,t)+k_0\frac{{\left|x\right|}^{r}t}{1+t^2}{\left|\xi \right|}^{\mu }-k_{1}t{\left|\xi \right|}^{\nu }.\end{array}$$

Let us check the conditions (H1)–(H4). Condition (H1) requires having $a^{-s}={\left|x\right|}^{-rs}\in L^1(B)$ for some $s\in \left(max\left\{\frac{N}{p},\frac{1}{p-1}\right\},+\infty \right)=(1,+\infty )$, which amounts to choosing $1<s<\frac{3}{r}$. Taking into account that $r\in (0,\frac{3}{2})$, condition (H1) is fulfilled for instance with $s=2$, a choice that we keep in the sequel.

Since $2r<3$, we have

$$\begin{array}{c}\hfill a^{-\frac{q}{p-q}}{b}^{\frac{p}{p-q}}={\left|x\right|}^{-2r}{(1-|x|}^2{)}^{3h}\in L^1(B),\end{array}$$

therefore, assumption (H2) is verified. For $s=2$, it holds $p_s=ps/(s+1)=2$, so we are in the situation of $N=3>p_s=2$, where $p_s^{*}=\frac{N{p}_s}{N-p_s}=6$, so ${(p_s^{*})}^{\prime }=\frac{6}{5}$ and $\frac{p_s}{{(p_s^{*})}^{\prime }}=\frac{5}{3}$.

We note that

$$\begin{array}{ccc}& & \left|f(x,t,\xi )\right|\le a_0{t}^2+b_0+k_0{\left|x\right|}^r{\left|\xi \right|}^{\mu }+k_1{\left|t\right|\left|\xi \right|}^{\nu }\hfill \\ & & \le a_0{t}^2+b_0+k_0{\left|\xi \right|}^{\mu }+\frac{k_1}{2}(t^2+{\left|\xi \right|}^{2\nu })\hfill \\ & & \le b_0+k_0+\frac{k_1}{2}+(a_0+\frac{k_1}{2})t^2+(k_0+\frac{k_1}{2}){\left|\xi \right|}^{\gamma }\hfill \end{array}$$

for a.e $x\in B$ and all $(t,\xi )\in \mathbb{R}\times {\mathbb{R}}^3$, where

$$\begin{array}{c}\hfill \gamma :=max\{\mu ,2\nu \}<\frac{5}{3}=\frac{p_s}{{(p_s^{*})}^{\prime }}.\end{array}$$

Therefore assumption (H3) is satisfied with $\sigma (x)=b_0+k_0+\frac{k_1}{2}$, $c_1=a_0+\frac{k_1}{2}$, $c_2=k_0+\frac{k_1}{2}$, $\alpha =2$ and $\beta =\gamma$. We also derive

$$\begin{array}{ccc}& & f(x,t,\xi )t=g(x,t)t+k_0\frac{{\left|x\right|}^r{t}^2}{1+t^2}{\left|\right|}^{\mu }-k_1{t}^2{\left|\xi \right|}^{\nu }\hfill \\ & & \le a_0{\left|t\right|}^3+b_0\left|t\right|+k_0{\left|x\right|}^r{\left|\xi \right|}^{\mu }\hfill \\ & & \le b_0+k_0+(a_0+b_0){\left|t\right|}^3+k_0{\left|x\right|}^r{\left|\xi \right|}^3\hfill \end{array}$$

for a.e $x\in B$ and all $(t,\xi )\in \mathbb{R}\times {\mathbb{R}}^3$. Assumption (H4) is verified with $\rho (x)=b_0+k_0$, $d_1=a_0+b_0$ and $d_2=k_0$ having been supposed that ${\lambda }_1^{-1}(a_0+b_0)+k_0<1$.

Since the assumptions (H1)-(H4) are satisfied, Theorems 1 and 2 can be applied to ensure the existence of a weak solution to problem (29) and of a generalized solution to problem (30). The weak solution to problem (29) can be approximated as described in Theorem 1.