1. INTRODUCTION
In a traditional supply chain, each member of the chain plans the level of inventory, as well as the number of orders and productions solely based on customer demand. This is mainly due to the fact that the information of each chain member is only related to its direct and immediate customers. Therefore, the productivity of the supply chain is not necessarily at the highest possible level (Chen et al., 2013;Choi et al., 2013). Accordingly, supply chain management (SCM) aiming at making the supply chain more productive and reducing overall costs has become a fundamental issue in the past few years and has been addressed in all industries, especially in the thermal energy industry and related technologies. The supply chain includes all activities from receiving raw materials to delivering the final product. In addition, these activities involve all construction stages, inventory control, distribution, warehousing and customer service. In general, SCM models all activities into one integrated process (Askarany et al., 2010).
2. INTRODUCTION OF THERMAL ENERGY SUPPLY CHAIN
The supply chain is a set of factors that generate added value in the economy. In the current global competition, diverse products must be available to customers based on their desire. Since customer demand for high quality and rapid services has increased pressure on industries, companies can no longer do all the work alone. In the existing competitive market, economic and production enterprises have found themselves in need of the management and monitoring of resources and related elements outside the organization in addition to dealing with their domestic resources and organization. Therefore, SCM has been identified as one of the fundamental aspects of electronic business implementation. SCM enables customers to receive reliable and quick services, along with highquality products at the minimum cost. As such, all activities such as supply and demand planning, material preparation, product planning and production, service, warehousing, inventory control, distribution and delivery to customers, which were formerly conducted at the company’s level, are now transferred to the supply chain level. A key issue in the supply chain is to manage and control all these activities coordinately (Ganeshan and Harrison, 1995;JanvierJames, 2012).
2.1 Two Echelon Supply Chain with Vendormanaged Inventory Policy
A twoechelon supply chain consists of two main components, namely the seller (vendor) and the buyer (consumer). In this system, which is known as the vendor managed inventory (VMI) system, the seller of goods plays the role of service provider and supporter and delivers goods to the buyer (consumer) without intermediaries. This issue is specifically important for marine industry products. Similar to any other supply chain, the goal of a VMIbased supply chain is to reduce the overall costs at various levels of the chain and increase profit using the mentioned policy. The managerial components of the marine industry supply chain are exhibited in Figure 1:
3. LITERATURE REVIEW
In a research, Diabat focused on the twoechelon supply chain using a VMI system. This scholar attempted to design a model in order to find an optimal sales amount. In the end, the model developed in the foregoing research was compared to the traditional method following defining the mixed solution technique (Diabat, 2014). Other studies have been conducted on a twoechelon supply chain, including a research by Nachiappan and Jawahar (2007). The model proposed by these scholars was a VMIbased twoechelon supply chain, and a heuristic algorithm was applied to solve the model. The VMI is one of the systems used for the twoechelon supply chain in various studies. Primarily, VMI is a process where the vendor creates orders for their customers based on the demand information that they receive from the customer. In other words, VMI is a model in which the vendor is responsible for meeting customer demand. Using this concept, Diabat showed that the vendor controls the inventory and makes the orders instead of the customer (Holmström, 1998). Meanwhile, the application of an inventory system that can help maintain, control and distribute goods is critical. In this regard, some of these systems include economic order quantity (EOQ) and economic production quantity (EPQ). In fact, these two inventory control systems support the main VMI system. Diabat (2014) as well as Nachiappan and Jawahar (2007) used the EOQ system in their twoechelon models. Nonetheless, the integration and coordination of supply chain components (e.g., taw material and semifinished goods purchasing, production, distribution and sales of the product) have a significant effect on the decrease of supply chain costs, a decrease of delivery time, and delays in product delivery from distribution centers to markets and sales centers, and consequently, customer satisfaction. In this regard, the VMI system is an important aspect of the supply chain (Zahedi et al., 2021). During the 60s, experts were able to reduce their inventory by studying the internal relationships between warehousing and transportation and their integration, the result of which was identified as distribution management. In the evolution path, the logistics concept was formed by adding the management of construction, procurement and orders to management. The current status (i.e., supply chain) is the result is a combination of different operational loops, the beginning and end of which include vendors and customers, respectively. The main goal of port logistics and supply chain management is to reduce uncertainty and risks in the supply chain. Conventional logistics systems are based on a paradigm that seeks to determine the optimal quantities and location available with the aim of reducing costs and providing excellent services (9). Recent studies, including a research by Pérez Rodríguez (2021), have considered new approaches for the design of a supply chain network, including considering the financial flow in the supply chain and using hybrid metaheuristic algorithms for their optimization.
4. METHOD
This paper designs and models a twoechelon supply chain, which is a type of supply chain formed of two main components, namely the seller (vendor) and the buyer (consumer). One of the fundamental issues in a twoechelon supply chain is VMI, where the vendor creates orders for their customers based on the demand information that they receive from the customer. In other words, the vendor is responsible for meeting customer demand in a VMI. In addition, the vendor controls the inventory and makes the order instead of the customer (Grammenos, 2013;Dong and Xu, 2002;Mastos et al., 2021). Moreover, the seller does all the work, including supplying the product and providing aftersales services, and delivers the product to the buyer (consumer) without intermediaries. In this research, we model and solve the problem of determining the optimal level of multiperiod sales of marine products in a twoechelon supply chain, aiming at maximizing the total profit from sales. In addition to the general information, there is a need for maximum product shelf life.
4.1 Model Premises
Similar to other supply chain models, modeling of a twoechelon supply chain requires some premises. In the present research, the model premises include:

 Customers have different and definite demands in various periods.

 There are certain upper and lower bounds for the sales amount of each seller in each period.

 Each inventory unit can be stored for a maximum of a certain number of periods.
