- Research
- Open Access

# A framework for end-of-life photovoltaics distribution routing optimization

- Qi Guo
^{1}Email author and - Hongjun Guo
^{2}

**29**:3

https://doi.org/10.1186/s42834-019-0005-8

© The Author(s) 2019

**Received:**8 August 2018**Accepted:**26 November 2018**Published:**19 March 2019

## Abstract

Transportation cost accounts for a significant portion among the total cost of photovoltaics (PV) recycling. The transportation cost can be significantly reduced with a well-planned vehicle routing. A generic mathematical framework was developed to generate the optimal distribution scheme for transporting retired PV. A static and a dynamic recycling algorithms were created as the optimization tool in the framework. As for the output, the framework generates the optimal distribution scheme, the transportation cost of the optimal scheme, as well as the amount of End-of-life PV that each PV recycling center receives. A case study was included to test the effectiveness of the proposed framework. The optimal costs generated by two algorithms were compared with the baseline cost. As a result, the cost can be reduced by 53% when compared with the baseline. In addition, the static algorithm can obtain a result with decent accuracy and low computational cost.

## Keywords

- Photovoltaic recycling
- Mathematical modeling
- Route optimization

## Introduction

In 2015, 50.6 GW new solar photovoltaics (PV) were installed worldwide. A 613 GW grid-connected solar power was forecasted by the end of 2019 all over the world [1]. The amount of PV installed in the United States in 2016 achieved 14.7 GW, doubling the amount of the year 2015. With the 39% renewable electric occupancy, solar energy has become the most used renewable electric energy in the United States in 2016. As the largest installer, in 2016, PV installation in California accounted for 35% of the United States. The cumulative PV installation in CA has researched 17 GW [2].

Due to containing toxic and value materials, all on site PVs need to be recycled after they retire in about 25 yr. [3]. However, recycling is less economical than disposing at present. Cost effectiveness has been emphasized as the greatest challenge while performing PV recycling [4].

Reducing the transportation cost is an effective way to reduce the cost of PV recycling. The transportation cost has been concluded as the leading cost while simulating the PV recycling process in Germany [5]. While planning EoL PV material recovery in New York, the transportation cost was also concluded as an indispensable parameter in the recycling framework [6]. The proportion of the transportation cost among the total cost was also emphasized while performing a financial analysis of *end-of-life* (EoL) PV recycling in Italy [7]. Therefore, in order to make the recycling affordable and profitable, the transportation cost needs to be reduced [8].

Implementing mathematical optimization is an effective way to reduce the transportation cost. Classical mathematical optimization models of the vehicle routing for the facility location decision has been summarized [9]. A waste collection optimization was performed by minimizing the route distance. By optimizing, the average waste collection efficiency was improved by 37%, and the fuel cost was reduced by 48% [10]. An optimization on solid waste collection and transport was performed by using the ArcGIS Network Analyst tool. As a result, a 48% gas saving was achieved on the optimal transport scenario [11].

As for the optimization algorithm, Mixed Integer Linear Programming (MILP) and Greedy Algorithm (GA) have been used broadly on mathematical models. In order to solve a EoL vehicle recovery problem in Turkey economically, a MILP optimization model was developed to minimize the costs on opening facilities, recovery processes, and transportation [12]. MILP was also applied on vehicle routing problems to maximize the on-time possibility with the shortest path. Some advanced GAs were developed on top of the original GA to solve specific problems. A modified GA was employed to solve the stochastic control problem on the networked storage operation area [13].

Though many mathematical optimization researches have been done, there is not an algorithm that can be directly used to minimize the transportation cost in PV recycling. The paper focuses on developing a generic mathematical optimization framework to facilitate EoL PV recycling to reduce the transportation cost. Two recycling algorithms were proposed in the framework. A case study was included to show the efficacy of the framework.

## Materials and methods

A mathematical modeling framework was developed to generate the optimal distribution scheme by minimizing the transportation cost. In the framework, three parameters are included: 1) Photovoltaics Installation Site (PVIS); 2) Photovoltaics Recycling Center (PVRC); and 3) Transportation Company (TC). This is typically a *combinatorial optimization problem*: each PVIS can be taken by any of TC, and sent to any PVRC. As a result, the optimal distribution scheme is presented by showing which PVIS is transported by which TC, and goes to which PVRC. The maximum annual capacities of each PVRC and TC are constraints to affect the distribution decision. Besides, the total transportation cost of the optimal distribution scheme can be calculated.

Three types of input data are required in the framework. These data are the transportation cost charged by TC, the distance between each PVIS, each PVRC, and each TC, and the size of each PVIS and the annual capacity of each PVRC and TC. First, in the cost charge aspect, the total transportation cost consists of the overhead charge, the mileage charge, and the weight charge. The overhead charge ($/trip) is a fixed expense charged by each trip. The mileage charge is the product of the unit mileage charge ($/km) and the mileage (km) travelled. The weight charge is the product of the unit weight cost ($/kg) and the weight of the material (kg) transported. The total cost of each transportation service is the sum of the overhead charge, the mileage charge, and the weight charge. The second input is the distance matrices consisting of a set of distance combination among PVIS, TC, and PVRC. With the latitude and the longitude of each PVIS, TC, and PVRC, the distance combination between two of three locations can be calculated and expressed in matrix forms. The third input is the size of each PVIS, TC, and PVRC. The weight of material of each PVIS affects the weight cost. The annual capacity of each TC and PVRC as a physical constraint affects the distribution result.

Two types of optimization algorithms were used in the modeling framework: a static optimization algorithm, and a dynamic optimization algorithm.

### Static optimization algorithm

The static optimization algorithm is a type of the GA. GA is an algorithm of finding the optimal solution by using a heuristic, which orders inputs by a seemingly logical way. In theory, as long as all inputs are ordered in the special way, the GA makes the same decisions that the optimal solution does. In practice, the optimal solution obtained by the GA is most likely a local optimum. However, the local optimum solution is sometimes very close to the global optimum with the reduction of significant amount of running time.

In the model, all PVIS can be sorted by their sizes from the largest to the smallest. The explanation is that the larger PVIS are the ones that would save the most money by choosing the cheapest transportation option. After sorting the size of PVIS, the largest PVIS has the first priority to be recycled by the first TC that offers the cheapest price. The less size PVIS has, the later that PVIS will be considered by the GA. As a result, for PVIS with larger sizes, the cost on transportation of those is minimal. Though the saving of each PVIS is small, the total saving can accumulate by considering some thousands of PVIS. The annual capacities of PVRC and TC were included in the model as constraints.

### Dynamic optimization algorithm

The dynamic optimization algorithm solves the problem by using the MILP. Eq. (1) is the objective function to minimize the total transportation cost. Again, the total transportation cost is the sum of the overhead charge, the mileage charge, and the weight charge. Each type of the charge was described in the previous section. Eq. (2) is to ensure all installed PVIS are sent as a whole to one of the PVRC. The linear inequalities in Eq. (3) restrict that the total weight that each PVRC receives is less than its maximum capacity. The linear inequalities in Eq. (4) ensure that the total weight of each TC transports is less than its maximum transportation capacity. The model iterates until it identifies the minimum total transportation cost.

*i*,

*j*, and

*k*represent PVIS, PVRC, and TC, respectively.

*o*

_{k},

*m*

_{k}, and

*w*

_{k}is the overhead charge, the unit mileage charge, and the unit weight charge that TC

*k*offeres, respectively.

*n*

_{i}is the number of required trips to transport PVIS

*i*.

*D*

_{ki}

*D*

_{ij}, and

*D*

_{jk}are the transportation distances between corresponding locations.

*W*

_{i}is the weight of material of PVIS

*i*.

*Y*

_{ijk}is a binary decision factor: if PVIS

*i*is transported to PVRC

*j*by TC

*k*, then

*Y*

_{ijk}= 1, otherwise,

*Y*

_{ijk}= 0.

*C*

_{j}and

*C*

_{k}is the maximum annual recycling capacity of PVRC

*j*and TC

*k*, respectively.

### Experiment

A case study in California is considered to show the efficacy of the proposed framework. The State of California was selected because of its relatively long history and the significant amount of PV installations compared to other states. In the case study, three recycling periods and several PVRC locations were assumed based on preliminary studies. In those studies, three recycling periods were recommended for recycling the EoL PV between the years 1984 and 2015 based on the amount of EoL PV as well as the increasing trend in each year. The location as well as the annual capacity of PVRC in each period was recommended as well: two-2.6 kt annual capacity PVRC in Period I, and addition 6–8 kt in Period II and Period III, respectively. The above assumptions were utilized for implementing the optimization model.

Considering the gigantic amount of PVIS in California, the timing of collection should be considered while operating and managing the recycling process in reality. The framework assumed all EoL PVIS have been disassembled and are ready to load on the truck. In addition, the following assumptions are made: the maximum load of a truck is 7.5 t, and the approximate weight of material of every 1 MW capacity PV modules is equivalent to 75 t [14]. With the above assumptions, each truck can take at most 100 kW PV modules. Considering the significant percentage (about 40%, in Fig. 1) of PV modules that are over 100 kW, the number of trips cannot be ignored.

Further reasonable assumptions are made to simplify the optimization process. Due to the gigantic number of PVIS (up to 160,000 PVIS per year) together with the combinational delivering options between PVIS, PVRC, and TC, reasonable simplifications need to make to satisfy the computational requirement on both software and hardware. Thus, another two reasonable assumptions are made. The first is the unit cost charge similarity assumption. In reality, the cost charge by each TC tends to be similar due to the competitive market. Therefore, the unit cost charge of each TC, such as $/trip, $/km, and $/kg, are assumed to be the same. The second is the dead mileage charge exemption assumption. In fact, TC usually do not charge on dead mileage. The dead mileage is the mileage generation when vehicles are driven on the way to or back from the worksite without taking the real workload. In the case study, the dead mileage is the mileage when trucks driving from TC to PVIS, and the mileage when trucks driving from PVRC to TC. Those dead mileage charges are reasonable to be exempted in the case study by considering the reality. With the cost charge similarity assumption and the dead mileage charge exemption assumption, TC is no longer a variable in the optimization process. Subsequently, the transportation cost function can be simplified from Eq. (1) to Eq. (5). The simplified transportation cost function consists of the overhead charge, *o*, the unit mileage charge *m*, the unit weight charge, *w*, the travel distance between the PVIS and the PVRC, *D*, the weight of material of the PVIS, *W*, and the number of trips, *n*. By using the simplified cost function, variables in optimization processes can be reduced by at least 50%.

With the simplification, the dynamic optimization algorithm is simplified below.

The simplified form can improve the feasibility for handling large numbers of PVIS and fasten the computational speed of MILP. In addition, rather than MILP, Linear Program (LP) is recommended and utilized for solving the particular California case. It is concluded that LP is able to handle more variables, compute and converge much faster, and conduct the total transportation cost with the error of less than 1% of the total transportation cost conducted by the MILP.

While performing the optimization, reasonable unit charges were assumed. The real market price of unit charges was not used as a main focus in the paper because those prices vary from time by time. All the assumed prices were for initiating the process and showing the efficacy of the PV recycling framework. Thus, the following unit costs were utilized: the overhead cost is $10/trip, the unit mileage cost is $0.06/km, and the unit weight cost is $0.008/kg. The framework can be employed in any multi-facility recycling problem anywhere as long as the unit cost rates are available.

## Results and discussion

Total cost and multiple of cost of PVIS size and number of PVIS combination

PVIS size (kW) | Number of trips | Overhead cost ($) | Mileage cost ($) | Weight cost ($) | Grand total ($) | Number of PVIS | Total cost ($) | Cost multiplier |
---|---|---|---|---|---|---|---|---|

1 | 1 | 10 | 6 | 1 | 17 | 100 | 1660 | 22x |

10 | 1 | 10 | 6 | 6 | 22 | 10 | 220 | 3x |

100 | 1 | 10 | 6 | 60 | 76 | 1 | 76 | 1 |

3 | 1 | 10 | 6 | 2 | 18 | 100 | 1780 | 8x |

30 | 1 | 10 | 6 | 18 | 34 | 10 | 340 | 1x |

300 | 3 | 30 | 18 | 180 | 228 | 1 | 228 | 1 |

5 | 1 | 10 | 6 | 3 | 19 | 100 | 1900 | 5x |

50 | 1 | 10 | 6 | 30 | 46 | 10 | 460 | 1x |

500 | 5 | 50 | 30 | 300 | 380 | 1 | 380 | 1 |

Saving of GA’ and LP of each year

Period | PVRC recycling capacity (MW/yr) | Installed year | Amount of installed PV (MW) | GA’ saving (10 | LP saving (10 | Computing time ratio (LP/GA’) |
---|---|---|---|---|---|---|

I | 70 | 84–01 | 19 | 0 | 0 | 11 |

2002 | 29 | 0 | 0 | 15 | ||

2003 | 40 | 0 | 0 | 14 | ||

2004 | 63 | 20 | 20 | 15 | ||

2005 | 65 | 10 | 10 | 15 | ||

II | 700 | 2006 | 113 | 0 | 0 | 70 |

2007 | 287 | 2 | 2 | 152 | ||

2008 | 346 | 5 | 5 | 148 | ||

2009 | 418 | 8 | 9 | 235 | ||

2010 | 603 | 44 | 45 | 311 | ||

III | 1300 | 2011 | 643 | 2 | 2 | 869 |

2012 | 953 | 26 | 27 | 1187 | ||

2013 | 671 | 3 | 3 | 1506 | ||

2014 | 866 | 20 | 22 | 2544 | ||

2015 | 1216 | 473 | 491 | 3731 |

From Table 2, though the saving is as high as $491,000 in the year 2015, the percentage saving is only 15%, which is not significant. The reason of the insignificance results from the category of the charges in the current cost function. The current cost function is shown in Eq. (5), and is denoted as Cost Structure 1. Cost Structure 1 is composed of the overhead charge, the weight cost, and the mileage charge. In Cost Structure 1, the overhead charge ($/trip) and the weight cost are not optimized and remain the same based on the unit cost charge similarity assumption. The weight cost is identical because the unit weight cost ($/kg) is identical from the unit cost charge similarity assumption. Therefore, the mileage charge is the only parameter that determines the total cost saving. The next step is creating another cost structure, Cost Structure 2, which is only in terms of the mileage charge as shown in Eq. (9). The transportation cost in Cost Structure 2 only consists of the unit mileage cost, *m*, the travel distance, *D*, and the number of trips, *n*. Assume the unit mileage cost is $0.2/km.

Percentage saving of in each cost structure of each year. (Unit: %)

Period | Installed year | Cost structure 1 | Cost structure 2 | ||
---|---|---|---|---|---|

GA’ saving | LP saving | GA’ saving | LP saving | ||

I | 84–01 | 0 | 0 | 0 | 0 |

2002 | 0 | 0 | 0 | 0 | |

2003 | 0 | 0 | 0 | 0 | |

2004 | 9 | 9 | 24 | 25 | |

2005 | 5 | 6 | 16 | 16 | |

II | 2006 | 0 | 0 | 0 | 0 |

2007 | 0 | 0 | 2 | 3 | |

2008 | 1 | 1 | 6 | 6 | |

2009 | 1 | 1 | 7 | 8 | |

2010 | 4 | 5 | 25 | 26 | |

III | 2011 | 0 | 0 | 1 | 1 |

2012 | 2 | 2 | 15 | 16 | |

2013 | 0 | 0 | 2 | 2 | |

2014 | 1 | 1 | 8 | 8 | |

2015 | 14 | 15 | 51 | 53 |

## Conclusions

A generic mathematical framework was constructed to generate the optimal distribution scheme for transporting EoL PVIS. A static algorithm (GA) and a dynamic optimization algorithm (LP) were employed in the framework. The following conclusions were made from the simulation result based on the PV installed case in CA. First, based on the size distribution of PV installed in CA, the focus of the optimization should be on small PVIS rather than large ones. The most saving on transportation cost is obtained by optimizing transport distribution on small PVIS. Second, GA is a reliable algorithm to perform the optimization on transportation cost with decent accuracy, low computational cost, and no optimization solver required for the CA case. Third, the transportation cost saving increases with the annual capacity of PVRC getting tighter. With the appropriate optimization process and the specific cost structure applied, the transportation cost savings can be expected at 53% compared with not doing optimization. When the annual recycling capacity of PVRC is twice or larger than the amount of the recycled PV demand, there will not be transportation savings, and thus no optimization is needed. The developed framework has a wide applicability and can be applied on any PV recycling case for reducing the transportation cost.

## Declarations

### Authors’ contributions

Both authors read and approved the final manuscript.

### Competing interests

The authors declare that they have no competing interests.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- SPE. Global Market Outlook for Solar Power 2016-2020. Brussels: SolarPower Europe; 2016.Google Scholar
- Perea A, Honeyman C, Kann S, Mond A, Shiao MJ, Jones J, et al. U.S. Solar Market Insight: 2016 Year in Review. Edinburgh and Washington, DC: Wood Mackenzie, Limited and Solar Energy Industries Association; 2017.Google Scholar
- Choi JK, Fthenakis V. Design and optimization of photovoltaics recycling infrastructure. Environ Sci Technol. 2010;44:8678–83.View ArticleGoogle Scholar
- Hosenuzzaman M, Rahim NA, Selvaraj J, Hasanuzzaman M, Malek ABMA, Nahar A. Global prospects, progress, policies, and environmental impact of solar photovoltaic power generation. Renew Sust Energ Rev. 2015;41:284–97.View ArticleGoogle Scholar
- Choi JK, Fthenakis V. Crystalline silicon photovoltaic recycling planning: macro and micro perspectives. J Clean Prod. 2014;66:443–9.View ArticleGoogle Scholar
- Goe M, Gaustad G, Tomaszewski B. System tradeoffs in siting a solar photovoltaic material recovery infrastructure. J Environ Manag. 2015;160:154–66.View ArticleGoogle Scholar
- Cucchiella F, D'Adamo I, Rosa P. End-of-life of used photovoltaic modules: a financial analysis. Renew Sust Energ Rev. 2015;47:552–61.View ArticleGoogle Scholar
- Göllei A, Görbe P, Magyar A. Modeling and optimization of electrical vehicle batteries in complex clean energy systems. J Clean Prod. 2012;34:138–45.View ArticleGoogle Scholar
- Daskin MS, Snyder LV, Berger RT. Facility location in supply chain design. In: Langevin A, Riopel D, editors. Logistics systems: design and optimization. Boston: Springer; 2005. p. 39–65.View ArticleGoogle Scholar
- Akhtar M, Hannan MA, Begum RA, Basri H, Scavino E. Backtracking search algorithm in CVRP models for efficient solid waste collection and route optimization. Waste Manag. 2017;61:117–28.View ArticleGoogle Scholar
- Kallel A, Serbaji MM, Zairi M. Using GIS-based tools for the optimization of solid waste collection and transport: case study of Sfax City, Tunisia. J Eng. 2016;2016:1–7.View ArticleGoogle Scholar
- Demirel E, Demirel N, Gökçen H. A mixed integer linear programming model to optimize reverse logistics activities of end-of-life vehicles in Turkey. J Clean Prod. 2016;112:2101–13.View ArticleGoogle Scholar
- Qin JJ, Chow YL, Yang JY, Rajagopal R. Distributed online modified greedy algorithm for networked storage operation under uncertainty. IEEE T Smart Grid. 2016;7:1106–18.Google Scholar
- Choi JK, Fthenakis V. Economic feasibility of recycling photovoltaic modules. J Ind Ecol. 2010;14:947–64.View ArticleGoogle Scholar