Processing math: 100%
Research article

Nonlocal fractional p()-Kirchhoff systems with variable-order: Two and three solutions

  • Received: 19 July 2021 Accepted: 23 September 2021 Published: 26 September 2021
  • MSC : 35J91, 35A15, 35R11, 35J67

  • In this article, we consider the following nonlocal fractional Kirchhoff-type elliptic systems

    {M1(RN×RN|η(x)η(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)     dxdy+Ω|η|¯p(x)¯p(x)dx)(Δs()p()η|η|¯p(x)η)=λFη(x,η,ξ)+μGη(x,η,ξ),xΩ,M2(RN×RN|ξ(x)ξ(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)     dxdy+Ω|ξ|¯p(x)¯p(x)dx)(Δs()p()ξ|ξ|¯p(x)ξ)=λFξ(x,η,ξ)+μGξ(x,η,ξ),xΩ,η=ξ=0,xRNΩ,

    where M1(t),M2(t) are the models of Kirchhoff coefficient, Ω is a bounded smooth domain in RN, (Δ)s()p() is a fractional Laplace operator, λ,μ are two real parameters, F,G are continuous differentiable functions, whose partial derivatives are Fη,Fξ,Gη,Gξ. With the help of direct variational methods, we study the existence of solutions for nonlocal fractional p()-Kirchhoff systems with variable-order, and obtain at least two and three weak solutions based on Bonanno's and Ricceri's critical points theorem. The outstanding feature is the case that the Palais-Smale condition is not requested. The major difficulties and innovations are nonlocal Kirchhoff functions with the presence of the Laplace operator involving two variable parameters.

    Citation: Weichun Bu, Tianqing An, Guoju Ye, Yating Guo. Nonlocal fractional p()-Kirchhoff systems with variable-order: Two and three solutions[J]. AIMS Mathematics, 2021, 6(12): 13797-13823. doi: 10.3934/math.2021801

    Related Papers:

    [1] Zehra Yucedag . Variational approach for a Steklov problem involving nonstandard growth conditions. AIMS Mathematics, 2023, 8(3): 5352-5368. doi: 10.3934/math.2023269
    [2] Armin Hadjian, Juan J. Nieto . Existence of solutions of Dirichlet problems for one dimensional fractional equations. AIMS Mathematics, 2022, 7(4): 6034-6049. doi: 10.3934/math.2022336
    [3] Yun-Ho Kim . Multiple solutions to Kirchhoff-Schrödinger equations involving the $ p(\cdot) $-Laplace-type operator. AIMS Mathematics, 2023, 8(4): 9461-9482. doi: 10.3934/math.2023477
    [4] 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
    [5] Khaled Kefi, Mohammed M. Al-Shomrani . On multiple solutions for an elliptic problem involving Leray–Lions operator, Hardy potential and indefinite weight with mixed boundary conditions. AIMS Mathematics, 2025, 10(3): 5444-5455. doi: 10.3934/math.2025251
    [6] Zhongqian Wang, Xuejun Zhang, Mingliang Song . Three nonnegative solutions for Sturm-Liouville BVP and application to the complete Sturm-Liouville equations. AIMS Mathematics, 2023, 8(3): 6543-6558. doi: 10.3934/math.2023330
    [7] Adnan Khaliq, Mujeeb ur Rehman . Fixed point theorem combined with variational methods for a class of nonlinear impulsive fractional problems with derivative dependence. AIMS Mathematics, 2021, 6(2): 1943-1953. doi: 10.3934/math.2021118
    [8] Shaomin Wang, Cunji Yang, Guozhi Cha . On the variational principle and applications for a class of damped vibration systems with a small forcing term. AIMS Mathematics, 2023, 8(9): 22162-22177. doi: 10.3934/math.20231129
    [9] Khaled Kefi, Abdeljabbar Ghanmi, Abdelhakim Sahbani, Mohammed M. Al-Shomrani . Infinitely many solutions for a critical $ p(x) $-Kirchhoff equation with Steklov boundary value. AIMS Mathematics, 2024, 9(10): 28361-28378. doi: 10.3934/math.20241376
    [10] Zhongqian Wang, Dan Liu, Mingliang Song . Existence of three periodic solutions for a quasilinear periodic boundary value problem. AIMS Mathematics, 2020, 5(6): 6061-6072. doi: 10.3934/math.2020389
  • In this article, we consider the following nonlocal fractional Kirchhoff-type elliptic systems

    {M1(RN×RN|η(x)η(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)     dxdy+Ω|η|¯p(x)¯p(x)dx)(Δs()p()η|η|¯p(x)η)=λFη(x,η,ξ)+μGη(x,η,ξ),xΩ,M2(RN×RN|ξ(x)ξ(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)     dxdy+Ω|ξ|¯p(x)¯p(x)dx)(Δs()p()ξ|ξ|¯p(x)ξ)=λFξ(x,η,ξ)+μGξ(x,η,ξ),xΩ,η=ξ=0,xRNΩ,

    where M1(t),M2(t) are the models of Kirchhoff coefficient, Ω is a bounded smooth domain in RN, (Δ)s()p() is a fractional Laplace operator, λ,μ are two real parameters, F,G are continuous differentiable functions, whose partial derivatives are Fη,Fξ,Gη,Gξ. With the help of direct variational methods, we study the existence of solutions for nonlocal fractional p()-Kirchhoff systems with variable-order, and obtain at least two and three weak solutions based on Bonanno's and Ricceri's critical points theorem. The outstanding feature is the case that the Palais-Smale condition is not requested. The major difficulties and innovations are nonlocal Kirchhoff functions with the presence of the Laplace operator involving two variable parameters.



    The study of the existence and multiplicity solutions for nonlocal elliptic systems boundary value problems and variational problems has attracted intense research interests for several decades. In this paper, we investigate the existence of nontrivial solutions for a class of p()-Kirchhoff systems

    {M1(RN×RN|η(x)η(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)     dxdy+Ω|η|¯p(x)¯p(x)dx)(Δs()p()η|η|¯p(x)η)=λFη(x,η,ξ)+μGη(x,η,ξ),xΩ,M2(RN×RN|ξ(x)ξ(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)     dxdy+Ω|ξ|¯p(x)¯p(x)dx)(Δs()p()ξ|ξ|¯p(x)ξ)=λFξ(x,η,ξ)+μGξ(x,η,ξ),xΩ,η=ξ=0,xRNΩ, (1.1)

    where M1(t),M2(t):[0,+)[0,+) are the models of Kirchhoff coefficient, Ω is a bounded smooth domain in RN, λ,μ are two real parameters. F,G are differentiable and measurable functions in R2 for all xΩ, Fη,Fξ are the partial derivatives of F, Gη,Gξ are the partial derivatives of G, whose assumptions will be introduced later.

    The fractional p()-Laplace operator (Δ)s()p() with variable s()-order is defined as

    (Δ)s()p()υ(x):=P.V.RN|υ(x)υ(y)|p(x,y)2(υ(x)υ(y))|xy|N+p(x,y)s(x,y)dy, for allxRN, (1.2)

    where vC0(RN) and P.V. stands for the Cauchy principal value. (Δ)s()p() is a nonlocal operator of elliptic type, which is connected with the Sobolev space of variable exponent. Concerning this kind of operator problems, here we just list a few pieces of literatures, see [1,2,3,4]. Especially, Biswas et al. [5] firstly proved a continuous embedding result and Hardy-Littlewood-Sobolev-type result, and then the existence and multiplicity of solutions were obtained by variational approaches. When s() constant and p() constant, (Δ)s()p() in (1.2) reduce to the usual fractional Laplace operator (Δ)sp, see [6,7,8] for the essential knowledge.

    Throughout this paper, s(),p()C+(D) are two continuous functions that the following assumptions are satisfied.

    (S): s():¯Ω×¯Ω (0,1) is symmetric, namely, s(x,y)=s(y,x) for any  (x,y)¯Ω×¯Ω  with ¯s(x)=s(x,x);

    (P): p():¯Ω×¯Ω (1,+) is symmetric, namely, p(x,y)=p(y,x) for any (x,y)¯Ω ×¯Ω with ¯p(x)=p(x,x).

    Kirchhoff in [9] introduced the following Kirchhoff equation

    ρ2ξ(x)t2(p0h+E2LL0|ξ(x)t|2dx)2ξ(x)x2=0, (1.3)

    where ρ,p0,h,E,L with physical meaning are constants. A characteristic of Eq (1.3) is the fact that it contains a nonlocal item p0h+E2LL0|ξ(x)t|2dx, and then this type of equation is called nonlocal problem. From then on, the existence, multiplicity, uniqueness, and regularity of solutions for various Kirchhoff-type equations have been studied extensively, such as, see [10,11,12,13,14] for further details.

    The continuous Kirchhoff terms Mi(t):R+0R+,(i=1,2) are strictly increasing functions, which the following conditions are satisfied.

    (M): There exist mi=mi(ι)>0 and Mi=Mi(ι)>0,(i=1,2) for any ι>0 such that

    MiMi(t)mifor all t>ι,(i=1,2)

    and put

    ˜Mi(t)=t0Mi(ς)dςfor all tR+0,(i=1,2).

    In recent years, a multitude of scholars has devoted themselves to the study of Kirchhoff-type systems. When M1(t)=1 and M2(t)=1, Chen et al. in [15] consider the nontrivial solutions for the following elliptic systems.

    {(Δ)spη=λ|η|q2η+2αα+β|η|α2η|ξ|β,xΩ,(Δ)spξ=λ|ξ|q2ξ+2βα+β|ξ|β2ξ|η|α,xΩ,η=ξ=0,xRNΩ, (1.4)

    by utilizing Nehari manifold method and Fibering maps, they studied the existence of weak solutions for this kind of problem (1.4). Moreover, it has been applied in the local case s=1 in [16].

    In the famous literature [17], the three critical points theorem was established by Ricceri. Starting from this paper, Marano and Motreanu in [18] extended the result of Ricceri to non-differentiable functionals. Subsequently, Fan and Deng in [19] studied the version of Ricceri's result including variables exponents. Ricceri's result in [20] has been successfully applied to Sobolev spaces W1,p0(Ω), and then at least three solutions are obtained. Furthermore, Bonanno in [21] established the existence of two intervals of positive real parameters λ for which the functional ΦλJ has three critical points, and applied the result to obtain two critical points.

    By using three critical points theorem, Azroul et al. [22] discussed the fractional p-Laplace systems with bounded domain

    {M1([η]ps,p)(Δ)spη=λFη(x,η,ξ)+μGη(x,η,ξ),xΩ,M2([ξ]ps,p)(Δ)spξ=λFξ(x,η,ξ)+μGξ(x,η,ξ),xΩ,η=ξ=0,xRNΩ, (1.5)

    thus, the existence and multiplicity of solutions were obtained by Azroul et al. In addition, there are many scholars who have used different methods to study the existence of elliptic systems on bounded and unbounded regions, for instance, see [23,24,25] for details.

    With respect to the fractional p()-Laplace operators, Azroul et al. [26] dealt with the class of Kirchhoff type elliptic systems in nonlocal fractional Sobolev spaces with variable exponents and constant order

    {M(Is,p(x,y)(η))((Δ)sp(x,)η+|η|¯p(x)2η)=λFη(x,η,ξ)+μGη(x,η,ξ),inΩ,M(Is,q(x,y)(ξ))((Δ)sq(x,)ξ+|ξ|¯q(x)2ξ)=λFξ(x,η,ξ)+μGξ(x,η,ξ),inΩ,η=ξ=0,inRNΩ, (1.6)

    where

    Is,r(x,y)(ω)=Ω×Ω|ω(x)ω(y)|r(x,y)r(x,y)|xy|N+r(x,y)s     dxdy+Ω|ω(x)|¯r(x)¯r(x)dx.

    Based on the three critical points theorem introduced by Ricceri and on the theory of fractional Sobolev spaces with variable exponents, the existence of weak solutions for a nonlocal fractional elliptic system of (p(x,),q(x,))-Kirchhoff type with homogeneous Dirichlet boundary conditions was obtained. By using Ekeland's variational principle and dual fountain theorem, Bu et al. in [27] obtained some new existence and multiplicity of negative energy solutions for the fractional p()-Laplace operators with constant order without the Ambrosetti-Rabinowitz condition.

    Previous studies have shown that the fractional p()-Laplace operators with variable-order are much more complex and difficult than p-Laplace operators. The investigation of these problems has captured the attention of a host of scholars. For example, Wu et al. in [28] considered the fractional Kirchhoff systems with a bounded set Ω in RN, as follows:

    {M1(R2N|η(x)η(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)     dxdy)(Δ)s()p()η(x)=f(η,ξ)+a(x),xΩ,M2(R2N|ξ(x)ξ(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)     dxdy)(Δ)s()p()ξ(x)=g(η,ξ)+b(x),xΩ,η=ξ=0,xRNΩ, (1.7)

    by applying Ekeland variational principle, they obtained the existence of a solution for this class of problem.

    When μ=0, problem (1.1) reduces to the following fractional Kirchhoff-type elliptic systems

    {M1(RN×RN|η(x)η(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)     dxdy+Ω|η|¯p(x)¯p(x)dx)(Δs()p()η|η|¯p(x)η)=λFη(x,η,ξ),xΩ,M2(RN×RN|ξ(x)ξ(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)     dxdy+Ω|ξ|¯p(x)¯p(x)dx)(Δs()p()ξ|ξ|¯p(x)ξ)=λFξ(x,η,ξ),xΩ,η=ξ=0,xRNΩ, (1.8)

    Motivated by the above cited works, we take into account the nonlocal fractional Kirchhoff-type elliptic systems with variable-order. Our aims are to establish the existence of at least three solutions for problem (1.1) by utilizing Ricceri's result in [29] and obtain the existence of at least two solutions for problem (1.8) with the help of the multiple critical points theorem in [37]. The primary consideration of the paper is an extension of the results found in the literatures and our results are new to the Kirchhoff-type systems in some ways.

    For simplicity, Cj(j=1,2,...,N) are used in various places to denote distinct constants, i=1,2, and we denote

    C+(D):={H()C(D,R):1<H()H()H+()<+},

    where H() is a real-valued function and

    H():=minDH(),H+():=maxDH().

    F:Ω×R2R is a C1-function, whose partial derivatives are Fη,Fξ, which satisfy the following conditions.

    (F1): For some positive constant C, there exist α(x),β(x)C+(D) and 2+α++β+<p such that

    |Fs(x,s,t)|C|s|α(x)|t|β(x)+1,|Ft(x,s,t)|C|s|α(x)+1|t|β(x) for all(x,s,t)¯Ω×R×R;

    (F2): F(x,s,t)>0 for any(x,s,t)¯Ω×[1,+)×[1,+),andF(x,s,t)<0 for any(x,s,t)¯Ω×(0,1)×(0,1), F(x,0,0)=0 for a.e.x¯Ω.

    G:Ω×R2R is a C1-function, whose partial derivatives are Gη,Gξ, which satisfy assumptions, as follows:

    (G): sup|s|σ,|t|σ(|Gη(x,s,t)|+|Gξ(x,s,t)|)L1(Ω) for all σ>0.

    Definition 2. We say that (η,ξ)X0 is a (weak) solution of nonlocal Kirchhoff systems (1.1), if

    M1(RN×RN|η(x)η(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|η|¯p(x)¯p(x)dx)×(RN×RN|η(x)η(y)|p(x,y)2(η(x)η(y))(φ(x)φ(y))|xy|N+p(x,y)s(x,y)dxdy+Ω|η|¯p(x)ηφdx)+M2(RN×RN|ξ(x)ξ(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξ|¯p(x)¯p(x)dx)×(RN×RN|ξ(x)ξ(y)|p(x,y)2(ξ(x)ξ(y))(ψ(x)ψ(y))|xy|N+p(x,y)s(x,y)dxdy+Ω|ξ|¯p(x)ξψdx)=λΩFη(x,η,ξ)φdx+λΩFξ(x,η,ξ)ψdx+μΩGη(x,η,ξ)φdx+μΩGξ(x,η,ξ)ψdx, (1.9)

    for any (φ,ψ)X0, and we will introduce X0 in Section 2.

    Define the corresponding functional I:X0R associated with Kirchhoff systems (1.1), by

    I(η,ξ):=Φ(η,ξ)+λΨ(η,ξ)+μJ(η,ξ), (1.10)

    for all (η,ξ)X0, where

    J(η,ξ)=ΩG(x,η,ξ)dx,Ψ(η,ξ)=ΩF(x,η,ξ)dx, (1.11)
    Φ(η,ξ)=˜M1(δp()(η))+˜M2(δp()(ξ)), (1.12)
    δp()(υ)=RN×RN|υ(x)υ(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|υ|¯p(x)¯p(x)dx. (1.13)

    The functions Φ,Ψ,J: X0R are well defined, and we define their Gâteaux derivatives at (η,ξ)X0, by

    J(η,ξ),(φ,ψ)=ΩGη(x,η,ξ)φdxΩGξ(x,η,ξ)ψdx, (1.14)
    Ψ(η,ξ),(φ,ψ)=ΩFη(x,η,ξ)φdxΩFξ(x,η,ξ)ψdx, (1.15)
    Φ(η,ξ),(φ,ψ)=M1(δp()(η))η,φ+M2(δp()(ξ))ξ,ψ, (1.16)

    for all (φ,ψ)X0, where

    υ,ϕ=RN×RN|υ(x)υ(y)|p(x,y)2(υ(x)υ(y))(ϕ(x)ϕ(y))|xy|N+p(x,y)s(x,y)dxdy+Ω|υ|¯p(x)υϕdx, (1.17)

    for all (υ,ϕ)X0.

    Hence, (η,ξ)X0 is a (weak) solution of Kirchhoff systems (1.1) if and only if (η,ξ) is a critical point of the functional I, that is

    I(η,ξ)=Φ(η,ξ)+λΨ(η,ξ)+μJ(η,ξ)=0. (1.18)

    Definition 2. For s+p+<N, and denote by A: there exists a kind of functions F:Ω×R2R such that two Carathéodory functions Fη=Fη and Fξ=Fξ, satisfying

    sup(x,s,t)Ω×R×RFs(x,s,t)1+|s|ϑ(x)1+|t|ϑ(x)1<+andsup(x,s,t)Ω×R×RFt(x,s,t)1+|s|ϑ(x)1+|t|ϑ(x)1<+, (1.19)

    for any ϑ(x)[1,ps(x)).

    Now, let us show our results in this article.

    Theorem 1.1. For s(),p()C+(D) with s+p+<N and FA, assume that (S), (P), (M), (F1) and (F2) are satisfied. There exist three constants a,c1,c2>0 with 0<γ1<c1<c2 such that

    M+A(c1)<mH(c1,a),M+A(c2)<mH(c1,a). (1.20)

    Then, for any

    λ(M+γCϕH(c1,a),mγCϕmin{1A(c1),1A(c2)}),

    there exists a positive real number ρ such that the system (1.8) has at least two weak solutions wj=(ηj,ξj)X0(j=1,2) whose norms wj in X0 are less than some positive constant ρ.

    Theorem 1.2. For s(),p()C+(D) with s+p+<N and FA, assume that (S), (P), (M), (F1) and (F2) are satisfied. Then there exists an open interval Λ(0,+) and a positive real number ρ with the following property: For each λΛ and for two Carathéodory functions Gη, Gξ :Ω×R2R satisfying (G), there exists δ>0 such that for each μ[0,δ], problem (1.1) has at least three weak solutions wj=(ηj,ξj)X0(j=1,2,3) whose norms wj in X0 are less than some positive constant ρ.

    Remark 1.1. Existence results for the Kirchhoff-type elliptic systems with both boundary value problems and variational problems were obtained according to using critical points theorem by Ricceri and Bonanno, respectively, where the condition of Palais-Smale is not satisfied.

    Remark 1.2. The nonlocal Kirchhoff coefficient M1(t),M2(t) stir up some of the fundamental difficulties. To deal with these difficulties, we suppose that M1(t),M2(t) are strictly increasing functions, and then prove that the function Φ is a homeomorphism.

    The remaining of this article is organized as follows: Some fundamental results about the fractional Lebesgue spaces and Sobolve spaces are given in Section 2. In Section 3, in order to use critical point theory, we prove some technical lemmas. Theorem 1.1 and Theorem 1.2 are proved in Section 4. Finally, we make a conclusion in Section 5.

    To study Laplacian problems with variable exponents, we need to recall a slice of preliminary theories on generalized Lebesgue spaces Lϑ(x)(Ω) and give some necessary lemmas and propositions.

    For any ϑ(x)C+(D), the generalized Lebesgue spaces with variable exponents is defined by

    Lϑ(x)(Ω)={ξ|ξ:ΩRis a measurable function andΩ|ξ|ϑ(x)dx<}

    with respect to the norm

    ξϑ(x)=inf{χ>0:Ω|ξχ|ϑ(x)dx1},

    then, the spaces (Lϑ(x)(Ω),ϑ(x)) is a separable and reflexive Banach space, see [30,31].

    Let ϑ(x) be the conjugate exponent of ˜ϑ(x), namely

    1ϑ(x)+1˜ϑ(x)=1, for all xΩ.

    Lemma 2.1. (see [31]) Assume that ξLϑ(x)(Ω) and ηL˜ϑ(x)(Ω), then

    |Ωξηdx|(1ϑ+1˜ϑ)ξϑ(x)η˜ϑ(x)2ξϑ(x)η˜ϑ(x).

    Proposition 2.1. (see [30,32]) If we define

    ρϑ(x)(ξ)=Ω|ξ|ϑ(x)dx,

    then for all ξn,ξLϑ(x)(Ω), the following properties are possessed.

    (1)ξϑ(x)>1ξϑϑ(x)ρϑ(x)(ξ)ξϑ+ϑ(x),(2)ξϑ(x)<1ξϑ+ϑ(x)ρϑ(x)(ξ)ξϑϑ(x),(3)ξϑ(x)<1(resp.=1,>1)ρϑ(x)(ξ)<1(resp.=1,>1)),(4)ξϑ(x)0(resp.+)ρϑ(x)(ξn)0(resp.+),(5)limn|ξnξ|ϑ(x)=0limnρϑ(x)(ξnξ)=0.

    Remark 2.1. Note that for any function ϑ1(x),ϑ2(x)C+(¯Ω) and ϑ1(x)<ϑ2(x), there exists an embedding Lϑ2(x)(Ω)Lϑ1(x)(Ω) for any x¯Ω. Especially, when ϑ(x) constant, the results of Proposition 2.1 still hold.

    From now on, we briefly review a slice of essential lemmas and propositions about the Sobolev spaces, which will be used later. The readers are invited to consult [33,34,35] and the references therein.

    The fractional Sobolev spaces Ws(),p()(Ω) is defined as

    W=Ws(),p()(Ω):={ξL¯p()(Ω):Ω×Ω|ξ(x)ξ(y)|p(x,y)|xy|N+p(x,y)s(x,y)     dxdy<}

    and it can be endowed with the norm

    ξW=ξ¯p()+[ξ]s(),p() for allξW,

    where

    [ξ]s(),p()=inf{χ>0:Ω×Ω|ξ(x)ξ(y)|p(x,y)χp(x,y)|xy|N+p(x,y)s(x,y)     dxdy<1},

    then, the spaces (W,W) is a separable and reflexive Banach space, see [2,5] for a more detailed.

    We define the new fractional Sobolev spaces W concerning variable exponent and variable-order for some χ>0.

    W={ξ:RNR:ξ|ΩL¯p()(Ω):andQ|ξ(x)ξ(y)|p(x,y)χp(x,y)|xy|N+p(x,y)s(x,y)     dxdy<},

    where Q=R2N(Ωc×Ωc) and it can also be endowed with the norm

    ξW=ξ¯p()+[ξ]W for allξW,

    where

    [ξ]W=inf{χ>0:Q|ξ(x)ξ(y)|p(x,y)χp(x,y)|xy|N+p(x,y)s(x,y)     dxdy<1}.

    Remark 2.2. Note that the norm W is different from W for the reason that Ω×ΩQ and Ω×ΩQ.

    Let

    W0={ξW:ξ=0,a.e.inRNΩ}

    with the norm

    ξW0=inf{χ>0:Q|ξ(x)ξ(y)|p(x,y)χp(x,y)|xy|N+p(x,y)s(x,y)     dxdy=RN×RN|ξ(x)ξ(y)|p(x,y)χp(x,y)|xy|N+p(x,y)s(x,y)     dxdy<1}.

    W0 denotes the dual spaces of W0.

    In what follows, X0 will denote the Cartesian product of two fractional Sobolev spaces W0 and W0, i.e. X0=W0×W0. Defined the norm

    (η,ξ)X0=ηW0+ξW0,

    where ηW0,ξW0 is the norm of W0.

    Theorem 2.1. (see [2]) Let s(), p()C+(D) satisfy (S) and (P), with N>p(x,y)s(x,y) for all (x,y)¯ΩׯΩ. Let ϕ(x)C+(D) satisfy

    1<ϕ=minx¯Ωϕ(x)ϕ(x)<ps(x)=N¯p(x)N¯p(x)¯s(x),for any x¯Ω,

    where ¯p(x)=p(x,x) and ¯s(x)=s(x,x). Then, there exists a constant Cϕ=Cϕ(N,s,p,ϕ,Ω)>0 such that

    ξϕ()CϕξW0

    for any ξW0. Moreover, the embedding W0Lϕ()(Ω) is compact.

    Proposition 2.2. (see [36]) If we define

    ρs()p()(ξ)=RN×RN|ξ(x)ξ(y)|p(x,y)|xy|N+p(x,y)s(x,y)     dxdy+Ω|ξ|¯p(x)dx,

    then for all ξn,ξW0, the following properties hold.

    (1)ξW0<1(resp.=1,>1)ρs()p()(ξ)<1(resp.=1,>1),(2)ξW0<1ξp+W0ρs()p()(ξ)ξpW0,(3)ξW0>1ξpW0ρs()p()(ξ)ξp+W0,(4)limnξnW0=0(resp.+)limnρs()p()(ξn)=0(resp.+),(5)limnξnξW0=0limnρs()p()(ξnξ)=0.

    In order to prove our main results, we present two multiple critical theorems: The first ensures the existence of two critical points, while the second establishes the existence of three critical points due to Ricceri.

    Theorem 2.2. (see [37]) Let X be a reflexive real Banach space, and let Φ,Ψ:XR be two sequentially weakly lower semicontinuous functions. Assume that Φ is (strongly) continuous and satisfies limηΦ(η)=+. Assume also that there exist two constants r1 and r2 such that

    (i)infXΦ<r1<r2,(ii)φ1(r1)<φ2(r1,r2),(iii)φ1(r2)<φ2(r1,r2),

    where

    φ1(ri)=infηΦ1 (,ri)Ψ(η)infη¯Φ1 (,ri)ωΨ(η)riΦ(η),φ2(r1,r2)=infηΦ1 (,r1)supξΦ1 [r1,r2)Ψ(η)Ψ(ξ)Φ(ξ)Φ(η),

    for i=1,2 and ¯Φ1(,ri)ω is the closure of Φ1(,ri) in the weak topology. Then, for each

    λ(1φ2(r1,r2),min{1φ1(r1),1φ2(r2)}),

    the functional Φ+λΨ has two local minima which lie in Φ1(,r1) and Φ1[r1,r2), respectively.

    Theorem 2.3. (see [29]) Let X be a reflexive real Banach space Φ:XR is a continuously Gâteaux differentiable and sequentially weakly lower semicontinuous functional whose Gâteaux derivative admits a continuous inverse on X and Φ is bounded on each bounded subset of X; Ψ:XR is a continuously Gâteaux differentiable functional whose Gâteaux derivative is compact, IR an interval. Assume that

    limη(Φ(η)+λΨ(η))=+,

    for λI, and that there exists rR and η0,η1X such that

    Φ(η0)<r<Φ(η1),infηΦ1((,r])Ψ(η)>(Φ(η1)r)Ψ(η0)+(rΦ(η0))Ψ(η1)Φ(η1)Φ(η0).

    Then there exists an open interval ΛI and a positive real number ρ with the following property: For every λΛ and every C1-functional J:XR with compact derivative, there exists δ>0 such that for each μ[0,δ] the equation

    Φ(η)+λΨ(η)+μJ(η)=0

    has at least three solutions in X whose norms are less than δ.

    In this section, in order to use critical point theory for Kirchhoff systems (1.1), we need the following crucial lemmas, which will play an important role in the proof of our results.

    Lemma 3.1. Assume that the functions G,FA, then J,ΨC1(X0,R) and their derivatives are defined as (1.14) and (1.15) for all (φ,ψ)X0. Moreover J,Ψ:X0R is compact.

    Proof. Here we follow the approach in [22], for completeness, we give the proof process. Suppose that s+p+<N, we prove that J(η,ξ) is continuous operator. Set {(ηn,ξn)}X0 with (ηn,ξn)(η,ξ) strongly in X0. According to Theorem 2.1, we obtain

    (ηn,ξn)(η,ξ) in Lϑ(x)(Ω)×Lϑ(x)(Ω).

    So, for a subsequence denoted by {(ηn,ξn)}, there exist functions ˜η,˜ξLϑ(x)(Ω) such that

    ηnη,ξnξa.e. in Ω,
    |ηn|˜η,|ξn|˜ξa.e. in Ω,

    for all nN.

    Fix (˜η,˜ξ)X0 with (˜η,˜ξ)X01. Since GA, combining Lemma 2.1 with Theorem 2.1, we obtain

    |J(ηn,ξn)J(η,ξ),(˜η,˜ξ)||Ω[Gη(x,ηn,ξn)dxGη(x,η,ξ)]˜ηdx+Ω[Gξ(x,ηn,ξn)dxGξ(x,η,ξ)]˜ξdx|C3Gη(x,ηn,ξn)Gη(x,η,ξ)Lϑ(x)(Ω)˜ηLϑ(x)(Ω)+C3Gξ(x,ηn,ξn)Gξ(x,η,ξ)Lϑ(x)(Ω)˜ξLϑ(x)(Ω)C4Gη(x,ηn,ξn)Gη(x,η,ξ)Lϑ(x)(Ω)˜ηW0+C4Gξ(x,ηn,ξn)Gξ(x,η,ξ)Lϑ(x)(Ω)˜ξW0.C4(Gη(x,ηn,ξn)Gη(x,η,ξ)Lϑ(x)(Ω)+Gξ(x,ηn,ξn)Gξ(x,η,ξ)Lϑ(x)(Ω))(˜η,˜ξ)X0.

    Consequently, for (˜η,˜ξ)X01, we get

    J(ηn,ξn)J(η,ξ),(˜η,˜ξ)W0C4Gη(x,ηn,ξn)Gξ(x,η,ξ)Lϑ(x)(Ω)+C4Gξ(x,ηn,ξn)Gξ(x,η,ξ)Lϑ(x)(Ω).

    According to Definition 2, we deduce

    Gη(x,ηn,ξn)Gη(x,η,ξ)0asna.e.xΩ,Gξ(x,ηn,ξn)Gξ(x,η,ξ)0asna.e.xΩ

    and

    |Gη(x,ηn,ξn)|C5(1+|¯η|ϑ(x)1+|¯ξ|ϑ(x)1),|Gξ(x,ηn,ξn)|C5(1+|¯η|ϑ(x)1+|¯ξ|ϑ(x)1),

    using the dominate convergence theorem, we have

    Gη(x,ηn,ξn)Gη(x,η,ξ))Lϑ(x)(Ω)0asna.e.xΩ,Gξ(x,ηn,ξn)Gξ(x,η,ξ))Lϑ(x)(Ω)0asna.e.xΩ.

    So, this prove that J(η,ξ) is continuous operator.

    We show the operator J(η,ξ) is compact. Let {(ηn,ξn)} be a bounded subsequence in X0. Arguing in the same way as above, we obtain that the sequence {J(ηn,ξn)} converges strongly, therefor the operator J(η,ξ) is compact.

    Similarly, we also deduce that Ψ(η,ξ) is continuous and compact.

    Lemma 3.2. Assume that (M) is satisfied. Then

    (i) & Φ is sequentially weakly lower semicontinuous and bounded on each bounded subset,

    (ii) & Φ:X0X0 is a strictly monotone and continuous operator,

    (iii) & Φ:X0X0 is a homeomorphism.

    Proof. (i) Since ˜Mi(t)>mi>0, ˜Mi(t) are increasing function on [0,+). Argue in a similar way from [38,Lemma 2.4], the operators ηδp()(η) and ξδp()(ξ) are strictly monotone. From [39,Proposition 25.10], δp()(η) and δp()(ξ) are strictly convex. Set {(ηn,ξn)}X0 be a subsequence such that

    ηnη,ξnξinW0.

    Based on the convexity of δp()(η) and δp()(ξ), we obtain

    δp()(ηn)δp()(η)δp()(η),ηnη (3.1)

    and

    δp()(ξn)δp()(ξ)δp()(ξ),ξnξ. (3.2)

    Thus, we have

    δp()(η)lim infnδp()(ηn) (3.3)

    and

    δp()(ξ)lim infnδp()(ξn), (3.4)

    namely, the operators ηδp()(η) and ξδp()(ξ) are sequentially weakly lower semicontinuous.

    On the other side, since t˜Mi(t) are continuous and monotonous functions, we obtain

    lim infnΦ(ηn,ξn)=lim infn˜M1(δp()(ηn))+lim infn˜M2(δp()(ξn))˜M1(lim infnδp()(ηn))+˜M2(lim infnδp()(ξn))˜M1(δp()(η))+˜M2(δp()(ξ))=Φ(η,ξ). (3.5)

    Consequently, the operator Φ is sequentially weakly lower semicontinuous.

    Right now, we claim that Φ is bounded on each bounded subset of X0. Set {(ηn,ξn)}X0 be a bounded subsequence. By Proposition 2.2, there exist constants C6,C7>0 so that

    δp()(η)C6andδp()(ξ)C7.

    Since ˜M1 and ˜M2 is monotone, we have

    Φ(ηn,ξn)=˜M1(δp()(ηn))+˜M2(δp()(ξn))˜M1(C6)+˜M2(C7). (3.6)

    Hence, the operator Φ is bounded.

    (ii) Let (M) is satisfied, then ΦC1(X0,R), and its derivatives are defined by (1.16), we have

    Φ(η,ξ),(φ,ψ)=Φη(η),φ+Φξ(ξ),ψ, (3.7)

    where

    Φη(η),φ=M1(δp()(η))η,φand Φξ(ξ),ψ=M2(δp()(ξ))ξ,ψ, (3.8)

    for all (η,ξ),(φ,ψ)X0. Therefore, Φ(η,ξ)X0, where X0 denotes the dual spaces of X0.

    In the first place, we show that Φ:X0X0 is a strictly monotone operator. Since Φη(η) and Φξ(ξ) are strictly monotone operators (see [40], Theorem 2.1), Φ is a strictly monotone operator.

    Next, we claim that Φ:X0X0 is a continuous operator. Set {(ηn,ξn)}X0 be a sequence, which converges strongly to (η,ξ) in X0. So, for a subsequence denoted by {(ηn,ξn)}, we suppose

    ηnη,ξnξa.e. in Ω.

    Then, the sequence

    {|ηn(x)ηn(y)|p(x,y)2(ηn(x)ηn(y))|xy|(N+p(x,y)s(x,y))/p(x,y)} (3.9)

    and

    {|ξn(x)ξn(y)|p(x,y)2(ξn(x)ξn(y))|xy|(N+p(x,y)s(x,y))/p(x,y)} (3.10)

    are bounded in Lp(RN×RN), and we have

    |ηn(x)ηn(y)|p(x,y)2(ηn(x)ηn(y))|xy|(N+p(x,y)s(x,y))/p(x,y)|η(x)η(y)|p(x,y)2(η(x)η(y))|xy|(N+p(x,y)s(x,y))/p(x,y) (3.11)

    and

    |ξn(x)ξn(y)|p(x,y)2(ξn(x)ξn(y))|xy|(N+p(x,y)s(x,y))/p(x,y)|ξ(x)ξ(y)|p(x,y)2(ξ(x)ξ(y))|xy|(N+p(x,y)s(x,y))/p(x,y). (3.12)

    According to the Brezis-Lieb lemma [41], combining (3.11) with (3.12) implies

    RN×RN||ηn(x)ηn(y)|p(x,y)2(ηn(x)ηn(y))|xy|(N+p(x,y)s(x,y))/p(x,y)|η(x)η(y)|p(x,y)2(η(x)η(y))|xy|(N+p(x,y)s(x,y))/p(x,y)|p(x,y)dxdy=limnRN×RN(|ηn(x)ηn(y)|p(x,y)|xy|N+p(x,y)s(x,y)|η(x)η(y)|p(x,y)|xy|N+p(x,y)s(x,y))dxdy (3.13)

    and

    RN×RN||ξn(x)ξn(y)|p(x,y)2(ξn(x)ξn(y))|xy|(N+p(x,y)s(x,y))/p(x,y)|ξ(x)ξ(y)|p(x,y)2(ξ(x)ξ(y))|xy|(N+p(x,y)s(x,y))/p(x,y)|p(x,y)dxdy=limnRN×RN(|ξn(x)ξn(y)|p(x,y)|xy|N+p(x,y)s(x,y)|ξ(x)ξ(y)|p(x,y)|xy|N+p(x,y)s(x,y))dxdy. (3.14)

    Based on the fact that ηnη, ξnξ strongly in X0 yields

    limnRN×RN(|ηn(x)ηn(y)|p(x,y)2(ηn(x)ηn(y))|xy|(N+p(x,y)s(x,y))/p(x,y)|η(x)η(y)|p(x,y)2(η(x)η(y))|xy|(N+p(x,y)s(x,y))/p(x,y))dxdy=0 (3.15)

    and

    limnRN×RN(|ξn(x)ξn(y)|p(x,y)2(ξn(x)ξn(y))|xy|(N+p(x,y)s(x,y))/p(x,y)|ξ(x)ξ(y)|p(x,y)2(ξ(x)ξ(y))|xy|(N+p(x,y)s(x,y))/p(x,y))dxdy=0. (3.16)

    Furthermore, M1(t) and M2(t) are continuous functions, which imply

    M1(δp()(ηn))M1(δp()(η)) (3.17)

    and

    M2(δp()(ξn))M2(δp()(ξ)). (3.18)

    With the help of Lemma 2.1, we get

    |Ω|ηn|¯p(x)2ηn(ηnη)dx|Ω|ηn|¯p(x)1|ηnη|dx||ηn|¯p(x)1|¯p(x)¯p(x)1|ηnη|¯p(x)0,

    as n, and then

    limnΩ|ηn|¯p(x)2ηn(ηnη)dx=0. (3.19)

    Arguing in the similar way as above, we have

    limnΩ|ξn|¯p(x)2ξn(ξnξ)dx=0. (3.20)

    Thus, from (3.15)–(3.20), we obtain

    Φ(ηn,ξn)Φ(η,ξ)X0=sup(φ,ψ)X0,(φ,ψ)X01|Φ(ηn,ξn),(φ,ψ)Φ(η,ξ),(φ,ψ)|0. (3.21)

    (iii) Since the operator Φ is a strictly monotone, it follows that Φ is injective. Set (η,ξ)X0 be such that (η,ξ)X0>1, from (M) and Proposition 2.2, we have

    Φ(η,ξ),(η,ξ)(η,ξ)X0=M1(δp()(η))η,η+M2(δp()(ξ))ξ,ξ(η,ξ)X0m1ρs()p()(η)+m2ρs()p()(ξ)(η,ξ)X0min(m1,m2)min(ηp+W0,ηpW0)+min(ξp+W0,ξpW0)ηW0+ξW0. (3.22)

    Thus, (3.22) implies that

    lim(η,ξ)Φ(η,ξ),(η,ξ)(η,ξ)X0=+. (3.23)

    Consequently, Φ is coercive operator, thanks to the Minty-Browder Theorem (see [39], Theorem 26A), Φ is a surjection. Due to its monotonicity, Φ is an injection. So, (Φ)1 exists.

    Let us first prove that Φ satisfies property:

    (Q):if(ηn,ξn)(η,ξ)andΦ(ηn,ξn)Φ(η,ξ),then(ηn,ξn)(η,ξ).

    Indeed, set (ηn,ξn)(η,ξ) in X0, Φ(ηn,ξn)Φ(η,ξ) in X0 and Theorem 2.1, which implies

    ηnη,ξnξa.e.xΩ.

    Based on Fatou's Lemma yields

    lim infn(RN×RN|ηn(x)ηn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ηn|¯p(x)dx)RN×RN|η(x)η(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|η|¯p(x)dx (3.24)

    and

    lim infn(RN×RN|ξn(x)ξn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξn|¯p(x)dx)RN×RN|ξ(x)ξ(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξ|¯p(x)dx. (3.25)

    Applying Young's inequality, there exist C8,C9>0 such that

    on=Φηn,ηnη+Φξn,ξnξ=M1(δp()(ηn))(RN×RN|ηn(x)ηn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ηn|¯p(x)dxRN×RN|ηn(x)ηn(y)|p(x,y)2(ηn(x)ηn(y))(η(x)η(y))|xy|N+p(x,y)s(x,y)dxdyΩ|ηn|¯p(x)2ηnηdx)+M2(δp()(ξn))(RN×RN|ξn(x)ξn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξn|¯p(x)dxRN×RN|ξn(x)ξn(y)|p(x,y)2(ξn(x)ξn(y))(ξ(x)ξ(y))|xy|N+p(x,y)s(x,y)dxdyΩ|ξn|¯p(x)2ξnξdx)m1(RN×RN|ηn(x)ηn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ηn|¯p(x)dxRN×RN|ηn(x)ηn(y)|p(x,y)1(η(x)η(y))|xy|N+p(x,y)s(x,y)dxdyΩ|ηn|¯p(x)1ηdx)+m2(RN×RN|ξn(x)ξn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξn|¯p(x)dxRN×RN|ξn(x)ξn(y)|p(x,y)1(ξ(x)ξ(y))|xy|N+p(x,y)s(x,y)dxdyΩ|ξn|¯p(x)1ξdx)m1(C8RN×RN|ηn(x)ηn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+C8Ω|ηn|¯p(x)dxC9RN×RN|η(x)η(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdyC9Ω|η|¯p(x)dx)+m2(C8RN×RN|ξn(x)ξn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+C8Ω|ξn|¯p(x)dxC9RN×RN|ξ(x)ξ(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdyC9Ω|ξ|¯p(x)dx). (3.26)

    The passage to the liminf implies in the above inequality, we get

    0m1[lim infn(RN×RN|ηn(x)ηn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ηn|¯p(x)dx)RN×RN|η(x)η(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdyΩ|η|¯p(x)dx]+m2[lim infn(RN×RN|ξn(x)ξn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξn|¯p(x)dx)RN×RN|ξ(x)ξ(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdyΩ|ξ|¯p(x)dx]. (3.27)

    This and combining (3.24) with (3.25), we obtain

    lim infn(RN×RN|ηn(x)ηn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ηn|¯p(x)dx)=RN×RN|η(x)η(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|η|¯p(x)dx (3.28)

    and

    lim infn(RN×RN|ξn(x)ξn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξn|¯p(x)dx)=RN×RN|ξ(x)ξ(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξ|¯p(x)dx. (3.29)

    Then, for a subsequence, we have

    limn(RN×RN|ηn(x)ηn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ηn|¯p(x)dx)=RN×RN|η(x)η(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|η|¯p(x)dx (3.30)

    and

    limn(RN×RN|ξn(x)ξn(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξn|¯p(x)dx)=RN×RN|ξ(x)ξ(y)|p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξ|¯p(x)dx, (3.31)

    that is

    limnρs()p()(ηn)=ρs()p()(η),limnρs()p()(ξn)=ρs()p()(ξ). (3.32)

    Since ηnη, ξnξ in X0, Proposition 2.2 imply that

    limnρs()p()(ηnη)=0,limnρs()p()(ξnξ)=0, (3.33)

    that is

    limnηnηW0=0,limnξnξW0=0. (3.34)

    Consequently

    (ηn,ξn)(η,ξ))=ηnηW0+ξnξW00. (3.35)

    Next, we prove that (Φ)1 is continuous operator. Let {(fn,gn)}X0 such that fnf0 in X0 and gng0 in X0. Set (ηn,ξn),(η,ξ)X0 such that

    Φ1η(fn)=ηn,Φ1η(f0)=η,Φ1ξ(gn)=ξn,Φ1ξ(g0)=ξ. (3.36)

    According to coercivity of Φ, we conclude that (ηn,ξn) is bounded. Then, up to subsequence (ηn,ξn)(~η0,~ξ0), which implies

    limnΦ(ηn,ξn)Φ(η,ξ),(ηn,ξn)(~η0,~ξ0)=limn(fn,gn)(f0,g0),(ηn,ξn)(~η0,~ξ0)=0. (3.37)

    Combining the property of (Q) and the continuity of Φ, we obtain

    (ηn,ξn)(~η0,~ξ0)inX0,Φ(ηn,ξn)Φ(~η0,~ξ0)=Φ(η,ξ)inX0. (3.38)

    Since Φ is an injection, we deduce that (η,ξ)=(~η0,~ξ0).

    In this subsection, we firstly prove Theorem 1.1 by applying Theorem 2.2.

    Proof of Theorem 1.1. Let X=X0, Ψ and Φ are given as (1.11) and (1.12), respectively. Note that Ψ is a compact derivative from Lemma 3.1, Lemma 3.2 ensures that Φ is a weakly lower semicontinuous and bounded operator in X0, and Φ admits a continuous inverse operator Φ:X0X0.

    Let (η,ξ)=(0,0), then Φ(0,0)=Ψ(0,0)=0. There exist a point x0Ω and pick two positive constants R2,R1(R2>R1) such that B(x0,R2)Ω. Set a,c be positive constants and define the function ω(x) by

    ω(x)={0,x¯ΩB(x0,R2),aR2R1{R2[ΣNi=1(xixi0)]1/2},xB(x0,R2)B(x0,R1),a,xB(x0,R1), (4.1)

    where B(x0,R) stands for the open ball in RN of radius R centered at x0, then (ω(x),ω(x))X0. Denote

    A(c)=Ωsup(s,t)R×R:|s|p()+|t|p()cF(x,s,t),K(a)=B(x0,R2)B(x0,R1)F(x,ω(x),ω(x))+B(x0,R1)F(x,a,a),H(c,a)=K(a)A(c),m=min{m1p+,m2p+},M+=max{M1p,M2p}. (4.2)

    Under condition (M) and by a simple computation, we obtain

    m(min{η(x)p+W0+ξ(x)p+W0,η(x)pW0+ξ(x)pW0})Φ(η,ξ)M+(max{η(x)p+W0+ξ(x)p+W0,η(x)pW0+ξ(x)pW0}). (4.3)

    When ηW0,ξW0 bounded (ξW0,ηW0 bounded) and ηW0,ξW0, we have

    (η,ξ)X0,

    this implies that

    lim(η,ξ)X0Φ(η,ξ)=+. (4.4)

    Fix γ such that 0<γ1<c1<c2 and set

    r1=mc1γCϕ,r2=M+c2γCϕ. (4.5)

    By virtue of (4.1), set η0=ξ0=ω for (η0,ξ0)X0 and

    η0W0=ωW0,ξ0W0=ωW0. (4.6)

    Consequently, we have

    r1<Φ(η0,ξ0)<r2. (4.7)

    Thus, (4.7) implies that

    φ2(r1,r2)=inf(η,ξ)Φ1(,r1)sup(η1,ξ1)Φ1[r1,r2)Ψ(η,ξ)Ψ(η1,ξ1)Φ(η1,ξ1)Φ(η,ξ)inf(η,ξ)Φ1(,r1)Ψ(η,ξ)Ψ(η0,ξ0)Φ(η0,ξ0)Φ(η,ξ). (4.8)

    According to (F2), (1.20) and (4.1), we obtain

    ΩF(x,η0,ξ0)dx=K(a)>H(c1,a)>M+mA(c1)>A(c1)=Ωsup(s,t)R×R:|s|p()+|t|p()c1F(x,s,t)dx. (4.9)

    For each (η,ξ)X0 with Φ(η,ξ)r1, we conclude

    |η(x)|p()+|ξ(x)|p()Cϕ(min{η(x)p+W0+ξ(x)p+W0,η(x)pW0+ξ(x)pW0})Cϕr1m=c1γc1, for allxΩ. (4.10)

    Fix γ0 such that 0<max{η0p+W0+ξ0p+W0,η0pW0+ξ0pW0}<γ0<γ. Thus, combining with (4.9) and (4.10), we have

    Ψ(η,ξ)Ψ(η0,ξ0)Φ(η0,ξ0)Φ(η,ξ)=ΩF(x,η0,ξ0)dxΩF(x,η,ξ)dxΦ(η0,ξ0)Φ(η,ξ)ΩF(x,η0,ξ0)dxΩsup(η,ξ)R×R:|η|p()+|ξ|p()c1F(x,η,ξ)dxΦ(η0,ξ0)Φ(η,ξ)ΩF(x,η0,ξ0)dxΩsup(η,ξ)R×R:|η|p()+|ξ|p()c1F(x,η,ξ)dxΦ(η0,ξ0)ΩF(x,η0,ξ0)dxΩsup(η,ξ)R×R:|η|p()+|ξ|p()c1F(x,η,ξ)dxM+(η0p()W0+ξ0p()W0)>CϕM+γH(c1,a). (4.11)

    From (4.8) and (4.11), we deduce

    φ2(r1,r2)>CϕM+γH(c1,a). (4.12)

    Similarly, for each (η,ξ)X0 such that Φ(η,ξ)r, we get

    |η(x)|p()+|ξ(x)|p()Cϕ(min{η(x)p+W0+ξ(x)p+W0,η(x)pW0+ξ(x)pW0})Cϕrm, for allxΩ. (4.13)

    According to Φ being sequentially weakly lower semicontinuous, then ¯Φ1(,r)ω=Φ1(,r]. Consequently, we have

    inf(η,ξ)Φ1(,r)Ψ(η,ξ)inf(η,ξ)¯Φ1(,r)ωΨ(η,ξ)rΦ(η,ξ)Ψ(0,0)inf(η,ξ)¯Φ1(,r)ωΨ(η,ξ)rΦ(0,0)inf(η,ξ)¯Φ1(,r)ωΨ(η,ξ)rΩsup(η,ξ)R×R:|η|p()+|ξ|p()CϕrmF(x,η,ξ)dxr, (4.14)

    this implies

    φi(ri)A(ci)ri=CϕmciγA(ci)<CϕM+γH(c1,a). (4.15)

    Thus, by virtue of (4.12) and (4.15), we obtain

    φi(ri)<CϕM+γH(c1,a)<φ2(r1,r2). (4.16)

    Therefore, the conditions of Theorem 2.2 are satisfied. Thus, the functional Φ+λΨ has two local minima (ηi,ξi)X0, which lie in Φ1(,ri), respectively. Since I=Φ+λΨC1(X0,R),(ηi,ξi)X0 are the two solutions of the following equation

    Φ(ηi,ξi)+λΨ(ηi,ξi)=0. (4.17)

    That is, (ηi,ξi)X0 are two weak solutions of the nonlocal fractional Kirchhoff-type elliptic systems (1.8). Since Φ(ηi,ξi)<ri, combining Theorem 2.1, we obtain

    |ηi(x)|p()+|ξi(x)|p()Cϕ(min{ηi(x)p+W0+ξi(x)p+W0,ηi(x)pW0+ξi(x)pW0})Cϕr2mc2, for allxΩ, (4.18)

    which implies that there is a real number ρ>0 so that wj=(ηj,ξj)X0 and wj<ρ.

    In what follows, we prove Theorem 1.2 Theorem by using Theorem 2.3.

    Proof of Theorem 1.2. Let X=X0, Ψ and Φ are given by (1.11) and (1.12), respectively. Note that Ψ is a compact derivative from Lemma 3.1, Lemma 3.2 ensures that Φ is a weakly lower semicontinuous and bounded operator in X0, and Φ admits a continuous inverse operator Φ:X0X0. Moreover,

    lim(η,ξ)x0Φ(η,ξ)+λΨ(η,ξ)=+ (4.19)

    for any λ(0,+). Indeed,

    Φ(η,ξ)=˜M1(RN×RN|η(x)η(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|η|¯p(x)¯p(x)dx)+˜M2(RN×RN|ξ(x)ξ(y)|p(x,y)p(x,y)|xy|N+p(x,y)s(x,y)dxdy+Ω|ξ|¯p(x)¯p(x)dx)=˜M1(δp()(η))+˜M2(δp()(ξ))m1p+ρs()p()(η)+m2p+ρs()p()(ξ)m1p+min(ηpW0,ηp+W0)+m2p+min(ξpW0,ξp+W0). (4.20)

    By virtue of (F1), we obtain

    |F(x,η,ξ)|C10|η|α(x)+1|ξ|β(x)+1 for all(x,η,ξ)¯Ω×R×R, (4.21)

    consequently,

    Ψ(η,ξ)=ΩF(x,η,ξ)dxC10Ω|η|α(x)+1|ξ|β(x)+1dxC10|Ω|max(η1+α,η1+α+)max(ξ1+β,ξ1+β+). (4.22)

    Since W0 are continuously embedded in C(¯Ω), there exists a constant C11 such that

    Ψ(η,ξ)C10C11|Ω|max(η1+αW0,η1+α+W0)max(ξ1+βW0,ξ1+β+W0). (4.23)

    When ηW0,ξW0 bounded (ξW0,ηW0 bounded) and ηW0,ξW0, we have

    (η,ξ)X0. (4.24)

    According to 2+α++β+<p, there exist p1<p such that 1+α+p1+1+β+p1=1. Therefore, from (4.23) and Young's inequality, we deduce

    Ψ(η,ξ)C10C11|Ω|(1+α+p1ηp1W0+1+β+p1max(ηp1W0,1)). (4.25)

    Hence, for λ>0 and p1<p, the combination of (4.20), (4.24) and (4.25) implies that

    lim(η,ξ)X0Φ(η,ξ)+λΨ(η,ξ)=+. (4.26)

    In view of (F2), we choose δ>1 such that

    F(x,s,t)>0 (4.27)

    for alls,t>δ,xΩ, and then

    F(x,s,t)>0=F(x,0,0)>F(x,τ1,τ2) (4.28)

    for allτ1,τ2[0,1). Set a0,b1 be two positive real numbers such that a0<min(1,Cϕ) with Cϕ defined by Theorem 2.1 and b1>δ with bp1|Ω|>1. From (4.28), we have

    ΩF(x,b1,b1)>0=F(x,0,0)>Ωsup0s,ta0F(x,s,t). (4.29)

    Let

    r:=min(m1p+(a0Cϕ)p+,m2p+(a0Cϕ)p+). (4.30)

    Choosing (η0(x),ξ0(x))=(0,0),(η1(x),ξ1(x))=(b1,b1) for allxΩ, such that we obtain

    Φ(η0,ξ0)=Ψ(η0,ξ0)=(0,0) (4.31)

    and

    Φ(η1,ξ1)=˜M1(Ω|b1|¯p(x)¯p(x)dx)+˜M2(Ω|b1|¯p(x)¯p(x)dx)m1(Ω|b1|¯p(x)¯p(x)dx)+m2(Ω|b1|¯p(x)¯p(x)dx)|Ω|(m1bp1p++m2bp1p+)m1p++m2p+r. (4.32)

    Thus, the combination of (4.31) and (4.32) implies that

    Φ(η0,ξ0)<r<Φ(η1,ξ1). (4.33)

    On the other hand

    (Φ(η1,ξ1)r)Ψ(η0,ξ0)+(rΦ(η0,ξ0))Ψ(η1,ξ1)Φ(η1,ξ1)Ψ(η0,ξ0)=rΩF(x,b1,b1)dx˜M1(Ω|b1|¯p(x)¯p(x)dx)+˜M2(Ω|b1|¯p(x)¯p(x)dx). (4.34)

    Let (η,ξ)X0 such that Φ(η,ξ)r, we get

    Φ(η,v)=˜M1(δp()(η))+˜M2(δp()(ξ))min{m1,m2}(1p+ρs()p()(η)+1p+ρs()p()(ξ)), (4.35)

    which implies that

    ρs()p()(η)rp+m1<1,ρs()p()(ξ)rp+m2<1. (4.36)

    According to Proposition 2.2, we derive

    ηw01,ξw01. (4.37)

    Combining with (4.36), we deduce

    ηw0(rp+m1)1p+,ξw0(rp+m2)1p+. (4.38)

    Therefore

    |η(x)|Cϕ(rp+m1)1p+a0,|ξ(x)|Cϕ(rp+m2)1p+a0 (4.39)

    for all xΩ. It follows from (4.34) that

    inf(η,ξ)Φ1((,r])Ψ(η,ξ)=supΦ(η,ξ)r(Ψ(η,ξ))sup{(η,ξ)X0:|η(x)|,|ξ(x)|a0,xΩ}ΩF(x,η,ξ)dxΩsup0s,taF(x,s,t)dx0(Φ(η1,ξ1)r)Ψ(η0,ξ0)+(rΦ(η0,ξ0))Ψ(η1,ξ1)Φ(η1,ξ1)Ψ(η0,ξ0). (4.40)

    Consequently, the conditions of Theorem 2.3 are satisfied. For every compact interval Λ(0,+) and GA, we fix λΛ and put

    J(η,ξ)=ΩG(x,η,ξ)dx,

    for all (η,ξ)X0. Then, J is a compact derivative. Therefore, there exists a positive number δ such that for every μ[0,δ], (ηi,ξi)X0 are three solutions of the following equation

    Φ(ηi,ξi)+λΨ(ηi,ξi)+μJ(ηi,ξi)=0.

    That is, (ηi,ξi)X0 are three weak solutions of the nonlocal fractional Kirchhoff-type elliptic systems (1.1).

    In this article, we consider a kind of p()-fractional Kirchhoff systems. Under some reasonable assumptions, we obtain two and three solutions based on Bonanno's multiple critical points theorems and Ricceri's three critical points theorem, where the condition of Palais-Smale is not requested. Several recent results of literatures are extended and improved.

    This work is supported by the Postgraduate Research & Practice Innovation Program of Jiangsu Province (KYCX21$ _{-} 54), the Natural Science Foundation of Jiangsu Province (BK20180500), National Key Research and Development Program of China (2018YFC1508100), Special Soft Science Project of Technological Innovation in Hubei Province (2019ADC146), and Natural Science Foundation of China (11701595).

    The authors declare that they have no competing interests.



    [1] J. B. Zuo, A. Fiscella, A. Bahrouni, Existence and multiplicity results for p()&q() fractional Choquard problems with variable order, Complex Var. Elliptic Equ., (2020), 1–17.
    [2] J. B. Zuo, T. Q. An, A. Fiscella, A critical Kirchhoff-type problem driven by a p()-fractional Laplace operator with variable s()-order, Math. Methods Appl. Sci., 44 (2020), 1071–1085.
    [3] J. B. Zuo, L. B. Yang, S. H. Liang, A variable-order fractional p()-Kirchhoff type problemin RN, Math. Methods Appl. Sci., 44 (2020), 3872–3889
    [4] Y. Cheng, B. Ge, R. Agarwal, Variable-order fractional sobolev spaces and nonlinear elliptic equations with variable exponents, J. Math. Phy., 61 (2020), 071507. doi: 10.1063/5.0004341
    [5] R. Biswas, S. Tiwari, Variable order nonlocal Choquard problem with variable exponents, Complex Var. Elliptic Equ., (2020), 853–875.
    [6] P. Pucci, M. Q. Xiang, B. L. Zhang, Multiple solutions for nonhomogeneous Schrödinger-Kirchhoff type equations involving the fractional p-Laplacian in RN, Calc. Var. Partial Dif. Equ., 54 (2015), 2785–2806. doi: 10.1007/s00526-015-0883-5
    [7] M. Q. Xiang, B. L. Zhang, Degenerate Kirchhoff problems involving the fractional p-Laplacian without the (AR) condition, Complex Var. Elliptic Equ., 60 (2015), 1277–1287. doi: 10.1080/17476933.2015.1005612
    [8] E. Di Nezz, G. Palatucci, E. Valdinoci, Hitchhiker's guide to the fractional Sobolev spaces, Bull. Math. Sci., 136 (2012), 521–573. doi: 10.1016/j.bulsci.2011.12.004
    [9] G. Kirchhoff, Vorlesungen über mathematische Physik: Mechanik, Leipzig: Druck Und Verlag Von B. G. Teubner, 1876.
    [10] G. M. Bisci, L. Vilasi, On a fractional degenerate Kirchhoff-type problem, Commun. Contemp. Math., 19 (2017), 1550088. doi: 10.1142/S0219199715500881
    [11] C. E. T. Ledesma, Multiplicity result for non-homogeneous fractional Schrödinger-Kirchhoff-type equations in RN, Adv. Nonlinear Anal., 7 (2018), 247–257. doi: 10.1515/anona-2015-0096
    [12] X. L. Fan, Q. H. Zhang, Existence of solutions for p(x)-Laplacian Dirichlet problem, Nonliear Anal., 52 (2003), 1843–1852. doi: 10.1016/S0362-546X(02)00150-5
    [13] N. Nyamoradi, Existence of three solutions for Kirchhoff nonlocal operators of elliptic type, Math. Commun., 18 (2013), 489–502.
    [14] N. Nyamoradi, N. T. Chung, Existence of solutions to nonlocal Kirchhoff equations of elliptic type via genus theory, Electron. J. Differ. Equ., 2014 (2014), 1–12. doi: 10.1186/1687-1847-2014-1
    [15] W. J. Chen, S. B. Deng, The Nehari manifold for a fractional p-Laplacian system involving concave-convex nonlinearities, Nonlinear Anal., 27 (2016), 80–92. doi: 10.1016/j.nonrwa.2015.07.009
    [16] T. S. Hsu, Multiple positive solutions for a critical quasilinear elliptic system with concave-convex nonlinearities, Nonlinear Anal., 71 (2009), 2688–2698. doi: 10.1016/j.na.2009.01.110
    [17] B. Ricceri, On a three critical points theorem, Arch. Math., 75 (2000), 220–226. doi: 10.1007/s000130050496
    [18] S. A. Marano, D. Motreanu, On a three critical points theorem for non differentiable functions and applications to nonlinear boundary value problems, Nonlinear Anal., 48 (2002), 37–52. doi: 10.1016/S0362-546X(00)00171-1
    [19] X. L. Fan, S. G. Deng, Remarks on Ricceri's variational principle and applications to the p(x)-Laplacian equations, Nonlinear Anal., 67 (2007), 3064–3075. doi: 10.1016/j.na.2006.09.060
    [20] G. Bonanno, R. Livrea, Multiplicity theorems for the Dirichlet problem involving the p-Laplacian, Nonlinear Anal., 54 (2003), 1–7. doi: 10.1016/S0362-546X(03)00027-0
    [21] G. Bonanno, A critical points theorem and nonlinear differential problems, J. Global Optim., 28 (2004), 249–258. doi: 10.1023/B:JOGO.0000026447.51988.f6
    [22] E. Azroul, A. Benkirane, A. Boumazourh, M. Srati, Three solutions for a nonlocal fractional p-Kirchhoff type elliptic system, Appl. Anal., (2019), 1–18.
    [23] F. J. S. A. Corrêa, R. G. Nascimento, On a nonlocal elliptic system of p-Kirchhoff-type under Neumann boundary condition, Math. Comput. Model., 49 (2009), 598–604. doi: 10.1016/j.mcm.2008.03.013
    [24] F. Alessio, P. Patrizia, B. L. Zhang, p-fractional hardy-schrdinger-kirchhoff systems with critical nonlinearities, Adv. Nonlinear Anal., 8 (2018), 1–21.
    [25] J. H. Chen, X. J. Huang, C. X. Zhu, Existence of multiple solutions for nonhomogeneous schrödinger-kirchhoff system involving the fractional p-laplacian with sign-changing potential-sciencedirect, Comput. Math. Appl., 77 (2019), 2725–2739. doi: 10.1016/j.camwa.2019.01.004
    [26] E. Azroul, A. Boumazourh, Three solutions for a fractional (p(x,),q(x,))-Kirchhoff type elliptic system, J. Nonlinear Funct. Anal., 40 (2020), 1–19.
    [27] W. C. Bu, T. Q. An, G. J. Ye, S. Taarabti, Negative energy solutions for a new fractional p(x)-Kirchhoff problem without the (AR) condition, J. Funct. Space, 2021 (2021), 8888078.
    [28] Y. Wu, Z. H. Qiao, M. K. Hamdani, B. Y. Kou, L. B. Yang, A class of variable-order fractional p()-Kirchhoff-type systems, J. Funct. Space, 2021 (2021), 5558074.
    [29] B. Ricceri, A three critical points theorem revisited, Nonlinear Anal., 70 (2009), 3084–3089. doi: 10.1016/j.na.2008.04.010
    [30] X. L. Fan, D. Zhao, On the spaces Lp(x)(Ω) and Wk,p(x)(Ω), J. Math. Anal. Appl., 263 (2001), 424–446. doi: 10.1006/jmaa.2000.7617
    [31] O. Kováčik, J. Rákosník, On spaces Lp(x)(Ω) and W1,p(x)(Ω), Czech. Math. J., 41 (1991), 592–618. doi: 10.21136/CMJ.1991.102493
    [32] L. Diening, P. Harjulehto, P. Hästö, M. Ružička, Lebesgue and Sobolev spaces with variable exponents, Lecture Notes in Mathematics, Heidelberg: Springer, 2011.
    [33] S. Bahrouni, H. Ounaies, Strauss and lions type theorems for the fractional sobolev spaces with variable exponent and applications to nonlocal Kirchhoff-Choquard problem, Mediterr. J. Math., 18 (2021), 46. doi: 10.1007/s00009-020-01661-w
    [34] R. Biswas, S. Tiwari, On a class of Kirchhoff-Choquard equations involving variable-order fractional p()-Laplacian and without Ambrosetti-Rabinowitz type condition, arXiv. Available from: https://arXiv preprint arXiv:2005.09221.
    [35] D. E. Edmunds, J. Rákosník, Sobolev embeddings with variable exponent, Studia Math., 143 (2000), 267–293. doi: 10.4064/sm-143-3-267-293
    [36] R. Biswas, S. Tiwari, Multiplicity and uniform estimate for a class of variable order fractional p(x)-Laplacian problems with concave-convex nonlinearities, arXiv. Available from: https://arXiv.org/abs/1810.12960.
    [37] G. Bonanno, Multiple critical points theorems without the Palais-Smale condition, J. Math. Anal. Appl., 299 (2004), 600–614. doi: 10.1016/j.jmaa.2004.06.034
    [38] A. Bahrouni, V. Ţ. D. Rǎdulescu, On a new fractional Sobolev space and applications to nonlocal variational problems with variable exponent, Discret. Contin. Dyn. Syst., 11 (2018), 379–389.
    [39] E. Zeidler, Nonlinear functional analysis and applications, In: Nonlinear monotone operators, Springer-Verlag, New York, 1990.
    [40] G. W. Dai, R. Y. Ma, Solutions for a p(x)-Kirchhoff type equation with Neumann boundary data, Nonlinear Anal., 12 (2011), 2666–2680. doi: 10.1016/j.nonrwa.2011.03.013
    [41] B. Barrios, E. Colorado, A. de Pablo, U. Sánchez, On some critical problems for the fractional Laplacian operator, J Differ. Equ., 252 (2012), 6133–6162. doi: 10.1016/j.jde.2012.02.023
  • This article has been cited by:

    1. Mdi Begum Jeelani, Abeer S. Alnahdi, Mohammed A. Almalahi, Mohammed S. Abdo, Hanan A. Wahash, Nadiyah Hussain Alharthi, Tianqing An, Qualitative Analyses of Fractional Integrodifferential Equations with a Variable Order under the Mittag-Leffler Power Law, 2022, 2022, 2314-8888, 1, 10.1155/2022/6387351
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2386) PDF downloads(144) Cited by(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog