2

The Optimizer

Workshop
Fast track to page 9.

In this example, a simple distillation column to separate Tetrahydrofuran (THF) from Toluene is simulated. The object of the exercise is to select the product specifications such that profit is maximized. A special tool in HYSYS, the Optimizer, will be used to find the optimum operating conditions.

Learning Objectives
Once you have completed this section, you will be able to:

Use the Optimizer tool in HYSYS to optimize flowsheets Use the Spreadsheet to perform calculations

Prerequisites
Before beginning this section you need to be able to:

Add Streams and Operations Model columns in HYSYS

Process Overview

Column Overview

The Optimizer

Building the Simulation

Defining the Simulation Basis
For this case, you will be using the Wilson Activity Model. The components are: Tetrahydrofuran and Toluene. 1. On the Binary Coeffs tab, check that the UNIFAC VLE estimation method is chosen, and press the Unknowns Only button to estimate the missing Interaction Parameters. Enter the Simulation Environment.

2.

Binary interaction parameters are used to correlate lab data with a thermodynamic model. When lab data is not available, you can estimate the parameters with HYSYS UNIFAC estimation method.

Adding the Feed Stream

Add a material stream with the following values: In This Cell...
Conditions Name Temperature Pressure Mass Flow Composition - Mass Fraction THF Toluene 0.44 0.56 Feed 10C (50F) 140 kPa (20 psia) 3700 kg/h (8160 lb/hr)

Enter...

The Optimizer

Add the Distillation Column

Add a Distillation Column with the following values: In This Cell...
Connections Column Name No. of Stages Feed Condenser Type Bottoms Liquid Ovhd Liquid Condenser Energy Stream Reboiler Energy Stream Pressure Delta P, Condenser and Reboiler Condenser Reboiler 0 kPa 103 kPa (15 psia) 107 kPa (15.5 psia) T-100 10 Feed, Stage 5 Total Toluene THF Cond Q Reb Q

Enter...

Make the following specifications: In This Cell...

Specs Reflux Ratio Distillate Rate 2 (Estimate) 1500 kg/h (3305 lb/hr) (Estimate)

Enter...

The Optimizer

Create two new Column Component Fraction specifications as shown below:

These two specifications should be Active. The DOF for the column should now be 0. The column should now solve automatically. If it does not, press the Run button to start the solver.

Save your case!

The Optimizer

Changing the Column Tolerances

In order for the optimizer to work properly, we have to tighten the tolerances of the column solver. There are two tolerances that the column must meet before it can be considered as solved.

The Equilibrium Error Tolerance The Heat/Spec Error Tolerance

We want to set the value for both of these tolerances at 1e-6. Note that this will increase the solving time of the column, but this column solves very fast anyway, so the tighter tolerances are acceptable in this case. To change the tolerances, follow these steps: 1. 2. Access the Solver page of the Parameters tab. Enter 1e-6 for both the tolerance values. These cells are located in the Solving Options group box.

This completes the changing of the tolerances. We are now ready to begin to optimize the column.

The Optimizer

Adding the Optimizer

In todays highly competitive market place, with stringent plant emissions controls and increased competition, a plants survivability is often determined by its ability to remain competitive. Optimization methods are now available that give the process engineer the necessary tools to perform on-going process improvement studies. Optimization studies lead directly to improved plant performance, efficient plant operation and finally to increased profitability. Typically, optimization studies involve an economic model composed of a Profit Function and operating constraints. HYSYS contains a multi-variable Steady State Optimizer. Once your Flowsheet has been built and a converged solution has been obtained, you can use the Optimizer to find the operating conditions which minimize or maximize an Objective Function. The Optimizer owns its own Spreadsheet for defining the Objective Functions as well as any constraint expressions to be used. This allows you to construct Objective Functions which maximize profit, minimize utilities or minimize exchanger UA.

Open the case OptOptimizer.hsc

Only user-specified process variables can be used as Primary Variables

Restrictions on the Optimizer only available for SteadyState calculations it cannot be used in Templates.

Primary Variables - these are flowsheet variables whose values are manipulated in order to minimize (or maximize) the objective function. You set the upper and lower bounds for the primary variables, which are used to set the search range. Objective Function - this is the function which is to be minimized or maximized. The function has to be defined within the Spreadsheet. This allows the user a great deal of flexibility in defining the function. Constraint Functions - inequality and equality functions are defined in the Spreadsheet. In solving the Objective Function, the Optimizer must also meet any constraints that are defined by the user.

10

The Optimizer

In our column example, the Profit Function/Cost Function (sometimes referred to as the Objective Function, Performance Criterion or Performance Index) requires calculating a net profit for the column. The Profit Function is a function of the revenue generated from desired products THF and Toluene, within the limit of product purity constraints. While increased profits is directly linked to increased production of the desired product, plant profitability is generally offset by operating costs. In our column, operating costs are associated with the column utility requirements for the reboiler and condenser. In general a profit function is given by: PF =

F pMp FfMf OC
p f

where: PF = Profit Function/time = Product Flows * Product Values = Feedstock Flows * Feedstock Costs = Operating Costs

FpMp FfMf
OC

Profit = (THF Product + Toluene Product) Heating Cost Cooling Cost Feed Cost Typically, the economic model includes operating constraints (equality or inequality constraints or equations). The operating constraints comprise the model of the process or equipment. In this case, we have no process constraints.

Which variable can we change to affect the Revenue?________ To invoke the Optimizer, select Optimizer under Simulation in the Menu Bar, or press <F5>.

10

The Optimizer

11

Variables tab
When you invoke the Optimizer for the first time, you are placed on the Variables tab. On the Variables tab you define the Adjusted (Primary) Variables to be used in the optimization. In this case, our Primary Variables are the purity of our products, THF and Toluene. 1. Press the Add button to add the first variable, T-100, Spec Value, THF Purity Spec.

2. 3.

Set the Upper Bound at 0.9, and the Lower Bound at 0.90. Add the second variable, T-100, Spec Value, Toluene Purity Spec, with the same bounds as above.

11

12

The Optimizer

Defining the Spreadsheet

The Optimizer has its own Spreadsheet for defining the Objective and Constraint functions. Primary Variables may be imported and functions defined within the Optimizer Spreadsheet, which possesses the same capabilities as the Main Flowsheet Spreadsheet 1. 2. Press the Spreadsheet button on the Optimizer view to open the Spreadsheet. Move to the Spreadsheet tab.

The Spreadsheet is an operation and thus the Spreadsheet cells get updated when Flowsheet variables change

Importing and Exporting Variables

You may import virtually any variable in the simulation into the Spreadsheet and you can export a cells value to any specific field in your simulation.

Object Inspection - object inspect (secondary mouse button) the cell which you want to Import into, or Export from. From the Menu that appears, select Import Variable or Export Formula Result. Then, using the Variable Navigator, select the variable you wish to import or export. Connections page tab - select the Add Import or Add Export button. Then, using the Variable Navigator, select the variable you wish to import or export. Drag n Drop - using the secondary mouse button, click the variable value (from the WorkBook or Property View) you wish to import, and drag it to the desired location in the Spreadsheet. If you are exporting the variable, drag it from the Spreadsheet to the exported location.

12

The Optimizer

13

Adding Formulas
Complex mathematical formulas can be created, using syntax that is similar to conventional Spreadsheets. Arithmetic, logarithmic and trigonometric functions can be performed in the Spreadsheet. All common functions must be preceded by a + symbol. Special Functions must be preceded by the @ symbol. Some of the functions available are:

Press the Function Help button to view the Available Spreadsheet Functions and Expressions.

Addition (+): +A1+A2 Subtraction (-): +A1-A2 Multiplication (*): +A1*A2 Division (/): +A1/A2 Power (^): +A1^3 Absolute Value (@ABS):@ABS(A1) Square Root (@SQRT):@SQRT(A1) Natural Log (@ln):@ln(A1) Exponential (@exp):@exp(A1)

The following variables need to be imported into the Spreadsheet. Text entries are added to the spreadsheet by typing them in the appropriate cell. Cell...
B1 B2 B3 B4 B5 B6

Object...
Cond Q Reb Q THF THF Toluene Toluene

Variable...
Heat Flow Heat Flow Mass Flow Comp Mass Frac THF Mass Flow Comp Mass Frac Toluene Mass Flow

B7

Feed

The quickest way of importing variables is to right-click on the desired cell. Select Import Variable, and use the variable navigator to locate the desired variable. Use this method to import the variables above now.

13

14

The Optimizer

Enter the following constants on the spreadsheet. The comments can be added, if desired. Cell
D1 D2 D7

Value (SI Units)

0.471 ($/kWh) 0.737 ($/kWh) 0.05 ($/kg)

Value (Field)
138 ($/MMBtu) 216 ($/MMBtu) 0.024 ($/lb)

Comment
Cooling Cost Heating Cost Feed Cost

Enter the following product prices: Equation (SI Units)

+0.333*b4^3 ($/kg)

Cell
D4

Equation (Field)
+0.151*b4^3 ($/lb)

Description
THF Price (corrected for purity) Toluene Price (corrected for purity)

D6

+0.163*b6^3 ($/kg)

+0.074*b6^3 ($/lb)

Note: the prices of the products decrease as the impurities increase. The Objective Function is placed in Cell D8. The equation is:
You can change the Variable Type to Unitless for dollar value variables.

+(b3*d4+b5*d6)-b7*d7-(b1*d1+b2*d2)/3600 Pay special attention to the units in this equation. If using field units, you may have to add an additional term to this equation to convert between hours and days. What is the value of cell D8, the Profit? __________

14

The Optimizer

15

Functions tab
The Functions tab contains two groups, the Objective Function and the Constraint Functions. However, in this example we do not have constraint functions. 1. In the Cell area of the Objective Function group, specify the Spreadsheet cell that defines the Objective Function. Use the drop down menu in the Edit Bar to select the appropriate cell. The Current Value of the Objective Function will be provided. Select the Maximize radio button.

2.

15

16

The Optimizer

Parameters tab
The Parameters tab is used for selecting the Optimization Scheme.

For more information on the Optimization Schemes, refer to the manual section 17.2 or the on-line Help.

Box - Handles inequality constraints but not equality constraints. It generally requires a large number of iterations to converge on the solution. SQP - Sequential Quadratic Programming, handles inequality and equality constraints. Considered by many to be the most efficient method for minimization. Mixed - Handles inequality constraints only. It is a combination of the Box and SQP methods. It starts the minimization with the Box method using a very loose convergence tolerance. After convergence, the SQP method is used to locate the final solution. Fletcher Reeves - Does not handle constraints. Efficient method for general minimization. Quasi-Newton - Does not handle constraints. Similar method to Fletcher Reeves.
1. 2. Select the Mixed method as the Scheme. Use the defaults for Tolerance and Number of Iterations.

Monitor tab
The Monitor tab displays the values of the Objective Function and Primary Variables during the Optimizer calculations. 1. Move to the Monitor tab and press the Start button to begin the optimization.

Save your case!

16

The Optimizer

17

Analysing the Results

Once the Optimization is complete, examine the results and fill in the following table: Base Case
THF mass flow 1650.8 kg/h (3639.4 lb/hr) 0.95 2049.2 kg/h (4517.7 lb/hr) 0.95 9.926e5 kJ/h (9.408e5 Btu/hr) 1.5980e6 kJ/h (1.5146e6 Btu/hr) 106.3 $/hr

Optimized Case

THF purity Toluene mass flow

Toluene purity Cond duty

Reb duty

Profit

17

18

The Optimizer

Exercise 1
We are going to introduce a constraint on the liquid volume flow of the stream THF and examine how it affects the results, profit and products purities. First, the case must be set back to its state before the optimization. 1. 2. 3. Go to the column subflowsheet. On the Design tab, Monitor page, enter the initial value of 0.95 for the THF and Toluene purity specifications. Re-run the column.

The constraint on the liquid flow is that it must not exceed 1.85 m3/h (65.5 ft3/hr). The THF Liquid Flow must first be imported into the spreadsheet (use an empty cell), and the constraint value must be written in the spreadsheet. Remember to write the comments next to the values so that they can be understood. 1. 2. 3. On the Functions tab of the Optimizer press the Add button. Type in the reference of the cell where the variable constraint is located. In the Cond column, use the scroll down arrow to find the less than sign. Type in the reference of the cell where the value of the constraint is written. Go to the Monitor tab and start the Optimizer.

4.
The constraint values are positive if inequality constraints are satisfied and negative if inequality constraints are not satisfied.

18

The Optimizer

19

What differences, if any, does the constraint make to the first optimized solution? Fill in the following table: Optimized Case Optimized Case 2

Base Case
THF Mass Flow THF Purity Toluene Mass Flow Toluene Purity Cond Duty Reb Duty Profit 1650.8 kg/h 0.95 2049.2 kg/h

0.95 9.926e5 kJ/h 1.5980e6 kJ/h 106.3 $

19

20

The Optimizer

Exercise 2
Here, we are going to model the reboiler with a steam-heated shell and tube heat exchanger. The heat exchanger will be modelled with 115 psia steam and the maximum flow of steam available to the reboiler is limited to 840 kg/h (1850 lb/hr). (Note that there are a few approaches that can be taken with the heat exchanger modelling so two of the possibilities are listed). The cost of the 115 psia steam is 0.682 $/kWh ($200/MMBtu).
Use the Mixed optimization scheme for this exercise.

1. 2. 3.

Add Water as a component. This must be done in the Basis Environment. Open the Columns Property viewer; on the Monitor page, reset the THF and Toluene purity specifications to 0.95. Add a new internal stream (on the Flowsheet tab), ToReb with the following attributes: Enter...
Liquid Total 10 Yes

In This Cell...
Type Net/Total Stage Export

4. 5.

Run the column. Return to the Main Environment (notice that the "To Reb" stream is now shown on the PFD. Add a heat exchanger. The stream ToReb is the shell side feed and steam is on the tube side. Set the Shell and Tube side pressure drops to 0. Specify the Steam inlet conditions at Vf = 1.0; P = 790 kPa (115 psia); Flow = 770 kg/h (1700 lb/hr). The outlet steam is at its bubble point (Vf = 0.0). Add a Duty Spec to the Heat Exchanger but do not specify a value (pass is Overall). What is the Heat Exchanger Duty? __________

6. 7.

8.

20

The Optimizer

21

9.

Add an Adjust operation. The Adjusted Variable is the Steam flow rate and the Target Variable is the Spec Calc Value (Duty Spec).

10.

Export the RebQ Heat Flow value (you will have to make this Exportable first) from the Optimizer spreadsheet into the Target value for the Adjust operation. A cell on the spreadsheet can not be simultaneously imported and exported. Copy the value of cell B2 to another empty cell by entering the formula "+B2". The value in this new cell can then be exported to the Target Value of the Adjust operation. Change the Heating Cost value to 0.682 $/kWh ($200/MMBtu) (Cell D2 on the spreadsheet) and Start the Adjust. What is the Exchanger Duty and what is the Steam mass flow rate? __________ & __________

11.

12. 13.

Add the Steam Mass Flow (must be less than 840 kg/h (1850 lb/ hr)) constraint to the Optimizer. Start the Optimizer.

21