Fuzzy Multiobjective Reliability Optimization Problem of Industrial Systems Using Particle Swarm Optimization

Thepresentwork investigates the reliability optimization problemof the repairable industrial systems by utilizing uncertain, limited, and imprecise data. In many practical situations where reliability enhancement is involved, the decision making is complicated because of the presence of several mutually conflicting objectives. Moreover, data collected or available for the systems are vague, ambiguous, qualitative, and imprecise in nature due to various practical constraints and hence create some difficulties in optimizing the design problems. To handle these problems, this work presents an interactive method for solving the fuzzy multiobjective optimization decision-making problem, which can be used for the optimization decisionmaking of the reliability with two or more objectives. Based on the preference of the decision makers toward the objectives, fuzzy multi-objective optimization problem is converted into crisp optimization problem and then solved with evolutionary algorithm.The proposed approach has been applied to the decomposition unit of a urea fertilizer plant situated in the northern part of India producing 1500–2000 metric tons per day.


Introduction
Reliability in general can be defined as the ability of a system to perform its required functions under stated conditions for a specified period of time.Reliability technology is an important phenomenon and is widely used for increasing the efficiency, risk analysis, production availability studies, and design of industrial systems.The industrial systems run continuously and suffer failure over a period of time which can be brought back in to service by proper repair or replacement.Consequently, it may be extremely difficult, if it is not possible to construct accurate and complete mathematical model for the system in order to access the reliability because of inadequate knowledge about the basic failure events [1,2].In highly competitive industrial market, the concept of failure analysis is an unavoidable fact in complex industrial systems.Reliability of such systems not only depends on the reliability of each element of these systems but also depends on occurrence of sequence of failures.
Generally, the design problems are always stated in precise mathematical forms.It must be recognized that many practical problems encountered by designers and decision makers would take place in an environment in which the statements might be vague or imprecise.Usually, it is difficult to describe the goals and constraints of such optimization problems by crisp relations through equations and/or descriptions.In such situations, the traditional reliability theory, based on probabilistic and binary state assumptions, does not always provide useful information to the practitioners due to the limitation of being able to handle only quantitative information [3][4][5].An alternative for this fuzzy set theory [6] can build a model to represent a subjective estimation of possible effect of the given values on the problem and permit the incorporation of vagueness in the conventional set theory that can be used to deal with uncertainty quantitatively.Due to incomplete and uncertain input information, mathematical models of such problems are developed in fuzzy environment, and the optimization problem under consideration becomes a fuzzy programming problem.The fuzzy set-based optimization was firstly introduced by Bellman and Zadeh [7] in their seminal paper on decision making in a fuzzy environment, in which the concepts of fuzzy constraint, fuzzy objective, and fuzzy decision were introduced.After this pioneering work, these concepts were subsequently profusely used and applied by many investigators.In most of the practical design situation, the presence of several conflicting objectives which are nonlinear and ambiguous complicates the reliability apportionment.For instance, a designer is required to minimize the system cost while simultaneously maximizing the system reliability.Therefore, multiobjective functions become an important aspect in the reliability design of the engineering systems, and hence various researchers [8][9][10][11][12][13][14][15] used different methods to solve reliability optimization problem in fuzzy environment.
In multiobjective optimization problems (MOOPs), it is difficult or rarely possible to find an optimal solution for all the objectives which simultaneously optimize the problem in fuzzy environment.For handling such types of situations, one usually tries to search for a solution which is as close to the decision makers (DMs) expectations as possible.For this reason, problem is solved is an interactive manner in which DM is initially asked to specify his or her preferences towards the objectives.Based on these preferences, the problem is solved, and the DM is provided with a possible solution.If the DM is satisfied with this solution, the problem ends there; otherwise, he or she is asked to modify his or her preferences in the light of the earlier obtained results.This iterative procedure is continued till a satisfactory solution is achieved which is close to DM's expectations.
The present work is an extension of the work earlier done by Sharma and Garg [16] and Kumar [17], in which the cost factor was not considered in mathematical modeling.In this study, a conflicting multiobjective nonlinear programming problem is considered in fuzzy environment where we maximize the reliability and minimize the cost of the system.A conflicting nature between the objectives is resolved with the help of the fuzzy technique.Also, the intention is to use compensatory operator for aggregation of the different fuzzy goals and a robust global optimization technique, namely particle swarm optimization (PSO), for solution of the resultant single objective optimization problem thus formulated with the choice of the DM/system expert regarding the priority amongst the objectives.The technique is explained through a case study of decomposition unit of a fertilizer plant, a complex repairable industrial system.The rest of the chapter is organized as follows.Section 2 describes the general definitions related to multiobjective optimization problem, while corresponding reliability optimization model is given in Section 3. The solution procedure for solving the multiobjective reliability optimization problem is described in Section 4. A system description for illustrating the approach has been given in Section 5, while their corresponding results and analysis are done in Section 6.Finally, some concrete conclusions are drawn and are presented in Section 7.

A General Multiobjective Nonlinear Programming Problem
A multiobjective nonlinear programming problem is to find the solution vector or vector of design variables  = ( 1 ,  2 , . . .,   )  that optimizes a vector of objective functions over the feasible design space X.The problem is modeled as follows (WLOG we assume the minimization): Minimize:  () where  1 (),   In general, there exist a number of Pareto optimal solutions to a multiobjective optimization problems.Due to multiobjectives, the selection of such objectives clearly depends on the problem under study and the decision maker (DM) criteria.Thus, the designer must select a compromise or a satisfying solution from the Pareto optimal solution set according to his or her preference.

Reliability Optimization Model
Reliability is one of the vital attributes of performance in arriving at the optimal design of a system because it directly and significantly influences the system's performance.In practical, the problem of system reliability may be formed as a typical nonlinear programming problem with nonlinear costfunctions in fuzzy environment.

Formulation of System Reliability
Model.Let us consider that the reliability problem of a system consists of  components.Each component has reliability   for the ith components for  = 1, 2, . . ., .Then, the system reliability is written in the form of reliability of each component as for series system (1 −   ) for parallel system or combination of series and parallel system. ( According to Aggarwal and Gupta [18], the cost of reliability is monotonically increasing function of reliability and hence based on the fact that the ith components reliability cost is   (  ).Therefore, the system cost is given by In reliability optimization problems, it is often required to minimize or maximize several objectives subject to several constraints.For instance, a designer is required to minimize the system cost while simultaneously maximizing the system reliability.Therefore, multiobjective functions become an important aspect in the reliability design of the engineering systems.Hence, the suitable form of optimization model of series system reliability problem by considering the system reliability and cost as objective is Maximize: or combination of series and parallel system Minimize: for  = 1, 2, . . ., . (5)

System Reliability Optimization Model in Fuzzy Environment.
The necessary features of the cost versus maintainability function are equivalent to the cost versus reliability function as given by Aggarwal and Gupta [18].It is very complicated decision-making process to determine the reliability components in fuzzy objective as well as constraint goal.It involves many uncertain factors and becomes a nonstochastic vague decision-making process.Therefore, the reliability allocation model ( 5) can be represented by fuzzy nonlinear programming to make the model more flexible and adoptable to the human decision process.Therefore, in fuzzy environment, the optimization problem (5) becomes The symbol Max denotes a relaxed or fuzzy version of "Max."

Interactive Methods for Solving Multiobjective Reliability Optimization Problems
Multiobjective optimization is intensively used in engineering applications for simultaneously optimizing the collection of objective functions systematically.In general, reliability optimization problem is solved with the assumption that the coefficients or cost of components is specified in a precise way.In real life, there are many diverse situations due to uncertainty in judgments, lack of evidence, and so forth such that there are incompleteness and unreliability of input information, and hence it is not possible to get relevant precise data for the reliability system.These types of impreciseness in data are well handled with the help of defining their membership functions by fuzzy set theory instead of representing by random variable.Thus, due to these, the concerned optimization problems are modeled in fuzzy environment, and hence the problem under consideration becomes a fuzzy programming problem.But in a multiobjective optimization problem, an optimal solution which simultaneously optimizes all the objectives, and that too when the problem is modeled in a fuzzy environment, is rarely possible.In such situations, one usually tries to search for the best possible solution in the presence of imprecise information which is as close to the DM's expectations as possible.Search of such a satisfying solution requires solving the multiobjective fuzzy optimization problem iteratively in an interactive manner, wherein the DM is initially asked to specify his or her preferences and expectations.Based on these preferences, the problem is solved, and the DM is provided with a possible solution.If the DM is satisfied with this solution, the problem ends there; otherwise, he or she is asked to modify his or her preferences in the light of the earlier obtained results.This iterative procedure is continued till a satisfactory solution is achieved which is close to DM's expectations.The detail of the computational procedure for solving the fuzzy multiobjective reliability optimization problem in an iterative way is described in the following steps.
Step 1 (find the ideal and anti-ideal values of each objective function).
where  * 1 and  * 2 are the ideal solutions of the objective functions corresponding to minimization and maximization, respectively.
Step 2 (establishing the fuzzy goals towards the objective functions).In the traditional optimization, the design feasibility is considered as either satisfied or violated.For many engineering applications, the transition from infeasibility to feasibility is not obvious, because of not only the vague information in the design constraints but also the factors that can affect the design scenario, such as designer's knowledge, manufacture precision, and material properties.For this reason, they are modeled in such a way that the transition from infeasible state to feasible state is smooth and gradual with subjectivity.To incorporate the DM's vague idea about in which region of the objective the optimum should be, the degree of fuzziness is used.Essentially, the approach is to translate the functions to a common scale [0, 1] by the means of mathematical transformations, combine them using the geometric mean, and optimize the overall metric.In this paper, the problem is fuzzified with the help of linear membership functions.Let R and C be the fuzzy region of satisfaction of system reliability (  ) and system cost (  ), respectively, and let    and    be their corresponding membership functions.Then, the fuzzy objective stated by a designer can be quantified by eliciting corresponding linear membership functions, using the minimal and maximal feasible values of each objective as obtained during Step 1, and is defined as.
Step 3 (equivalent single optimization problem).Using the achieved objectives' membership functions of   and   , it is very important to choose the aggregation operator.Even since Zadeh [6] suggested that ∧, that is, min be used for the intersection of fuzzy sets as The two objectives    and    in (10) are equally important.This is not true in a real-life situation; that is, DMs sometimes do not pay equal attention to these two objectives.Owing to this (10), is modified, according to the importance of the objective, by Huang [11] as where  1 and  2 ∈ [0, 1] are called the objective weights.Thus, using the achieved objectives' membership functions and DM/system expert preferences in the form of weights problem is formulated as a single objective optimization problem and is given as where ∧ indicates the intersection or min operator,   represents the tth objective weight, in the form of the importance of the objective functions, suggested by DM,  is the vector of decision variables, and    and    are the lower and upper bounds of decision vector   , respectively.The obtained optimization problem is solved with the particle swarm optimization algorithm which has been described in Section 4.1.
Step 4 (adjusting the preference parameters).If the DM is satisfied by the solution obtained in Step 3, then the approach stops successfully.Otherwise, the key preference parameters, that is, decision maker's preferences regarding the relative importance of each objective function ( = [ 1 ,  2 ]), can be altered to meet the DM's choice, and the method again goes back to Step 3. The process is repeated until DM is satisfied.We are just showing one run of the approach here as we assume that in this problem DM is satisfied by the results obtained in Step 3.

4.1.
A Survey of PSO Algorithm.The particle swarm optimization (PSO) algorithm was firstly proposed by Kennedy and Eberhart [19,20] and has deserved some attention during the last years in the global optimization field.PSO algorithm is based on the population of agents or particles and tries to simulate its social behavior in optimal exploration of problem space.The PSO algorithm is inspired by social behavior of bird flocking, animal hording, or fish schooling.In PSO, the potential solutions, called particles, fly through the problem space by following the current optimum particles.PSO is initialized with a group of random particles (solutions) and then searches for optima by updating generations.During every iteration, each particle is updated by following two "best" values.The first one is the position vector of the best solution (fitness), this particle has achieved so far.The fitness value is also stored.This position is called pbest.Another "best" position that is tracked by the particle swarm optimizer is the best position, obtained so far by any particle in the population.This best position is the current global best and is called gbest.After finding the two best values, the particle updates its velocity and position according to ( 13) and ( 14), respectively, where V   is the velocity of ith particle at the kth iteration and    is currently the solution (or position) of the ith particle.ud and Ud are random numbers generated uniformly between 0 and 1.  1 is the self-confidence (cognitive) factor, and  2 is the swarm confidence (social) factor.Finally,  is the inertia factor that takes linearly decreasing values downward according to a predefined number of iterations.The first term in (13) represents the effect of the inertia of the particle, the second term represents the particle memory influence, and the third term represents the swarm (society) influence.The algorithm for the PSO can be summarized as follows.
(1) Initialize the swarm   , the position of particles is randomly initialized within the hypercube of feasible space.
(3) Evaluate the performance  of each particle, using its current position   .
(6) Update the velocity of the particle according to (13).
(7) Move each particle to a new position using (14).
(9) Go to Step 3, and repeat until convergence.

Illustrative Example
To illustrate, a fertilizer plant situated in the northern part of India and producing approximately 1500-2000 metric tons per day has been considered as a main system [17].The brief description of the system is given below.

System Description.
The fertilizer plants are large, complex, and repairable engineering unit which is a combination of ammonia and urea plant.The urea plant is composed of synthesis, decomposition, crystallization, and prilling system, arranged in predetermined configuration [16,17].In this process, the ammonia and CO 2 enter the urea synthesis reactor.The reactants from urea synthesis reactor enter the urea decomposer in which urea is separated from reactants.These are further sent to urea crystallizer in which the urea solution is concentrated and crystallized.The urea crystals are separated by centrifuge and conveyed pneumatically to the top of urea prilling system.In this system, urea crystals are melted and sprayed through distributers and fall down in urea prilling tower against the ascending air allowing getting prilled on the way.The prilled urea is collected at the bottom of urea prilling, system and sent to bagging section.Among the various functional units in the plant such as urea synthesis, urea decomposition, urea crystallization, urea prilling and urea recovery, urea decomposition is one of the most important and vital functional processes which is the subject of our discussion.

Decomposition System.
The gas-liquid mixture (urea, NH 3 , CO 2 , Biuret) flows from a reactor at 126 ∘ C into the upper part of a high-pressure decomposer, where the flashed gases are separated.The liquid falls through a sieve plate, which comes into contact with high-temperature gas available from the reboiler and the falling film heater.The same process is repeated in a low-pressure absorber.In the reboiler, the liquid is further heated to 151 ∘ C with medium-pressure steam, so that the remaining ammonia and carbonate are released as gases.The solution is then further heated to 165 ∘ C in the falling film heater, which reduces the Biuret formation and hydrolysis of urea.The overhead gases from the highpressure decomposer go to the high-pressure absorber cooler.The liquid flows to the top of the low-pressure absorber and is cooled in a heat exchanger.Additional flashing of the solution takes place in the upper part of the low-pressure absorber to reduce the solution pressure from 17.5 to 2.5 kg/cm 2 .The low-pressure absorber has four sieve trays and a packed bed.
In the packed bed, the remaining ammonia is stripped off by CO 2 gas.The overhead gases go to the low-pressure absorber cooler, in which the pressure is controlled at 2.2 kg/cm 2 .Most of the excess ammonia and carbonate are separated from the solution flowing to the gas separator.The gas separator has two parts: (a) the upper part is at 105 ∘ C and 0.3 kg/cm 2 ; here the remaining small amounts of ammonia and CO 2 are recovered by reducing the pressure; the sensible heat of solution is enough to vaporize these gases; (b) the lower part has a packed section at 110 ∘ C and atmospheric pressure.Air containing a small amount of ammonia and CO 2 is fed off from the gas absorber by an off-gas blower to remove the remaining small amounts of ammonia and CO 2 present on the solution.Off-gases from the lower and upper parts are mixed and fed to the off-gas condenser.The urea solution, concentrated to 70-75%, is fed to a crystallizer.In brief, the various subsystems and the components associated with them are defined as below.(i) Subsystem   has two units.Unit  1 is called the reboiler for the high-pressure absorber, and unit  2 is called the falling film heater for the low-pressure absorber.Failure of  1 or  2 causes complete failure of the system.(ii) Subsystem   has two units in series.Unit  1 is called high-pressure absorber, and  2 is called lowpressure absorber (the component contained of sieve trays and packed bed for tripping off the remaining ammonia).Failure of either unit causes the failure of whole system.(iii) Subsystem , the gas separator, has one unit only (used to separate the gases obtained from pressure absorbers.The solution is fed to crystallization unit for further processing) arranged in series with  1 and  2 .
(iv) Subsystem   has two units in series that is  1 and  2 where  1 is low-pressure heat exchanger, and  2 is high-pressure heat exchanger with standby unit (the heat exchangers are used to recover the heat of the gases).Failure of both at a time will cause complete failure of the system.
The schematic diagram of the system is shown in Figure 1.
The data related to main components of the system, in the form of failure rate (  's) and repair time (  's), is collected from the historical or available records of the industry and is integrated with expertise of maintenance personnel which is given in Table 1 [16].

Reliability Optimization Problem.
The multiobjective reliability optimization problem for the considered system is formulated by taking systems reliability and cost as an objective and crisp failure rates (  's) and repair times (  's) as decision variables under the considered ±15% uncertainty level or support towards the data as Maximize   = exp (−  ) .

Computation Results
This section turns to the description and analysis of the results obtained by the optimization tests.
6.1.Parameter Setting.In all algorithms, the values of the common parameters such as population size and total evaluation number are chosen to be the same.Population size and the maximum evaluation number are taken as 20 × , where  is dimension of the problem and 1500, respectively, for the function.The method has been implemented in MATLAB (MathWorks); and in order to eliminate stochastic discrepancy, 30 independent runs have been made that involve 30 different initial trial solutions.The termination criterion has been set either limited to a maximum number of generations or to the order of relative error equal to 10 −6 , whichever is achieved first.The other randomly specified parameters of algorithms are given below.
6.1.1.GA Settings.In our experiment, real-coded genetic algorithm is utilized to find optimal values.Roulette wheel selection criterion is employed to choose better-fitted chromosomes.One-point crossover with the rate of 0.9 and random point mutation with the rate of 0.01 are used in the present analysis for the reproduction of new solution.
6.1.2.PSO Settings.Except common parameters (population number and maximum evaluation number), cognitive ( 1 ) and social ( 2 ) components are constants that can be used to change the weighting between personal and population experience, respectively.In our experiments, cognitive and the social components were both set to 1.49.Inertia weight (), which determines how the previous velocity of the particle influences the velocity in the next iteration, was defined as the linearly decreases from initial weight  1 = 0.9 to final weight  2 = 0.4 with the relation  =  2 + (iter max − iter)( 1 −  2 )/iter max where iter max is the maximum number of iteration, and iter is used iteration number [21].

Results and Discussions.
The ideal values of the objective functions are given as    = 0.98046407,    = 0.98554720,    = 1039.731722,and    = 1065.984110.Based on these ideal values, a fuzzy region of satisfaction corresponding to the objective functions is to be constructed as follows: Using these constructed membership functions and their weight vector as suggested by DM/system experts corresponding to the two objective functions, the equivalent crisp optimization problem is formulated as follows: PSO has been used to solve (18) with the initial preference of the weight towards the objectives as 1 = [1,1].Since in the first iteration DM does not want the preference towards the objectives; that is, they pay equal attention towards each objective.Hence, the weight vector is taken as 1 for each objective.Results corresponding to the first iteration are shown in Table 2, which is (  ,   ) = (0.98359477, 1049.92107) with membership value (   ,    ) = (61.59%,61.89%).Not satisfied with this outcome or willing to know other possible optimal solutions, keeping this result in view, let DM be decided to give more importance on reliability objective than cost objective and give a preference weight vector 1 corresponding to reliability and 0.5 to cost; that is, 2 = [1, 0.5].In other words, DM wants to pay two times more attention towards the reliability objective than cost objective.Outcome of this iteration is (  ,   ) = (0.98451928, 1055.59802) with membership value (   ,    ) = (79.77%,39.56%).Process is repeated till DM is fully satisfied.In this way, DM obtained the different possible solution for different satisfaction levels.Process is stopped after V iteration (it may continue further till DM is satisfied).Outcome of iteration V is (  ,   ) = (0.98255695, 1044.58702) with membership value (   ,    ) = (41.17%,81.50%).This result shows that 41.17% achievement for reliability and 81.50% for cost of respective fuzzy goals.A comparison of the results listed in the table with the results obtained by GA.It has been seen from the table that the results computed by PSO are better in terms of preferences given by the DM towards the objectives as compared to GA results.Thus, for different preferences suggested by DMs, optimum values of systems' reliability and cost are achieved.The optimum design parameters of design variables corresponding to optimum values are also summarized in Table 2. Based on these decision variables, the system analyst or decision makers may plan the schedule for proper maintenance in regular interval of time.

Conclusion
This paper reports the multiobjective optimization problem of the urea decomposition unit of a fertilizer plant under the fuzzy environment where we maximize the reliability and minimize the cost of the system.A proposed approach suggests a way for achieving the optimum performance of the system with the help of experts/DMs elicitation.A mutual conflicting nature of the objectives is resolved with the help of fuzzy set theory after constructing the fuzzy region of satisfaction by taking linear membership functions.Since reliability decision is usually made in the earliest stage of system design and the information at this stage is incomplete and imprecise, it is necessary to rely on the experience of DMs and experts.Based on their preferences in the form of weights towards the objectives along with the achieved objective membership functions, FMOOP is reformulated as a single objective optimization problem and then solves with the PSO and GA iteratively.The one run of the approach has been shown here as we assume that in the problem DM is satisfied by the obtained results.Through this approach, a decision support system has been developed which helps the plant maintenance personnel in deciding his future strategy to gain optimum performance of the system.The decision variable corresponding to the main components of each subsystem/unit of a plant is reported which may be targeted so that optimum system performance could be achieved by using the discussed approach.
The present work done may be extended for a system in hybrid configuration, for example, series, parallel, seriesparallel, and so forth.The presented methodology will be further extended and improved using other optimization tools/algorithms and artificial neural network will be used to handle the complex nature of the systems.As a constant, failure rate model has been taken in the analysis, so authors are engaged in the approach for further analysis such as resource allocation, facility planning and management, inventory control, network analysis, and job shop scheduling for an arbitrary failure rate model instead of constant failure rate model.

Table 2 :
Solution of the reliability optimization problem of decomposition unit.