In addition to the mentioned items, the model’s parameters include seller demand, seller capacity, inventory costs and order submission, cost of starting production in each period, and cost of delivery of products from the vendor to the seller.
4.2 Modeling
According to the mentioned premises, the modeling parameters are presented as follows:

I_{jt} : Level of inventory for the jth buyer at the end of the tth period

d_{jt} : The demand of the jth buyer in the tth period

a_{jt} : The yintercept of the pricedemand curve for the jth customer

b_{jt} : The slope of the pricedemand curve for the jth customer

C : Seller capacity

Hb_{jt} : Inventory cost for the jth buyer independently in the tth period

H_{st} : Inventory cost for the seller independently in the tth period

Q_{jt} : Economic order amount in the tth period

Sb_{jt} : Cost of establishment for the jth buyer in each order in the tth period independently

S_{st} : Cost of establishment for seller in each order in the tth period independently

W_{jt} : Cost of the contract between the jth buyer and seller in the tth period

y_{jtmin} : Minimum expected sales for the jth buyer in tth period

y_{itmax} : Maximum expected sales for the jth buyer in the tth period

θ_{jt} : Flow cost per unit, from the seller to the jth buyer in the tth period

v_{jt} : Cost of transportation resources per unit, from seller to the jth buyer in the tth period

δ : Manufacturing cost per unit

τ_{max} : Maximum product storage time
After solving the model, the amount of goods delivery from the vendor to each seller in each period must be determined, which is presented by y_{jt}, which shows the amount of product delivery from the vendor to the jth seller in the tth period.
4.3 Objective Function
As mentioned earlier in the literature review section, the objective function of the problem is to maximize the profit from the sale of products to customers. Seller costs are based on costs related to production, distribution, orders, and inventory maintenance. The objective function of the problem is calculated from the difference between the net profit and the sellerrelated costs. The price in each period will be a function of the amount of demand in that period. If the jth customer demand in the tth period is assumed to be equal to y_{jt}, the cost of the problem will be estimated using the equation below:
where a_{jt} and b_{jt} show the pricedemand yintercept and curve slope for the jth costumer, respectively. According to Equation (1), the selling price of goods to customers in different time periods can be different. Therefore, the gross profit from sales is equal to the product of the number of sales P(y_{jt}) multiplied by the sale price y_{jt} :
 Production Costs
The production costs are obtained by multiplying the cost of manufacturing each product by the sales amount (demand for goods), and this value is expressed in the form of δy_{it} based on the parameters defined in the problem.
 Distribution Costs
The distribution costs are obtained by multiplying the flow cost θ_{jt} y_{jt} by the cost of transportation resources v_{jt}y_{jt}. Accordingly, the distribution costs are estimated, as follows:
 Inventory Ordering and Maintenance Costs
