Next Article in Journal
A Study on the Effects of Digital Finance on Green Low-Carbon Circular Development Based on Machine Learning Models
Next Article in Special Issue
Using a Node–Child Matrix to Address the Quickest Path Problem in Multistate Flow Networks under Transmission Cost Constraints
Previous Article in Journal
On Erdélyi–Kober Fractional Operator and Quadratic Integral Equations in Orlicz Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Design Problem for Multi-Source Multi-Sink Flow Networks Based on Genetic Algorithm Approach

by
Sahbi Boubaker
1,*,
Noha Hamdy Radwan
2,
Moatamad Refaat Hassan
2,
Faisal S. Alsubaei
3,
Ahmed Younes
4,5 and
Hameda A. Sennary
2
1
Department of Computer & Network Engineering, College of Computer Science and Engineering, University of Jeddah, Jeddah 21959, Saudi Arabia
2
Department of Mathematics and Computer Science, Faculty of Science, Aswan University, Aswan 81528, Egypt
3
Department of Cybersecurity, College of Computer Science and Engineering, University of Jeddah, Jeddah 23218, Saudi Arabia
4
Department of Computer Science, College of Applied Studies and Community Service, Imam Abdulrahman bin Faisal University, Dammam 34212, Saudi Arabia
5
Department of Computer Science, Faculty of Computer and Artificial Intelligence, Sohag University, Sohag 82524, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(18), 3902; https://doi.org/10.3390/math11183902
Submission received: 12 August 2023 / Revised: 7 September 2023 / Accepted: 8 September 2023 / Published: 13 September 2023
(This article belongs to the Special Issue Graph Theory and Network Theory)

Abstract

:
Robust design problems in flow networks involve determining the optimal capacity assignments that enable the network to operate effectively even in the case of events’ occurrence such as arcs or nodes’ failures. Multi-source multi-sink flow networks (MMSFNs) are frequent in many real-life systems such as computer and telecommunication, logistics and supply-chain, and urban traffic. Although numerous studies on the design of MMSFNs have been conducted, the robust design problem for multi-source multi-sink stochastic-flow networks (MMSFNs) remains unexplored. To contribute to this field, this study addresses the robust design problem for MMSFNs using an approach of two steps. First, the problem is mathematically formulated as an optimization problem and second, a sub-optimal solution is proposed based on a genetic algorithm (GA) involving two components. The first component, an outer genetic algorithm, is employed to search the optimal capacity assigned to the network components with minimum sum. The second component, an inner genetic algorithm, is used to find the optimal flow vectors that maximize the system’s reliability. Through extensive experimentation on three different networks with different topologies, the proposed solution has been found to be efficient.

1. Introduction

Multi-source multi-sink stochastic-flow networks (MMSFNs) play a vital role in real-life systems, such as telecommunication, computer, logistics, supply-chain, and transportation networks. The optimization of such networks has been a popular topic for applied research for the last few years. The main concept behind these kinds of networks is related to optimizing both the network capacity and reliability under several constraints such as budget. As an example, Hassan [1] presented a strategy to maximize the reliability of the capacity vector for MMSFNs, considering an assignment budget. The strategy was divided into two parts: the first part searches for the preferable components that can be specified in the network with minimum cost, while the second part searches for the flow arc with maximum reliability of the capacity arc for the assigned components. Hassan [2] solved the problem of system reliability optimization for MMSFNs, subject to transmission budget constraints, by utilizing genetic algorithms to achieve maximum system reliability through obtaining the best optimal set of lower boundary points. Forghani-elahabad and Kagan [3] presented an algorithm to calculate the exact value of reliability in terms of minimal number of paths in MMSFNs. Network reliability is defined as the probability that a gathering of sinks is fed by a gathering of sources [4]. It refers to the ability of a network to fulfil its purpose within the specified timeframe, despite environmental factors, as determined by the connection between a source and a target [5].
Hamid et al. [6] employed a genetic algorithm (GA) to assign optimal capacities with minimum total capacities and maximum reliability. The authors presented a GA to optimize the reliability of MMSFN systems while solving flow allocation problems [7]. Liu et al. [8] used the GA to find the optimal allocation pattern (vector and node) by the implementation of the equations in the algorithm. The formulas sum up the flow quantity in each minimal pass when the node or arc belongs to this minimal pass. Chen [9] introduced an algorithm to achieve the optimal double-resource assignment for the robust design problem in multi-state computer networks. Atzori and Raccis [10] used a genetic algorithm to assign capacity for multicast services.
From the previous selected studies, it can be noticed that the concept of robustness was not considered when solving the MMSFN problem. For this reason, the current study will concentrate on solving these kinds of problems under constraints and robustness considerations.
Robust Design is defined as a powerful mechanism for improving the reliability of a system at a low cost in a short time [11,12]. It is a process used to find the best configuration for a product design in order to fulfil the end-user’s needs while considering feasibility requirements [13,14]. Robust Design (RD) is defined as an engineering approach that concentrates on the development/improvement of new/existing products (physical products or services), processes, and specific equipment. The robustness of a design refers to the insensitivity to different disturbances and variations coming from the surrounding environment [15]. Robust design, based on the concept of building quality into products or processes, is considered one of the most important systems engineering design concepts for improving quality and the optimization process. In this direction, Cho and Shin [16] presented new robust design empirical and optimization models for time-oriented data by developing a set of methods following three separate phases: modeling, estimation, and optimization.
The robust design problem in a capacitated flow network aims to achieve the minimum capacity allocated to each arc so that the network continues surviving even under the arc’s failure and external disturbances that may affect its operation [17]. Hassan and Abdou [16] conducted a study to evaluate the MMSFNs’ reliability under time constraints. The studies by [18,19] solved the robust design problem for single and two-commodity flow networks with node failure by using a genetic algorithm. Chen and Lin [20] proposed a solution to find feasible flows that meet flow requirements while minimizing the maximum occurring cost among all demand realizations. López-Prado [21] defined the robustness not only in terms of the algorithm’s ability to supply “good” solutions to problem instances of different sizes and numerical characteristics but also in terms of its ability to keep working well even when constraints are added or removed. A successful robust design should ensure a compromise between feasibility and optimality.
Reliability evaluation is the process of determining whether an existing system has achieved a particular level of operational reliability. For example, the reliability evaluation of electric power systems is a fundamental and vital problem in the planning, design, and operation of power systems. The reliability assessment in a conventional power system is the probability that it implements its functions properly, without any failure within a stipulated period of time, when it is subject to normal operating conditions [22]. In transport networks, reliability is a useful tool to represent its operational goodness [23]. Network reliability is important to measure the degree of stability and quality of infrastructure assembly of a network [24].
In order to carefully cope with the problems usually faced during the design phase of a network, this study’s main contributions are the formulation of the robust design problem for MMSFNs and its sub-optimal solution using efficient mathematical tools and proposing suitable algorithms. The solution of the problem is built around the concept of robustness defined as the ability of a network to operate even under uncertainties and failures. In addition, an approach based on genetic algorithms (GAs) is developed to solve the problem in MMSFNs by minimizing the sum of arcs capacities and searching for the best set of flow vectors that check the conditions and achieve maximum reliability. The proposed approach is illustrated through three networks having different topologies to show its efficiency. To the best of the authors’ knowledge, this study is the first that considers the concept of robustness in the design problem of MMSFNs.
The remaining parts of this paper are organized as follows: in Section 2, materials and methods including the structural analysis, problem formulation, the proposed approach for solving the problem, and the whole algorithm are presented in Section 2. The results and discussion are provided in Section 3. Finally, Section 4 concludes the paper.

