The present work is dedicated to a study that focuses on solving space-fractional advection-diffusion equations (SFADEs) using the Galerkin method. Through our analysis, we demonstrate the effectiveness of this approach in solving the considered equations. After introducing the Chebyshev cardinal functions (CCFs), the Caputo fractional derivative (CFD) was represented based on these bases as an operational matrix. Applying the Galerkin method reduces the desired equation to a system of algebraic equations. We have proved that the method converges analytically. By solving some numerical examples, we have demonstrated that the proposed method is effective and yields superior outcomes compared to existing methods for addressing this problem.
Citation: Haifa Bin Jebreen, Hongzhou Wang. On the effective method for the space-fractional advection-diffusion equation by the Galerkin method[J]. AIMS Mathematics, 2024, 9(9): 24143-24162. doi: 10.3934/math.20241173
[1] | Alexander Fekete . Urban and Rural Landslide Hazard and Exposure Mapping Using Landsat and Corona Satellite Imagery for Tehran and the Alborz Mountains, Iran. AIMS Geosciences, 2017, 3(1): 37-66. doi: 10.3934/geosci.2017.1.37 |
[2] | Roberto Gianardi, Marina Bisson, Lisa Beccaro, Riccardo De Ritis, Vincenzo Sepe, Laura Colini, Cristiano Tolomei, Luca Cocchi, Claudia Spinetti . High-resolution susceptibility mapping of seismically induced landslides on Ischia island: the 2017 earthquake case study. AIMS Geosciences, 2024, 10(3): 573-597. doi: 10.3934/geosci.2024030 |
[3] | Cem Kıncal, Zhenhong Li, Jane Drummond, Peng Liu, Trevor Hoey, Jan-Peter Muller . Landslide Susceptibility Mapping Using GIS-based Vector Grid File (VGF) Validating with InSAR Techniques: Three Gorges, Yangtze River (China). AIMS Geosciences, 2017, 3(1): 116-141. doi: 10.3934/geosci.2017.1.116 |
[4] | Chih-Wei Chuang, Hao-Yu Huang, Chun-Wei Tseng . Assessment of vegetation dynamic and its effects in a large-scale landslide in Central Taiwan with multitemporal Landsat images. AIMS Geosciences, 2025, 11(2): 318-342. doi: 10.3934/geosci.2025014 |
[5] | Raissa Mazova, Kirill Polyakov, Dmitry Smirnov, Ivan Smirnov, Andrey Kurkin, Jorge Van Den Bosch F . Pre-eruption landslide as a possible cause of early tsunami waves from Tongo volcano on January 15, 2022. AIMS Geosciences, 2025, 11(2): 274-297. doi: 10.3934/geosci.2025012 |
[6] | Tonje Eide Helle, Per Aagaard, Steinar Nordal, Michael Long, Sara Bazin . A geochemical, mineralogical and geotechnical characterization of the low plastic, highly sensitive glaciomarine clay at Dragvoll, Norway. AIMS Geosciences, 2019, 5(4): 704-722. doi: 10.3934/geosci.2019.4.704 |
[7] | Ariane Locat, Pascal Locat, Hubert Michaud, Kevin Hébert, Serge Leroueil, Denis Demers . Geotechnical characterization of the Saint-Jude clay, Quebec, Canada. AIMS Geosciences, 2019, 5(2): 273-302. doi: 10.3934/geosci.2019.2.273 |
[8] | Eduardo Antonio Rossello . Effects of the Guando Slump on the tectonic interpretation of the Boquerón Thrust in the Guando oilfield [Tolima, Colombia]. AIMS Geosciences, 2025, 11(2): 370-386. doi: 10.3934/geosci.2025016 |
[9] | Leila Maria Ramos, Thiago Bazzan, Mariana Ferreira Benessiuti Motta, George de Paula Bernardes, Heraldo Luiz Giacheti . Landslide susceptibility mapping based on rainfall scenarios: a case study from Sao Paulo in Brazil. AIMS Geosciences, 2022, 8(3): 438-451. doi: 10.3934/geosci.2022024 |
[10] | M.C. Raghucharan, Surendra Nadh Somala . Stochastic Extended Simulation (EXSIM) of Mw 7.0 Kumamoto-Shi earthquake on 15 April 2016 in the Southwest of Japan using the SCEC Broadband Platform (BBP). AIMS Geosciences, 2018, 4(2): 144-165. doi: 10.3934/geosci.2018.2.144 |
The present work is dedicated to a study that focuses on solving space-fractional advection-diffusion equations (SFADEs) using the Galerkin method. Through our analysis, we demonstrate the effectiveness of this approach in solving the considered equations. After introducing the Chebyshev cardinal functions (CCFs), the Caputo fractional derivative (CFD) was represented based on these bases as an operational matrix. Applying the Galerkin method reduces the desired equation to a system of algebraic equations. We have proved that the method converges analytically. By solving some numerical examples, we have demonstrated that the proposed method is effective and yields superior outcomes compared to existing methods for addressing this problem.
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)|x−y|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)|x−y|N+p(x,y)s(x,y) dxdy+∫Ω|ξ|¯p(x)¯p(x)dx)(Δs(⋅)p(⋅)ξ−|ξ|¯p(x)ξ)=λFξ(x,η,ξ)+μGξ(x,η,ξ),x∈Ω,η=ξ=0,x∈RN∖Ω, | (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))|x−y|N+p(x,y)s(x,y)dy, for allx∈RN, | (1.2) |
where v∈C∞0(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+E2L∫L0|∂ξ(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+E2L∫L0|∂ξ(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+0→R+,(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
Mi≥Mi(t)≥mifor all t>ι,(i=1,2) |
and put
˜Mi(t)=∫t0Mi(ς)dςfor all t∈R+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η=λ|η|q−2η+2αα+β|η|α−2η|ξ|β,x∈Ω,(−Δ)spξ=λ|ξ|q−2ξ+2βα+β|ξ|β−2ξ|η|α,x∈Ω,η=ξ=0,x∈RN∖Ω, | (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,x∈RN∖Ω, | (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)|x−y|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)|x−y|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)|x−y|N+p(x,y)s(x,y) dxdy)(−Δ)s(⋅)p(⋅)ξ(x)=g(η,ξ)+b(x),x∈Ω,η=ξ=0,x∈RN∖Ω, | (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)|x−y|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)|x−y|N+p(x,y)s(x,y) dxdy+∫Ω|ξ|¯p(x)¯p(x)dx)(Δs(⋅)p(⋅)ξ−|ξ|¯p(x)ξ)=λFξ(x,η,ξ),x∈Ω,η=ξ=0,x∈RN∖Ω, | (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:Ω×R2→R 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:Ω×R2→R 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)|x−y|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))|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|η|¯p(x)ηφdx)+M2(∫RN×RN|ξ(x)−ξ(y)|p(x,y)p(x,y)|x−y|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))|x−y|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:X0→R 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)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|υ|¯p(x)¯p(x)dx. | (1.13) |
The functions Φ,Ψ,J: X0→R 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))|x−y|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:Ω×R2→R 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,p∗s(x)).
Now, let us show our results in this article.
Theorem 1.1. For s(⋅),p(⋅)∈C+(D) with s+p+<N and F∈A, 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)<m−H(c1,a),M+A(c2)<m−H(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 F∈A, 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ξ :Ω×R2→R 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)dx≤1}, |
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)=0⇔limn→∞ρϑ(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)|x−y|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)|x−y|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′={ξ:RN→R:ξ|Ω∈L¯p(⋅)(Ω):and∫Q|ξ(x)−ξ(y)|p(x,y)χp(x,y)|x−y|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)|x−y|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)|x−y|N+p(x,y)s(x,y) dxdy=∫RN×RN|ξ(x)−ξ(y)|p(x,y)χp(x,y)|x−y|N+p(x,y)s(x,y) dxdy<1}. |
W∗0 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)<p∗s(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 W0↪Lϕ(⋅)(Ω) is compact.
Proposition 2.2. (see [36]) If we define
ρs(⋅)p(⋅)(ξ)=∫RN×RN|ξ(x)−ξ(y)|p(x,y)|x−y|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(⋅)(ξ)≤‖ξ‖p−W0,(3)‖ξ‖W0>1⇒‖ξ‖p−W0≤ρs(⋅)p(⋅)(ξ)≤‖ξ‖p+W0,(4)limn→∞‖ξn‖W0=0(resp.→+∞)⇔limn→∞ρs(⋅)p(⋅)(ξn)=0(resp.→+∞),(5)limn→∞‖ξn−ξ‖W0=0⇔limn→∞ρ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 Φ,Ψ:X→R 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 Φ:X→R 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; Ψ:X→R is a continuously Gâteaux differentiable functional whose Gâteaux derivative is compact, I⊆R an interval. Assume that
lim‖η‖→∞(Φ(η)+λΨ(η))=+∞, |
for λ∈I, and that there exists r∈R and η0,η1∈X 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:X→R 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,F∈A, then J,Ψ∈C1(X0,R) and their derivatives are defined as (1.14) and (1.15) for all (φ,ψ)∈X0. Moreover J′,Ψ′:X0→R 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 n∈N.
Fix (˜η,˜ξ)∈X0 with ‖(˜η,˜ξ)‖X0≤1. Since G∈A, combining Lemma 2.1 with Theorem 2.1, we obtain
|⟨J′(ηn,ξn)−J′(η,ξ),(˜η,˜ξ)⟩|≤|∫Ω[Gη(x,ηn,ξn)dx−Gη(x,η,ξ)]˜ηdx+∫Ω[Gξ(x,ηn,ξn)dx−Gξ(x,η,ξ)]˜ξdx|≤C3‖Gη(x,ηn,ξn)−Gη(x,η,ξ)‖Lϑ′(x)(Ω)‖˜η‖Lϑ(x)(Ω)+C3‖Gξ(x,ηn,ξn)−Gξ(x,η,ξ)‖Lϑ′(x)(Ω)‖˜ξ‖Lϑ(x)(Ω)≤C4‖Gη(x,ηn,ξn)−Gη(x,η,ξ)‖Lϑ′(x)(Ω)‖˜η‖W0+C4‖Gξ(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 ‖(˜η,˜ξ)‖X0≤1, we get
‖⟨J′(ηn,ξn)−J′(η,ξ),(˜η,˜ξ)⟩‖W∗0≤C4‖Gη(x,ηn,ξn)−Gξ(x,η,ξ)‖Lϑ′(x)(Ω)+C4‖Gξ(x,ηn,ξn)−Gξ(x,η,ξ)‖Lϑ′(x)(Ω). |
According to Definition 2, we deduce
Gη(x,ηn,ξn)−Gη(x,η,ξ)→0asn→∞a.e.x∈Ω,Gξ(x,ηn,ξn)−Gξ(x,η,ξ)→0asn→∞a.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)(Ω)→0asn→∞a.e.x∈Ω,‖Gξ(x,ηn,ξn)−Gξ(x,η,ξ))‖Lϑ′(x)(Ω)→0asn→∞a.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) & Φ′:X0→X∗0 is a strictly monotone and continuous operator,
(iii) & Φ′:X0→X∗0 is a homeomorphism.
Proof. (i) Since ˜M′i(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, Φ′(η,ξ)∈X∗0, where X∗0 denotes the dual spaces of X0.
In the first place, we show that Φ′:X0→X∗0 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 Φ′:X0→X∗0 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))|x−y|(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))|x−y|(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))|x−y|(N+p(x,y)s(x,y))/p′(x,y)→|η(x)−η(y)|p(x,y)−2(η(x)−η(y))|x−y|(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))|x−y|(N+p(x,y)s(x,y))/p′(x,y)→|ξ(x)−ξ(y)|p(x,y)−2(ξ(x)−ξ(y))|x−y|(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))|x−y|(N+p(x,y)s(x,y))/p′(x,y)−|η(x)−η(y)|p(x,y)−2(η(x)−η(y))|x−y|(N+p(x,y)s(x,y))/p′(x,y)|p′(x,y)dxdy=limn→∞∫RN×RN(|ηn(x)−ηn(y)|p(x,y)|x−y|N+p(x,y)s(x,y)−|η(x)−η(y)|p(x,y)|x−y|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))|x−y|(N+p(x,y)s(x,y))/p′(x,y)−|ξ(x)−ξ(y)|p(x,y)−2(ξ(x)−ξ(y))|x−y|(N+p(x,y)s(x,y))/p′(x,y)|p′(x,y)dxdy=limn→∞∫RN×RN(|ξn(x)−ξn(y)|p(x,y)|x−y|N+p(x,y)s(x,y)−|ξ(x)−ξ(y)|p(x,y)|x−y|N+p(x,y)s(x,y))dxdy. | (3.14) |
Based on the fact that ηn→η, ξn→ξ strongly in X0 yields
limn→∞∫RN×RN(|ηn(x)−ηn(y)|p(x,y)−2(ηn(x)−ηn(y))|x−y|(N+p(x,y)s(x,y))/p′(x,y)−|η(x)−η(y)|p(x,y)−2(η(x)−η(y))|x−y|(N+p(x,y)s(x,y))/p′(x,y))dxdy=0 | (3.15) |
and
limn→∞∫RN×RN(|ξn(x)−ξn(y)|p(x,y)−2(ξn(x)−ξn(y))|x−y|(N+p(x,y)s(x,y))/p′(x,y)−|ξ(x)−ξ(y)|p(x,y)−2(ξ(x)−ξ(y))|x−y|(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)−Φ′(η,ξ)‖X∗0=sup(φ,ψ)∈X0,‖(φ,ψ)‖X0≤1|⟨Φ′(η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(⋅)(ξ))⟨ξ,ξ⟩‖(η,ξ)‖X0≥m1ρs(⋅)p(⋅)(η)+m2ρs(⋅)p(⋅)(ξ)‖(η,ξ)‖X0≥min(m1,m2)min(‖η‖p+W0,‖η‖p−W0)+min(‖ξ‖p+W0,‖ξ‖p−W0)‖η‖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 X∗0 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)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ηn|¯p(x)dx)≥∫RN×RN|η(x)−η(y)|p(x,y)|x−y|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)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ξn|¯p(x)dx)≥∫RN×RN|ξ(x)−ξ(y)|p(x,y)|x−y|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)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ηn|¯p(x)dx−∫RN×RN|ηn(x)−ηn(y)|p(x,y)−2(ηn(x)−ηn(y))(η(x)−η(y))|x−y|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)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ξn|¯p(x)dx−∫RN×RN|ξn(x)−ξn(y)|p(x,y)−2(ξn(x)−ξn(y))(ξ(x)−ξ(y))|x−y|N+p(x,y)s(x,y)dxdy−∫Ω|ξn|¯p(x)−2ξnξdx)≥m1(∫RN×RN|ηn(x)−ηn(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ηn|¯p(x)dx−∫RN×RN|ηn(x)−ηn(y)|p(x,y)−1(η(x)−η(y))|x−y|N+p(x,y)s(x,y)dxdy−∫Ω|ηn|¯p(x)−1ηdx)+m2(∫RN×RN|ξn(x)−ξn(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ξn|¯p(x)dx−∫RN×RN|ξn(x)−ξn(y)|p(x,y)−1(ξ(x)−ξ(y))|x−y|N+p(x,y)s(x,y)dxdy−∫Ω|ξn|¯p(x)−1ξdx)≥m1(C8∫RN×RN|ηn(x)−ηn(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy+C8∫Ω|ηn|¯p(x)dx−C9∫RN×RN|η(x)−η(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy−C9∫Ω|η|¯p(x)dx)+m2(C8∫RN×RN|ξn(x)−ξn(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy+C8∫Ω|ξn|¯p(x)dx−C9∫RN×RN|ξ(x)−ξ(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy−C9∫Ω|ξ|¯p(x)dx). | (3.26) |
The passage to the liminf implies in the above inequality, we get
0≥m1[lim infn→∞(∫RN×RN|ηn(x)−ηn(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ηn|¯p(x)dx)−∫RN×RN|η(x)−η(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy−∫Ω|η|¯p(x)dx]+m2[lim infn→∞(∫RN×RN|ξn(x)−ξn(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ξn|¯p(x)dx)−∫RN×RN|ξ(x)−ξ(y)|p(x,y)|x−y|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)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ηn|¯p(x)dx)=∫RN×RN|η(x)−η(y)|p(x,y)|x−y|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)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ξn|¯p(x)dx)=∫RN×RN|ξ(x)−ξ(y)|p(x,y)|x−y|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)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ηn|¯p(x)dx)=∫RN×RN|η(x)−η(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|η|¯p(x)dx | (3.30) |
and
limn→∞(∫RN×RN|ξn(x)−ξn(y)|p(x,y)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|ξn|¯p(x)dx)=∫RN×RN|ξ(x)−ξ(y)|p(x,y)|x−y|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−ξ‖W0→0. | (3.35) |
Next, we prove that (Φ′)−1 is continuous operator. Let {(fn,gn)}⊂X∗0 such that fn→f0 in X∗0 and gn→g0 in X∗0. 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)=Φ′(η,ξ)inX∗0. | (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 Φ′:X∗0→X0.
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),aR2−R1{R2−[ΣNi=1(xi−xi0)]1/2},x∈B(x0,R2)∖B(x0,R1),a,x∈B(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)‖p−W0+‖ξ(x)‖p−W0})≤Φ(η,ξ)≤M+(max{‖η(x)‖p+W0+‖ξ(x)‖p+W0,‖η(x)‖p−W0+‖ξ(x)‖p−W0}). | (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=m−c1γCϕ,r2=M+c2γCϕ. | (4.5) |
By virtue of (4.1), set η0=ξ0=ω for (η0,ξ0)∈X0 and
‖η0‖W0=‖ω‖W0,‖ξ0‖W0=‖ω‖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+m−A(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)‖p−W0+‖ξ(x)‖p−W0})≤Cϕr1m−=c1γ≤c1, for allx∈Ω. | (4.10) |
Fix γ0 such that 0<max{‖η0‖p+W0+‖ξ0‖p+W0,‖η0‖p−W0+‖ξ0‖p−W0}<γ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+(‖η0‖p(⋅)W0+‖ξ0‖p(⋅)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)‖p−W0+‖ξ(x)‖p−W0})≤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ϕrm−F(x,η,ξ)dxr, | (4.14) |
this implies
φi(ri)≤A(ci)ri=Cϕm−ciγ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)‖p−W0+‖ξi(x)‖p−W0})≤Cϕr2m−≤c2, 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 Φ′:X∗0→X0. Moreover,
lim‖(η,ξ)‖x0→∞Φ(η,ξ)+λΨ(η,ξ)=+∞ | (4.19) |
for any λ∈(0,+∞). Indeed,
Φ(η,ξ)=˜M1(∫RN×RN|η(x)−η(y)|p(x,y)p(x,y)|x−y|N+p(x,y)s(x,y)dxdy+∫Ω|η|¯p(x)¯p(x)dx)+˜M2(∫RN×RN|ξ(x)−ξ(y)|p(x,y)p(x,y)|x−y|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(‖η‖p−W0,‖η‖p+W0)+m2p+min(‖ξ‖p−W0,‖ξ‖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,η,ξ)dx≥−C10∫Ω|η|α(x)+1|ξ|β(x)+1dx≥−C10|Ω|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 bp−1|Ω|>1. From (4.28), we have
∫ΩF(x,b1,b1)>0=F(x,0,0)>∫Ωsup0≤s,t≤a0F(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)≥|Ω|(m1bp−1p++m2bp−1p+)≥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
‖η‖w0≤1,‖ξ‖w0≤1. | (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≤∫Ωsup0≤s,t≤aF(x,s,t)dx≤0≤−(Φ(η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 G∈A, 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. Liu, F. Geng, An explanation on four new definitions of fractional operators, Acta Math. Sci., 44 (2024), 1271–1279. https://doi.org/10.1007/s10473-024-0405-7 doi: 10.1007/s10473-024-0405-7
![]() |
[2] |
M. Arif, F. Ali, I. Khan, K. S. Nisar, A time fractional model with non-singular kernel the generalized couette flow of couple stress nanofluid, IEEE Access, 8 (2020), 77378–77395. https://doi.org/10.1109/ACCESS.2020.2982028 doi: 10.1109/ACCESS.2020.2982028
![]() |
[3] |
A. Chang, H. Sun, C. Zheng, B. Lu, C. Lu, R. Ma, et al., A time fractional convection-diffusion equation to model gas transport through heterogeneous soil and gas reservoirs, Physica A, 502 (2018), 356–369. https://doi.org/10.1016/j.physa.2018.02.080 doi: 10.1016/j.physa.2018.02.080
![]() |
[4] |
J. A. Tenreiro Machado, M. F. Silva, R. S. Barbosa, I. S. Jesus, C. M. Reis, M. G. Marcos, et al., Some applications of fractional calculus in engineering, Math. Probl. Eng., 2010 (2010), 1–34. https://doi.org/10.1155/2010/639801 doi: 10.1155/2010/639801
![]() |
[5] | F. Mainardi, Fractional calculus and waves in linear viscoelasticity, Imperial College Press, 2010. https://doi.org/10.1142/p614 |
[6] |
M. Asadzadeh, B. N. Saray, On a multiwavelet spectral element method for integral equation of a generalized Cauchy problem, BIT Numer. Math., 62 (2022), 1383–1416. https://doi.org/10.1007/s10543-022-00915-1 doi: 10.1007/s10543-022-00915-1
![]() |
[7] |
A. K. Gupta, S. Saha Ray, Wavelet methods for solving fractional order differential equations, Math. Probl. Eng., 2014 (2014), 140453. https://doi.org/10.1155/2014/140453 doi: 10.1155/2014/140453
![]() |
[8] |
L. Shi, B. N. Saray, F. Soleymani, Sparse wavelet Galerkin method: Application for fractional Pantograph problem, J. Comput. Appl. Math., 451 (2024), 116081. https://doi.org/10.1016/j.cam.2024.116081 doi: 10.1016/j.cam.2024.116081
![]() |
[9] |
P. Thanh Toan, T. N. Vo, M. Razzaghi, Taylor wavelet method for fractional delay differential equations, Eng. Comput., 37 (2021), 231–240. https://doi.org/10.1007/s00366-019-00818-w doi: 10.1007/s00366-019-00818-w
![]() |
[10] |
V. Daftardar-Gejji, A. Jafari, Adomian decomposition: A tool for solving a system of fractional differential equations, J. Math. Anal. Appl., 301 (2005), 508–518. https://doi.org/10.1016/j.jmaa.2004.07.039 doi: 10.1016/j.jmaa.2004.07.039
![]() |
[11] |
B. Benkerrouche, D. Baleanu, M. S. Souid, A. Hakem, Boundary value problem for nonlinear fractional differential equations of variable order via Kuratowski MNC technique, Adv. Differ. Equ., 2021 (2021), 365. https://doi.org/10.1186/s13662-021-03520-8 doi: 10.1186/s13662-021-03520-8
![]() |
[12] |
M. Lakestani, M. Dehghan, The use of Chebyshev cardinal functions for the solution of a partial differential equation with an unknown time-dependent coefficient subject to an extra measurement, J. Comput. Appl. Math., 235 (2010), 669–678. https://doi.org/10.1016/j.cam.2010.06.020 doi: 10.1016/j.cam.2010.06.020
![]() |
[13] |
G. J. Fix, J. P. Roop, Least squares finite element solution of a fractional order two-point boundary value problem, Comput. Math. Appl., 48 (2004), 1017–1033. https://doi.org/10.1016/j.camwa.2004.10.003 doi: 10.1016/j.camwa.2004.10.003
![]() |
[14] |
S. Maji, S. Natesan, Adaptive-grid technique for the numerical solution of a class of fractional boundary-value-problems, Comput. Methods Differ. Equ., 12 (2010), 338–349. https://doi.org/10.22034/CMDE.2023.55266.2296 doi: 10.22034/CMDE.2023.55266.2296
![]() |
[15] |
R. Garrappa, On some explicit Adams multistep methods for fractional differential equations, J. Comput. Appl. Math., 229 (2009), 392–399. https://doi.org/10.1016/j.cam.2008.04.004 doi: 10.1016/j.cam.2008.04.004
![]() |
[16] |
Y. L. Zhao, P. Zhu, X. M. Gu, X. Zhao, H. Y. Jian, An implicit integration factor method for a kind of spatial fractional diffusion equations, J. Phys. Conf. Ser., 1324 (2019), 012030. https://doi.org/10.1088/1742-6596/1324/1/012030 doi: 10.1088/1742-6596/1324/1/012030
![]() |
[17] |
Z. Lin, D. Wang, D. Qi, L. Deng, A Petrov-Galerkin finite element-meshfree formulation for multi-dimensional fractional diffusion equations, Comput. Mech., 66 (2020), 323–350. https://doi.org/10.1007/s00466-020-01853-x doi: 10.1007/s00466-020-01853-x
![]() |
[18] |
H. Y. Jian, T. Z. Huang, X. L. Zhao, Y. L. Zhao, A fast second-order accurate difference schemes for time distributed-order and Riesz space fractional diffusion equations, J. Appl. Anal. Comput., 9 (2019), 1359–1392. https://doi.org/10.11948/2156-907X.20180247 doi: 10.11948/2156-907X.20180247
![]() |
[19] |
P. Biler, W. A. Woyczynski, Global and exploding solutions for nonlocal quadratic evolution problems, SIAM J. Appl. Math., 59 (1998), 845–869. https://doi.org/10.1137/S003613999631344 doi: 10.1137/S003613999631344
![]() |
[20] |
W. Y. Tian, W. Deng, Y. Wu, Polynomial spectral collocation method for space fractional advection-diffusion equation, Numer. Meth. Partial Differ. Equ., 30 (2014), 514–535. https://doi.org/10.1002/num.21822 doi: 10.1002/num.21822
![]() |
[21] |
R. Metzler, J. Klafter, The random walks guide to anomalous diffusion: A fractional dynamics approach, Phys. Rep., 339 (2000), 1–77. https://doi.org/10.1016/S0370-1573(00)00070-3 doi: 10.1016/S0370-1573(00)00070-3
![]() |
[22] |
A. Pablo, F. Quirós, A. Rodríguez, J. L. Vázquez, A fractional porous medium equation, Adv. Math., 226 (2011), 1378–1409. https://doi.org/10.1016/j.aim.2010.07.017 doi: 10.1016/j.aim.2010.07.017
![]() |
[23] |
Y. Zheng, C. Li, Z. Zhao, A note on the finite element method for the space-fractional advection diffusion equation, Comput. Math. Appl., 59 (2010), 1718–1726. https://doi.org/10.1016/j.camwa.2009.08.071 doi: 10.1016/j.camwa.2009.08.071
![]() |
[24] |
H. Hejazi, T. Moroney, F. Liu, Stability and convergence of a finite volume method for the space fractional advection-dispersion equation, J. Comput. Appl. Math., 255 (2014), 684–697. https://doi.org/10.1016/j.cam.2013.06.039 doi: 10.1016/j.cam.2013.06.039
![]() |
[25] |
A. Jannelli, M. Ruggieri, M. P. Speciale, Numerical solutions of space-fractional advection-diffusion equations with nonlinear source term, Appl. Numer. Math., 155 (2020), 93–102. https://doi.org/10.1016/j.apnum.2020.01.016 doi: 10.1016/j.apnum.2020.01.016
![]() |
[26] |
H. Y. Jian, T. Z. Huang, X. M. Gu, Y. L. Zhao, Compact implicit integration factor method for two-dimensional space-fractional advection-diffusion-reaction equations, J. Phys. Conf. Ser., 1592 (2020), 012048. https://doi.org/10.1088/1742-6596/1592/1/012048 doi: 10.1088/1742-6596/1592/1/012048
![]() |
[27] |
H. Y. Jian, T. Z. Huang, A. Ostermann, X. M. Gu, Y. L. Zhao, Fast IIF-WENO method on non-uniform meshes for nonlinear space-fractional convection-diffusion-reaction equations, J. Sci. Comput., 89 (2021), 13. https://doi.org/10.1007/s10915-021-01622-9 doi: 10.1007/s10915-021-01622-9
![]() |
[28] | J. P. Boyd, Chebyshev and fourier spectral methods, Dover Publications, 2001. |
[29] |
A. Afarideh, F. D. Saei, M. Lakestani, B. N. Saray, Pseudospectral method for solving fractional Sturm-Liouville problem using Chebyshev cardinal functions, Phys. Scr., 96 (2021), 125267. https://doi.org/10.1088/1402-4896/ac3c59 doi: 10.1088/1402-4896/ac3c59
![]() |
[30] |
A. Afarideh, F. D. Saei, B. N. Saray, Eigenvalue problem with fractional differential operator: Chebyshev cardinal spectral method, J. Math. Model., 11 (2021), 343–355. https://doi.org/10.22124/JMM.2023.24239.2169 doi: 10.22124/JMM.2023.24239.2169
![]() |
[31] |
F. Tchier, I. Dassios, F. Tawfiq, L. Ragoub, On the approximate solution of partial integro-differential equations using the pseudospectral method based on Chebyshev cardinal functions, Mathematics, 9 (2021), 286. https://doi.org/10.3390/math9030286 doi: 10.3390/math9030286
![]() |
[32] |
M. Shahriari, B. N. Saray, B. Mohammadalipour, S. Saeidian, Pseudospectral method for solving the fractional one-dimensional Dirac operator using Chebyshev cardinal functions, Phys. Scr., 98 (2023), 055205. https://doi.org/10.1088/1402-4896/acc7d3 doi: 10.1088/1402-4896/acc7d3
![]() |
[33] | K. Sayevand, H. Arab, An efficient extension of the Chebyshev cardinal functions for differential equations with coordinate derivatives of non-integer order, Comput. Methods Differ. Equ., 6 (2018), 339–352. https://doi.org/20.1001.1.23453982.2018.6.3.6.5 |
[34] | A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Amsterdam: Elsevier, 2006. |
[35] |
Y. Saad, M. H. Schultz, GMRES: A generalized minimal residual method for solving nonsymmetric linear systems, SIAM J. Sci. Stat. Comput., 7 (1986), 856–869. https://doi.org/10.1137/090705 doi: 10.1137/090705
![]() |
[36] | G. Dahlquist, A. Björck, Numerical methods, Englewood Cliffs: Prentice-Hall, 1974. |
[37] |
B. N. Saray, M. Lakestani, M. Dehghan, On the sparse multiscale representation of 2‐D Burgers equations by an efficient algorithm based on multiwavelets, Numer. Math. Partial Differ. Equ., 39 (2023), 1938–1961. https://doi.org/10.1002/num.22795 doi: 10.1002/num.22795
![]() |
[38] |
G. Pang, W. Chen, Z. Fu, Space-fractional advection-dispersion equations by the Kansa method, J. Comput. Phys., 293 (2015), 280–296. https://doi.org/10.1016/j.jcp.2014.07.020 doi: 10.1016/j.jcp.2014.07.020
![]() |
[39] |
M. M. Meerschaert, H. P. Scheffler, C. Tadjeran, Finite difference methods for two-dimensional fractional dispersion equation, J. Comput. Phys., 211 (2006), 249–261. https://doi.org/10.1016/j.jcp.2005.05.017 doi: 10.1016/j.jcp.2005.05.017
![]() |
1. | Vincenzo Barrile, Luigi Bibbò, Giuliana Bilotta, Giuseppe M. Meduri, Emanuela Genovese, Geomatics Innovation and Simulation for Landslide Risk Management: The Use of Cellular Automata and Random Forest Automation, 2024, 14, 2076-3417, 11853, 10.3390/app142411853 |