Skip to main content

Compatibility conditions for the existence of weak solutions to a singular elliptic equation

Abstract

This paper deals with the existence of positive solutions to the singular elliptic boundary value problem involving p-Laplace operator

$$-\operatorname{div}\bigl(|\nabla u|^{p-2}\nabla u\bigr)= \frac{h(x)}{u^{\alpha}}+k(x)u^{\beta},\ x\in\Omega;\ \ u(x)>0,\ x\in\Omega;\ \ u(x)=0,\ x \in \partial\Omega; $$

where \(\Omega\subset\mathbb{R}^{N}\) (\(N\geq1\)) is a bounded domain with smooth boundary Ω, \(h\in L^{1}(\Omega)\), \(h(x)>0\) almost everywhere in Ω, \(k\in L^{\infty}(\Omega)\) is nonnegative, \(p>2\), \(\alpha>1\) and \(\beta\in(0,p-1)\). A compatibility condition on the couple \((h(x),\alpha)\) is given for the problem to have at least one solution. More precisely, it is shown that the problem admits a solution if and only if there exists \(u_{0}\in H_{0}^{1}(\Omega)\) such that \(\int_{\Omega}hu_{0}^{1-\alpha}\,\mathrm{d}x<\infty\).

1 Introduction

In this paper, we consider the following quasilinear elliptic equation:

$$ \begin{cases} -\Delta_{p} u=\frac{h(x)}{u^{\alpha}}+k(x)u^{\beta},&x\in\Omega,\\ u(x)=0,&x\in\partial\Omega,\\ u(x)>0,&x\in\Omega, \end{cases} $$
(1)

where Ω is a bounded domain in \(\mathbb{R}^{N}\) (\(N\geq1\)) with smooth boundary Ω, \(\Delta_{p} u=\operatorname{div}(|\nabla u|^{p-2}\nabla u)\) is the standard p-Laplace operator with \(p>2\), \(h\in L^{1}(\Omega)\), \(h(x)>0\) almost everywhere in Ω, \(k\in L^{\infty}(\Omega)\) is nonnegative, \(\alpha>1\) and \(\beta\in(0,p-1)\).

In the past few decades much attention has been devoted to nonlinear elliptic equations with singularities because of their wide applications in applied sciences, for example, non-Newtonian fluids, boundary layer phenomena for viscous fluids, chemical heterogenous catalysts, etc. (see, for example, [1, 2]). When \(p=2\) and the singularity is weak (i.e., \(0<\alpha<1\)), the existence, uniqueness and multiplicity of positive solutions have been established (see [38] and the references therein). By a result of Shi and Yao [5] one knows that if \(h\in C^{2,\gamma}(\overline{\Omega})\), \(h(x)\geq0\), then (1) with \(p=2\), \(0<\alpha<1\) and \(k(x)\equiv\lambda\) has one and only one solution \(u\in C^{2,\gamma }(\Omega)\cap C^{1,1-\alpha}(\overline{\Omega})\) for all \(\lambda\geq0\). However, when the singularity is strong (i.e., \(\alpha>1\)), various difficulties arise and the situation becomes more complicated.

A famous result obtained by Lazer and Mckenna (see [9]) showed that if \(h\in C^{\gamma}(\overline{\Omega})\), \(h(x)>0\) for all \(x\in \overline{\Omega}\), then (1) with \(p=2\), \(k(x)\equiv0\) and \(0<\alpha<\infty\) admits a unique solution \(u\in C^{2,\gamma}(\Omega)\cap C(\overline{\Omega})\). Furthermore, it was also shown in [9] that the solution u is not in \(C^{1}(\overline{\Omega})\) if \(\alpha>1\) and \(u\in H^{1}_{0}(\Omega)\) if and only if \(\alpha<3\). However, when \(h(x)\) does not have a positive lower bound on \(\overline{\Omega}\), whether or not there is an \(H_{0}^{1}(\Omega)\) solution of (1) may be determined by both \(h(x)\) and the parameter α. In particular, when \(h(x)\) behaves as \(\operatorname{dist}^{\gamma}(x,\partial\Omega)\) with \(\gamma\in\mathbb{R}\) (i.e., there exist \(c,C>0\) such that \(c \operatorname{dist}^{\gamma}(x,\partial \Omega)\leq h(x)\leq C\operatorname{dist}^{\gamma}(x,\partial\Omega)\) in Ω), (1) with \(p=2\) and \(k(x)\equiv0\) has a solution \(u\in H_{0}^{1}(\Omega)\) if and only if \(\alpha-2\gamma<3\). Therefore, it is reasonable to conjecture that (1) admits a solution in \(H_{0}^{1}(\Omega)\) if and only if the couple \((h(x),\alpha)\) satisfies some ‘compatibility condition’. A positive answer to this conjecture was given by Sun in [10], where she showed that (1) with \(p=2\) admits at least one solution in \(H_{0}^{1}(\Omega)\) if and only if there exists \(u_{0}\in H_{0}^{1}(\Omega)\) such that

$$ \int_{\Omega}h(x)|u_{0}|^{1-\alpha} \,\mathrm{d}x<\infty. $$
(CD)

Note that the necessity of (CD) is obvious since if (1) has an \(H_{0}^{1}(\Omega)\) solution u, then by choosing \(\varphi=u\) in the definition of weak solution of (1) (see the definition below) it follows that \(\int_{\Omega}h(x)|u|^{1-\alpha}\,\mathrm{d}x=\int_{\Omega}|\nabla u|^{2}\,\mathrm{d}x-\int_{\Omega}k(x)|u|^{1+\beta}\,\mathrm{d}x<\infty\). An additional significant paper is the paper by Crandall et al. [11], where the existence of solutions to the more general problem

$$Lu=g(x,u) \quad \mbox{in }\Omega, \qquad u=0 \quad\mbox{on }\partial\Omega $$

was studied. Here L is a linear second order elliptic operator which satisfies the maximum principle and g is positive and becomes singular uniformly in x as \(u\rightarrow0\). Their techniques are also based on the use of sub-supersolution theorems. Another paper worth mentioning is the one by Boccardo and Orsina [12], in which they proved for \(p=2\), \(k(x)\equiv0\) and general \(\alpha>0\) that (1) has only an \(H^{1}_{\mathrm{loc}}(\Omega)\)-solution when \(h(x)\in L^{1}(\Omega)\), \(h(x)\geq\not\equiv0\). However, whether or not the solutions are \(H_{0}^{1}(\Omega)\) functions was not answered there.

When \(p\neq2\), the corresponding results are much fewer. By using a sub-supersolution approach and a mountain pass theorem, Giacomoni et al. [13] proved, among other things, that when \(h(x)\equiv\lambda\), \(k(x)\equiv1\), \(0<\alpha<1\) and \(p-1<\beta <p^{*}-1\) (\(p^{*}\) is the critical Sobolev exponent of p), problem (1) has multiple weak solutions (depending on certain value of the parameter λ). Later, Loc and Schmitt [14] investigated the following singular problem:

$$ \begin{cases} -\Delta_{p} u=a(x)g(u)+\lambda h(u),&x\in\Omega,\\ u(x)=0,&x\in\partial\Omega,\\ u(x)>0,&x\in\Omega, \end{cases} $$
(2)

where \(p>1\), \(g(u)\) is a non-monotone singular term, \(a\in L^{\infty}(\Omega)\), \(a(x)\geq1\) for almost every \(x\in\Omega\), λ is a nonnegative parameter and \(h(u)\) is continuous. They first established a sub-supersolution theorem and then proved, under some structural conditions, that the problem had at least one weak solution in \(W^{1,p}_{0}(\Omega)\) by constructing ordered sub- and super-solutions.

It is our purpose to investigate problem (1) with a general positive function \(h\in L^{1}(\Omega)\). We shall show that the compatibility condition (CD) on the couple \((h(x),\alpha)\) is also optimal for the existence of weak solutions to problem (1). By weak solutions we mean that \(u\in W_{0}^{1,p}(\Omega)\) satisfying \(u(x)>0\) a.e. in Ω and

$$ \int_{\Omega}\bigl[|\nabla u|^{p-2}\nabla u\cdot\nabla\varphi -h(x)u^{-\alpha}\varphi-k(x)u^{\beta}\varphi \bigr]\,\mathrm{d}x=0,\quad \forall \varphi\in W_{0}^{1,p}( \Omega). $$
(3)

Since \(h(x)\) does not have a positive lower bound on \(\overline{\Omega}\) in general, the sub- and super-solution method used in [14] cannot be applied in our paper to prove the existence of weak solutions to (1). Instead, we will deal with (1) in variational framework and introduce the following singular energy functional (\(\alpha>1\)):

$$I(u)=\frac{1}{p}\int_{\Omega}|\nabla u|^{p}\, \mathrm{d}x+\frac{1}{\alpha -1}\int_{\Omega}h(x)|u|^{1-\alpha}\, \mathrm{d}x-\frac{1}{1+\beta}\int_{\Omega}k(x)|u|^{1+\beta} \,\mathrm{d}x. $$

As was pointed out in [10], the main difficulty that arises in such problems is the absence of integrability of \(u^{-\alpha}\) for \(u\in W^{1,p}_{0}(\Omega)\), and all the inequalities related to \(u\in W^{1,p}_{0}(\Omega)\) will not be of much help. Therefore, one will encounter many difficulties in trying to find \(W^{1,p}_{0}\)-solutions by standard variational methods. The first one is that there is a sharp contrast between the case \(0<\alpha<1\), for which the energy functional is continuous, and the case \(\alpha>1\), for which the energy functional is singular. In order to deal with this difficulty, we use variational techniques in the spirit of the work by Sun et al. [7] and consider the functional \(I(u)\) under appropriate constraints to restore integrability. Another difficulty arises from the nonlinearity of the p-Laplace operator, since generally one cannot deduce from \(u_{n}\rightharpoonup u\) in \(W_{0}^{1,p}(\Omega)\) the convergence \(|\nabla u_{n}|^{p-2}\nabla u_{n}\rightharpoonup |\nabla u|^{p-2}\nabla u\) in \(L^{\frac{p}{p-1}}(\Omega,\mathbb{R}^{N})\), which forces us to derive more properties of the functional \(I(u)\) to deal with the weak convergence. The main results of this paper will be stated and proved in Section 2.

2 Main results

Very little can be done about \(I(u)\) unless we confine ourselves to some subset of \(W_{0}^{1,p}(\Omega)\) because of the singularity of \(I(u)\), which results from that fact \(\alpha>1\). To overcome this difficulty, we first define the set of the specific choice of constraints

$$N_{1}= \biggl\{ u\in W_{0}^{1,p}(\Omega);\int _{\Omega}\bigl[|\nabla u|^{p}-h(x)|u|^{1-\alpha}-k(x)|u|^{1+\beta} \bigr]\,\mathrm{d}x\geq0 \biggr\} , $$

and

$$N_{2}= \biggl\{ u\in W_{0}^{1,p}(\Omega);\int _{\Omega}\bigl[|\nabla u|^{p}-h(x)|u|^{1-\alpha}-k(x)|u|^{1+\beta} \bigr]\,\mathrm{d}x=0 \biggr\} . $$

Typically, \(N_{2}\) is not a closed set anymore for \(\alpha>1\) (certainly not weakly closed). Moreover, we will equip \(W_{0}^{1,p}(\Omega)\) with the norm \(\|u\|_{1,p;\Omega}=(\int_{\Omega}|\nabla u|^{p} \,\mathrm{d}x)^{\frac{1}{p}}\) due to Poincare’s inequality. The subscript will be omitted when no confusion arises.

Under the above hypothesis, our main results of this paper read as follows.

Theorem 2.1

Let \(\Omega\subset\mathbb{R}^{N} \) (\(N\geq1\)) be a bounded domain with smooth boundary Ω, \(k\in L^{\infty}(\Omega)\) be a nonnegative function, \(h\in L^{1}(\Omega)\), \(h(x)>0\) a.e. in Ω (not necessarily with a positive lower bound and possibly unbounded in Ω), and let \(p>2\), \(\alpha>1\), \(\beta\in(0,p-1)\). Then problem (1) admits a unique \(W_{0}^{1,p}\)-solution if and only if there exists \(u_{0}\in W_{0}^{1,p}(\Omega)\) such that

$$ \int_{\Omega}h(x)|u_{0}|^{1-\alpha} \,\mathrm{d}x<+\infty. $$
(4)

Proof

Necessity. Suppose that \(u\in W_{0}^{1,p}(\Omega)\) is a solution of (1), by choosing \(\varphi=u\) in (3), we then have

