Efficient Routing Protection Algorithm Based on Optimized Network Topology

Haijun Geng,Zikun Jin,Jiangyuan Yao,Han Zhang,Zhiguo Hu,Bo Yang,Yingije Guo,Wei Wang,Qidong Zhang and Guoao Duan

1School of Computer and Information Technology,Shanxi University,Shanxi,030006,China

2School of Automation and Software Engineering,Shanxi University,Shanxi,030006,China

3School of Computer Science and Technology,Hainan University,Haikou,570228,China

4Department of Computer Science&Technology,Tsinghua University,Beijing,100084,China

5Department of Urban&Regional Planning San José State University,San José,CA,95192,USA

6Key Laboratory of Embedded System and Service Computing(Tongji University),Ministry of Education,Shanghai,201804,China

7Institute of Fundamental and Frontier Sciences,University of Electronic Science and Technology of China,Chengdu,610054,China

8Beijinng University of Technology,Beijing-Dublin International College,Beijing University of Technology,Beijing,100124,China

Abstract: Network failures are unavoidable and occur frequently.When the network fails,intra-domain routing protocols deploying on the Internet need to undergo a long convergence process.During this period, a large number of messages are discarded, which results in a decline in the user experience and severely affects the quality of service of Internet Service Providers(ISP).Therefore, improving the availability of intra-domain routing is a trending research question to be solved.Industry usually employs routing protection algorithms to improve intra-domain routing availability.However, existing routing protection schemes compute as many backup paths as possible to reduce message loss due to network failures, which increases the cost of the network and impedes the methods deployed in practice.To address the issues,this study proposes an efficient routing protection algorithm based on optimized network topology(ERPBONT).ERPBONT adopts the optimized network topology to calculate a backup path with the minimum path coincidence degree with the shortest path for all source purposes.Firstly,the backup path with the minimum path coincidence with the shortest path is described as an integer programming problem.Then the simulated annealing algorithm ERPBONT is used to find the optimal solution.Finally,the algorithm is tested on the simulated topology and the real topology.The experimental results show that ERPBONT effectively reduces the path coincidence between the shortest path and the backup path, and significantly improves the routing availability.

Keywords: Routing protection algorithm; network availability; network failure;shortest path;backup path

1 Introduction

With the development of the Internet,the range of supported applications has shown significant changes.Initially,the Internet mainly supported some non-real-time applications,such as e-mail,file transfer, etc.Nowadays, a large amount of real-time business data [1,2] is widely spread over the Internet, such as Voice over Internet Protocol (VoIP), online stock trading, remote surgery, video streaming,instant messaging,and so forth.These new applications put forward higher requirements for routing availability [3,4].Thus, routing availability will directly affect the property security and even users’life quality.Relevant studies have shown that network failures frequently occur[5,6],and are unavoidable.During the failure repair period,the routing protocol needs to go through a period of the convergence process.In this process,there will be routing black holes,routing loops[7-9],which lead to network interruption,packet loss,and greatly reduce routing availability.Previous studies[10]demonstrated that when an OC-48 link is disconnected for 10 s, 3 million 1KB packets will be lost.In real-time applications,it is usually required to repair the failures in milliseconds.For current intradomain routing protocols,e.g.,OSPF and IS-IS,it takes several or even tens of seconds to complete convergence,which fails to meet the requirements of real-time applications for routing availability[11].Therefore,how to improve the availability of intra-domain routing[12]has become an urgent Internet research problem.More research are focused on the problem of how to improve the network’s ability to respond to failures quickly.

Current algorithms aiming to improve the availability of intra-domain routing can be divided into two categories[13],passive recovery algorithms and routing protection algorithms.Among them,the passive recovery algorithm[14]mainly investigates how to speed up the route convergence when the network fails,and shorten the network interruption time as much as possible.However,when the links in a network are frequently disconnected,the algorithm may cause routing instability.According to the relevant rules,the routing protection algorithm[15]calculates the backup routes in advance.When the network fails,the backup routes calculated in advance are used to forward packets,which minimizes the loss of messages and shortens the network service interruption time.In comparison,the routing protection algorithm is more favored by the academic community.