2. Materials and Methods

2.1. Notations

  • Network notations:
npNo. of paths
ncNo. of arcs
vNo. of nodes
σNo. of sources
θNo. of sinks
AThe set of arcs {a_b│1 ≤ bnc}
SThe set of source node {s_1,s_2,…,s_σ}
TThe set of sink node {t_1,t_2,…,t_θ}
MPMinimal path
MPi, j, kThe kth minimal path from source i to sink j
mNo. of resources.
sdw,jThe demand for resource w at sink node tj
DDThe demand DD = {sdw,jǀ1 ≤ wm, 1≤ jθ}
srw,iThe maximum quantity of resource w that source node si can supply
RRThe resource RR = {srw,iǀ1≤ wm, 1 ≤ iσ}
MiMaximum capacity of an arc i
FFlow vector defined as F = (f_1,1,1,1,f_1,1,2,1 …, f_(i,j,k_(i,j,) 1), …, f_(i,j,k_(i,j),m), …,f_(σ,θ,k_(σ,θ),m)), where f_(i,j,k,w) represents the flow quantity of resource w on MP_(i,j,k)
NFLength of the flow vector, NF = np × m
RSThe system reliability
  • Parameters of the outer GA:
NP1Number of chromosomes
NG1Number of genes equals to nc
MG1The maximum number of generations
Cr1The crossover rate
Mu1The mutation rate
  • Parameters of the inner GA:
NP2Number of chromosomes
NG2Number of genes equals to NF
MG2The maximum number of generations
Cr2The crossover rate
Mu2The mutation rate
XCapacity vector, X = (x_1,x_2,…,x_neq)
RXThe reliability of the capacity vectors

2.2. Assumptions

  • The flow conservation law must apply to the flow network.
  • The arcs have capacities that are statistically independent.
  • The capacity of an arc is an integer-valued random variable, which takes values 0 < 1 < 2 <…< Me according to a given distribution.
  • The flow along a path does not exceed its maximum capacity of that path.

2.3. The Robust Design Problem for MMSFNs: Definitions

The analysis of a network structure is important to assign the optimal capacity to network edges to guarantee that the network will continue surviving even under edges’ failures. In the following subsections, the progressive steps for getting a robust design for MMSFNs as well as the main concepts used during the design will be presented.
  • Coverage set
Suppose that Γ i = { M P i , j , k | a i M P i , j , k } and Γ j = { M P i , j , l | a j M P i , j , l } , the arc a j belongs to the coverage set Q i of a i if and only if Γ j Γ i , [15].
  • Structural impact
The structural impact ( S i ) for a i is given by:
S i = Q i n c  
where   Q i = { a j | Γ j Γ i } . The edge a i has the strongest impact if S i = 1.
  • Critical edge
The arc a i is said to be critical if S i = 1. The network reliability R s is zero if and only if a i has zero capacity [15].
  • The capacity assignment
In a single-source single-sink network, the maximum capacity Mi of a i ranges from 0 to the demand value. However, in MMSFNs, there are multiple demands requested by the sinks. As a result, the following expression is proposed here to set M i values.
Let E i = { s d w , j | a i M P i , j } , then 0 M i E i + δ , δ 0 , N , N is a positive integer.
If a i is a critical edge, then M i should not be less than E i .
  • Probabilities for each edge
The probability for the current capacity γ ( 0 γ M i ) for a i , Pr { γ } :
P r γ = M i γ r i γ 1 r i M i γ
where r i is the probability (availability) of a i [15].

2.4. MMSFN Problem Formulation

The problem is divided into two subproblems, as detailed in the next subsections. The first subproblem is to determine the optimal capacities to be assigned to the network’s arcs. The second one is to find the optimal feasible solutions for the flow vector to evaluate the system’s reliability.

2.4.1. First Subproblem

Let us consider M = ( M 1 , M 2 , , M n c ) as assigned capacities to the set of arcs ( a 1 , a 2 , , a n c ) . The mathematical formulation of the subproblem is provided as follows:
M i n i m i z e   S
M a x i m i z e   R S
where S = i n c M i refers to the sum of assigned capacities and Rs is the corresponding system reliability.

2.4.2. Second Subproblem

Searching the flow vectors that succeed to achieve Equations (5)–(7) (discussed in, [1,2,7,8]).
i = 1 σ k = 1 k i , j f i , j , k , w = s d w , j , w = 1 , , m ; j = 1 , , θ
j = 1 θ k = 1 k i , j f i , j , k , w s r w , i , w = 1 , , m ; i = 1 , , σ
x l M l l = 1 , , n c
where,
x l = i = 1 σ j = 1 θ k = 1 k i , j w = 1 m f i , j , k , w a l ϵ M P i , j , k
The reliability of the capacity vector is given by:
R X = l = 1 n c P r { x x l }

2.5. The Proposed Approach

The proposed approach to solve the design problem is based on a two genetic algorithms (GAs) approach. The first GA (outer) is used to assign the optimal capacities to the arcs and the second GA (inner) is employed to search for the flow vectors that satisfy the conditions in (3)–(5) while maximizing the network reliability. The proposed approach (as required in a GA solution) should include the representation of chromosomes, initial population generation, fitness function, crossover process, and mutation process. Interested readers can refer to [18,25] for details about the GA and its structure, components, and hyperparameters.

2.5.1. The Outer GA

Representation

The representation chromosome is used to characterize each chromosome in the GA. The chromosome M is exemplified by a series of length (nc), where (nc) refers to the number of arcs ( M 1 , M 2 , , M n c ).

Initial Population

The initial population is the first step in the genetic algorithm (GA) and consists of a set of random potential solutions. Each solution is symbolized by a chromosome. Typically, this initial population is generated randomly to produce a range of possible solutions [26]. Although random initialization may allow diversification of solutions over the search space, other initialization methods such as pre-defined known solutions can be used.

Fitness Function

The fitness function is the heart of a genetic algorithm. This function makes it possible to evaluate a given individual solution and determine how well it satisfies the optimization criterion that the algorithm is developed for. The sum of the assigned capacities ( S ) for the ith individual is selected to be the fitness function of this candidate solution, i.e., f i t ɨ = S . To normalize the fitness function for each individual i, the calculated value is divided by the total sum of all fitnesses ( i f i t ɨ ) in the population. The following Algorithm 1 outlines how to calculate the fitness for the ith solution in the population:
Algorithm 1: Fitness evaluation
f i t ɨ = S
  f o r   ɨ = 1 : N P 1
   N o r m a l i z e _ f i t ( ɨ ) = f i t ( ɨ ) / i f i t ɨ
  E n d   f o r

Selection

The selection process is one of the important steps in genetic algorithms, used to choose individuals based on their fitness value. Chromosomes with higher fitness values are more likely to be selected for reproduction, while those with lower values have a lower chance of being selected. In other words, the probability of a chromosome being chosen for reproduction is proportional to its fitness value [27].
This study employs the roulette wheel selection method, a widely adopted technique in genetic algorithms to identify promising individuals for the crossover and mutation process (Algorithm 2). In the roulette wheel selection, like all other selection methods, each prospective solution is assigned a solution through the fitness function. This measure of physical fitness determines the likelihood of an individual’s selection. Solutions with superior fitness have a higher probability of being selected, while inferior solutions may still have the opportunity to endure through the selection process. This characteristic may incorporate certain elements that could prove advantageous in subsequent recombination processes. In this study, the roulette wheel mechanism employed for parent selection relies on the cumulative sum of the respective individual’s fitness values.
Algorithm 2: Roulette wheel algorithm
Begin
f o r   j = 1   t o   N P 1
f o r   i = j   t o   N P 1
  C S u m j = C S u m j + N o r m a l i z e _ f i t i      //CSum stands for the commutative sum
E n d   f o r   i
E n d   f o r   j
Generate number random ƙ   ϵ   [ 0 , 1 ]
f o r   j = 1   t o   N P 1
i f   ƙ > C S u m ( j )
p a r e n t   1 = j 1
e n d   i f
e n d   f o r   j
Generate number random ƙ ϵ [ 0 , 1 ]
f o r   j = 1   t o   N P 1
i f   ƙ > C S u m ( j )
p a r e n t   2 = j 1
e n d   i f
e n d   f o r   j
E n d

Crossover

Crossover is a genetic operator that combines two chromosomes to produce a new chromosome. The idea behind crossover is that the new chromosome may inherit the best characteristics from each of the parent chromosomes.
In this study, a one-point crossover operator that randomly chooses one crossover point is used. The process involves copying everything before this point from the first parent and everything after the crossover points from the second parent, as follows [28].
M C = [ M p 1 ( j ) ] j = 1 α + [ M p 2 ( j ) ] j = α + 1 n c M D = [ M p 2 ( j ) ] j = 1 α + [ M p 1 ( j ) ] j = α + 1 n c
where M C and M D are the children and M p 1 and M p 2 parents selected based on Mc1 value and α is the random cut-point.

Mutation

Mutation represents the final mechanism of evolution, wherein one or multiple genes undergo genetic variation and are subsequently transmitted to offspring for potential adaptation. The genetic algorithm (GA) typically employs a low mutation rate to maintain its functionality since high rates of mutation may lead to a less efficient and more rudimentary random search process (Algorithm 3). The introduction of the mutation operator in the algorithm adds an additional level of randomness, thereby preserving the diversity of the population. The utilization of the genetic algorithm (GA) technique is fruitful in circumventing the emergence of analogous solutions and escalating the likelihood of evading local solutions [29].
Algorithm 3: Genetic algorithm
Generate a random number rm [0,1]
if rm M u 1 then
{
  for i = 1 to nc, do 
   M D ( i ) = β ,   β { 0 , 1 , , M i } .
End for
}

2.5.2. The Inner GA

The inner GA is used to locate the optimal set of lower boundary points that maximize the system reliability. The flowing formulation represents the problem of searching for the optimal lower boundaries to maximize the system reliability.
Find the optimal set of X s . t . R S is maximized
The equations needed to solve flow vectors searching problem are given in Section 2.4.2.

Representation

The chromosome F is represented by a series of length (NF), where (NF) refers to the number resulting from multiplying a number to the minimal path (np) and the number of resources (m) ( f 1,1 , 1,1 , f 1,1 , 2,1 , …, f i , j , k i , j , 1 , …, f i , j , k i , j , m , …, f σ , θ , k σ , θ , m ).

Fitness Function

The following Algorithm 4 shows how to calculate the fitness for each solution (i) in the population, where X i corresponds to F i , and each F i satisfies Equations (5)–(7) presented in Section 2.4.2.
Algorithm 4: Fitness function
Begin
Calculate the fitness value for X i = R X i
f o r   a l l   e l e m e n t s   i n   t h e   c u r r e n t   p o p u l a t i o n
 The N o r m a l i z e d   f i t n e s s   f o r   e a c h   e l e m e n t   i =
         t h e   f i t n e s s   v a l u e   f o r   X i / S u m   o f   a l l   f i t n e s s  
    E n d   f o r
 End

Selection

We used the roulette wheel mechanism in this work to select two parents (Algorithm 5).
Algorithm 5: Roulette wheel algorithm
Begin
f o r   j = 1   t o   N P 2
f o r   i = j   t o   N P 2
C S u m j = C S u m j + N o r m a l i z e _ f i t ( i )
E n d   f o r   i
E n d   f o r   j
Generate number random ƙ   ϵ   [ 0 , 1 ]
f o r   j = 1   t o   N P 2
i f   ƙ > C S u m ( j )
p a r e n t   1 = j 1
e n d   i f
e n d   f o r   j
Generate number random ƙ   ϵ   [ 0 , 1 ]
f o r   j = 1   t o   N P 2
i f   ƙ > C S u m ( j )
p a r e n t   2 = j 1
e n d   i f
e n d   f o r   j
E n d

Crossover

One-point crossover is used to generate new offspring for F as follows.
F C = [ F p 1 ( j ) ] j = 1 α + [ F p 2 ( j ) ] j = α + 1 N F F D = [ F p 2 ( j ) ] j = 1 α + [ F p 1 ( j ) ] j = α + 1 N F
where F C and F D are the new vectors generated by pairing up F p 1 and F p 2 based on Mc2 value and α is the random cut-point.

Mutation

A proposed a mutation mechanism to improve the flow vector F based on Mu2 is provided as follows (Algorithm 6):
Algorithm 6: Mutation mechanism
Generate a random number rm [0,1]
if rm M u 2 then
{
  for i = 1 to NF, do
    F C ( i ) = β ,   β { 0 , 1 , I , M i n c } , M i n c the maximum value can be assigned to F(i), [2].
End for
}

Evaluating Rs

If X 1 , X 2 , , X N P 2 represent the generated set of capacity vectors that correspond to the flow vectors, then by removing the non-minimal ones in X 1 , X 2 , , X N P 2 , we obtain all lower boundary points if the network is cyclic (as discussed in, [30]). Next, if X 1 , X 2 , , X l are all lower boundary points, then the system reliability R S is calculated by Equation (12).
R S = p r   i = 1 l Z Z     X i
p r Z = p r z 1 . p r z 2 . . p r { z n e q } . Then, we use the recursive sum of disjoint products (RSDP) procedure presented in [31].
when   T M 1 = p r Z X 1 and   T M i = p r Z X i p r j = 1 i 1 Z X j , i , f o r   i 2
R s = p r   i = 1 l Z Z     X i = i = 1 l T M i
If Z = z 1 , z 2 , , z e , , z n e q
p r ( Z ) = e = 1 n e q p r z e

2.6. The Whole Algorithm of the Proposed Approach

Algorithm 7 shows the whole algorithm of the proposed approach.
Algorithm 7: The whole algorithm of the proposed approach
Start
Input the network information, such as minimal path, demand, resources, MG1, NP1, NG1, Cr1, Mu1, MG2, NP2, NG2, Cr2, Mu2.
Start outer GAGenerate the initial population randomly with size N P 1 .
Evaluate the initial population.
While g M G 1 , do
While N P 1 , do
                Select two chromosomes.
        Generate new offspring after applying crossover and mutation.
              = + 1
end do
Evaluate the current population.
Save the best solution to M.
    g = g + 1