$$\int_{\Omega}h(x)|u|^{1-\alpha}\,\mathrm{d}x=\int _{\Omega}\bigl[|\nabla u|^{p}-k(x)|u|^{1+\beta} \bigr]\,\mathrm{d}x<+\infty. $$

Sufficiency. The key to solving (1) relies on a natural interpolation between the two constraints \(N_{1}\) and \(N_{2}\). For the convenience of the reader, we divide this part into three steps.

Step 1. Properties of \(N_{1}\) and \(N_{2}\) .

Note that, for any \(u\in W_{0}^{1,p}(\Omega)\) with \(\int_{\Omega}h(x)|u|^{1-\alpha}\,\mathrm{d}x<\infty\), there exists a rectilinear curve \(t\rightarrow tu\) in the positive semi-axis through \(N_{2}\) at \(t=t(u)\) satisfying \(tu\in N_{1}\) for all \(t\geq t(u)\), which implies that \(N_{1}\) is unbounded. Therefore, we know from (4) that \(N_{2}\) and \(N_{1}\) are nonempty since \(t(u_{0})u_{0}\in N_{2}\subset N_{1}\). Moreover, along this curve,

$$ I(tu)\geq I\bigl(t(u)u\bigr),\quad \forall t>0. $$
(5)

The closedness of \(N_{1}\) follows easily from Fatou’s lemma. However, since \(\int_{\Omega}h(x)|u|^{1-\alpha}\,\mathrm{d}x\) is not continuous even on \(N_{2}\), \(N_{2}\) is not anymore a closed set in \(W_{0}^{1,p}(\Omega)\).

Finally, unbounded set \(N_{1}\) lies in the exterior of \(W_{0}^{1,p}(\Omega )\) (i.e., it stays away from a ball centered at 0). Suppose on the contrary that there is a sequence \((u_{n})\subset N_{1}\) with \(u_{n}\rightarrow0\) in \(W_{0}^{1,p}(\Omega)\). Then, by using the reversed Hölder’s inequality, it follows that

$$\begin{aligned} \biggl(\int_{\Omega}h^{\frac{1}{\alpha}}(x)\,\mathrm{d}x \biggr)^{\alpha}\biggl( \int_{\Omega}|u_{n}|\, \mathrm{d}x \biggr)^{1-\alpha} \leq&\int_{\Omega}h(x)|u_{n}|^{1-\alpha}\,\mathrm{d}x \\ \leq&\|u_{n}\|^{p}-\int_{\Omega}k(x)|u_{n}|^{1+\beta}\,\mathrm{d}x\rightarrow0. \end{aligned}$$

Consequently, \(\int_{\Omega}|u_{n}|\,\mathrm{d}x\rightarrow\infty\), since \(\alpha>1\), which is clearly impossible. Therefore, there is a positive constant \(c_{1}\) such that \(\|u\|\geq c_{1}\) for any \(u\in N_{1}\).

Step 2. Properties of the minimizing sequence \(\{u_{n}\}\) .

Now we consider \(\inf_{N_{1}}I\). By Ekeland’s variational principle [15], there exists a corresponding minimizing sequence \(\{u_{n}\}\in N_{1}\) satisfying:

  1. (i)

    \(I(u_{n})\leq \inf_{N_{1}}I+\frac{1}{n} \);

  2. (ii)

    \(I(u_{n})\leq I(\omega)+\frac{1}{n}\|u_{n}-\omega\|\), \(\forall \omega \in N_{1}\).

We may assume that \(u_{n}\geq0\) in Ω since \(I(|u_{n}|)=I(u)\). Furthermore, the fact \(h(x)>0\) a.e. in Ω, \(\alpha>1\) and the following inequality \(\int_{\Omega}h(x)|u_{n}|^{1-\alpha}\,\mathrm{d}x \leq\|u_{n}\|^{p}-\int_{\Omega}k(x)|u_{n}|^{1+\beta}\,\mathrm{d}x\) (since \(u_{n}\in N_{1}\)) guarantee that \(u_{n}(x)>0\) a.e. in Ω. Moreover, one can observe from \(\alpha>1\) and \(k\in L^{\infty}(\Omega)\) that

$$\begin{aligned} I(u) =&\frac{1}{p}\int_{\Omega}|\nabla u|^{p}\, \mathrm{d}x+\frac{1}{\alpha -1}\int_{\Omega}h(x)|u|^{1-\alpha}\, \mathrm{d}x-\frac{1}{1+\beta}\int_{\Omega}k(x)|u|^{1+\beta} \,\mathrm{d}x \\ \geq&\frac{1}{p}\|u\|^{p}-C\|u\|^{\beta+1}_{\beta+1} \end{aligned}$$

for some constant \(C>0\), which, due to \(\beta< p-1\), implies the coercivity of \(I(u)\) on \(N_{1}\) and then the existence of a constant \(c_{2}>0\) such that \(\|u_{n}\|\leq c_{2}\). Therefore, there exist a subsequence of \(\{u_{n}\}\), which is still denoted by \(\{u_{n}\}\), and \(u^{*}\in W_{0}^{1,p}(\Omega)\) such that \(u_{n}\rightharpoonup u^{*}\) weakly in \(W_{0}^{1,p}(\Omega)\), strongly in \(L^{p}(\Omega)\), and pointwise a.e. in Ω. We will try to show that \(u^{*}\in N_{2}\) by analyzing the minimizing sequence \(\{u_{n}\}\). First, by Fatou’s lemma we know that

$$\int_{\Omega}h(x)\bigl|u^{*}\bigr|^{1-\alpha}\,\mathrm{d}x\leq\liminf _{n\rightarrow\infty}\int_{\Omega}h(x)|u_{n}|^{1-\alpha} \,\mathrm{d}x<\infty, $$

which implies that \(u^{*}>0\) a.e. in Ω since \(\alpha>1\). Furthermore, we have

