In this work, the ADI-FDTD method with fourth-order accuracy in time for the 2-D Maxwell's equations without sources and charges is proposed. We mainly focus on energy analysis of the proposed ADI-FDTD method. By using the energy method, we derive the numerical energy identity of the ADI-FDTD method and show that the ADI-FDTD method is approximately energy-preserving. In comparison with the energy in theory, the numerical one has two perturbation terms and can be used in computation in order to keep it approximately energy-preserving. Numerical experiments are given to show the performance of the proposed ADI-FDTD method which confirm the theoretical results.
Citation: Li Zhang, Maohua Ran, Hanyue Zhang. Energy analysis of the ADI-FDTD method with fourth-order accuracy in time for Maxwell's equations[J]. AIMS Mathematics, 2023, 8(1): 264-284. doi: 10.3934/math.2023012
[1] | Aimin Xu . Some identities involving derangement polynomials and r-Bell polynomials. AIMS Mathematics, 2024, 9(1): 2051-2062. doi: 10.3934/math.2024102 |
[2] | Feng Qi, Da-Wei Niu, Bai-Ni Guo . Simplifying coefficients in differential equations associated with higher order Bernoulli numbers of the second kind. AIMS Mathematics, 2019, 4(2): 170-175. doi: 10.3934/math.2019.2.170 |
[3] | Taekyun Kim, Dae San Kim, Hyunseok Lee, Lee-Chae Jang . A note on degenerate derangement polynomials and numbers. AIMS Mathematics, 2021, 6(6): 6469-6481. doi: 10.3934/math.2021380 |
[4] | Waseem A. Khan, Abdulghani Muhyi, Rifaqat Ali, Khaled Ahmad Hassan Alzobydi, Manoj Singh, Praveen Agarwal . A new family of degenerate poly-Bernoulli polynomials of the second kind with its certain related properties. AIMS Mathematics, 2021, 6(11): 12680-12697. doi: 10.3934/math.2021731 |
[5] | Hye Kyung Kim . Note on r-central Lah numbers and r-central Lah-Bell numbers. AIMS Mathematics, 2022, 7(2): 2929-2939. doi: 10.3934/math.2022161 |
[6] | Bai-Ni Guo, Dongkyu Lim, Feng Qi . Series expansions of powers of arcsine, closed forms for special values of Bell polynomials, and series representations of generalized logsine functions. AIMS Mathematics, 2021, 6(7): 7494-7517. doi: 10.3934/math.2021438 |
[7] | Jongkyum Kwon, Patcharee Wongsason, Yunjae Kim, Dojin Kim . Representations of modified type 2 degenerate poly-Bernoulli polynomials. AIMS Mathematics, 2022, 7(6): 11443-11463. doi: 10.3934/math.2022638 |
[8] | Rabeb Sidaoui, Abdulghani Muhyi, Khaled Aldwoah, Ayman Alahmade, Mohammed Rabih, Amer Alsulami, Khidir Mohamed . A family of trivariate Gould-Hopper-Bell-Apostol-Type polynomials: Construction and properties. AIMS Mathematics, 2025, 10(6): 13251-13277. doi: 10.3934/math.2025594 |
[9] | Dojin Kim, Patcharee Wongsason, Jongkyum Kwon . Type 2 degenerate modified poly-Bernoulli polynomials arising from the degenerate poly-exponential functions. AIMS Mathematics, 2022, 7(6): 9716-9730. doi: 10.3934/math.2022541 |
[10] | Hye Kyung Kim, Dmitry V. Dolgy . Degenerate Catalan-Daehee numbers and polynomials of order r arising from degenerate umbral calculus. AIMS Mathematics, 2022, 7(3): 3845-3865. doi: 10.3934/math.2022213 |
In this work, the ADI-FDTD method with fourth-order accuracy in time for the 2-D Maxwell's equations without sources and charges is proposed. We mainly focus on energy analysis of the proposed ADI-FDTD method. By using the energy method, we derive the numerical energy identity of the ADI-FDTD method and show that the ADI-FDTD method is approximately energy-preserving. In comparison with the energy in theory, the numerical one has two perturbation terms and can be used in computation in order to keep it approximately energy-preserving. Numerical experiments are given to show the performance of the proposed ADI-FDTD method which confirm the theoretical results.
Fixed point theory is one of the most important and significant research fields in nonlinear and functional analysis since it provides some of the most useful tools to solve many problems in applied sciences and engineering, such as dynamical systems, game theory, optimization theory, the existence of solutions to integrals, and differential and matrix equations. One of the most useful results in fixed point theory is the Banach-Picard-Caccioppoli contraction principle [1,2]. It states that a mapping T on a complete metric space X onto itself has a unique fixed point provided that T is a contractive mapping, i.e., there exists a constant k∈[0,1) such that
d(Tx,Ty)≤kd(x,y) |
for all x, y∈X. Due to its importance, the result has since been studied and generalized.
Broadly, there are two manners in the generalization of fixed point theorems. One is by weakening the conditions of the contraction and the other way is by generalizing the underlying space (which is usually done by weakening the triangle inequality).
In 1962, Rakotch [3] generalized the contraction condition by defining a family of functions such that the contraction principle was still valid when the constant k was replaced by a function with suitable properties. After that, in 1969, Boyd and Wong [4] improved the result of Rakotch by introducing in the right hand member of the contraction an upper semicontinuous function φ(⋅) such that φ(t)<t, for any t>0 and, in the same year, Meir and Keeler [5] generalized the latter result by introducing a condition of weakly uniformly strict contraction. Inspired by the contraction principle, in 1968, the concept of the Kannan contraction [6] was introduced. Note that this type of contraction may not always be continuous, thus the independence of the Kannan fixed point theorem of the contraction principle. Also, note that a Kannan contraction characterizes the completeness of a metric space, as was shown in [7]. In 1972, the Chatterjea's fixed point theorem [8] was stated, while in 1971, Reich [9] generalized Banach's fixed point theorem and Kannan's fixed point theorem by putting together the terms from the right-hand side of both contractions. In the same manner, Hardy and Rogers [10] generalized, in 1973, the results of Reich and Chatterjea and one year later, in 1974, Ćirić [11] improved the Hardy-Rogers contraction. Also, one can consider different types of rational contractions, for example, in [12], researchers used an ICS function and some generalized weak contractions of Boyd-Wong-type to generalize the results from [13]. For a more perspicuous and comprehensive view of the contractive mappings that admit a unique fixed point, we address the reader to Rhoades [14] and the references therein. In 1982, Istrăţescu [15] inaugurated the class of convex contractive mappings by introducing several convexity conditions, thus bringing out new generalizations for the contraction principle. The convex contractions, too, were later subjected to generalization. For example, in [16], Miandaragh et al. defined the notions of generalized convex contractions and generalized convex contractions of order 2 and proved fixed point theorems regarding these type of mappings.
In recent years, several generalizations of classical metric spaces have been given. For instance, Bakhtin [17] and Czerwik [18] introduced the concept of b-metric spaces. With the emergence of this space, a myriad of novel results has concurrently surfaced. For example, in [19], Ali et al. defined on b-metric spaces the notions of Hardy-Rogers-type (F-α)-contractions and Hardy-Rogers-type (F-α∗)-contractions and then established fixed point theorems for these contractions. In [20], Kamran et al. introduced on the same ambiental space the concepts of Feng-Liu-type (F-α)-contractions and Feng-Liu-type (F-α∗)-contractions and proved fixed point results regarding these contractions. In [21], Shatanawi et al. used a contraction condition by means of a comparison function to prove a result regarding a unique common fixed point of two mappings.
In 2014, Kirk and Shahzad [22] defined the notion of strong b-metric spaces. In 2000, Branciari [23] defined the framework of generalized metric spaces (also known as rectangular metric spaces) and generalized metric spaces of order ν, and in the same year, Hitzler and Seda [24] introduced the notion of dislocated metric spaces. In 2014, Khojasteh et al. [25] defined the concept of θ-metric spaces. For related further results, including the metrization of such spaces, see [26]. In 2015 and 2018, Jleli and Samet introduced the generalized metric spaces [27] and F-metric spaces [28], respectively. In 2017, Kamran et al. [29] defined the concept of extended b-metric spaces, in this way generalizing the b-metric spaces. For fixed point results in the setting of this space, the reader can consult, for example, some of the following articles and references therein. In [30], Samreen et al. came up with a generalization of some of the main results from [31,32,33]. In [34], Alqahtani et al. proved fixed point theorems for two mappings that form an (α,β)-orbital-cyclic-admissible pair and obtained corollaries for (α,β)-orbital-cyclic-admissible mappings and α-orbital-admissible mappings. In [35], Abdeljawad et al. defined the concepts of Θe-contractions and Hardy-Rogers-type Θ-contractions and proved fixed point theorems for each one of them in the setting of extended b-metric spaces. Also, in [36], Shatanawi et al. introduced the notion of α−ψ-contractive mappings and proved a fixed point result for such functions. In [37], Alqahtani et al. proved some fixed point theorems for an orbitally continuous self-map T on a T-orbitally complete extended b-metric space. In [38], Mitrović et al. proved the fixed point theorems of Reich [9] and Nadler [39] in the setting of extended b-metric spaces. Huang et al. [40] and Alqahtani et al. [41] determined fixed point results for some rational type contractions and in [42], Kiran et al. generalized the Hardy-Rogers fixed point theorem [10] in the setting of extended b-metric spaces and also proved some theorems for multi-valued mappings.
This paper is organized as follows: In Section 2, preliminary concepts and notions are recalled, such as b-metric spaces, extended b-metric spaces, and the convex contractive mappings used in the main results. Also, a new type of convex contractive mapping is presented—the Ćirić-convex contraction. In Section 3, some fixed point theorems are formulated and proven in the setting of extended b-metric spaces for the contractive mappings defined in the previous section.
First, recall the definition of a b-metric space.
Definition 2.1. [17,18] Let X be a nonempty set and s≥1 be a given real number. A function d:X×X→[0,∞) is called a b-metric provided that, for all x, y, z∈X,
● d(x,y)=0 if and only if x=y,
● d(x,y)=d(y,x),
● d(x,z)≤s[d(x,y)+d(y,z)].
A pair (X,d) is called a b-metric space of constant s.
Second, recall the definition of an extended b-metric space.
Definition 2.2. [29] Let X be a nonempty set and θ: X×X→[1,∞). A function dθ: X×X→[0,∞) is called an extended b-metric if for all x, y, z∈X, it satisfies:
● dθ(x,y)=0 if and only if x=y,
● dθ(x,y)=dθ(y,x),
● dθ(x,z)≤θ(x,z)[dθ(x,y)+dθ(y,z)].
A pair (X,dθ) is called an extended b-metric space.
It is obvious that a b-metric space is a particular case of an extended b-metric space by taking θ(x,y)=s. The following example shows that the class of extended b-metric spaces is larger than the class of b-metric spaces, in the sense of inclusion.
Example 2.1. [43] Let X=[−1,1] and θ: X×X→[1,∞) be defined by
θ(x,y)=1+x2+y2x2+y2, |
if x2+y2>0 and θ(0,0)=1. Define dθ: X×X→[0,∞),
dθ(x,y)={0, if and only if x=y,1x2, if xy=0andx2+y2≠0,1x2y2, if 0≠x≠y≠0. |
Thus, dθ defines an extended b-metric on X, therefore (X,dθ) is an extended b-metric space. Note that (X,dθ) is not a b-metric space. To prove this, consider x, y∈[−1,1]∖{0} such that x≠y. We have that
dθ(x,y)dθ(x,0)+dθ(0,y)=1x2y21x2+1y2=1x2+y2. |
Note that
sup{1x2+y2:x,y∈[−1,1]∖{0},x≠y}=+∞. |
Therefore, it is impossible to find s≥1 such that
dθ(x,y)≤s[dθ(x,0)+dθ(0,y)]. |
Onwards, we recollect the concepts of convergence, Cauchy sequence, and completeness in an extended b-metric space.
Definition 2.3. [29] Let (X,dθ) be an extended b-metric space. Then a sequence {xn}n in X is said to be:
● Convergent if and only if there exists x∈X such that dθ(xn,x)→0 as n→∞, and we write limn→∞xn=x,
● Cauchy if and only if dθ(xn,xm)→0 as n, m→∞.
The extended b-metric space (X,dθ) is complete if every Cauchy sequence converges in X. We note that the extended b-metric dθ is not a continuous function in general.
A useful result in proving fixed point theorems in the setting of extended b-metric spaces is presented here.
Lemma 2.1. [30] Let (X,dθ) be an extended b-metric space. Then every convergent sequence has a unique limit.
Proof. Consider a convergent sequence {xn}n of X and presume there exist u, v∈X such that
limn→∞xn=uandlimn→∞xn=v. |
Then,
dθ(u,v)≤θ(u,v)[dθ(u,xn)+dθ(xn,v)]. |
As θ(⋅,⋅) is finite, by taking the limit when n→∞ in the previous inequality, we obtain that dθ(u,v)=0, thus u=v. Therefore, the limit of {xn}n is unique.
Remark 2.1. If we define θ(⋅,⋅) as θ: X×X→[1,∞], then the uniqueness of a convergent sequence would not yield from Definition 2.2.
Henceforth, consider X to be a nonempty set and (X,d) to be an extended b-metric space.
A usual property of functions used in proving fixed point results for convex contractive mappings is stated next.
Definition 2.4. [44] A mapping T: X→X is called orbitally continuous if
limi→∞Tnix=z |
implies
limi→∞T(Tnix)=Tz, |
where Tn denotes the n-fold composition of T with itself.
In the following, we present the types of convex contractive mappings that will be used in the main results.
Definition 2.5. [15] A mapping T: X→X is said to be a convex contraction of order 2 if there exist a, b∈[0,1) with a+b<1 such that, for all x, y∈X, the following inequality holds:
d(T2x,T2y)≤ad(Tx,Ty)+bd(x,y). |
Note that this class of mappings contains the class of contractive mappings in the sense of [1,2].
Definition 2.6. [15] A mapping T: X→X is said to be a two-sided convex contraction if there exist a1, a2, b1, b2∈[0,1) with a1+a2+b1+b2<1 such that, for all x, y∈X, the following inequality holds:
d(T2x,T2y)≤a1d(x,Tx)+a2d(Tx,T2x)+b1d(y,Ty)+b2d(Ty,T2y). |
Definition 2.7. [15] A mapping T: X→X is said to be a convex contraction of type 2 if there exist constants c0, c1, a1, a2, b1, b2∈[0,1) with c0+c1+a1+a2+b1+b2<1 such that, for all x, y∈X, the following inequality holds:
d(T2x,T2y)≤c0d(x,y)+c1d(Tx,Ty)+a1d(x,Tx)+a2d(Tx,T2x)+b1d(y,Ty)+b2d(Ty,T2y). |
Definition 2.8. [15] A mapping T: X→X is said to be a convex contraction of order k≥2 if there exist a0, a1, …, ak−1∈[0,1) with a0+a1+…+ak−1<1 such that, for all x, y∈X, the following inequality holds:
d(Tkx,Tky)≤a0d(x,y)+a1d(Tx,Ty)+…+ak−1d(Tk−1x,Tk−1y). |
Another type of convex contractive mapping is presented next.
Definition 2.9. A mapping T: X→X is said to be a Ćirić-convex contraction if there exists h∈[0,1) such that, for all x, y∈X, the following inequality holds:
d(T2x,T2y)≤hmax{d(x,y),d(Tx,Ty),d(x,Tx),d(Tx,T2x),d(y,Ty),d(Ty,T2y)}. |
A valuable lemma that will be used throughout the proofs of the theorems in this section is given below.
Lemma 3.1. Let (X,d) be an extended b-metric space, let T: X→X be a mapping and the sequence {xn}n of Picard iterations based on an initial point x0∈X, i.e., xn=Tnx0, for all n≥0. If there exists λ∈[0,1) such that
d(xn,xn+1)≤√λn−1max{d(x0,x1),d(x1,x2)} |
for all n≥0, and there exist
β<1√λ |
and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, then {xn}n is a Cauchy sequence.
Proof. Set dn=d(xn,xn+1), for all n≥0 and M=max{d0,d1}. Consider n, p, j∈N with n≥n0, p≥1 and j∈{n,n+1,…,n+p−1}. Obviously, β√λ<1.
The following estimations will justify that {xn}n is a Cauchy sequence:
d(xn,xn+p)≤θ(xn,xn+p)d(xn,xn+1)+θ(xn,xn+p)d(xn+1,xn+p)≤θ(xn,xn+p)dn+θ(xn,xn+p)θ(xn+1,xn+p)d(xn+1,xn+2)+θ(xn,xn+p)θ(xn+1,xn+p)d(xn+2,xn+p)≤…≤θ(xn,xn+p)dn+θ(xn,xn+p)θ(xn+1,xn+p)dn+1+…+θ(xn,xn+p)θ(xn+1,xn+p)⋅⋅⋅θ(xn+p−1,xn+p)dn+p−1≤β√λn−1M+β2√λnM+…+βp√λn+p−2M≤β√λn−1M[1+β√λ+…+βp−1√λp−1]=β√λn−1M1−(β√λ)p1−β√λ≤β√λn−1M11−β√λ→0, |
when n→∞.
The bounding condition imposed on the function θ(⋅,⋅) was dictated by the necessity to ensure that the Picard sequence under consideration meets the Cauchy sequence definition. Note that, in the previous lemma, if we replace the condition
j∏i=nθ(xi,xn+p)≤βj−n+1 |
with the condition
limn,m→∞θ(xn,xm)≤βorlimn,m→∞θ(xn,xm)<1√h, |
the conclusion would still be valid.
Another lemma, which will prove to be useful, is presented as follows. For the forthcoming result, which is self-evident and straightforward, we shall opt to omit the proof.
Lemma 3.2. Let (X,d) be a complete extended b-metric space and let T: X→X be an orbitally continuous mapping. If the Picard iterations sequence based on an initial point x0∈X, i.e., xn=Tnx0, for all n≥0, is a Cauchy sequence, then T has a fixed point.
From now on, we present fixed point theorems for the previously introduced convex contractive mappings.
Theorem 3.1. Let (X,d) be a complete extended b-metric space and let T: X→X be an orbitally continuous convex contraction of type 2. Suppose there exist x0∈X,
β<√1−b2c0+c1+a1+a2+b1 |
and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, where xn=Tnx0, n≥0.
Then T has a unique fixed point.
Proof. {xn}n is the Picard iteration sequence based on the initial point x0.
Set dn=d(xn,xn+1), for all n≥0,
λ=c0+c1+a1+a2+b11−b2andM=max{d0,d1}. |
Affirmation 1. dn≤√λn−1M, for all n≥0.
The statement will be proven by complete induction on n. For n=0 and n=1, the inequality is obvious. For n=2, it follows that
d2=d(x2,x3)=d(T2x0,T2x1)≤c0d(x0,x1)+c1d(x1,x2)+a1d(x0,x1)+a2d(x1,x2)+b1d(x1,x2)+b2d(x2,x3), |
then
d2(1−b2)≤(c0+a1)d0+(c1+a2+b1)d1≤(c0+a1)M+(c1+a2+b1)M=(c0+c1+a1+a2+b1)M, |
thus
d2≤c0+c1+a1+a2+b11−b2M=λM<√λM=√λ2−1M. |
Therefore, the base step is verified. Henceforth, consider k≥1 such that
dk≤√λk−1Manddk−1≤√λk−2M. |
Considering these inequalities, we present the inductive step:
dk+1=d(xk+1,xk+2)=d(T2xk−1,T2xk)≤c0dk−1+c1dk+a1dk−1+a2dk+b1dk+b2dk−1, |
then
dk+1(1−b2)≤(c0+a1)dk−1+(c1+a2+b1)dk≤(c0+a1)√λk−2M+(c1+a2+b1)√λk−1M≤(c0+a1)√λk−2M+(c1+a2+b1)√λk−2M=√λk−2M(c0+a1+c1+a2+b1), |
thus
dk+1≤√λk−2Mc0+c1+a1+a2+b11−b2=√λk−2Mλ=√λ(k+1)−1M. |
Therefore, by complete induction on n, we conclude that dn≤√λn−1M, for all n≥0.
Affirmation 2. T has a fixed point.
By using Lemma 3.1, we conclude that the sequence {xn}n is Cauchy. Also, by making use of Lemma 3.2, we get that T has a fixed point u∈X.
Affirmation 3. T has a unique fixed point.
Assume that there exists v∈X such that v≠u and Tv=v. Then,
d(u,v)=d(T2u,T2v)≤c0d(u,v)+c1d(Tu,Tv)+a1d(u,Tu)+a2d(Tu,T2u)+b1d(v,Tv)+b2d(Tv,T2v)=(c0+c1)d(u,v)<d(u,v) |
is a contradiction.
Consequently, u is the only fixed point of T.
By setting a1=a2=b1=b2=0 and c0=c1=0 in Theorem 3.1, we obtain the corresponding results for convex contractions of order 2 and two-sided convex contractions, respectively.
Corollary 3.1. Let (X,d) be a complete extended b-metric space and let T: X→X be an orbitally continuous convex contraction of order 2. If there exist x0∈X,
β<1√a+b, |
n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, where xn=Tnx0, n≥0, then T has a unique fixed point.
Corollary 3.2. Let (X,d) be a complete extended b-metric space and let T: X→X be an orbitally continuous two-sided convex contraction. Suppose there exist x0∈X,
β<√1−b2a1+a2+b1 |
and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, where xn=Tnx0, n≥0.
Then T has a fixed point that is unique.
Theorem 3.2. Let (X,d) be a complete extended b-metric space and let T: X→X be an orbitally continuous Ćirić-convex contraction. Presume there exist x0∈X,
β<1√h |
and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, where xn=Tnx0, n≥0.
Then T has a unique fixed point.
Proof. Consider {xn}n as the Picard iteration sequence based on the initial point x0.
If there exists n0∈N such that xn0=xn0+1, then xn0 is a fixed point of T. Without loss of generality, we may assume that xn≠xn+1 for any n∈N. Set dn=d(xn,xn+1), for all n≥0 and M=max{d0,d1}.
Affirmation 1. dn≤√hn−1M, for all n≥0.
The statement will be proven by complete induction on n. For n=0 and n=1, the inequality is obvious. For n=2, it follows that:
d2=d(x2,x3)=d(T2x0,T2x1)≤hmax{d0,d1,d0,d1,d1,d2}=hmax{d0,d1,d2}=hmax{M,d2}. |
If max{M,d2}=d2, then
d2≤hd2<d2, |
is a contradiction.
Thus, max{M,d2}=M. Eventually,
d2≤hM≤√hM=√h2−1M. |
Therefore, the base step is verified. Henceforth, consider k≥1 such that
dk≤√hk−1Manddk−1≤√hk−2M. |
Considering these inequalities, we present the inductive step:
dk+1=d(xk+1,xk+2)=d(T2xk−1,T2xk)≤max{dk−1,dk,dk+1}. |
If max{dk−1,dk,dk+1}=dk+1, then
dk+1≤hdk+1<dk+1, |
which is a contradiction.
Therefore,
max{dk−1,dk,dk+1}=max{dk−1,dk}. |
Then,
dk+1≤hmax{dk−1,dk}≤hmax{√hk−2M,√hk−1M}=h√hk−2M=√h(k+1)−1M. |
Thus, by complete induction on n, we conclude that dn≤√hn−1M, for all n≥0.
Affirmation 2. T has a fixed point.
By using Lemma 3.1 with λ=h, {xn}n is a Cauchy sequence. Now, employing Lemma 3.2, T has a fixed point u∈X.
Affirmation 3. T has a unique fixed point.
Assume that there exists v∈X such that v≠u and Tv=v. Then,
d(u,v)=d(T2u,T2v)≤hmax{d(u,v),d(Tu,Tv),d(u,Tu),d(Tu,T2u),d(v,Tv),d(Tv,T2v)}=hd(u,v)<d(u,v), |
which is a contradiction.
Consequently, u is the only fixed point of T.
In order to prove the next theorem, we first need to prove the following result.
Lemma 3.3. Let (X,d) be an extended b-metric space and let T: X→X be a mapping and the sequence {xn}n of Picard iterations based on an initial point x0∈X, i.e., xn=Tnx0, for all n≥0. If there exist λ∈[0,1) and an integer k≥2 such that
d(xn,xn+1)≤k√λn−kM |
for all n≥0, where
M=max{d(x0,x1),d(x1,x2),…,d(xk−1,xk)} |
and there exist β<1k√λ and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, then {xn}n is a Cauchy sequence.
Proof. Set dn=d(xn,xn+1), for all n≥0 and M=max{d0,d1}. Consider n, p, j∈N with n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}. Obviously, βk√λ<1.
The following estimations will justify that {xn}n is a Cauchy sequence.
d(xn,xn+p)≤θ(xn,xn+p)d(xn,xn+1)+θ(xn,xn+p)d(xn+1,xn+p)≤θ(xn,xn+p)dn+θ(xn,xn+p)θ(xn+1,xn+p)d(xn+1,xn+2)+θ(xn,xn+p)θ(xn+1,xn+p)d(xn+2,xn+p)≤…≤θ(xn,xn+p)dn+θ(xn,xn+p)θ(xn+1,xn+p)dn+1+…+θ(xn,xn+p)θ(xn+1,xn+p)⋅⋅⋅θ(xn+p−1,xn+p)dn+p−1≤βk√λn−kM+β2k√λn−k+1M+…+βpk√λn−k+p−1M≤βk√λn−kM[1+βk√λ+…+βp−1k√λp−1]=βk√λn−kM1−(βk√λ)p1−βk√λ≤βk√λn−kM11−βk√λ→0, |
when n→∞.
The following result is an extension of Corollary 3.1.
Theorem 3.3. Let (X,d) be a complete extended b-metric space, let k be an integer such that k≥2, and let T: X→X be an orbitally continuous convex contraction of order k. Suppose there exist x0∈X,
β<1k√a0+a1+…+ak−1 |
and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, where xn=Tnx0, n≥0.
Then T has a unique fixed point.
Proof. Let x0∈X. Construct the sequence {xn}n of Picard iterations based on the initial point x0. Set dn=d(xn,xn+1), for all n≥0. Set
λ=a0+a1+…+ak−1andM=max{d0,d1,…,dk−1}. |
Affirmation 1. dn≤k√λn−kM, for all n≥0.
The statement will be proven by complete induction on n.
For n∈{0,1,…,k−1}, the inequality is obvious. For n=k, it follows that
dk=d(xk,xk+1)=d(Tkx0,Tkx1)≤a0d(x0,x1)+a1d(x1,x2)+…+ak−1d(xk−1,xk)=a0d0+a1d1+…+ak−1dk−1≤a0M+a1M+…+ak−1M=λM≤M=λ0M=k√λk−kM. |
Thus, the base step is verified. Henceforth, consider s≥k such that dj≤k√λj−kM, for all j≤s. Considering these inequalities, we present the inductive step:
ds+1=d(xs+1,xs+2)=d(Tkxs−k+1,Tkxs−k+2)≤a0ds−k+1+a1ds−k+2+…+ak−1ds≤a0k√λs−2k+1M+a1k√λs−2k+2M+…+ak−1k√λs−kM≤a0k√λs−2k+1M+a1k√λs−2k+1M+…+ak−1k√λs−2k+1M=Mk√λs−2k+1[a0+a1+…ak−1]=Mk√λs−2k+1λ=k√λ(s+1)−kM. |
Affirmation 2. T has a fixed point.
By using Lemma 3.3, {xn}n is a Cauchy sequence. Now, employing Lemma 3.2, T has a fixed point u∈X.
Affirmation 3. T has a unique fixed point.
Assume that there exists v∈X such that v≠u and Tv=v. Then,
d(u,v)=d(Tku,Tkv)≤a0d(u,v)+a1d(Tu,Tv)+…+ak−1d(Tk−1u,Tk−1v)=a0d(u,v)+a1d(u,v)+…+ak−1d(u,v)=(a0+a1+…+ak−1)d(u,v)<d(u,v), |
which is a contradiction.
Consequently, u is the only fixed point of T.
By taking c0=a1=b1=0 and a2=b2 in Theorem 3.1, we get the following power contraction version of the Reich fixed point theorem for extended b-metric spaces.
Corollary 3.3. Let (X,dθ) be a complete extended b-metric space and let T: X→X be an orbitally continuous mapping. Suppose there exist a, b∈[0,1) with a+2b<1 such that the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+b[dθ(Tx,T2x)+dθ(Ty,T2y)] |
holds for any x, y∈X. Also, presume there exist x0∈X, β<√1−ba+b, and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, where xn=Tnx0, n≥0. Then T has a unique fixed point.
By setting c0=a1=a2=b1=b2=0 in Theorem 3.1, we obtain the following extension of the contraction principle in the framework of extended b-metric spaces.
Corollary 3.4. Let (X,dθ) be a complete extended b-metric space and let T: X→X be an orbitally continuous mapping. If there exists k∈[0,1) such that the inequality
dθ(T2x,T2y)≤kdθ(Tx,Ty) |
holds for any x, y∈X, and if there exist x0∈X, β<1√k, and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, where xn=Tnx0, n≥0, then T has a unique fixed point.
If we pick c1=a1=a2=b1=b2=0 in Theorem 3.1, we get a power contraction version of the contraction principle on extended b-metric spaces.
Corollary 3.5. Let (X,dθ) be a complete extended b-metric space and let T: X→X be an orbitally continuous mapping. Presume there exists k∈[0,1) such that the inequality
dθ(T2x,T2y)≤kdθ(x,y) |
holds for any x, y∈X, and suppose there exist x0∈X, β<1√k, and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, where xn=Tnx0, n≥0. Then, T has a fixed point, which is unique.
By choosing c0=c1=a1=b1=0 and a2=b2 in Theorem 3.1, we obtain a power contraction version of the Kannan fixed point theorem in the setting of extended b-metric spaces.
Corollary 3.6. Let (X,dθ) be a complete extended b-metric space and let T: X→X be an orbitally continuous mapping. If there exists k∈[0,12) such that the inequality
dθ(T2x,T2y)≤k[dθ(Tx,T2x)+dθ(Ty,T2y)] |
holds for any x, y∈X, and if there exist x0∈X, β<√1−kk, and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, where xn=Tnx0, n≥0, then T has a unique fixed point.
By taking θ(⋅,⋅)≡s≥1, we obtain the following analogous result for Theorem 3.1 in the setting of b-metric spaces.
Corollary 3.7. Let (X,d) be a complete b-metric space and let T: X→X be an orbitally continuous convex contraction of type 2. If
s<√1−b2c0+c1+a1+a2+b1, |
then T has a unique fixed point.
Also, we obtain the corresponding results for Corollaries 3.1 and 3.2 for b-metric spaces.
Corollary 3.8. Let (X,d) be a complete b-metric space and let T: X→X be an orbitally continuous two-sided convex contraction. If
s<√1−b2a1+a2+b1, |
then T has a fixed point that is unique.
Corollary 3.9. Let (X,d) be a complete b-metric space and let T: X→X be an orbitally continuous convex contraction of order 2. Suppose that
s<1√a+b. |
Then T has a unique fixed point.
Consequently, we also acquire the analogous results for Theorems 3.2 and 3.3 in the setting of b-metric spaces.
Corollary 3.10. Let (X,d) be a complete b-metric space and let T: X→X be an orbitally continuous Ćirić-convex contraction. Presume
s<1√h. |
Then T has a unique fixed point.
Corollary 3.11. Let (X,d) be a complete extended b-metric space, let k be an integer such that k≥2, and let T: X→X be an orbitally continuous convex contraction of order k. Suppose that
s<1k√a0+a1+…+ak−1. |
Then T has a unique fixed point.
By taking s=1, we obtain the main results from [15] and the analogue of Theorem 3.2 in the setting of metric spaces, stated next.
Corollary 3.12. Let (X,d) be a complete metric space and let T: X→X be a Ćirić-convex contraction. If T is orbitally continuous, then T has a unique fixed point.
Now, we provide an application of our proven results by presenting an example of a mapping which fulfills the conditions of a convex contractive mapping of order 2 so that it has a unique fixed point, following the work in [34]. Similar examples can be found for all the other types of convex contractive mappings presented in this paper.
Example 3.1. Let X=[0,1]. Define θ: X×X→[1,∞),
θ(x,y)={x+y+1x+y, if (x,y)≠(0,0),1.25, if (x,y)=(0,0). |
Define dθ: X×X→[0,∞),
dθ(x,y)={0, if and only if x=y,1x, if xy=0,andx2+y2≠0,1xy, if 0≠x≠y≠0. |
Define T: X→X,
Tx={2x, if x∈[0,14],0, if x∈(14,1]. |
We shall prove that the mapping T satisfies the conditions stated in Theorem 3.1, thus T has a unique fixed point.
Affirmation 1. (X,dθ) is an extended b-metric space and is not a b-metric space.
Justification 1. As it was shown in [34], dθ is an extended b-metric, thus (X,dθ) is an extended b-metric space.
Due to an argument similar to the one used in Example 2.1, we conclude that (X,dθ) is not a b-metric space.
Affirmation 2. T is a convex contraction of order 2.
Justification 2. Note that
T2x={4x, if x∈[0,18],0, if x∈(18,1]. |
Observe, by induction, that
Tkx={2kx, if x∈[0,12k+1],0, if x∈(12k+1,1]. |
To prove that T is a convex contraction of order 2, we have to find constants a, b∈[0,1) with a+b<1 such that the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
holds for any x, y∈[0,1]. We will prove that a=12 and b=18 satisfy the contractive inequality.
Let x, y∈X. We have to consider different cases.
Case 1. If x=y, note that Tx=Ty and T2x=T2y, thus
dθ(x,y)=dθ(Tx,Ty)=dθ(T2x,T2y)=0, |
and then, the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
becomes 0≤0, thus the contractive inequality holds for any a and b.
Next, we have to consider the cases where x=0 or y=0. Without loss of generality, due to the symmetry of dθ, we can assume that y=0.
Case 2. If x∈(0,18] and y=0, note that Tx=2x, T2x=4x, Ty=0, and T2y=0. Thus, the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
becomes
14x≤a2x+bx, |
or, equivalently,
1≤2a+4b, |
which holds for a=12 and b=18.
Case 3. If x∈(18,14] and y=0, note that Tx=2x, T2x=0, Ty=0, and T2y=0. Thus, the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
becomes
0≤a2x+bx, |
which holds for any a and b∈[0,1).
Case 4. If x∈(14,1] and y=0, note that Tx=0, T2x=0, Ty=0, and T2y=0. Thus, the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
becomes 0≤0, thus the contractive inequality holds for any a and b.
Henceforth, assume that 0≠x≠y≠0.
Case 5. If x, y∈(0,18], note that Tx=2x, T2x=4x, Ty=2y, and T2y=4y. Thus, the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
becomes
116xy≤a4xy+bxy, |
or, equivalently,
1≤4a+16b, |
which holds for a=12 and b=18.
Case 6. If x, y∈(18,14], note that Tx=2x, T2x=0, Ty=2y, and T2y=0. Thus, the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
becomes
0≤a4xy+bxy, |
thus the contractive inequality holds for any a and b.
Case 7. If x, y∈(14,1], note that Tx=0, T2x=0, Ty=0, and T2y=0. Thus, the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
becomes
0≤0+bxy, |
which holds for any b∈[0,1).
Due to the symmetry of dθ, it is enough to consider only three more cases.
Case 8. If x∈(0,18] and y∈(18,14], note that Tx=2x, T2x=4x, Ty=2y, and T2y=0. Thus, the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
becomes
14x≤a4xy+bxy, |
or, equivalently,
y≤a+4b. |
Note that with a=12 and b=18, we get that
y≤1=a+4b. |
Case 9. If x∈(0,18] and y∈(14,1], note that Tx=2x, T2x=4x, Ty=0, and T2y=0. Thus, the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
becomes
14x≤a2x+bxy, |
or, equivalently,
y≤2ay+4b. |
Considering a=12 and b=18, we get that
y≤y+12=2ay+4b. |
Case 10. If x∈(18,14] and y∈(14,1], note that Tx=2x, T2x=0, Ty=0, and T2y=0. Thus, the inequality
dθ(T2x,T2y)≤adθ(Tx,Ty)+bdθ(x,y) |
becomes
0≤a2x+bxy, |
which holds for any a, b∈[0,1).
Consequently, considering the cases presented above, we can conclude that T is a convex contraction of order 2 with a=12 and b=18.
Affirmation 3. There exist
β<1√a+b |
and n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤βj−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, where xn=Tnx0, for all n≥0.
Justification 3. Let p≥1.
1√a+b=1√12+18=2√105. |
Let
β=1.26<2√105. |
To prove that there exists n0∈N such that the inequality
j∏i=nθ(xi,xn+p)≤1.26j−n+1 |
holds for all n≥n0, p≥1, and j∈{n,n+1,…,n+p−1}, it is enough to prove that there exists n0∈N such that θ(xk,xn+p)<1.26 holds for all k∈{n,n+1,…,n+p−1}, k≥n0.
Let k∈{n,n+1,…,n+p−1}. If x0=0, Tnx0=0, for all n∈N. Therefore,
θ(xk,xn+p)=θ(0,0)=1.25≤1.26. |
If x0∈(0,1], let q∈N be the largest number such that x0∈(0,12q]. Set n0=q+1. Therefore, xk=Tkx0=0, for any k≥n0. As k≤n+p, obviously, xn+p=0 also. Therefore,
θ(xk,xn+p)=θ(0,0)=1.25≤1.26 |
for any k≥n0.
Therefore, all the conditions from Theorem 3.1 are satisfied.
Remark that 0 is a fixed point of T. We can conclude that 0 is the only fixed point of T.
This research article delved into the exploration of convex contractions within the realm of extended b-metric spaces, with a particular emphasis on the existence and uniqueness of fixed points under various convexity conditions. The introduction of the Ćirić-convex contraction enriches the fixed point theory by broadening the understanding of contractive mappings.
The findings underscore the versatility and potency of convex contractions in generalized metric spaces, offering novel insights and extending the boundaries of traditional fixed point theorems. This work not only contributes to the existing literature by providing generalizations and new perspectives on convex contractions, but also lays the groundwork for future research in the setting of generalized extended b-metric spaces, such as new extended b-metric spaces [45], controlled metric-type spaces [46] and double controlled metric-type spaces [47], etc.
The author declares he has not used Artificial Intelligence (AI) tools in the creation of this article.
This research was funded by the National University of Science and Technology Politehnica Bucharest [PUB Art] and the APC was funded by [PUB Art].
The author declares that he has no competing interests.
[1] |
K. S. Yee, Numerical solution of inition of initial boundary value problems involving Maxwell's equations in isotropic media, IEEE T. Antenn. Propag., 14 (1966), 302–307. http://dx.doi.org/10.1109/TAP.1966.1138693 doi: 10.1109/TAP.1966.1138693
![]() |
[2] |
T. Namiki, A new FDTD algorithm based on alternating direction implicit method, IEEE T. Microw. Theory, 47 (1999), 2003–2007. http://dx.doi.org/10.1109/22.795075 doi: 10.1109/22.795075
![]() |
[3] |
P. Ciarlet, J. Zou, Fully discrete finite element approaches for time-dependent Maxwell's equations, Numer. Math., 82 (1999), 193–219. http://dx.doi.org/10.1007/s002110050417 doi: 10.1007/s002110050417
![]() |
[4] |
L. Mu, J. Wang, X. Ye, S. Zhang, A weak Galerkin finite element method for the Maxwell's equations, J. Sci. Comput., 65 (2016), 363–386. http://dx.doi.org/10.1007/s10915-014-9964-4 doi: 10.1007/s10915-014-9964-4
![]() |
[5] |
A. Taflove, M. E. Brodwin, Numerical solution of steady-state electromagnetic scattering problems using the time-dependent Maxwell's equations, IEEE T. Microw. Theory, 23 (1975), 623–630. http://dx.doi.org/10.1109/tmtt.1975.1128640 doi: 10.1109/tmtt.1975.1128640
![]() |
[6] |
J. E. Dendy, G. Fairweather, Alternating-direction Galerkin methods for parabolic and hyperbolic problems on rectangular polygons, SIAM J. Numer. Anal., 12 (1975), 144–163. http://dx.doi.org/10.2307/2156287 doi: 10.2307/2156287
![]() |
[7] |
R. Holland, Implicit three-dimensinal finite differencing of Maxwell's equations, IEEE T. Nucl. Sci., 31 (1984), 1322–1326. http://dx.doi.org/10.1109/tns.1984.4333504 doi: 10.1109/tns.1984.4333504
![]() |
[8] | R. Leis, Initial boundary value problems in mathematical physics, New York, Wiley, 1986. http://dx.doi.org/10.1007/978-3-663-10649-4 |
[9] |
H. Yoshida, Construction of higher order symplectic integrators, Phys. Lett. A, 150 (1990), 262–268. http://dx.doi.org/10.1016/0375-9601(90)90092-3 doi: 10.1016/0375-9601(90)90092-3
![]() |
[10] |
S. Li, L. Vu-Quoc, Finite difference calculus invariant structure of a class of algorithms for the nonlinear Klein-Gordon equation, SIAM J. Numer. Anal., 32 (1995), 1839–1875. http://dx.doi.org/10.2307/2158531 doi: 10.2307/2158531
![]() |
[11] | A. Taflove, S. Hagness, Computational electrodynamics: The finite-difference time-domain method, 2Eds., Boston, Artech House, 2000. http://dx.doi.org/10.1016/0021-9169(96)80449-1 |
[12] |
F. Zheng, Z. Chen, J. Zhang, Toward the development of a three-dimensional unconditionally stable finite-difference time-domain method, IEEE T. Microw. Theory, 48 (2000), 1550–1558. http://dx.doi.org/10.1109/22.869007 doi: 10.1109/22.869007
![]() |
[13] |
T. Namiki, K. Ito, Investigation of numerical errors of the two-dimensional ADI-FDTD method, IEEE T. Microw. Theory, 48 (2000), 1950–1956. http://dx.doi.org/10.1109/22.883876 doi: 10.1109/22.883876
![]() |
[14] |
F. Zheng, Z. Chen, Numerical dispersion analysis of the unconditionally stable 3-D ADI-FDTD method, IEEE T. Microw. Theory, 49 (2001), 1006–1009. http://dx.doi.org/10.1109/22.920165 doi: 10.1109/22.920165
![]() |
[15] |
S. D. Gedney, G. Liu, J. A. Roden, A. Zhu, Perfectly matched layer media with CFS for an unconditional stable ADI-FDTD method, IEEE T. Antenn. Propag., 49 (2001), 1554–1559. http://dx.doi.org/10.1109/8.964091 doi: 10.1109/8.964091
![]() |
[16] |
J. Douglas, S. Kim, Improved accuracy for locally one-dimensional methods for parabolic equations, Math. Mod. Method. Appl. S., 11 (2001), 1563–1579. http://dx.doi.org/10.1142/s0218202501001471 doi: 10.1142/s0218202501001471
![]() |
[17] |
A. P. Zhao, Analysis of the numerical dispersion of the 2-D alternating-direction implicit FDTD method, IEEE T. Microw. Theory, 50 (2002), 1156–1164. http://dx.doi.org/10.1109/22.993419 doi: 10.1109/22.993419
![]() |
[18] |
Z. Q. Xie, C. H. Chan, B. Zhang, An explict four-order staggered finite-difference time-domain method for Maxwell's equations, J. Comput. Appl. Math., 147 (2002), 75–98. http://dx.doi.org/10.1016/S0377-0427(02)00394-1 doi: 10.1016/S0377-0427(02)00394-1
![]() |
[19] |
L. Gao, B. Zhang, D. Liang, The splitting finite-difference time-domain methods for Maxwell's equations in two dimensions, J. Comput. Appl. Math., 205 (2007), 207–230. http://dx.doi.org/10.1016/j.cam.2006.04.051 doi: 10.1016/j.cam.2006.04.051
![]() |
[20] |
B. Li, W. Sun, Unconditional convergence and optimal error estimates of a Galerkin-mixed FEM for incompressible miscible flow in porous media, SIAM J. Numer. Anal., 51 (2013), 1959–1977. http://dx.doi.org/10.1137/120871821 doi: 10.1137/120871821
![]() |
[21] |
J. Wang, A new error analysis of Crank-Nicolson Galerkin FEMs for a generalized nonlinear Schrödinger equation, J. Sci. Comput., 60 (2014), 390–407. http://dx.doi.org/10.1007/s10915-013-9799-4 doi: 10.1007/s10915-013-9799-4
![]() |
[22] |
D. Li, J. Wang, Unconditionally optimal error analysis of Crank-Nicolson Galerkin FEMs for a strongly nonlinear parabolic system, J. Sci. Comput., 72 (2017), 892–915. http://dx.doi.org/10.1007/s10915-017-0381-3 doi: 10.1007/s10915-017-0381-3
![]() |
[23] |
E. L. Tan, Y. H. Ding, ADI-FDTD Method with fourth order accuracy in time, IEEE T. Microw. Theory, 18 (2008), 296–298. http://dx.doi.org/10.1109/lmwc.2008.922099 doi: 10.1109/lmwc.2008.922099
![]() |
[24] |
W. Chen, X. Li, D. Liang, Energy-conserved splitting FDTD methods for Maxwell's equations, Numer. Math., 108 (2008), 445–485. http://dx.doi.org/10.1007/s00211-007-0123-9 doi: 10.1007/s00211-007-0123-9
![]() |
[25] |
W. Chen, X. Li, D. Liang, Symmetric energy-conserved splitting FDTD scheme for the Maxwell's equations, Commun. Comput. Phys., 6 (2009), 804–825. http://dx.doi.org/10.4208/cicp.2009.v6.p804 doi: 10.4208/cicp.2009.v6.p804
![]() |
[26] |
W. Li, D. Liang, Symmetric energy-conserved S-FDTD scheme for two-dimensional Maxwell's equations in negative index metamaterials, J. Sci. Comput., 69 (2016), 696–735. http://dx.doi.org/10.1007/s10915-016-0214-9 doi: 10.1007/s10915-016-0214-9
![]() |
[27] |
M. H. Carpenter, D. Gottlieb, S. Abarbanel, Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes, J. Comput. Phys., 111 (1994), 220–236. http://dx.doi.org/10.1006/jcph.1994.1057 doi: 10.1006/jcph.1994.1057
![]() |
[28] |
D. Appel¨o, V. A. Bokil, Y. Cheng, F. Li, Energy stable SBP-FDTD methods for Maxwell-Duffing models in nonlinear photonics, IEEE J. Multiscale Mu., 4 (2019), 329–336. http://dx.doi.org/10.1109/jmmct.2019.2959587 doi: 10.1109/jmmct.2019.2959587
![]() |
[29] |
N. Sharan, P. T. Brady, D. Livescu, Time stability of strong boundary conditions in finite-difference schemes for hyperbolic systems, SIAM J. Numer. Anal., 60 (2022), 1331–1362. http://dx.doi.org/10.1137/21M1419957 doi: 10.1137/21M1419957
![]() |
[30] |
D. Liang, Q. Yuan, The spatial fourth-order energy-conserved S-FDTD scheme for Maxwell's equations, J. Comput. Phys., 243 (2013), 344–364. http://dx.doi.org/10.1016/j.jcp.2013.02.040 doi: 10.1016/j.jcp.2013.02.040
![]() |
[31] |
L. Gao, X. Li, W. Chen, New energy identities and super convergence analysis of the energy conserved splitting FDTD methods for 3D Maxwell's equations, Math. Meth. Appl. Sci., 36 (2013), 440–455. http://dx.doi.org/10.1002/mma.2605 doi: 10.1002/mma.2605
![]() |
[32] |
L. Gao, M. Cao, R. Shi, H. Guo, Energy conservation and super convergence analysis of the EC-S-FDTD schemes for Maxwell's equations with periodic boundaries, Numer. Meth. Part. D. E., 35 (2019), 1562–1587. http://dx.doi.org/10.1002/num.22364 doi: 10.1002/num.22364
![]() |
[33] |
J. Xie, D. Liang, Z. Zhang, Energy-preserving local mesh-refined splitting FDTD schemes for two dimensional Maxwell's equations, J. Comput. Phys., 425 (2021), 109896. http://dx.doi.org/10.1016/j.jcp.2020.109896 doi: 10.1016/j.jcp.2020.109896
![]() |
[34] |
D. Wang, A. Xiao, W. Yang, A linearly implicit conservative difference scheme for the space fractional coupled nonlinear Schrödinger equations, J. Comput. Phys., 272 (2014), 644–655. http://dx.doi.org/10.1016/j.jcp.2014.04.047 doi: 10.1016/j.jcp.2014.04.047
![]() |
[35] |
P. Wang, C. Huang, An energy conservative difference scheme for the nonlinear fractional Schrödinger equations, J. Comput. Phys., 293 (2015), 238–251. http://dx.doi.org/10.1016/j.jcp.2014.03.037 doi: 10.1016/j.jcp.2014.03.037
![]() |
[36] |
M. Ran, C. Zhang, A conservative difference scheme for solving the strongly coupled nonlinear fractional Schrödinger equations, Commun. Nonlinear Sci., 41 (2016), 64–83. http://dx.doi.org/10.1016/j.cnsns.2016.04.026 doi: 10.1016/j.cnsns.2016.04.026
![]() |
[37] |
M. Ran, C. Zhang, A linearly implicit conservative scheme for the fractional nonlinear Schrödinger equation with wave operator, Int. J. Comput. Math., 93 (2016), 1103–1118. http://dx.doi.org/10.1080/00207160.2015.1016924 doi: 10.1080/00207160.2015.1016924
![]() |
[38] |
H. Li, Y. Wang, Q. Sheng, An energy-preserving Crank-Nicolson Galerkin method for Hamiltonian partial differential equations, Numer. Meth. Part. D. E., 32 (2016), 1485–1504. http://dx.doi.org/10.1002/num.22062 doi: 10.1002/num.22062
![]() |
1. | Rashad A. Al-Jawfi, Abdulghani Muhyi, Wadia Faid Hassan Al-shameri, On Generalized Class of Bell Polynomials Associated with Geometric Applications, 2024, 13, 2075-1680, 73, 10.3390/axioms13020073 |