end do
Report the best solution found.
Start inner GA
    Generate the initial population randomly with size N P 2 .
    Evaluate the initial population.
    While g M G 2 , do
While N P 2 , do
    Select two chromosomes.
     Generate new offspring after applying crossover and mutation.
        = + 1
       end do
 Evaluate the current population.
   Save the best solution to F.
  g = g + 1
 end do
 Report the optimal lowers found.
 End inner GA Calculate the system reliability Rs.
 End outer GA
End

3. Results and Discussion

In the following subsections, the presented approach has been applied to three networks. The valid values of the GA parameters were as follows, MG1 = 100, NP1 = 20, Cr1 = 0.95, Mu1 = 0.05, MG2 = 100, NP2 = 10, 15, and 20. The value of Cr2 = 0.95, and that of Mu2 = 0.05. The proposed approach was implemented in MATLAB 2017a environment with hardware specification that comprises a core i7 Intel processor with 8 GB of RAM.

3.1. Case1: Network with Three Sources and Two Sinks

This network consists of 10 arcs as shown in Figure 1; it has 7 MPs; MP1,1,1 = {a1, a7}, MP1,1,2 = {a2, a9}, MP1,2,1 = {a1, a8}, MP2,1,1 = {a3, a9}, MP2,2,1 = {a4, a10}, MP3,1,1 = {a5, a9} and MP3,2,1 = {a6, a10}. Where nc = NG1 = 10, NF = NG2 = 21, Resources: R = (r1,1, r1,2, r1,3, r2,1, r2,2, r2,3, r3,1, r3,2, r3,3) = (5, 2, 3, 5, 3, 2, 2, 2, 3). Demand: D = (d1,1, d1,2, d2,1, d2,2, d3,1, d3,2) = (3, 1, 2, 2, 1, 3). Table 1 shows the best minimum S achieves maximum R s when NP2 = 10, 15, and 20, and Table 2 clarifies F and R X .

3.2. Case2: Network with Two Sources and Two Sinks

This network consists of 14 arcs as shown in Figure 2, it has 11 MPs; MP1,1,1 = {a1, a5}, MP1,1,2 = {a1, a6, a9}, MP1,1,3 = {a2, a7, a9}, MP1,2,1 = {a1, a6, a14}, MP1,2,2 = {a2, a7, a14}, MP2,1,1 = {a3, a7, a9}, MP2,1,2 = {a4, a8, a13, a9}, MP2,2,1 = {a3, a7, a14}, MP2,2,2 = {a4, a8, a13, a14}, MP2,2,3 = {a4, a8, a10}, and MP2,2,4 = {a4, a11, a12}. Where nc = NG1 = 14, NF = NG2=22, Resources: R = (r1,1, r1,2, r2,1, r2,2) = (15, 17, 10, 13). Demand: D = (d1,1, d1,2, d2,1, d2,2) = (7, 8, 5, 8). Table 3 shows the best minimum S achieves maximum R s when NP2 = 10, 15, and 20, and Table 4 clarifies F and R X .

3.3. Case3: Network with Two Sources and Three Sinks

This network consists of 11 arcs as shown in Figure 3, it has 13 MPs; MP1,1,1 = {a1, a7}, MP1,1,2 = {a2, a5, a7}, MP1,2,1 = {a1, a8}, MP1,2,2 = {a2, a9}, MP1,2,3 = {a2, a5, a8}, MP1,3,1 = {a2, a10}, MP2,1,1 = {a3, a5, a7}, MP2,1,2 = {a4, a6, a5, a7}, MP2,2,1 = {a3, a9}, MP2,2,2 = {a4, a6, a9}, MP2,3,1 = {a3, a11}, MP2,3,2 = {a4, a6, a10}, and MP2,3,3 ={a4, a11}. Where nc = NG1 = 11, NF = NG2 = 26, Resources: R = (r1,1, r1,2, r2,1, r2,2) = (10, 19, 14, 19). Demand: D = (d1,1, d1,2, d1,3, d2,1, d2,2, d2,3) = (3, 2, 2, 2, 3, 3). Table 5 shows the best minimum S achieves maximum R s when NP2 = 10, 15, and 20, and Table 6 clarifies F and   R X .
The main goal of the Robust Design Problem (RDP) for MMSFN is to find the optimal set of capacities assigned to arcs, ensuring network survival even under arc failures, i.e., R S > 0 . Initially, we studied a network with three sources and two sinks. The best values for S and R S were 36 and 0.973585, respectively, with NP2(NF) = 10. Next, in a network with two sources and two sinks, the best values for S and R S were found to be, respectively, 169 and 0.996641 with NP2(NF) = 10. Finally, in a network with two sources and three sinks, the best values for S and R S were 111 and 0.999923, respectively, with NP2(NF) = 15. The NP2 for the inner GA (corresponding to NF, number of flow vectors) was selected to be 10, 15, and 20 for all studied cases to find the optimal or near-optimal value for the system reliability with NG2 equal to 100. The number of flow vectors for MMFNs cannot be anticipated. Therefore, we searched for the optimal flows that maximize the system reliability. According to [1,2,7], the obtained maximum value for NF number does not exceed 20 as summarized in Table 7. Furthermore, we follow up the reliability values of capacity vectors R X for each example, [1], o identify the best flow strategy with maximum R X .

4. Conclusions

This paper presented a structural analysis of MMSFN to discuss and formulate the Robust Design Problem for this type of SFN. Additionally, it successfully solves the RDP for MMSFN using a two-stage GA approach. We divided the problem into two subproblems: first, searching for the optimal capacity assigned to network components with a minimum sum, and second, searching for the optimal lower vectors to achieve the maximum system reliability value. The proposed GA-based approach was applied to a group of networks, and it yielded satisfactory results. No comparisons were made as the RDP for MMSFNs has not been discussed previously.