$$\begin{aligned} \inf_{N_{1}} I =&\lim_{n\rightarrow\infty}I(u_{n}) \\ =&\liminf_{n\rightarrow\infty} \biggl[\frac{1}{p}\int _{\Omega}|\nabla u_{n}|^{p}\,\mathrm{d}x+ \frac{1}{\alpha -1}\int_{\Omega}h(x)u_{n}^{1-\alpha} \,\mathrm{d}x-\frac{1}{1+\beta}\int_{\Omega}k(x)u_{n}^{1+\beta} \,\mathrm{d}x \biggr] \\ \geq&\liminf_{n\rightarrow\infty} \biggl[\frac{1}{p}\int _{\Omega}|\nabla u_{n}|^{p} \,\mathrm{d}x \biggr]+\liminf_{n\rightarrow \infty} \biggl[\frac{1}{\alpha-1}\int _{\Omega}h(x)u_{n}^{1-\alpha}\,\mathrm{d}x \biggr] \\ &{}-\frac{1}{1+\beta}\int_{\Omega}k(x)u^{*1+\beta}\,\mathrm{d}x \\ \geq&\frac{1}{p}\int_{\Omega}\bigl|\nabla u^{*}\bigr|^{p} \,\mathrm{d}x +\frac{1}{\alpha-1}\int_{\Omega}h(x)u^{*1-\alpha} \,\mathrm{d}x-\frac{1}{1+\beta} \int_{\Omega}k(x)u^{*1+\beta} \,\mathrm{d}x \\ =&I\bigl(u^{*}\bigr)\geq I\bigl(t\bigl(u^{*}\bigr)u^{*}\bigr) \\ \geq&\inf_{N_{2}}I\geq\inf_{N_{1}} I. \end{aligned}$$
(6)

Thus, the above inequalities are actually equalities. Particularly, it follows that

$$\liminf_{n\rightarrow\infty}\int_{\Omega}|\nabla u_{n}|^{p} \,\mathrm{d}x =\int_{\Omega}\bigl| \nabla u^{*}\bigr|^{p} \,\mathrm{d}x, $$

and there exists a subsequence of \(\{u_{n}\}\) (still denoted by \(\{u_{n}\} \)) such that

$$\lim_{n\rightarrow\infty}\int_{\Omega}|\nabla u_{n}|^{p} \,\mathrm{d}x =\int_{\Omega}\bigl| \nabla u^{*}\bigr|^{p} \,\mathrm{d}x. $$

This together with the weak convergence of \(u_{n}\rightharpoonup u^{*}\) in \(W_{0}^{1,p}(\Omega)\) implies \(u_{n}\rightarrow u^{*}\) strongly in \(W_{0}^{1,p}(\Omega)\). Consequently, we have

$$ |\nabla u_{n}|^{p-2}\nabla u_{n} \rightharpoonup\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*} \quad\mbox{in } L^{ \frac{p}{p-1}}\bigl(\Omega,\mathbb{R}^{N}\bigr). $$
(s)

We shall distinguish two cases according to \(u_{n}\) belonging to \(N_{1}\setminus N_{2}\) or \(N_{2}\) to show that \(u^{*}\in N_{2}\).

Case 1. Suppose that \(\{u_{n}\}\subset N_{1}\setminus N_{2}\) for all n large. Fix \(\varphi\in W_{0}^{1,p}(\Omega)\), \(\varphi\geq0\) and n by now. Since \(u_{n}\in N_{1}\setminus N_{2}\), \(h(x)\) is nonnegative in Ω and \(\alpha>1\), we have

$$\int_{\Omega}h(x) (u_{n}+t\varphi)^{1-\alpha}\, \mathrm{d}x\leq \int_{\Omega}h(x)u_{n}^{1-\alpha}\, \mathrm{d}x<\|u_{n}\|^{p}-\int_{\Omega}k(x)u_{n}^{1+\beta }\,\mathrm{d}x, \quad\forall t\geq0. $$

Subsequently, we can choose \(t>0\) sufficiently small such that

$$\int_{\Omega}h(x) (u_{n}+t\varphi)^{1-\gamma}\, \mathrm{d}x<\|u_{n}+t\varphi\|^{p}-\int_{\Omega}k(x) (u_{n}+t\varphi)^{1+\beta}\,\mathrm{d}x, $$

which implies that \(u_{n}+t\varphi\in N_{1}\). Furthermore, by choosing \(\omega=u_{n}+t\varphi\) in (ii), we obtain

$$\begin{aligned} &\frac{1}{n}\|t\varphi\|+\frac{1}{p}\bigl(\|u_{n}+t \varphi\|^{p}-\|u_{n}\| ^{p}\bigr)- \frac{1}{1+\beta}\int_{\Omega}k(x)\bigl[(u_{n}+t \varphi)^{1+\beta }-u_{n}^{1+\beta}\bigr]\,\mathrm{d}x \\ &\quad\geq\frac{1}{1-\alpha}\int_{\Omega}h(x) \bigl[(u_{n}+t\varphi)^{1-\alpha }-u_{n}^{1-\gamma} \bigr]\,\mathrm{d}x. \end{aligned}$$

Dividing both sides of the above inequality by \(t>0\), passing to the liminf as \(t\rightarrow0\) and using Fatou’s lemma, we have

$$\begin{aligned} &\frac{\|\varphi\|}{n}+\int_{\Omega}|\nabla u_{n}|^{p-2} \nabla u_{n}\cdot \nabla\varphi\,\mathrm{d}x-\int_{\Omega}k(x)u_{n}^{\beta}\varphi\,\mathrm{d}x \\ &\quad\geq\int_{\Omega}\liminf_{t\rightarrow0} \frac{h(x)[(u_{n}+t\varphi )^{1-\alpha}-u_{n}^{1-\alpha}]}{t(1-\alpha)}\,\mathrm{d}x \\ &\quad=\int_{\Omega}h(x)u_{n}^{-\alpha}\varphi \,\mathrm{d}x. \end{aligned}$$

By letting \(n\rightarrow\infty\), using Fatou’s lemma again and noticing (s), we arrive at

$$ \int_{\Omega}\bigl[\bigl|\nabla u^{*}\bigr|^{p-2} \nabla u^{*}\cdot\nabla\varphi -h(x)u^{*-\alpha}\varphi \bigr] \,\mathrm{d}x\geq \int_{\Omega}k(x)u^{*\beta }\varphi\,\mathrm{d}x, \quad \forall \varphi\in W_{0}^{1,p}(\Omega) , \varphi\geq0. $$
(7)

Choosing \(\varphi=u^{*}\) in (7) yields

$$ \int_{\Omega}\bigl[\bigl|\nabla u^{*}\bigr|^{p}-h(x)u^{*1-\alpha} \bigr]\,\mathrm{d}x\geq \int_{\Omega}k(x)u^{*1+\beta}\, \mathrm{d}x, $$
(8)

which implies that \(u^{*}\in N_{1}\). Furthermore, we know from (6) that

$$ u^{*}\in N_{2} \quad\bigl(\textit{i.e. } t\bigl(u^{*} \bigr)=1\bigr). $$
(9)

Case 2. There exists a subsequence of \(\{u_{n}\}\), still denoted by \(\{ u_{n}\}\), that belongs to \(N_{2}\). Let \(\varphi\in W_{0}^{1,p}(\Omega)\), \(\varphi\geq0\) and n be fixed again. Then, for all \(t\geq0\), \(\int_{\Omega}h(x)(u_{n}+t\varphi)^{1-\alpha}\,\mathrm{d}x \leq\int_{\Omega}h(x)u_{n}^{1-\alpha}\,\mathrm{d}x=\|u_{n}\|^{p}<+\infty\) since \(\alpha >1\), which ensures the existence of a unique positive number corresponding to \(u_{n}+t\varphi\), denoted by \(f_{n,\varphi}(t)\), such that \(f_{n,\varphi}(t)(u_{n}+t\varphi)\in N_{2}\), that is, \(f_{n,\varphi}(t)\) satisfies

$$f^{p-1-\beta}_{n,\varphi}(t)\|u_{n}+t\varphi \|^{p}-f^{-\alpha-\beta }_{n,\varphi}(t)\int_{\Omega}h(x) (u_{n}+t\varphi)^{1-\alpha}\,\mathrm{d}x =\int _{\Omega}k(x) (u_{n}+t\varphi)^{1+\beta}\, \mathrm{d}x. $$

It is easy to show that \(f_{n,\varphi}(t)\) is continuous for all \(t\geq 0\) by the dominated convergence theorem, and \(f_{n,\varphi}(0)=1\) since \(u_{n}\in N_{2}\). Next, we shall show that \(f_{n,\varphi}(t)\) has uniform behavior at 0 with respect to n, i.e., \(|f_{n,\varphi}'(0)|\leq C\) for suitable \(C>0\) independent of n. For this purpose, we assume henceforth that

$$f_{n,\varphi}'(0)=\lim_{t\rightarrow0} \frac{f_{n,\varphi}(t)-1}{t}\in [-\infty,+\infty] . $$

If the above limit does not exist, we can choose a positive sequence \(\{ t_{k}\}\) with \(t_{k}\rightarrow0\) (instead of \(t\rightarrow0\)) in such a way that \(q_{n} :=\lim_{k\rightarrow\infty}\frac{f_{n,\varphi }(t_{k})-1}{n} \in[-\infty, +\infty]\), and replace \(f_{n,\varphi}'(0)\) by \(q_{n}\). We first derive the upper bound of \(f_{n,\varphi}'(0)\). From \(f_{n,\varphi}(t)(u_{n}+t\varphi)\), \(u_{n} \in N_{2}\), we have

$$\begin{aligned}& 0=f^{p}_{n,\varphi}(t)\|u_{n}+t \varphi\|^{p}-f^{1-\alpha}_{n,\varphi}(t)\int_{\Omega}h(x) (u_{n}+t\varphi)^{1-\alpha}\,\mathrm{d}x -f^{1+\beta}_{n,\varphi}(t)\int_{\Omega}k(x) (u_{n}+t\varphi )^{1+\beta}\,\mathrm{d}x, \\& 0=\|u_{n}\|^{p}-\int_{\Omega}h(x)u_{n}^{1-\alpha}\,\mathrm{d}x-\int_{\Omega}k(x)u_{n}^{1+\beta}\,\mathrm{d}x. \end{aligned}$$

Combining the above two equalities and noticing \(\alpha>1\), we obtain

$$\begin{aligned} 0 \geq& \biggl\{ p\bigl[f_{n,\varphi}(0)+o(1)\bigr]^{p-1} \|u_{n}+t\varphi\|^{p} \\ &{}-(1-\alpha)\bigl[f_{n,\varphi}(0)+o(1)\bigr]^{-\alpha}\int _{\Omega}h(x) (u_{n}+t\varphi)^{1-\alpha}\, \mathrm{d}x \\ &{}-(1+\beta)\bigl[f_{n,\varphi}(0)+o(1)\bigr]^{\beta}\int _{\Omega}k(x) (u_{n}+t\varphi )^{1+\beta}\, \mathrm{d}x \biggr\} \bigl[f_{n,\varphi}(t)-f_{n,\varphi}(0) \bigr] \\ &{}+\|u_{n}+t\varphi\|^{p}-\|u_{n} \|^{p}-\int_{\Omega}k(x)\bigl[(u_{n}+t \varphi)^{1+\beta }-u_{n}^{1+\beta}\bigr]\,\mathrm{d}x. \end{aligned}$$

Dividing by \(t>0\) on both sides of the above inequalities, letting \(t\rightarrow0\) and noticing the continuity of \(f_{n,\varphi}(t)\), we arrive at

$$\begin{aligned} 0 \geq& \biggl[p\|u_{n}\|^{p}+(\alpha-1)\int _{\Omega}h(x)u_{n}^{1-\alpha}\,\mathrm{d}x -(1+ \beta)\int_{\Omega}k(x)u_{n}^{1+\beta}\, \mathrm{d}x \biggr] f_{n,\varphi}'(0) \\ &{}+p\int_{\Omega}|\nabla u_{n}|^{p-2}\nabla u_{n}\cdot\nabla\varphi\,\mathrm{d}x -(1+\beta)\int _{\Omega}k(x)u_{n}^{\beta}\varphi\,\mathrm{d}x. \end{aligned}$$

Thus, from \(\{u_{n}\}\subset N_{2}\) we know that

$$\begin{aligned} 0 \geq& \biggl[(\alpha-1)\int_{\Omega}h(x)u_{n}^{1-\alpha} \,\mathrm{d}x+(p-1-\beta) \int_{\Omega}k(x)u_{n}^{1+\beta} \,\mathrm{d}x \biggr]f_{n,\varphi}'(0) \\ &{}+p\int_{\Omega}|\nabla u_{n}|^{p-2}\nabla u_{n}\cdot\nabla\varphi\,\mathrm{d}x -(1+\beta)\int _{\Omega}k(x)u_{n}^{\beta}\varphi\,\mathrm{d}x. \end{aligned}$$

Since \(\{u_{n}\}\subset N_{2}\) (\(\subset N_{1}\)) is bounded and \(N_{1}\) stays away from some ball centered at 0, the above inequality necessarily implies that

$$ f_{n,\varphi}'(0)\in[-\infty, +\infty) \quad \mbox{and}\quad f_{n,\varphi}'(0)\leq c_{3} \quad \mbox{uniformly for all large } n $$
(10)

for some suitable constant \(c_{3}>0\).

Next we shall derive the lower bound for \(f_{n,\varphi}'(0)\). By applying (ii) again, one gets

$$\begin{aligned} &\frac{1}{n}\bigl|f_{n,\varphi}(t)-1\bigr|\|u_{n}\|+\frac{1}{n}t f_{n,\varphi}(t) \|\varphi\| \\ &\quad\geq\frac{1}{n}\bigl\| u_{n}-f_{n,\varphi}(t) (u_{n}+t\varphi)\bigr\| \\ &\quad\geq I(u_{n})-I\bigl(f_{n,\varphi}(t) (u_{n}+t \varphi)\bigr), \end{aligned}$$

which is equivalent to

$$\begin{aligned} &\frac{1}{n}t f_{n,\varphi}(t)\|\varphi\|+\biggl(\frac{1}{p}+ \frac{1}{\alpha-1}\biggr) \bigl[\|u_{n}+t\varphi\|^{p}- \|u_{n}\|^{p} \bigr] \\ &\quad\geq \biggl\{ -\biggl(1+\frac{p}{\alpha-1}\biggr)\bigl[f_{n,\varphi}(0)+o(1) \bigr]^{p-1}\| u_{n}+t\varphi\|^{p} \\ &\qquad{}+\biggl(1+\frac{1+\beta}{\alpha-1}\biggr)\bigl[f_{n,\varphi}(0)+o(1) \bigr]^{\beta}\int_{\Omega}k(x) (u_{n}+t \varphi)^{1+\beta}\,\mathrm{d}x \\ &\qquad{}-\frac{\|u_{n}\|}{n}\operatorname{sgn} \bigl[f_{n,\varphi}(t)-1 \bigr] \biggr\} \bigl[f_{n,\varphi}(t)-1 \bigr] \\ &\qquad{}+\biggl(\frac{1}{1+\beta}+\frac{1}{\alpha_{1}}\biggr)\int _{\Omega}k(x)\bigl[(u_{n}+t\varphi )^{1+\beta}-u^{1+\beta}_{n}\bigr]\,\mathrm{d}x. \end{aligned}$$

Noticing

$$\begin{aligned} &{-}\biggl(1+\frac{p}{\alpha-1}\biggr)\bigl[f_{n,\varphi}(0)+o(1) \bigr]^{p-1}\|u_{n}+t\varphi\| ^{p} \\ &\qquad{}+\biggl(1+\frac{1+\beta}{\alpha-1}\biggr)\bigl[f_{n,\varphi}(0)+o(1) \bigr]^{\beta}\int_{\Omega}k(x) (u_{n}+t \varphi)^{1+\beta}\,\mathrm{d}x \\ &\quad\rightarrow-\biggl(1+\frac{p}{\alpha-1}\biggr)\|u_{n} \|^{p}+\biggl(1+\frac{1+\beta}{\alpha -1}\biggr)\int_{\Omega}k(x)u_{n}^{1+\beta}\,\mathrm{d}x \quad(t\rightarrow0) \\ &\quad=-\frac{p-1-\beta}{\alpha-1}\|u_{n}\|^{p}-\biggl(1+ \frac{1+\beta}{\alpha-1}\biggr)\int_{\Omega}h(x)u_{n}^{1-\alpha} \,\mathrm{d}x \quad(\mbox{since } u_{n}\in N_{2}) \\ &\quad\leq-\frac{p-1-\beta}{\alpha-1}\|u_{n}\|^{p}\leq- \frac{p-1-\beta}{\alpha-1}c_{1}^{p}, \end{aligned}$$

we see that \(f_{n,\varphi}'(0)\) cannot diverge to −∞ as \(n\rightarrow\infty\), that is,

$$ f_{n,\varphi}'(0)\in(-\infty, +\infty)\quad \mbox{and}\quad f_{n,\varphi}'(0)\geq c_{4}\quad \mbox{uniformly in all large } n $$
(11)

for some suitable real number \(c_{4}\).

Putting together (10) and (11), we find that

$$ f_{n,\varphi}'(0)\in(-\infty, +\infty)\quad \mbox{and}\quad \bigl|f_{n,\varphi}'(0)\bigr|\leq C \quad\mbox{uniformly in all large } n. $$
(12)

Based on (12), we are now in the position to locate \(u^{*}\) in Case 2. Let \(\varphi\in W_{0}^{1,p}(\Omega)\), \(\varphi\geq0\) and n be fixed. By (ii) we have

$$\begin{aligned} &\frac{\|u_{n}\|}{n}\bigl|f_{n,\varphi}(t)-1\bigr|+\frac{\|\varphi\| }{n}f_{n,\varphi}(t)t \\ &\quad\geq\frac{1}{n}\bigl\| u_{n}-f_{n,\varphi}(t) (u_{n}+t\varphi)\bigr\| \\ &\quad\geq I(u_{n})-I\bigl(f_{n,\varphi}(t) (u_{n}+t \varphi)\bigr) \\ &\quad= \biggl\{ -\bigl[f_{n,\varphi}(0)+o(1)\bigr]^{p-1} \|u_{n}+t\varphi\|^{p}+\bigl[f_{n,\varphi }(0)+o(1) \bigr]^{-\alpha}\int_{\Omega}h(x) (u_{n}+t \varphi)^{1-\alpha}\,\mathrm{d}x \\ &\qquad{}+\bigl[f_{n,\varphi}(0)+o(1)\bigr]^{\beta}\int _{\Omega}k(x) (u_{n}+t\varphi )^{1+\beta}\, \mathrm{d}x \biggr\} \bigl(f_{n,\varphi}(t)-1\bigr) \\ &\qquad{}-\frac{1}{p}\bigl(\|u_{n}+t\varphi\|^{p}- \|u_{n}\|^{p}\bigr)+\frac{1}{1-\alpha}\int _{\Omega}h(x)\bigl[(u_{n}+t\varphi)^{1-\alpha}-u_{n}^{1-\alpha} \bigr]\,\mathrm{d}x \\ &\qquad{}+\frac{1}{1+\beta}\int_{\Omega}k(x) \bigl[(u_{n}+t\varphi)^{1+\beta }-u_{n}^{1+\beta} \bigr]\,\mathrm{d}x. \end{aligned}$$

Thus, dividing by \(t>0\) on both sides of the inequality, passing to the liminf as \(t\rightarrow0\) and noticing that \(u_{n}\in N_{2}\), we get

$$\begin{aligned} &\frac{1}{n} \bigl[\bigl|f_{n,\varphi}'(0)\bigr| \|u_{n}\|+\|\varphi\| \bigr] \\ &\quad\geq \biggl[-\|u_{n}\|^{p}+\int_{\Omega}h(x)u_{n}^{1-\alpha}\,\mathrm{d}x+\int_{\Omega}k(x)u_{n}^{1+\beta}\,\mathrm{d}x \biggr]f_{n,\varphi}'(0) \\ &\qquad{}-\int_{\Omega}|\nabla u_{n}|^{p-2} \nabla u_{n}\cdot\nabla\varphi\,\mathrm{d}x +\int_{\Omega}h(x)u_{n}^{-\alpha}\varphi\,\mathrm{d}x+\int_{\Omega}k(x)u_{n}^{\beta}\varphi\,\mathrm{d}x \\ &\quad=-\int_{\Omega}|\nabla u_{n}|^{p-2} \nabla u_{n}\cdot\nabla\varphi\,\mathrm{d}x +\int_{\Omega}h(x)u_{n}^{-\alpha}\varphi\,\mathrm{d}x+\int_{\Omega}k(x)u_{n}^{\beta}\varphi\,\mathrm{d}x. \end{aligned}$$
(13)

Since \(|f_{n,\varphi}'(0)|\leq C\) uniformly in all n large, we know that \(h(x)u_{n}^{-\alpha}\varphi\) is integrable in Ω. Furthermore, letting \(n\rightarrow\infty\) in (13) and using Fatou’s lemma again, we obtain

$$\begin{aligned} &\int_{\Omega}\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*} \cdot\nabla\varphi\,\mathrm{d}x\geq \int_{\Omega}h(x)u^{*-\alpha}\varphi\,\mathrm{d}x+\int_{\Omega}k(x)u^{*\beta }\varphi\,\mathrm{d}x, \\ &\quad\forall \varphi\in W_{0}^{1,p}(\Omega), \varphi \geq0. \end{aligned}$$
(14)

Taking \(\varphi=u^{*}\) in (14), we get \(u^{*}\in N_{1}\). By the same argument as in Case 1, we see also

$$ u^{*}\in N_{2}. $$
(15)

Thus, from (7), (9), (14) and (15), we conclude that in either case \(u^{*}\in N_{2}\) and

$$\begin{aligned} &\int_{\Omega}\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*} \cdot\nabla\varphi\,\mathrm{d}x\geq \int_{\Omega}h(x)u^{*-\alpha}\varphi\,\mathrm{d}x+\int_{\Omega}k(x)u^{*\beta }\varphi\,\mathrm{d}x, \\ &\quad\forall \varphi\in W_{0}^{1,p}(\Omega), \varphi \geq0. \end{aligned}$$
(16)

The strict positivity \(u^{*}(x)>0\), \(\forall x\in\Omega\) follows from the maximum principle [16].

Step 3. Existence of weak solutions.

Finally, we will prove that \(u^{*}\in W_{0}^{1,p}(\Omega)\) is a solution to (1). To show this, for arbitrary \(\phi\in W_{0}^{1,p}(\Omega)\) and \(\varepsilon >0\), set \(\Omega^{+}_{\varepsilon}=\{x\in\Omega: u^{*}(x)+\varepsilon\phi (x)\geq0\}\) and \(\Omega^{-}_{\varepsilon}=\{x\in\Omega: u^{*}(x)+\varepsilon\phi(x)<0\} \). Taking \(\varphi(x)=(u^{*}(x)+\varepsilon\phi(x))^{+}\in W_{0}^{1,p}(\Omega)\) in (16) to obtain

$$\begin{aligned} 0 \leq&\int_{\Omega}\bigl[\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*} \cdot\nabla\varphi \,\mathrm{d}x-h(x)u^{*-\alpha}\varphi-k(x)u^{*\beta} \varphi \bigr]\,\mathrm{d}x \\ =&\int_{\Omega^{+}_{\varepsilon}} \bigl[\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*}\cdot \nabla\bigl(u^{*}+\varepsilon\phi\bigr) -h(x)u^{*-\alpha}\bigl(u^{*}+\varepsilon \phi\bigr)-k(x)u^{*\beta}\bigl(u^{*}+\varepsilon \phi\bigr) \bigr]\, \mathrm{d}x \\ =&\int_{\Omega}\bigl[\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*}\cdot \nabla \bigl(u^{*}+\varepsilon\phi\bigr) -h(x)u^{*-\alpha}\bigl(u^{*}+\varepsilon \phi\bigr)-k(x)u^{*\beta}\bigl(u^{*}+\varepsilon \phi\bigr) \bigr]\, \mathrm{d}x \\ &{}-\int_{\Omega^{-}_{\varepsilon}} \bigl[\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*} \cdot \nabla\bigl(u^{*}+\varepsilon\phi\bigr) -h(x)u^{*-\alpha}\bigl(u^{*}+ \varepsilon\phi\bigr)-k(x)u^{*\beta}\bigl(u^{*}+\varepsilon \phi\bigr) \bigr] \,\mathrm{d}x \\ =&\bigl\| u^{*}\bigr\| ^{p}-\int_{\Omega}h(x)u^{*1-\alpha}\, \mathrm{d}x-\int_{\Omega}k(x)u^{*1+\beta}\,\mathrm{d}x \\ &{}+\varepsilon\int_{\Omega}\bigl[\bigl|\nabla u^{*}\bigr|^{p-2} \nabla u^{*}\nabla\phi- h(x)u^{*-\alpha}\phi-k(x)u^{*\beta}\phi \bigr]\, \mathrm{d}x \\ &{}-\int_{\Omega^{-}_{\varepsilon}} \bigl[\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*} \nabla \bigl(u^{*}+\varepsilon\phi\bigr)-h(x)u^{*-\alpha}\bigl(u^{*}+\varepsilon \phi \bigr)-k(x)u^{*\beta}\bigl(u^{*}+\varepsilon\phi\bigr) \bigr]\, \mathrm{d}x \\ =&\varepsilon\int_{\Omega}\bigl[\bigl|\nabla u^{*}\bigr|^{p-2} \nabla u^{*}\cdot\nabla \phi- h(x)u^{*-\alpha}\phi-k(x)u^{*\beta}\phi \bigr]\,\mathrm{d}x \\ &{}-\int_{\Omega^{-}_{\varepsilon}} \bigl[\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*} \nabla \bigl(u^{*}+\varepsilon\phi\bigr)-h(x)u^{*-\alpha}\bigl(u^{*}+\varepsilon \phi \bigr)-k(x)u^{*\beta}\bigl(u^{*}+\varepsilon\phi\bigr) \bigr]\, \mathrm{d}x \\ \leq&\varepsilon\int_{\Omega}\bigl[\bigl|\nabla u^{*}\bigr|^{p-2} \nabla u^{*}\nabla\phi- h(x)u^{*-\alpha}\phi-k(x)u^{*\beta}\phi \bigr]\, \mathrm{d}x-\varepsilon\int_{\Omega ^{-}_{\varepsilon}}\bigl|\nabla u^{*}\bigr|^{p-2} \nabla u^{*}\nabla\phi\,\mathrm{d}x. \end{aligned}$$

Since the measure of \(\Omega^{-}_{\varepsilon}\) tends to zero as \(\varepsilon\rightarrow0\), dividing the above expression by \(\varepsilon>0\) and letting \(\varepsilon\rightarrow0\), we obtain

$$\int_{\Omega}\bigl[\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*}\nabla \phi- h(x)u^{*-\alpha}\phi-k(x)u^{*\beta}\phi \bigr]\,\mathrm{d}x\geq0, \quad \forall \phi\in W_{0}^{1,p}(\Omega). $$

Replacing ϕ by −ϕ we obtain the reserved inequality, and hence \(u^{*}\) actually satisfies the following:

$$\int_{\Omega}\bigl[\bigl|\nabla u^{*}\bigr|^{p-2}\nabla u^{*}\nabla \phi- h(x)u^{*-\alpha}\phi-k(x)u^{*\beta}\phi \bigr]\,\mathrm{d}x=0, \quad \forall \phi\in W_{0}^{1,p}(\Omega), $$

which means that \(u^{*}\) is a \(W_{0}^{1,p}\)-solution to (1). The proof is complete. □

References

  1. Callegari, A, Nachman, A: Some singular nonlinear differential equations arising in boundary layer theory. J. Math. Anal. Appl. 64, 96-105 (1978)

    Article  MATH  MathSciNet  Google Scholar 

  2. Callegari, A, Nachman, A: A nonlinear singular boundary value problem in the theory of pseudoplastic fluids. SIAM J. Appl. Math. 38, 275-281 (1980)

    Article  MATH  MathSciNet  Google Scholar 

  3. Diaz, JI, Hernández, J, Rakotoson, JM: On very weak positive solutions to some semilinear elliptic problems with simultaneous singular nonlinear and spatial dependence terms. Milan J. Math. 79, 233-245 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  4. Hernández, J, Mancebo, FJ, Vega, JM: Positive solutions for singular nonlinear elliptic equations. Proc. R. Soc. Edinb. A 137, 41-62 (2007)

    Article  MATH  MathSciNet  Google Scholar 

  5. Shi, JP, Yao, MX: On a singular semilinear elliptic problem. Proc. R. Soc. Edinb. A 128, 1389-1491 (1998)

    Article  MATH  MathSciNet  Google Scholar 

  6. Sun, YJ, Wu, SP: An exact estimate result for a class of singular equations with critical exponents. J. Funct. Anal. 260, 1257-1284 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  7. Sun, YJ, Wu, SP, Long, YM: Combined effects of singular and superlinear nonlinearities in some singular boundary value problems. J. Differ. Equ. 176, 511-531 (2001)

    Article  MATH  MathSciNet  Google Scholar 

  8. Yang, HT: Multiplicity and asymptotic behavior of positive solutions for a singular semilinear elliptic problem. J. Differ. Equ. 189, 487-512 (2003)

    Article  MATH  MathSciNet  Google Scholar 

  9. Lazer, AC, McKenna, PJ: On a singular nonlinear elliptic boundary value problems. Proc. Am. Math. Soc. 111, 720-730 (1991)

    Article  MathSciNet  MATH  Google Scholar 

  10. Sun, YJ: Compatibility phenomena in singular problems. Proc. R. Soc. Edinb. A 143, 1321-1430 (2013)

    Article  MATH  MathSciNet  Google Scholar 

  11. Crandall, MG, Rabinowitz, PH, Tartar, L: On a Dirichlet problem with a singular nonlinearity. Commun. Partial Differ. Equ. 2, 193-222 (1977)

    Article  MATH  MathSciNet  Google Scholar 

  12. Boccardo, L, Orsina, L: Semilinear elliptic equations with singular nonlinearities. Calc. Var. Partial Differ. Equ. 37, 363-380 (2010)

    Article  MATH  MathSciNet  Google Scholar 

  13. Giacomoni, J, Schindler, I, Takác̆, P: Sobolev versus Hölder local minimizers and existence of multiple solutions for a singular quasilinear equation. Ann. Sc. Norm. Super. Pisa, Cl. Sci. 6, 117-158 (2007)

    MATH  MathSciNet  Google Scholar 

  14. Loc, NH, Schmitt, K: Boundary value problems for singular elliptic equations. Rocky Mt. J. Math. 41, 555-572 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  15. Aubin, JP, Ekeland, I: Applied Nonlinear Analysis. Pure and Applied Mathematics. Wiley, New York (1984)

    MATH  Google Scholar 

  16. Vazquez, JL: A strong maximum principle for some quasilinear elliptic equations. Appl. Math. Optim. 12, 191-202 (1984)

    Article  MATH  MathSciNet  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the referees for their valuable comments and suggestions which improved the original manuscript. The authors also express their sincere gratitude to Professor Wenjie Gao for his enthusiastic guidance and constant encouragement. The second author is supported by NSFC (11271154, 11401252) and by Fundamental Research Funds of Jilin University (450060501179).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yuzhu Han.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to the manuscript and read and approved the final manuscript.

Rights and permissions

Open Access This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cong, S., Han, Y. Compatibility conditions for the existence of weak solutions to a singular elliptic equation. Bound Value Probl 2015, 27 (2015). https://doi.org/10.1186/s13661-015-0285-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0285-9

Keywords