
The non-Fourier heat transfer in heterogeneous media is crucial for material science and biomedical engineering. The optimized Schwarz waveform relaxation (OSWR) method is an efficient approach for solving such problems due to its divide-and-conquer strategy. Despite the wave-type nature of non-Fourier heat transfer, the short phase-lag time leads to more parabolic-like behavior. To address this, in the OSWR method, we employed Robin boundary conditions to transmit information along the interface. Using Fourier analysis, we derived and rigorously optimized the convergence factors of the OSWR algorithm with scaled Robin and Robin-Robin transmission conditions. The resulting optimized transmission parameters were provided in explicit form for direct application in the OSWR algorithm, along with corresponding convergence factor estimates. Interestingly, the results show that a larger heterogeneity contrast actually accelerates the convergence, rather than deteriorating it. Furthermore, the OSWR algorithm with the Robin-Robin condition exhibits mesh-independent convergence asymptotically. However, the presence of the phase-lag time is found to slow down the convergence of the OSWR algorithm. These theoretical findings were validated through numerical experiments.
Citation: Feng Hu, Yingxiang Xu. Optimized Schwarz waveform relaxation for heterogeneous Cattaneo-Vernotte non-Fourier heat transfer[J]. AIMS Mathematics, 2025, 10(3): 7370-7395. doi: 10.3934/math.2025338
[1] | Zhi-Yu Shi, Jia-Bao Liu . Topological indices of linear crossed phenylenes with respect to their Laplacian and normalized Laplacian spectrum. AIMS Mathematics, 2024, 9(3): 5431-5450. doi: 10.3934/math.2024262 |
[2] | Jahfar T K, Chithra A V . Central vertex join and central edge join of two graphs. AIMS Mathematics, 2020, 5(6): 7214-7233. doi: 10.3934/math.2020461 |
[3] | Zenan Du, Lihua You, Hechao Liu, Yufei Huang . The Sombor index and coindex of two-trees. AIMS Mathematics, 2023, 8(8): 18982-18994. doi: 10.3934/math.2023967 |
[4] | Jia-Bao Liu, Kang Wang . The multiplicative degree-Kirchhoff index and complexity of a class of linear networks. AIMS Mathematics, 2024, 9(3): 7111-7130. doi: 10.3934/math.2024347 |
[5] | Akbar Ali, Sadia Noureen, Akhlaq A. Bhatti, Abeer M. Albalahi . On optimal molecular trees with respect to Sombor indices. AIMS Mathematics, 2023, 8(3): 5369-5390. doi: 10.3934/math.2023270 |
[6] | Kun Wang, Wenjie Ning, Yuheng Song . Extremal values of the modified Sombor index in trees. AIMS Mathematics, 2025, 10(5): 12092-12103. doi: 10.3934/math.2025548 |
[7] | Zhen Lin, Ting Zhou, Xiaojing Wang, Lianying Miao . The general Albertson irregularity index of graphs. AIMS Mathematics, 2022, 7(1): 25-38. doi: 10.3934/math.2022002 |
[8] | Zhibin Du, Ayu Ameliatul Shahilah Ahmad Jamri, Roslan Hasni, Doost Ali Mojdeh . Maximal first Zagreb index of trees with given Roman domination number. AIMS Mathematics, 2022, 7(7): 11801-11812. doi: 10.3934/math.2022658 |
[9] | Tariq A. Alraqad, Igor Ž. Milovanović, Hicham Saber, Akbar Ali, Jaya P. Mazorodze, Adel A. Attiya . Minimum atom-bond sum-connectivity index of trees with a fixed order and/or number of pendent vertices. AIMS Mathematics, 2024, 9(2): 3707-3721. doi: 10.3934/math.2024182 |
[10] | Zhenhua Su, Zikai Tang, Hanyuan Deng . Higher-order Randić index and isomorphism of double starlike trees. AIMS Mathematics, 2023, 8(12): 31186-31197. doi: 10.3934/math.20231596 |
The non-Fourier heat transfer in heterogeneous media is crucial for material science and biomedical engineering. The optimized Schwarz waveform relaxation (OSWR) method is an efficient approach for solving such problems due to its divide-and-conquer strategy. Despite the wave-type nature of non-Fourier heat transfer, the short phase-lag time leads to more parabolic-like behavior. To address this, in the OSWR method, we employed Robin boundary conditions to transmit information along the interface. Using Fourier analysis, we derived and rigorously optimized the convergence factors of the OSWR algorithm with scaled Robin and Robin-Robin transmission conditions. The resulting optimized transmission parameters were provided in explicit form for direct application in the OSWR algorithm, along with corresponding convergence factor estimates. Interestingly, the results show that a larger heterogeneity contrast actually accelerates the convergence, rather than deteriorating it. Furthermore, the OSWR algorithm with the Robin-Robin condition exhibits mesh-independent convergence asymptotically. However, the presence of the phase-lag time is found to slow down the convergence of the OSWR algorithm. These theoretical findings were validated through numerical experiments.
Throughout this article, we only consider simple, undirected, and finite graphs and assume that all graphs are connected. Suppose G is a graph with the vertex set V(G)={v1,v2,⋅⋅⋅,vn} and the edge set E(G)={e1,e2,⋅⋅⋅,em}. The adjacency matrix A(G) is a 0−1n×n matrix indexed by the vertices of G and defined by aij=1 if and only if vsvt∈EG. For more notation, one refer to [1].
The Laplacian matrix of graph G is defined as L(G)=D(G)−A(G), and assume that the eigenvalues of L(G) are labeled 0≤μ1<μ2<⋅⋅⋅<μn.
(L(G))st={1,s=t;−1,s≠t,vs∼vt;0,otherwise. | (1.1) |
The normalized Laplacian matrix is given by
(L(G))st={1,s=t;−1√dsdt,s≠t,vs∼vt;0,otherwise. | (1.2) |
The distance, dij=dG(vs,vt), between vertices us and ut of G is the length of the shortest us, ut-path in G. The Wiener index [2,3] is the sum of the distances of any two vertices in the graph G, that is
W(G)=∑s<tdst. |
In 1947, the distance-based invariant first appeared in chemistry [2,3] and began to the applied to mathematics 30 years later [5]. Nowadays, the Wiener index is widely used in mathematics [6,7,8] and chemistry [9,10,11].
In a simple graph G, the degree, dij=dG(vi), of a vertex vi is the number of edges at vi. The Gutman index [12] of the simple graph G is expressed by
Gut(G)=∑s<tdsdtdst. |
Klein and Randic initially outlined the concepts associated with the resistance distance of the graph. Assume that each edge is replaced by a unit resistor, and we use rst to denote the resistance distance between two vertices s and t. Similar to the Wiener index, the Kirchhoff index [13,14] of graph G is expressed as the sum of the resistance distances between each two vertices, that is
Kf(G)=∑s<trst. |
In 2007, Chen and Zhang [15] defined the multiplicative degree-Kirchhoff index [16,17], that is
Kf∗(G)=∑s<tdsdtrst. |
Phenyl is a conjugated hydrocarbon, and L6,4,4n denote a linear chain, containing n hexagons and 2n−1 squares, as shown in Figure 1.
With the rapid changes of the times, organic chemistry has also developed rapidly, which has led to a growing interest in polycyclic aromatic compounds.
In 1985, the computational method and procedure of the matrix decomposition theorem were proposed by Yang [18]. This led to the solution of some problems in graph networks and allowed the unprecedented development of self-homogeneous linear hydrocarbon chains. For example, in 2021, X.L. Ma [20] got the normalized Laplacian spectrum of linear phenylene, and the linear phenylene containing it has n hexagons and n−1 squares. L. Lan [21] explored linear phenylene with n hexagons and n squares. Umar Ali [22] analyzed the strong prism of a graph G, which is the strong product of the complete graph of order 2, and X.Y. Geng [23] obtained the Laplacian spectrum of L6,4,4n, which contains n hexagons and 2n−1 squares. J.B. Liu [24] derived the Kirchhoff index and complexity of On, which denoting linear octagonal-quadrilateral networks. C. Liu [25] got the Laplacian spectrum and Kirchhoff index of Ln, and Ln has t hexagons and 3t+1 quadrangles.
Inspired by these recent works, we try to investigate the Laplacians and the normalized Laplaceians for graph transformations on phenyl dicyclobutadieno derivatives.
The various sections of this article are as follows: In Section 2, we proposed some concepts and lemmas and used them in subsequent articles. In Section 3 and Section 4, we acquired the Laplacian matrix and the normalized Laplacian matrix, then we made sure the Kirchoff index, the multiplicative degree-Kirchoff index, and the complexity of Ln. In Section 5, we obtained conclusions based on the calculations in this paper.
In this article, graph Ln and graph L6,4,4n are portrayed in Figure 1. Define the characteristic polynomial of matrix U of order n as PU(x)=det(xI−U).
It is easy to understand that π=(1,1′)(2,2′)⋯((4n),(4n)′) is an automorphism. Let V1={u1,u2,⋯,u4n+1,v1,⋯,v4n},V2={u′1,u′2,⋯,u′4n,v′1,⋯,v′4n+1}, |V(Ln)|=8n and |E(Ln)|=19n−4. Thus, the (normalized) Laplacians matrix can be expressed in the form of a block matrix, that is
L(Ln)=(LV0V0LV0V1LV0V2LV1V0LV1V1LV1V2LV2V0LV2V1LV2V2), |
where LVsVt and LVsVt are a submatrix consisting of rows corresponding to the vertices in Vs and columns corresponding to the vertices in Vt, s,t=0,1,2. Let
Q=(It0001√2I4n1√2I4n01√2I4n−1√2I4n), |
then
QL(LG)Q′=(LA(G)00LS(G)),QL(LG)Q′=(LA(G)00LS(G)), |
and Q′ is the transposition of Q.
LA=LV1V1+LV1V2, LS=LV1V1−LV1V2, LA=LV1V1+LV1V2, LS=LV1V1−LV1V2
Lemma 2.1. [20] If G is a graph and suppose that LA(G),LS(G),LA(G), and LS(G) are determined as above, then
ϑL(Ln)(y)=θLA(G)(y)θLS(G)(y),ϑL(Ln)(y)=θLA(G)(y)θLS(G)(y). |
Lemma 2.2. [26] With the extensive study of the Kirchhoff index, Gutman and Mohar proposed an algorithm based on the relation between the Kirchhoff index and the Laplacian eigenvalues, namely
Kf(G)=nn∑t=21ξt, |
and the eigenvalues of L(G) are 0=ξ1<ξ2≥⋯≥ξn(n≥2).
Lemma 2.3. [14] Suppose that the eigenvalues of L(G) are ξ1≤ξ2≤⋯≤ξn, then its multiplicative degree-Kirchhoff index can be denoted by
Kf∗(G)=2mn∑t=21ξt. |
Lemma 2.4. [1] The number of spanning trees in G can also be called the complexity of G. If G is a graph with |VG|=n and |EG|=m. Let λi(i=2,3,⋯,n) be the eigenvalues of L(G). Then the complexity of G is
2mτ(G)=n∏i=1di⋅n∏i=2λi. |
In this section, the main objective is to find out the Kirchhoff index of Ln. Then, combining the definition of the Laplacian matrix and Eq (1.1), we can write these block matrices as follows:
LV1V1=(3−1−14−1−15−1−15−1−15−1−14−1⋱−15−1−14−1−15−1−13)(4n)×(4n), |
and
LV1V2=(−1−1−10−1−1−1−1−1−1−1−1−1−1−10−1⋱−1−1−1−10−1−1−1−1−1−1)(4n)×(4n). |
Hence,
LA=(2−2−24−2−24−2−24−2−24−2−24−2⋱−24−2−24−2−24−2−22)(4n)×(4n), |
and
LS=diag(4,4,6,6,6,4,⋅⋅⋅,6,4,6,4). |
Assume that 0≤α1<α2≤α3≤⋯≤α4n are the roots of PLA(x)=0, and 0≤β1≤β2≤β3≤⋯≤β4n are the roots of PLS(x)=0. By Lemma 2.2, we immediately have
Kf(L2n)=2(4n)(4n∑i=21αi+4n∑j=11βj). | (3.1) |
Obviously, ∑4nj=11βj can be obtained according to LS.
4n∑j=11βj=16×(3n−2)+14×(n+2)=9n+212. | (3.2) |
Next, we focus on computing ∑4ni=11αi. Let
PLA(x)=det(xI−LA)=x(x4n−1+a1x4n−2+⋯+a4n−2x+a4n−1),a4n−1≠0. |
Based on Vieta's theorem of PLA(x), we can exactly get the following equation:
4n∑i=21αi=(−1)4n−2a4n−2(−1)4n−1a4n−1. | (3.3) |
For the sake of convenience, let Ms be used to express the s−th order principal minors of matrix A, and ms=detMs is recorded. We can get m1=2,m2=4,m3=8.
And
ms=4ms−1−4ms−2,4≤s≤4n, |
by further induction, we have
ms=2s. |
In this way, we can get two theorems.
Theorem 3.1. (−1)4n−1a4n−1=(4n)24n−1.
Proof. Due to the sum of all the principal minors of order 4n−1 of LA is (−1)4n−1a4n−1,
(−1)4n−1a4n−1=4n∑s=1detLA[s]=4n∑s=1det(Ms−100U4n−s)=4n∑s=1detMs−1⋅detU4n−s, |
where
Ms−1=(l11−2⋯0−2l22⋯0⋮⋮⋱⋮00⋯ls−1,s−1)(s−1)×(s−1), |
U4n−s=(ls+1,s+1⋯⋯0⋮⋱⋮⋮⋮⋯l4n−1,4n−1⋮0⋯−2l4n,4n)(4n−s)×(4n−s). |
Let m0=1,detU0=1, because of the symmetry of matrix LA, then detU4n−s=detM4n−s. Hence
(−1)4n−1a4n−1=4n∑s=1detms−1⋅detm4n−s=(4n)24n−1, |
as desired.
Theorem 3.2. (−1)4n−2a4n−2=(4n−1)(4n)(4n+1)24n−13.
Proof. Since the (−1)4n−2a4n−2 is the tatal of all the principal minors of order 4n−2 of LA, we have
(−1)4n−1a4n−1=∑1≤s<t≤4ndetLA[s,t]=∑1≤s<t≤4ndetMs−1detNt−s−1detU4n−t, |
where
LA[s,t]=(Mp−1000Nt−s−1000U4n−t),1≤s<t≤4n |
and
Nt−s−1=|4−2−24−2−24−2⋱−24−2−24−2−24|(t−s−1)=(t−s)2t−s−1. |
Therefore, we can have
(−1)4n−2a4n−2=∑1≤s<t≤4ndetMs−1⋅detNt−s−1⋅detU4n−t=∑1≤s<t≤4n(t−s)2t−s−1⋅detms−1⋅m4n−t=(4n−1)(4n)(4n+1)24n−13. |
The proof is over.
From the results of Theorem 3.1 and Theorem 3.2, we can get
4n∑i=21αi=(−1)4n−2a4n−2(−1)4n−1a4n−1=16n2−112, | (3.4) |
where the eigenvalues of LA are 0≤α1<α2≤α3≤⋯≤α4n.
Lemma 3.3. Suppose L6,4,4n be the dicyclobutadieno derivative of phenylenes and the graph Ln be obtained from the transformation of the graph L6,4,4n.
Kf(Ln)=32n3+18n2+2n3. |
Proof. Substituting Eqs (3.5) and (3.6) into (3.4), the Kirchhoff index of Ln can be expressed
Kf(Ln)=2(4n)(4n∑i=21αi+4n∑j=11βj)=(8n)(9n+212+(4n+1)(4n−1)12)=32n3+18n2+2n3. |
The result is as desired.
The Kirchhoff index of Ln from L1 to L12 is shown in Table 1.
G | Kf(G) | G | Kf(G) | G | Kf(G) |
L1 | 17.3 | L5 | 1486.7 | L9 | 8268.0 |
L2 | 110.7 | L6 | 2524.0 | L10 | 24457.3 |
L3 | 344.0 | L7 | 3957.3 | L11 | 30454.7 |
L4 | 781.3 | L8 | 5850.7 | L12 | 37360.0 |
Next, we will further consider the Wiener index of Ln.
Lemma 3.4. Let L6,4,4n be the dicyclobutadieno derivative of [n]phenylenes, and the graph Ln be obtained from the transformation of the graph L6,4,4n, then
limn→∞Kf(Ln)W(Ln)=14. |
Proof. Consider dst for all vertices. For the sake of convenience, we divide the vertices of the graph into the following five categories:
Case 1. Vertex 1 of Ln:
g1(i)=1+2(4n−1∑k=1k). |
Case 2. Vertex 4j−3(j=1,2,⋯,n) of Ln, i=4j−3:
g2(i)=1+2(4n−1∑k=1k+4n−i∑k=1k). |
Case 3. Vertex 4j−2(j=1,2,⋯,n) of Ln, i=4j−2:
g3(i)=1+2(4n−1∑k=1k+4n−i∑k=1k). |
Case 4. Vertex 4j−1(j=1,2,⋯,n) of Ln, i=4j−1:
g4(i)=1+2(4n−1∑k=1k+4n−i∑k=1k). |
Case 5. Vertex 4j(j=1,2,⋯,n) of Ln, i=4j:
g5(i)=1+2(4n−1∑k=1k+4n−i∑k=1k). |
Hence, we have
W(Ln)=4g1(i)+2∑i=4j−3g2(i)+2∑i=4j−2g3(i)+2∑i=4j−1g4(i)+2∑i=4jg5(i)2=4(1+2∑4n−1k=1k)+2∑nj=1[1+2(∑4j−kk=1k+∑4n−4j+2k=1k)]2+2∑nj=1[2+2(∑4j−3k=1k+∑4n−4j+2k=1k)]+2∑nj=1[1+2(∑4j−2k=1k+∑4n−4j+1k=1k)]2+2∑n−1j=1[1+2(∑4j−1k=1k+∑4n−4jk=1k)]2=128n3+48n2−5n+33. |
Consider the above results of the Kirchhoff index and the Wiener index. We can get following equation when n tends to infinity:
limn→∞Kf(Ln)W(Ln)=14. |
The result is as desired.
In this section, we use the eigenvalues of the normalized Laplacian matrix to determine the multiplicative degree-Kirchhoff index of Ln. Besides, we calculate the complexity of Ln. Then
LV1V1=(1−1√12−1√121−1√20−1√201−15−151−15−151−1√20−1√201−1√20⋱−151−1√20−1√201−1√20−1√201−1√15−1√151)(4n)×(4n), |
and
LV1V2=(−13−1√12−1√120−1√20−1√20−15−15−15−15−15−15−15−1√20−1√200−1√20⋱−15−15−1√20−1√200−1√20−1√20−15−1√15−1√15−13)(4n)×(4n). |
Therefore,
LA=(23−1√3−1√31−1√5−1√545−25−2545−25−2545−1√5−1√51−1√5⋱−2545−1√5−1√51−1√5−1√545−2√15−2√1523)(4n)×(4n), |
and
LS=diag(43,1,65,65,65,⋅⋅⋅,65,1,65,43). |
Assume that the roots of PLA(x)=0 are 0≤ξ1<ξ2≤ξ3≤⋯≤α4n, and 0≤γ1≤γ2≤γ3≤⋯≤γ4n are the roots of PLS(x)=0. By Lemma 2.3, we can get
Kf∗(L2n)=2(19n−4)(4n∑i=21ξi+4n∑j=11γj). |
Since Ls is a diagonal matrix. Obviously, its diagonal elements 1,43 and 65 correspond to the eigenvalues of Ls, respectively. Then, it can be clearly obtained as
4n∑i=11γi=21n−16. | (4.1) |
Let
PLA(x)=det(xI−LA)=x4n+b1x4n−1+⋯+b4n−1x,b4n−1≠0, |
1ξ2,1ξ3,⋅⋅⋅,1ξ4n+1are the roots of the following equation
b4n−1x4n−1+b4n−2x4n−2+⋅⋅⋅+b1x+1=0. |
Based on the Vieta' s theorem of PLA(x), we can get
4n∑i=21ξi=(−1)4n−2b4n−2(−1)4n−1b4n−1. |
Similarly, we can get another two interesting facts.
Theorem 4.1. (−1)4n−1b4n−1=259(38n−8)(4125)n.
Proof. Let sp=detFp, then we have s1=23,s2=13,s3=215,s4=475, and
{s4p=45s4p−1−425s4p−2;s4p+1=45s4p−425s4p−1;s4p+2=s4p+1−15s4p;s4p+3=45s4p+2−15s4p+1. |
After further simplification, the transformation form of the above formula is obtained.
{s4p=53⋅(4125)p,1≤p≤n;s4p+1=23⋅(4125)i,0≤p≤n−1;s4p+2=13⋅(4125)i,0≤p≤n−1;s4p+3=115⋅(4125)i,0≤p≤n−1. |
Similarly, we have t1=23,t2=415,t3=215,t4=475, and
{t4p=25t4p−1−25t4p−2;t4p+1=45t4p−425t4p−1;t4p+2=45t4p+1−425t4p;t4p+3=t4p+2−15t4p+1. |
Therefore, the transformation form of the above formula is obtained.
{t4p−4=53⋅(4125)p,1≤p≤n;t4p−3=23⋅(4125)p,0≤p≤n−1;t4p−2=415⋅(4125)p,0≤p≤n−1;t4p−1=215⋅(4125)p,0≤p≤n−1. |
Since the (−1)4n−1b4n−1 is the total of all the principal minors of order 4n−1 of LA, we have
(−1)4n−1b4n−1=4n∑i=2detNLA[i]+s4n+t4n=145(38n−8)(4125)n. |
The proof of Theorem 4.1 completed.
Theorem 4.2. (−1)4n−2b4n−2=13240(14520n3+4599n2−1496n+3)(4125)n.
Proof. We observe that the sum of all the principal minors of order 4n in LA is (−1)4n−2b4n−2, then
(−1)4n−2b4n−2=∑1≤s<t≤4ndetLA[s,t]⋅fs−1⋅f′4n−t. | (4.2) |
By Eq (4.8), we know that the result of detLA[s,t] will change with the values of s and t. Then we can get the following twenty cases:
Case 1. i=4s, j=4t, 1≤s<t≤n,
detφ=|45−1√5−2√51−15−1545−2√5⋱−2545−25−25−45−1√5−151−1√5−25−45|(4t−4s−1)=10(t−s)(4125)t−s. |
Case 2. i=4s, j=4t+1, 1≤s≤t≤n−1,
detφ=|45−1√5−1√51−15−1545−2√5⋱−2545−15−151−1√5−1√545−2√5−25−45|(4t−4s)=(4t−4s+1)(4125)t−s. |
Case 3. i=4s, j=4t+2, 1≤s≤t≤n−1,
detφ=|45−1√5−1√51−15−1545−2√5⋱−1√51−1√5−1√545−25−2545−25−25−45|(4t−4s+1)=45(2t−2s+1)(4125)t−s. |
Case 4. i=4s, j=4t+3, 1≤s≤t≤n−1,
detφ=|45−1√5−1√51−15−1545−2√5⋱−1√545−25−2545−25−2545−1√5−1√51|(4t−4s+2)=15(4t−4s+3)(4125)t−s. |
Case 5. i≡0, j=4n, 1≤s≤t,
detφ=|45−1√5−2√51−15−1√545−25⋱−2545−25−2545−1√5−1√51−1√5−2√545|(4n−4s−1)=10(n−s)(4125)n−s. |
Case 6. i=4s+1, j=4t, 1≤s<t≤n,
detφ=|1−1√5−1√545−2√5−2√545−2√5⋱−2545−25−2545−1√5−1√51−1√5−1√545|(4t−4s−2)=254(4t−4s−1)(4125)t−s. |
Case 7. i=4s+1, j=4t+1, 1≤s<t≤n−1,
detφ=|1−1√5−1√545−2√5−2√545−2√5⋱−2545−1√5−1√51−1√5−1√545−25−2545|(4t−4s−1)=10(t−s)(4125)t−s. |
Case 8. i=4s+1, j=4t+2, 1≤s<t≤n−1,
detφ=|1−1√5−1√545−2√5−2√545−2√5⋱−251−1√5−1√545−25−2545−25−2545|(4t−4s)=(4t−4s+1)(4125)t−s. |
Case 9. i=4s+1, j=4t+3, 0≤s≤t≤n−1,
detφ=|1−1√5−1√545−2√5−2√545−2√5⋱−1√545−25−2545−25−2545−1√5−1√51|(4t−4s+1)=(2t−2s+1)(4125)t−s. |
Case 10. i≡1, j=4n+1, 0≤s≤n,
detφ=|1−1√5−1√545−2√5−2√545−2√5⋱−2545−25−2545−1√5−1√51−1√5−1√545|(4n−4s−2)=254(4n−4s−1)(4125)n−s. |
Case 11. i=4s+2, j=4t, 0≤s<t≤n,
detφ=|45−25−2545−25−2545−25⋱−2545−25−2545−1√5−1√51−1√5−1√545|(4t−4s−3)=25(2t−2s−1)(4125)t−s. |
Case 12. i=4s+2, j=4t+1, 0≤s<t≤n−1,
detφ=|45−25−2545−25−2545−25⋱−1√51−1√5−1√545−25−2545−25−2545|(4t−4s−2)=5(4t−4s−1)(4125)t−s. |
Case 13. i=4s+2, j=4t+2, 0≤s<t≤n−1,
detφ=|45−25−2545−25−2545−25⋱−1√51−1√5−1√545−25−2545−25−2545|(4t−4s−1)=8(t−s)(4125)t−s. |
Case 14. i=4s+2, j=4t+3, 0≤s≤t≤n−1,
detφ=|45−25−2545−25−2545−25⋱−1√545−25−2545−25−2545−1√5−1√51|(4t−4s)=(4t−4s+1)(4125)t−s. |
Case 15. i≡2, j=4t+3, 0≤s≤n−1,
detφ=|45−25−2545−25−2545−25⋱−2545−25−2545−1√5−1√51−1√5−1√545|(4n−4s−3)=25(2n−2s−1)(4125)n−s. |
Case 16. i=4s+3, j=4t, 0≤s<t≤n,
detφ=|45−25−2545−1√5−1√51−1√5⋱−2545−25−2545−1√5−1√51−1√5−1√545|(4t−4s−4)=1254(4t−4s−3)(4125)t−s. |
Case 17. i=4s+3, j=4t+1, 0≤s<t≤n−1,
detφ=|45−25−2545−1√5−1√51−1√5⋱−2545−1√5−1√51−1√5−1√545−1√5−1√545|(4t−4s−3)=25(2t−2s−1)(4125)t−s. |
Case 18. i=4s+3, j=4t+2, 0≤s<t≤n−1,
detφ=|45−25−2545−1√5−1√51−1√5⋱−1√51−1√5−1√545−25−2545−25−2545|(4t−4s−3)=253(4t−4s−1)(4125)t−s. |
Case 19. i=4s+3, j=4t+3, 0≤s<t≤n−1,
detφ=|45−25−2545−1√5−1√51−1√5⋱−1√545−25−2545−25−2545−1√5−1√51|(4t−4s−1)=10(t−s)(4125)t−s. |
Case 20. i≡3, j=4t, 0≤s≤n−1,
detφ=|45−25−2545−1√5−1√51−1√5⋱−2545−25−251−1√5−1√51−1√5−1√545|(4n−4s−4)=1254(4n−4s−3)(4125)n−s. |
Therefore, we can get
(−1)4n−2b4n−2=E1+E2+E3+E4, |
where
E1=118(227n3+347n2−574n+4)(4125)n−1. |
E2=172(908n3+3431n2+523n)(4125)n. |
E3=145(454n3+1375n2−1079n)(4125)n. |
E4=181(92n3+561n2−611n)(4125)n−1. |
Hence
(−1)4n−2b4n−2=E1+E2+E3+E4=13240(14520n3+4599n2−1496n+4)(4125)n. |
The proof of Theorem 4.2 completed.
Let 0≤ξ1≤ξ2≤ξ3≤⋯≤ξ3n+2 are the eigenvalues of LA. We can get the following exact equation:
4n∑i=21ξi=(−1)4n−2b4n−2(−1)4n−1b4n−1=172(14520n3+4599n2−1496n+838n−8). |
Theorem 4.3. Set L6,4,4n be the derivative [n]pheylenes, and the expression of the multiplicative degree-Kirchhoff index is
Kf∗(Ln)=(29040n3+8996n2−3198n+8144). |
Proof. Together with Eq.(4.7) and Theorems 4.1 and 4.2, one can get
Kf∗(Ln)=2(19n−4)(4n∑i=21ξi+4n∑i=11γi)=2(19n−4)[172(14520n3+4599n2−1496n+838n−8)+21n−16]=(29040n3+8996n2−3198n+8144). |
The result as desired.
The multiplicative degree-Kirchhoff indices of Ln from L1 to L12, see Table 2.
G | Kf∗(G) | G | Kf∗(G) | G | Kf∗(G) |
L1 | 241.98 | L5 | 26659.15 | L9 | 151875.4 |
L2 | 1818.86 | L6 | 45675.81 | L10 | 207691.9 |
L3 | 5940.68 | L7 | 72077.4 | L11 | 275733.2 |
L4 | 13817.44 | L8 | 107073.9 | L12 | 357209.6 |
Then we want to calculate the Gutman index of Ln.
Theorem 4.4. Suppose that L6,4,4n is the dicyclobutadieno derivative of [n]phenylenes and the graph Ln is obtained from the transformation of the graph L6,4,4n, then
limn→∞Kf∗(Ln)Gut(Ln)=14. |
Proof. Consider dij for all vertices. We divide the vertices of Ln into the following four categories.
Case 1. Vertex 4i−2(i=1,2,⋯,n) of Ln:
f4i−2=103n(56n2−24n+37). |
Case 2. Vertex 4i−1(i=1,2,⋯,n) of Ln:
f4i−1=103n(152n2−48n+29). |
Case 3. Vertex 4i(i=1,2,⋯,n) of Ln:
f4i=103n(140n2−48n+43). |
Case 4. Vertex 4i−3(i=1,2,⋯,n) of Ln:
f4i−3=103n(136n2−6n+71). |
According to Eq.(1.3), the Gutman index of Ln is
Gut(Ln)=f4i+f4i−1+f4i−2+f4i−32=10n3(242n2−63n+61). |
Therefore, combining with Kf∗(Ln) and Gut(Ln), we have
limn→∞Kf∗(Ln)Gut(Ln)=14. |
The result as desired.
Finally, we want to know the complexity of Ln.
Theorem 4.5. For the graph Ln, we have
τ(Ln)=23n+2⋅33n−2. |
Proof. Based on Lemma 2.4, we can get
8n∏i=1di4n∏i=2αi4n∏j=1βi=2(19n−4)⋅τ(Ln). |
Note that
8n∏i=1di=34⋅42n⋅56n−4. |
4n∏i=2αi=259⋅(38n−8)⋅(4125)n. |
4n∏j=1βi=(43)2⋅(65)3n−2. |
Hence,
τ(Ln)=23n+2⋅33n−2. |
The proof is over.
Thus, we can get the complexity of Ln from W1 to W8 which are listed in Table 3.
G | τ(G) | G | τ(G) |
W1 | 96 | W5 | 208971104256 |
W2 | 20736 | W6 | 45137758519296 |
W3 | 4478976 | W7 | 9749755840167936 |
W4 | 967458816 | W8 | 2105947261476274176 |
In this paper, the linear chain network with n hexagons and 2n−1 squares is considered. We have devoted ourselves to calculating the (multiplicative degree) Kirchhoff index, Wiener index, Gutman index, and complexity. In the meantime, we deduced that the ratio of the (multiplicative degree) Kirchhoff index to the (Gutman) Wiener index is nearly a quarter when n tends to infinity. These rules also apply to some other graphs.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This research was funded by the Anhui Provincial Natural Science Research Major Project (No. 2022AH040317) and the Anhui Provincial 2023 Action Project for Cultivating Young and Middle Aged Teachers in Universities (No. DTR2023095).
No potential conflicts of interest were reported by the authors.
[1] |
A. Narasimhan, S. Sadasivam, Non-Fourier bio heat transfer modelling of thermal damage during retinal laser irradiation, Int. J. Heat Mass Tran., 60 (2013), 591–597. https://doi.org/10.1016/j.ijheatmasstransfer.2013.01.010 doi: 10.1016/j.ijheatmasstransfer.2013.01.010
![]() |
[2] |
J. Zhou, Y. Zhang, J. Chen, Non-Fourier heat conduction effect on laser-induced thermal damage in biological tissues, Numer. Heat Tr. A-Appl., 54 (2008), 1–19. https://doi.org/10.1080/10407780802025911 doi: 10.1080/10407780802025911
![]() |
[3] |
T. Xue, X. Zhang, K. K. Tamma, Investigation of thermal interfacial problems involving non-locality in space and time, Int. Commun. Heat Mass, 99 (2018), 37–42. https://doi.org/10.1016/j.icheatmasstransfer.2018.10.008 doi: 10.1016/j.icheatmasstransfer.2018.10.008
![]() |
[4] |
Z. Y. Guo, Y. S. Xu, Non-Fourier heat conduction in IC chip, J. Electron. Packag. Sep., 117 (1995), 174–177. https://doi.org/10.1115/1.2792088 doi: 10.1115/1.2792088
![]() |
[5] |
H. D. Wang, B. Y. Cao, Z. Y. Guo, Non-Fourier heat conduction in carbon nanotubes, J. Heat Transfer, 134 (2012), 051004. https://doi.org/10.1115/1.4005634 doi: 10.1115/1.4005634
![]() |
[6] | C. Cattaneo, A form of heat-conduction equations which eliminates the paradox of instantaneous propagation, Comptes Rendus, 247 (1958), 431. |
[7] | P. Vernotte, Some possible complications in the phenomena of thermal conduction, Compte Rendus, 252 (1961), 2190–2191. |
[8] |
M. Ozisik, D. Y. Tzou, On the wave theory in heat conduction, J. Heat Transfer, 116 (1994), 526–535. https://doi.org/10.1115/1.2910903 doi: 10.1115/1.2910903
![]() |
[9] |
J. Ordonez-Miranda, J. Alvarado-Gil, Thermal wave oscillations and thermal relaxation time determination in a hyperbolic heat transport model, Int. J. Therm. Sci., 48 (2009), 2053–2062. https://doi.org/10.1016/j.ijthermalsci.2009.03.008 doi: 10.1016/j.ijthermalsci.2009.03.008
![]() |
[10] |
A. Vedavarz, S. Kumar, M. K. Moallemi, Significance of non-Fourier heat waves in conduction, J. Heat Transfer, 116 (1994), 221–224. https://doi.org/10.1115/1.2910859 doi: 10.1115/1.2910859
![]() |
[11] | D. Y. Tzou, J. C. Dowell, Computational techniques for microscale heat transfer, In: Handbook of Numerical Heat Transfer, 2000. https://doi.org/10.1002/9780470172599.ch20 |
[12] |
K. C. Liu, P. J. Cheng, Finite propagation of heat transfer in a multilayer tissue, J. Thermophys. Heat Tr., 22 (2008), 775–782. https://doi.org/10.2514/1.37267 doi: 10.2514/1.37267
![]() |
[13] |
C. Barman, P. Rath, A. Bhattacharya, A non-Fourier bioheat transfer model for cryosurgery of tumor tissue with minimum collateral damage, Comput. Meth. Prog. Bio., 200 (2021), 105857. https://doi.org/10.1016/j.cmpb.2020.105857 doi: 10.1016/j.cmpb.2020.105857
![]() |
[14] |
P. Namakshenas, A. Mojra, Microstructure-based non-Fourier heat transfer modeling of HIFU treatment for thyroid cancer, Comput. Meth. Prog. Bio., 197 (2020), 105698. https://doi.org/10.1016/j.cmpb.2020.105698 doi: 10.1016/j.cmpb.2020.105698
![]() |
[15] |
Z. S. Deng, J. Liu, Non-Fourier heat conduction effect on prediction of temperature transients and thermal stress in skin cryopreservation, J. Therm. Stresses, 26 (2003), 779–798. https://doi.org/10.1080/01495730390219377 doi: 10.1080/01495730390219377
![]() |
[16] |
A. K. Kheibari, M. Jafari, M. B. Nazari, Propagation of heat wave in composite cylinder using Cattaneo-Vernotte theory, Int. J. Heat Mass Transfer, 160 (2020), 120208. https://doi.org/10.1016/j.ijheatmasstransfer.2020.120208 doi: 10.1016/j.ijheatmasstransfer.2020.120208
![]() |
[17] | A. Fehér, R. Kovacs, Novel evaluation method for non-Fourier effects in heat pulse experiments, 2021. https://doi.org/10.48550/arXiv.2101.01123 |
[18] | A. Fehér, D. Markovics, T. Fodor, R. Kovacs, Size effects and non-Fourier thermal behaviour in rocks, ISRM EUROCK, 2020. |
[19] |
P. Duhamel, A new finite integral transform pair for hyperbolic conduction problems in heterogeneous media, Int. J. Heat Mass Transfer, 44 (2001), 3307-3320. https://doi.org/10.1016/S0017-9310(00)00360-4 doi: 10.1016/S0017-9310(00)00360-4
![]() |
[20] |
S. Singh, Z. Li, A high order compact scheme for a thermal wave model of bio-heat transfer with an interface, Numer. Math. Theory Me., 11 (2018), 321–337. https://doi.org/10.4208/nmtma.OA-2017-0048 doi: 10.4208/nmtma.OA-2017-0048
![]() |
[21] |
M. Asif, F. Bilal, R. Bilal, N. Haider, S. A. M. Abdelmohsenc, S. M. Eldind, An efficient algorithm for the numerical solution of telegraph interface model with discontinuous coefficients via Haar wavelets, Alex. Eng. J., 72 (2023), 275–285. https://doi.org/10.1016/j.aej.2023.03.074 doi: 10.1016/j.aej.2023.03.074
![]() |
[22] |
T. Ahmod, J. Dutta, Finite element method for hyperbolic heat conduction model with discontinuous coefficients in one dimension, Proc. Math. Sci., 132 (2022), 6. https://doi.org/10.1007/s12044-021-00646-3 doi: 10.1007/s12044-021-00646-3
![]() |
[23] |
B. Deka, J. Dutta, Finite element methods for non-Fourier thermal wave model of bio heat transfer with an interface, J. Appl. Math. Comput., 62 (2020), 701–724. https://doi.org/10.1007/s12190-019-01304-8 doi: 10.1007/s12190-019-01304-8
![]() |
[24] |
B. Deka, J. Dutta, Convergence of finite element methods for hyperbolic heat conduction model with an interface, Comput. Math. Appl., 79 (2020), 3139–3159. https://doi.org/10.1016/j.camwa.2020.01.013 doi: 10.1016/j.camwa.2020.01.013
![]() |
[25] |
J. Dutta, B. Deka, Optimal a priori error estimates for the finite element approximation of dual-phase-lag bio heat model in heterogeneous medium, J. Sci. Comput., 87 (2021), 58. https://doi.org/10.1007/s10915-021-01460-9 doi: 10.1007/s10915-021-01460-9
![]() |
[26] |
S. Han, Finite volume solution of 2-D hyperbolic conduction in a heterogeneous medium, Numer. Heat Tr. A-Appl., 70 (2016), 723–737. https://doi.org/10.1080/10407782.2016.1193347 doi: 10.1080/10407782.2016.1193347
![]() |
[27] |
H. Sauerland, T. P. Fries, The stable XFEM for two-phase flows, Comput. Fluids, 87 (2013), 41–49. https://doi.org/10.1016/j.compfluid.2012.10.017 doi: 10.1016/j.compfluid.2012.10.017
![]() |
[28] |
B. Ayuso de Dios, M. Holst, Y. Zhu, L. Zikatanov, Multilevel preconditioners for discontinuous Galerkin approximations of elliptic problems with jump coefficients, Math. Comp., 83 (2014), 1083–1120. https://doi.org/10.1090/S0025-5718-2013-02760-3 doi: 10.1090/S0025-5718-2013-02760-3
![]() |
[29] |
M. J. Gander, Optimized Schwarz methods, SIAM J. Numer. Anal., 44 (2006), 699–731. https://doi.org/10.1137/S0036142903425409 doi: 10.1137/S0036142903425409
![]() |
[30] |
Y. Maday, F. Magoules, Optimized Schwarz methods without overlap for highly heterogeneous media, Comput. Method. Appl. M., 196 (2007), 1541–1553. https://doi.org/10.1016/j.cma.2005.05.059 doi: 10.1016/j.cma.2005.05.059
![]() |
[31] |
M. J. Gander, O. Dnbois, Optimized Schwarz methods for a diffusion problem with discontinuous coefficient, Numer. Algor., 69 (2015), 109–144. https://doi.org/10.1007/s11075-014-9884-2 doi: 10.1007/s11075-014-9884-2
![]() |
[32] |
M. J. Gander, T. Vanzan, Heterogeneous optimized Schwarz methods for second order elliptic PDEs, SIAM J. Sci. Comput., 41 (2019), A2329–A2354. https://doi.org/10.1137/18M122114X doi: 10.1137/18M122114X
![]() |
[33] |
M. J. Gander, S. B. Lunowa, Ch. Rohde, Non-overlapping Schwarz Waveform-Relaxation for Nonlinear Advection-Diffusion Equations, SIAM J. Sci. Comput., 45 (2023), A49–A73. https://doi.org/10.1137/21M1415005 doi: 10.1137/21M1415005
![]() |
[34] | M. J. Gander, L. Halpern, M. Kern, A Schwarz waveform relaxation method for advection-diffusion-reaction problems with discontinuous coefficients and non-matching grids, In: Domain Decomposition Methods in Science and Engineering, Berlin: Springer, 2007,283–290. https://doi.org/10.1007/978-3-540-34469-8_33 |
[35] | L. Halpern, C. Japhet, Discontinuous Galerkin and nonconforming in time optimized Schwarz waveform relaxation for heterogeneous problems, In: Domain Decomposition Methods in Science and Engineering XVII, Berlin: Springer, 2008,211–219. https://doi.org/10.1007/978-3-540-75199-1_23 |
[36] |
L. Halpern, C. Japhet, J. Szeftel, Optimized Schwarz waveform relaxation and discontinuous Galerkin time stepping for heterogeneous problems, SIAM J. Numer. Anal., 50 (2012), 2588–2611. https://doi.org/10.1137/120865033 doi: 10.1137/120865033
![]() |
[37] |
M. D. Al-Khaleel, M. J. Gander, and P. M. Kumbhar, Optimized Schwarz waveform relaxation methods for the telegrapher equation, SIAM J. Sci. Comput., 46 (2024), A3528–A3551. https://doi.org/10.1137/24M1642962 doi: 10.1137/24M1642962
![]() |
[38] | M. J. Gander, J. Lin, S. L. Wu, X. Yue, T. Zhou, Parareal: parallel-in-time algorithms based on the diagonalization technique, 2020. https://doi.org/10.48550/arXiv.2005.09158 |
[39] |
M. J. Gander, S. Vandewalle, Analysis of the Parareal time-parallel time integration method, SIAM J. Sci. Comput., 29 (2007), 556–578. https://doi.org/10.1137/05064607X doi: 10.1137/05064607X
![]() |
[40] | G. Ciaramella, M. J. Gander, I. Mazzieri, Unmapped tent pitching schemes by waveform relaxation, In: Decomposition Methods in Science and Engineering XXVII, Cham: Springer, 2024,455–462. https://doi.org/10.1007/978-3-031-50769-4_54 |
[41] |
M. J. Gander, L. Halpern, Optimized Schwarz waveform relaxation methods for advection-reaction-diffusion problems, SIAM J. Numer. Anal., 45 (2007), 666–697. https://doi.org/10.1137/050642137 doi: 10.1137/050642137
![]() |
[42] |
D. Bennequin, M. J. Gander, L. Gouarin, L. Halpern, Optimized Schwarz waveform relaxation for advection-reaction-diffusion equations in two dimensions, Numer. Math., 134 (2016), 513–567. https://doi.org/10.1007/s00211-015-0784-8 doi: 10.1007/s00211-015-0784-8
![]() |
[43] |
V. Martin, An optimized Schwarz waveform relaxation method for the unsteady convection-diffusion equation in two dimensions, Appl. Numer. Math., 52 (2005), 401–428. https://doi.org/10.1016/j.apnum.2004.08.022 doi: 10.1016/j.apnum.2004.08.022
![]() |
[44] |
Y. Xu, The influence of domain truncation on the performance of optimized Schwarz methods, Electron. T. Numer. Ana., 49 (2018), 182–209. https://doi.org/10.1553/etna_vol49s182 doi: 10.1553/etna_vol49s182
![]() |
[45] |
W. B. Dong, H. S. Tang, Convergence analysis of Schwarz waveform relaxation method to compute coupled advection-diffusion-reaction equations, Math. Comput. Simul., 218 (2024), 462–481. https://doi.org/10.1016/j.matcom.2023.11.026 doi: 10.1016/j.matcom.2023.11.026
![]() |
[46] |
M. J. Gander, V. Martin, Why Fourier mode analysis in time is different when studying Schwarz waveform relaxation, J. Comput. Phys., 491 (2023), 112316. https://doi.org/10.1016/j.jcp.2023.112316 doi: 10.1016/j.jcp.2023.112316
![]() |
G | Kf(G) | G | Kf(G) | G | Kf(G) |
L1 | 17.3 | L5 | 1486.7 | L9 | 8268.0 |
L2 | 110.7 | L6 | 2524.0 | L10 | 24457.3 |
L3 | 344.0 | L7 | 3957.3 | L11 | 30454.7 |
L4 | 781.3 | L8 | 5850.7 | L12 | 37360.0 |
G | Kf∗(G) | G | Kf∗(G) | G | Kf∗(G) |
L1 | 241.98 | L5 | 26659.15 | L9 | 151875.4 |
L2 | 1818.86 | L6 | 45675.81 | L10 | 207691.9 |
L3 | 5940.68 | L7 | 72077.4 | L11 | 275733.2 |
L4 | 13817.44 | L8 | 107073.9 | L12 | 357209.6 |
G | τ(G) | G | τ(G) |
W1 | 96 | W5 | 208971104256 |
W2 | 20736 | W6 | 45137758519296 |
W3 | 4478976 | W7 | 9749755840167936 |
W4 | 967458816 | W8 | 2105947261476274176 |
G | Kf(G) | G | Kf(G) | G | Kf(G) |
L1 | 17.3 | L5 | 1486.7 | L9 | 8268.0 |
L2 | 110.7 | L6 | 2524.0 | L10 | 24457.3 |
L3 | 344.0 | L7 | 3957.3 | L11 | 30454.7 |
L4 | 781.3 | L8 | 5850.7 | L12 | 37360.0 |
G | Kf∗(G) | G | Kf∗(G) | G | Kf∗(G) |
L1 | 241.98 | L5 | 26659.15 | L9 | 151875.4 |
L2 | 1818.86 | L6 | 45675.81 | L10 | 207691.9 |
L3 | 5940.68 | L7 | 72077.4 | L11 | 275733.2 |
L4 | 13817.44 | L8 | 107073.9 | L12 | 357209.6 |
G | τ(G) | G | τ(G) |
W1 | 96 | W5 | 208971104256 |
W2 | 20736 | W6 | 45137758519296 |
W3 | 4478976 | W7 | 9749755840167936 |
W4 | 967458816 | W8 | 2105947261476274176 |