In this study, the EOQ is used to calculate the total inventory ordering and maintenance cost. In this system, the inventory ordering cost is obtained using the equation below:
By placing the related amounts in Equation (4) for ordering cost, we will have:
To calculate the maintenance cost in the same way:
Finally, and based on Diabat’s model (6), the total cost of ordering and maintaining inventory can be calculated by the equation below after replacing Q_{jt} with the optimal amount obtained from the economic order model $EOQ={\left[\frac{2({S}_{st}+{S}_{bjt}){y}_{jt}}{({H}_{st}+{H}_{bjt})}\right]}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{1ex}{$2$}\right.}$.
The objective function (8) minimizes profit from the sale of perishable items at different times. In addition, given the use of coefficients related to the flow cost and the use of the inventory system of the amount of economic order, the optimization problem is of nonlinear type. It is notable that the v_{jt} coefficient includes costs such as transportation costs, human march costs, and administrative costs in each unit. Ultimately, the mentioned coefficient is calculated at 0.5 following reviewing various articles (Dubey et al., 2021;Pahlevan et al., 2021).
4.4 Constraints
The problem of determining the optimal level of multiperiod sales of marine products in a twoechelon supply chain encompasses various constraints, which are classified below:
 Capacity Constraints
In constraints (9), the sales amount related to different periods is considered lower or equal to the sales capacity, as shown below:
 Sales Minimization and Maximization Constraints
Given the diversity of items and different customer demands at various periods, a lower and upper bound is considered for selling different products in the present problem. This is clearly shown by constraints (10):
 Demand Balance Constraints
Constraints (11) express that the total inventory of the previous period and sales of the current period must be equal to the amount of demand and inventory in the present period:
 Symbol Constraints
Constraints (12) show that the values related to sales, inventory and demand must be positive in each period:
5. RESULTS
According to the objective function and constraints of the proposed model, the problem is of constrained optimization type, which is solved by two methods: A) using classic techniques and B) using smart or evolutionary methods. In this paper, PSO metaheuristic methods and CPSO evolutionary method are used to solve the proposed model and finally, the solutions are compared with the method solved by GAMS software. The problems are assumed to have five customers and five periods, and the demand of customers in the five different periods is predetermined. The final solution of the problem for sales amount in 20 problems is presented in Table 2 based on the problem’s parameters. In addition, the results of the comparison of different methods are shown in Table 1.
According to Figure 2, all methods are between the lower and upper bound. In addition, GAMS has the lowest difference with the lower bound in all problems since GAMS accurately solves mathematical models and the lowest amount of objective function is obtained by using this method. Among the metaheuristic methods, the CPSO has the lowest amount in different problems. In addition, the PSO algorithm has very little difference with the upper bound, which shows the weakness of the algorithm. In order to more accurately examine the solution methods, a problem is presented in more detail. The complete solution to the data set of problem 2 is shown in Table 2:
6. CONCLUSION
The critical role of SCM in society is perhaps overlooked and less emphasized. In general, SCM can help people survive by improving the healthcare conditions and protecting people against overindulgence and weather conditions. People trust supply chains to provide necessities such as food and water as well as medicines and health care. The supply chain is also critical to supplying electricity to homes and businesses and provides the energy needed for light, heat, air conditioning and refrigeration. In addition, SCM improves the quality of life by improving employment, thereby laying the foundation for economic growth and improvement of living standards. This leads to many employment opportunities since supply chain experts design and control all supply chains of a community. Moreover, they manage the inventory, warehousing, packaging and logistics. Furthermore, a common feature of most poor nations is the lack of a developed supply chain. Communities with strong and welldeveloped supply chain infrastructures, such as large rail networks, interurban highway systems and a range of airports and ports, can trade goods at lower costs, allowing consumers to buy more products, thus providing economic growth. The present study evaluated a VMIbased twoechelon model for optimal sales management in a supply chain, which encompassed two seller and buyer sections and its objective function related to the sales amount was optimized. In addition, the solution obtained from the accurate method was compared to the solutions achieved from metaheuristic techniques. As observed, the accurate techniques yielded more optimal solutions. However, metaheuristic algorithms could generate a nearoptimal solution in a much shorter period. It is recommended that hybrid algorithms (e.g., simulated annealing and genetics algorithms) be used to evaluate and compare the mentioned solutions.