Abstract

Evaluating the reliability of a network with multiple sources to multiple sinks is a critical issue from the perspective of quality management. Due to the unrealistic definition of paths of network models in previous literature, existing models are not appropriate for real-world computer networks such as the Taiwan Advanced Research and Education Network (TWAREN). This paper proposes a modified stochastic-flow network model to evaluate the network reliability of a practical computer network with multiple sources where data is transmitted through several light paths (LPs). Network reliability is defined as being the probability of delivering a specified amount of data from the sources to the sink. It is taken as a performance index to measure the service level of TWAREN. This paper studies the network reliability of the international portion of TWAREN from two sources (Taipei and Hsinchu) to one sink (New York) that goes through a submarine and land surface cable between Taiwan and the United States.

1. Introduction

The issue of the QoS [1] of networks has been studied in the past decades. QoS is an important element of understanding the efficiency of real-world computer networks. It refers to the ability to provide a predictable, consistent data transfer service and the ability to satisfy customersโ€™ application needs while maximizing the use of network resources, especially a network reliability analysis. One of the traditional issues in this area of network reliability research is known as the source-sink (๐‘ -๐‘ก) network reliability problem [2โ€“16], which some articles refer to as two-terminal network reliability (TTNR) [14, 15]. In TTNR analyses, it is interesting to compute the network reliability in relation to the connecting paths between two specific network nodes, usually the source-sink (๐‘ -๐‘ก). Generally speaking, people are interested in obtaining the probability that the source connects the sink. Some researchers extend the study of TTNR to the ๐‘˜-terminal network reliability (KTNR) problem [17, 18], which contains at least one path from the source node to other ๐‘˜ nodes. Besides TTNR and KTNR, there is an all-terminal network reliability (ATNR) (also called overall or uniform network reliability), which is calculated from the probability that each and every node in the network is connected to each other [19, 20]. In a binary-state flow network, the capacity of each arc has two levels 0 and a positive integer. System with various states is called a stochastic model [21โ€“23]. For a more realistic system, the arc should have several possible states/capacities, and such a network is named a stochastic-flow network (or multistate network). The previous problems, TTNR, KTNR, and ATNR, are discussed for binary-state flow networks. However, this paper addresses the evaluation of the network reliability of a stochastic-flow network with multiple sources.

The Taiwan Advanced Research and Education Network (TWAREN) [24] is Taiwanโ€™s academic research network that mainly provides network communication services for Taiwanโ€™s research and academic society. It also offers a tunnel between Taiwan and the United States to connect the global research network through a land surface line and the Asia Pacific regionโ€™s submarine cable. Since TWARENโ€™s resources (i.e., bandwidth) are limited, it is a critical issue to find a technique to optimize its utility. Using efficient evaluation tools to understand TWARENโ€™s performance to improve its infrastructure is one of the major tasks of Taiwanโ€™s National High Performance Computing Center (NCHC). To measure TWARENโ€™s capability, network analysis is a useful tool. For a practical computer network, the transmission media (physical lines such as fiber optics or coaxial cables) may be modeled as arcs, while transmission facilities (switches or routers) may be modeled as nodes. In particular, the capacity of each arc should be stochastic due to the possibility of failure, partial failure, or maintenance. Thus, the computer network characterized by such arcs also has stochastic capacities and it is a typical stochastic-flow network [2โ€“13, 25, 26]. Network reliability evaluation of a stochastic-flow network has been studied as a performance index for decades [2โ€“13, 25, 26]. Most of these studies examined the network reliability from source node ๐‘  to sink node ๐‘ก in terms of minimal paths (MPs), in which an MP is a path with proper subsets which are no longer paths [2โ€“4, 6โ€“9]. This implies that an MP is a set that connects an (๐‘ , ๐‘ก) pair, here not limited to one (๐‘ , ๐‘ก) pair, without any surplus arcs from the perspective of the network topology.

Those previous studies assume that data can be sent through all possible MP from ๐‘  to ๐‘ก according to the network topology, where each MP is composed of some physical lines (arcs). However, in a real computer system, data can only be sent through some unique light paths (LPs) between specific node pairs, where an LP is a virtual tunnel between two end-to-end nodes which combined by some segments (i.e., arcs or lines) and nodes; however, an MP is a path that connects a specific source and a specific sink, while an LP can be a link between any two nodes (not limited to source and sink pair). That is, data may be transmitted from source node ๐‘  to sink node ๐‘ก via more than one LP. In particular, any segments that LP goes through cannot be divided during transmission. Therefore, the previous studies [2โ€“4, 6โ€“9] based on MP to transmit data are not appropriate for TWAREN. In TWAREN, each LP is composed of a set of light path segments (LPSs) linking two nodes. In particular, each physical line can be divided into several LPS, and each LPS belongs to only one LP. Since TWAREN involves the light path, which cannot be divided through any part of its nodes or arcs during transmission, this kind of network model is different than the MP concept described in [2โ€“4, 6โ€“9]. Therefore, we implement a minimal light path (MLP) concept to find all LPs to evaluate TWARENโ€™s network reliability. In this paper, the MLP is defined as a series of nodes and LPSs, from source node to sink node, which contains no cycle.

A revised stochastic-flow network with multiple sources is constructed to model the TWAREN in terms of LP. The difference of single source and multiple sources is that previous one only dedicates on the network reliability between one source and one sink. But in real world, system may transfer the real time data to the sink that exceed the total capacity of LPS which are beside the single source node. That is, we need to transfer data from at least two source nodes, where the data from different source might influence each other, the theory that developed in traditional one source and one sink not applicable here. In generally, we have to transfer the real-time data from multiple sources to one sink to handle the practical worldsโ€™ data transmission. Therefore we consider multiple sources and implement the new technique in this paper to realize the operation of real system instead of single source. Besides, we need to deal with the assignment of multiple sources and the flow conflict on the intersectional arcs. The two-source case is firstly addressed for convenience. A general case with multiple sources can be extended by the proposed algorithm. Then we can evaluate the network reliability for the international part of TWAREN whose tunnel mainly connects to the global academic research network, especially the Internet2 Network [27]. Taiwanโ€™s largest network service provider (NSP), Chunghwa Telecom (CHT) [28], integrates those NSPs that the lines pass through to organize the whole portion of TWARENโ€™s international infrastructure in two areas: on the land surface of both Taiwan and the United States and in the under-sea areas of the Asia Pacific, including the Japan-US submarine cable that disconnected when it was hit by the earthquake and tsunami in Japan on March 11th 2011. Nakagawa [29] has mentioned the influence of that earthquake regarding reliability, so we study this disasterโ€™s effect as well. In fact, when a line breaks, the NSP of these pass-through lines will offer serviceable lines as backups; therefore they offer some degree of the network reliability. However, in this study, we only concentrate on the portion that includes the regular lines to determine the factors that affect TWARENโ€™s network reliability, as the NCHCโ€™s prime task, aside from improving TWARENโ€™s overall performance, is to anticipate major factors which could fail the regular lines. The issue of the network reliability of the backup line [30โ€“34] has not been considered yet.

This paper mainly emphasizes the network reliability that the network can send specified units of data from two source nodes (Taipei city and Hsinchu city) to a single sink node (New York) through TWARENโ€™s light path. The remainder of this paper is organized as follows. The TWAREN network is introduced in Section 2. The research scope, problem formulation, the concept of the minimal light path and the evaluation technique, recursive sum of disjoint of products, (RSDP [9]) are all described in Section 3. Network reliability of TWAREN is evaluated in Section 4. A summary and conclusion are presented in Section 5.

2. TWAREN Network

2.1. Introduction to TWAREN

TWAREN has been funded by the National Science Council of Taiwan since 1998 and was built by the NCHC. Construction was completed at the end of 2003, and service and operation started in 2004. Today, more than 100 academic and research institutions connect with TWAREN in Taiwan and this number is increasing continuously. As well, since 2005, over 1,000 elementary schools and junior and senior high schools have been using TWARENโ€™s internal backbone. TWAREN provides network infrastructure for general use but is also an integrated platform for network research. For instance, TWAREN was instrumental in developing applications and network technology such as IPv6, MPLS, VoIP, e-learning, multicast, multimedia, and performance measurement and has supported GRID computing applications such as e-Learning Grid, Medical Grid, and EcoGrid. As promoting Taiwan as an international R&D center is one of NCHCโ€™s objectives, a stable and reliable TWAREN is the foundation to achieve this goal.

Many countries fund national research and education network (NREN) infrastructure. TWAREN, Taiwanโ€™s NREN, connects to the international research community through global advanced networks, specifically the Internet2 Network [27] of the United States, the major NREN in the world. Therefore, network reliability analyses of TWAREN will help to continuously improve its infrastructure so it can continue to cooperate and connect globally.

2.2. TWARENโ€™s Light Path

TWAREN is network that connects to the world-wide research network through light path international tunnel. TWARENโ€™s physical topology is an optical infrastructure and its virtual topology is constructed by connecting light paths and routers. A light path is a tunnel between two sites connected by various cables and is an end-to-end, preallocated optical network resource, according to usersโ€™ needs. It allows signals to be delivered sequentially without jitters and congestion. Each light path is generally a 155โ€‰Mbps~10โ€‰Gbps dedicated channel that transports various applications.

Figure 1 is the light path international infrastructure that TWAREN leases from CHT, including major sites located at Taipei and Hsinchu in Taiwan, and Los Angeles, Chicago, and New York in the United States. This infrastructure contains the land surface and submarine cable between these cities. Each light path is denoted by LP๐‘– where ๐‘– is the light path number, ๐‘–=1,2,โ€ฆ,๐‘™ with ๐‘™ being the number of light path.

Most of these city sites connect to each other with 2.5โ€‰Gbps physical line connections, divided into four light path channels at 622โ€‰Mb bandwidths. The research scope of this paper is to study the network reliability of the transmission from two sources (Taipei city and Hsinchu city) to the sink node (New York) by means of the light path tunnel.

3. Problem Description and Model Formulation

3.1. Problem Description

This paper describes how the probability that a specified amount of data can be sent from Taipei and Hsinchu to New York via TWAREN is measured. This is referred to as network reliability. Also, Figure 1 is transformed into Figure 2 which is constructed by the light path segments and nodes.

3.2. Some Definitions

As Figure 2 shows, those cities or site devices defined as nodes are denoted by ๐‘›๐‘˜, where ๐‘˜=1,2,โ€ฆ,๐‘ with ๐‘ being the number of nodes. For example, Taipei City is ๐‘›1 and TP-1 is ๐‘›2. We denote each LPSs as ๐‘™๐‘–,๐‘— where ๐‘™๐‘–,๐‘—โˆˆLP๐‘– means the ๐‘—th segment in LP๐‘– (๐‘—=1,2,โ€ฆ,๐‘Ÿ๐‘– with ๐‘Ÿ๐‘– being the number of LPS in LP๐‘–). For example, in Figure 2, LP1 is a tunnel from Taipei (๐‘›1) to Chicago (๐‘›8), which is combined with three LPS ๐‘™1,1,๐‘™1,2, and ๐‘™1,3, and goes through two nodes ๐‘›2 (TP-1) and ๐‘›6 (San Francisco). Its connection sequence is ๐‘›1โ†”๐‘™1,1โ†”๐‘›2โ†”๐‘™1,2โ†”๐‘›6โ†”๐‘™1,3โ†”๐‘›8. The capacity of each LP is 622โ€‰Mb, and each LP is combined by four 155โ€‰Mb channels. As each channel is regarded as one unit, there are 4 units for each LPs.

The physical line (PL) is the actual optical cable where the LP is located and used for data transmission. For example, LPS ๐‘™1,3,๐‘™4,4,๐‘™11,1 is combined in one PL from San Francisco to Chicago, as shown as PL ๐‘ƒ10 in Figure 3. The capacity of each PL is 2.5โ€‰G and is divided into four 622โ€‰Mbโ€‰LP.

The capacity state of an LPS is the same as a PL either when connected or disconnected. Each LPS has two capacity states: 0 units (0โ€‰G) and 4 units (622โ€‰Mb with four 155โ€‰Mbโ€‰LP), respectively. That is, once the PL fails, all the LPSs that are located in this PL also fail. Those LPSs located in the same PL have the same disconnection probability (or, conversely, the same connection probability). For example, LPS ๐‘™1,3,๐‘™4,4,๐‘™11,1 located in one PL ๐‘ƒ10 have the same disconnection probability.

3.3. Model Formulation

The stochastic-flow network evaluation technology developed in [3] is a method that is not suitable to be applied to TWAREN in Figure 2. There are some differences in this problem, since each LP๐‘– is combined with LPS ๐‘™๐‘–,๐‘—, which cannot be divided through any nodes. To create an easier expression, we re-sort all LPSs as ๐‘Ž1,๐‘Ž2,โ€ฆ,๐‘Ž๐‘›, where ๐‘› is the total number of LPS, instead of ๐‘™๐‘–,๐‘—. Let ๐บ=(๐ด,๐‘,๐‘€) be a stochastic flow network where ๐ด={๐‘Ž๐‘–โˆฃ1โ‰ค๐‘–โ‰ค๐‘›} is the set of LPS, ๐‘ is the set of nodes, and ๐‘€=(๐‘€1,๐‘€2,โ€ฆ,๐‘€๐‘›) with ๐‘€๐‘– (an integer) being the maximum capacity of each LPS ๐‘Ž๐‘–. Such a ๐บ is assumed to further satisfy the following assumptions.(1)Each node is perfectly reliable.(2)The capacity of each LPS is stochastic with a given probability distribution according to historical data.(3)The capacities of different LPS are statistically independent.

Let Taipei be the first source node denoted by ๐‘ 1, and let Hsinchu be the second source node denoted by ๐‘ 2. Then let ๐‘†={๐‘ 1,๐‘ 2}. A minimal light path (MLP) is a series of LPSs from a source node to a sink node, which contains no cycle. In particular, any segment used by LP๐‘– cannot be divided during transmission in LP๐‘–. That is, each LPS belongs to only one LP. Suppose ๐‘š๐‘™1,๐‘š๐‘™2,โ€ฆ,๐‘š๐‘™๐‘Ÿ are all MLPs from ๐‘ 1 to ๐‘ก and ๐‘š๐‘™๐‘Ÿ+1,๐‘š๐‘™๐‘Ÿ+2,โ€ฆ,๐‘š๐‘™๐‘ž are all MLPs from ๐‘ 2 to ๐‘ก. Then, the stochastic flow network can be described by the capacity vector ๐‘‹=(๐‘ฅ1,๐‘ฅ2,โ€ฆ,๐‘ฅ๐‘›) and the flow vector ๐น=(๐‘“1,๐‘“2,โ€ฆ,๐‘“๐‘ž) where ๐‘ฅ๐‘– denotes the current capacity of ๐‘Ž๐‘–, and ๐‘“๐‘— denotes the current flow on ๐‘š๐‘™๐‘—. The following constraint shows that the flow through ๐‘Ž๐‘– cannot exceed the maximum capacity of ๐‘Ž๐‘–:๐‘ž๎“๐‘—=1๎€ฝ๐‘“๐‘—โˆฃ๐‘Ž๐‘–โˆˆ๐‘š๐‘™๐‘—๎€พโ‰ค๐‘€๐‘–.(3.1) Let the total demand to New York be ๐‘. Then demand set ๐ท๐‘={(๐‘‘1,๐‘‘2)โˆฃ(๐‘‘1+๐‘‘2)=๐‘} where ๐‘‘1 and ๐‘‘2 are the demand from Taipei and Hsinchu to New York, respectively. To meet the demand pair (๐‘‘1,๐‘‘2), the flow vector ๐น=(๐‘“1,๐‘“2,โ€ฆ,๐‘“๐‘ž) has to satisfy๐‘Ÿ๎“๐‘—=1๐‘“๐‘—=๐‘‘1,๐‘ž๎“๐‘—=๐‘Ÿ+1๐‘“๐‘—=๐‘‘2.(3.2) For convenience, let ๐…(๐‘‘1,๐‘‘2) = {๐นโˆฃ๐น satisfy constraints (3.1) and (3.2)}. For each ๐นโˆˆ๐…(๐‘‘1,๐‘‘2), the corresponding capacity vector ๐‘‹๐น=(๐‘ฅ1,๐‘ฅ2,โ€ฆ,๐‘ฅ๐‘›) is generated via๐‘ฅ๐‘–=๐‘ž๎“๐‘—=1๎€ฝ๐‘“๐‘—โˆฃ๐‘Ž๐‘–โˆˆ๐‘š๐‘™๐‘—๎€พ,๐‘ฅ๐‘–=1,2,โ€ฆ,๐‘›.(3.3) Let ฮฉ(๐‘‘1,๐‘‘2)={๐‘‹๐นโˆฃ๐นโˆˆ๐…(๐‘‘1,๐‘‘2)} be such capacity vectors, and let ฮฉ(๐‘‘1,๐‘‘2),min={๐‘‹โˆฃ๐‘‹ be โ‰ค with respect to in ฮฉ(๐‘‘1,๐‘‘2)} (where ๐‘Œโ‰ค๐‘‹ if and only if ๐‘ฆ๐‘–โ‰ค๐‘ฅ๐‘– for each ๐‘–=1,2,โ€ฆ,๐‘› and ๐‘Œ<๐‘‹, if and only if ๐‘Œโ‰ค๐‘‹ and ๐‘ฆ๐‘–<๐‘ฅ๐‘– for at least one ๐‘–). For convenience, each ๐‘‹โˆˆฮฉ(๐‘‘1,๐‘‘2),min is named a (๐‘‘1,๐‘‘2)-MLP in this paper. Suppose all MLPs have been precomputed. All (๐‘‘1,๐‘‘2)-MLP can be derived by the following steps.

Step 1. Do the following steps for each (๐‘‘1,๐‘‘2)โˆˆ๐ท๐‘.

Step 2. Find all feasible solutions ๐น=(๐‘“1,๐‘“2,โ€ฆ,๐‘“๐‘ž) of the constraints (3.1) and (3.2).

Step 3. Transform each ๐น into ๐‘‹๐น=(๐‘ฅ1,๐‘ฅ2,โ€ฆ,๐‘ฅ๐‘›) via (3.3) to get ฮฉ(๐‘‘1,๐‘‘2).

Step 4. Remove the nonminimal ones in ฮฉ(๐‘‘1,๐‘‘2) to obtain ฮฉ(๐‘‘1,๐‘‘2),min, that is, (๐‘‘1,๐‘‘2)-MLP.

Step 5. Next (๐‘‘1,๐‘‘2).

Step 6. End.

3.4. Network Reliability Evaluation

Network reliability ๐‘…๐ท๐‘ is the probability that the system can transmit ๐‘ units of data to the sink, that is, ๐‘…๐ท๐‘=โˆ‘(๐‘‘1,๐‘‘2)โˆˆ๐ท๐‘Pr{๐‘Œโˆˆฮฉ(๐‘‘1,๐‘‘2)โˆฃ๐‘Œโ‰ฅ๐‘‹ for a (๐‘‘1,๐‘‘2)-MLP๐‘‹}. If {๐‘‹1,๐‘‹2,โ€ฆ,๐‘‹โ„Ž} is the set of minimal capacity vectors capable of satisfying any (๐‘‘1,๐‘‘2)โˆˆ๐ท๐‘, then network reliability ๐‘…๐ท๐‘ is๐‘…๐ท๐‘๎ƒฏ=Prโ„Ž๎š๐‘ฃ=1๐‘„๐‘ฃ๎ƒฐ,(3.4) where ๐‘„๐‘ฃ={๐‘‹โˆฃ๐‘‹โ‰ฅ๐‘‹๐‘ฃ}, ๐‘ฃ=1,2,โ€ฆ,โ„Ž. Several methods such as the RSDP algorithm (Algorithm 1) [9], the inclusion-exclusion method (IEM) [10, 25], the disjoint-event method (DEM) [35], and state-space decomposition (SSD) [11, 12] may be applied to compute ๐‘…๐ท๐‘. The IEM [10, 25] principle is a simple way to calculate network reliability, which basically is similar to the theorem in traditional probability theory that is recursively plus (inclusion) and minus (exclusion) the intersection portion, but easily results in memory overload as there are lots of input data. SSDs [12] are based upon the decomposition method, in which the state space is decomposed into three sets of states: acceptable (๐ด) sets, nonacceptable (๐‘) sets, and unspecified (๐‘ˆ) sets, which recursively decompose the ๐‘ˆ sets into smaller ๐ด, ๐‘, and ๐‘ˆ sets to get the whole system reliability in terms of the summation of the reliability of all A sets. Aven [12] proved that somehow SSD has much better performance than IEM [10, 25]. Zuo et al. [9] implemented a new technique RSDP; it calculates one recordโ€™s reliability first and then continuously and, respectively, handles another single record that is minus the intersection portion with previous records that those reliability already been calculated, which quite different than the IEM that recursively plus and minus the intersection portions for all records. It has been proved by Zuo et al. [9] that RSDP has better efficiency than SSD [12] and easier than IEM [10, 25]. Therefore, recently most network reliability evaluation articles apply the RSDP to assess the related issue. It calculates the probability of a union with ๐‘Ÿ vectors in terms of the probabilities unions with (๐‘Ÿโˆ’1) vectors or less by using a special maximum operator [9] โ€œโŠ•โ€, which is defined as๐‘‹1,2=๐‘‹1โŠ•๐‘‹2โ‰ก๎€ท๎€ท๐‘ฅmax1๐‘–,๐‘ฅ2๐‘–๎€ธ๎€ธ,for๐‘–=1,2,โ€ฆ,๐‘›.(3.5) For example, if ๐‘‹1 = (2, 2, 1, 1, 0) and ๐‘‹2 = (3, 0, 1, 0, 1), ๐‘‹1,2=๐‘‹1โŠ•๐‘‹2 = (max(2,3), max(2,0), max(1,1), max(1,0), max(0,1)) = (3, 2, 1, 1, 1). The RSDP algorithm is presented as follows.

//Calculate the network reliability for all
function = RSDP
//Input vectors and connection probability of each LPS
โ€ƒfor i = 1โ€‰:โ€‰h
โ€ƒโ€ƒif i == 1
โ€ƒโ€ƒโ€ƒ = Pr( );
โ€ƒโ€ƒelse
โ€ƒโ€ƒโ€ƒTemp_R_1 = Pr ;
โ€ƒโ€ƒโ€ƒIf i == 2
โ€ƒโ€ƒโ€ƒโ€ƒTemp_R_2 = Pr(X โ‰ฅ max( )); //max( ) = ( โŠ• )
โ€ƒโ€ƒโ€ƒelse
โ€ƒโ€ƒโ€ƒโ€ƒfor j = 1โ€‰:โ€‰iโ€“1
โ€ƒโ€ƒโ€ƒโ€ƒโ€ƒ = max( ); //max( ) = ( โŠ• )
โ€ƒโ€ƒโ€ƒโ€ƒend
โ€ƒโ€ƒโ€ƒโ€ƒTemp_R_2 = RSDP ;
โ€ƒโ€ƒโ€ƒend
โ€ƒโ€ƒend
โ€ƒโ€ƒ = + (Temp_R_1) โˆ’ (Temp_R_2);

4. Case Study: TWAREN between Taiwan and the US through the Light Path

4.1. Level of Demand and MLP from Taipei and Hsinchu to New York

To calculate TWARENโ€™s network reliability from Taipei and Hsinchu to New York, there must be a reasonable demand level. For each arcโ€™s capacity, each LP occupies a bandwidth 622โ€‰Mb, and each 622โ€‰Mb bandwidth has four 155โ€‰Mb channels. We regard each 155โ€‰Mb as one unit. Therefore, there are four units in each 622โ€‰Mbโ€‰LP channel.

Let the total demand be ๐‘=20 units, that is, 5 ร— 622โ€‰Mb = 3,110โ€‰Mb. For ฮฉ(๐‘‘1,๐‘‘2),min the demand set ๐ท20={(20,0),(16,4),(12,8),(8,12),(4,16),(0,20)}, we try to evaluate ๐‘…๐ท20=โˆ‘(๐‘‘1,๐‘‘2)โˆˆ๐ท20Pr{๐‘Œโˆˆฮฉ(๐‘‘1,๐‘‘2)โˆฃ๐‘Œโ‰ฅ๐‘‹for๐‘Ž(๐‘‘1,๐‘‘2)-MLP ๐‘‹}. In these cases, there are 10 MLPs from ๐‘›1 (Taipei) to ๐‘›9 (New York) as Table 1(a) and 10 MLPs from ๐‘›13 (Hsinchu) to ๐‘›9 (New York) as shown in Table 1(b).

4.2. Probability of All LPSs Breaking

To compute the connection probability of each PL, we use the disconnection data from 2008 through 2011. The longest duration of every break for each physical line during the 168 hours of every week is used to determine the disconnection probability of each line. For example, as the physical line ๐‘ƒ10 from San Francisco to Chicago broke for 403 minutes on 2010/5/25, its connection probability is (168 ร— 60 โˆ’ 403)/(168 ร— 60) = 0.90. Therefore, its disconnection probability is (1 โˆ’ 0.9) = 0.1. All the LPSs ๐‘™1,3,๐‘™4,4 and ๐‘™11,1 located in this physical line ๐‘ƒ10 have the same disconnection probability of 0.1.

Table 2 shows all LPSsโ€™ connection probability after screening all physical linesโ€™ disconnection records and selecting the longest broken time for each. These breaks include disabled card devices, circuit failures, and breaks from March 11, 2011 Japanese earthquake and tsunami that caused the physical submarine line ๐‘ƒ8 to break. This line uses a submarine cable connection between TP-3 and Los Angeles. Artificial devices, short circuits, and natural disasters simultaneously influence TWARENโ€™s network reliability from Taipei and Hsinchu to New York. Since each failure of a node device has been included and recorded in the physical lineโ€™s disconnection record, each node is supposed to be perfect with a reliability of 1. For computational convenience, as described in Section 3.3, we converted LPS ๐‘™๐‘–,๐‘— by using ๐‘Ž๐‘– and the probability of ๐‘Ž๐‘–, as Table 3 shows.

4.3. Network Reliability Computation

When line breaks occur, the suppliers of these pass-through physical lines provide all serviceable lines as backup lines, therefore increasing the network reliability. In this study, we do not discuss the backup lines and concentrate only on the regular lines to determine those factors that affect their network reliability. Firstly, we focus on the demand set ๐ท20={(20,0),(16,4),(12,8),(8,12),(4,16),(0,20)}, given all MLPs in Tables 1(a) and 1(b) and by using the algorithm in Section 3.3 as follows to obtain ฮฉ(๐‘‘1,๐‘‘2),min.

Step 1. Do the following steps for (4,16) โˆˆ๐ท20 (since there is no solution for ฮฉ(0,20),min in this example, we only demonstrate ฮฉ(4,16),min here).

Step 2. Find all feasible solutions ๐น that satisfy constraints (4.1): ๐‘“1+๐‘“2+๐‘“9โ‰ค๐‘€1๐‘“=4,1+๐‘“2+๐‘“9โ‰ค๐‘€2โ‹ฎ๐‘“=4,17+๐‘“18+๐‘“19โ‰ค๐‘€32๐‘“=4,17+๐‘“18+๐‘“19โ‰ค๐‘€33๐‘“=4,1+๐‘“2+โ‹ฏ+๐‘“10=๐‘‘1๐‘“=4,11+๐‘“12+โ‹ฏ+๐‘“20=๐‘‘2=16.(4.1)
In this step, each ๐‘“๐‘– has two values, say 0 and 4, standing for the two capacity states of failure or success. From this, we obtain 4 flow vectors as shown in Table 4(a) (column 1).

Step 3. Transform each ๐น into LPS ๐‘‹ to get ฮฉ4,16 by (4.2).
For ๐น1=(0,0,0,0,4,0,0,0,0,0,0,4,4,0,4,0,4,0,0,0), the capacity vector ๐‘‹1 is transformed by ๐‘ฅ1=๐‘“1+๐‘“2+๐‘“9๐‘ฅ=0+0+0=0,2=๐‘“1+๐‘“2+๐‘“9โ‹ฎ๐‘ฅ=0+0+0=0,32=๐‘“17+๐‘“18+๐‘“19๐‘ฅ=4+0+0=4,33=๐‘“17+๐‘“18+๐‘“19=4+0+0=4.(4.2)
Thus, ๐‘‹1=(0,0,0,0,0,0,0,4,4,4,4,0,0,0,0,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4). Similarly, we obtain 4 capacity vectors as shown in Table 4(a) (column 2).

Step 4. The non-minimal ones in ฮฉ(4,16) are removed to obtain ฮฉ(4,16),min, that is, (4,16)-MLP as shown in Table 4(a) (column 3).

When repeating the previous steps, we can also obtain ฮฉ(8,12),min (resp., ฮฉ(12,8),min and ฮฉ(16,4),min) in Table 4(b) (resp. Tables 4(c) and 4(d)). In terms of RSDP [9], we calculate the network reliability ๐‘…๐ท20= โˆ‘(๐‘‘1,๐‘‘2)โˆˆ๐ท20Pr{๐‘Œโˆˆฮฉ(๐‘‘1,๐‘‘2)โˆฃ๐‘Œโ‰ฅ๐‘‹ for a (๐‘‘1,๐‘‘2)-MLP ๐‘‹} = 0.4140. Similarly, ๐‘…๐ท16 = 0.8195, ๐‘…๐ท12= 0.9707, ๐‘…๐ท8 = 0.9976, and ๐‘…๐ท4= 0.9999 can be evaluated, respectively. The network reliability can be observed to decrease as the total demand increases, as shown in Figure 4.

In regard to QoS, this is only a concern when there are insufficient networks resources. When there are enough resources and demand is low, for instance, as above with ๐ท4, there are still plenty of resources to handle other transmission requests, so the network reliability is quite high. On the other hand, if demand is high, say above set ๐ท20, the network reliability will be low, since there are not enough resources to handle other data transmissions. To maintain the network reliability, it is important to avoid full transmission loads or increase line capacity. Depending on the results of our analysis, we may decide to allocate more economic resources to TWAREN to maximize future network utilities.

5. Summary and Conclusion

Instead of the classical TTNR, KTNR, and ATNR analysis of a binary-state flow network, this paper evaluates the network reliability of a stochastic-flow network with multiple sources. It also designs an MLP-based network reliability evaluation technique for the international LP portion of TWARENโ€™s academic and research network. This portion contains the domestic land surface line and the Asia Pacific submarine cables which connect to the global academic research network, including the Internet2 Network [27]. Since the LP cannot be divided through any of its nodes or LPSs during transmission, MLP is a new concept to evaluate the network reliability in an LP environment. MLP is used to discuss the flow assignment and to evaluate the network reliability. This research contributes by making real TWAREN data available to be analyzed in a stochastic-flow network model. By using the MLP analysis technique, we will know how to continuously adjust TWARENโ€™s infrastructure to achieve higher network reliability. In this study, we concentrate on the portion of the network that includes regular lines and does not include backup cables yet. This allows us to determine those factors that influence the dedicated regular linesโ€™ network reliability. We also consider the effects of the earthquake that hit Japan on March 11, 2011. All factors are studied, including artificial, machine, and cable failures and natural disasters that simultaneously influence TWARENโ€™s network reliability from the two source nodes, Taipei and Hsinchu, to the single sink node, New York. In addition, the MLP network reliability technique used in the multiple sources case will enable us to increase the efficiency of TWAREN and help us to learn how to improve its network infrastructure and performance in the near future. Subsequently, further study may be undertaken on the network reliability of TWARENโ€™s multisource to multisink (terminal) issue.

Acknowledgments

This work was supported in part by the National Science Council, Taiwan, Republic of China, under Grant no. NSC 98-2221-E-011-051-MY3. This work was supported in part by the National High Performance Computing Center, Taiwan, Republic of China.