Location Decision Making and Transportation Route Planning Considering Fuel Consumption

: This study presents the Location Routing Problem (LRP) for which we have created a model for the integration of locating facilities and vehicle routing decisions to solve the problem. The case study is the Palm Oil Collection Center, which is also important for the supply chain system. A mathematical model was made to minimize the total cost of a facility-opening cost, ﬁxed cost of vehicle uses and fuel consumption cost. The fuel consumption cost relies on the distance and road conditions, in case of poor physical condition of a road, and its width, which can be a ﬀ ected the speed of the vehicle as well as the used fuel. Thus, we propose an Adaptive Large Neighborhood Search (ALNS) based on heuristic for solving the LRP. The ALNS method was tested with three datasets of samples divided into small, medium and large problems. Then, the results were compared with the results from the exact method by the Lingo program. The computational study indicated that the ALNS algorithm was competitive to the results of the Lingo for all instance sizes. Moreover, the ALNS was more e ﬀ ective than the exact method; approximately 99% in terms of processing time. We extended this approach to solve the case study, which was considered to be the largest problem, and the ALNS algorithm was e ﬃ cient with acceptable solutions and short processing time. Therefore, the proposed method provided an e ﬀ ective solution to manage location routing decision of the palm oil collection center. bounds for medium-sized and large-sized. ALNS on the other hand, was able to find an acceptable solution within a short processing time and it was faster than Lingo program. Thus, we could summarize that ALNS was the optimal method for solving large problems. We have also implemented ALNS of the case study problem, which was larger than the instances. The ALNS produced the solution in an average processing time of 20.46 min, while the Lingo consumed over 120 hours, but was unable to find the solution. Thus, ALNS was efficient to solve the case study problem. Furthermore, we extended a comparison between the model of fuel consumption and distance minimization with respect to fuel cost. The result showed that this model could reduce 7.72% of fuel costs as compared to the classical LRP, which is distance minimization.


Introduction
Logistics management and operations cost management are the major driver of economic growth in many countries. According to the World Bank in 2016, the logistic capability index (LPI) of the G20 countries is at the highest level. Germany is the country that has the highest level of LPI at 4.23. As for Southeast Asian countries, Singapore is ranked 5 th , followed by Malaysia ranked 32 nd with an LPI score at 4.14 and 3.43, respectively [1]. Thailand is ranked 45 th with an LPI score of 3.26, which shows that the logistics performance of Thailand is at a low level. Therefore, it is influenced by the development of logistics in various sectors, especially in the agricultural sector, which is considered the main sector of Thailand.
Road transport is the preferred mode of transport, and is thus widely used for transportation of agricultural products in Thailand. This mode of transport relies on energy from fuel, resulting in high transportation costs. In 2017, total fuel consumption of road transport increased by 0.37%; compared to the previous year in terms of diesel and gasoline usage, this rate increased by 2.68% and 3.43%, effects of greenhouse gas (GHG) emissions from the transportation section in Korea from 1990 to 2013. The results indicate that economic growth, population growth and the transportation sector are the factors that can cause GHG emissions expansion. This study recommended some strategies for GHG deduction to achieve the national goal for transportation sector. Many countries in Asia have a number of REDD (Reducing Emissions from Deforestation and Degradation) projects that focus on reducing emissions from deforestation and forest degradation all over the world. Even though the benefits of REDD projects are clearly interesting, and environmental responsibility has become an important aspect of development both at the micro and macroeconomic level, these carbon deduction strategies may be hard to attain, because it depends on an area or country's specific circumstances [6].
Along with concerns about energy crisis and environmental impact, more researchers pay attention to the energy saving logistics. Since the logistics industry is a significant source of energy consumption and carbon emission, it is becoming more important to relieve the situation in logistics operation. Chen et al. [7] found that the ordinary navigation system generally provides the shortest path for drivers as maintained by geographic maps; however, the provided path may be slower because of traffic congestion. This study introduced a navigation system to determine a fuel saving transportation route by collecting and analysis real time traffic data. Bektaş and Laporte [4] used an exact method by using the CPLEX program to solve the Pollution-Routing Problem (PRP), a special case of the classical Vehicle Routing Problem (VRP) with the objective function that considers for the sum of GHG emissions, travelling times, travelling cost and fuel. Xiao et al. [8] proposed a Simulated Annealing (SA) heuristic with a hybrid exchange rule to solve the Capacitated Vehicle Routing Problem (CVRP) with the objective of minimizing fuel consumption. Furthermore, Wang and Li [9] introduced a hybrid heuristic algorithm with two phases to solve the low carbon emissions model of Location Routing Problem with heterogeneous fleet, simultaneous pickup-delivery and time windows. They combined the Genetic Algorithm (GA) and Variable Neighborhood Search (VNS) into the method. From the literature, the meta-heuristics method is quite popular and effective to solve the energy and pollution minimization model.
The classical LRP consist of finding facility locations, assigning customers and determining vehicle routes to transport product. The LRP is classified as a NP-Hard problem and usually resolve by clustering-based, iterative or hierarchical meta-heuristics (Nagy and Salhi) [10]. There are extensive studies conducted on LRPs and the majority of research concerned with the distributed problem. Alumur and Kara [11] proposed a new model for the hazardous waste location-routing problem to manage waste hazardous. The model is to find the solution of problems: the treatment centers location and suitable technologies, disposal centers location, hazardous waste transportation route to which of the suitable treatment technologies, waste residues transportation route to disposal centers. Karaoglan et al. [12] studied an LRP with simultaneous pickup and delivery (LRPSPD), which is included decisions for finding depot locations and planning transportation routes. The model minimized overall cost with as a special feature pickup and delivery customer demands, which must be executed with same vehicle while minimizing the overall cost. In a contrasting problem, Albareda-Sambola et al. [13] studied the multi-period LRP, in which travelling costs were combined together with location costs to plan the operating facility pattern through a time horizon. Moreover, depots were able to be closed or opened but only in a subgroup of time periods. Samanlioglu [14] also proposed the multi-objective location routing model which consisted of three objectives; (1) total cost minimization, which included total travelling cost of waste residues and hazardous materials plus fixed cost of the treatment construction, recycling and disposal center, (2) total transportation risk minimization associate with the population exposure throughout transportation routes of waste residues and hazardous materials and (3) total risk minimization associate with the population located around the disposal center and treatment center. Vincent and Lin [15] introduced the open location-routing problem (OLRP), which is an extension of the LRP. The special characteristic of OLRP is that after servicing all customers, vehicles do not return to the depot. The objective of this model is to minimize the total cost, consisting of transportation costs, depot opening costs and vehicle fixed cost. Moshref-Javadi and Lee [16] studied the Latency Location Routing Problem (LLRP), of which the objective was different from the traditional LRP model. LLRP focused on minimize latency and waiting time of recipients in the post disaster area. Furthermore, the depot location and distribution routes needed to be decided. Moreover, Schiffer and Walther [17] presented an Electric Vehicles (EVs) location routing problem. EVs had a limited driving distance, thus, the location of charging station and travelling routes have to be designed. The objectives of this model were to minimize travelling distance, number of vehicles used and number of charging stations. They also proposed contrasting recharging alternatives owing to real world restrictions.
From the logistic viewpoint, most of the studies formulated LRP in the context of the distribution problem, such as supply products or commodities to customers and when dealing with reverse flow. Thus, it is in the field of hazardous waste management. Zhao and Verter [18] provided an analytical framework of used oil location routing problem to find the solution of these problems: the used oil storage location, treatment center and disposal center location, the capacity level of used oil for these centers and the used oil transportation routes in collection network. They presented a model for the LRP problem to minimize the total cost and an environmental risk. In addition, Vidović at al. [19] proposed the 2 Echelon Location Routing Problem (2ELRP) model of non-hazardous recyclables collection, which aimed to maximize the profit calculated as income acquired from the collected recyclables. However, there were only few studies in reverse flows of non-hazardous products same palm oil problem.
The traditional LRP formulates an objective function to minimize the total cost, which consists of the following three aspects; (1) location opening cost, (2) the fixed cost associated with vehicle uses and (3) the travelling cost which these objectives have been widely studied. As mentioned above, the risk and time minimization are taken into the account while the total transportation cost reduction is still popular [11][12][13][14][15]20]. However, there were few papers that studied fuel consumption minimization.
Most studies on LRP proposed meta-heuristic method as the solution approach to solve the problem. Since the LRP is an NP-hard problem, the model is unable to find optimal solutions for medium and large size problem. Therefore, meta-heuristics seem to be the suitable method to find the optimal solutions in appropriated processing time. Prodhon and Prins [21] compared the recent meta-heuristics on LRP, which considered two indicators to assess the performance for each method: (1) the percentage gap in either the best known solutions from the previous studies or the lower bounds and (2) the calculation time in seconds. The result showed that meta-heuristics is effective to solve LRP problem.
An adaptive large neighborhood search (ALNS) was first proposed by Ropke and Pisinger [22]. Since then, ALNS has been successfully used to solve various combinatorial optimizations. Yu et al. [23] applied ALNS to solve a robust gate assignment problem; the results were compared with the benchmark algorithm and showed the competitiveness of the proposed ALNS in solving the considered problem. Carvalho and Santos [24] also used ALNS to solve the electronic circuits design problem. Their objective was to connect its nets. The proposed algorithm was tested on seven datasets from the literature. The proposed ALNS improved the best known result of some instances. Furthermore, He et al. [25] applied ALNS to the satellite scheduling problem. They used the destroy and repair concept to design the algorithm and combine it with the excellent perturbation method. The ALNS performance is more efficient than competing previous studies of satellite scheduling problem.
Many researches have been successful by using ALNS to solve transportation problem. Riberiro and Laporte [26] presented ALNS heuristic for the cumulative capacitated vehicle routing problem (CCVRP). Their objective was not to minimize the total routing cost, but to minimize the total arrival times at customers. The ALNS was tested with state-of-the-art instances and compared with memetic algorithms in related studies. Azi et al. [27] proposed ALNS for solving the Vehicle Routing Problem with Multiple routes (VRPM). VRPM is the extension of VRP, but the vehicles are allowed to travel multiple routes during a working day. The multi-objective model was presented; the number of customers was maximized and the transportation distance was minimized. The results show the benefit of ALNS for this approach. Emec et al. [28] developed ALNS for solving an E-grocery Delivery Routing Problem by presenting new destroy and repair operators to select, insert and allocate vendor. The results showed that the proposed solution is effective to obtain high quality solutions rapidly. Li et al. [29] also proposed ALNS to solve the pickup and delivery problem with time windows. The algorithm performance show that the ALNS remarkably outperforms the competitor, both in solution quality and processing time. Alinaghian and Shokouhi [30] combined ALNS and Variable Neighborhood Search (VNS) to solve the model for multi-depot multi-compartment vehicle routing problem. The computational results show good performance of the proposed hybrid algorithm. Sirirak and Pitakaso [31] also developed ALNS to manage marketplace location decision making and tourism route planning in the Chiang Rai province, Thailand. Six destroy methods and five repair methods were applied to solve the tourism routing problem in order to find the best travelling route that tourist can visit all attractive places. The result showed that the proposed ALNS provides an effective solution for the tourism route problem. The previous research demonstrated the effectiveness of ALNS method which often used for solving transportation problems. Although the solution obtained by ALNS is not an exact solution, the result is statistically accepted. Therefore, the ALNS was chosen to solve the LRP problem for this case study.
Reviewing literature led to the following LRP observation: (a) for the single objective LRP model, most studies aimed to minimize total cost, which were comprised of opening cost, fixed cost of vehicle uses and travelling cost and (b) a few studies deal with the reverse flows of non-hazardous products; most of them deal with the recyclable product.
This paper includes the following aspect considering the contributions of this paper: (1) the real-world problem, which is the special case of the LRP that has been introduced, (2) the objective of the model is to minimize the fuel consumption cost which relies on the road and vehicle conditions and (3) the characteristic of problem is the palm oil collection as the reverse flows logistic system of non-hazardous product.

