Column

Figure 1. Distillation Column

Authors: Scott Smith ^[2015] and Harry Poppick ^[2015]

Stewards: Jian Gong and Fengqi You

This article highlights the use of HYSYS for simulations of chemical distillation. The two primary distillation simulation utilities are the shortcut column, which is useful for quick and cursory simulations, and the full column utility, for more involved simulations.

HYSYS Column Types

The table's caption

Column heading 1	Column heading 2	Column heading 3
Row heading 1	Cell2	Cell 3
Row heading A	Cell B	Cell C

Distillation Column HYSYS Simulation

Shortcut Distillation Column

The shortcut distillation utility within HYSYS allows users to quickly model a distillation separation of a stream with relatively little known information. Once the user attaches a fully defined inlet stream to the shortcut distillation column, the essential information that must be specified are as follows:

• The vapor/liquid quality of the top product
• The light key component and composition in the bottoms product
• The heavy key component and composition in the distillate (top) product
• The condenser pressure
• The reboiler pressure
• The external reflux ratio

With this information, the utility will run and, if possible, converge on a solution that will meet the key component specifications. The critical information that the utility will calculate are as follows:

• Fully specified distillate product stream (temperature, pressure, flow rate, composition)
• Fully specified bottom product stream (temperature, pressure, flow rate, composition)
• Condenser operating temperature and heat duty
• Reboiler operating temperature and heat duty
• Minimum number of trays
• Actual number of trays
• Optimum feed tray location

The primary purpose of the shortcut distillation is to allow for users to perform a first pass estimate on the performance requirements of a given system for a specified reflux ratio. The values gathered from the shortcut distillation utility can be used to assist the user when setting up the complete distillation utility in HYSYS, or when attempting to find reasonable column properties to allow the system to converge (e.g. number of trays, feed location, etc.)

Tutorial

This example will simulate the separation of benzene, toluene and p-xylene from some upstream process.

Initialize

Figure 2. Adding Components in HYSYS

1. Start Aspen HYSYS v.8 and click “new” on the main screen to open a new simulation window.
2. Once on the component list screen, select “add” near the bottom of the screen
3. Search the database for the desired components using the search bar, and add them to the component list by clicking “add.” In this case, we wish to add Benzene, Toluene and p-Xylene.
4. On the left hand “properties” menu, click the “fluid packages” tab and select “add” to add a new package.
5. Select the desired fluid package. In this case we are using Peng-Robinson
6. In the bottom, left-hand corner of the screen, click “simulation” to enter the simulation environment.

Setup Streams

Figure 3. Add Streams

1. Start by clicking the blue arrow at the top of the “palette” window and then clicking on the main flowsheet to produce a new material stream. This will create our feed stream. It is a good idea to place your other streams as well: one for distillate product and one for bottoms product.
2. Add two energy streams that will be the condenser and reboiler duties.
3. Double click on the newly created feed stream to pull up the worksheet for the stream. Please note that we are specifying this feed stream to demonstrate how the column would behaved with a fully specified inlet stream. The temperature, pressure, flow rate and composition of the feed stream can be specified by the user.

Figure 4. Define Feed

1. We will make the following changes within this worksheet to specify the feed stream:
2. Enter “feed” into the stream name.
3. Set the temperature to 20 degrees C.
4. Set the pressure to 1 atm. (Be sure to select the “atm” unit from the dropdown when you begin to enter this value)
5. Set the mass flow to 500 kg/hr

Figure 5. Define Compositions

1. On the left hand column within the worksheet, click “composition” to enter the composition tab.
2. Click “edit” on the composition tab to change the composition of the feed stream
3. Enter .25 to Toluene, .35 to Benzene and .40 to p-Xylene.

The worksheet should indicate that the stream is fully specified with a green “ok” indicator at the bottom. Close the worksheet.

Fully Specify Shortcut Column

Figure 6. Connect the inlet and outlet streams

1. Create the shortcut column by clicking on the "columns" tab within the palette and selecting "shortcut column"
2. Click on the simulation environment to place the column
3. The shortcut column should currently appear red to represent that it is not currently specified.
4. Double click on the column to pull up the column worksheet
5. Under each drop down, add the appropriate stream: E.g. Under "Distillate" add the distillate product stream that was defined earlier.
6. Additionally, be sure to specify the quality of the top product stream: whether it is liquid or vapor phase. In this case, select "liquid."

Figure 7. Set the key specifications

1. On the left hand side of the worksheet, select "Parameters."
2. This will move you to the parameters section of the column window where column specifications can be added.
3. At this point it is important to identify which components are most volatile. In this case, benzene is most volatile, xylene is the least volatile and toluene rests in between the two. One rule of thumb is that more volatile components tend to be more valuable, and so in this example we will be purifying benzene in the distillate stream. In order to do this we will specify that benzene is the light key in the bottoms and that toluene is the heavy key in the distillate. This arrangement makes sense since it will force the heavier xylene and toluene products to drop out in the bottoms.
4. Set the Light Key in Bottoms to Benzene and the mole fraction to .001 and set the Heavy Key in Distillate to Toluene and the mole fraction to .001.
5. Set the condenser and reboiler pressure to the approximate values that they would be. In this case, we are assuming that our column is sufficiently small that atmospheric pressure is a fair assumption for both values. Enter 1 atm for both of these values.

Figure 8. Viewing the Worksheet

1. Finally set the reflux ratio to a value that is higher than the minimum reflux ratio that HYSYS calculates. In this case the minimum reflux ratio was calculated to be 2.119 and so it is best to input around 2.5. This value can be tweaked later to see how performance is changed for different reflux ratios.
2. At this point it is worthwhile to see what outputs the column is producing, since this is one of the primary uses of the shortcut distillation column.
3. Click on the "worksheet" tab at the top of the window. This will bring you to your stream information where you can view the temperatures and flow rates of your outgoing streams.
4. Now click on the "compositions" tab on the left side to view the compositions of the outgoing stream. As expected, the benzene is leaving in the distillate stream at 99.9% purity with a .1% toluene impurity, while all of the xylene is leaving out in the bottoms.

Figure 9. View the Performance Tab

1. At this point it is valuable to record the calculated minimum number of trays, actual number of trays, as well as the condenser/reboiler temperatures.
2. To view these values click on the "performance" tab
3. You can observe here that the minimum number of trays calculated was around 15.7, the actual number of trays used was 35.7 the optimal feed stage was around 20.5, the condenser temperature was found to be around 80 C, while the reboiler was calculated to be around 123 C.

Standard Distillation Column

The standard distillation column, found under the "Columns" tab on the equipment palette, may be used to perform a more rigorous simulation of your process than that offered by the shortcut column. In addition, the shortcut column may be used as a "first pass" analysis and then these specification estimations can used as the starting point for the standard column.

Tutorial

In this example, we will be setting up a basic distillation column in HYSYS to simulate the separation of benzene and toluene.

1. Start a new case in HYSYS by clicking on the button on the top menu bar.
2. On the Components screen, add Benzene and Toluene.
3. On the Fluid Package screen, you have a wide range of options with which to simulate vapor-liquid equilibrium. For this example, choose the Peng-Robinson fluid package and enter the simulation environment.

Figure 10. Component Selection in HYSYS

Figure 12. Fluid Package Selection in HYSYS

We are using the simple distillation column and thus we only need a single feed stream into the distillation column as well as a distillate stream and a bottoms stream. We will also need 2 energy streams for the condenser and reboiler.

1. Add the material and energy streams to the simulation environment. Naming them appropriately will make it easer to keep everything straight as we move forward with the tutorial.
2. Select the distillation column from the palette and place anywhere in the simulation environment.
3. Switch to “Attach Mode” and attach the mass and energy streams to their appropriate locations on the distillation column. NOTE: You could also complete this step by double-clicking on the column itself and selecting each respective stream from the drop down menus.
4. Double click on the column and select the condenser type. For this example, select “total condenser” and click next
5. This next page is where the pressures of the column are specified. For this example, we will assume that the column is run at 1 atm with no pressure drops within the column. Therefore, specify reboiler pressure and condenser pressure to 101.3 kPa and the condenser pressure drop to 0 kPa.
6. The next page has fields where temperature estimates can be made for the condenser, top stage, and reboiler. As these are optional fields, simply click next.
7. The final screen of the distillation column setup is where the reflux ratio and distillate flow rates should be specified. For this example, specify the distillate flow rate as 1000 kgmol/hr and the reflux ratio as 1.5.
8. Click “Done” at the bottom of the screen and the simulation page should be back in full view.

Figure 11. Simulation View of Distillation Column in HYSYS

We are now ready to specify the flow rates and compositions of the streams coming in and out of the column.

1. Double click on the feed stream coming into the column. A specifications page should come up.
2. This page allows you specify a number of parameters pertaining to the feed stream, including temperature, pressure, and compositions. Specify the temperature as 37C and pressure as 101.3 kPa (1 atm). In addition, specify the molar flow as 3500 kgmole/h.
3. Now click on the composition tab on the left side of the window. Here specify, the benzene as 70 mass % and the toluene as 30 mass %. The stream should be fully specified. Close the window to return to the simulation.

Figure 13. Column Configuration Page

The column and feed stream should now be fully specified. It is now time to run the simulation.

1. If the solver is not already activated, click the green traffic light at the top of the simulation window.
2. Now double click on the distillation column to bring up the main column window. A feed stage and total number of stages should already be specified.
3. Click the “Run” button of the window to start the simulation. The small rectangle to the right of this button should light up green and say “converged” when it finds a solution.

Figure 14. Outlet compositions in HYSYS

Manipulating Column Specifications

Information taken from Towler Ch. 17

Parameters that directly impact the design and operation of the distillation column include:

• Reflux ratio
• Feed point location
• Column Pressure
• Number of equilibrium stages/real stages

Assuming that we have already specified a desired specification for either the distillate or bottoms stream, we can manipulate the above parameters when running our simulation within HYSYS.

Reflux Ratio

The number of stages required for separation is directly affected by the reflux ratio that can be defined as follows:

$R={\frac {Flow\ returned\ as\ reflux}{Flow\ of\ top\ product\ taken\ off}}$

Increasing the reflux ratio will reduce the total number of stages required but will increase the utility costs for the column. Therefore, when simulating a distillation column, special care should be taken to strike a balance between capital and operating costs.

When starting a new distillation column simulation within HYSYS, it is a good idea to calculate a minimum reflux ratio to advise the specification of the reflux ratio. The minimum reflux ratio can be calculated through use of a McCabe-Thiele diagram. According to Towler, a reflux ratio of 1.15 times the minimum reflux ratio should be used as a first approximation. Once this ratio is inputted into the simulation environment, it can be further manipulated to fine-tune the optimum reflux ratio for the system.

Feed Point Location

The column simulation environment within HYSYS allows for manipulation of the feed point location. The general rule for designating feed point location is to have the feed enter the column at the stage which best matches the composition of the feed. There can be some uncertainty when it comes to finding the precise feed point location. Instead, a trial-and-error method may be useful when using HYSYS to determine which feed stage results in the desired distillate and bottoms conditions.

Column Pressure

When specifying the distillation column within the HYSYS simulation environment, you must specify the operating pressure as well as any pressure drops. If you are running a HYSYS simulation to get a first-pass approximation of compositions and flow rates, it is a typical practice to assume the pressure in the column is constant. If you are running a more intensive and rigorous simulation, you can approximate the pressure drop per tray using the following equation:

$\Delta P_{\text{tray}}=2\rho _{L}gh_{w}$

where $\rho _{L}$ is the liquid density, g is gravitational acceleration, and $h_{w}$ is the weir height.