As the earliest routing protection algorithm applied to the Internet, Equal Cost Multiple Paths(ECMP)proposed that if there are multiple paths from the source node to the destination node with the same minimum cost, they can be used as a backup path.Although the algorithm is simple to implement and easy to deploy,it requires the backup path to have the same minimum cost.Therefore,the ECMP algorithm does not contribute much to routing availability.Research [16] was developed to configure the shortest path with the same cost by adjusting the link weights,but this problem was proved to be an NP-Hard problem.

In response to the problems of ECMP, the Internet Engineering Task Force (IETF)released a framework for fast rerouting.On the basis of this framework,Loop-free alternates(LFA),Not-Viabased routing protection scheme[17]and tunnel-based routing protection scheme[18]are proposed.In all routing protection schemes,LFA has received close attention from the industry for its simplicity,and has been deployed and supported by router manufacturers such as Huawei and H3C.Although LFA is simple and easy to deploy,it has a fatal disadvantage—LFA cannot protect all possible single failure scenarios in the network.In order to solve the problems, the author used the theoretical knowledge of graph theory to analyze the problem of LFA failure coverage in the paper [19].The routing availability of LFA is increased by adjusting the cost of the link in the network,but this method does not necessarily guarantee to deal with all single failure situations.For this reason, the author theoretically analyzed in-depth relationship between LFA’s routing availability and network topology[20],and improved the single-failure protection situation of LFA by adding links.

Reference [21] was proposed to use the routing deflection algorithm to improve the routing availability,that is,first use the loop-free rule to calculate all the optional next hops from the source node to the destination node,and then use the label technology to realize the flexible forwarding of the message.Although this algorithm can improve routing availability,its implementation is complicated and expensive,making it difficult to actually deploy.

Multiple Routing Configurations (MRC)[22] proposes to save multiple configuration diagrams for each router.Each configuration includes all nodes and links.By adjusting link weights,some nodes and links are protected in this configuration.Finally,a configuration diagram for all possible single failures is constructed.When the message encounters a failure in the forwarding process, the preconfigured route can be used to forward the message.However,in real networks,this algorithm needs to consume a lot of computing resources and storage overhead.

Failure Carrying Packet(FCP)[23]proposes to store the failure information encountered in the forwarding process of the message in the header of the message.When the message reaches a node,the node first detects the failure information field in the header of the message,constructs a new topology according to the field,and then re-calculates the shortest path using the new topology,so as to realize the loop-free forwarding of messages.The major advantage of this algorithm is that it eliminates the routing convergence process.However,it has high computational complexity and large changes to the routing protocol,so it is not easy to deploy in practice.

In order to solve the problem that the existing routing protection methods have low failure protection ration and can not be deployed in practice, this paper proposes an Efficient Routing Protection Algorithm Based on Optimized Network Topology(ERPBONT).

The contributions of the paper can be summarized as follows:

1.ERPBONT is easy to implement and deploy in the real networks.

2.ERPBONT can effectively solve the problem of low failure recovery ratio caused by network failures.

3.ERPBONT can significantly reduce the path coincidence degree,and has the ability to quickly recover from the failed network.

2 Network Model and Problem Description

2.1 Network Model

The network topology can be represented by an undirected connected graph [24-26]G=(V,E,W).WhereV,E,Wrespectively represents the set of nodes, the set of links, and the set of costs of links in the network,is the complement ofE.For link(u,v) ∈Ein the network,w(u,v)represents the cost of the link.Currently,the intra-domain routing protocols deployed on the Internet,such as OSPF and IS-IS, takeG=(V,E,W) as input and use the Dijkstra algorithm to calculate the shortest path of all source-destination pairs.When the network is in a normal state, all sourcedestination pairs use the shortest path to transmit packets,but when a failure occurs in the network,the packets passing through the failure will be discarded.Therefore, the routing protection method is widely used in the industry to minimize the message loss caused by network failure.The routing protection method is to calculate some backup paths for all source-destination pairs in the network.When the shortest path fails,the backup path can be used to forward packets,thereby greatly reducing the packet loss ratio caused by network failures.However, the existing routing protection algorithm calculates multiple backup paths for all source-destination pairs in order to increase the probability of the message being successfully forwarded,which not only increases the calculation overhead,but is also not easy to be deployed in actual networks.In order to solve the above-mentioned problems of existing routing protection algorithms, this paper calculates only one backup path for all sourcedestination pairs,and uses path coincidence to measure the difference between the backup path and the shortest path.