Problem Statement
The case study consisted of 80 palm oil farms. These farms can be selected as the collection center. We aimed to choose the suitable location and set transportation route as shown in Figure 1. The results showed that the proposed solution is effective to obtain high quality solutions rapidly. Li et al. [29] also proposed ALNS to solve the pickup and delivery problem with time windows. The algorithm performance show that the ALNS remarkably outperforms the competitor, both in solution quality and processing time. Alinaghian and Shokouhi [30] combined ALNS and Variable Neighborhood Search (VNS) to solve the model for multi-depot multi-compartment vehicle routing problem. The computational results show good performance of the proposed hybrid algorithm. Sirirak and Pitakaso [31] also developed ALNS to manage marketplace location decision making and tourism route planning in the Chiang Rai province, Thailand. Six destroy methods and five repair methods were applied to solve the tourism routing problem in order to find the best travelling route that tourist can visit all attractive places. The result showed that the proposed ALNS provides an effective solution for the tourism route problem. The previous research demonstrated the effectiveness of ALNS method which often used for solving transportation problems. Although the solution obtained by ALNS is not an exact solution, the result is statistically accepted. Therefore, the ALNS was chosen to solve the LRP problem for this case study.
Reviewing literature led to the following LRP observation: (a) for the single objective LRP model, most studies aimed to minimize total cost, which were comprised of opening cost, fixed cost of vehicle uses and travelling cost and (b) a few studies deal with the reverse flows of non-hazardous products; most of them deal with the recyclable product.
This paper includes the following aspect considering the contributions of this paper: (1) the realworld problem, which is the special case of the LRP that has been introduced, (2) the objective of the model is to minimize the fuel consumption cost which relies on the road and vehicle conditions and (3) the characteristic of problem is the palm oil collection as the reverse flows logistic system of nonhazardous product.

