Evaluation of Jet Flooding in Distillation Column Oleﬁns Plant on Naphtha to LPG Feed Substitution

: The naphtha cracking process is the most commonly used technology for the production of ethylene, propylene, mixed C4s (including 1,3-butadiene and other C 4 components), and pygas (pyrolysis gasoline, a mixture of benzene, toluene, and xylene), all of which are oleﬁns. The cracking furnace and distillation columns are the primary operational units. The raw material is cracked and undergoes reactions in the cracking furnaces, while the distillation columns are responsible for separating the products. Raw material costs account for 80% of production costs. There is also the possibility of using LPG as a less expensive alternative to some of the naphtha. However, changing the raw material would affect the operability of the distillation columns and inﬂuence the yield on the cracking side. To determine the optimal naphtha substitution for LPG without causing hydraulic problems (such as jet ﬂooding) in the distillation columns, analysis using simulation tools must be conducted. A reliability model is being developed to simulate the substitution of naphtha with other feed stocks by comparing simulation results with data from the actual plant. The LPG ﬂow is a variable that is freely adjusted to substitute for naphtha. Simulation tools can be used to assess the effects of economically advantageous naphtha substitution for LPG without compromising plant operability. The optimum naphtha substitution rate is 21.14% from the base case, resulting in jet ﬂooding occurring at Propylene Fractionator No. 2. By implementing this substitution, the beneﬁts that can be obtained amount to USD 22,772.02 per hour.


Introduction
Ethylene, propylene, crude C4 (1,3-butadiene), and pygas (benzene, toluene, xylene) are the fundamental building blocks of the petrochemical industry.Naphtha cracking is the process that produces olefins most frequently (Kirk, 1988) [1].Up to 80% of the operational capacity of olefin plants is contributed by naphtha price.Olefins plants around the world offer the potential to replace naphtha feed into LPG due to ample LPG thanks to the shale oil production trend, LPG price decline, and LPG plenty.(Argus, 2018) [2,3].
However, switching from naphtha to LPG feed carries certain risks.The cracking furnace, where the feed undergoes steam cracking reactions, is a critical unit in olefin plants (Schaschke, 2014) [4].The cracking gas is then separated into ethylene, propylene, mixed C4 (including butadiene and other C4 compounds), and pygas (a mixture of benzene, toluene, and xylene streams), which are all pure olefins, in the distillation column.The process flow diagram is shown in Figure 1.

Figure 1.
Olefins plant (PT Chandra Asri, 2021), as indicated in the block diagram product portfolio of olefins plant under study [5].
Olefins plants utilizing LPG feed are expected to experience an increase in propylene and propane production (Argus, 2018) [2].As naphtha is commonly assumed to be the primary feed for most olefin plants, including the plant discussed in this article, it is crucial to consider the effects of LPG cracking on column hydraulics (Fakhroleslam, 2020) [6].
According to Kister (1992) [7], there are six key hydraulic parameters that significantly impact the operability of distillation columns: a. Jet flooding: This occurs when a distillation column is operated with a high vapor load, causing liquid to be transported over the trays.b.Down comer back flood: when a distillation column is run with a substantial liquid and vapor load, liquid accumulates in the top half of the column.c.Down comer choke flood: excessive vapor flow results in liquid buildup in the trays' down comers.d.Weeping: low vapor load results in liquid draining from the tray channel.e. Excessive pressure drops: the pressure drop across each tray should not exceed the design limit.f.
Turndown ratio: it is essential to feed the distillation column within its operational capacity to maintain acceptable efficiency.
The olefins plant under investigation consists of nine distillation columns, with the following configurations (Lummus, 2012) [8]: a. Demethanizer: This nine-bed distillation column operates at a pressure of 5.8 kg/cm 2  and extracts heavier C2 compounds from methane.The top product is methane, with gauge and bottom side temperatures of −53 °C and −131 °C, respectively.b.Deethanizer: A distillation column with 177 sieve trays is used to separate C2 compounds (the top product) from C3 compounds and heavier substances (the bottom product).It operates at a pressure of 21.3 kg/cm  Olefins plants utilizing LPG feed are expected to experience an increase in propylene and propane production (Argus, 2018) [2].As naphtha is commonly assumed to be the primary feed for most olefin plants, including the plant discussed in this article, it is crucial to consider the effects of LPG cracking on column hydraulics (Fakhroleslam, 2020) [6].
According to Kister (1992) [7], there are six key hydraulic parameters that significantly impact the operability of distillation columns: a.
Jet flooding: This occurs when a distillation column is operated with a high vapor load, causing liquid to be transported over the trays.b.
Down comer back flood: when a distillation column is run with a substantial liquid and vapor load, liquid accumulates in the top half of the column.c.Down comer choke flood: excessive vapor flow results in liquid buildup in the trays' down comers.d.
Weeping: low vapor load results in liquid draining from the tray channel.e.
Excessive pressure drops: the pressure drop across each tray should not exceed the design limit.f.
Turndown ratio: it is essential to feed the distillation column within its operational capacity to maintain acceptable efficiency.
The olefins plant under investigation consists of nine distillation columns, with the following configurations (Lummus, 2012) [8]: a.
Demethanizer: This nine-bed distillation column operates at a pressure of 5.8 kg/cm 2 and extracts heavier C2 compounds from methane.The top product is methane, with gauge and bottom side temperatures of −53 A model will be constructed for each distillation column based on actual data and simulated effects on column hydraulics.The objective is to determine the optimal substitution of naphtha with LPG without causing any adverse effects on the distillation column.
When evaluating the performance of distillation columns, the American Institute of Chemical Engineers (AIChE) Equipment Testing Procedures Committee considers the following parameters: ensuring that column performance meets vendor guarantees, identifying capacity bottlenecks, troubleshooting performance issues, determining the operating range of the column, defining optimum operating conditions, developing basic data and correlations for new designs, and calibrating computer simulations for use in optimization, bottleneck analysis, and design studies (CEP AIChE, 2013) [9].Kister (1992) [3,6,10] explained that flooding can occur due to various mechanisms, including spray entrainment flooding, froth entrainment flooding, downcomer malfunctions, and defects in large diameter columns.Further details on these mechanisms can be found in the cited literature.For sieve trays, the entrainment flooding point can be predicted using the method developed by Kister and Haas [3,6,10], shown in Equation (1).This method has been shown to accurately reproduce a large database of measured flood points within a 15 percent margin.
where, d h = hole diameter, mm; σ = surface tension, mN/m (dyn/cm); ρ G , ρ L = vapor and liquid densities, kg/m 3 ; TS = tray spacing, mm; h ct = clear liquid height at froth to spray transition, mm; h ct is obtained from the Equation (2): h ct, H2O derived from Equation (3) In Equations ( 3) and ( 4), Q L = m 3 liquid down flow/(h,m weir length) and A f = fractional hole based on active bubbling area; for example, derived from Equation (5).
The simulation tool used in this study is ASPEN HYSYS V.12, which is capable of simulating the thermodynamic characteristics and equilibrium of distillation column separations (Hanley, 2016) [11].The selection of properties for ASPEN HYSYS simulation is crucial to obtain reliable results (Luyben, 2014) [12].In this investigation, a fluid program was utilized to replicate the thermodynamic model, with guidance from (Yadav, 2020) [13].
Here, the user must select the method based on the component type or process type.If one does not have an idea of selecting the appropriate method, there is an option called method assistant on the page.Method assistant gives the suitable method based on the component properties or specific area of application.
When developing a robust model, it is important to carefully choose the dependent and independent variables (Shinskey, 1991) [14].In this study, the independent variable is the LPG feed flow, while the dependent variables are the naphtha composition and LPG composition.The simulation model should be evaluated by adjusting tray efficiency to minimize the percentage error, and plant data should be used to validate the simulation against actual conditions (Loshchev, 2010) [15].

Data Input
The simulation is based on the following data: a.
Distillation column data sheets, which provide information about the geometry of each column.b.
Actual operating conditions of each distillation column, including pressure, temperature, flow rate, and composition.These operating conditions are derived from actual plant data.c.
Naphtha composition, which serves as the basis for the study.
The naphtha feed composition used in the study is presented in Table 1.Table 1 shows the component composition in the naphtha feed based on lab sampling results.The components are categorized based on the number of carbons present in each compound.The composition is given in weight percentages (% wt).d.Naphtha flow rate basis for study.
To produce 100 t/h of ethylene product, the naphtha flow rate used as the basis for the study, according to actual plant conditions, is 252 t/h.This flow rate is achieved by operating 7 furnaces with naphtha feed.

e.
Feed composition of LPG feed for substitute naphtha use for study.
The LPG feed composition used in the study is presented in Table 2. Table 2 shows the corresponding chemical compounds and their respective composition in the LPG feed, based on laboratory sampling results.f.
LPG flow rate basis for study.
The LPG flow rate basis for the study, depending on the number of furnaces run with LPG feed as a substitution for naphtha feed, is shown in Table 3.The process plant under investigation comprises seven furnaces, each with a maximum capacity of 36 t/h.The percentage substitution of LPG furnaces for naphtha furnaces can be calculated using Equation ( 6): Table 3 shows the flow rates of naphtha and LPG feed based on the number of furnaces running with LPG feed as a substitution for naphtha feed.The flow rates are given in tons per hour (t/h).The last column represents the percentage of LPG substitution to naphtha.

Theory
The Kister and Haas method, along with the Fair's correlation method, are used to estimate the flood capacity of the distillation column.These methods take into account various factors such as column geometry, fluid properties, and operating conditions to predict the flood point.The net area of the column is considered in the flood calculation.
Kister and Haas gave a correlation which is said to reproduce a large database of measured flood points to within ±15 percent.CSB, flood is based on the net area.
By using Equation ( 2) and obtaining the necessary parameters, such as the net area and hydrostatic head Equation (5), h ct was subsequently obtained from Equation (3) and can be applied to calculate the percentage flood for each distillation column.This percentage indicates the extent to which the column is operating at or near its flood capacity, as calculated in Equation (7).It is important to note that the accuracy of these predictions may vary, and it is recommended to validate the results against actual plant data to ensure reliability.

Fluid Package Selection
The selection of a suitable fluid package in ASPEN HYSYS plays a crucial role in accurately simulating the thermodynamic behavior and equilibrium of the system under investigation.The fluid package, or thermodynamic model, serves as the foundation for calculating properties such as phase behavior, heat capacity, enthalpy, and vapor-liquid equilibrium.
In the context of this journal publication, the fluid package selection process was of particular importance.It involved replicating the thermodynamic model using a fluid program, guided by the work of Yadav (2020) [13].The choice of fluid package depends on the specific requirements of the components involved in the separation process or the particular application area.
To assist users in selecting an appropriate fluid package, ASPEN HYSYS provides a method assistant tool.This tool aids users in identifying the most suitable method based on component properties or the specific application area.By utilizing the method assistant, users without prior knowledge or expertise in fluid package selection can make informed choices to ensure reliable and accurate simulation results.
By addressing the significance of fluid package selection and highlighting the availability of the method assistant tool, this journal publication emphasizes the importance of choosing the appropriate thermodynamic model to enhance the reliability and accuracy of the simulation results in the context of the distillation cascade model being studied.A summary of fluid package selection from the method assistant tool for each distillation column under study is shown in Table 4.The selection of the Fluid Package in ASPEN HYSYS for each column was based on the fluid properties outlined in Table 4 and the guidelines provided by Hysys [11].

Model Development
Based on the provided assumptions, the model development for the distillation cascade in the olefins plant aims to achieve the separation of components to obtain pure products.Here is a summary of the assumptions [16,17]: The equilibrium conditions are only valid at bubble point temperature.Once the bubble point temperature is reached, the equations are switched to the equilibrium state.
We assumed that all the vapor entering the stage is condensed until the temperature on each stage reaches the equilibrium state.
Overall, the objective of the distillation cascade model is to accurately simulate the separation process in the olefins plant, allowing for the production of pure products by effectively separating the components.[18] The simulation block under study is depicted in Figure 2, illustrating the configuration of the distillation column in ASPEN HYSYS could be accessed in Supplementary Materials.culated based on the sieve tray characteristics.The temperature of the timated based on the tower bottom pressure.(f) Stream configuration: The distillation column is assumed to have only tillate, and bottoms streams, without any other side streams.(g) Uniform composition: The liquid hold-up in the reboiler, reflux drum tray of the column is assumed to be well-mixed with uniform compos (h) Negligible dynamics: The dynamics of the piping, reboiler, and cond sidered negligible, implying that there are no significant time lags in t por phase dynamics are neglected as they are much faster compared phase.(i) Constant liquid hold-ups: The liquid hold-ups are assumed to be co tray, as well as in the reboiler and reflux drums.(j) Adiabatic column: The column is assumed to be adiabatic, neglecting a from components.Decay heat is also neglected in this case.(k) Actual feed composition: The feed composition to the column is obtain plant operation data, ensuring that the model reflects the real-world c (l) The equations for the non-equilibrium state describe liquid and heat The equilibrium conditions are only valid at bubble point temperature.ble point temperature is reached, the equations are switched to the equ We assumed that all the vapor entering the stage is condensed until th on each stage reaches the equilibrium state.
Overall, the objective of the distillation cascade model is to accuratel separation process in the olefins plant, allowing for the production of pur effectively separating the components.[18] The simulation block under stu in Figure 2, illustrating the configuration of the distillation column in A could be accessed in Supplementary Materials.The data input stages for the simulation block are illustrated in Figure 3.
The data input stages for the simulation block are illustrated in Figure 3.

Model Validation
The simulation of distillation columns using the rate-based approach involves solving a large and highly non-linear system of algebraic equations [19].In this research work, the equations are implemented in the ASPEN HYSYS simulation software, which is a process simulator tool.One advantage of using ASPEN HYSYS is the availability of the AS-PEN Properties Database, which provides reliable thermodynamic properties for the calculations.
To solve the complex system of equations, good starting values are required.In this work, these values are generated through simulations using a model with lower complexity, specifically the equilibrium stage model.The column can be divided into segments, and different correlations for the process hydraulics are implemented in each segment based on the column internal modeling approach derived from the plant data sheet.
To ensure the simulation results match the actual plant results, calculated parameters for the composition of the top and bottom components are compared.The comparison is performed using the Absolute Percent Error (% error) calculation, as shown in Equation (8).The tray efficiency in the ASPEN HYSYS simulator software is adjusted iteratively to minimize the Absolute Percent Error (AAPE, %).
Absolute Percent Error (AAPE, %) = 100 x By iteratively adjusting the tray efficiency and minimizing the Absolute Percent Error, the simulation results can be improved to better match the actual plant results, thereby enhancing the accuracy and reliability of the simulation model.

Sensitivity Analysis Substitution of Naphtha to LPG into % Flooding
To determine the maximum percentage of naphtha substitution with LPG without exceeding 100% flooding, a trial-and-error approach can be employed using the simulation model.The simulation model can be run for different percentages of LPG substitution in incremental steps, and the corresponding flooding percentages can be calculated.Table 3 can be used to guide the trial process by specifying the number of furnaces run with LPG feed and the corresponding flow rates of naphtha and LPG.The simulation model will calculate the flooding percentage for each case.
It is important to note that the simulation should be limited to ensure that the flooding percentage does not exceed 100%.Flooding occurs when the liquid flow rate in the

Model Validation
The simulation of distillation columns using the rate-based approach involves solving a large and highly non-linear system of algebraic equations [19].In this research work, the equations are implemented in the ASPEN HYSYS simulation software, which is a process simulator tool.One advantage of using ASPEN HYSYS is the availability of the ASPEN Properties Database, which provides reliable thermodynamic properties for the calculations.
To solve the complex system of equations, good starting values are required.In this work, these values are generated through simulations using a model with lower complexity, specifically the equilibrium stage model.The column can be divided into segments, and different correlations for the process hydraulics are implemented in each segment based on the column internal modeling approach derived from the plant data sheet.
To ensure the simulation results match the actual plant results, calculated parameters for the composition of the top and bottom components are compared.The comparison is performed using the Absolute Percent Error (% error) calculation, as shown in Equation ( 8).The tray efficiency in the ASPEN HYSYS simulator software is adjusted iteratively to minimize the Absolute Percent Error (AAPE, %).

Absolute Percent Error
By iteratively adjusting the tray efficiency and minimizing the Absolute Percent Error, the simulation results can be improved to better match the actual plant results, thereby enhancing the accuracy and reliability of the simulation model.

Sensitivity Analysis Substitution of Naphtha to LPG into % Flooding
To determine the maximum percentage of naphtha substitution with LPG without exceeding 100% flooding, a trial-and-error approach can be employed using the simulation model.The simulation model can be run for different percentages of LPG substitution in incremental steps, and the corresponding flooding percentages can be calculated.Table 3 can be used to guide the trial process by specifying the number of furnaces run with LPG feed and the corresponding flow rates of naphtha and LPG.The simulation model will calculate the flooding percentage for each case.
It is important to note that the simulation should be limited to ensure that the flooding percentage does not exceed 100%.Flooding occurs when the liquid flow rate in the column exceeds its capacity to separate the components effectively, resulting in reduced separation efficiency and potential operational issues.Exceeding 100% flooding can lead to severe flooding conditions and loss of column performance.
Therefore, during the trial process, the LPG substitution percentage should be monitored to ensure it does not exceed the maximum allowable value that would result in 100% flooding.By carefully selecting the substitution percentages and evaluating the corresponding flooding percentages, an optimal point can be determined where the desired naphtha substitution is achieved without exceeding the flooding limit [6,8].

Economic Calculation
To calculate the benefit of feed substitution from naphtha to LPG, the Equation ( 9) can be used.The benefit is determined by considering the product flow rate, product prices, feed flow rates, feed prices, and utilities cost.
Here's how the calculation can be performed using the provided data from Table 4.

Determine the Relevant Flow Rates
Product flow rate: This is the flow rate of the desired product (e.g., ethylene, propylene, mixed C4, or pygas) obtained from the simulation model.
Naphtha flow rate: This is the initial flow rate of naphtha before substitution, which can be obtained from the plant data or the simulation model.
LPG Flow rate: This is the flow rate of LPG obtained from the simulation model.

Determine the Relevant Prices
Product price: This is the price per ton of the desired product (e.g., ethylene, propylene, mixed C4, or pygas) obtained from the cost reference in Table 4.
Naphtha price: This is the price per ton of naphtha obtained from the cost reference in Table 5. LPG Price: This is the price per ton of LPG obtained from the cost reference in Table 5.

Determine the Utilities Cost
Utilities Cost: This is the cost per ton of product for utilities (e.g., energy, water, etc.) required for the operation of the distillation column, obtained from the cost reference in Table 5.
By plugging in the corresponding values into Equation ( 8), the benefit of feed substitution from naphtha to LPG can be calculated.The result will represent the economic advantage or savings achieved by substituting naphtha with LPG as the feed in the distillation column.

Simulation Matching with Actual Plant Data
Table 5 provides a summary of the assessment of simulation results in ASPEN HYSYS compared to actual plant data, including the efficiencies achieved in each distillation tower.The fluid package used in ASPEN HYSYS was selected based on the guidelines for properties specific to each column under study.The column efficiency was determined through trial and error, following the rule of thumb for distillation column operation.The average absolute percent error (AAPE) was calculated to evaluate the goodness of fit between the simulation and plant data.
Analyzing the efficiencies obtained in each distillation tower (as shown in Table 6), it is evident that they vary across the columns under study.The Demethanizer column achieved an efficiency of 72.87%, while the Deethanizer column had an efficiency of 74.31%.Similarly, the Ethylene Fractionator exhibited an efficiency of 71.92%.It is important to note that the efficiencies were determined through a trial-and-error process, taking into account the specific operational characteristics of each column in order to achieve a minimum error percentage in the simulated column compared with the actual data.These efficiency values highlight the effectiveness of the separation process within the distillation cascade.Table 6.Summary of assessment simulation result in ASPEN HYSYS compare with actual plant data, Fluid Package in ASPEN HYSYS selected based on guideline of properties in each column under study [11], column efficiency as trial result by following rule of thumb in distillation column operating [13], and the average absolute percent error (AAPE) to assess the goodness of fit.In terms of validation, the AAPE values provide an assessment of the agreement between the simulation results and the actual plant data.In this study, the AAPE values for all columns were found to be less than 2%, indicating a good fit between the simulation results and the experimental data.This suggests that the simulation model accurately predicts the behavior of the distillation columns and provides reliable results for the study [20,21].
It is crucial to consider the level of validation achievable through the experimental data, taking into account the obtained error.The AAPE values below 2% indicate a high level of agreement between the simulation and plant data.However, it is important to acknowledge the limitations and uncertainties associated with experimental measurements.Factors such as measurement accuracy, instrumentation, and process variations may contribute to the observed error.Nevertheless, the small magnitude of the AAPE values suggests a favorable level of validation for the simulation results [21].
By discussing the different efficiencies achieved in each distillation tower and considering the level of validation achievable through the experimental data, this journal publication provides a comprehensive evaluation of the simulation results and highlights the reliability of the findings.
Based on the calculation result shown in Table 6, it can be concluded that the model gave a satisfactory result due to AAPE (%) ≤ 2% [6,9].

Sensitivity Anaysis
After the reliable model was developed and tested, sensitivity analysis was conducted to determine the correlation between % flooding based on Equations ( 1)-( 6) with % naphtha substitution to LPG in Table 3.
In the sensitivity analysis, certain variables were held constant while others were subject to variation.The fixed variables include the operational parameters such as column dimensions, feed composition, and operating conditions.The variable of interest is the % of naphtha substitution to LPG, which was systematically varied in order to observe its impact on % flooding.
The R-squared value or root-mean-square values (RMS) measures how close the data are to the fitted regression line.An R-squared value of 1 is an ideal approximation, and 0 indicates that the model explains none of the data points [22].
Based on the calculation results shown in Figures 4-7, all graphs have R-squared values greater than 0.7, indicating satisfactory correlation [10,11].The assessment result of the limitation of % naphtha substitution in each distillation column under study is presented in Table 7.
Based on the calculation result shown in Table 6, it can be concluded that the model gave a satisfactory result due to AAPE (%) ≤ 2% [6,9].

Sensitivity Anaysis
After the reliable model was developed and tested, sensitivity analysis was conducted to determine the correlation between % flooding based on Equations ( 1)-( 6) with % naphtha substitution to LPG in Table 3.
In the sensitivity analysis, certain variables were held constant while others were subject to variation.The fixed variables include the operational parameters such as column dimensions, feed composition, and operating conditions.The variable of interest is the % of naphtha substitution to LPG, which was systematically varied in order to observe its impact on % flooding.
The R-squared value or root-mean-square values (RMS) measures how close the data are to the fitted regression line.An R-squared value of 1 is an ideal approximation, and 0 indicates that the model explains none of the data points [22].
Based on the calculation results shown in Figures 4-7, all graphs have R-squared values greater than 0.7, indicating satisfactory correlation [10,11].The assessment result of the limitation of % naphtha substitution in each distillation column under study is presented in Table 7.     3. The solid line correspondence to regression result and R-squared value is shown in the graph.
Figures 4-7 illustrate the correlation between the percentage of naphtha substitution to LPG (X-axis) and the calculated % flooding (Y-axis) to determine the maximum limit of naphtha substitution percentage that results in 100% jet flooding, acting as a limitation for each column under study.The maximum percentage of naphtha substitution to LPG and the RMS for each column are presented in Table 7. Figures 4-7 could be accessed in Supplementary Materials.
Table 7. Summary of the goodness of fit using the correlation (Equations ( 1)-( 6)) for % flooding as compared to % naphtha substitution with LPG.The goodness of fit is quantified using root-meansquared error (RMS).In each distillation column under study, the maximum amount of % naphtha substitution with LPG when % flooding hit 100% is also presented.3.

Column
The solid line correspondence to regression result and R-squared value is shown in the graph.
Table 7. Summary of the goodness of fit using the correlation (Equations ( 1)-( 6)) for % flooding as compared to % naphtha substitution with LPG.The goodness of fit is quantified using root-meansquared error (RMS).In each distillation column under study, the maximum amount of % naphtha substitution with LPG when % flooding hit 100% is also presented.illustrate the correlation between the percentage of naphtha substitution to LPG (X-axis) and the calculated % flooding (Y-axis) to determine the maximum limit of naphtha substitution percentage that results in 100% jet flooding, acting as a limitation for each column under study.The maximum percentage of naphtha substitution to LPG and the RMS for each column are presented in Table 7. Figures 4-7 could be accessed in Supplementary Materials.

Column
Based on the summary in Table 7, jet flooding limitation in distillation columns under study occurred in Propylene Fractionator No. 2, with maximum naphtha substitution 21.14%.This value, correlated with the norm in the cracking of LPG, will cause a higher yield of propylene [2,3,5] which mainly adds load in Propylene Fractionator No. 2.
The sensitivity analysis allows us to observe how changes in % naphtha substitution to LPG influence the occurrence of flooding in each distillation column.By systematically varying the % naphtha substitution while keeping other variables constant, we can determine the impact and limitations of naphtha substitution on the distillation process.

Economical Evaluation
Economical evaluation being conducted following Equation (8) with data as per Tables 4 and 8.The calculated optimum economic benefit, considering limitation on % flooding from Equations ( 1)-( 6) and Table 8, and using the cost reference from Table 4 by applied calculation based on Equation ( 8), the resulting optimum benefit is USD 22,771.02per hour from the reference basis of 100% naphtha feed consumption.

Summary and Conclusions
In the present study, ASPEN HYSYS was utilized to conduct sensitivity analysis and assess the ability of the validated model in predicting the flooding phenomenon in distillation columns when substituting naphtha feed with LPG.ASPEN HYSYS is a simulation-based tool that incorporates thermodynamic properties' packages, allowing for the selection of appropriate packages based on the system's physical conditions.
By employing modeling, optimization, and prediction techniques within ASPEN HYSYS, the operating conditions and limitations of distillation columns-particularly regarding jet flooding-could be predicted without risking the actual plant conditions.This approach provides a safe and efficient way to explore different scenarios and evaluate the impact of modifications on column internals.
To further enhance the accuracy of the model and ensure its alignment with actual plant conditions, future efforts will be directed towards improving the model's performance, and conducting thorough assessments of the modification effects on column internals.These ongoing developments will contribute to advancing the understanding and optimization of distillation column operations.

Figure 1 .
Figure 1.Olefins plant (PT Chandra Asri, 2021), as indicated in the block diagram product portfolio of olefins plant under study [5].
(a) Ideal solution: Both vapor and liquid phases are assumed to be ideal solutions, and constant average relative volatilities are considered.This assumption simplifies the thermodynamic calculations involved in the separation process.(b) Condenser and operating conditions: The column top temperature is estimated based on the assumed condenser operating pressure.The feed and reflux streams are assumed to be at their dew point and bubble point, respectively.A total condenser is assumed for the reflux condenser.(c) Constant molal overflow: The assumption of constant molal overflow simplifies the calculation of liquid and vapor flows in the column.(d) Stage efficiency: The stages in the column are initially assumed to be 100% efficient with respect to mass transfer.However, adjustments may be made to match the simulation results with actual plant conditions.(e)Pressure drop and reboiler temperature: The pressure drop through the tower is calculated based on the sieve tray characteristics.The temperature of the reboiler is estimated based on the tower bottom pressure.(f) Stream configuration: The distillation column is assumed to have only the feed, distillate, and bottoms streams, without any other side streams.(g) Uniform composition: The liquid hold-up in the reboiler, reflux drum, and on each tray of the column is assumed to be well-mixed with uniform composition.(h) Negligible dynamics: The dynamics of the piping, reboiler, and condenser are considered negligible, implying that there are no significant time lags in the system.Vapor phase dynamics are neglected as they are much faster compared to the liquid phase.(i) Constant liquid hold-ups: The liquid hold-ups are assumed to be constant on each tray, as well as in the reboiler and reflux drums.(j) Adiabatic column: The column is assumed to be adiabatic, neglecting any heat release from components.Decay heat is also neglected in this case.(k) Actual feed composition: The feed composition to the column is obtained from actual plant operation data, ensuring that the model reflects the real-world conditions.(l) The equations for the non-equilibrium state describe liquid and heat accumulation.

Figure 2 .
Figure 2. Simulation block of distillation column under study in ASPEN HYSYS.

Figure 2 .
Figure 2. Simulation block of distillation column under study in ASPEN HYSYS.

Figure 3 .
Figure 3. Input stages for simulation block under study in ASPEN HYSYS.

Figure 3 .
Figure 3. Input stages for simulation block under study in ASPEN HYSYS.

Figure 4 .
Figure 4. Parity plot of % flooding based on Equations (1)-(6) for (a) Demethanizer, and (b) Deethanizer column with % naphtha substitution to LPG in Table3.The solid line correspondence to regression result and R-squared value is shown in the graph.

Figure 4 .Figure 5 .Figure 5 .
Figure 4. Parity plot of % flooding based on Equations (1)-(6) for (a) Demethanizer, and (b) Deethanizer column with % naphtha substitution to LPG in Table3.The solid line correspondence to regression result and R-squared value is shown in the graph.

Figure 4 .Figure 5 .Figure 6 .Figure 6 .
Figure 4. Parity plot of % flooding based on Equations (1)-(6) for (a) Demethanizer, and (b) Deethanizer column with % naphtha substitution to LPG in Table3.The solid line correspondence to regression result and R-squared value is shown in the graph.

Figure 6 .Figure 6 .Figure 7 .Figure 7 .
Figure 6.Parity plot of % flooding based on Equations (1)-(6) for (a) Depropanizer No. 2, and (b) Propylene Fractionator No. 1 column with % naphtha substitution to LPG in Table 3.The solid line correspondence to regression result and R-squared value is shown in the graph.
2 , with top side and bottom side temperatures of −23 °C and 66 °C, respectively.c.Ethylene Fractionator: With 137 sieve trays, this distillation column separates ethylene and ethane.It operates at pressures of 16.48 kg/cm 2 , with top side and bottom side temperatures of −35 °C and −12 °C, respectively.
d. Depropanizer No. 1: This distillation column, equipped with 48 sieve trays, separates C3 compounds (the top product) from heavier C4 compounds (the bottom product).It operates at a pressure of 16.7 kg/cm 2 , with a top side temperature of 44 °C and a bottom side temperature of 82 °C.
. Propylene Fractionator No. 1: With 55 valve trays, this distillation column separates C 3 compounds from propane.It operates at a pressure of 19.7 kg/cm 2 , with top side and bottom side temperatures of 50 • C and 58 • C, respectively.g.Propylene Fractionator No. 2: Propylene and propane are separated using 149 sieve trays in this distillation column, which operates at pressures of 19.2 kg/cm 2 and temperatures of 46 • C (top side) and 50 • C (bottom side).h.Propylene Fractionator No. 3: Propylene and propane are separated using 210 sieve trays in this distillation column, which operates at pressures of 18.3 kg/cm 2 and temperatures of 46 • C (top side) and 58 • C (bottom side).This column can be interchanged with Propylene Fractionator No. 1 and No. 2. i. Debutanizer: This distillation column, equipped with 34 valve trays, separates C4 compounds (the top product) from C5 compounds and heavier substances (the bottom product).It operates at a pressure of 4.34 kg/cm 2 , with top side and bottom side temperatures of 47 • C and 116 • C, respectively.
• C and −131 • C, respectively.Depropanizer No. 2: Using 30 sieve trays, this distillation column separates C 3 compounds (the top product) from heavier C4 compounds (the bottom product).It operates at a pressure of 6.1 kg/cm 2 , with top side and bottom side temperatures of 38.2 • C and 82 • C, respectively.f

Table 1 .
Feed composition of LPG feed component correspondence for chemical compound in naphtha feed, based on lab sampling results.

Table 2 .
Feed composition of LPG feed component correspondence for chemical compound in LPG feed, based on lab sampling results.

Table 3 .
Flow rate LPG feed as naphtha feed substitution based on furnace arrangement.

Table 4 .
[11]ary of column under study, Fluid Package in ASPEN HYSYS selected based on guideline of properties in each column under study[11].

Maximum Percentage Naphtha Substitution to LPG (%) RMS
[2,3,5]n the summary in Table7, jet flooding limitation in distillation columns under study occurred in Propylene Fractionator No. 2, with maximum naphtha substitution 21.14%.This value, correlated with the norm in the cracking of LPG, will cause a higher yield of propylene[2,3,5]which mainly adds load in Propylene Fractionator No. 2.The sensitivity analysis allows us to observe how changes in % naphtha substitution to LPG influence the occurrence of flooding in each distillation column.By systematically Figure 7. Parity plot of % flooding based on Equations (1)-(6) for (a) Propylene Fractionator No. 2, (b) Propylene Fractionator No. 3, and (c) Debutanizer with % naphtha substitution to LPG in Table

Table 8 .
Data input for economic evaluation as required for calculating Equation (8); data calculated from validated model with limitation feed flow rate % LPG from naphtha based on minimization result of Table7.