The following details the path coincidence degree.Given a network topologyG=(V,E,W),for any source-destination pairsanddin the network,SP(s,d)represents the set of edges in the shortest path between the pair of nodessanddin network topologyG=(V,E,W).The optimized network topology isG′=(V,E′,W′),BP(s,d)represents the set of edges in the shortest path between node pairssanddin network topologyG′=(V,E′,W′).

The path coincidence degree between the pair of nodes can be expressed as:

The path coincidence degree between all node pairs in the network can be expressed as:

The network model defined above is explained by a specific example.The figure on the left in Fig.1 shows a network topology with four nodes and four links,E= {(a,b),(a,d),(b,c),(c,d)},E= {(a,c),(b,d)},w(a,b) =w(b,a) = 3.The shortest pathSP(a,c) =(a,b,c)between node pairsa-cand the shortest pathSP(b,d) =(b,c,d)between node pairsb-d.The figure on the right in Fig.1 shows the optimized network topology,the backup path between node paira-cisBP(a,c)=(a,c),CD(a,c) == 0.Similarly, the backup path between node pairsb-disBP(b,d) =(b,d),CD(b,d)==0.

Figure 1:Original network topology(left)and optimized network topology(right)

2.2 Problem Description

This paper aims to calculate two paths with the minimum coincidence degree for all sourcedestination pair in the network, which are the shortest path and backup path.Given a network topologyG=(V,E,W), the Dijkstra algorithm can be used to calculate the shortest path of all source-destination pairs.How to calculate the backup path is the key issue that this article needs to solve.

Given a network topologyG=(V,E,W), there are many ways to calculate the backup path between all pairs of nodes.In order to be easy to deploy in the network,we choose the straightforward calculation method.That is, by expanding the original network topology structure, a new network topology structure is obtained, and the shortest path of the source-destination pair calculated in the new network topology structure is the backup path between the node pair.The problem can be expressed as:

Given a network topologyG=(V,E,W), how to construct an extended topologyG′=(V,E′,W′) of the topologyG=(V,E,W), and use the extended topologyG′=(V,E′,W′) to calculate the backup path,so as to minimize the overlap of the network path.According to the original network topology, the method of constructing the extended network topology mainly includes: (1)Adjust the link weight(2)Virtually delete the link(3)Increase the link(4)Increase the nodes in the network.Method(1)and method(2)are relatively simple.Method(2)does not actual delete links in the network,but saves these links when calculating the backup path.Method(3)needs to add a certain number of links in the network.If the increased number of links is controlled within a certain range,this method will not introduce large additional overhead.Method(4)requires an additional increase in the number of nodes,and the deployment overhead is relatively large.In summary,this article adopts fusion methods(2)and(3)to construct the extended network topology.Mathematically,

Eq.(1)is the objective function to minimize the path coincidence between the shortest path and the backup path.Eq.(2)indicates that the number of added links is the minimum value between 0.1 times and 5 of the number of links in the original network,whereALis the set of added links.Eq.(3)shows that the cost of adding links satisfies the triangle rule,that is,the cost of links in the network is the minimum of all path costs.Eq.(4)indicates that the number of deleted links is the minimum value between 0.1 times and 5 of the number of links in the original network,whereDLis the set of deleted links.Eq.(5)indicates that the number of links in the extended network topology is equal to the number of links in the original network topology plus the increased number of links,and then minus the number of deleted links.Eqs.(6)and(7)show that the cost of links in the network is symmetrical.Eqs.(8)and(9)show that the link cost of the network used in this paper is symmetrical.

3 Algorithm

The scientific problem solved in this paper is an NP-Hard problem, and it is difficult to obtain the optimal solution in limited computing time.Therefore, in order to reduce the deployment cost,this paper uses the heuristic method ERPBONT to calculate the approximate optimal solution.The algorithm are described in detail below.

3.1 ERPBONT Algorithm

ERPBONT is a two-stage algorithm,which is composed of adding link algorithm AddLink and Deleting Link algorithm DeleteLink.AddLink and DeleteLink expand the original network topology by adding and deleting links in the network,respectively.Through AddLink and DeleteLink,a better extended network topology can be calculated,thereby reducing the coincidence of network paths.

Algorithm 1 introduces the execution process of ERPBONT in detail.The input of the algorithm is the original network topology, the shortest path set.The output of the algorithm is the extended network topology and the backup path.First use the AddLink algorithm to calculate the final added link setAL?L, and then use the DeleteLink algorithm to calculate the deleted link setDL?L.Finally,the extended network topologyG′=(V,E′,W′)is obtained,andE′=E-DL+AL.According to the extended network topology,the backup pathBPof all source-destination pairs in the network is calculated.

Algorithm 1: ERPBONT Input: G =(V,E,W),T0 Output: BP 1: Begin 2: AL ←AddLink(G,T0)3: DL ←DeleteLink(G ∪AL,T0)4: G′=(V,E+AL-DL,W′)5: BP=Dijkstra(G′)6: End

3.2 Add Link Algorithm

Algorithm 2:AddLink Input:G =(V,E,W),T0 Output : AL 1: T =T0 2: G′=G 3: AL ←?4: original =current ←C(G,G′)5: While current ≥0 and T ＞0 and|AL|＜5 6:(u,v)←argmin E(C(G,G′))7:E′←E ∪(u,v)8:G′=(V,E′,W′)9:(m,n)∈E-(u,v)10:current ←C(G,G′)11:If original ＜current or T ＞random(T0)Then 12:current ←C(G,G′)E ←(Continued)

Algorithm 2:Continued 13:AL ←AL ∪(u,v)14:else 15:E ∪(u,v)16:EndIf 17:T =T -1 18: EndWhile 19: Return AL E ←

The AddLink algorithm needs to solve the following three problems:(1)How many links to be added(2)Which links to be added(3)How to set the cost of adding links.The solutions to the three problems will be described below.

As for the problem (1), adding too many links will change the original network topology and increase the cost of service providers.On the contrary, adding too few links will not significantly improve the performance of the network.This paper uses a compromise method to solve this problem,increasing the number of links to |AL| = min(5,|E|*0.1) increase the number of links will not exceed five.

As for the problem(2),the links added by the topologyG=(V,E,W)can be selected from the set.The complexity of calculating the added link is,and the computational complexity is high,so this paper uses a simulated annealing algorithm to select the added link.

As for the problem (3), if the link cost is set too large, the shortest path between all nodes in the network will not pass through the link,and the value of increasing the link will be lost.For any additional link(u,v) ∈, the cost of the link isw(u,v) =cost(u,v)-min(W), Wherecost(u,v)represents the cost between the shortest path between nodeuand nodev, min(W) represents the minimum cost of all links in the network.

The input of the AddLink algorithm is the original network topology,the shortest path set and the initial temperature,where the initial temperature indicates the number of iterations of the algorithm,and the output is the increased link set.The variableTis set toT0,the extended network topology is initialized toG,the setALof added links is initialized to an empty set,and the original records the path coincidence degree between all pairs of nodes in the network(lines 1-4).The algorithm requires an iterative process to calculate the set of increased links.When the current path coincidence degree is equal to 0 or the system temperature is equal to 0 or the number of increase links is 5,the program execution is completed.In each iteration process,select a link(u,v)from the setto join the network,at this timeC(G,G′)has the minimum value.Update the parameters of the extended network topology(lines 7-8),delete the link(u,v)from the set,and calculate the current network coincidence degree(line 10).Whenoriginal ＜currentor the temperature in the system at this time is greater than the temperature generated by the random function, the link(u,v) is added to the setAL(lines 12-14);Otherwise,the operation of adding links(u,v)will be canceled(line 15 of the algorithm),and finally return to the set of adding linksAL(line 19).

Algorithm 3:DeleteLink Input: G =(V,E,W),AL,T0 Output : DL 1: T ←T0 2: G′←G U AL 3: DL ←?8: L ←E 5: original =current ←C(G,G′)6 : While current ≥0 and T ＞0 7:(u,v)←arg min(m,n)∈E (C(G,G′))8 :E′←L-(u,v)9 :G′=(V,E′)10 :If isconnected(G′)Then 11:current ←C(G,G′)12:else 13:E′←E′∪(u,v)14:continue 15:EndIf 16:If original ＜current or T ＞random(T0)Then 17:original ←C(G,G′)18:L ←L-(u,v)19:DL ←DL ∪(u,v)20:else 21:E′←E′∪(u,v)22:EndIf 23:T =T -1 24: EndWhile 25: Return AL

3.3 DeleteLink Algorithm

The Deletelink algorithm needs to solve the following two problems: (1)how many links are deleted and(2)which links are deleted.The solutions to the two problems will be described below.

For problem(1),the number of links added by the DeleteLink algorithm is|DL|=min(5,|E|*0.1).In order to keep the number of edges in the network unchanged,the number of links deleted by the DeleteLink algorithm is|DL|.

For problem(2),the deleted link can be selected from setDL?E.The complexity of calculating the deleted link is,which is high.Therefore,this paper uses a simulated annealing algorithm to select the deleted link.

The input of the DeleteLink algorithm is the original network topology, the outputALof AddLink and the initial temperature,where the initial temperature indicates the number of iterations of the algorithm,and the output is the deleted link set.Set the variable t toT0,initialize the extended network topology to the union ofGandAL, initialize the setDLof deleted links to an empty set,initialize the setLtoE, and original records the path coincidence degree between all node pairs in the network(lines 1-5).The algorithm needs to calculate the set of deleted links through an iterative process.When the current path coincidence degree is equal to 0 or the system temperature is equal to 0 or the number of deleted links is 5,the program execution is completed.In each iteration process,a link(u,v) is selected from the setLto join the network, andC(G,G′) has the minimum value at this time.If the networkG′is a connected graph at this time,update the parameters of the extended network topology,delete the link(u,v)from the setE,and calculate the current network coincidence degree(lines 7-11).If the networkG′is not a connected graph,the previous operation of updatingG′is canceled,and the next iteration process is executed.Whenoriginal ＜currentor the temperature in the system is greater than the temperature generated by the random function at this time,the link(u,v)is added to the setDL(lines 16-19);Otherwise,the operation of adding links(u,v)will be canceled(line 21),and finally return to the set of adding linksAL(line 25).

3.4 Algorithm Example

The following will use the form of a table to show the entire process of ERPBONT algorithm calculation on the topology of Fig.2.The left picture of Fig.2 represents the original topology,and the right picture is the topology optimized by the ERPBONT algorithm.The yellow edge represents the newly added link, and the red edge represents the virtual deleted link.The table only shows the change in the shortest path of the update link.A blank item indicates that the shortest path has not changed in the current state and is the same as the shortest path shown in the second column.The original topology generates links(a,c)and(b,d)under the action of the AddLink algorithm.At this time,the shortest path of the node toa-c,b-d,c-a,d-bchanges,as shown in the third column.It is continuing to calculate the previous topology by the DeleteLink algorithm,the shortest path of the node toa-b,b-a,c-d,d-chas also changed,as shown in the fourth column.According to the path coincidence degree formula,the path coincidence degree of the backup path calculated by the optimized topology is 25%.It can be concluded from Tab.1 that the public link of the backup path obtained by the ERPBONT algorithm is greatly reduced,and the path coincidence degree decreases significantly.

Figure 2:Algorithm example topology

Table 1: The calculation steps and results of the ERPBONT algorithm on the topology of Fig.2

Table 1:Continued

3.5 Algorithm Complexity Analysis

Theorem 1:The time complexity of the AddLink algorithm is:

Proof:It takes min(T0,cd1)cycles to find the link set that needs to be added.In each cycle,it is necessary to find the structureG′through the for loop,so as to calculate the shortest path and path coincidence degree of all nodes ofG′.At the end of the for loop,the shortest path and path coincidence degree calculation is performed.The time complexity of this algorithm isTherefore,the time complexity of AddLink is min

Theorem 2:The time complexity of the DeleteLink algorithm is:

Proof:The same as the AddLink proof process.

Theorem 3:The time complexity of the ERPBONT algorithm is:

Proof: ERPBONT mainly includes five steps: (1)To construct an undirected connected graph,the time complexity is(2)AddLink and DeleteLink,the time complexity is 2×min(T0,cd1)×(3)The shortest path algorithm and the path coincidence degree algorithm, the time complexity isTherefore, the time complexity of ERPBONT is: 2×min

4 Experiment Analysis

This section uses experiments to evaluate the path coincidence degree of ERPBONT, AddLink and DeleteLink.The following first introduces the network topology used in the experiment,and then describes the performance of different algorithms.

4.1 Network Topology

This paper uses three types of network topologies for experiments, the real network topology,the network topology published by Rocketfuel and the network topology generated by Brite.The real network topology includes Abilene,NJLATA,TORONTO and USLD[27].The topologies published by Rocketfuel [27,28] are 1755 (Ebone), 3257 (Tiscali), and 3967 (Exodus).When using Brite to generate the topology, it is assumed that the link weights are symmetrical.Brite’s model is set to Waxman,the number of nodes is 50 to 1000,the values of alpha and beta are 0.15 and 0.2,respectively,and the average degree of nodes is set between 2 and 10.The mode is set to the router,the node position obeys heavy-tail distribution,the size of the link bandwidth is between 10 and 1024,and the link cost is the inverse of the bandwidth.

The real topology used in the experiment is shown in Tab.2:

Table 2: Real topology

The Rocketfuel topology used in the experiment is shown in Tab.3:

Table 3: Rocketfuel topology

The Brite topology used in the experiment is shown in Tab.4:

Table 4: Brite topology

Table 4:Continued

4.2 Path Coincidence Degree

Figs.3 and 4 show the path coincidence degree corresponding to different algorithms in the real network topology and the measured network topology.As shown in Figs.3 and 4,in the real topology and measurement network topology,the ERPBONT algorithm’s performance is considerably better than the AddLink algorithm and the DeleteLink algorithm.Compared with the two algorithms,the path coincidence degree of the ERPBONT algorithm is reduced by 23.6%and 6.9%,respectively,and on the Rocketfule topology, the path coincidence degree of the ERPBONT algorithm is reduced by 17.8%and 3.7%,respectively.

Figure 3:Comparison of path coincidence degree of three algorithms in real topology

Fig.5 shows how the path coincidence degree of different algorithms changes with the size of the network topology when the average degree of nodes is 4 in the network.Fig.6 shows the relationship between the path coincidence degree of different algorithms and the average degree of network nodes when the network topology size is 200.It can be seen from Fig.5 that ERPBONT reduces the path overlap by 18.9%and 5.8%,respectively,compared with the AddLink algorithm and the DeleteLink algorithm.It can be seen from Fig.6 that the path coincidence degree of ERPBONT is 9.3%and 2.1%lower than that of the AddLink algorithm and the DeleteLink algorithm,respectively.

Figure 4:Comparison of path coincidence degree of three algorithms in Rocketfuel topology

Figure 5:Comparison of path coincidence degree of three algorithms in Brite degree 4 topology

Figure 6:Comparison of path coincidence degree of three algorithms in Brite node 200 topology

5 Conclusion and Future Work

The redundancy problem of backup path calculated by the existing routing protection methods is serious.The use of path coincidence degree as a measurement index is not considered in the selection of backup path, resulting in a large number of duplicate links between backup path and best path,which greatly reduces the routing availability.In order to solve the above problem,this paper proposes a routing protection method ERPBONT based on optimizing network topology.ERPBONT can be employed in the following scenarios, such as Voice over Internet Protocol (VoIP), online stock trading, remote surgery, video streaming, instant messaging, and so forth.ERPBONT takes the path coincidence degree as the index to select the backup path, calculates only one backup path for all source-destination pairs.ERPBONT has the minimum path coincidence degree between the backup path and the original path,which greatly improves the routing availability.According to the simulation experiment,it can be seen that the path coincidence degree of the backup path calculated by ERPBONT algorithm under the real topology and the simulated topology is significantly reduced.

At present,ERPBONT algorithm only stays in the stage of calculating the backup path to improve the routing availability.In the future, we will test the algorithm in a real network, and apply the algorithm ERPBONT to SDN network.We also consider how to improve ERPBONT to deal with the single node failures and concurrent failures scenarios in the network.

Funding Statement:This work is supported by the Hainan Provincial Natural Science Foundation of China (620RC562), the Natural Science Foundation of Shanxi Province (Grant Nos.20210302123444,20210302123455), the China University industry university research innovation fund(No.2021FNA02009),the Open Project Program of the Key Laboratory of Embedded System and Service Computing of Ministry of Education (Tongji University)ESSCKF 2021-04, the National Natural Science Foundation of China (Grant Nos.61702315, 61802092), the Applied Basic Research Plan of Shanxi Province(No.201901D211168),the Program of Hainan Association for Science and Technology Plans to Youth R&D Innovation (QCXM201910), and the Key R&D Program(International Science and Technology Cooperation Project)of Shanxi Province China(No.201903D421003).

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the study.