Author Contributions

Conceptualization, N.H.R., M.R.H. and H.A.S.; methodology, N.H.R., M.R.H. and H.A.S.; software, validation, N.H.R., A.Y., M.R.H. and H.A.S.; formal analysis, S.B. and F.S.A.; investigation, S.B.; writing—original draft preparation, N.H.R., H.A.S. and S.B.; writing—review and editing, A.Y., M.R.H., S.B. and F.S.A.; visualization, A.Y. and M.R.H.; supervision, S.B. and M.R.H.; funding acquisition, S.B. and F.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia—project number MoE-IF-UJ-22-4220772-4.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia, for funding this research work through the project number MoE-IF-UJ-22-4220772-4.

Conflicts of Interest

The authors declare no conflict of interest in this work.

References

  1. Hassan, M.R. Maximizing reliability of the capacity vector for multi-source multi-sink stochastic-flow networks subject to an assignment budget. J. Ind. Manag. Optim. 2021, 17, 1253–1267. [Google Scholar] [CrossRef]
  2. Hassan, M.R. System Reliability Optimization of Multi-Source Multi-Sink Stochastic Flow Networks with Budget Constraint. Int. J. Reliab. Qual. Saf. Eng. 2021, 28, 2150025. [Google Scholar] [CrossRef]
  3. Forghani-Elahabad, M.; Kagan, N. Reliability evaluation of a stochastic-flow network in terms of minimal paths with budget constraint. IISE Trans. 2018, 51, 547–558. [Google Scholar] [CrossRef]
  4. Bobbio, A.; Terruggia, R.; Ciancamerla, E.; Minichino, M. Reliability analysis of multi-source multi-sink critical interacting systems. In Proceedings of the 2011 3rd International Workshop on Dependable Control of Discrete Systems, Saarbruecken, Germany, 15–17 June 2011; pp. 127–132. [Google Scholar] [CrossRef]
  5. Sahinoglu, M.; Rice, B. Network reliability evaluation. In WIREs Computational Statistics; John Wiley & Sons: Hoboken, NJ, USA, 2010; Volume 2, pp. 189–211. [Google Scholar]
  6. Hamed, A.Y.; Alkinani, M.H.; Hassan, M.R. A Genetic Algorithm to Solve Capacity Assignment Problem in a Flow Network. Comput. Mater. Contin. 2020, 64, 1579–1586. [Google Scholar] [CrossRef]
  7. Hassan, M. Solving Flow Allocation Problems and Optimizing System Reliability of Multi-Source Multi-Sink Stochastic Flow Network. Int. Arab. J. Inf. Technol. (IAJIT) 2016, 13, 477–483. [Google Scholar]
  8. Liu, Q.; Zhang, H.; Ma, X.; Zhao, Q. Genetic Algorithm-based Study on Flow Allocation in a Multicommodity Stochastic-flow Network with Unreliable Nodes. In Proceedings of the Eighth ACIS International Conference on Software Engineering, Artificial Intelligence, Networking, and Parallel/Distributed Computing (SNPD 2007), Qingdao, China, 30 July 2007–1 August 2007; Volume 1, pp. 576–581. [Google Scholar] [CrossRef]
  9. Chen, S.-G. Optimal double-resource assignment for the robust design problem in multistate computer networks. Appl. Math. Model. 2014, 38, 263–277. [Google Scholar] [CrossRef]
  10. Atzori, L.; Raccis, A. Network Capacity Assignment for Multicast Services Using Genetic Algorithms. IEEE Commun. Lett. 2004, 8, 403–405. [Google Scholar] [CrossRef]
  11. Yang, G. Life Cycle Reliability Engineering; John Wiley & Sons: Hoboken, NJ, USA, 2007. [Google Scholar] [CrossRef]
  12. Büsing, C.; Koster, A.M.C.A.; Schmitz, S. Robust minimum cost flow problem under consistent flow constraints. Ann. Oper. Res. 2021, 312, 691–722. [Google Scholar] [CrossRef]
  13. Fowlkes, W.Y.; Creveling, C.M. Engineering Methods for Robust Product Design: Using Taguchi Methods in Technology and Product Development; Addison-Wesley Publishing Company: Boston, MA, USA, 1995. [Google Scholar]
  14. Chabrier, A.; Danna, E.; Le Pape, C.; Perron, L. Solving a Network Design Problem. Ann. Oper. Res. 2004, 130, 217–239. [Google Scholar] [CrossRef]
  15. Yong, K. Robust Design and Reliability. Trans. Nanjing Univ. Aeronaut. Astronaut. 1998, 15, 9–14. [Google Scholar]
  16. Cho, B.R.; Shin, S. Quality Improvement and Robust Design Methods to a Pharmaceutical Research and Development. Math. Probl. Eng. 2012, 2012, 193246. [Google Scholar] [CrossRef]
  17. Chen, S.-G. An optimal capacity assignment for the robust design problem in capacitated flow networks. Appl. Math. Model. 2012, 36, 5272–5282. [Google Scholar] [CrossRef]
  18. Hamdy, N.; Hassan, M.R.; Hussein, M.E. A genetic algorithm to solve the robust design problem for a Flow Network with Node Failure. Trans. Networks Commun. 2020, 8, 1–10. [Google Scholar] [CrossRef]
  19. Radwan, N.H.; Hassan, M.R.; Hussein, M.E. Solving the Robust Design Problem for a Two-Commodity Flow Network with Node Failure. Am. J. Eng. Appl. Sci. 2020, 13, 837–845. [Google Scholar] [CrossRef]
  20. Chen, S.-G.; Lin, Y.-K. An Approximate Algorithm for the Robust Design in a Stochastic-Flow Network. Commun. Stat.—Theory Methods 2010, 39, 2440–2454. [Google Scholar] [CrossRef]
  21. López-Prado, J.L.; Vélez, J.I.; Garcia-Llinás, G.A. Reliability Evaluation in Distribution Networks with Microgrids: Review and Classification of the Literature. Energies 2020, 13, 6189. [Google Scholar] [CrossRef]
  22. Al-Shaalan, A.M. Reliability Evaluation of Power Systems. In Reliability and Maintenance-An Overview of Cases; IntechOpen: London, UK, 2020. [Google Scholar]
  23. Niu, Y.-F.; Xu, X.-Z.; He, C.; Ding, D.; Liu, Z.-Z. Capacity Reliability Calculation and Sensitivity Analysis for a Stochastic Transport Network. IEEE Access 2020, 8, 133161–133169. [Google Scholar] [CrossRef]
  24. Satitsatian, S.; Kapur, K.C. An algorithm for lower reliability bounds of multistate two-terminal networks. IEEE Trans. Reliab. 2006, 55, 199–206. [Google Scholar] [CrossRef]
  25. Hassan, M.R.; Abdou, H. Multi-Source Multi-Sink Stochastic-Flow Networks Reliability under Time Constraints. Indian J. Sci. Technol. 2019, 12, 22. [Google Scholar] [CrossRef]
  26. Abed, M.H.; Tang, A.Y. Hybridizing Genetic Algorithm and Record-to-Record Travel Algorithm for Solving Uncapacitated Examination Timetabling Problem. electronic J. Comput. Sci. Inf. Technol. 2013, 4, 25–31. [Google Scholar]
  27. Sharma, M. Role and Working of Genetic Algorithm in Computer Science. Int. J. Comput. Appl. Inf. Technol. 2013, 2, 27–32. [Google Scholar]
  28. Malik, A. A Study of Genetic Algorithm and Crossover Techniques. Int. J. Comput. Sci. Mob. Comput. 2019, 8, 335–344. [Google Scholar]
  29. Mirjalili, S. Studies in Computational Intelligence. In Evolutionary Algorithms and Neural Networks; Springer: Berlin/Heidelberg, Germany, 2019; Volume 780. [Google Scholar]
  30. Lin, Y.-K. A simple algorithm for reliability evaluation of a stochastic-flow network with node failure. Comput. Oper. Res. 2001, 28, 1277–1285. [Google Scholar] [CrossRef]
  31. Zuo, M.J.; Tian, Z.; Huang, H.-Z. An efficient method for reliability evaluation of multistate networks given all minimal path vectors. IIE Trans. 2007, 39, 811–817. [Google Scholar] [CrossRef]
