Citation: Yilei Tang, Dongmei Xiao, Weinian Zhang, Di Zhu. Dynamics of epidemic models with asymptomatic infection and seasonal succession[J]. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1407-1424. doi: 10.3934/mbe.2017073
[1] | Lulu Tao, Rui He, Sihua Liang, Rui Niu . Existence and multiplicity of solutions for critical Choquard-Kirchhoff type equations with variable growth. AIMS Mathematics, 2023, 8(2): 3026-3048. doi: 10.3934/math.2023156 |
[2] | Zehra Yucedag . Variational approach for a Steklov problem involving nonstandard growth conditions. AIMS Mathematics, 2023, 8(3): 5352-5368. doi: 10.3934/math.2023269 |
[3] | Yun-Ho Kim . Multiple solutions to Kirchhoff-Schrödinger equations involving the p(⋅)-Laplace-type operator. AIMS Mathematics, 2023, 8(4): 9461-9482. doi: 10.3934/math.2023477 |
[4] | Deke Wu, Hongmin Suo, Linyan Peng, Guaiqi Tian, Changmu Chu . Existence and multiplicity of positive solutions for a class of Kirchhoff type problems with singularity and critical exponents. AIMS Mathematics, 2022, 7(5): 7909-7935. doi: 10.3934/math.2022443 |
[5] | Wangjin Yao . Variational approach to non-instantaneous impulsive differential equations with p-Laplacian operator. AIMS Mathematics, 2022, 7(9): 17269-17285. doi: 10.3934/math.2022951 |
[6] | Zhenluo Lou, Jian Zhang . On general Kirchhoff type equations with steep potential well and critical growth in R2. AIMS Mathematics, 2024, 9(8): 21433-21454. doi: 10.3934/math.20241041 |
[7] | Fugeng Zeng, Peng Shi, Min Jiang . Global existence and finite time blow-up for a class of fractional p-Laplacian Kirchhoff type equations with logarithmic nonlinearity. AIMS Mathematics, 2021, 6(3): 2559-2578. doi: 10.3934/math.2021155 |
[8] | Najla Alghamdi, Abdeljabbar Ghanmi . Multiple solutions for a singular fractional Kirchhoff problem with variable exponents. AIMS Mathematics, 2025, 10(1): 826-838. doi: 10.3934/math.2025039 |
[9] | Liuyang Shao, Haibo Chen, Yicheng Pang, Yingmin Wang . Multiplicity of nontrivial solutions for a class of fractional Kirchhoff equations. AIMS Mathematics, 2024, 9(2): 4135-4160. doi: 10.3934/math.2024203 |
[10] | Nanbo Chen, Honghong Liang, Xiaochun Liu . On Kirchhoff type problems with singular nonlinearity in closed manifolds. AIMS Mathematics, 2024, 9(8): 21397-21413. doi: 10.3934/math.20241039 |
We study the following Neumann problem of Kirchhoff type equation with critical growth
{−(a+b∫Ω|∇u|2dx)Δu+u=Q(x)|u|4u+λP(x)|u|q−2u,in Ω,∂u∂v=0,on ∂Ω, | (1.1) |
where Ω ⊂ R3 is a bounded domain with a smooth boundary, a,b>0, 1<q<2, λ>0 is a real parameter. We assume that Q(x) and P(x) satisfy the following conditions:
(Q1) Q(x)∈C(ˉΩ) is a sign-changing;
(Q2) there exists xM∈Ω such that QM=Q(xM)>0 and
|Q(x)−QM|=o(|x−xM|)asx→xM; |
(Q3) there exists 0∈∂Ω such that Qm=Q(0)>0 and
|Q(x)−Qm|=o(|x|)asx→0; |
(P1) P(x) is positive continuous on ˉΩ and P(x0)=maxx∈ˉΩP(x);
(P2) there exist σ>0, R>0 and 3−q<β<6−q2 such that P(x)≥σ|x−y|−β for |x−y|≤R, where y is xM∈Ω or 0∈∂Ω.
In recent years, the following Dirichlet problem of Kirchhoff type equation has been studied extensively by many researchers
{−(a+b∫Ω|∇u|2dx)Δu=f(x,u),in Ω,u=0,on ∂Ω, | (1.2) |
which is related to the stationary analogue of the equation
utt−(a+b∫Ω|∇u|2dx)Δu=f(x,u) | (1.3) |
proposed by Kirchhoff in [13] as an extension of the classical D'Alembert's wave equation for free vibrations of elastic strings. Kirchhoff's model takes into account the changes in length of the string produced by transverse vibrations. In (1.2) and (1.3), u denotes the displacement, b is the initial tension and f(x,u) stands for the external force, while a is related to the intrinsic properties of the string (such as Young's modulus). We have to point out that such nonlocal problems appear in other fields like biological systems, such as population density, where u describes a process which depends on the average of itself (see Alves et al. [2]). After the pioneer work of Lions [18], where a functional analysis approach was proposed. The Kirchhoff type Eq (1.2) with critical growth began to call attention of researchers, we can see [1,9,14,17,23,24,28,30] and so on.
Recently, the following Kirchhoff type equation has been well studied by various authors
{−(a+b∫R3|∇u|2dx)Δu+V(x)u=f(x,u),inR3,u>0,u∈H1(R3). | (1.4) |
There has been much research regarding the concentration behavior of the positive solutions of (1.4), we can see [10,11,12,25,33]. Many papers studied the existence of ground state solutions of (1.4), for example [5,8,15,16,21,22,24]. In addition, the authors established the existence of sign-changing solutions of (1.4) in [20,31]. In papers [27,32] proved the existence and multiplicity of nontrivial solutions of (1.4) by using mountain pass theorem.
In particular, Chabrowski in [6] studied the solvability of the Neumann problem
{−Δu=Q(x)|u|2∗−2u+λf(x,u),in Ω,∂u∂v=0,on ∂Ω, |
where Ω ⊂ RN is a smooth bounded domain, 2∗=2NN−2(N≥3) is the critical Sobolev exponent, λ>0 is a parameter. Assume that Q(x)∈C(¯Ω) is a sign-changing function and ∫ΩQ(x)dx<0, under the condition of f(x,u). Using the space decomposition H1(Ω)=span1⊕V, where V={v∈H1(Ω):∫Ωvdx=0}, the author obtained the existence of two distinct solutions by the variational method.
In [14], Lei et al. considered the following Kirchhoff type equation with critical exponent
{−(a+b∫Ω|∇u|2dx)Δu=u5+λuq−1|x|β,in Ω,u=0,on ∂Ω, |
where Ω ⊂ R3 is a smooth bounded domain, a,b>0, 1<q<2, λ>0 is a parameter. They obtained the existence of a positive ground state solution for 0≤β<2 and two positive solutions for 3−q≤β<2 by the Nehari manifold method.
In [34], Zhang obtained the existence and multiplicity of nontrivial solutions of the following equation
{−(a+b∫Ω|∇u|2dx)Δu+u=λ|u|q−2u+f(x,u)+Q(x)u5,in Ω,∂u∂v=0,on ∂Ω, | (1.5) |
where Ω is an open bounded domain in R3, a,b>0, 1<q<2, λ≥0 is a parameter, f(x,u) and Q(x) are positive continuous functions satisfying some additional assumptions. Moreover, f(x,u)∼|u|p−2u with 4<p<6.
Comparing with the above mentioned papers, our results are different and extend the above results to some extent. Specially, motivated by [34], we suppose Q(x) changes sign on Ω and f(x,u)≡0 for (1.5). Since (1.1) is critical growth, which leads to the cause of the lack of compactness of the embedding H1(Ω)↪L6(Ω), we overcome this difficulty by using P.Lions concentration compactness principle [19]. Moreover, note that Q(x) changes sign on Ω, how to estimate the level of the mountain pass is another difficulty.
We define the energy functional corresponding to problem (1.1) by
Iλ(u)=12‖u‖2+b4(∫Ω|∇u|2dx)2−16∫ΩQ(x)|u|6dx−λq∫ΩP(x)|u|qdx. |
A weak solution of problem (1.1) is a function u∈H1(Ω) and for all φ∈H1(Ω) such that
∫Ω(a∇u∇φ+uφ)dx+b∫Ω|∇u|2dx∫Ω∇u∇φdx=∫ΩQ(x)|u|4uφdx+λ∫ΩP(x)|u|q−2uφdx. |
Our main results are the following:
Theorem 1.1. Assume that 1<q<2 and Q(x) changes sign on Ω. Then there exists Λ0>0 such that for every λ∈(0,Λ0), problem (1.1) has at least one nontrivial solution.
Theorem 1.2. Assume that 1<q<2, 3−q<β<6−q2 and Q(x) changes sign on Ω, there exists Λ∗>0 such that for all λ∈(0,Λ∗). Then problem (1.1) has at least two nontrivial solutions.
Throughout this paper, we make use of the following notations:
● The space H1(Ω) is equipped with the norm ‖u‖2H1(Ω)=∫Ω(|∇u|2+u2)dx, the norm in Lp(Ω) is denoted by ‖⋅‖p.
● Define ‖u‖2=∫Ω(a|∇u|2+u2)dx for u∈H1(Ω). Note that ‖⋅‖ is an equivalent norm on H1(Ω) with the standard norm.
● Let D1,2(R3) is the completion of C∞0(R3) with respect to the norm ‖u‖2D1,2(R3)=∫R3|∇u|2dx.
● 0<QM=maxx∈ˉΩQ(x), 0<Qm=maxx∈∂ΩQ(x).
● Ω+={x∈Ω:Q(x)>0} and Ω−={x∈Ω:Q(x)<0}.
● C,C1,C2,… denote various positive constants, which may vary from line to line.
● We denote by Sρ (respectively, Bρ) the sphere (respectively, the closed ball) of center zero and radius ρ, i.e. Sρ={u∈H1(Ω):‖u‖=ρ}, Bρ={u∈H1(Ω):‖u‖≤ρ}.
● Let S be the best constant for Sobolev embedding H1(Ω)↪L6(Ω), namely
S=infu∈H1(Ω)∖{0}∫Ω(a|∇u|2+u2)dx(∫Ω|u|6dx)1/3. |
● Let S0 be the best constant for Sobolev embedding D1,2(R3)↪L6(R3), namely
S0=infu∈D1,2(R3)∖{0}∫R3|∇u|2dx(∫R3|u|6dx)1/3. |
In this section, we firstly show that the functional Iλ(u) has a mountain pass geometry.
Lemma 2.1. There exist constants r,ρ,Λ0>0 such that the functional Iλ satisfies the following conditions for each λ∈(0,Λ0):
(i) Iλ|u∈Sρ≥r>0; infu∈BρIλ(u)<0.
(ii) There exists e∈H1(Ω) with ‖e‖>ρ such that Iλ(e)<0.
Proof. (i) From (P1), by the H¨older inequality and the Sobolev inequality, for all u∈H1(Ω) one has
∫ΩP(x)|u|qdx≤P(x0)∫Ω|u|qdx≤P(x0)|Ω|6−q6S−q2‖u‖q, | (2.1) |
and there exists a constant C>0, we get
|∫ΩQ(x)|u|6dx|≤C∫Ω|u|6dx≤CS−3‖u‖6. | (2.2) |
Hence, combining (2.1) and (2.2), we have the following estimate
Iλ(u)=12‖u‖2+b4(∫Ω|∇u|2dx)2−16∫ΩQ(x)|u|6dx−λq∫ΩP(x)|u|qdx≥12‖u‖2−C6∫Ω|u|6dx−λqP(x0)|Ω|6−q6S−q2‖u‖q≥‖u‖q(12‖u‖2−q−C6S−3‖u‖6−q−λqP(x0)|Ω|6−q6S−q2). |
Set h(t)=12t2−q−C6S−3t6−q for t>0, then there exists a constant ρ=(3(2−q)S3C(6−q))14>0 such that maxt>0h(t)=h(ρ)>0. Letting Λ0=qSq2P(x0)|Ω|6−q6h(ρ), there exists a constant r>0 such that Iλ|u∈Sρ≥r for every λ∈(0,Λ0). Moreover, for all u∈H1(Ω)∖{0}, we have
limt→0+Iλ(tu)tq=−λq∫ΩP(x)|u|qdx<0. |
So we obtain Iλ(tu)<0 for every u≠0 and t small enough. Therefore, for ‖u‖ small enough, one has
m≜ |
(\mathrm{ii}) Let v\in H^1(\Omega) be such that supp v\subset\Omega^+ , v\not\equiv 0 and t > 0 , we have
\begin{equation*} I_{\lambda}(tv) = \frac{t^2}{2}\|v\|^2+\frac{bt^4}{4}\left(\int_{\Omega}|\nabla v|^2dx\right)^{2}-\frac{t^6}{6}\int_{\Omega}Q(x)|v|^{6}dx-\frac{\lambda t^q}{q}\int_{\Omega}P(x)|v|^qdx\\ \rightarrow -\infty \end{equation*} |
as t\rightarrow\infty , which implies that I_{\lambda}(tv) < 0 for t > 0 large enough. Therefore, we can find e\in H^1(\Omega) with \|e\| > \rho such that I_{\lambda}(e) < 0 . The proof is complete.
Denote
\begin{equation*} \begin{cases} \Theta_1 = \frac{abS_0^3}{4Q_M}+\frac{b^3S_0^6}{24Q_M^2}+\frac{aS_0\sqrt{b^2S_0^4+4aS_0Q_M}}{6Q_M}+\frac{b^2S_0^4\sqrt{b^2S_0^4+4aS_0Q_M}}{24Q_M^2}, \\ \Theta_2 = \frac{abS_0^3}{16Q_m}+\frac{b^3S_0^6}{384Q_m^2}+\frac{aS_0\sqrt{b^2S_0^4+16aS_0Q_m}}{24Q_m}+\frac{b^2S_0^4\sqrt{b^2S_0^4+16aS_0Q_m}}{384Q_m^2}. \end{cases} \end{equation*} |
Then we have the following compactness result.
Lemma 2.2. Suppose that 1 < q < 2 . Then the functional I_{\lambda} satisfies the (PS)_{c_{\lambda}} condition for every c_{\lambda} < c_* = min \{\Theta_1-D\lambda^{\frac{2}{2-q}}, \Theta_2-D\lambda^{\frac{2}{2-q}}\} , where D = \frac{2-q}{3q}(\frac{6-q}{4}P(x_0)S^{-\frac{q}{2}}|\Omega|^{\frac{6-q}{6}})^{\frac{2}{2-q}} .
Proof. Let \{u_n\}\subset H^1(\Omega) be a (PS)_{c_\lambda} sequence for
\begin{equation} I_{\lambda}(u_n)\rightarrow c_{\lambda}\; \mathrm{and}\; I'_{\lambda}(u_n)\rightarrow 0\; \mathrm{as}\; n\rightarrow\infty. \end{equation} | (2.3) |
It follows from (2.1), (2.3) and the H \ddot{\mathrm{o}} lder inequality that
\begin{equation*} \begin{aligned} c_{\lambda}+1+o{(\|u_n\|)}&\geq I_{\lambda}(u_n)-\frac{1}{6}\langle I'_{\lambda}(u_n), u_n\rangle\\ &\geq \frac{1}{3}\|u_n\|^2+\frac{b}{12}\left(\int_{\Omega}|\nabla u_n|^2dx\right)^{2}\\ &\quad- \lambda\left(\frac{1}{q}-\frac{1}{6}\right)P(x_0)S^{-\frac{q}{2}}|\Omega|^{\frac{6-q}{6}}\|u_n\|^q\\ &\geq \frac{1}{3}\|u_n\|^2-\frac{\lambda(6-q)}{6q}P(x_0)S^{-\frac{q}{2}}|\Omega|^{\frac{6-q}{6}}\|u_n\|^q. \end{aligned} \end{equation*} |
Therefore \{u_n\} is bounded in H^1(\Omega) for all 1 < q < 2 . Thus, we may assume up to a subsequence, still denoted by \{u_n\} , there exists u\in H^{1}(\Omega) such that
\begin{align} \begin{cases} u_n\rightharpoonup u, \; \; \; \mathrm{weakly\; in}\; H^{1}(\Omega), \\ u_n\rightarrow u, \; \; \; \mathrm{strongly\; in}\; L^{p}(\Omega)\; (1\leq p \lt 6), \\ u_n(x)\rightarrow u(x), \; \; \; \mathrm{a.e.\; in}\; \Omega, \end{cases} \end{align} | (2.4) |
as n\rightarrow\infty . Next, we prove that u_n\rightarrow u strongly in H^1(\Omega) . By using the concentration compactness principle (see [19]), there exist some at most countable index set J , \delta_{x_j} is the Dirac mass at x_j\subset \bar{\Omega} and positive numbers \{\nu_j\} , \{\mu_j\} , j\in J , such that
\begin{equation*} \left. \begin{array}{rcl} & |u_n|^6dx\rightharpoonup d\nu = |u|^6dx+\sum\limits_{j\in J}\nu_j\delta_{x_j}, \\ & |\nabla u_n|^2dx\rightharpoonup d\mu \geq |\nabla u|^2dx+\sum\limits_{j\in J}\mu_{j}\delta_{x_j}. \end{array} \right. \end{equation*} |
Moreover, numbers \nu_{j} and \mu_{j} satisfy the following inequalities
\begin{equation} \left. \begin{array}{rcl} & S_0\nu_{j}^\frac{1}{3}\leq \mu_{j}\; \; \mathrm{if}\; x_{j}\in \Omega, \\ & \frac{S_0}{2^\frac{2}{3}}\nu_{j}^\frac{1}{3}\leq \mu_{j}\; \; \; \mathrm{if}\; x_{j}\in \partial \Omega. \end{array} \right. \end{equation} | (2.5) |
For \varepsilon > 0 , let \phi_{\varepsilon, j}(x) be a smooth cut-off function centered at x_j such that 0\leq\phi_{\varepsilon, j}\leq 1, |\nabla \phi_{\varepsilon, j}|\leq\frac{2}{\varepsilon} , and
\begin{equation*} \phi_{\varepsilon, j}(x) = \begin{cases} 1, &\text{ in $B(x_j, \frac{\varepsilon}{2})\cap\bar{\Omega}$}, \\ 0, &\text{ in $\Omega \setminus B(x_j, \varepsilon)$.} \end{cases} \end{equation*} |
There exists a constant C > 0 such that
\begin{equation*} \lim\limits_{\varepsilon\rightarrow 0}\lim\limits_{n\rightarrow \infty}\int_{\Omega}P(x)|u_n|^q\phi_{\varepsilon, j}dx\leq P(x_0)\lim\limits_{\varepsilon\rightarrow 0}\lim\limits_{n\rightarrow \infty}\int_{B(x_j, \varepsilon)}|u_n|^{q}dx = 0. \end{equation*} |
Since |\nabla \phi_{\varepsilon, j}|\leq\frac{2}{\varepsilon} , by using the H \ddot{\mathrm{o}} lder inequality and L^2(\Omega) -convergence of \{u_n\} , we have
\begin{equation*} \begin{aligned} &\quad\lim\limits_{\varepsilon\rightarrow 0}\lim\limits_{n\rightarrow \infty}\left(a+b\int_{\Omega}|\nabla u_n|^2dx\right)\int_{\Omega}\langle\nabla u_n, \nabla \phi_{\varepsilon, j}\rangle u_n dx\\ &\leq C \lim\limits_{\varepsilon\rightarrow 0}\lim\limits_{n\rightarrow \infty} \left(\int_{\Omega}|\nabla u_n|^2dx\right)^{\frac{1}{2}}\left(\int_{\Omega}|u_n|^2|\nabla \phi_{\varepsilon, j}|^2dx\right)^{\frac{1}{2}}\\ &\leq C \lim\limits_{\varepsilon\rightarrow 0} \left(\int_{B(x_j, \varepsilon)}|u|^6dx\right)^{\frac{1}{6}}\left(\int_{B(x_j, \varepsilon)}|\nabla \phi_{\varepsilon, j}|^3dx\right)^{\frac{1}{3}}\\ &\leq C \lim\limits_{\varepsilon\rightarrow 0} \left(\int_{B(x_j, \varepsilon)}|u|^6dx\right)^{\frac{1}{6}}\left(\int_{B(x_j, \varepsilon)}\left(\frac{2}{\varepsilon}\right)^3dx\right)^{\frac{1}{3}}\\ &\leq C_1 \lim\limits_{\varepsilon\rightarrow 0} \left(\int_{B(x_j, \varepsilon)}|u|^6dx\right)^{\frac{1}{6}}\\ & = 0, \end{aligned} \end{equation*} |
where C_1 > 0 , and we also derive that
\begin{equation*} \lim\limits_{\varepsilon\rightarrow 0}\lim\limits_{n\rightarrow\infty} \int_{\Omega}|\nabla u_n|^2 \phi_{\varepsilon, j}dx \geq\lim\limits_{\varepsilon\rightarrow 0}\int_{\Omega}|\nabla u|^2\phi_{\varepsilon, j}dx+\mu_j = \mu_j, \end{equation*} |
\begin{equation*} \lim\limits_{\varepsilon\rightarrow 0}\lim\limits_{n\rightarrow\infty} \int_{\Omega}Q(x)|u_n|^6 \phi_{\varepsilon, j}dx = \lim\limits_{\varepsilon\rightarrow 0}\int_{\Omega}Q(x)|u|^6\phi_{\varepsilon, j}dx+Q(x_j)\nu_j = Q(x_j)\nu_j, \end{equation*} |
\begin{equation*} \lim\limits_{\varepsilon\rightarrow 0}\lim\limits_{n\rightarrow\infty} \int_{\Omega}u_n^2 \phi_{\varepsilon, j}dx = \lim\limits_{\varepsilon\rightarrow 0}\int_{\Omega}u^2\phi_{\varepsilon, j}dx \leq\lim\limits_{\varepsilon\rightarrow 0}\int_{B(x_j, \varepsilon)}u^2dx = 0. \end{equation*} |
Noting that u_n\phi_{\varepsilon, j} is bounded in H^1(\Omega) uniformly for n , taking the test function \varphi = u_n\phi_{\varepsilon, j} in (2.3), from the above information, one has
\begin{equation*} \begin{aligned} 0& = \lim\limits_{\varepsilon\rightarrow 0}\lim\limits_{n\rightarrow \infty}\langle I'_{\lambda}(u_n), u_n\phi_{\varepsilon, j}\rangle\\ & = \lim\limits_{\varepsilon\rightarrow 0}\lim\limits_{n\rightarrow \infty}\bigg\{\left(a+b\int_{\Omega}|\nabla u_n|^2dx\right)\int_{\Omega}\langle\nabla u_n, \nabla(u_n\phi_{\varepsilon, j})\rangle dx+\int_{\Omega}u_n^2 \phi_{\varepsilon, j}dx\\ &\quad- \int_{\Omega}Q(x)|u_n|^6\phi_{\varepsilon, j}dx -\lambda\int_{\Omega}P(x)|u_n|^q\phi_{\varepsilon, j}dx \bigg\}\\ & = \lim\limits_{\varepsilon\rightarrow 0}\lim\limits_{n\rightarrow \infty}\bigg\{\left(a+b\int_{\Omega}|\nabla u_n|^2dx\right)\int_{\Omega}\left(|\nabla u_n|^2 \phi_{\varepsilon, j}+\langle\nabla u_n, \nabla \phi_{\varepsilon, j}\rangle u_n\right)dx\\ &\quad-\int_{\Omega}Q(x)|u_n|^6\phi_{\varepsilon, j}dx\bigg\}\\ &\geq \lim\limits_{\varepsilon\rightarrow 0}\bigg\{\left(a+b\int_{\Omega}|\nabla u|^2dx+b\mu_{j}\right)\left(\int_{\Omega}|\nabla u|^2 \phi_{\varepsilon, j}dx+\mu_{j}\right)\\ &\quad- \int_{\Omega}Q(x)|u|^6\phi_{\varepsilon, j}dx-Q(x_{j})\nu_{j}\bigg\}\\ &\geq \left(a+b\mu_{j}\right)\mu_{j}-Q(x_{j})\nu_{j}, \end{aligned} \end{equation*} |
so that
Q(x_{j})\nu_{j}\geq(a+b\mu_{j})\mu_{j}, |
which shows that \{u_n\} can only concentrate at points x_{j} where Q(x_{j}) > 0 . If \nu_{j} > 0 , by (2.5) we get
\begin{equation} \left. \begin{array}{rcl} & \nu_{j}^\frac{1}{3}\geq\frac{bS_0^2+\sqrt{b^2S_0^4+4aS_0Q_M}}{2Q_M}\; \; \mathrm{if}\; x_{j}\in \Omega, \\ & \nu_{j}^\frac{1}{3}\geq\frac{bS_0^2+\sqrt{b^2S_0^4+16aS_0Q_m}}{2^\frac{7}{3}Q_m}\; \; \mathrm{if}\; x_{j}\in \partial \Omega. \end{array} \right. \end{equation} | (2.6) |
From (2.5) and (2.6), we have
\begin{equation} \left. \begin{array}{rcl} & \mu_{j}\geq\frac{bS_0^3+\sqrt{b^2S_0^6+4aS_0^3Q_M}}{2Q_M}\; \; \mathrm{if}\; x_{j}\in \Omega, \\ & \mu_{j}\geq\frac{bS_0^3+\sqrt{b^2S_0^6+16aS_0^3Q_m}}{8Q_m}\; \; \mathrm{if}\; x_{j}\in \partial \Omega. \end{array} \right. \end{equation} | (2.7) |
To proceed further we show that (2.7) is impossible. To obtain a contradiction assume that there exists j_0\in J such that \mu_{j_0}\geq\frac{bS_0^3+\sqrt{b^2S_0^6+4aS_0^3Q_M}}{2Q_M} and x_{j_0}\in \Omega . By (2.1), (2.3) and (2.4), one has
\begin{equation*} \begin{aligned} c_\lambda& = \lim\limits_{n\rightarrow \infty}\bigg\{I_{\lambda}(u_n)-\frac{1}{6}\langle I'_{\lambda}(u_n), u_n\rangle\bigg\}\\ & = \lim\limits_{n\rightarrow \infty}\bigg\{\frac{a}{3}\int_{\Omega}|\nabla u_n|^2dx+\frac{b}{12}\left(\int_{\Omega}|\nabla u_n|^2dx\right)^2\\ &\quad+ \frac{1}{3}\int_{\Omega}u_n^2dx-\lambda\frac{6-q}{6q}\int_{\Omega}P(x)|u_n|^qdx\bigg\}\\ &\geq \frac{a}{3}\left(\int_{\Omega}|\nabla u|^2dx+\sum\limits_{j\in J}\mu_j\right)+\frac{b}{12}\left(\int_{\Omega}|\nabla u|^2dx+\sum\limits_{j\in J}\mu_j\right)^2\\ &\quad+ \frac{1}{3}\int_{\Omega}u^2dx-\lambda\frac{6-q}{6q}P(x_0)S^{-\frac{q}{2}}|\Omega|^{\frac{6-q}{6}}\|u\|^q\\ &\geq \frac{a}{3}\mu_{j_0}+\frac{b}{12}\mu^2_{j_0}+\frac{1}{3}\|u\|^2-\lambda\frac{6-q}{6q}P(x_0)S^{-\frac{q}{2}}|\Omega|^{\frac{6-q}{6}}\|u\|^q. \end{aligned} \end{equation*} |
Set
\begin{equation*} g(t) = \frac{1}{3}t^2-\lambda\frac{6-q}{6q}P(x_0)S^{-\frac{q}{2}}|\Omega|^{\frac{6-q}{6}}t^{q}, \; \; t \gt 0, \end{equation*} |
then
\begin{equation*} g'(t) = \frac{2}{3}t-\lambda \frac{6-q}{6}P(x_0)S^{-\frac{q}{2}}|\Omega|^{\frac{6-q}{6}}t^{q-1} = 0, \end{equation*} |
we can deduce that \min_{t\geq0}g(t) attains at t_0 > 0 and
\begin{equation*} t_0 = \left(\lambda\frac{6-q}{4}P(x_0)S^{-\frac{q}{2}}|\Omega|^{\frac{6-q}{6}}\right)^{\frac{1}{2-q}}. \end{equation*} |
Consequently, we obtain
\begin{equation*} \begin{aligned} c_\lambda&\geq \frac{abS_0^3}{4Q_M}+\frac{b^3S_0^6}{24Q_M^2}+\frac{aS_0\sqrt{b^2S_0^4+4aS_0Q_M}}{6Q_M}\\ &\quad+\frac{b^2S_0^4\sqrt{b^2S_0^4+4aS_0Q_M}}{24Q_M^2}- D\lambda^{\frac{2}{2-q}}\\ & = \Theta_1-D\lambda^{\frac{2}{2-q}}, \end{aligned} \end{equation*} |
where D = \frac{2-q}{3q}\left(\frac{6-q}{4}P(x_0)S^{-\frac{q}{2}}|\Omega|^{\frac{6-q}{6}}\right)^{\frac{2}{2-q}} . If \mu_{j_0}\geq\frac{bS_0^3+\sqrt{b^2S_0^6+16aS_0^3Q_m}}{8Q_m} and x_{j_0}\in \partial\Omega , then, by the similar calculation, we also get
\begin{equation*} \begin{aligned} c_\lambda&\geq \frac{abS_0^3}{16Q_m}+\frac{b^3S_0^6}{384Q_m^2}+\frac{aS_0\sqrt{b^2S_0^4+16aS_0Q_m}}{24Q_m}\\ &\quad+\frac{b^2S_0^4\sqrt{b^2S_0^4+16aS_0Q_m}}{384Q_m^2}-D\lambda^{\frac{2}{2-q}}\\ & = \Theta_2-D\lambda^{\frac{2}{2-q}}. \end{aligned} \end{equation*} |
Let c_* = \min\{\Theta_1-D\lambda^{\frac{2}{2-q}}, \Theta_2-D\lambda^{\frac{2}{2-q}}\} , from the above information, we deduce that c_\lambda\geq c_* . It contradicts our assumption, so it indicates that \nu_j = \mu_j = 0 for every j\in J , which implies that
\begin{equation} \int_{\Omega}|u_n|^6dx\rightarrow\int_{\Omega}|u|^6dx \end{equation} | (2.8) |
as n\rightarrow\infty . Now, we may assume that \int_{\Omega}|\nabla u_n|^2dx\rightarrow A^2 and \int_{\Omega}|\nabla u|^2dx\leq A^2 , by (2.3), (2.4) and (2.8), one has
\begin{equation*} \begin{aligned} 0& = \lim\limits_{n\rightarrow \infty}\langle I_\lambda'(u_n), u_n-u\rangle\\ & = \lim\limits_{n\rightarrow\infty}\bigg[\left(a+b\int_{\Omega}|\nabla u_n|^2dx\right)\left(\int_{\Omega}|\nabla u_n|^2dx-\int_{\Omega}\nabla u_n \nabla udx\right)\\ &\quad+\int_{\Omega}u_n(u_n-u)dx-\int_{\Omega}Q(x)|u_n|^5(u_n-u)dx-\lambda\int_{\Omega}P(x)|u_n|^{q-1}(u_n-u)dx\bigg]\\ & = \left(a+bA^2\right)\left(A^2-\int_{\Omega}|\nabla u|^2dx\right). \end{aligned} \end{equation*} |
Then, we obtain that u_n\rightarrow u in H^1(\Omega) . The proof is complete.
As well known, the function
\begin{equation*} U_{\varepsilon, y}(x) = \frac{(3\varepsilon^2)^{\frac{1}{4}}}{(\varepsilon^2+|x-y|^2)^{\frac{1}{2}}}, \; \mathrm{for\; any}\; \varepsilon \gt 0, \end{equation*} |
satisfies
\begin{equation*} -\Delta U_{\varepsilon, y} = U_{\varepsilon, y}^{5}\; \; \mathrm{in}\; \mathbb{R}^3, \end{equation*} |
and
\begin{equation*} \int_{\mathbb{R}^3}|\nabla U_{\varepsilon, y}|^2dx = \int_{\mathbb{R}^3}|U_{\varepsilon, y}|^6dx = S_0^{\frac{3}{2}}. \end{equation*} |
Let \phi\in C^1(\mathbb{R}^3) such that \phi(x) = 1 on B(x_M, \frac{R}{2}) , \phi(x) = 0 on \mathbb{R}^3-B(x_M, R) and 0\leq\phi(x)\leq 1 on \mathbb{R}^3 , we set v_\varepsilon(x) = \phi(x)U_{\varepsilon, x_M}(x) . We may assume that Q(x) > 0 on B(x_M, R) for some R > 0 such that B(x_M, R)\subset \Omega . From [4], we have
\begin{equation} \begin{cases} \|\nabla v_\varepsilon\|_2^2 = S_0^\frac{3}{2}+O(\varepsilon), \\ \|v_\varepsilon\|_6^6 = S_0^\frac{3}{2}+O(\varepsilon^{3}), \\ \|v_\varepsilon\|_2^2 = O(\varepsilon), \\ \|v_\varepsilon\|^2 = aS_0^{\frac{3}{2}}+O(\varepsilon). \end{cases} \end{equation} | (2.9) |
Moreover, by [28], we get
\begin{equation} \begin{cases} \|\nabla v_\varepsilon\|_2^4\leq S_0^3+O(\varepsilon), \\ \|\nabla v_\varepsilon\|_2^8\leq S_0^6+O(\varepsilon), \\ \|\nabla v_\varepsilon\|_2^{12}\leq S_0^9+O(\varepsilon). \end{cases} \end{equation} | (2.10) |
Then we have the following Lemma.
Lemma 2.3. Suppose that 1 < q < 2 , 3-q < \beta < \frac{6-q}{2} , Q_M > 4Q_m , (Q_1) and (Q_2) , then \sup_{t\geq0}I_{\lambda}(tv_\varepsilon) < \Theta_1-D\lambda^{\frac{2}{2-q}}.
Proof. By Lemma 2.1, one has I_{\lambda}(tv_{\varepsilon})\rightarrow-\infty as t\rightarrow\infty and I_{\lambda}(tv_{\varepsilon}) < 0 as t\rightarrow 0 , then there exists t_{\varepsilon} > 0 such that I_{\lambda}(t_{\varepsilon} v_{\varepsilon}) = \sup_{t > 0}I_{\lambda}(tv_{\varepsilon})\geq r > 0 . We can assume that there exist positive constants t_1, t_2 > 0 and 0 < t_1 < t_{\varepsilon} < t_2 < +\infty . Let I_{\lambda}(t_{\varepsilon}v_{\varepsilon}) = \beta(t_{\varepsilon}v_{\varepsilon})-\lambda\psi(t_{\varepsilon}v_{\varepsilon}) , where
\begin{equation*} \beta(t_{\varepsilon}v_{\varepsilon}) = \frac{ t_{\varepsilon}^2}{2}\| v_{\varepsilon}\|^2+\frac{b t_{\varepsilon}^4}{4}\|\nabla v_{\varepsilon}\|_2^{4}-\frac{t_{\varepsilon}^6}{6}\int_{\Omega}Q(x)|v_{\varepsilon}|^{6}dx, \end{equation*} |
and
\begin{equation*} \psi(t_{\varepsilon}v_{\varepsilon}) = \frac{ t_{\varepsilon}^q}{q}\int_{\Omega}P(x)|v_{\varepsilon}|^qdx. \end{equation*} |
Now, we set
\begin{equation*} h(t) = \frac{t^2}{2}\| v_{\varepsilon}\|^2+\frac{b t^4}{4}\|\nabla v_{\varepsilon}\|_2^{4}-\frac{t^6 }{6}\int_{\Omega}Q(x)|v_{\varepsilon}|^{6}dx. \end{equation*} |
It is clear that \lim_{t\rightarrow0}h(t) = 0 and \lim_{t\rightarrow\infty}h(t) = -\infty . Therefore there exists T_{1} > 0 such that h(T_{1}) = \max_{t\geq0}h(t) , that is
\begin{equation*} h'(t)|_{T_{1}} = T_{1}\| v_{\varepsilon}\|^2+bT_{1}^3\|\nabla v_{\varepsilon}\|_2^{4}- T_{1}^5 \int_{\Omega}Q(x)|v_{\varepsilon}|^{6}dx = 0, \end{equation*} |
from which we have
\begin{equation} \| v_{\varepsilon}\|^2+bT_{1}^2\|\nabla v_{\varepsilon}\|_2^{4} = T_{1}^4 \int_{\Omega}Q(x)|v_{\varepsilon}|^{6}dx. \end{equation} | (2.11) |
By (2.11) we obtain
\begin{equation*} T_{1}^2 = \frac{b\|\nabla v_{\varepsilon}\|_2^4+\sqrt{b^2\|\nabla v_{\varepsilon}\|_2^8+4\| v_{\varepsilon}\|^2\int_{\Omega}Q(x)|v_{\varepsilon}|^{6}dx}}{2 \int_{\Omega}Q(x)|v_{\varepsilon}|^{6}dx}. \end{equation*} |
In addition, by (Q_2) , for all \eta > 0 , there exists \rho > 0 such that |Q(x)-Q_M| < \eta|x-x_M| for 0 < |x-x_M| < \rho , for \varepsilon > 0 small enough, we have
\begin{align*} \left|\int_{\Omega}Q(x)v_\varepsilon^6dx-\int_{\Omega}Q_Mv_\varepsilon^6dx\right| &\leq \int_{\Omega}|Q(x)-Q_M|v_\varepsilon^6dx\\ & \lt \int_{B(x_M, \rho)}\eta|x-x_M|\frac{(3\varepsilon^2)^\frac{3}{2}}{(\varepsilon^2+|x-x_M|^2)^3}dx\\ &\quad+ C\int_{\Omega\setminus{B(x_M, \rho)}}\frac{(3\varepsilon^2)^\frac{3}{2}}{(\varepsilon^2+|x-x_M|^2)^3}dx\\ &\leq C \eta\varepsilon^3\int_{0}^{\rho}\frac{r^3}{(\varepsilon^2+r^2)^{3}}dr+C\varepsilon^3\int_\rho^{R}\frac{r^2}{(\varepsilon^2+r^2)^{3}}dr\\ &\leq C \eta\varepsilon\int_{0}^{\rho/\varepsilon}\frac{t^3}{(1+t^2)^3}dt+C\int_{\rho/\varepsilon}^{R/\varepsilon}\frac{t^2}{(1+t^2)^3}dt\\ &\leq C_1 \eta\varepsilon+C_2\varepsilon^3, \end{align*} |
where C_1, C_2 > 0 (independent of \eta , \varepsilon ). From this we derive that
\begin{equation} \limsup\limits_{\varepsilon\rightarrow 0}\frac{|\int_{\Omega}Q(x)v_\varepsilon^6dx-\int_{\Omega}Q_Mv_\varepsilon^6dx|}{\varepsilon}\leq C_1\eta. \end{equation} | (2.12) |
Then from the arbitrariness of \eta > 0 , by (2.9) and (2.12), one has
\begin{equation} \int_{\Omega}Q(x)|v_\varepsilon|^6dx = Q_M\int_{\Omega}|v_\varepsilon|^6dx+o(\varepsilon) = Q_MS_0^{\frac{3}{2}}+o(\varepsilon). \end{equation} | (2.13) |
Hence, it follows from (2.9), (2.10) and (2.13) that
\begin{align*} \beta(t_{\varepsilon}v_{\varepsilon})&\leq h(T_{1})\\ & = T_{1}^2\left(\frac{1}{3}\| v_{\varepsilon}\|^2+\frac{b T_{1}^2}{12}\|\nabla v_{\varepsilon}\|_2^{4}\right)\\ & = \frac{b\|\nabla v_\varepsilon\|_2^4\|v_\varepsilon\|^2}{4\int_{\Omega}Q(x)|v_\varepsilon|^6dx}+\frac{b^3\|\nabla v_\varepsilon\|_2^{12}}{24(\int_{\Omega}Q(x)|v_\varepsilon|^6dx)^2}\\ &\quad+ \frac{\|v_\varepsilon\|^2\sqrt{b^2\|\nabla v_\varepsilon\|_2^8+4\|v_\varepsilon\|^2\int_{\Omega}Q(x)|v_\varepsilon|^6dx}}{6\int_{\Omega}Q(x)|v_\varepsilon|^6dx}\\ &\quad+ \frac{b^2\|\nabla v_\varepsilon\|_2^8\sqrt{b^2\|\nabla v_\varepsilon\|_2^8+4\|v_\varepsilon\|^2\int_{\Omega}Q(x)|v_\varepsilon|^6dx}}{24(\int_{\Omega}Q(x)|v_\varepsilon|^6dx)^2}\\ &\leq \frac{b(S_0^3+O(\varepsilon))(aS_0^\frac{3}{2}+O(\varepsilon))}{4(Q_MS_0^\frac{3}{2}+o(\varepsilon))} +\frac{b^3(S_0^9+O(\varepsilon))}{24(Q_MS_0^\frac{3}{2}+o(\varepsilon))^2}\\ &\quad+ \frac{(aS_0^\frac{3}{2}+O(\varepsilon))\sqrt{b^2(S_0^6+O(\varepsilon)) +4(aS_0^\frac{3}{2}+O(\varepsilon))(Q_MS_0^\frac{3}{2}+o(\varepsilon))}}{6(Q_MS_0^\frac{3}{2}+o(\varepsilon))}\\ &\quad+ \frac{b^2(S_0^6+O(\varepsilon))\sqrt{b^2(S_0^6+O(\varepsilon)) +4(aS_0^\frac{3}{2}+O(\varepsilon))(Q_MS_0^\frac{3}{2}+o(\varepsilon))}}{24(Q_MS_0^\frac{3}{2}+o(\varepsilon))^2}\\ &\leq \frac{abS_0^3}{4Q_M}+\frac{b^3S_0^6}{24Q_M^2}+\frac{aS_0\sqrt{b^2S_0^4+4aS_0Q_M}}{6Q_M}\\ &\quad+ \frac{b^2S_0^4\sqrt{b^2S_0^4+4aS_0Q_M}}{24Q_M^2}+C_3\varepsilon\\ & = \Theta_1+C_3\varepsilon, \end{align*} |
where the constant C_3 > 0 . According to the definition of v_\varepsilon , from [29], for \frac{R}{2} > \varepsilon > 0 , there holds
\begin{align} \psi(t_{\varepsilon}v_{\varepsilon}) &\geq \frac{1}{q}3^{\frac{q}{4}}t_1^q\int_{B(x_M, \frac{R}{2})} \frac{\sigma\varepsilon^\frac{q}{2}}{(\varepsilon^2+|x-x_M|^2)^{\frac{q}{2}}|x-x_M|^\beta}dx \\ &\geq C\varepsilon^\frac{q}{2}\int_0^{R/2}\frac{r^2}{(\varepsilon^2+r^2)^{\frac{q}{2}}r^\beta}dr \\ & = C\varepsilon^{\frac{6-q}{2}-\beta}\int_0^{R/2\varepsilon}\frac{t^2}{(1+t^2)^{\frac{q}{2}}t^\beta}dt \\ &\geq C\varepsilon^{\frac{6-q}{2}-\beta}\int_0^1t^{2-\beta}dt \\ & = C_4 \varepsilon^{\frac{6-q}{2}-\beta}, \end{align} | (2.14) |
where C_4 > 0 (independent of \varepsilon, \lambda ). Consequently, from the above information, we obtain
\begin{equation*} \begin{aligned} I_{\lambda}(t_{\varepsilon}v_{\varepsilon})& = \beta(t_{\varepsilon}v_{\varepsilon})-\lambda\psi(t_{\varepsilon}v_{\varepsilon})\\ &\leq \Theta_1+C_3\varepsilon-C_4\lambda \varepsilon^{\frac{6-q}{2}-\beta}\\ & \lt \Theta_1-D\lambda^{\frac{2}{2-q}}. \end{aligned} \end{equation*} |
Here we have used the fact that \beta > 3-q and let \varepsilon = \lambda^{\frac{2}{2-q}} , 0 < \lambda < \Lambda_1 = \min\{1, (\frac{C_3+D}{C_4})^{\frac{2-q}{6-2q-2\beta}}\} , then
\begin{equation} \begin{aligned} C_3\varepsilon-C_4\lambda \varepsilon^{{\frac{6-q}{2}-\beta}}& = C_3\lambda^{\frac{2}{2-q}}-C_4\lambda^{\frac{8-2q-2\beta}{2-q}}\\ & = \lambda^{\frac{2}{2-q}}(C_3-C_4\lambda^{\frac{6-2q-2\beta}{2-q}})\\ & \lt -D\lambda^{\frac{2}{2-q}}. \end{aligned} \end{equation} | (2.15) |
The proof is complete.
We assume that 0\in\partial\Omega and Q_m = Q(0) . Let \varphi\in C^1(\mathbb{R}^3) such that \varphi(x) = 1 on B(0, \frac{R}{2}) , \varphi(x) = 0 on \mathbb{R}^3-B(0, R) and 0\leq\varphi(x)\leq 1 on \mathbb{R}^3 , we set u_\varepsilon(x) = \varphi(x)U_\varepsilon(x) , the radius R is chosen so that Q(x) > 0 on B(0, R)\cap \Omega . If H(0) denotes the mean curvature of the boundary at 0 , then the following estimates hold (see [6] or [26])
\begin{equation} \begin{cases} \|u_\varepsilon\|_2^2 = O(\varepsilon), \\ \frac{\|\nabla u_{\varepsilon}\|_2^2}{\|u_{\varepsilon}\|_6^2}\leq\frac{S_0}{2^\frac{2}{3}}-A_3H(0)\varepsilon \log\frac{1}{\varepsilon}+O(\varepsilon), \end{cases} \end{equation} | (2.16) |
where A_3 > 0 is a constant. Then we have the following lemma.
Lemma 2.4. Suppose that 1 < q < 2 , 3-q < \beta < \frac{6-q}{2} , Q_M\leq4Q_m , H(0) > 0 , Q is positive somewhere on \partial\Omega , (Q_1) and (Q_3) , then \sup_{t\geq0}I_{\lambda}(tu_\varepsilon) < \Theta_2-D\lambda^{\frac{2}{2-q}}.
Proof. Similar to the proof of Lemma 2.3, we also have by Lemma 2.1, there exists t_{\varepsilon} > 0 such that I_{\lambda}(t_{\varepsilon} u_{\varepsilon}) = \sup_{t > 0}I_{\lambda}(tu_{\varepsilon})\geq r > 0 . We can assume that there exist positive constants t_1, t_2 > 0 such that 0 < t_1 < t_{\varepsilon} < t_2 < +\infty . Let I_{\lambda}(t_{\varepsilon}u_{\varepsilon}) = A(t_{\varepsilon}u_{\varepsilon})-\lambda B(t_{\varepsilon}u_{\varepsilon}) , where
\begin{equation*} A(t_{\varepsilon}u_{\varepsilon}) = \frac{ t_{\varepsilon}^2}{2}\| u_{\varepsilon}\|^2+\frac{b t_{\varepsilon}^4}{4}\|\nabla u_{\varepsilon}\|_2^{4}-\frac{t_{\varepsilon}^6}{6}\int_{\Omega}Q(x)|u_{\varepsilon}|^{6}dx, \end{equation*} |
and
\begin{equation*} B(t_{\varepsilon}u_{\varepsilon}) = \frac{ t_{\varepsilon}^q}{q}\int_{\Omega}P(x)|u_{\varepsilon}|^qdx. \end{equation*} |
Now, we set
\begin{equation*} f(t) = \frac{t^2}{2}\| u_{\varepsilon}\|^2+\frac{b t^4}{4}\|\nabla u_{\varepsilon}\|_2^{4}-\frac{t^6}{6}\int_{\Omega}Q(x)|u_{\varepsilon}|^{6}dx. \end{equation*} |
Therefore, it is easy to see that there exists T_{2} > 0 such that f(T_{2}) = \max_{f\geq0}f(t) , that is
\begin{equation} f'(t)|_{T_{2}} = T_{2}\| u_{\varepsilon}\|^2+bT_{2}^3\|\nabla u_{\varepsilon}\|_2^{4}- T_{2}^5 \int_{\Omega}Q(x)|u_{\varepsilon}|^{6}dx = 0. \end{equation} | (2.17) |
From (2.17) we obtain
\begin{equation*} T_{2}^2 = \frac{b\|\nabla u_{\varepsilon}\|_2^4+\sqrt{b^2\|\nabla u_{\varepsilon}\|_2^8+4\| u_{\varepsilon}\|^2 \int_{\Omega}Q(x)|u_{\varepsilon}|^{6}dx}}{2\int_{\Omega}Q(x)|u_{\varepsilon}|^{6}dx}. \end{equation*} |
By the assumption (Q_3) , we have the expansion formula
\begin{equation} \int_{\Omega}Q(x)|u_\varepsilon|^6dx = Q_m\int_{\Omega}|u_\varepsilon|^6dx+o(\varepsilon). \end{equation} | (2.18) |
Hence, combining (2.16) and (2.18), there exists C_5 > 0 , such that
\begin{align*} A(t_{\varepsilon}u_{\varepsilon})&\leq f(T_{2})\\ & = T_{2}^2\left(\frac{1}{3}\| u_{\varepsilon}\|^2+\frac{b T_{2}^2}{12}\|\nabla u_{\varepsilon}\|_2^4\right)\\ & = \frac{b\|\nabla u_{\varepsilon}\|_2^4\|u_{\varepsilon}\|^2}{4\int_{\Omega}Q(x)|u_\varepsilon|^6dx} +\frac{b^3\|\nabla u_\varepsilon\|_2^{12}}{24(\int_{\Omega}Q(x)|u_{\varepsilon}|^6dx)^2}\\ &\quad+ \frac{\|u_\varepsilon\|^2\sqrt{b^2\|\nabla u_\varepsilon\|_2^8 +4\|u_\varepsilon\|^2\int_{\Omega}Q(x)|u_\varepsilon|^6dx}}{6\int_{\Omega}Q(x)|u_\varepsilon|^6dx}\\ &\quad+ \frac{b^2\|\nabla u_\varepsilon\|_2^8\sqrt{b^2\|\nabla u_\varepsilon\|_2^8 +4\|u_\varepsilon\|^2\int_{\Omega}Q(x)| u_\varepsilon|^6dx}}{24(\int_{\Omega}Q(x)|u_\varepsilon|^6dx)^2}\\ &\leq \frac{ab}{4Q_m}\left(\frac{\|\nabla u_{\varepsilon}\|_2^6}{\int_{\Omega}|u_\varepsilon|^6dx}+O(\varepsilon)\right)+\frac{b^3}{24Q_m^2}\left(\frac{\|\nabla u_\varepsilon\|_2^{12}}{(\int_{\Omega}|u_\varepsilon|^6dx)^2}+O(\varepsilon)\right)\\ &\quad+ \frac{a}{6Q_m}\left(\frac{\|\nabla u_\varepsilon\|_2^2}{(\int_{\Omega}|u_\varepsilon|^6dx)^\frac{1}{3}}\sqrt{\frac{b^2\|\nabla u_\varepsilon\|_2^8}{(\int_{\Omega}|u_\varepsilon|^6dx)^\frac{4}{3}}+\frac{4aQ_m\|\nabla u_\varepsilon\|_2^2}{(\int_{\Omega}|u_\varepsilon|^6dx)^\frac{1}{3}}}+O(\varepsilon)\right)\\ &\quad+ \frac{b^2}{24Q_m^2}\left(\frac{\|\nabla u_\varepsilon\|_2^8}{(\int_{\Omega}|u_\varepsilon|^6dx)^\frac{4}{3}}\sqrt{\frac{b^2\|\nabla u_\varepsilon\|_2^8}{(\int_{\Omega}|u_\varepsilon|^6dx)^\frac{4}{3}}+\frac{4aQ_m\|\nabla u_\varepsilon\|_2^2}{(\int_{\Omega}|u_\varepsilon|^6dx)^\frac{1}{3}}}+O(\varepsilon)\right)\\ &\leq \frac{abS_0^3}{16Q_m}+\frac{b^3S_0^6}{384Q_m^2}+\frac{aS_0\sqrt{b^2S_0^4+16aS_0Q_m}}{24Q_m}\\ &\quad+ \frac{b^2S_0^4\sqrt{b^2S_0^4+16aS_0Q_m}}{384Q_m^2}+C_5\varepsilon\\ & = \Theta_2+C_5\varepsilon. \end{align*} |
Consequently, by (2.14) and (2.15), similarly, there exists \Lambda_2 > 0 such that 0 < \lambda < \Lambda_2 , we get
\begin{equation*} \begin{aligned} I_{\lambda}(t_{\varepsilon}u_{\varepsilon})& = A(t_{\varepsilon}u_{\varepsilon})-\lambda B(t_{\varepsilon}u_{\varepsilon})\\ &\leq \Theta_2+C_5\varepsilon-C_6\lambda \varepsilon^{\frac{6-q}{2}-\beta}\\ & \lt \Theta_2-D\lambda^{\frac{2}{2-q}}. \end{aligned} \end{equation*} |
where C_6 > 0 (independent of \varepsilon, \lambda ). The proof is complete.
Theorem 2.5. Assume that 0 < \lambda < \Lambda_0 ( \Lambda_0 is as in Lemma 2.1) and 1 < q < 2 . Then problem (1.1) has a nontrivial solution u_\lambda with I_\lambda(u_\lambda) < 0 .
Proof. It follows from Lemma 2.1 that
\begin{equation*} m\triangleq\inf\limits_{u\in \overline{B_\rho(0)}}I_\lambda(u) \lt 0. \end{equation*} |
By the Ekeland variational principle [7], there exists a minimizing sequence \{u_n\}\subset {\overline{B_\rho(0)}} such that
\begin{equation*} I_\lambda(u_n)\leq \inf\limits_{u\in{\overline{B_\rho(0)}}}I_\lambda(u)+\frac{1}{n}, \; \; I_\lambda(v)\geq I_\lambda(u_n)-\frac{1}{n}\|v-u_n\|, \; \; v\in{\overline{B_\rho(0)}}. \end{equation*} |
Therefore, there holds I_\lambda(u_n)\rightarrow m and I_\lambda'(u_n)\rightarrow 0 . Since \{u_n\} is a bounded sequence and {\overline{B_\rho(0)}} is a closed convex set, we may assume up to a subsequence, still denoted by \{u_n\} , there exists u_{\lambda}\in{\overline{B_\rho(0)}}\subset H^{1}(\Omega) such that
\begin{align*} \begin{cases} u_n\rightharpoonup u_{\lambda}, \; \; \; \mathrm{weakly\; in}\; H^{1}(\Omega), \\ u_n\rightarrow u_{\lambda}, \; \; \; \mathrm{strongly\; in}\; L^{p}(\Omega), \; 1\leq p \lt 6, \\ u_n(x)\rightarrow u_{\lambda}(x), \; \; \; \mathrm{a.e.\; in}\; \Omega. \end{cases} \end{align*} |
By the lower semi-continuity of the norm with respect to weak convergence, one has
\begin{equation*} \begin{aligned} m&\geq\liminf\limits_{n\rightarrow \infty}\left[I_\lambda(u_n)-\frac{1}{6}\langle I_\lambda'(u_n), u_n\rangle\right]\\ & = \liminf\limits_{n\rightarrow \infty}\bigg[\frac{1}{3}\int_\Omega \left(a|\nabla u_n|^2+u_n^2\right)dx +\frac{b}{12}\left(\int_\Omega|\nabla u_n|^2dx\right)^2\\ &\quad+ \lambda\left(\frac{1}{6} -\frac{1}{q}\right)\int_\Omega P(x)|u_n|^{q}dx\bigg]\\ &\geq \frac{1}{3}\int_\Omega \left(a|\nabla u_\lambda|^2+u_\lambda^2\right)dx +\frac{b}{12}\left(\int_\Omega|\nabla u_\lambda|^2dx\right)^2\\ &\quad+ \lambda\left(\frac{1}{6} -\frac{1}{q}\right)\int_\Omega P(x)|u_\lambda|^{q}dx\\ & = I_\lambda(u_\lambda)-\frac{1}{6}\langle I_\lambda'(u_\lambda), u_\lambda\rangle = I_\lambda(u_\lambda) = m. \end{aligned} \end{equation*} |
Thus I_{\lambda}(u_\lambda) = m < 0 , by m < 0 < c_\lambda and Lemma 2.2, we can see that \nabla u_n\rightarrow \nabla u_\lambda in L^2(\Omega) and u_\lambda\not\equiv0 . Therefore, we obtain that u_\lambda is a weak solution of problem (1.1). Since I_\lambda(|u_\lambda|) = I_\lambda(u_\lambda) , which suggests that u_\lambda\geq0 , then u_\lambda is a nontrivial solution to problem (1.1). That is, the proof of Theorem 1.1 is complete.
Theorem 2.6. Assume that 0 < \lambda < \Lambda_{*} (\Lambda_{*} = \min\{\Lambda_0, \Lambda_1, \Lambda_2\}) , 1 < q < 2 and 3-q < \beta < \frac{6-q}{2} . Then the problem (1.1) has a nontrivial solution u_{1}\in H^1(\Omega) such that I_{\lambda}(u_{1}) > 0 .
Proof. Applying the mountain pass lemma [3] and Lemma 2.2, there exists a sequence \{u_n\}\subset H^1(\Omega) such that
\begin{equation*} I_{\lambda}(u_n)\rightarrow c_\lambda \gt 0\; \; \mathrm{and}\; \; I'_\lambda(u_n)\rightarrow0\; \mathrm{as}\; n\rightarrow \infty, \end{equation*} |
where
c_\lambda = \inf\limits_{\gamma\in\Gamma}\max\limits_{t\in[0, 1]}I_{\lambda}(\gamma(t)), |
and
\Gamma = \left\{\gamma\in C([0, 1], H^{1}(\Omega)): \gamma(0) = 0, \gamma(1) = e\right\}. |
According to Lemma 2.2, we know that \{u_n\}\subset H^1(\Omega) has a convergent subsequence, still denoted by \{u_n\} , such that u_n\rightarrow u_{1} in H^1(\Omega) as n\rightarrow\infty ,
\begin{equation*} I_{\lambda}(u_{1}) = \lim\limits_{n\rightarrow\infty}I_{\lambda}(u_n) = c_\lambda \gt r \gt 0, \end{equation*} |
which implies that u_{1}\not\equiv0 . Therefore, from the continuity of I'_\lambda , we obtain that u_{1} is a nontrivial solution of problem (1.1) with I_{\lambda}(u_{1}) > 0 . Combining the above facts with Theorem 2.5 the proof of Theorem 1.2 is complete.
In this paper, we consider a class of Kirchhoff type equations with Neumann conditions and critical growth. Under suitable assumptions on Q(x) and P(x) , using the variational method and the concentration compactness principle, we proved the existence and multiplicity of nontrivial solutions.
This research was supported by the National Natural Science Foundation of China (Grant Nos. 11661021 and 11861021). Authors are grateful to the referees for their very constructive comments and valuable suggestions.
The authors declare no conflict of interest in this paper.
[1] | [ M. E. Alexander,S. M. Moghadas, Bifurcation analysis of SIRS epidemic model with generalized incidence, SIAM J. Appl. Math., 65 (2005): 1794-1816. |
[2] | [ R. M. Anderson and R. M. May, Infectious Diseases of Humans: Dynamics and Control, Oxford University Press, New York, 1991. |
[3] | [ J. Arino,F. Brauer,P. van den Driessche,J. Watmoughd,J. Wu, A model for influenza with vaccination and antiviral treatment, J. Theor. Biol., 253 (2008): 118-130. |
[4] | [ L. Bourouiba,A. Teslya,J. Wu, Highly pathogenic avian influenza outbreak mitigated by seasonal low pathogenic strains: insights from dynamic modeling, J. Theor. Biol., 271 (2011): 181-201. |
[5] | [ F. Brauer and C. Castillo-Chávez, Mathematical Models in Population Biology and Epidemiology, Springer-Verlag, New York, 2001. |
[6] | [ J. M. Cushing, A juvenile-adult model with periodic vital rates, J. Math. Biol., 53 (2006): 520-539. |
[7] | [ O. Diekmann,J. A. P. Heesterbeek,J. A. J. Metz, On the definition and the computation of the basic reproduction ratio R_0 in themodels for infectious disease in heterogeneous populations, J. Math. Biol., 28 (1990): 365-382. |
[8] | [ K. Dietz, The incidence of infectious diseases under the influence of seasonal fluctuations, Lecture Notes in Biomath, 11 (1976): 1-15, Berlin-Heidelberg-New York: Springer. |
[9] | [ D. J. D. Earn,P. Rohani,B. M. Bolker,B. T. Grenfell, A simple model for complex dynamical transitions in epidemics, Science, 287 (2000): 667-670. |
[10] | [ Z. Guo,L. Huang,X. Zou, Impact of discontinuous treatments on disease dynamics in an SIR epidemic model, Math. Biosci. Eng., 9 (2012): 97-110. |
[11] | [ H. W. Hethcote, The mathematics of infectious diseases, SIAM Rev., 42 (2000): 599-653. |
[12] | [ Y. Hsieh,J. Liu,Y. Tzeng,J. Wu, Impact of visitors and hospital staff on nosocomial transmission and spread to community, J. Theor. Biol., 356 (2014): 20-29. |
[13] | [ W. O. Kermack,A. G. McKendrick, Contributions to the mathematical theory of epidemics, Proc. R. Soc. London A: Math., 115 (1927): 700-721. |
[14] | [ I. M. Longini,M. E. Halloran,A. Nizam,Y. Yang, Containing pandemic influenza with antiviral agents, Am. J. Epidemiol., 159 (2004): 623-633. |
[15] | [ R. Olinky,A. Huppert,L. Stone, Seasonal dynamics and thresholds governing recurrent epidemics, J. Math. Biol., 56 (2008): 827-839. |
[16] | [ I. B. Schwartz, Small amplitude, long periodic out breaks in seasonally driven epidemics, J. Math. Biol., 30 (1992): 473-491. |
[17] | [ H. L. Smith, Multiple stable subharmonics for a periodic epidemic model, J. Math. Biol., 17 (1983): 179-190. |
[18] | [ H. L. Smith and P. Waltman, The Theory of the Chemostat: Dynamics of Microbial Competition, Cambridge Univ. Press, 1995. |
[19] | [ L. Stone,R. Olinky,A. Huppert, Seasonal dynamics of recurrent epidemics, Nature, 446 (2007): 533-536. |
[20] | [ Y. Tang,D. Huang,S. Ruan,W. Zhang, Coexistence of limit cycles and homoclinic loops in a SIRS model with a nonlinear incidence rate, SIAM J. Appl. Math., 69 (2008): 621-639. |
[21] | [ S. Towers,Z. Feng, Social contact patterns and control strategies for influenza in the elderly, Math. Biosci., 240 (2012): 241-249. |
[22] | [ S. Towers,K. Vogt Geisse,Y. Zheng,Z. Feng, Antiviral treatment for pandemic influenza: Assessing potential repercussions using a seasonally forced SIR model, J. Theor. Biol., 289 (2011): 259-268. |
[23] | [ P. Van den Driessche,J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002): 29-48. |
[24] | [ W. Wang,S. Ruan, Bifurcation in an epidemic model with constant removal rate of the infectives, J. Math. Anal. Appl., 291 (2004): 775-793. |
[25] | [ W. Wang,X. Zhao, Threshold dynamics for compartmental epidemic models in periodic environments, J. Dyn. Diff. Equ., 20 (2008): 699-717. |
[26] | [ D. Xiao, Dynamics and bifurcations on a class of population model with seasonal constant-yield harvesting, Discrete Contin. Dynam. Syst. -B, 21 (2016): 699-719. |
[27] | [ D. Xiao,S. Ruan, Global analysis of an epidemic model with nonmonotone incidence rate, Math. Biosci., 208 (2007): 419-429. |
[28] | [ F. Zhang,X. Zhao, A periodic epidemic model in a patchy environment, J. Math. Anal. Appl., 325 (2007): 496-516. |
[29] | [ X. Zhao, Dynamical Systems in Population Biology, Springer-Verlag, New York, 2003. |
1. | Ying Zhou, Jun Lei, Yue Wang, Zonghong Xiong, Positive solutions of a Kirchhoff–Schrödinger--Newton system with critical nonlocal term, 2022, 14173875, 1, 10.14232/ejqtde.2022.1.50 | |
2. | Deke Wu, Hongmin Suo, Jun Lei, Multiple Positive Solutions for Kirchhoff-Type Problems Involving Supercritical and Critical Terms, 2024, 23, 1575-5460, 10.1007/s12346-024-00999-w | |
3. | Jiaqing Hu, Anmin Mao, Multiple solutions to nonlocal Neumann boundary problem with sign‐changing coefficients, 2024, 0170-4214, 10.1002/mma.10419 | |
4. | A. Ahmed, Mohamed Saad Bouh Elemine Vall, MULTIPLICITY OF WEAK SOLUTIONS FOR A (P(X), Q(X))-KIRCHHOFF EQUATION WITH NEUMANN BOUNDARY CONDITIONS, 2024, 14, 2156-907X, 2441, 10.11948/20230449 | |
5. | Ahmed Ahmed, Mohamed Saad Bouh Elemine Vall, Weak solutions in anisotropic (α→(z),β→(z))-Laplacian Kirchhoff models, 2025, 1747-6933, 1, 10.1080/17476933.2025.2500342 |