Loading [MathJax]/jax/output/SVG/jax.js
Research article

A nonlinear correction finite volume scheme preserving maximum principle for diffusion equations with anisotropic and discontinuous coefficient

  • A nonlinear correction finite volume scheme preserving the discrete maximum principle (DMP) was presented to solve diffusion equations with anisotropic and discontinuous coefficients. It is well-known that existing cell-centered finite volume schemes for the diffusion problem with the general discontinuous coefficient often impose severe restrictions on the mesh-cell geometry to maintain the DMP. We proposed a nonlinear method for modifying the flux to obtain a new scheme which eliminated the requirement for nonnegative interpolation coefficients at the midpoint of cell-edge unknowns while still preserving the DMP. That is, our new scheme satisfied the DMP unconditionally and can be applied to the diffusion problem with the discontinuous coefficient on arbitrary distorted meshes. We then provided a priori estimation under a coercivity assumption and proved that the scheme satisfied the DMP. Numerical results were presented to demonstrate that our scheme can handle diffusion equations with anisotropic and discontinuous coefficients, satisfied the DMP, and, in some cases, outperformed existing schemes which preserved the DMP in terms of accuracy.

    Citation: Yao Yu, Guanyu Xue. A nonlinear correction finite volume scheme preserving maximum principle for diffusion equations with anisotropic and discontinuous coefficient[J]. Electronic Research Archive, 2025, 33(3): 1589-1609. doi: 10.3934/era.2025075

    Related Papers:

    [1] Aziz Belmiloudi . Time-varying delays in electrophysiological wave propagation along cardiac tissue and minimax control problems associated with uncertain bidomain type models. AIMS Mathematics, 2019, 4(3): 928-983. doi: 10.3934/math.2019.3.928
    [2] Zuliang Lu, Xiankui Wu, Fei Huang, Fei Cai, Chunjuan Hou, Yin Yang . Convergence and quasi-optimality based on an adaptive finite element method for the bilinear optimal control problem. AIMS Mathematics, 2021, 6(9): 9510-9535. doi: 10.3934/math.2021553
    [3] Zahra Pirouzeh, Mohammad Hadi Noori Skandari, Kamele Nassiri Pirbazari, Stanford Shateyi . A pseudo-spectral approach for optimal control problems of variable-order fractional integro-differential equations. AIMS Mathematics, 2024, 9(9): 23692-23710. doi: 10.3934/math.20241151
    [4] Xin Yi, Rong Liu . An age-dependent hybrid system for optimal contraception control of vermin. AIMS Mathematics, 2025, 10(2): 2619-2633. doi: 10.3934/math.2025122
    [5] Yuanyuan Cheng, Yuan Li . A novel event-triggered constrained control for nonlinear discrete-time systems. AIMS Mathematics, 2023, 8(9): 20530-20545. doi: 10.3934/math.20231046
    [6] Woocheol Choi, Young-Pil Choi . A sharp error analysis for the DG method of optimal control problems. AIMS Mathematics, 2022, 7(5): 9117-9155. doi: 10.3934/math.2022506
    [7] Xiang Wu, Yuzhou Hou, Kanjian Zhang . Optimal feedback control for a class of fed-batch fermentation processes using switched dynamical system approach. AIMS Mathematics, 2022, 7(5): 9206-9231. doi: 10.3934/math.2022510
    [8] Tainian Zhang, Zhixue Luo . Optimal harvesting for a periodic competing system with size structure in a polluted environment. AIMS Mathematics, 2022, 7(8): 14696-14717. doi: 10.3934/math.2022808
    [9] Qian Li, Zhenghong Jin, Linyan Qiao, Aichun Du, Gang Liu . Distributed optimization of nonlinear singularly perturbed multi-agent systems via a small-gain approach and sliding mode control. AIMS Mathematics, 2024, 9(8): 20865-20886. doi: 10.3934/math.20241015
    [10] Asaf Khan, Gul Zaman, Roman Ullah, Nawazish Naveed . Optimal control strategies for a heroin epidemic model with age-dependent susceptibility and recovery-age. AIMS Mathematics, 2021, 6(2): 1377-1394. doi: 10.3934/math.2021086
  • A nonlinear correction finite volume scheme preserving the discrete maximum principle (DMP) was presented to solve diffusion equations with anisotropic and discontinuous coefficients. It is well-known that existing cell-centered finite volume schemes for the diffusion problem with the general discontinuous coefficient often impose severe restrictions on the mesh-cell geometry to maintain the DMP. We proposed a nonlinear method for modifying the flux to obtain a new scheme which eliminated the requirement for nonnegative interpolation coefficients at the midpoint of cell-edge unknowns while still preserving the DMP. That is, our new scheme satisfied the DMP unconditionally and can be applied to the diffusion problem with the discontinuous coefficient on arbitrary distorted meshes. We then provided a priori estimation under a coercivity assumption and proved that the scheme satisfied the DMP. Numerical results were presented to demonstrate that our scheme can handle diffusion equations with anisotropic and discontinuous coefficients, satisfied the DMP, and, in some cases, outperformed existing schemes which preserved the DMP in terms of accuracy.



    A century ago, the first metric fixed-point theorem was published by Banach [1]. In fact, before Banach, some famous mathematicians such as Picard and Liouville had used the fixed point approach to solve certain differential equations, more precisely, initial value problems. Inspired by their results, Banach considered it as a separated and independent result in the framework of the nonlinear functional analysis and point-set topology. The statement and the proof of an outstanding work of Banach, also known as contraction mapping principle, can be considered as an art piece: Each contraction, in the setting of a complete metric space, possesses a unique fixed point. Metric fixed point theory has been appreciated and investigated by several researchers. These researchers have different reasons and motivations to study this theory. The most important reason why the researchers find worthful to work and investigate metric fixed point theory is the natural and strong connection of the theoretical result in nonlinear functional analysis with applied sciences. If we look at it with the chronological aspect, we note that the fixed point theory was born as a tool to solve certain differential equations. Banach liberated the theory from being a tool in applied mathematics to an independent work of nonlinear functional analysis. On Picard and Liouville's side, it is a tool to solve the initial value problem. On the other side, from Banach's point of view, the fixed point theory is an independent research topic that has enormous application potential on almost all qualitative sciences, including applied mathematics. Secondly, Banach fixed point theorem not only guarantee the solution (the existence of a fixed point) but also indicate how we reach the mentioned solution (how to find the fixed point). Finally, we need to underline that almost all real world problems can be transferred to a fixed point problem, easily.

    With this motivation, several generalizations and extensions of Banach's fixed point theory have been released by introducing new contractions or by changing the structure of the studied abstract space. Among, we shall mention only a few of them that set up the skeleton of the contraction dealt with it. Historically, the first contraction we shall focus on it is the Meir-Keeler contraction [2]. Roughly speaking, the Meir-Keeler contraction can be considered as a uniform contraction. The second contraction that we dealt with is Jaggi contraction [3]. The interesting part of Jaggi's contraction is the following: Jaggi's contraction is one of the first of its kind that involves some rational expression. The last one is called as an interpolative contraction [4]. In the interpolative contraction the terms are used exponentially instead of standard usage of them.

    In this paper, we shall introduce a new contraction, hybrid Jaggi-Meir-Keeler type contraction, as a unification and generalization of the Meir-Keeler's contraction, the Jaggi's contraction and interpolative contraction in the setting of a complete metric space. We propose certain assumptions to guarantee the existence of a fixed point for such mappings. In addition, we express some example to indicate the validity of the derived results.

    Before going into details, we would like to reach a consensus by explaining the concepts and notations: Throughout the paper, we presume the sets, we deal with, are non-empty. The letter N presents the set of positive integers. Further, we assume that the pair (X,d) is a complete metric space. This notation is required in each of the following theorems, definitions, lemma and so on. We shall use the pair (X,d) everywhere without repeating that it is a complete metric space.

    In what follows we recall the notion of the uniform contraction which is also known as Meir-Keeler contraction:

    Definition 1.1. [2] A mapping f:(X,d)(X,d) is said to be a Meir-Keeler contraction on X, if for every E>0, there exists δ>0 such that

    Ed(x,y)<E+δimpliesd(fx,fy)<E, (1.1)

    for every x,yX.

    Theorem 1.1. [2] Any Meir-Keeler contraction f:(X,d)(X,d) possesses a unique fixed point.

    Very recently, Bisht and Rakočević [5] suggested the following extension of the uniform contraction:

    Theorem 1.2. [5] Suppose a mapping f:(X,d)(X,d) fulfills the following statements:

    (1) for a given E>0 there exists a δ(E)>0 such that

    E<M(x,y)<E+δ(E)impliesd(fx,fy)E;

    (2) d(fx,fy)<M(x,y), whenever M(x,y)>0;

    for any x,yX, where

    M(x,y)=max{d(x,y),αd(x,fx)+(1α)d(y,fy),(1α)d(x,fx)+αd(y,fy),β[d(x,fy)+d(y,fx)]2},

    with 0<α<1,0β<1.Then, f has a unique fixed point uX and fnxu for each xX.

    On the other hand, in 2018, the idea of interpolative contraction was consider to revisit the well-known Kannan's fixed point theorem [6]:

    Definition 1.2. [4] A mapping f:(X,d)(X,d) is said to be an interpolative Kannan type contraction on X if there exist κ[0,1) and γ(0,1) such that

    d(fx,fy)κ[d(x,fx)]γ[d(y,fy)]1γ, (1.2)

    for every x,yXFix(f), where Fix(f)={xX|fx=x}.

    Theorem 1.3. [4] Any interpolative Kannan-contraction mapping f:(X,d)(X,d) possesses a fixed point.

    For more interpolative contractions results, we refer to [7,8,9,10,11] and related references therein.

    Definition 1.3. A mapping f:(X,d)(X,d) is called a Jaggi type hybrid contraction if there is ψΨ so that

    d(fx,fy)ψ(Jsf(x,y)), (1.3)

    for all distinct x,yX where p0 and σi0,i=1,2,3,4, such that σ1+σ2=1 and

    Jsf(x,y)={[σ1(d(x,fx)d(y,fy)d(x,y))s+σ2(d(x,y))s]1/p,ifp>0,x,yX,xy(d(x,fx))σ1(d(y,fy))σ2,ifp=0,x,yXFf(X), (1.4)

    where Ff(X)={zX:fz=z}.

    Theorem 1.4. A continuous mapping f:(X,d)(X,d) possesses a fixed point x if it forms a Jaggi-type hybrid contraction.In particular, for any x0X, the sequence {fnx0} converges to x.

    Definition 1.4. [12] Let α:X×X[0,+) be a mapping, where X. A self-mapping f:(X,d)(X,d) is called triangular α-orbital admissible and denote as fTαX if

    α(x,fx)1impliesα(fx,f2x)1,

    and

    α(x,y)1,andα(y,fy)1,impliesα(x,fy)1

    for all x,yX.

    This concept, was used by many authors, in order to prove variant fixed point results (see, for instance [13,14,15,16,17,18,19] and the corresponding references therein).

    Lemma 1.1. [12] Assume that fTαX. If there exists x0X such that α(x0,fx0)1, then α(xm,xk)1, for all m,nN, where the sequence {xk} is defined by xk+1=xk.

    The following condition is frequently considered to avoid the continuity of the mappings involved.

    (R) if the sequence {xn} in X is such that for each nN,

    α(xn,xn+1)1andlimn+xn=xx,

    then there exists a subsequence {xn(j)} of {xn} such that

    α(xn(j),x)1, for each jN.

    We start this section by introducing the new contraction, namely, hybrid Jaggi-Meir-Keeler type contraction.

    Consider the mapping f:(X,d)(X,d) and the set of fixed point, Ff(X)={zX:fz=z}. We define the crucial expression Rsf as follows:

    Rsf(x,y)={[β1(d(x,fx)d(y,fy)d(x,y))s+β2(d(x,y))s+β3(d(x,fy)+d(y,fx)4)s]1/s,fors>0,x,yX,xy(d(x,fx))β1(d(y,fy))β2(d(x,fy)+d(y,fx)4)β3,fors=0,x,yX, (2.1)

    where p1 and βi0, i=1,2,3 are such that β1+β2+β3=1.

    Definition 2.1. Assume that fTαX. We say that f:(X,d)(X,d) is an α-hybrid Jaggi-Meir-Keeler type contraction on X, if for all distinct x,yX we have:

    (a1) for given E>0, there exists δ>0 such that

    E<max{d(x,y),Rsf(x,y)}<E+δα(x,y)d(fx,fy)E; (2.2)

    (a2)

    α(x,y)d(fx,fy)<max{d(x,y),Rsf(x,y)}. (2.3)

    Theorem 2.1. Any continuous α-hybrid Jaggi-Meir-Keeler type contraction f:(X,d)(X,d) provide a fixed point if there exists x0X, such that α(x0,fx0)1 and α(x0,f2x0)1.

    Proof. Let x0X be an arbitrary, but fixed point. We form the sequence {xm}, as follows:

    xm=fxm1=fmx0,

    for all mN and assume that d(xm,xm+1)>0, for all nN{0}. Indeed, if for some l0N{0} we have d(xl0,xl0+1)=0, it follows that xl0=xl0+1=fxl0. Therefore, xl0 is a fixed point of the mapping f and the proof is closed.

    Since, by assumption, the mapping f is triangular α-orbital admissible, it follows that

    α(x0,fx0)1α(x1,x2)=α(fx0,f2x0)1...
    α(xn,xn+1)1, (2.4)

    for every nN.

    We shall study two cases; these are s>0 and s=0.

    Case (A). For the case s>0, letting x=xn1 and y=xn=fxn1 in (a2), we get

    d(xn,xn+1)α(xn1,xn)d(fxn1,fxn)<max{d(xn1,xn),Rsf(xn1,xn)}, (2.5)

    where

    Rsf(xn1,xn)=[β1(d(xn1,fxn1)d(xn,fxn)d(xn1,xn))s+β2(d(xn1,xn))s++β3(d(xn1,fxn)+d(xn,fxn1)4)s]1/s=[β1(d(xn1,xn)d(xn,xn+1)d(xn1,xn))s+β2(d(xn1,xn))s++β3(d(xn1,xn+1)+d(xn,xn)4)s]1/s[β1(d(xn,xn+1))s+β2(d(xn1,xn))s++β3(d(xn1,xn)+d(xn,xn+1)4)s]1/s.

    If we can find n0N such that d(xn0,xn0+1)d(xn01,xn0), we have

    Rsf(xn01,xn0)[β1(d(xn0,xn0+1))s+β2(d(xn0,xn0+1))s++β3(d(xn0,xn0+1))s]1/s=d(xn0,xn0+1)(β1+β2+β3)1/s=d(xn0,xn0+1).

    Then, max{d(xn0,xn0+1),Rsf(xn01,xn0)}=d(xn0,xn0+1), and using (2.4), respectively (2.5) we get

    d(xn0,xn0+1)α(xn01,xn0)d(fxn+01,fxn0)<max{d(xn0,xn0+1),Rsf(xn01,xn0)}d(xn0,xn0+1),

    which is a contradiction. Therefore, d(xn,xn+1)<d(xn1,xn) for all nN and (2.5) becomes

    d(xn,xn+1)<d(xn1,xn),

    for all nN. Consequently, there exists b0 such that limn+d(xn1,xn)=b. If b>0, we have

    d(xm,xm+1)b>0,

    for any mN. On the one hand, since (2.2) holds for every given E>0, it is possible to choose E=b and let δ>0 be such that (2.2) is satisfied. On the other hand, since, also, limn+max{d(xn1,xn),Rsf(xn1,xn)}=b, there exists m0N such that

    b<max{d(xm01,xm0),Rsf(xm01,xm0)}<b+δ.

    Thus, by (2.2), together with (2.4) we obtain

    d(xm0,xm0+1)α(xm0,xm0+1)d(fxm01,fxm0)<b,

    which is a contradiction. Therefore,

    limn+d(xn,xn+1)=b=0. (2.6)

    We claim now that {xn} is a Cauchy sequence. Let E>0 be fixed and we can choose that δ=min{δ(E),E,1}. Thus, from (2.6) it follows that there exists j0N such that

    d(xn,xn+1)<δ2, (2.7)

    for all nj0. Now, we consider the set

    A={xl|lj0,d(xl,xj0)<E+δ2}. (2.8)

    We claim that fyA whenever y=xlA. Indeed, in case of l=j0, we have fxl=fxj0=xj0+1, and taking (2.7) into account, we get

    d(xj0,xj0+1)<δ2<E+δ2. (2.9)

    Thus, we will assume that l>j0, and we distinguish two cases, namely:

    Case 1. Suppose that

    E<d(xl,xj0)<E+δ2. (2.10)

    We have

    Rsf(xl,xj0)=[β1(d(xl,fxl)d(xj0,fxj0)d(xl,xj0))s+β2(d(xl,xj0))s+β3(d(xl,fxj0)+d(xj0,fxl)4)s]1/s=[β1(d(xl,xl+1)d(xj0,xj0+1)d(xl,xj0))s+β2(d(xl,xj0))s++β3(d(xl,xj0+1)+d(xj0,xl+1)4)s]1/s[β1(d(xl,xl+1)d(xj0,xj0+1)d(xl,xj0))s+β2(d(xl,xj0))s++β3(d(xl,xj0)+d(xj0,xj0+1)+d(xl,xj0)+d(xl,xl+1)4)s]1/s<[β1(d(xl,xl+1))s+β2(d(xl,xj0))s+β3(2d(xl,xj0)+d(xj0,xj0+1)+d(xl,xl+1)4)s]1/s<[β1(δ2)s+β2(E+δ2)s+β3(E2+δ4+δ4)s]1/s(β1+β2+β3)1/s(E+δ2)E+δ.

    In this case,

    E<d(xl,xj0)max{d(xl,xj0),Rsf(xl,xj0)}<max{(E+δ2),(E+δ)}=(E+δ),

    which implies by (a1) that

    α(xl,xj0)d(fxl,fxj0)E. (2.11)

    But, taking into account that the mapping f is triangular α-orbital admissible, together with (2.4) we have

    α(xn,xn+1)1 and α(xn+1,fxn+1)1 implies α(xn,xn+2)1,

    and recursively we get that

    α(xn,xl)1, (2.12)

    for all n,lN. Therefore, from (2.11) and (2.12), we have

    d(xl+1,xj0+1)=d(fxl,fxj0)E. (2.13)

    Now, by the triangle inequality together with (2.7)and (2.13) we get

    d(xl+1,xj0)d(xl+1,xj0+1)+d(xj0+1,xj0)<(E+δ2),

    which means that, indeed fxl=xl+1A.

    Case 2. Suppose that

    d(xl,xj0)E. (2.14)

    Thus,

    d(fxl,xj0)d(fxl,fxj0)+d(fxj0,xj0)α(xl,xj0)d(fxl,fxj0)+d(xj0+1,xj0)<max{d(xl,xj0),Rsf(xl,xj0)}+d(xj0+1,xj0), (2.15)

    where

    Rsf(xl,xj0)=[β1(d(xl,xl+1)d(xj0,xj0+1)d(xl,xj0))s+β2(d(xl,xj0))s++β3(d(xl,xj0)+d(xj0,xj0+1)+d(xl,xj0)+d(xl,xl+1)4)s]1/s.

    We must consider two subcases

    (2a). d(xl,xj0)d(xj0,xj0+1). Then,

    Rsf(xl,xj0)[β1(d(xl,xl+1))s+β2(d(xl,xj0))s++β3(2d(xl,xj0)+d(xj0,xj0+1)+d(xl,xl+1)4)s]1/s<[β1(δ2)s+β2(E))s+β3(2E+2δ24))s]1/s<[β1+β2+β3]1/s(E2+δ4).

    But, since δ=min{δ,E,1}, we get

    Rsf(xl,xj0)<3E4,

    and then

    d(fxl,xj0)<max{d(xl,xj0),Rsf(xl,xj0)}+d(xj0+1,xj0)<max{E,3E4}+δ2=(E+δ2),

    which shows that fxlA.

    (2b). d(xl,xj0)<d(xj0,xj0+1). Then,

    d(fxl,xj0)xl+1,xl)+d(xl,xj0)<δ2+δ2<E+δ2.

    Consequently, choosing some m,nN such that m>n>j0, we can write

    d(xm,xn)d(xm,xj0)+d(xj0,xn)<2(E+δ2)<4E,

    which leads us to

    limm,n+d(xm,xn)=0.

    Therefore, {xm} is a Cauchy sequence in a complete metric space. Thus, we can find a point uX such that limm+xm=u. Moreover, since the mapping f is continuous we have

    u=limm+fm+1x0=limm+f(fmx0)=f(limm+fmx0)=fu,

    that is, u is a fixed point of f.

    Case (B). For the case s=0, letting x=xn1 and y=xn=fxn1 in (2.2), we get

    d(xn,xn+1)α(xn1,xn)d(fxn1,fxn)<max{d(xn1,xn),Rf(xn1,xn)}, (2.16)

    where

    Rf(xn1,xn)=[d(xn1,fxn1)]β1[d(xn,fxn)]β2[d(xn1,fxn)+d(xn,fxn14]β3=[d(xn1,xn)]β1[d(xn,xn+1)]β2[d(xn1,xn+1)+d(xn,xn4]β3=[d(xn1,xn)]β1[d(xn,xn+1)]β2[d(xn1,xn+1)+d(xn,xn4]β3[d(xn1,xn)]β1[d(xn,xn+1)]β2[d(xn1,xn)+d(xn,xn+1)4]β3=[d(xn1,xn)]β1[d(xn,xn+1)]β2[d(xn1,xn)+d(xn,xn+1)4]β3

    Thus, by (2.3) and taking (2.4) into account we have

    d(xn,xn+1)α(xn1,xn)d(fxn1,fxn)<max{d(xn1,xn),Rf(xn1,xn)}.

    Now, if there exists n0N such that d(xn0,xn0+1)d(xn01,xn0), we get

    d(xn0,xn0+1)<max{d(xn01,xn0),Rf(xn01,xn0)}max{d(xn01,xn0),d(xn0,xn0+1)}<d(xn0,xn0+1),

    which is a contradiction. Therefore, d(xn,xn+1)<d(xn1,xn) for all nN, that is, the sequence {xn} decreasing and moreover, converges to some b0. Moreover, since

    Rf(xn1,xn)=[d(xn1,xn)]β1[d(xn,xn+1)]β2[d(xn1,xn+1)4]β3,

    we get that

    limn+max{d(xn1,xn),Rf(xn1,xn)}=b.

    If we suppose that b>0, then, 0<b<d(xn1,xn) and we can find δ>0 such that

    b<max{d(xn1,xn),Rf(xn1,xn)}<b+δ.

    In this way, taking E=b, we get

    b=E<max{d(xn1,xn),Rf(xn1,xn)}<E+δ,

    which implies (by (a1)) that

    d(xn1,xn)α(xn1,xn)d(fxn1,fxn))E=b,

    which is a contradiction. We thus proved that

    limmd(xn1,xn)=0. (2.17)

    We claim now, that the sequence {xn} is Cauchy. Firstly, we remark that, since d(xn1,xn)=0, there exists j0N, such that

    d(xn1,xn)<δ2, (2.18)

    for any nj0, where δ=min{δ,E,1}. Reasoning by induction, we will prove that the following relation

    d(xj0,xj0+m)<E+δ2 (2.19)

    holds, for any mN. Indeed, in case of m=1,

    d(xj0,xj0+1)<δ2<E+δ2,

    so, (2.19) is true. Now, supposing that (2.19) holds for some l, we shall show that it holds for l+1. We have

    Rf(xj0,xj0+l)=(d(xj0,fxj0))β1(xj0+l,fxj0+l)s(d(xj0,fxj0+l)+d(xj0+l,fxj0)4)β3=(d(xj0,xj0+1))β1(xj0+l,xj0+l+1)s(d(xj0,xj0+l+1)+d(xj0+l,xj0+1)4)β3(d(xj0,xj0+1))β1(d(xj0+l,xj0+l+1))s(d(xj0,xj0+l)+d(xj0+l,xj0+l+1)+d(xj0+l,xj0)+d(xj0,xj0+1)4)β3<(δ2)β1+β2((E2+δ4)+δ4)β3(E+δ2). (2.20)

    As in the Case (A), if d(xj0,xj0+l)>E, by (a2), and keeping in mind the above inequalities, we get

    E<d(xj0,xj0+l)max{d(xj0,xj0+l),Rf(xj0,xj0+l)}<max{δ2,(E+δ2)}=E+δ implies α(xj0,xj0+l)d(fxj0,fxj0+l)E.

    But, since using (2.12), it follows that

    d(xj0+1,xj0+l+1)=d(fxj0,fxj0+l)E,

    and then, by (b3) we get

    d(xj0,xj0+l+1)d(xj0,xj0+1)+d(xj0+1,xj0+l+1)<δ2+E<E+δ2.

    Therefore, (2.19) holds for (l+1). In the opposite situation, if d(xj0,xj0+l)E, again by the triangle inequality, we obtain

    d(xj0,xj0+l+1)d(xj0,xj0+1)+d(xj0+1,xj0+l+1)d(xj0,xj0+1)+α(xj0,xj0+l)d(fxj0,fxj0+l)<δ2+max{d(xj0,xj0+l),Rf(xj0,xl)}<δ2+max{E,E2+δ4}=δ2+E.

    Consequently, the induction is completed. Therefore, {xm} is a Cauchy sequence in a complete metric space. Thus, there exists uX such that fu=u.

    In the above Theorem, the continuity condition of the mapping f can be replace by the continuity of f2.

    Theorem 2.2. Suppose that f:(X,d)(X,d) forms an α-hybrid Jaggi-Meir-Keeler type contraction such that f2 is continuous.Then, f has a fixed point, provided that there exists x0X, such that α(x0,fx0)1.

    Proof. Let x0X such that α(x0,fx0)1 and the sequence {xn}, where xn=fxn1, for any nN. Thus, from Theorem 2.3 we know that this is a convergent sequence. Letting u=limn+xn, we claim that u=fu.

    Since the mapping f2 is supposed to be continuous,

    f2u=limn+f2xn=u.

    Assuming on the contrary, that ufu, we have

    Rsf(u,fu)={[β1(d(u,fu)d(fu,f2u)d(u,fu))s+β2(d(u,fu))s+β3(d(u,f2u)+d(fu,fu)4)s]1/s,fors>0(d(u,fu))β1(d(fu,f2u))β2(d(u,f2u)+d(fu,fu)4)β3,fors=0[16pt]={[β1(d(u,fu)d(fu,u)d(u,fu))s+β2(d(u,fu))s+β3(d(u,u)+d(fu,fu)4)s]1/s,fors>0(d(u,fu))β1(d(fu,u))β2(d(u,u)+d(fu,fu)4)β3,fors=0[16pt]={[β1(d(fu,u))s+β2(d(u,fu))s]1/s,fors>00,fors=0={[β1+β2]1/sd(u,fu)),fors>00,fors=0

    Example 2.1. Let X=[0,+), d:X×X[0,+), d(x,y)=|xy|, and the mapping f:XX, where

    f={12,ifx[0,1]16,ifx>1.

    We can easily observe that f is discontinuous at the point x=1, but f2 is a continuous mapping. Let also the function α:X×X[0,+),

    α(x,y)={x2+y2+1,ifx,y[0,1]ln(x+y)+1,ifx,y(1,+)1,ifx=56,y=760, otherwise ,

    and we choose β1=14,β2=12,β1=14 and s=2. The mapping f is triangular α-orbital admissible and satisfies (a2) in Definition 2.1 for any x,y[0,1], respectively for x,y(1,+). Taking into account the definition of the function α, we have more to check the case x=56, y=76. We have

    Rf(56,76)=[14(d(56,f56)d(76,f76)d(56,76))2+12(d(56,76))2+14(d(56,f76)+d(76,f56)4)2]1/2=14+1219+19=13.

    Therefore,

    α(56,76)d(f56,f76)=d(f56,f76)=d(12,16)=13<13=max{d(56,76),Rf(56,76)}.

    Moreover, since the mapping f satisfies condition (a1) for

    δ(E)={1E,forE<11,forE1,

    it follows that the assumptions of Theorem 2.3 are satisfied, and u=12 ia a fixed point of the mapping f.

    Theorem 2.3. If to the hypotheses of Theorem we add the following assumption

    α(u,v)1for anyu,vFf(X),

    then the mapping f admits an unique fixed point.

    Proof. Let uX be a fixed point of f. Supposing on the contrary, that we can find vX such that fu=uv=fv, we have

    (i) For s>0,

    Rsf=[β1(d(u,fu)d(v,fv)d(u,v))s+β2(d(u,v))s+β3(d(u,fv)+d(v,fu)4)s]1/s=[β2(d(u,v))s+β3(d(u,v)2)s]1/s(β2+β3)1/sd(u,v)d(u,v).

    Thus, taking x=u and y=v in (2.3) we get

    d(u,v)α(u,v)d(fu,fv)<max{d(u,v),Rf(u,v)}=d(u,v),

    which is a contradiction.

    (ii) For s=0,

    d(u,v)α(u,v)d(fu,fv)<max{d(u,v),Rf(u,v)}=max{d(u,v),(d(u,fu))β1(d(v,fv))β2(d(u,fv)+d(v,fu)4)β3}=d(u,v),

    which is a contradiction.

    Consequently, if there exists a fixed point of the mapping f, under the assumptions of the theorem, this is unique.

    Example 2.2. Let the set X=[1,+), d:X×X[0,+), d(x,y)=|xy|, and the mapping f:XX, where

    fx={x2+1,ifx[1,0)1,ifx[0,1]1x,ifx>1.

    Let also α:X×X[0,+) defined as follows

    α(x,y)={34,ifx,y[1,0)x2+y2+1,ifx,y[0,1]1,ifx[1,0),y[0,1]0, otherwise .

    It is easy to check that, with these chooses, f is a continuous triangular α-orbital admissible mapping and also, it follows that the mapping f satisfies the conditions (a2) from Definition (2.1). Moreover, f satisfies the condition (a1), considering δ(E)=1E in case of E<1 and δ(E)=1 for E1. Consequently, f satisfies the conditions of Theorem 2.3 and has a unique fixed point, u=0.

    In particular, for the case s=0, the continuity assumption of the mapping f can be replace by the condition (R).

    Theorem 2.4. We presume thatf:(X,d)(X,d)TαX and fulfills

    (ai) for given E>0, there exists δ>0 such that

    E<O(x,y)<E+δimpliesα(x,y)d(fx,fy)E, (2.21)

    with

    O(x,y)=max{d(x,y),(d(x,fx))β1(d(y,fy))β2(d(x,fy)+d(y,fx)4)β3},

    for all x,yX, where βi0, i=1,2,3 so that β1+β2+β3=1;

    (aii)

    α(x,y)d(fx,fy)<O(x,y). (2.22)

    The mapping f has a unique fixed point provided that:

    (α1) there exists x0X such that α(x0,fx0)1;

    (α3) α(u,v)1 for any u,vFf(X);

    (α2) if the sequence {xn} in X is such that for each nN

    α(xn,xn+1)1andlimn+xn=xX,

    then there exists a subsequence {xn(j)} of {xn} such that

    α(xn(j),x)1,for eachjN.

    Proof. Let x0X such that α(x0,fx0)1. Then, we know (following the proof of Theorem 2.3) that the sequence {xn}, with xn=fnx0 is convergent; let u=limn+xn. On the other hand, from (α2), we can find a subsequence {xn(j)} of {xn} such that

    α(xn(j),u)1, for each jN.

    Since we can suppose that d(xn(j)+1,fu)>0, from (aii) we have

    d((xn(j)+1,fu))α(xn(j),u)d(fxn(j),fu)<O(x,y)=max{d(xn(j),u),(d(xn(j),fxn(j)))β1(d(u,fu))β2(d(xn(j),fu)+d(u,fxn(j))4)β3}=max{d(xn(j),u),(d(xn(j),xn(j)+1))β1(d(u,fu))β2(d(xn(j),fu)+d(u,xn(j)+1)4)β3}.

    Letting n+ in the above inequality, we get d(u,fu)=0. Thus, fu=u.

    To proof the uniqueness, we consider that we can find another fixed point of f. From (aii), we have

    d(u,v)=d(fu,fv)α(u,v)d(fu,fv)<O(u,v)=d(u,v)<d(u,v),

    which is a contradiction. Therefore, u=v.

    Considering α(x,y)=1 in the above theorems, we can easily obtain the following result.

    Definition 2.2. A mapping f:(X,d)(X,d) is called hybrid Jaggi-Meir-Keeler type contraction on X if for all distinct x,yX we have:

    (a1) for given E>0, there exists δ>0 such that

    E<max{d(x,y),Rsf(x,y)}<E+δd(fx,fy)E; (2.23)

    (a2) whenever Rf(x,y)>0,

    d(fx,fy)<max{d(x,y),Rsf(x,y)}. (2.24)

    Corollary 2.1. Any hybrid Jaggi-Meir-Keeler type contraction f:(X,d)(X,d) possesses a unique fixed point provided that f is continuous or f2 is continuous.

    The authors declare that they have no conflicts of interest.



    [1] F. Hermeline, A finite volume method for the approximation of diffusion operators on distorted meshes, J. Comput. Phys., 160 (2000), 481–499. https://doi.org/10.1006/jcph.2000.6466 doi: 10.1006/jcph.2000.6466
    [2] J. M. Nordbotten, I. Aavatsmark, Monotonicity conditons for control volume methods on uniform parallelogram grids in homogeneous media, Comput. Geosci., 9 (2005), 61–72. https://doi.org/10.1007/s10596-005-5665-2 doi: 10.1007/s10596-005-5665-2
    [3] E. Bertolazzi, G. Manzini, A second-order maximum principle preserving finite volume method for steady convection-diffusion problems, SIAM J. Numer. Anal., 43 (2005), 2172–2199. https://doi.org/10.1137/040607071 doi: 10.1137/040607071
    [4] J. Droniou, Finite volume schemes for diffusion equations: Introduction to and review of modern methods, Math. Models Methods Appl. Sci., 24 (2014), 1575–1619. https://doi.org/10.1142/S0218202514400041 doi: 10.1142/S0218202514400041
    [5] C. L. Potier, Finite volume monotone scheme for highly anisotropic diffusion operators on unstructured triangular meshes, C. R. Math., 341 (2005) 787–792. https://doi.org/10.1016/j.crma.2005.10.010
    [6] K. Lipnikov. D. Svyatskiy, Y. Vassilevski, Monotone finite volume method for diffusion equations on unstructured triangular and shape-regular polygonal meshes, J. Comput. Phys., 227 (2007), 492–512. https://doi.org/10.1016/j.jcp.2007.08.008 doi: 10.1016/j.jcp.2007.08.008
    [7] Z. Sheng, J. Yue, G. Yuan, Monotone finite volume schemes of nonequilibrium radiation diffusion equations on distorted meshes, SIAM J. Sci. Comput., 31 (2009), 2915–2934. https://doi.org/10.1137/080721558 doi: 10.1137/080721558
    [8] K. Lipnikov, D. Svyatskiy, Y. Vassilevski, A monotone finite volume method for advection-diffusion equations on unstructured polygonal meshes, J. Comput. Phys., 229 (2010), 4017–4032. https://doi.org/10.1016/j.jcp.2010.01.035 doi: 10.1016/j.jcp.2010.01.035
    [9] K. Lipnikov, G. Manzini, D. Svyatskiy, Analysis of the monotonicity condition in the mimetic finite difference method for elliptic problems, J. Comput. Phys., 230 (2011), 2620–2642. https://doi.org/10.1016/j.jcp.2010.12.039 doi: 10.1016/j.jcp.2010.12.039
    [10] J. Camier, F. Hermeline, A monotone nonlinear finite volume method for approximating diffusion operators on general meshes, Int. J. Numer. Methods Eng., 107 (2016), 496–519. https://doi.org/10.1002/nme.5184 doi: 10.1002/nme.5184
    [11] C. L. Potier, A nonlinear finite volume scheme satisfying maximum and minimum principles for diffusion operators, Int. J. Finite Vol., 2 (2009).
    [12] Z. Sheng, G. Yuan, The finite volume scheme preserving extremum principle for diffusion equations on polygonal meshes, J. Comput. Phys., 230 (2011), 2588–2604. https://doi.org/10.1016/j.jcp.2010.12.037 doi: 10.1016/j.jcp.2010.12.037
    [13] Z. Sheng, G. Yuan, Construction of nonlinear weighted method for finite volume schemes preserving maximum principle, SIAM J. Sci. Comput., 40 (2018), A607–A628. https://doi.org/10.1137/16m1098000 doi: 10.1137/16m1098000
    [14] Y. Yu, G. Yuan, Z. Sheng, The finite volume scheme preserving maximum principle for diffusion equations with discontinuous coefficient, Comput. Math. Appl., 79 (2020), 2168–2188. https://doi.org/10.1016/j.camwa.2019.10.012 doi: 10.1016/j.camwa.2019.10.012
    [15] S. Su, J. Wu, A symmetric and coercive finite volume scheme preserving the discrete maximum principle for anisotropic diffusion equations on star-shaped polygonal meshes, Appl. Numer. Math., 198 (2024), 217–235. https://doi.org/10.1016/J.APNUM.2024.01.008 doi: 10.1016/J.APNUM.2024.01.008
    [16] J. Droniou, C. L. Potier, Construction and convergence study of schemes preserving the elliptic local maximum principle, SIAM J. Numer. Anal., 49 (2011), 459–490. https://doi.org/10.1137/090770849 doi: 10.1137/090770849
    [17] Z. Sheng, G. Yuan, Analysis of nonlinear scheme preserving maximum principle for anisotropic diffusion equation on distorted meshes, Sci. China Math., 65 (2022), 2379–2396. https://doi.org/10.1007/S11425-021-1931-3 doi: 10.1007/S11425-021-1931-3
    [18] S. Korotov, M. Krizek. P. Neittaanmäki, Weakened acute type condition for tetrahedral triangulations and the discrete maximum principle, Math. Comput., 70 (2000), 107–119. https://doi.org/10.1090/s0025-5718-00-01270-9 doi: 10.1090/s0025-5718-00-01270-9
    [19] E. Burman, A. Ern, Discrete maximum principle for Galerkin approximatios of the Laplace operator on arbitrary meshes, C. R. Math., 338 (2004), 641–646. https://doi.org/10.1016/j.crma.2004.02.010 doi: 10.1016/j.crma.2004.02.010
    [20] C. Lu, W. Huang, J. Qiu, Maximum principle in linear finite element approximations of anisotropic diffusion-convection-reaction problems, Numer. Math., 127 (2014), 515–537. https://doi.org/10.1007/s00211-013-0595-8 doi: 10.1007/s00211-013-0595-8
    [21] W. Huang, Y. Wang, Discrete maximum principle for the weak Galerkin method for anisotropic diffusion problems, Commum. Comput. Phys., 18 (2015), 65–90. https://doi.org/10.4208/cicp.180914.121214a doi: 10.4208/cicp.180914.121214a
    [22] C. Cances, M. Cathala, C. L. Potier, Monotone corrections for generic cell-centered finite volume approximations of anisotropic diffusion equations, Numer. Math., 125 (2013), 387–417. https://doi.org/10.1007/s00211-013-0545-5 doi: 10.1007/s00211-013-0545-5
    [23] H. Yang, B. Yu, Y. Li, Monotonicity correction for second order element finite volume methods of anisotropic diffusion problems, J. Comput. Phys., 449 (2022), 110759. https://doi.org/10.1016/J.JCP.2021.110759 doi: 10.1016/J.JCP.2021.110759
    [24] Z. Sheng, G. Yuan, A nonlinear scheme preserving maximum principle for heterogeneous anisotropic diffusion equation, J. Comput. Appl. Math., 436 (2024). https://doi.org/10.1016/J.CAM.2023.115438
    [25] Z. Sheng, G. Yuan, A nine point scheme for the approximation of diffusion operators on distorted quadrilateral meshes, SIAM J. Sci. Comput., 30 (2008), 1341–1361. https://doi.org/10.1137/060665853 doi: 10.1137/060665853
    [26] Z. Sheng, G. Yuan, A new nonlinear finite volume scheme preserving positivity for diffusion equations, J. Comput. Phys., 315 (2016), 182–193. https://doi.org/10.1016/j.jcp.2016.03.053 doi: 10.1016/j.jcp.2016.03.053
  • This article has been cited by:

    1. Sana Hadj Amor, Ameni Remadi, Self similarity sets via fixed point theory with lack of convexity, 2023, 37, 0354-5180, 10055, 10.2298/FIL2329055A
    2. Vo Tri, Continuous dependence on parameters of differential inclusion using new techniques of fixed point theory, 2023, 37, 0354-5180, 5469, 10.2298/FIL2316469T
    3. Mohammed Shehu Shagari, Manzuma Mustapha, Hala H. Taha, Sarah Aljohani, Nabil Mlaiki, On Combinational Contractions with Applications, 2025, 24058440, e41905, 10.1016/j.heliyon.2025.e41905
    4. Mohammed Shehu Shagari, Faryad Ali, Monairah Alansari, Akbar Azam, New views on RLC-electric circuit models via combinational contractions, 2025, 2025, 1687-2770, 10.1186/s13661-025-02068-w
    5. Sirajo Yahaya, Mohammed Shagari, Ibrahim Fulatan, Fixed points of bilateral multivalued contractions, 2024, 38, 0354-5180, 2835, 10.2298/FIL2408835Y
  • Reader Comments
  • © 2025 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(390) PDF downloads(21) Cited by(0)

Figures and Tables

Figures(15)  /  Tables(6)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog