Shuffled Complex Evolution-Based Performance Enhancement and Analysis of Cascade Liquefaction Process for Large-Scale LNG Production

: Among all large-scale natural gas (NG) liquefaction processes, the mixed fluid cascade (MFC) process is recognized as a best-alternative option for the LNG production, mainly due its competitive performance. However, from a thermodynamic point of view, the MFC process is still far from its potential maximum energy efficiency due to non-optimal execution of design variables. Therefore, the energy efficiency enhancement of the MFC process remains an ongoing issue. The design optimization after fixing the main configuration of the process is one of the most economic, but challenging exercises during the design stages. In this study, shuffled complex evolution (SCE) is studied to find the optimal design of the MFC process corresponding to minimal energy consumption in refrigeration cycles. The MFC process is simulated using Aspen Hysys ® v10 and then coupled with the SCE approach, which is coded in MATLAB ® 2019a. The refrigerant composition and operating pressures for each cycle of the MFC process were optimized considering the approach temperature inside the LNG heat exchanger as a constraint. The resulting optimal MFC process saved 19.76% overall compression power and reduced the exergy destruction up to 28.76%. The thermodynamic efficiency (figure of merit) of the SCE-optimized process was 25% higher than that of the published base case. Furthermore, the optimization results also imply that there is a trade-off between the thermodynamic performance improvement and the computational cost (no. of iterations). In conclusion, SCE exhibited potential to improve the performance of highly nonlinear and complex processes such as LNG processes.


Introduction
In recent years, there has been growing interest in the consumption of natural gas [1] as markets continue to expand [2,3]. Global consumption of natural gas increased by 1.7% in 2015, representing 23.8% of primary energy consumption [4]. The origin of this growing interest is the negative environmental effects of oil and coal [5,6], which has resulted in the implementation of the Kyoto Protocol [7], stringent environmental regulations, enhanced energy security and an increase in the proven reserves of natural gas [1]. Global gas reserves are steadily increasing without any immediate threat of depletion. There was an increase in proven reserves from 120 trillion cubic feet (TCF) in 1995 to 186.1 TCF in 2014 [8]. The burning of oil and coal is associated with greenhouse emissions and consequently, climate change and global warming. On the other hand, natural gas has higher energy conversion efficiencies for power generation [9], in addition to lower carbon and oxide emissions compared to oil and coal. However, gas reserves are located remotely, and a significant portion of known reserves is stranded far from market. Moreover, reserves are not distributed globally. All these aforementioned factors play a key role in the introduction and development of liquefied natural gas (LNG) technologies, which have transformed this resource from a regional to international commodity [10]. Therefore, the current trend in terms of LNG production is expected to continue over the next decades.
Liquefaction consumes approximately 30 to 35% of the total required energy for LNG value chain [11]. This high energy consumption makes LNG production an energy-intensive step of LNG value chain [12]. The energy requirement for LNG production depends not only on the liquefaction technology, but also on the site conditions [13] and mixed refrigerant combination with optimal composition [14]. Many liquefaction technologies have been considered for land-based LNG plants that vary in complexity, capacity and efficiency [15][16][17]. These technologies also differ on the basis of their refrigerant type and involved refrigerant cycles [18]. Among mixed refrigerant based processes, propane pre-cooled mixed refrigerant (C3MR), dual mixed refrigerant (DMR), single mixed refrigerant (SMR) and cascade cycles are popular choices [19][20][21][22]. Among all LNG technologies, cascade refrigeration cycles are quite popular for high capacity LNG plants. Generally, two or more refrigeration cycles are used in a cascade system with either a single component (pure) or mixed refrigerant in each loop. It is worth mentioning that increasing the number of refrigeration cycles improves the thermal efficiency of the process and enhances plant performance. However, there is an associated increase in capital cost [23].
Mixed fluid cascade-LNG (MFC-LNG) process consists of three refrigeration cycles and each consists of a different composition and components of refrigerants to formulate a mixed refrigerant. The use of mixed refrigerants with different compositions in pre-cooling, liquefaction and subcooling sections renders the MFC process as a potential candidate for low-cost LNG production, mainly due to the relative high exergy efficiency, which ultimately leads to the low energy consumption. In the evidence of pretext, Vatani et al. [23], compared the performance of the MFC process with the SMR, DMR and C3MR processes and found that MFC process has the lowest energy consumption with the maximum exergetic efficiency, i.e., 51.82%. However, from thermodynamic point of view, the MFC process is still far away from its maximum achievable energy efficiency due to complex and highly nonlinear thermodynamic interactions between design variables, constraints (approach temperature inside multistream heat exchangers) and energy efficiency. Therefore, the performance enhancement of the MFC process by searching best optimal design variables is an ongoing issue, so far. Thus the optimization of mixed refrigerant LNG processes in general [24] and MFC process in particular, is one of the challenging tasks during the design stages.
A review of the literature revealed that there have been optimization studies that focus on achieving optimum operating conditions (temperature and pressure) and mixed refrigerant composition for the LNG processes. For instance, Qyyum et al. [16] presented a comprehensive review focusing on the design optimization of LNG processes. However, the design optimization studies of MFC process are not widely reported in the open literature, mainly due to its high dimensionality and complexity. Ding et al. [25], simulated the MFC process using Hysys and optimized through genetic algorithm (GA) to find the operating conditions and refrigerant composition against minimal compression power. Nawaz et al. [26], used Coggins algorithm to optimize the Hysys-simulated MFC process and demonstrated approximately 11% savings in compression power in comparison to the MFC process optimized by Ding et al. [25]. Ghorbani et al. [27] integrated the MFC process with natural gas liquids (NGL) recovery and nitrogen rejection units (NRU). They simulated the integrated process in Aspen Hysys and then optimized using GA to minimize the compression power. Mehrpooya et al. [28] replaced the precooling compression refrigeration cycles with absorption refrigeration cycle and reduced the compression power up to 30%. In a similar study, Ghorbani et al. [29], integrated NGL recovery with LNG process while replacing absorption refrigeration system in lieu of precooling stage of MFC to lower required energy consumption. Lin et al. [30] designed and optimized the two-stage cascade processes for pressurized LNG production. They conducted simulation-based optimization to find the optimal design corresponding to minimal energy consumption. Recently, Brodal et al. [31] modified the MFC process by replacing the Joule-Thomson (JT) valves with liquid expanders same as did by Qyyum et al. [32] for the SMR process. Moreover, Brodal et al. [31] investigated the impact of atmospheric temperature on the performance of MFC process and optimized the process based on temperature difference at the outlet of heat exchangers. Kamalinejad et al. [33] applied mixed integer nonlinear programming model for the synthesis of refrigeration systems and the procedure was applied on a cascade process. The final design depicted 3.3% reduction in compression shaftwork when compared to the base case.
Since the optimization of LNG processes has presented a significant reduction in exergy destruction and enhanced the energy efficiency. Most of the studies used GA and Hysys-built-in optimizers to find the optimal design variables for the MFC process. Notably there is no algorithm that could provide global optimal solution with decisively conclusion. However, a most widely used approach to evaluate optimization algorithms is just to try them with same basis and highlight its potential superiority for the particular problem based on the observation results. In this context, this study evaluates the shuffled complex evolution (SCE) technique to find the optimal design variables corresponding to minimal compression power. The SCE algorithm has hitherto not been employed for the design optimization of any complex LNG process such as MFC. The SCE is built on the synthesis of four basic proven concepts for global optimization problems including the concepts of random and deterministic approach, clustering, evolution of point complex and competitive evolution. Thus, the efficiency of global convergence is improved [34]. To evaluate and enhance the performance of the MFC process, a well-known commercial simulator Aspen Hysys ® v10 is used. Hysys-simulated MFC process is linked with MATLAB ® 2019a to establish optimization environment. The performance enhancement is analyzed in terms of composite curves, exergy destruction and figure of merit.

Methodology: Shuffled Complex Evolution Approach
The Shuffled Complex Evolution (SCE) (originally proposed by Duan et al. [35]) is an efficient approach for obtaining the global optimal solution of high dimensional and complex optimization problems. The SCE is a search methodology that exploits the strengths of controlled random search (CRS) algorithms in combination with the emerging concepts of competitive evolution and complex shuffling. This technique can effectively use the response surface information to guide the search using a deterministic strategy. The inclusion of random elements is helpful in achieving the flexibility and robustness of the algorithm. In the SCE method, the systematic complex evolution technique ensures that a comparatively robust search is performed. The surfaces of rough, non-convex and insensitive objective functions can be readily used in this approach because of its complex structure. The approach is relatively unaltered by the minor local minima that appear in the way to global solution without the requirement of further derivation. Thus, the efficiency of global convergence is improved. The working flowchart of SCE is shown in Figure 1. The SCE approach achieves the optimal solution in the following steps: Step 1: Initialize p and m; such as p ≥ 1 and m ≥ n+1 where p and m represent the number of complexes and the number of points in each complex, respectively and n is the number of design variables of the MFC process. Generate a sample size S p m   Step 2: A set S of samples 1 2 3 , , ,..., s x x x x in the feasible space n    is generated (feasible space means that sample set of decision variables is generated within the constraints). The corresponding value of objective function fi (energy consumption in our case) is computed at each x sample in the absence of prior distribution and by using uniform sampling distribution.
Step 3: All the samples in set S are sorted in ascending order on the basis of their energy consumption. They can be stored in an array such that x1 represents a sample with minimum energy consumption in the MFC process.
Step 4: Array D is portioned in to p complexes A 1 to A p such that each complex contains m samples, as shown in Equation (1). Here complex or community or population is a subset of array D that contains S feasible samples of decision variables of MFC process and their corresponding energy consumption. It means, in the MFC case, all samples of decision variables in array D are divided into p complexes in such a way that each complex will have few samples of decision variables.
Step 5: Each complex 1,..., k A where k p  evolves according to the competitive complex evolution (CCE) algorithm illustrated separately in Figure 2.
Step 6: . This process is called shuffling of complexes. D is subsequently sorted on the basis of higher energy consumption by the MFC process.
Step 7: At the end, a check is performed, and the process is terminated if the criteria set for convergence is satisfied. Otherwise, the algorithm returns to step 4. In our case, the number of iterations is the termination criteria.
For the evolution of each complex in step 5, a competitive complex evolution (CCE) algorithm, as demonstrated in Figure 2, is used. The CCE working procedure can be explained in the steps as follows: Step 1: The initial values of q, α and β are set such that: 2 , q m   1, 1. and     Step 2: Each complex k A is assigned a weight in terms of a triangular probability distribution, i.e., The point 1 k x has the highest probability, 1  Step 4: In this step, offspring are generated in the following way: (i) B and L are sorted such that q points are arranged in ascending order with respect to their respective function values and centroid g is thereafter computed using the expression given in Equation (2): is computed in the so-called reflection step.
(iii) A mutation is performed in this step such that if r lies within Ω, energy consumption by the MFC process is computed and Step (iv) is performed; else, the smallest hypercube n H   containing k A is calculated and fz is computed by randomly generating a point z within H. Finally, we set is computed followed by fc.
Step ( i ) through ( vi ) are repeated α times, where α is a user-defined parameter that has a value of ≥1.0.
Step 5: Decision variables in B are replaced in A k using the original positions in L. A k is subsequently sorted in order of the increasing energy consumption value.
Step 6: Step 2 through Step 3 are repeated β times where β is a user-defined parameter that has a value greater than or equal to one and determines the number of offspring that should be generated.

Process Simulation and Description
For the specific feed NG (conditions and composition are listed in Table 1), the MFC-LNG process was simulated using Aspen Hysys ® V10. The fundamental assumptions considered while simulating the process are as follows: (i) Peng-Robinson with EOS enthalpy/entropy calculation option was selected for thermodynamic properties calculations.
(ii) The adiabatic efficiency for each compressor was fixed at 80%.
(iii) The minimum temperature approach and pressure drop in LNG exchangers was kept constant at 3 °C and 1 bar, respectively.
(iv) Interstage cooling medium was water, and the water-cooler outlet temperature was 40 °C.
(v) The pressure drop of the water coolers was 0.30 bar.
(vi) Furthermore, it was also assumed that the NG feed composition and conditions remain constant. A basic process flow diagram of the MFC-LNG process, which consists of three compression refrigeration loops, namely, pre-cooling, liquefaction and sub-cooling, is shown in Figure 3. Each cycle entails a different composition of refrigerants to formulate a mixed refrigerant (MR) for that particular cycle. Different MR compositions perform different duties by cooling NG to −25, −100 and −155 °C in pre-cooling, liquefaction and sub-cooling cycles, respectively. NG passes through the LNG-exchangers and exchanges its heat with the respective MR. During this exchange of heat, the temperature of the NG drops to −147.1 °C and it condenses. To lower the pressure of this liquefied NG to the storage pressure, a Joule-Thomson valve (JTV-111) is used, whereby its temperature decreases further as spin-off from reduced pressure. The low-pressure LNG is passed through a phase separator (V-111) and flash gas free LNG is obtained as the final product.
Pre-cooling MR (MR-100) comprises of ethylene (ethane), propane and n-butane. The refrigerant is separated into liquid and vapor phases by passing through a phase separator (V-100). The vapor phase MR passes through a single-stage compression (K-100) and is condensed using a water cooler (E-100). The high-pressure condensed refrigerant is mixed with high-pressure liquid refrigerants in a mixer and cooled using a heat exchanger (E-111). The high-pressure pre-cooling MR (MR-100) after exchanging heat in LNG-100 is expanded using a Joule-Thomson valve (JTV-100). Once the MR is expanded, it vaporizes and its pressure (and consequently temperature) decreases to perform its precooling duty. The low-pressure cold refrigerant passes through the LNG-100 again and absorbs heat from the NG, MR-101 and MR-102. The vaporized, hot and low-pressure MR (MR-100) subsequently returns to the phase separator in order to accomplish the cycle.
Liquefaction MR (MR-101) is composed of ethane as the major constituent along with methane and propane. The vapor phase MR is separated from the liquid phase using a phase separator (V-101). It is subsequently compressed, condensed and mixed with the liquid phase and passed through LNG-100 and LNG-101 before being passed through JTV-101. After passing through JTV-101, the low-pressure MR serves as the cold mixed refrigerant in LNG-101 and is ultimately converted to lowpressure and warm mixed refrigerant after an exchange of its cold energy.
The sub-cooling MR (MR-102) consists of low boiling temperature refrigerants, namely methane, ethane and nitrogen with methane as the major and nitrogen as the minor constituent. After the detachment of vapor from a liquid phase in a separator, it is compressed, condensed and mixed again with the liquid phase before passing through all LNG exchangers. Finally, the MR is passed through JTV-102 to lower its pressure (and consequently temperature) to facilitate its sub-cooling duty. The low-pressure, vaporized MR enters LNG-102 to exchange cold energy with NG and liquefies it before moving to the phase separator as a low-pressure super-heated MR to complete the cycle.

Simulation-Optimization Environment
The key design variables that govern the performance of MFC process include temperature, flow rate and condensation and evaporation pressure of the MR for the pre-cooling, liquefaction and subcooling cycles. Given that these key design variables have an effect on the overall process performance, their optimization to achieve specific energy consumption at lowest level is of utmost importance. Since, objective of this optimization is to minimize energy consumption for natural gas liquefaction, therefore, this task is selected as objective function.
Given that the flow rates of MR, their pressures and temperature significantly affect the rate of energy consumption, are chosen as decision variables for this optimization problem. Lower and upper limits of these variables are tabulated in Table 2. The optimization is further constrained by the cryogenic heat exchanger MITA value to 3.0 °C. In addition, the liquid entry into the compressor inlet stream must be restricted during normal operation of the compressors.

Results and Discussion
MFC has several design variables and many constraints. Nonlinear thermodynamic interactions among decision variables cause a notable increase in complexity and thus, rigorous optimization with the objective of selecting optimal values for the variables is required. In this study, an SCE optimization algorithm was used to address the optimization problem. SCE-optimized values of the design variables corresponding to a different number of iterations and their overall impact on the objective function is presented in Table 3. Furthermore, the temperature and pressure corresponding to each stream of base case, SC-optimized with 300 iterations and 700 iterations are listed in Table 4. Table 3, the methane, ethane and n-butane have the highest contribution (flowrates) in the MR for subcooling, liquefaction and precooling purposes, respectively, for both the SCEoptimized process and the base case. However, a comparison of the total MR flow rate of the base case with that of the SCE-optimized case corresponding to different iterations in all the three regimes shows that the SCE-optimized flow rates change in an unpredictable manner, which is indicative of the interdependence of multiple variables in a complex manner. In addition, the pressure ratio across the three compressors in three different regimes of the base case MFC-LNG process is much higher than that of the optimized case. A lower compression ratio results in a reduction of the power requirement of a compressor, as confirmed by the specific compression power. The specific compression power for the base case is 0.334 kWh/kg-NG compared to 0.284 and 0.268 kWh/kg-NG corresponding to 300 and 700 iterations, respectively. It is worth mentioning that the SCE-optimized processes corresponding to 700 iterations can reduce the specific compression power by more than 19% compared to the base case. This reduced specific energy consumption creates an opportunity for energy saving and ultimately, the reduction of the operational cost of LNG production. It was also observed that specific compression power decreased with increasing number of iterations, but at computational time expense and consequently, computational cost. The computational cost for 700 iterations is around 3-fold higher than that for 300 iterations. Hence, there is a trade-off between the minimization of objective function (specific compression power) and the computational cost.

As seen in
The efficiency of heat exchange systems is considered as a key variable in determining operational cost of LNG production due to its effect on the required refrigeration. Composite curve analysis is widely used for this purpose (to describe the efficiency level). There were many studies on the efficiency of heat exchange systems based on the analysis of composite curves [36][37][38]. The gap between the temperature difference composite curves (TDCC) represents the impact of the MR components on the performance of cryogenic heat exchangers in terms of the heat transfer efficiency, while the temperature-heat flow composite curves (THCC) indicates the work lost due to entropy generation. Generally, the area/gap between composite curves is directly related to the irreversibility of the process and any change in irreversibility has an inverse relationship with the heat transfer efficiency level, due to the work lost. The larger the gap, the lower the efficiency and vice versa.
To evaluate the effect of SCE approach on the efficiency of cryogenic heat exchangers, TDCC and THCC of the base case and SCE-optimized processes are shown in Figures 4-6 for LNG-100, LNG-101 and LNG-102, respectively.

Streams
Base Case MFC_300 MFC_700

Pre-Cooling Cycle Composite Curves
Pre-cooling regime composite curves corresponding to LNG-100 are shown in Figure 4. A larger gap and minimum approach temperature larger than 5 °C are evident from Figure 4a. In addition to the larger gaps between composite curves, the curves started deviating from each other as well, at the hot end. The former represents an optimization opportunity to reduce the gap by tuning the individual refrigerant flow rates using the new constraint of 3 °C MITA. A close approach temperature between the hot and cold streams leads to efficient cooling of NG.
TDCC and THCC for SCE-optimized process corresponding to 300 and 700 iterations are shown in Figure 4c, e and Figure 4d, f, respectively. It is obvious from the comparison of THCC and TDCC that the difference between optimized curves is lesser than that of the base case. In addition to this reduced gap, the curves come close to each other again after some deviation. A comparable gap for TDCC; and the reduced gap in center and re-joining the curves at hot end of THCC suggest that the optimized process has an efficient heat flow, and ultimately improved the thermodynamic performance. Given that the generation of low temperatures costs more energy, a minimum gap between THCC leads to enhanced exergy efficiency especially at low-temperature ends, as confirmed by Figure 4.

Sub-Cooling Cycle Composite Curves
TDCC and THCC for the base case and SCE process for the sub-cooling regime are shown in Figure 6. The regime temperature ranges from −160 °C to −100 °C with a minimum approach temperature at the hot end. A larger gap between TDCC and THCC of the base case suggest an optimization opportunity for energy saving. After optimization, the gap between the TDCC and THCC was reduced in comparison to the base case and the curves come closer to each other. In addition, a minimum approach temperature and efficient heat transfer occurred at the ends of the exchanger. For these exceptions with a minimum approach temperature, a major part of the length of the exchanger depicted non-optimal execution. Consequently, THCC exhibits higher gaps in composite curves, which lead to exergy destruction.
It is quite apparent from Figures 4-6 that the gap between the optimized process composite curves is smaller than that of the base case, and the curves approach each other after optimization. In addition, the smallest temperature difference is observed for the pre-cooling regime and the largest for the sub-cooling regime. Given that the temperature difference plays a major role in minimizing irreversibility, the pre-cooling regime is the major contributor in minimizing the irreversibility of the entire process, and thus the optimized process becomes more energy efficient than the base case.

Thermodynamic Performance: Exergy Analysis and Figure of Merit
Exergy analysis and figure of merit (FOM) are pertinent attributes when describing and comparing the thermodynamic performance of liquefaction and refrigeration systems [39]. Exergy analysis of a system reveals its work potential, considering environmental parameters as the reference state. The analysis provides an opportunity not only to evaluate the thermodynamic performance of a cryogenic system, but also identifies possible improvement in the performance of each inescapable equipment.
In general, exergy analysis of a system is conducted to highlight the destruction that occur in the entire system and in each inescapable equipment of the process in particular. In this study, the exergy destruction of each amenity that is intrinsic to the MFC-LNG process is calculated based on exergy balance equations adopted from Venkatarathnam and Timmerhaus [21] and are provided in Table 5. The exergy destruction results associated with each of the MFC-LNG equipment are detailed in Table  6. In addition, this integrated exergy destruction for each amenity in the SCE-optimized MFC-LNG process are depicted in Figure 7. Table 5. Expressions for exergy destruction calculations in different equipment associated with the MFC-LNG processes [21].

Equipment
Exergy Destruction (kW) Exergy analysis of the base case and SCE-optimized process corresponding to distinct iterations revealed that the optimized process operates at a higher efficiency in comparison to the base case, as shown in Table 6. In addition, there is an improvement in the exergetic efficiency with an increase of the number of iterations. Overall process exergy destruction was minimized to around 22.4% and 28.8% corresponding to 300 and 700 iterations, respectively, when compared to the base case. The exergy results also revealed that the exergy destruction corresponding to different individual equipment's intrinsic to the MFC-LNG process are lower for the optimized process. Given that the SCE improvement potential continues to increase with the increasing number of iterations, the SCE-optimized system corresponding to 300 iterations has 22.39% less exergy destruction compared to the base case. Increasing the number of iterations increases the computational time and consequently the overall cost but results in a reduction of the exergy destruction. These findings confirm that there is a trade-off between exergetic efficiency and computational cost.
A percentage wise breakdown of the exergy destruction by each physical equipment associated with the MFC-LNG process is depicted in Figure 7. The main sources of exergy destruction are air coolers, compressors and LNG-exchangers. During the process of heat exchange, entropy generation could have different origins including heat transfer between hot and cold fluids, the exchanger and its surrounding and fluid movement. The temperature difference between the hot and cold fluid is another source of exergy destruction in the process of heat exchange [40].
The FOM is considered to be an important index in describing the thermodynamic efficiency of a liquefaction process because it provides insight into the overall exergy efficiency [41]. It characterizes the thermodynamic efficiency of a liquefaction process using the Carnot refrigeration cycle (ideal liquefaction process) in terms of the work required to achieve liquefaction. Thus, the FOM of a liquefaction process is defined as the ratio of the ideal work Wi to the actual required work Wr for liquefaction, which can be described by Equation (3): The thermodynamic efficiency can be described in terms of FOM, exergy efficiency, percent Carnot and thermodynamic limit [21,42,43]. Thermodynamically, an ideal liquefaction process gives rise to ideal work (minimum work), whereas actual work is used in actual and practical liquefaction processes. The minimum work required for liquefaction is the maximum work that can be achieved from a fluid when it follows reversible processes to reach equilibrium with its surroundings and can be expressed as follows: where ̇ and ̇ are the mass flow rates of the product LNG and feed NG, respectively. In contrast, ̇ and ̇ are the mass exergy rate of the product LNG and feed NG, respectively. When the ambient temperature is 32 °C, the ideal work required to convert 17.82 kg/h NG to LNG via the liquefaction process is 1.955 kW. Due to irreversibility in actual (real) liquefaction systems, the actual work needed for liquefaction (conversion of NG to LNG) is always more than the ideal work. The actual required work and FOM (thermodynamic efficiency) for the base case and SCE-optimized MFC-LNG process are given in Table 7. It is evident from Table 7 that the actual work requirement for the base case is 5.593 kW, including approximately 67% work lost (difference between actual work and ideal work [44]). Given that the minimum reversible work requirement is 1.955 kW, the thermodynamic efficiency of the base case is only 32.8%. Optimization of the refrigerant composition and their flow rates resulted in a significant reduction in the actual work requirement, which also reduced the work destruction and increased the thermodynamic efficiency. MFC_SCE-300 is around 6% thermodynamically efficient compared to the base case with 18% relative improvement in the thermodynamic efficiency. The results also revealed that the thermodynamic efficiency continued to increase with an increase in the number of iterations. The thermodynamic efficiency of MFC_SCE-700 increased to 41%, with 25% relative improvement in energy saving.

Conclusions
Highly interacting and nonlinear thermodynamic interactions in the MFC-LNG process are unavoidable due to the use of MR, which poses a challenge with respect to their optimization. To achieve near-optimal solutions, various optimization methods have been proposed and attempted. In this study, an evolutionary approach is adopted to optimize key decision variables that satisfy constraints on the approach temperature. Shuffled complex evolution (SCE) was coupled with Aspen Hysys and applied to the MFC-LNG process to minimize the energy requirements by enhancing its thermodynamic efficiency. The SCE-optimized MFC process corresponding to 700 iterations consumes 19.76% less energy as compared to the base case. The exergy destruction was reduced from 4.1304 to 2.9425 kW, which represented a 28.76% reduction. Overall, thermodynamic analysis (figure of merit) revealed that the SCE-optimized process that corresponded to 700 iterations has a 25% higher thermodynamic efficiency than the base case. Conclusively, the SCE approach has a potential to improve the performance of mixed refrigerant-based LNG as well as other high dimensional and complex processes through sole design optimization.

Conflicts of Interest:
The authors declare no conflict of interest.