Problem Statement
The case study consisted of 80 palm oil farms. These farms can be selected as the collection center. We aimed to choose the suitable location and set transportation route as shown in Figure 1. The objective was to minimize the total cost. According to fuel consumption cost, we categorized the roads into five types, which were classified by using the average driving speed and the fuel consumption rate that was adopted from Akararungruangkul and Kaewman [32]. Table 1 show the road type and fuel consumption rate. The objective was to minimize the total cost. According to fuel consumption cost, we categorized the roads into five types, which were classified by using the average driving speed and the fuel consumption rate that was adopted from Akararungruangkul and Kaewman [32]. Table 1 show the road type and fuel consumption rate. Because of road condition, the average speed of vehicles is different on each road type. The wider the road is, the faster a vehicle can travel. We collected the data of speed used for the road, then calculated the average speed. After that, road type and fuel used were converted by using average speed data. The fuel consumption rate the amount (liter) of fuel used after travelling 1 kilometer. A road type metric of 6 palm oil farms was used, as shown in Table 2. Table 2. Road type from farm to farm.
From the distance metric shown in Table 3, the fuel consumption rate was multiplied by the distance of farm location. Finally, we obtained the fuel consumption metric of the road for travelling between six farms as shown in Table 4. For example, the distance from farm 1 to farm 2 is 12 kilometer; the type of road is A and the fuel consumption rate is 0.112 liter/km. Thus, the fuel consumption was 12 (0.112) = 1.344 liter.  For the traditional LRP, the travelling cost would directly depend on the distance, but this is not the case for the current LRP model. The distance from farm 2 to farm 3 and the distance from farm 2 to farm 6 are 23 km and 19 km, respectively. On the other hand, the fuel consumptions from 2-3 and 2-6 are 2.070 liter and 2.128 liter, respectively, which shows that in some routes, the fuel consumption is not increased because of the distance, but rather caused by road conditions.
Another example: we consider a different routing in case farm 1 was selected to be the collection center. Therefore, the two possible routings were 1-4-2-6-5-3-1 and 1-4-3-2-5-6-1, which have a total distance of 156 km and 158 km and fuel consumptions of 15.210 liter and 14.944 liter, respectively. However, the traditional LRP model selects routing of 1-4-2-6-5-3-1 because it concerns with the distance minimization. On the contrary, fuel consumption model chooses routing 1-4-3-2-5-6-1, which found that the shorter path consumes more fuel than the longer path. Nevertheless, this study focused on fuel consumption.
Vehicle condition is also a factor of fuel consumption rate. For example, heavy duty trucks usually consume more fuel than small trucks. Moreover, the age of the trucks also affects fuel consumption. The older the trucks, the more fuel is needed due to the decrease of capability of engine. This paper assumed an attribution of the vehicle, then classified them into three types categorized by age and different fuel rate consumption as shown in Table 5.
Objective function Constraints k∈K i∈V The objective function (1) measures the total cost consist of the fixed cost of opening the collection center and the fuel consumption cost using for transportation and the fixed cost related with vehicle uses. Constraint (2) ensures that each customer is in exactly one route and each has only one predecessor in the route. Constraints (3) and (4) are capacity constraints related to the vehicles and collection center. Constraints (5) and (6) assure the continuity of each route and terminate at the depot where the route starts. Constraint (7) is a sub-tour elimination constraints which said that for any subset S of the set of customers J and for any route k, the number of arcs belonging to route k that connect the members of S, must not exceed the cardinality of S-1. Constraint (8) ensures that a customer must be assigned to a depot if there is a route connecting them. Finally, constraints (9), (10) and (11) specify the binary variables used in the formulation.

Solution Approach
We propose a heuristic method based on Adaptive Large Neighborhood Search (ALNS) as the solution to this problem. The ALNS framework was presented by Rope and Pisinger [22] to solve the variant of combinatorial optimization problem such as VRP and LRP. They developed the Large Neighborhood Search (LNS), which was introduced by Shaw [33] for the Capacitated Vehicle Routing Problem (CVRP). The procedure starts with an initial solution and gradually improves objective value. A destroy and repair operator are applied in each iteration. The destroy operator removes a small group of customers randomly, while the repair operator reinserts them to improve objective value. Furthermore, Ropke and Pisinger [22] proposed an ALNS to improve LNS, which allowed the ALNS to use several destroys and repair methods in the same search iteration. In each iteration, the operator would destroy part of the current solution and repair it in a different way to create a new solution. Destroy and repair operators are selected according to an adaptive probability mechanism, and in each iteration, the operator has a probability of which the weights are adjusted dynamically depending on how well it has performed in the past.
In this paper, however, we have adapted the ALNS algorithm of Lutz [34], in which the ALNS starts with initialization of solution then select destroy and repair operator applying to the solution.
Let D = {d i | i = 1, . . . , k} be the set of destroy operator and R = {r i | i = 1, . . . , l} be the set of repair operator. The initially equal weights of the operators are denoted by w(r i ) and w(d i ). During the runtime, these weights are adjusted periodically and an updated period contains p u iterations. The selection of these operators of each iteration is based on their weights. The structure of ALNS is shown in Algorithm 1.
if c(s) < c(s min ) then s min = s adjust the weight w and probability p of the operators return s min

Initial Solution
We constructed the initial solution with the following steps: Step 1: Sorting the candidate locations as the product quantity by descending order, then calculate the accumulation probability of each based on their product quantity. Next, randomly select the location of the collection center by using a roulette wheel. The location thus has more product quantity as well as more chance to select.
Step 2: Assigning a farm to collection center follow fuel consumption; the farm that has the lowest fuel cost can be selected first.
Step 3: Assigning to the next farm through sequence of fuel consumption. Meanwhile, the total amount of products needs to be checked to match with vehicle and collection center constraints. When the capacity has satisfied the constraint, the assignment can be stopped.
Step 4: If there are farms left over, return to step 1-3 to open the next collection center and assign tasks for the farms.
Step 5: If there are no other farms, the procedure is terminated.

Degree of Destruction
Degree of destruction (d) is the level of eradication that is applicable to the current solution.
For example, the current solution contains 30 farms and the degree of destruction is equivalent to 20%, which means we need to remove six farms (30 × 0.2) from the current solution. Then, set q = 6 and put them into the customer pool, and the repair operator puts them back and creates a new solution. We designed the selection degree of destruction randomly from the interval {10%, 15%, 20%, 25%, and 30%} in every iteration.

Destroy Operators
The following describe the destroy operators used in the algorithm.

Random Removal (RaR)
This operator is a simple method in which q farms are chosen randomly and removed from the solution, as shown in Figure 2.
Step 1: Randomly choose q farms from the current solution.
Step 2: Remove them from the current solution to the customer pool.
The following describe the destroy operators used in the algorithm.

Random Removal (RaR)
This operator is a simple method in which q farms are chosen randomly and removed from the solution, as shown in Figure 2.
Step 1: Randomly choose q farms from the current solution.
Step 2: Remove them from the current solution to the customer pool.
Random farm for removal (let q = 2) Farm removal from solution

Worst Removal (WoR)
The worst removal entails removing the worst value in order to improve the solution. This operator considers the highest cost in current solution, shown in Figure 3.
Step 1: Calculate the total fuel consumption of this routing by using data in Table 4

Worst Removal (WoR)
The worst removal entails removing the worst value in order to improve the solution. This operator considers the highest cost in current solution, shown in Figure 3.
Step 1: Calculate the total fuel consumption of this routing by using data in Table 4 Step 4: Ranking the farms by considering the different value by decreasing order. Thus, the list of farms who have to be removed are number 3,4,2,5 and 6 respectively and set to q = 2; thus farms number 3 and 4 will be removed.
The following describe the destroy operators used in the algorithm.

Random Removal (RaR)
This operator is a simple method in which q farms are chosen randomly and removed from the solution, as shown in Figure 2.
Step 1: Randomly choose q farms from the current solution.
Step 2: Remove them from the current solution to the customer pool.

Random farm for removal (let q = 2)
Farm removal from solution

Worst Removal (WoR)
The worst removal entails removing the worst value in order to improve the solution. This operator considers the highest cost in current solution, shown in Figure 3.
Step 1: Calculate the total fuel consumption of this routing by using data in Table 4

Related Removal (ReR)
This operator is similar to the one used by Hemmelmayr et al. [35]. A farm is randomly selected for removal from the solutions as shown in Figure 4.
Step 1: Randomly choose a farm (a seed farm), for instance number 2.
Step 2: Select q-1 farms that consume the least amount of fuel when traveling to the seed farm. For example, the fuel consumptions when traveling from farm number 2 to numbers 3,4,5 and 6 are 2.070, 2.548, 1.568 and 2.128, respectively (see Table 4). Select two farms (q − 1): numbers 5 and 3.
Step 3: Remove those farms from the current solution to the customer pool. Step 1: Among all the opened collection centers, randomly choose one and close it.
Step 2: All assigned farms are removed and put into the customer pool.
Step 3: Randomly choose another one which is not opened yet and open it.
This operator is also important for diversification. Figure 5 illustrate the method of this operator.

Repair Operators
The repair operators re-insert the customers from the pool to the solution.

Random Insertion (RaI)
The farms which are removed from the destroy method are inserted randomly in both position and order as shown in Figure 6.  Let the current solution be {1-3-5-2-4-6-1} and q = 3.
Step 1: Randomly choose a farm (a seed farm), for instance number 2.
Step 2: Select q − 1 farms that consume the least amount of fuel when traveling to the seed farm. For example, the fuel consumptions when traveling from farm number 2 to numbers 3,4,5 and 6 are 2.070, 2.548, 1.568 and 2.128, respectively (see Table 4). Select two farms (q − 1): numbers 5 and 3.
Step 3: Remove those farms from the current solution to the customer pool.

Collecting Centre Removal (CcR)
Step 1: Among all the opened collection centers, randomly choose one and close it.
Step 2: All assigned farms are removed and put into the customer pool.
Step 3: Randomly choose another one which is not opened yet and open it.
This operator is also important for diversification. Figure 5 illustrate the method of this operator. This operator is similar to the one used by Hemmelmayr et al. [35]. A farm is randomly selected for removal from the solutions as shown in Figure 4.
Step 1: Randomly choose a farm (a seed farm), for instance number 2.
Step 2: Select q-1 farms that consume the least amount of fuel when traveling to the seed farm. For example, the fuel consumptions when traveling from farm number 2 to numbers 3,4,5 and 6 are 2.070, 2.548, 1.568 and 2.128, respectively (see Table 4). Select two farms (q − 1): numbers 5 and 3.
Step 3: Remove those farms from the current solution to the customer pool.

Collecting Centre Removal (CcR)
Step 1: Among all the opened collection centers, randomly choose one and close it.
Step 2: All assigned farms are removed and put into the customer pool.
Step 3: Randomly choose another one which is not opened yet and open it.
This operator is also important for diversification. Figure 5 illustrate the method of this operator.

Repair Operators
The repair operators re-insert the customers from the pool to the solution.

Random Insertion (RaI)
The farms which are removed from the destroy method are inserted randomly in both position and order as shown in Figure 6.

Repair Operators
The repair operators re-insert the customers from the pool to the solution.

Random Insertion (RaI)
The farms which are removed from the destroy method are inserted randomly in both position and order as shown in Figure 6. This operator is a solution improvement to find the best position of insertion considering minimal fuel consumption as shown in Figure 7.

Best Insertion: (BeI)
This operator is a solution improvement to find the best position of insertion considering minimal fuel consumption as shown in Figure 7. minimal fuel consumption as shown in Figure 7.
Step 1: Attempt to insert farms in the routing that generate the lowest fuel consumption.
Step 2: Repeat step 1 until all the removed farms are inserted.

Adaptive Weight Adjustment and Acceptance Criteria
Step 1: Attempt to insert farms in the routing that generate the lowest fuel consumption.
Step 2: Repeat step 1 until all the removed farms are inserted.

Adaptive Weight Adjustment and Acceptance Criteria
Let w(h) is the weight of an operator h which is either a destroy or repair operator and p u is the update duration, i.e., the non-improved number of iteration. The amount of period time the operator h has been used throughout the p u iterations is called u h . The success s(h) of h is set to zero at the start of p u iterations. After h is used in an iteration, s(h) is improved by ω i if the new solution is of corresponding quality: The purpose of making the decision of the acceptance criteria of ALNS is whether to proceed with the previous solution (s) or with the new solution (s'); the purpose of all acceptance criteria is to accept solutions that improve the current situation. However, there are some criteria that allow them to accept non-improving solutions. There are different ways in previous studies which have attempted to solve this problem. Shaw [32] proposed the Greedy acceptance method. This method does not accept any worse solutions than the previous one. This might be the limitation of the search because the worse solutions are not fully considered, but often rejected at first sight. Ropke and Pisinger [22] used a simulated annealing (SA) as an acceptance criteria. This is one of the most popular acceptance method for the ALNS algorithm. Thus, we used the SA acceptance criterion in this paper, in which all improving solutions s' are accepted or rejected with probability p = exp(c(s) − c(s )/kT, where T is the temperature of the considered iteration and k is the constant used in the probability function to determine whether to accept a worse solution or not.

Computational Result
We performed the testing of mathematical model by Lingo program V.11 and the ALNS algorithm (CPU: Intel Core i5 up to 2.7 GHz, RAM: 4 GB). The three datasets were provided and tested which the datasets of sample consisted of small (n = 10), medium (n = 30) and large (n = 50). The parameters of the test instances were randomly generated using uniform distribution. The distances between the farm locations were randomly generated from 10 km to 50 km. The capacity of the collection center and the vehicles included three levels; 10,000 kg, 20,000 kg and 30,000 kg. The palm oil quantities were randomly generated from 500 kg to 6000 kg. The simulation was run for five times and the best solution was selected to be the representative of the algorithm. The stopping criterion was the number of iteration which was set to 500 iterations. Finally, we used the algorithm to solve the case study problem. The result from the Lingo program and the ALNS algorithm are shown in Table 6.

Three Datasets of Sample
According to Table 6, the five small problems shown a global optimal status on Lingo that was only 1 second different in Sm.2, Sm.3 and Sm.5 which could be concluded the Lingo program and the ALNS algorithm calculated the same results for this sample size. As a result, the small-sized sample contain 10 farms, so the problem was not complicated. The optimal solution for all small problems was the average total cost at 1,037,791 Baht. Moreover, there was no gap between Lingo program and ALNS algorithm, thus, the ALNS algorithm was effectively equal to the Lingo program for small-sized samples.
For medium-sized samples, the result from ALNS was a different total cost compared with the Lingo program. Although the ALNS result was slightly higher than the Lingo result, the Lingo program required a longer processing time than ALNS. The Lingo program took an average processing time of 20.03 h but ALNS took only 5.57 min; in other words, it was able to reduce the processing time by 99.77%. Therefore, the Lingo program consumed long processing time as the number of farms increase.
For the large-sized sample, the Lingo program could not find a global or feasible solution even though it had been taking more than 120 h of processing time. The Lingo program was able to find a lower bound, which was compared with the ALNS result. The total cost obtained by ALNS was slightly higher than Lingo result. However, the ALNS algorithm provided the acceptable result with a fast processing time. The ALNS algorithm had taken only 14.47 min of processing time which was able to reduce the computational time by 99.79%. Therefore, the Lingo program seemed ill-fitted for large samples.
We performed the analysis using statistical methods to compare the performance. The result of total cost and processing time of all problem sizes were tested normality with a 95% level of confidence as shown in Table 7 which all the p-values were greater than 0.05. Therefore, the data for all datasets were normal distribution and after that the data of all datasets were tested for the average difference using a paired t-test with a 95% level of confidence and the results were presented in Table 8.
From Table 8, the statistics showed a significant p-value greater than 0.05 of total cost for all datasets. It was insignificant difference between the Lingo program and ALNS algorithm with a 95% level of confidence. The processing time for the small-sizes sample had p-value greater than 0.05, denoted that insignificant difference between two method. However, the processing time for the medium-sized sample and the large-sized sample had a p-value of less than 0.05, which denoted a significant difference between the Lingo program and ALNS algorithm. Therefore, the ALNS algorithm was able to reduce the processing time while finding a total cost was not different compared to the Lingo's result. The ALNS algorithm was the optimal method for solving the location routing problem with a large number of farms.

The Case Study
We applied the ALNS algorithm for solving the case study problem, which was conducted to demonstrate how the propose model could be used for the palm oil collection location routing problem. The opened location, routing plan and total cost in each collection center are summarized in Table 9. There were 80 customers (farm groups) in the case study. The exact method could not find the solution in spite of taking processing time over 120 hours. On the other hand, ALNS took an average of 20.46 min and came along with the solution as showed in Table 10.
As mentioned in Section 3.1, the traditional LRP dealt with distance minimization; however, our model dealt with fuel consumption minimization. The conditions of the road affected to the fuel consumption; thus, considering distance was did not minimize fuel consumption. Table 11 compares the solutions of this study model to the solution from classical LRP; the distance was reduced, and it saved 7.72% of fuel cost. To verify the effectiveness of the proposed algorithm, we compare the ALNG algorithm with the Differential Evolution (DE) algorithm from the literature [32]. This study proposed DE to solve the problem and modified some points in the original DE, called MED. The computational results are shown in Table 12. The ALNS has a percentage difference of fuel usage 7.72%, more than that of DE, MDE-1 and MDE-2. MDE-3 has the best performance in term of the quality of solution with 15.05% of fuel saving. However, the size of problem may effect to the percentage difference. Our problem is smaller, thus the impact of percentage difference is lower than the bigger problem. In term of processing time, DE algorithm set stopping criteria by runtime limitation at 10 min but ALNS set that by 500 iterations. Moreover, both algorithms were tested by different computer specification. Therefore, we only present the performance of solution quality. Figure 8 showed the nature of ALNS during searching the solution of the case study which shows that ALNS continuously improved the solution from until stopping criteria were met. During 1-70 iterations, the solution was improved rapidly. After that, the solution was gradually improved until 250 iterations, after which it was stable.

DE-Fuel Used (L)
ALNS-Fuel Used (L)  Figure 8 showed the nature of ALNS during searching the solution of the case study which shows that ALNS continuously improved the solution from until stopping criteria were met. During 1-70 iterations, the solution was improved rapidly. After that, the solution was gradually improved until 250 iterations, after which it was stable.

Conclusion and Outlook
The purpose of this study was to solve the Location Routing Problem in the case study of a palm oil collection center. The mathematical model has been formulated to minimize the total cost consisted of opening cost, fixed cost of vehicle uses and fuel cost. We have presented an ALNS algorithm for this LRP problem, which used several operators that were simple and relied on few parameters.
Three datasets were tested (small, medium and large-sized samples) and the results were compared to the exact method that was computed by the Lingo program. The objective value from Lingo and ALNS of the small-sized sample were equal. In addition, the processing time of ALNS was slightly higher than the Lingo program. For the medium and large-sized sample, ALNS calculated a higher objective value than the Lingo program but the Lingo program consumed longer processing time. Although the Lingo program had a long processing time, it found only feasible and lower bounds for medium-sized and large-sized. ALNS on the other hand, was able to find an acceptable solution within a short processing time and it was faster than Lingo program. Thus, we could summarize that ALNS was the optimal method for solving large problems.
We have also implemented ALNS of the case study problem, which was larger than the instances. The ALNS produced the solution in an average processing time of 20.46 min, while the Lingo consumed over 120 hours, but was unable to find the solution. Thus, ALNS was efficient to solve the case study problem. Furthermore, we extended a comparison between the model of fuel consumption and distance minimization with respect to fuel cost. The result showed that this model could reduce 7.72% of fuel costs as compared to the classical LRP, which is distance minimization.

Conclusions and Outlook
The purpose of this study was to solve the Location Routing Problem in the case study of a palm oil collection center. The mathematical model has been formulated to minimize the total cost consisted of opening cost, fixed cost of vehicle uses and fuel cost. We have presented an ALNS algorithm for this LRP problem, which used several operators that were simple and relied on few parameters.
Three datasets were tested (small, medium and large-sized samples) and the results were compared to the exact method that was computed by the Lingo program. The objective value from Lingo and ALNS of the small-sized sample were equal. In addition, the processing time of ALNS was slightly higher than the Lingo program. For the medium and large-sized sample, ALNS calculated a higher objective value than the Lingo program but the Lingo program consumed longer processing time. Although the Lingo program had a long processing time, it found only feasible and lower bounds for medium-sized and large-sized. ALNS on the other hand, was able to find an acceptable solution within a short processing time and it was faster than Lingo program. Thus, we could summarize that ALNS was the optimal method for solving large problems.
We have also implemented ALNS of the case study problem, which was larger than the instances. The ALNS produced the solution in an average processing time of 20.46 min, while the Lingo consumed over 120 hours, but was unable to find the solution. Thus, ALNS was efficient to solve the case study problem. Furthermore, we extended a comparison between the model of fuel consumption and distance minimization with respect to fuel cost. The result showed that this model could reduce 7.72% of fuel costs as compared to the classical LRP, which is distance minimization.
The results provided by the ALNS algorithm were useful for palm oil logistic management. The suitable locations were selected to be collection center and arranged the routing plan to collect palm oil from farmers. The total costs of the system were minimized and the farmers would make more profit compared to the past. Palm oil was one of the sources of renewable energy and this research aimed to minimize fuel energy consumed as well. Therefore, this could be implied that this study considered as green logistic management work as well as environmental improvement.
Future research will focus on extending this approach to solve a larger LRP problem or another LRP variant such as 2E-LRP (2 Echelon Location Routing Problem); in the meantime, developing the ALNS is needed. Because we will solve bigger problems, the algorithm should be more efficient to compute both the high solution quality and a short processing time.
The proposed algorithm is not only used to solve the case study problem, but can also solve the similar problems. To develop the performance of the algorithm, a good software design is needed. After software development, the software can be suitable for related businesses. The license of the developed software can be customized and shared for worldwide usage. Therefore, this is not restricted to palm oil businesses, but is also an open opportunity for other fields. We can imply that the open innovation concept has been applied to the palm oil and related business. The main idea of open innovation is to look for new technologies and ideas outside of the firm for other ways to increase the efficiency and effectiveness of their innovation processes. An example of open innovation research studies can be revealed in Yun et al. [36] and Yun et al. [37].
Author Contributions: C.T. and R.P. collected the data, designed the algorithm and performed the experiment, K.S. analyzed the result, while S.K. review and editing.