Figure 1. Network with three sources and two sinks.
Figure 1. Network with three sources and two sinks.
Mathematics 11 03902 g001
Figure 2. Network with two sources and two sinks.
Figure 2. Network with two sources and two sinks.
Mathematics 11 03902 g002
Figure 3. Network with two sources and three sinks.
Figure 3. Network with two sources and three sinks.
Mathematics 11 03902 g003
Table 1. Case1: The best solution found while varying NP2 and its corresponding M, S, and Rs.
Table 1. Case1: The best solution found while varying NP2 and its corresponding M, S, and Rs.
NP2M S R S
10[4242543354]360.973585
15[4434115554]360.968533
20[5333314545]360.968419
Table 2. Case1: F and its corresponding RX for different NP2 (10, 15, and 20).
Table 2. Case1: F and its corresponding RX for different NP2 (10, 15, and 20).
F R X
1000010111101011100110.8749
0101000111101011010110.8749
1001100111100010011110.8745
1000001111110111100010.8407
0011100111101011000110.8222
1101000110110111001010.7904
1001000110100011111110.7904
1101101011000110010110.7899
1101010110100011001110.7594
1101000110000111111010.7413
1101010111001110100010.9012
1100000111111111000010.9012
1000101111110111000010.8855
1011000111101110100010.8855
1001111111100110000010.8832
1011100111101110000010.8832
0011101111100111000010.8832
0101011111100111000010.8120
1011100111001110100010.8102
1001100110111111000010.8082
1001110111101110000010.7976
1001011111100110100010.7976
1101100111011110000010.7976
1100100111011111000010.7958
1101000111011110100010.7299
1011101111100110000010.8760
1100101111000111010010.8760
1100101111010111000010.8760
1100001111000111110010.8313
1101101111000110010010.8313
1111000110101111000010.8313
1101100111001110010010.8301
0111000111101111000010.8301
1010101111100111000010.8301
1011100111001110100010.8292
1101011111100110000010.8292
1101100111011110000010.8222
1101000111001110110010.8222
0111100111001111000010.8222
1101011110100111000010.8222
1110100111001111000010.8210
1110100111001111000010.8210
1101001110000111110010.8200
1100011111100111000010.7834
1011000111101110100010.7834
Table 3. Case2: The best solution found while varying NP2 and its corresponding M, S, and Rs.
Table 3. Case2: The best solution found while varying NP2 and its corresponding M, S, and Rs.
NP2M S R S
10[15164174189111499191014]1690.996641
15[819171119101552013138815]1810.989280
20[196181815711161319941717]1890.999984
Table 4. Case2: F and its corresponding RX for different NP2 (10, 15, and 20).
Table 4. Case2: F and its corresponding RX for different NP2 (10, 15, and 20).
F R X
20100220122122012211220.9916
12221122211020200212110.9645
11211111121022212012120.8583
20122121222002002211120.8530
21102210120121012222210.7601
21211021211211012202210.7536
22222012220110022110210.6973
02221022122021022210200.6722
21220211222010012222100.6591
02212021112012102202220.6070
10202222211211111002220.9234
11122112122110222002200.9210
01212121122111122001220.9206
10201222201221220002220.8997
11002220222220210102220.8997
20002221221221121100220.8997
20110212222220222000120.8996
10102222221221210100220.8943
21001022222222121000220.8370
12202211212100220211210.8333
22101112122111122010220.8331
22112020212002220201220.8331
11222120220120221200120.8330
10222110221210211201220.8263
11212212112100212210120.5717
10222112201111022221110.9992
11212002211112222012110.9957
21211101120022122122200.9957
20201122122202112111110.9941
11202102202201012222210.9906
22012122011021201221120.9584
22221001220021220221200.9566
22022022111120220022110.9526
11201112121221022120210.9522
11202222021210200220220.9522
21100012212222221110120.9229
21211102200220222112110.9176
20111102111122122221200.8803
02202222111021202112020.7910
10102112011222212222100.7149
11022012212022221002210.7109
21001122102122122111120.7105
12211122221120022010120.7098
21110122222022220200200.7095
20212022211202010222200.6072
Table 5. Case3: The best solution found while varying NP2 and its corresponding M, S, and Rs.
Table 5. Case3: The best solution found while varying NP2 and its corresponding M, S, and Rs.
NP2M S R S
10[8 12144 15 1372 141210]1110.995669
15[12 6 12145912 6 12149]1110.999923
20[15 12 12159 4 1011 1445]1110.999403
Table 6. Case3: F and its corresponding RX for different NP2 (10, 15, and 20).
Table 6. Case3: F and its corresponding RX for different NP2 (10, 15, and 20).
F R X
011101100110101011000101110.9595
111000111111011001000100110.9595
101110110111011001001101000.9573
101100110111100110010110010.9573
111101011011001001100111000.9573
101101101011110001011100010.9212
101101011001100101101111000.8851
101000001111110111011101000.8851
001101101011101101011101000.8851
110100011101101011001100110.8851
110010000111101111011101000.9999
110101110001101010010111100.9998
111011000011110001110110010.9962
001101100001111001111111000.9962
110000000101111011110110110.9962
011100111001101001011011010.9962
101110010101100110011111000.9962
101111000010011011010101110.9960
101111110000001100011101110.9039
101110010101100111001110010.9038
011000010111101110110111000.9038
101111000101001101100011110.8988
101101000010011011110111010.8988
110001010110101110101101010.8988
011011011000101101000111110.8007
011111010001101010011100110.9994
111011100001100101011001110.9994
100111010100111000100111110.9988
111010010001100111010110110.9988
101101100110110000110111010.9988
101110010010100111011101010.9988
111110010010100011010101110.9988
001100011110111001100101110.9982
100111010110101110000111010.9606
101111011011100100010011010.9600
111000110010100111011101010.9600
101100110011011011010101100.9600
100100011111101110011100010.9600
011000111101111000010110110.9600
101111011010011000011101010.9600
011011000101111010011001110.9227
110011001010111001010111010.9227
101100110101110011000111100.9227
101111010000100101101111010.9227
011011001110101101000101110.9227
Table 7. The number of NF found by [1,2,7].
Table 7. The number of NF found by [1,2,7].
ApproachProblem Maximum   No .   of   N F
[1]Maximizing reliability of the capacity vector by searching the single optimal flow vector with minimum assignment cost 1
[2]Optimal system reliability under transmission budget10
[7]Optimal system reliability evaluation under DD 9
Proposed
approach
Optimal system reliability under maximum capacities (Robust Design Problem)20
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Boubaker, S.; Radwan, N.H.; Refaat Hassan, M.; Alsubaei, F.S.; Younes, A.; Sennary, H.A. Robust Design Problem for Multi-Source Multi-Sink Flow Networks Based on Genetic Algorithm Approach. Mathematics 2023, 11, 3902. https://doi.org/10.3390/math11183902

AMA Style

Boubaker S, Radwan NH, Refaat Hassan M, Alsubaei FS, Younes A, Sennary HA. Robust Design Problem for Multi-Source Multi-Sink Flow Networks Based on Genetic Algorithm Approach. Mathematics. 2023; 11(18):3902. https://doi.org/10.3390/math11183902

Chicago/Turabian Style

Boubaker, Sahbi, Noha Hamdy Radwan, Moatamad Refaat Hassan, Faisal S. Alsubaei, Ahmed Younes, and Hameda A. Sennary. 2023. "Robust Design Problem for Multi-Source Multi-Sink Flow Networks Based on Genetic Algorithm Approach" Mathematics 11, no. 18: 3902. https://doi.org/10.3390/math11183902

APA Style

Boubaker, S., Radwan, N. H., Refaat Hassan, M., Alsubaei, F. S., Younes, A., & Sennary, H. A. (2023). Robust Design Problem for Multi-Source Multi-Sink Flow Networks Based on Genetic Algorithm Approach. Mathematics, 11(18), 3902. https://doi.org/10.3390/math11183902

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop