Hysys 2020 PDF
Hysys 2020 PDF
Hysys 2020 PDF
Aspen HYSYS
Aspen HYSYS is a software program used for simulation of chemical process. It is a market-
leading process modeling tool for conceptual design, optimization, business planning, asset
management, and performance monitoring. Due to its accuracy it is widely used in industries such as
oil & gas production, gas processing, petroleum refining, and air separation industries. Study shows
the advantages of hysys and why it is preferred over other softwares. These lectures give
introduction of working of hysys in various chemical fields. Aspen HYSYS is a core element of
AspenTech’s aspenONE Engineering applications. Aspen HYSYS has established itself as a very
intuitive and easy to use process simulator in oil and gas refining industry. Users with little prior
knowledge of Aspen HYSYS can pick up and train themselves in its modelling capabilities. Some of
the very intuitive capabilities include a highly interactive process flow diagram for building and
navigating through large simulations. The program also provides a very flexible and easy to use
distillation column modelling environment. Additionally the interactive nature of HYSYS enables users
to build and use their models quickly and effectively. Aspen HYSYS offers a comprehensive
thermodynamics foundation for accurate calculation of physical properties, transport properties, and
phase behaviour for the oil & gas and refining industries. Comprehensive library of unit operation
models including distillation, reactors, heat transfer operation, rotating equipments, controllers and
logical operations in both the steady state and dynamic environments.
Simulation is an situation in which a particular set of conditions is created artificially in order to study
or experience something that could really exist in reality. It is the act if pretending that something is
real when it is not. A computer simulation is an attempt to model a real-life or hypothetical
situation on a computer so that it can be studied to see how the system works.
There are two modes of simulation Steady state mode and Dynamic mode.
Steady State Mode: Initially process simulation was used to simulate steady state processes. Steady
state models perform a mass and energy balance of a stationary process (a Process in an equilibrium
state) but any changes over time had to be ignored.
Dynamics Mode: Dynamic simulations require increased calculation time and are mathematically more
complex than a steady state simulation. It can be seen as a multiply repeated steady state simulation
(based on a fixed time step) with constantly changing parameters.
1
Chapter 1
This chapter begins by starting Aspen HYSYS and how to select the right components and fluid
package for simulation purposes. Knowing how to start Aspen HYSYS and get familiar with its desktop is
very important in this chapter. The second part is about how to enter and reenter the simulation
environment, and get familiar with simulation flowsheet. In this part, users will be informed some
important features of Aspen HYSYS. The last part is dealing with how to add and specify material streams
for simulation. Variables specification is one of the important steps that users need to understand when
dealing with Aspen HYSYS.
Learning Outcomes: At the end of this chapter, the user will be able to:
• Select Components
Mostly in this manual, and for reduction it will named Aspen Hysys with (Hysys)
2
3
Edit Menu
4
Hot Keys
5
6
1.1 Starting Aspen HYSYS
From “Start Menu” Select: All Program/Aspen Tech/ Process Modeling/ Aspen HYSYS/ Aspen HYSYS
V8.4
Figure 1-1
Before any simulation can occur, Aspen HYSYS needs to undergo an initial setup. During an initial setup,
the components and the fluids package that will be used will be selected.
Aspen HYSYS used the concept of the fluid package to contain all necessary information for
performing flash and physical property calculations. This approach allows you to define all information
(property package, components, hypothetical components, interaction parameters, reactions, tabular data,
etc.) inside a single entity.
• All associated information is defined in a single location, allowing for easy creation and modification of
the information.
• Fluid packages can be stored as completely defined entities for use in any simulation.
• Component lists can be stored out separately from the Fluid Packages as completely defined entities for
use in any simulation.
• Multiple Fluid Packages can be used in the same simulation. However, they are defined inside the
common Basis Manager.
The Simulation Basis Manager is property view that allows you to create and manipulate multiple fluid
packages or component lists in the simulation.
7
1.3 Creating A New Simulation
Select File/New/Case or press Crtl+N or click on the New Case to start a new case. In HYSYS,
your simulation is referred to as a “case”. This will open up the Simulation Basis Manager which is where
all of the components and their properties can be specified.
Before proceeding any further, save your file in an appropriate location. Select File/Save As and select
where to save the file. Do not save the file to the default location.
The first step in establishing the simulation basis is to set the chemical components which will be
present in your simulation.
1. To add components to the simulation, click on the Add button in the Simulation Basis Manager as in
Figure 1-2.
2. Clicking on Add (Figure 1-2) will bring up the Component List View which is a list of all the
components available in Aspen HYSYS.
Navigation pane
Add Component
Figure 1-2
8
3. Select the desired components for your simulation. You can search through the list of components in
one of three ways, Figure 1-3:
a. Simulation Name
b. Full Name
c. Formula
Type of Search
Figure 1-3
4. Once you have located the desired component, either double click on the component or click < Add
button to add it to the list of components for the simulation.
5. At the bottom of the components page, you can give your component list a name.
6. Once this is complete, click on under the Component List-1 OR click on the
in the Home Ribbon to select proper fluid package.
The fluid package is used to calculate the fluid/thermodynamic properties of the components and
mixtures in your simulation (such as enthalpy, entropy, density, vapour-liquid equilibrium … etc.).
Therefore, it is very important that you select the correct fluid package since this forms the basis for the
results returned by your simulation.
1. Select the Fluid Package by clicking on Add button as shown in figure 1-4.
9
Click for help of select the Fluid package
Figure 1-4
3. From the list of Property package Selection, select the desired Fluid package.
4. Once the desired model has been located, select it by clicking on it once (no need to double click). For
example, select Peng-Robinson property package for your simulation.
When faced with choosing a thermodynamic model, it is helpful to at least a logical procedure for
deciding which model to try first.
Elliott and Lira (1999) suggested a decision tree as shown in Figure 1-5.
The property packages available in HYSYS allow you to predict properties of mixtures ranging from
well defined light hydrocarbon systems to complex oil mixtures and highly non-ideal (non-electrolyte)
chemical systems. HYSYS provides enhanced equations of state (PR and PRSV) for rigorous treatment of
hydrocarbon systems; semiempirical and vapor pressure models for the heavier hydrocarbon systems;
steam correlations for accurate steam property predictions; and activity coefficient models for chemical
systems. All of these equations have their own inherent limitations and you are encouraged to become
more familiar with the application of each equation.
10
Figure 1-5
11
The following table lists some typical systems and recommended correlations.
For oil, gas and petrochemical applications, the Peng-Robinson EOS (PR) is generally the recommended
property package. For more details, please refer Aspen HYSYS Simulation Basis Manual.
If You can’t select the proper fluid package, it can be used Methods Assistant in the
Home Ribbon, then click next , then select Specify component type or select Specify process type. If
you choose Specify component type, according to the our example select Hydrocarbon system , from the
message appears, click No, the program tell you to Use Peng-Robenson, Lee-Kesler-Polcker or SRK
12
1.7 Enter Simulation Environment
You have now completed all necessary input to begin your simulation. Click on the Simulation button in
the left bottom , Figure 1-3.
Once you have specified the components and fluid package, and entered the simulation environment,
you will see the view as shown in Figure 1-6. Before proceeding, you should taking care of a few features
of this simulation window:
1. HYSYS, unlike the majority of other simulation packages, solves the flowsheet after each
addition/change to the flowsheet. This feature can be disabled by clicking the Solver Holding button (the
red light button ) located in the Ribbon of the Home menu (see Figure 1-6). If this button is
selected, then HYSYS will not solve the simulation and it will not provide any results. In order to allow
HYSYS to return results, the Solver Active button (the green light button ) must be selected.
2. Unlike most other process simulators, HYSYS is capable of solving for information both downstream
and upstream. Therefore, it is very important to pay close attention to your flowsheet specification to
ensure that you are not providing HYSYS with conflicting information. Otherwise, you will get an error
and the simulation will not solve.
Object Palette
Figure 1-6
13
Re-Entering the Simulation Basis Manager
When the basis of the simulation has to be changed, the Simulation Basis Manager needs to be re-
Object Palette
On the left hand side of Figure 1-6, you will notice a vertical toolbar. This is known as the Object Palette.
If for any reason this palette is not visible, got to the Flowsheet/Modify Ribbon and select or
press F4 to display the palette. From this toolbar that you will add streams and unit operations to your
simulation.
Material Streams are used to transport the material components from process units in the simulation. A
material stream can be added to the flowsheet in one of two ways:
1. Click on the blue arrow button on the Object Palette, Figure 1-6, Or
2. Pressing F11
Using any of the above methods will create a new material stream (a Blue arrow) on the flowsheet,
refer Figure 1-7. The HYSYS default names the stream in increasing numerical order (i.e. the first stream
created will be given the name “1”). This name can be modified at any time.
To enter information about the material stream, double click on the stream to show the window shown
in Figure 1-7. It is within this window that the user specifies the details regarding the material stream. For
material stream that will be used as an input, we need to specify four variables. Within HYSYS
environment, input material stream always have four degree of freedoms. Meaning, we need to supply
four information in order to fulfill the requirement for HYSYS to start its calculations.
14
Figure 1-7
From Figure 1-7, you will see the warning yellow message bar at the bottom of the window indicating
what information is needed (unknown compositions). Just follow what the message wants, for example,
the first thing that you need to supply is compositions. In order to specify the composition of the stream,
select the “Composition” option from this list to display the window in Figure 1-8. It is within this
window that the user specifies the composition of the stream. Note that only the components that you
specified in the simulation basis manager will appear in this list. You can specify the composition in many
different ways by clicking on the “Basis…” button. The HYSYS default is mole fractions, however the
user can also specify mass fractions, liquid volume fractions, or flows of each component. If the user is
specifying fractions, all fractions must add up to 1. Enter mole fraction of 1 in the methane section to
indicate 1 mole fraction of methane.
15
Tabs
pages
Figure 1-8
Next, the warning yellow message bar indicates that you need to specify the input temperature for this
stream. In order to specify the temperature of the stream, select the “Conditions” option from this list to
display the window in Figure 1-9. It is within this window that the user specifies the temperature of the
stream.
When entering the conditions for a stream, it is not necessary to enter the values in the default units
provided. When the user begins to enter a value in one of the cells, a drop down arrow appears in the units
box next to the cell. By clicking on this drop down arrow, the user can specify any unit for the
corresponding value and HYSYS will automatically convert the value to the default unit set. Enter the
temperature of 25 in the temperature section to indicate the temperature of 25oC. Next, the yellow warning
message bar indicates that you need to specify the input pressure for this stream. In the same window,
enter the pressure of 10 in the pressure section to indicate the pressure of 10 bar_g as shown in
Figure 1-10 , (Pabsolute = Pg + Patm).
16
Figure 1-9
Figure 1-10
17
Next, the last variable that you need to specify is flowrate. For this, you have two options either to
specify molar flowrate or mass flowrate. In the same window, enter the molar flowrate of 100 in the
molar flowrate section to indicate the flowrate of 100 kgmole/h as shown in Figure 1-10.
Once all of the stream information has been entered, HYSYS will calculate the remaining properties and
data provided it has enough information from the rest of the flowsheet. Once a stream has enough
information to be completely characterized, a green message bar appears at the bottom of the window
within the stream input view indicating that everything is “OK” (See Figure 1-10). Otherwise, the input
window will have a yellow message bar at the bottom of the window indicating what information is
missing.
Values shown in blue have been specified by the user and can be modified while values shown in black
have been calculated by HYSYS and can not be modified.
The following color code for material streams on the flowhseet indicates whether HYSYS has enough
information to completely characterize the stream:
In the first part of this chapter, we opened it with how to start HYSYS and get familiar with its desktop
environment. We also discussed how to select components that will be used in simulation. Selecting the
right fluid/thermodynamic package is very important and therefore we provided a flowchart that will assist
users to select the right thermodynamics models.
The second part of this chapter was about how to enter and re-enter the simulation environment, and get
familiar with simulation flowsheet. In this part, users are also informed some important features of
HYSYS.
The last part of this chapter was dealing with how to add and specify material streams for simulation.
Variables specification is one of the important steps that users need to understand when dealing with
HYSYS. When users wanted to specify streams especially materials, they need to specify at least four
variables in order to have HYSYS to calculate the remaining properties.
1.10 Problems
1.1. Create one materials stream that contains only water with following conditions:
Fluid Package: Peng-Robinson, Flowrate: 100 kgmole/h, Pressure: 1 atm, Vapor/Phase Fraction: 1.0
18
1.2. Repeat the above procedures by replacing pressure with temperature of 150oC.
1.3. With the same condition in (2), reduce the temperature to 70oC.
1.4. Create one new materials stream that contains only water with following conditions:
Fluid Package: Peng-Robinson, Flowrate: 100 kgmole/h, Pressure: 2 atm, Vapor/Phase Fraction:1.00
1.5. With the same condition in (4), increase the pressure to 5 atm.
1.6. With the same condition in (4), decrease the pressure to 0.5 atm.
____________________________________________________________________
19
Chapter 2
Equation of State
Problem: Find the specific volume of n-butane at 500 K and 18 atm using the following equation of state:
• Soave-Redlich-Kwong (SRK)
• Peng-Robinson (PR)
The first step in building any simulation is defining the fluid package. A brief review on how to define a
fluid package and install streams is described below. For a complete description, see the previous chapter
(Chapter 1: Starting with HYSYS).
Components: n-butane
2. Click the Simulation button when you are ready to start building the simulation.
Installing a Stream
There are several ways to create streams. (For complete description, see the previous chapter.)
Name 1
Temperature 500 K
Pressure 18 atm
Composition 1
20
Saving
File menu/ Select Save As/ Give the HYSYS file the name EOS SRK then press the OK button.
Figure 2-1
1. In the Home Ribbon select Workbook , from Navigation pane; click on Workbook, or click
Ctrl+W as shown in Figure 2-2.
Figure 2-2
21
2. Specific volume in HYSYS is defined as Molar Volume. From Figure 2-2, there is no Molar Volume
shown in the Workbook. In order to preview the value of Molar Volume, we have to add it to the
Workbook.
3. To add the Molar Volume or other variables, click on the icon on the left top. The setup window for
Workbook can be viewed as shown in Figure 2-3.
Figure 2-3
4. Under the Variables group, click the Add button at the right side of the window.
5. Window for you to select variables will appear as shown in Figure 2-4.
22
Figure 2-4
6. In the Variable column, scroll down until you find the Molar Volume and then click OK.
7. The Molar Volume value is presented in the Workbook as shown in Figure 2-5.
Figure 2-5
In this section, we will analyze the specific volume of n-butane when the temperature is changing. To
achieve this analysis, do the following:
1. In the Navigation pane click on Case Studies or In the Home menu click on icon as
shown in Figure 2-6.
23
OR
Either
Figure 2-6
2. Next, click Add button, then Set up view displays as shown in Figure 2-7.
Figure 2-7
3. Click Add button, the Variable Navigator will show as in the figure 2-8. In the Object column select
stream 1, and in the Variable column, select Molar Volume. Then, click OK button.
24
Figure 2-8
4. Repeat step 3 to insert Temperature. The new updated Case Study is shown in Figure 2-9.
Figure 2-9
Figure 2-9 shows the Temperature is the Independent Variable and Molar Volume is the dependent
variable.
5. Give the Low Bound, High Bound and Step Size of the Temperature (oC) as shown in the figure 2-10.
25
Figure 2-10
7. Click the Run button to Run the analysis. Once the analysis finished, click Plots to view the result as
shown in figure 2-11.
Figure 2-11
26
The results can be viewed as a Table by click on the Results button
1. Press the Properties button which is located left side, click on the Basis-1 folder in the
Navigation Pane.
2. This should take you to the Fluid Package window. In the Property package selection column, select
Peng Robinson, a message will appear, click No
3. Press the Simulation button located left side to return to the PFD.
Compare the result using two different fluid packages; Soave-Redlich-Kwong and Peng-Robinson.
___________________________________________________________________________
___________________________________________________________________________
You have solved a very simple problem to find the specific volume of a pure component sing Aspen
HYSYS. When you use Aspen HYSYS, the parameters are stored in a database, and the calculations are
pre-programmed. Your main concern is to use the graphical user interface (GUI) correctly.
In this chapter, you are able to preview the result using Workbook. Workbook is the most concise way to
display process information in a tabular format. The Workbook is designed for this purpose and extends
the concept to the entire simulation. In addition to displaying stream and general unit operation
information, the Workbook is also configured to display information about any object type (streams,
pipes, controllers, separators, etc.).
You are also should be to analyze the process property using Case Studies. The Case Study is used to
monitor the response of key process variables to changes in your steady state process. After the Case
Study solves, you can view the results in a plot.
27
Lastly, you are able to compare the result from two different equation of state, Peng-Robinson (PR) and
Soave-Redlich-Kwong (SRK).
2.2 Problems
1. Find the molar volume of ammonia gas at 56 atm and 450 K using Soave-Redlich-Kwong (SRK)
equation of state.
2. Find the molar volume of methanol gas at 100 atm and 300oC using Peng-Robinson (PR) equation of
state. Compare its molar volume when you are using Soave-Redlich-Kwong (SRK) equation of state.
3. Consider the following mixture going into a Water-Gas-Shift reactor to make hydrogen for the
hydrogen economy. CO, 630 kmol/h; H2O, 1130 kmol/h; CO2, 189 kmol/h; H2, 63 kmol/h. The gas is at
1 atm and 500 K. Compute the specific volume of this mixture using Soave-Redlich-Kwong (SRK)
equation of state.
4. Consider a mixture of 25 percent ammonia, and the rest nitrogen and hydrogen in a 1:3 ratio. The gas is
at 270 atm and 550 K. Use Peng-Robinson equation of state to compute the specific volume of this
mixture.
5. Consider the following mixture that is coming out of a methanol reactor: CO, 100 kmol/h; H2, 200
kmol/h; methanol, 100 kmol/h. The gas is at 100 atm and 300oC. Compute the specific volume using
Soave-Redlich-Kwong (SRK) equation of state and compare it with Peng-Robinson (PR) equation of
state.
28
Chapter 3
Pump
This chapter begins with a problem to find the pump outlet temperature when given the pump
efficiency. The user will operate a pump operation in HYSYS to model the pumping process. The user
will learn how to connect streams to unit operations such as pump. At the end of this chapter, the user will
determine the pump outlet temperature when given pump efficiency or vice versa.
The Pump operation is used to increase the pressure of an inlet liquid stream. Depending on the
information specified, the Pump calculates either an unknown pressure, temperature or pump efficiency.
Pumps are used to move liquids. The pump increases the pressure of the liquid. Water at 120oC and 3
bar is fed into a pump that has only 10% efficiency. The flowrate of the water is 100 kgmole/h and its
outlet pressure from the pump is 84 bar. Using Peng-Robinson equation of state as a fluid package,
determine the outlet temperature of the water.
Components: Water
3. Click the Simulation button when you are ready to start building the simulation.
Name Feed
Temperature 120oC
Pressure 3 bar
29
3.4 Adding Unit Operations
- Object Palette: From Object Palette; Double click on the Icon of the operation you want to add
- Workbook: Open Workbook (in the Home Menu or from Navigation Pane) and go to the UnitOps
page, then click the Add UnitOp button. The UnitOps view appears
1. From Figure 3-1, double-click on the Pump P-100 icon to open the pump window as shown in
Figure 3-2.
30
Figure 3-2: Pump Window Property
2. In the Inlet, scroll down to select Feed and in the Outlet select Product as shown in Figure 3-3.
3. From Figure 3-3, the warning red message bar at the bottom of the window indicating that we need an
energy stream. To create an energy stream for the pump, click to the space in the Energy, and type work.
This will create energy stream name work for the pump as shown in Figure 3-4.
31
Figure 3-4: Creating Energy Stream for Pump
4. Once a pump has enough information, a green bar appears at the bottom of the window indicating that
everything is “OK” (See Figure 3-4).
Default efficiency for the pump is 75%. "The efficiency of the pump is lower than 0.5"
3. In the Adiabatic Efficiency box on the parameters page, enter 10. The units should be in per cent as
shown in Figure 3-5.
32
Figure 3-5: Changing Pump Efficiency
4. After the efficiency is entered the streams of the pumps should be solved. Click on the Worksheet tab
to view the results.
3.7 Saving: Go to the File menu, Select Save As, Give the HYSYS file the name Pump then press the OK
button.
3.8 Discussion
This example shows that pumping liquid can increase their temperature. In this case, the pump was only
10% efficient and it caused 18°C in the temperature of the water. The less efficient a pump is, the greater
the increase in the temperature of the fluid being pumped. This arises because in a low efficient pump,
more energy is needed to pump the liquid to get the same outlet pressure of a more efficient pump. So the
extra energy gets transferred to the fluid.
In the first part of this chapter, we started with a problem to find the pump outlet temperature when
given the pump efficiency. Pump basically used to move liquids. In this chapter the user operated a pump
operation in HYSYS to model the pumping process. The user also been trained on how to connect streams
to unit operations such as pump.
At the end of this chapter, the user was trained to determine the pump outlet temperature when the pump
efficiency was given. On the other hand, when the outlet temperature was given, the pump efficiency can
be determined using HYSYS.
33
Chapter 4
Expander
This chapter begins with a problem to find the expander outlet temperature when given the expander
efficiency. The user will operate an expander operation in HYSYS to model the expansion process. At the
end of this chapter, the user will determine the expander outlet temperature when given expansion
efficiency or vice versa.
The Expander operation is used to decrease the pressure of a high pressure inlet gas stream to produce an
outlet stream with low pressure and high velocity. An expansion process involves converting the internal
energy of the gas to kinetic energy and finally to shaft work. The Expander calculates either a stream
property or an expansion efficiency.
There are several methods for the Expander to solve, depending on what information has been specified.
In general, the solution is a function of flow, pressure change, applied energy, and efficiency. The
Expander provides a great deal of flexibility with respect to what you can specify and what it then
calculates. You must ensure that you do not enable too many of the solution options or inconsistencies
may result.
The Expander operation is used to decrease the pressure of a high pressure inlet gas stream to
produce an outlet stream with low pressure and high velocity. A mixture of natural gas (methane, ethane
and propane) at 25oC and 20 bar is fed into an expander that has only 30% efficiency. The flowrate of the
natural gas is 100 kgmole/h and its outlet pressure from the compressor is 5 bar. Using Peng-Robinson
equation of state as a fluid package, determine the outlet temperature of the natural gas.
2. Click the Simulation button when you are ready to start building the simulation.
34
4.4 Adding an Expander
Name Expander
Outlet Out
Energy Work
Figure 4-1
35
3. Switch to the Parameters page. Change the Adiabatic Efficiency to 30%.
4. Go to the Worksheet tab. On the Conditions page, complete the page as shown in the following figure.
Specify the pressure for Out = 5 bar.
Figure 4-2
4.5 Discussion
This example shows that expansion gases can decrease their temperature. In this case, the expander
was only 30% efficient and it caused 31°C decrease in the temperature of the natural gas. The less
efficient a compressor is, the less the decrease in the temperature of the gases being expanded.
In the first part of this chapter, we started with a problem to find the expander outlet temperature when
given the expansion efficiency. The Expander operation is used to decrease the pressure of a high pressure
inlet gas stream to produce an outlet stream with low pressure and high velocity. In this chapter the user
operated an expander operation in HYSYS to model the expansion process.
36
Note: Bubble point calculations are performing by insert the Vapor/phase fraction = 0 in conditions of the
Material Stream, where if the temperature is given the bubble pressure will calculated (Figure 4-3), and if the
pressure is given the bubble temperature will calculated.
Dew point calculations are performing by insert the Vapor/phase = 1 in conditions of the Material Stream where if
the temperature is given the dew pressure will calculated, and if the pressure is given the dew temperature will
calculated.
Flash Calculations: If the temperature and pressure are given and the Vapor/Phase result between (0 – 1) that’s
mean the pressure is the flash pressure.
Figure 4-3
Ex: For the acetone(1)/ acetonitrile(2)/ nitromethane (3) system, we have the following Antoine
equations:
2940.46
ln P1sat = 14.5463− …(a)
T + 237.22
2945.47
ln P2sat = 14.2724 − …(b)
T + 224
2972.64
ln P3sat = 14.2043− …(c)
T + 209
Where T in oC and the vapor pressure are in kPa. Assuming that Raoult's law is approximate to this
system, calculate:
(a) P and yk, given that T = 80oC, x1 = 0.25, x2 = 0.35, and x3 = 0.4.
(b) P and xk, given that T = 70oC, y1 = 0.5, y2 = 0.3, and y3 = 0.2.
(c) T and yk, given that P = 80 kPa, x1 = 0.3, x2 = 0.45 and x3 = 0.25.
T and xk, given that P = 90 kPa, y1 = 0.6, y2 = 0.2, and y3 = 0.2.
37
Solution: Open new case
Vapour/ Phase Fraction = 0, T = 80oC and mole fractions: x1 = 0.25, x2 =0.35 , x3 = 0.4
Figure 4-4
The estimated pressure = 107.7 kPa, and the estimated composition will be as in Figure 4-5
Figure 4-5
38
Chapter 5
Compressor
This chapter begins with a problem to find the compressor outlet temperature when given the
compressor efficiency. The user will operate a compressor operation in HYSYS to model the compressing
process. At the end of this chapter, the user will determine the compressor outlet temperature when given
compressor efficiency or vice versa.
The Compressor operation is used to increase the pressure of an inlet gas stream. Depending on the
information specified, the Compressor calculates either a stream property (pressure or temperature) or
compression efficiency.
Compressors are used to move gases. The compressor increases the pressure of the gases. A mixture of
natural gas (C1, C2, C3, i-C4, n-C4, i-C5, n-C5, n-C6, C7+) at 100oC and 1 bar is fed into a compressor that
has only 30% efficiency. The flowrate of the natural gas is 100 kgmole/h and its outlet pressure from the
compressor is 5 bar. Using Peng-Robinson equation of state as a fluid package, determine the outlet
temperature of the natural gas.
Components: C1, C2, C3, i-C4, n-C4, i-C5, n-C5, n-C6, C7+
2. Component C7+ is not available in the Component Library. Therefore, you need to create this
component using Hypothetical.
Figure 5-1
39
2. From Method Menu choose Create and Edit Hypos, Figure 5-2.
Figure 5-2
A hypothetical component can be used to model non-library components, defined mixtures, undefined
mixtures, or solids. You will be using a hypothetical component to model the component in the gas
mixture heavier than hexane.
3. In the view for the hypo component you are creating, in the Component Name cell type C7+ and in the
Normal Boiling Point cell enter 110oC, Figure 5-3.
Figure 5-3
Since you do not know the structure of the hypothetical component and you are modeling a mixture, the
Structure Builder will not be used.
4. Click Estimate Unknown to estimate all the other properties and fully define the hypothetical
component.
40
8. You have now finished defining the fluid package. Click the Simulation button when you are ready to
start building the simulation.
1. From Object Palette; Double-click on the unit operation button of the Compressor.
2. On the Connections tab, enter the following information as shown in Figure 5-4:
41
Figure 5-4
3. Switch to the Parameters page. Change the Adiabatic Efficiency to 30%, Figure 5-5.
Figure 5-5
4. Go to the Worksheet tab. On the Conditions page, complete the page as shown in the figure 5-6. Set
the pressure for Comp_Out to be 500 kPa.
42
Figure 5-6
5.7 Discussion
This example shows that compressing gases can increase their temperature. In this case, the compressor
was only 30% efficient and it caused 165.3°C increase in the temperature of the natural gas. The less
efficient a compressor is, the greater the increase in the temperature of the gases being compressed.
In the first part of this chapter, we started with a problem to find the compressor outlet temperature when
given the compressor efficiency. Compressor basically used to move gases. In this chapter the user
operated a compressor operation in HYSYS to model the compressing process.
At the end of this chapter, the user was trained to determine the compressor outlet temperature when the
compressor efficiency was given.
43
Chapter 6
Heating & Cooling
This chapter begins with a problem to find the temperature of the cold stream passing through the heat
exchanger at the given stream conditions. In this chapter, HYSYS’s shell and tube heat exchanger will be
used to model the process. The heat exchanger performs two-sided energy and material balance
calculations. The heat exchanger is very flexible, and can solve for temperatures, pressures, heat flows
(including heat loss and heat leak), material streams flows, or UA.
In HYSYS, you can choose the Heat Exchanger Model for your analysis. Your choices include an End
Point analysis design model, an ideal (Ft=1) counter-current Weighted designed model, a steady state
rating method, Rigorous shell & tube model , and a dynamic rating method for use in dynamic
simulations. The dynamic rating method is available as either a Basic or Detailed model, and can also be
used in Steady State mode for Heat Exchanger rating.
1. Fixed tubes: This configuration has fixed the setting, straight pipes on both sides using welded stationary tube sheets. It is
low-cost construction. Its major disadvantages are maintenance since being fixed to the shell cannot be mechanically clean this
part of the tube; so it is very difficult to change a damaged tube.
2. U-tubes: In this configuration requires only a stationary tube sheet; however the cost is equal because it is expensive to
produce tubes that form. One advantage is that remove the tube bundle thus facilitating cleaning of the shell side. However,
cleaning the side of the tubes is very complicated because the U-shape is necessary to use cleaning equipment with moving
head. Such tubes should not be used for services with very dirty fluids.
3. Floating head: It is the most versatile and also the most expensive configuration. In this design, one of the head is relatively
fixed to the shell, while the other end is left free to “float” inside the shell. Thus thermal expansion joints without allowing and
likewise mechanically possible to clean both sides of the tubes. This design can be used for applications where both fluids (the
tubes and the shell) are dirty. To apply this design TEMA headers S and T are recommended.
Cold water at 20°C and 150 kPa is used to cool a hot stream of (ethanol & 1-Pentanol) in a shell and tube
heat exchanger. The inlet temperature and pressure of the hot stream is 65°C and 150 kPa, respectively. It
is required cool the hot stream to 35°C. If the flow rate of the hot and cold streams are 300 and 500 kg/h,
determine the temperature of the outlet cold stream passing through the exchanger.
1. Use HYSYS’ shell and tube heat exchanger to model the process.
2. Define the inlet and outlet conditions for the streams as given in the problem statement.
44
6.1.3 Defining the Simulation Basis
Property package: NRTL, click on the Binary Coeff tab, then click on Unknown only button
3. Click the Simulation button when you are ready to start building the simulation.
Temperature, oC 20 65
The heat exchanger performs two-sided energy and material balance calculations. The heat exchanger is
capable of solving the temperatures, pressures, heat flows (including heat loss and heat leak), material
stream flows, and UA.
45
Figure 6-1
3. Switch to the Parameters page. The pressure drops for the Shell and Tube sides, will be 30 kPa as
shown in the Figure 6-2.
Figure 6-2
46
4. Go to the Worksheet tab. On the Conditions page, Specify the temperature for Shell out = 35oC as in
figure 6-3.
Figure 6-3
Note: To achieve the design of heat exchanger successfully the pressure drop across shell and tube must
be > 0.
On the Design tab, Parameters page click on Size Rigorous Shell & Tube button, wait some minutes, from
new message appears, click Yes then click Rigorous Shell & Tube tab (Figure 6-4), then click on Model
Details button,
2
1
Click for Exchanger design
Figure 6-4
47
The design of Heat-Exchanger can see from figure 6-5.
Figure 6-5
For more results, click on Exchanger Geometry Folder in the Navigation Pane and can see the results in
each icon, figure 6-6.
Figure 6-6
48
For more accurate results: In the Process tab (Figure 6-5), set Fouling resistance = 0.001 hF/Btu
(0.00018 m2K/W) for hot and cold side.
To view full details of the design, click on the Model Details button at the bottom of the form.
Note: The standard diameter of the pipes in the Heat Exchangers are:
2. Simple Weighted Model: The use of this model for systems where heat curves exhibit a nonlinear, such as phase-change
systems for pure components on either side of the equipment recommended behavior. It is available only for counter-current
exchangers.
3. Simple Steady State Rating Model: It is similar to End point model but requires very detailed information on the
equipment geometry; if you have this information, we recommend running the calculation with this model.
49
Tubular Exchanger Manufacturers Association (TEMA)
50
Floating head
Fixed tubesheet
51
Removable U-bundle
Kettle reboiler
52
1 Shell pass , 2 tube passes
Shells in Series
Shells in parallel
53
4 Shell passes , 8 tube passes
54
Shell Baffle arrangement
55
6.2 Air Cooler
An Air Cooler will be used to cool a liquid process stream by air, Figure 6-7.
Figure 6-7
1. Create the following Fluid Package and add the Material stream
2. Specify the AIR COOLER connections and parameters as shown in the table below. The Parameters
page contains the fan rating and speed specifications.
56
3. From Worksheet tab enter the temperature of the product = 80oF
4. Using the following specifications, HYSYS is able to calculate the actual volumetric flow of air:
Number of Fans 1
5. At this point, the Air Cooler unit operation is fully specified. The exit stream conditions are displayed
in the Conditions page of the Worksheet tab.
6. The Air Cooler unit operation parameters are shown in the Results page of the Performance tab:
57
Air Cooler design
Design
Design Tab > Parameters page: click on the [Size Rigorous Air Cooler] button; wait some minutes to view
the results Figure 6-8.
Figure 6-8
58
For more accurate results: In the Process tab; specify the fouling resistance about 0.00018 m2 K/W in the
Hot Side, and 0.00035 m2K/W in cold side, then click on [Size Rigorous Air Cooler] another time
59
60
6.3 Fired Heater (Furnace)
It is performs energy and material balances in steady state or dynamic modes to model a direct
Fired Heater type furnace. This type of equipment requires a large amount of heat input. Heat is generated
by fuel combustion and transferred to process streams. A simplified schematic of a direct Fired Heater is
illustrated in the figure 6-9.
Figure 6-9
Note: The dynamic mode does not support multiple fuel streams, but uses one combined fuel/air stream. If
there are multiple fuel streams defined for steady state, when switching from steady state to dynamics a
new stream is created from the sum of the fuel streams and the air stream defined in steady state. When
switching back to steady state, one fuel and one air stream is created from the dynamic settings.
Ex: The feed stream (mole flows kmol/h: C1 = 300, C2 = 250, C3 = 200, i-C4 = 150, n-C4 = 150, i-C5 =
100, n-C5 = 100) at temperature of 50oC and pressure of 5 atm is heated to 300oC by fired heater and using
a methane gas as a fuel (40oC, 200 kPa) with 100% excess air at 40oC & 200kPa. What is the temperature
of flue gas and what is the amount of air required if the efficiency of furnace is 60%?
61
Solution:
Construct Air stream with molar composition (0.21 O2 & 0.79 N2) and with conditions as in problem
above
Figure 6-10
Enter to the parameters page and specify the Excess Air = 100% and Efficiency = 60%, Figure 6-11.
62
Figure 6-11
Figure 6-12
To study the effect of the molar flow of fuel required to heating the feed to the 300 oC, first delete the
temp of the product from worksheet and enter the molar flow rate of the fuel (e.g. 70 kmol/h). From
Navigation Pane OR from Home Ribbon click on Case study, click Add, and then click Add:
In the Object Column select Fuel with its variable = molar flow, and click Add
Similarly in the Object Column select Product with its variable = Temperature, and click Add
63
Specify the range of the Fuel molar flow: low bound = 50 kmol/h, high bound = 100 kmol & the step size
= 2 kmol/h as shown in the figure 6-13.
Figure 6-13
Click on Run and then click plot to see the results, Figure 6-14.
Figure 6-14
It is seen from the figure 6-14 that the temperature of the product at 300oC require molar flow of the
fuel = 63.4 kmol/h
64
Chapter 7
Flash Separator
This chapter begins with a problem to find the flowrate of the liquid and vapor outlet streams of the
flash separator. In steady state mode, the Separator divides the vessel contents into its constituent vapor
and liquid phases. The vapor and liquid in the vessel are allowed to reach equilibrium, before they are
separated.
A Flash Separator is performed to determine the product conditions and phases. The pressure at which
the flash is performed is the lowest feed pressure minus the pressure drop across the vessel. The enthalpy
is the combined feed enthalpy plus or minus the duty.
The Separator has the ability to back-calculate results. In addition to the standard application (completely
defined feed stream (s) being separated at the vessel pressure and enthalpy), the Separator can also use a
known product composition to determine the composition(s) of the other product stream(s), and by a
balance the feed composition.
7.1 Problem Statement (Gases mixture are separated by compression and cooling)
We have a stream containing 15% ethane, 20% propane, 60% i-butane and 5% n-butane at 50o F and
atmospheric pressure, and a flowrate of 100 lbmole/hr. This stream is to be compressed to 50 psia, and
then cooled to 32o F. The resulting vapour and liquid are to be separated as the two product streams. What
are the flowrates and compositions of these two streams?
Click the Simulation button when you are ready to start building the simulation.
Name Gas
Temperature 50oF
Pressure 1 atm
65
7.4 Adding a Compressor
Figure 7-1
3. Go to the Worksheet tab. At the Conditions page, specify the pressure for the Comp Gas = 50 psia.
The amount of duty (power) required can be seen by double click on the stream of the energy (Comp Q),
OR just place the mouse on it OR from Worksheet the amount of Comp Q (Duty) can be read it, OR from
the Parameters tab it can read the amount of Duty as in the figure below
66
7.5 Adding a Cooler
Figure 7-2
3. Switch to the Parameters page and set the pressure drop to 0 psia.
4. Switch to the Worsheet page and set the temperature of the Cool Gas to 32oF
67
Figure 7-3
3. Go to the Worksheet tab to preview the result as shown in the Figure 7-4 and Figure 7-5.
68
Figure 7-5, Stream’s compositions of the flash separator.
Figure 7-6 shows the process flow diagram of the compressor, cooler, and flash separator.
Figure 7-6
If the Cool Gas temperature is 10oF, what are the new flowrates and compositions of these two streams?
Create a case study to see the effect of the cool gas temperature on the vapor liquefied fraction molar or
the flow of the liquid from the separator.
69
Vessel Design
From the Navigation Pane, click on the Equipment Design Icon
70
From the Select Separator Windows, select V-100, click OK
From the Design tab, Sizing page, it can be select any specification from the Available Specifications group
and click on the Add Spec button to the list of the Active Specifications.
71
Miscellaneous Examples:
Ex.1: Steam Power Plant: Steam generated (basis = 1
Ex
QH
kg/h) in a power plant at a pressure of 8600 kPa and a
temperature of 500oC is fed to a turbine. Exhaust from 1
2
the turbine enters a condenser at 10 kPa where it is ∆P = 0
condensed to saturated liquid which then pumped to the
boiler. (a) Determine the thermal efficiency of a
practical cycle operating at these conditions if the
turbine efficiency and pump efficiency are both 75%. ∆P = 0
(b) If the rating of the power cycle of part (a) is 80000 4 3
kW, what is the steam rate and what are the heat transfer
QC
rates in the boiler and condenser? (note: Fluid Package
= NBS steam or Peng-Robinson)
Ex.2: Refregiration: A vapor-compression refrigeration system operates on the cycle of figure below.
The refrigerant is water. Given that the evaporation T = 277.15 K, the condensation T = 307.15 K,
η(compressor) = 0.76, and the refrigeration rate = 1200 kW, determine the circulation rate of the
refrigerant, the heat-transfer rate in the condenser, the power requirement, the coefficient of performance
of the cycle.
QC
∆P = 0
4 3
∆P = 0
1
2
QH Q
72
Ex.4: Claude Liquefaction Process:
Natural gas, assumed here to be pure methane, is liquefied in a Claude Process. Compression is to
60 bar. The expander and throttling exhaust to a pressure of 1 bar. Recycle methane at this pressue leaves
the exchanger system (point 15) at 295 K. Assume no heat leakes into the system from the surrounding, an
expander efficiency of 75%, and an expander exaust of saturated vapor. For a draw-off to the expander of
25% of the methane entering the exchanger system (x = 0.25), what fraction of z of the methane is
liquefied, and what is the temperature of the high-presure stream leaving the cooler? The temperture of the
stream 7 is – 75oC, the pressure of the stream 1 is 2 bar, the temperature and flow rate of stream 2 are
20oC and 100 kg/h respectively.
∆P=0
∆P=0 ∆P=0
73
Chapter 8
Logical Operations
The Set is an operation used to set the value of a specific Process Variable (PV) in relation to another
PV. The relationship is between the same PV in two like objects; for instance, the temperature of two
streams, or the UA of two exchangers.
The Set unit operation can be used in both Dynamic and Steady State mode.
The dependent, or target, variable is defined in terms of the independent, or source, variable according
to the following linear relation:
Y = MX + B
where:
M = multiplier
B = offset (intercept)
74
• Target Object: The stream or operation to which the dependent variable belongs. This is chosen by
clicking the Select Var button. This brings up the Variable Navigator property view.
• Target Variable. The type of variable you want to set, for example, temperature, pressure, and flow.
The available choices for Variable are dependent on the Object type (stream, heat exchanger, and so forth)
Your choice of Variable is automatically assigned to both the Target and Source object.
• Source Object. The stream or operation to which the independent variable belongs.
The Adjust operation varies the value of one stream variable (the independent variable) to meet a
required value or specification (the dependent variable) in another stream or operation.
Adjusted Variable
Target Variable
75
Target Value Group
Once the target object and variable are defined, there are three choices for how the target is to be satisfied:
• If the target variable is to meet a certain numerical value, select the User Supplied radio button (as
shown in the figure above), and enter the appropriate value in the Specified Target Value field.
• If the target variable is to meet the value (or the value plus an offset) of the same variable in another
stream or operation, select the Another Object radio button (as shown in the figure below), and select the
stream or operation of interest from the Matching Value Object drop-down list. If applicable, enter an
offset in the available field.
• If the target variable is to meet the value (or the value plus an offset) of the same variable specified in the
spreadsheet, select the SpreadSheetCell Object radio button (as shown in the figure below), and select the
cell that you want from the Matching Value Object drop-downlist. This allows the SpreadSheetCell to be
attached as an adjusted variable, and source to the target variable. You can also specify the offset in the
available field.
76
Parameters Tab
Once you have chosen the dependent and independent variables, the convergence criteria must be defined.
8.3 Balance
The Balance operation provides a general-purpose heat and material balance facility. The only
information required by the Balance is the names of the streams entering and leaving the operation. For
the General Balance, component ratios can also be specified. Since HYSYS permits streams to enter or
leave more than one operation, the Balance can be used in parallel with other units for overall material and
energy balances.
The Balance overrides the filtering of streams that HYSYS typically performs.
77
The Parameters tab contains two groups:
• Balance Type
• Ratio List
The Balance Type group contains a series of radio buttons, which allow you to choose the type of Balance
you want to use. The radio buttons are shown above.
This operation performs an overall mole balance on selected streams; no energy balance is made. It can
be used to provide material balance envelopes in the flowsheet or to transfer the flow and composition of a
process stream into a second stream.
• The direction of flow of the unknown is of no consequence. HYSYS calculates the molar flow of a feed
to the operation based on the known products, or vice versa.
Mass Flow
This operation performs an overall balance where only the mass flow is conserved. An application is the
modeling of reactors with no known stoichiometry, but for which analyses of all feeds and products are
available. If you specify the composition of all streams, and the flowrate for all but one of the attached
streams, the Mass Balance operation determines the flowrate of the unknown stream. This is a common
application in alkylation units, hydrotreaters, and other non-stoichiometric reactors.
• The flowrate must be specified for all but one of the streams. HYSYS determines the flow of that stream
by a mass balance.
• Energy, moles, and chemical species are not conserved. The Mass Balance operation determines the
equivalent masses of the components you have defined for the inlet and outlet streams of the operation.
78
Heat Flow
This operation performs an overall heat balance on selected streams. It can be used to provide heat
balance envelopes in the flowsheet or to transfer the enthalpy of a process stream into a second energy
stream.
• The composition and material flowrate must be specified for all material streams. The heat flow is not
passed to streams which do not have the composition and material flowrate specified, even if there is only
one unknown heat flow.
• The direction of flow for the unknown stream is of no consequence. HYSYS calculates the heat flow of a
feed to the operation based on the known products, or vice versa.
The most common application for this balance is to perform overall material (molar basis), and energy
balance calculations of selected process streams to either check for balances or force HYSYS to calculate
an unknown variable, such as a flowrate.
• The Mole and Heat Balance independently balance energy and material.
• The Mole and Heat Balance calculate ONE unknown based on a total energy balance, and ONE
unknown based on a total material balance.
• The operation is not directionally dependent for its calculations. Information can be determined about
either a feed or product stream.
• The balance remains a part of your flowsheet and as such defines a constraint; whenever any change is
made, the streams attached to the balance always balances with regard to material and energy. As such,
this constraint reduces by one the number of variables available for specification.
• Since the Mole and Heat Balance work on a molar basis, it should not be used in conjunction with a
reactor where chemical species are changing.
Similar to the Mass balance mode, this balance mode performs a balance on the overall mass flow. In
addition, however, energy is also conserved.
• Flow rate must be specified for all but one of the streams. HYSYS determines the flow of that stream by
a mass balance.
• Enthalpy must be specified for all but one of the streams. HYSYS determines the enthalpy of that stream
by a heat balance.
79
• Moles and chemical species are not conserved.
General Balance
The General Balance is capable of solving a greater scope of problems. It solves a set of n unknowns in
the n equations developed from the streams attached to the operation. This operation, because of the
method of solution, is extremely powerful in the types of problems that it can solve. Not only can it solve
unknown flows and compositions in the attached streams (either inlet or outlet can have unknowns), but
ratios can be established between components in streams. When the operation determines the solution, the
prescribed ratio between components are maintained.
• The General balance solves material and energy balances independently. An Energy Stream is an
acceptable inlet or outlet stream.
• The operation solves unknown flows or compositions, and can have ratios specified between
components in one of the streams.
Ratios
A Ratio, which is unique to the general Balance, is defined between two components in one of the
attached streams. Multiple ratios within a stream (for example 1:2 and 1:1.5) can be set with a single Ratio
on a mole, mass or liquid volume basis. Each individual ratio (1:2, 1:1, and 1:1.5), however uses a degree
of freedom.
To set a ratio:
1. On the Parameters tab of the Balance operation property view, select the General Balance radio button.
The Ratio List group should now be visible.
2. Click the Add Ratio button to access the Ratio property view.
80
• Name. The name of the Ratio.
• Ratio Type. Allows you to specify the Ratio Type: Mole, Mass, or Volume.
• Component/Ratio. Provides the relative compositions of two or more components. Other components in
the stream are calculated accordingly, and it is not necessary nor advantageous to include these in the
table. All ratios must be positive; non-integer values are acceptable.
81
Chapter 9
Piping Equipment
The Pipe Segment is used to simulate a wide variety of piping situations ranging from single or
multiphase plant piping with rigorous heat transfer estimation, to a large capacity looped pipeline
problems.
Ex: In the Figure 9-1, 700 kg/h of Water at 60oF and one atmosphere is being transferred with a 2 hp
pump that is 45% efficient. All the piping is 4 inch, schedule 40, & cast iron pipe except for the last
section, which is 2 inch schedule 40 steel pipe. The elbows are standard 90o and one reducer to connect
the 2 inch pipe.
Figure 9-1
82
Solution: (a)
1. Open New Case. Component: Water, Fluid package: Peng-Robenson. Enter Simulation
Environment
Figure 9-2
To enable the valve estimate the pressure drop, it must be given the size of the valve first.
4. Click on the Rating tab, the size of valve is determined by enter specific value of Cv.
So for current case the volumetric flow rate = 0.7014 m3/h (taken from worksheet)
83
Figure 9-3
Note the pressure drop is calculated to be 31.72 kPa. This is for 50% opening valve.
5. Add Pipe Segment by double click on Pipe Segment icon in the palette. Fill the data
required as in Figure 9-4 below:
84 Figure 9-4
6. Click on the Rating tab, Sizing page, Click on Append Segment button:
Enter the Length = 6 ft, Elevation change = - 6 ft, Figure 9-5.
Figure 9-5
Figure 9-6
85
9. From Nominal Diameters select 4.0 inch. In the Pipe Material Cell select Cast Iron from drop-
down list. Close the Pipe Info View, figure 9.6.
10. Click Append Segment button: from drop-down list of Pipe cell Select Elbow 90o Std, Inner
Diameter: 4.026 in. Material: Cast Iron
Figure 9.7
14. Select Heat Transfer Page: In the Heat Loss Cell: Enter 0 Btu/h.
Close the Pipe segment view.
86
Note: To estimate the heat transfer coefficient; click on the Estimate HTC button radio and enter
the Ambient temperature, then determine the parameters that required to estimate, figure 9.8
Figure 9.8
18. From Available Nominal Diameters Group select 4.00 and Click Specify button. In the Pipe
Material Cell select Cast Iron from drop-down list. Close the Pipe Info View.
19. Click Append Segment button: from drop-down list of Pipe cell Select Elbow 90o Std,
Inner Diameter: 4.026 in. Material: Cast Iron
20. Click Append Segment button: Length = 300 ft, Elevation change = 300 ft.
21. Click on View Segment button: In the Pipe Schedule Cell select Schedule 40 From drop-down
list.
87
22. In the Available Nominal Diameters Group select 4.00 and Click Specify button. In the Pipe
Material Cell select Cast Iron from drop-down list. Close the Pipe Info View.
23. Click Append Segment button: from drop-down list of Pipe cell Select Elbow 90o Std,
Inner Diameter: 4.026 in. Material: Cast Iron
24. Click Append Segment button: from drop-down list of Pipe cell Select: Coupling/Union
Inner Diameter: 2 in. Material: Cast Iron.
25. Click Append Segment button: Length = 375 ft, Elevation change = 0.
26. Click on View Segment button: In the Pipe Schedule Cell select Schedule 40 From drop-down
list.
In the Available Nominal Diameters Group select 2.00 and Click Specify button. In the Pipe Material
Cell select Cast Iron from drop-down list.
27. Click Append Segment button: from drop-down list of Pipe cell Select: Gate Valve: Half,
diameter = 2 in.
Close the Pipe Info View.
Note: the valve can be added to piping system separately from object palette as in the step 3 or from pipe
segment as in step 27.
Figure 9-9 shows the PFD of the piping system of schematic diagram of figure 9-1.
Figure 9-9
88
Figure 9-10
Note:
(1) The sizing of the pipe can be estimated from the Home Ribbon; Equipment Design → Pipe Sizing
Then click on the Select Stream button to choose the required stream, and determine the pressure drop to
estimate the diameter or the determine diameter required to estimate pressure drop.
89
(2) For more details of the pipes, in the Palette objects, Upstream group: Add the Aspen Hydraulic icon
to the flowsheet
Then click on
90
Standard pipe fittings
Standard Valves
91
PIPESYS
A pipeline must transport fluids over diverse topography and under varied conditions. Ideally this
would be done efficiently with a correctly sized pipeline that adequately accounts for pressure drop, heat
losses and includes the properly specified and sized inline facilities, such as compressors, heaters or
fittings.
Figure 9-11
Distance = Σ Runs
EX: The figure 9-12 shows the physical configuration of the performance of a small gas condensate
gathering system. The system consists of three wells distributed over an area of approximately 1.0 square
mile connected to a gas plant via a network of pipelines.
Figure 9-12
Field data shows that the wells are delivering the following rates:
92
Well MMSCFD
A 8.6
B 7.4
C 10.1
All three wells have the same composition. The residue of all the heavier components in the condensate
has a molecular weight of 122 and a density of 760 kg/m3. The characteristics of this component will be
accounted for by using the hypothetical component facility in HYSYS. The compositional analysis of the
gas condensate resulted in the following information:
Branch diameter
Branch 1 3”
Branch 2 3”
Branch 3 3”
Branch 4 4”
Branch 5 6”
A schedule 40 steel pipe is used throughout and all branches are buried at a depth of three feet. All pipes
are uninsulated. The following table summarizes the elevation data for each of the branches. The elevation
given for the pipe units is for the endpoint of the pipe (i.e., the downstream end).
The flow rate at each well is specified and is independent of the flow rate at each of the other wells. In
cases such as this, the system can be modeled with only one pressure drop determination per branch.
93
Simultaneous pressure and temperature calculations can be performed if the temperature at each wellhead
is also known.
The PFD generated by HYSYS for the completed case is shown in Figure 9-12.
Figure 9-13
A pressure of 1060 psia will be specified for Well A. PIPESYS will then determine the pressures
elsewhere in the network that are consistent with this specification.
The fluid temperatures at wells A, B and C are known and must be entered as fixed conditions.
PIPESYS will then perform an iterative pipeline calculation in branches where the upstream temperature
and downstream pressure are known. Temperatures of the blended fluids will be computed on a mass
basis downstream of the junctions of two or more streams.
Pure Components: C1, C2, C3, i-C4, n-C4, i-C5, n-C5, C6, N2, CO2, H2S
Property Package: Peng Robinson
Name C7+
Molecular weight 122
Ideal Liquid Density (lb/ft3) 47.45
94
6. Enter the compositional data for Well A as specified in the following table:
The composition of the streams Well A, Well B and Well C are identical
7. Use the Define from Other Stream... button on Well B’s stream view to define its composition
8. Open the Energy Streams tab on the Workbook view. Enter the names of five energy streams:
1. Add a PIPESYS Extension to the flowsheet (click F12, select Extension & double click on PIPESYS)
2. Complete the Connections tab as shown in the figure 9-14.
Figure 9-14
3. Go to the Elevation Profile tab.
4. In the Pipeline Origin Enter 0 ft into the Distance field and 2095 ft into the Elevation field figure 9-15
95
Note: A starting point for the profile must be specified at the top of the tab in the Pipeline Origin group,
using the Distance and Elevation input cells.
Figure 9-15
5. Add the first pipe unit for this extension on the Elevation Profile Page. The Pipe Unit view appears in
figure 9-16.
Figure 9-16
6. Go to the Dimensions tab.
7. Specify the Nominal diameter = 3 inches, and select Pipe Schedule = 40 as shown figure 9-17.
Figure 9-17
96
8. Go to the Heat Transfer tab:
Click on the Centre Line Depth cell and click the Default button. Leave all other parameters at their
default values as shown in figure 9-18. Close the Pipe Unit view.
Figure 9-18
9. Complete the Elevation Profile tab by entering 2110 ft for the Elevation parameter and 945 ft for the
Run parameter as in Figure 9-19.
Figure 9-19
10. Add the remaining two pipe units. Because all the pipe units for the extension have identical properties
to Pipe #1, you can use the Copy and Paste buttons as a time saving measure for adding the new
units.
97
11. Complete the elevation profile as shown by adding the Elevation and Run parameters for all units,
Figure 9-20.
Figure 9-20
12. Go to the Temperature Profile tab and enter 40°F in the Ambient Temperature field (in the Pipeline
Origin group) , Figure 9-21.
Figure 9-21
98
1. Enter the required information for the second extension as defined in the following table:
Name Branch 2
Distance
(ft) 0
Pipeline Origin
Elevation 2015
Elevation Profile (ft)
Output = FS4
On the Parameters page of the Mixer, select the Equalize All radio button in the Automatic Pressure
Assignment group.
99
Add the Third PIPESYS Extension
This third PIPESYS extension is used to represent Branch 4. It consists of a single pipe unit.
1. Enter the data for the third PIPESYS extension as defined in the following table.
Distance 0
Pipeline Origin (ft)
Elevation 2090
Elevation Profile (ft)
The following table contains the information required to complete the fourth PIPESYS extension:
100
Tab Object Specification
Name Branch 3
Distance 0
Pipeline Origin (ft)
Elevation 2125
(ft)
101
Add the Fifth PIPESYS Extension
In this segment, the total gas flows from Wells A, B and C are merged and the endpoint of the extension
is the gas plant.
Enter the information for the final extension as defined in the following table:
Distance 0
Pipeline
(ft)
Origin
Elevation 2090
(ft)
Elevation Profile
Pipeline Unit Pipe # 1
On the Parameters page of the Mixer, select the Equalize All radio button.
102
Stream Object Specification
Save your case, and view the results from Material streams tab in the workbook.
Figure 9-22
The adjust will be used to maintain a constant pressure of 1000 psia at the gas plant with the pressure at
Well A being the adjusted variable.
103
Object Specification
Method Secant
2. Click the Start button at the bottom of the Adjust property view to begin the adjust calculations.
(HYSYS may need several minutes to reach a solution.)
Add an Inline Compressor: The next step will be to add an inline compressor at the upstream end
of PIPESYS Branch 5.
2. Modify the pipe extension, so that the Compressor is placed before the Pipe. See the figure below for
placement of the pipe units:
104
2. Open the Inline Compressor view:
4. Use the following information to complete the tabs as specified on the Inline Compressor view.
Object Specification
The solution process may take several minutes to perform the iterative calculation for the PIPESYS
network and converge.
To compare the performance between two compressors, the same calculations are repeated using a 750 hp
compressor.
6. In the Brake Power group, change the 1000 hp in the Specified cell to 750 hp.
105
When HYSYS completes the iteration, the new pressures are:
Well Pressure (psia)
The results indicate that the 1000 hp compressor would increase production by less than 5%, over that
of the 750 hp compressor. It is therefore reasonable to conclude that adding compression to the system is
worthwhile since both compressors lower the wellhead pressures by a large amount, but the small increase
in production may not be enough to justify the choice of the 1000 hp compressor. For this example,
assume that economic and engineering considerations favour installing the 750 hp compressor.
106
Chapter 10
Reactors
This chapter begins with a problem to develop a model that represents the partial oxidation reaction of
methane to produce hydrogen. The partial oxidation method relies on the reaction of the methane with air
in order to produce carbon oxides and hydrogen. The user will learn how to add the conversion reactions
and reactions sets in HYSYS.
This reaction type does not require any thermodynamic knowledge. You must input the stoichiometry
and the conversion of the basis reactant. The specified conversion cannot exceed 100%. The reaction will
proceed until either the specified conversion has been reached or a limiting reagent has been exhausted.
Conversion reactions may not be grouped with any other form of reaction in a reaction set. However,
they may be grouped with other conversion reactions and ranked to operate either sequentially or
simultaneously. Lowest ranking occurs first (may start with either 0 or 1). Just as with single reactions,
simultaneous reactions cannot total over 100% conversion of the same basis.
Conversion reactions cannot be used with Plug Flow Reactors or CSTRs. In general, they should only be
used in Conversion Reactors.
The interest in production of hydrogen from hydrocarbons has grown significantly in the last decade.
Efficient production of hydrogen is an enabling technology, directly related to the fuel cell energy
conversion device. The conversion of fuels to hydrogen can be carried out by the partial oxidation. The
partial oxidation method relies on the reaction of the fuel for example methane with air in order to produce
carbon oxides and hydrogen.
1. The first step in simulating a hydrogen production is choosing an appropriate fluid package. Enter the
following values in the specified fluid package view:
Reactions in HYSYS are added in a manner very similar to the method used to add components to the
simulation:
107
1. Click on the Reactions folder in the Navigation Pane, then add Reaction Set by click on Add button,
Figure 10-1.
Figure 10-1
Figure 10-2
108
3. Choose Conversion as the type from the displayed list & click Add Reaction
4. Under the Active Reaction: Double click on Rxn-1; Figure 10-3 shows the windows of Conversion
Reaction.
Figure 10-3
4. Add components according to the reaction 1. Give minus sign for Stoich Coeff of reactants and
Positive sign for Stoich Coeff of products as in Figure 10-4.
Figure 10-4
5. In Base Component: select Methane , Rxn Phase: VapourPhase & Give the conversion of the reaction
in Co cell = 40%.
6. Add the second reaction by the similar way for reaction 1: enter the information as in Figure 10-5:
109
Figure 10-5
Conversion reactions can be grouped with other conversion reactions and ranked to operate either
sequentially or simultaneously. Lowest ranking occurs first (may start with either 0 or 1).
1. To make the reactions operate sequentially, in the Set-1, click the Ranking… button and enter the
information as shown in Figure 10-6:
Figure 10-6
After the reaction set has been created, it must be added to the current fluid package in order or HYSYS
to use them.
1. Click on Add to FP button, from the new windows click on Add Set to Fluid Package.
Once the reaction set is added to the Fluid Package, you can enter the Simulation and begin construction
of the simulation.
111
10.1.8 Adding the Conversion Reactor
1. From the Object Palette, click Column tab (Figure 10-9) and select the Conversion Reactor (Figure
10-10) and enter it into the PFD.
2. Name this reactor Oxidation Reactor and attach Methane and Air as feeds. Name the vapor outlet
Ox_Vap and even though the liquid product from this reactor will be zero, we still must name the stream.
Name the liquid product stream as Ox_Liq, Figure 10-11.
Figure 10-11
Note: If the Reaction is adiabatic the Energy stream must be empty and if the Reaction at specified
temperature, it must be add Energy stream. In this example the Reaction is adiabatic
In the Parameters Page it must enter the pressure drop if there is, in this example the ∆P = 0.
112
3. On the Reactions tab, Details page and from the Reaction Set select Set-1. This will automatically
connect the proper reactions to this reactor, figure 10-12.
Figure 10-12
4. Go to the Worksheet tab. On the Composition page, analyze the composition in the Ox_Vap stream.
Save your Case as Conversion. The final view of the PFD will like figure 10-13.
Figure 10-13
113
10.1.9 Review and Summary
In the first part of this chapter, we started with a problem to develop a model that represents the partial
oxidation reaction of methane to produce hydrogen. The partial oxidation method relies on the reaction of
the methane with air in order to produce carbon oxides and hydrogen. The user also learns how to add the
conversion reactions and reactions sets in HYSYS.
This chapter begins with a problem to develop a model that represents the water gas shift reaction. The
role of the WGS reaction is to increase the H2 yield and decrease the CO concentration to cell
requirements to prevent the anode being poisoned and the cell efficiency abruptly drops. The user will
learn how to add the equilibrium reactions and reactions sets in HYSYS.
The Equilibrium reactor is a vessel which models equilibrium reactions. The outlet streams of the
reactor are in a state of chemical and physical equilibrium. The reaction set which you attach to the
Equilibrium reactor can contain an unlimited number of equilibrium reactions, which are simultaneously
or sequentially solved. Neither the components nor the mixing process need to be ideal, since HYSYS can
compute the chemical activity of each component in the mixture based on mixture and pure component
fugacities.
You can also examine the actual conversion, the base component, the equilibrium constant, and the
reaction extent for each reaction in the selected reaction set. The conversion, the equilibrium constant and
the extent are all calculated based on the equilibrium reaction information which you provided when the
reaction set was created.
The new application of hydrogen as a raw material for fuel cells for mobile power sources (PEM fuel
cells) requires that the anode inlet gas have a CO concentration lower than 10-20 ppm. Otherwise, the
anode is poisoned and the cell efficiency abruptly drops.
Hence, if the hydrogen is produced from hydrocarbon or alcohol reforming, purification is required in
order to reduce the CO levels to cell requirements. The most technologically feasible purification train
consists of a water gas shift reaction (WGS). The reaction
CO + H2O → CO2 + H2
has been employed for 40 years in the industrial process for H2 production from liquid and gaseous
hydrocarbons. The role of the WGS reaction is to increase the H2 yield and decrease the CO
concentration, which is a poison for some catalysts used.
1. For this section, you will be using the saved case from the previous section (section 10.1: Conversion
Reaction) with one additional component, H2O.
114
3. Click the and from Navigation Pane click on Components List-1, Figure 10-14.
Figure 10-14
6. Add Reaction-Set by Click the Add button, then click Add Reaction, choose Equilibrium as the type
from the displayed list, then click Add Reaction, Figure 10-15.
Figure 10-15
8. From the Components column add the components of the following reaction with its stoichiometry
115
Figure 10-16
Click Keq tab and see the source of the Keq, Figure 10-17.
Figure 10-17
In this example the Keq depends on the Gibbs Free Energy which is taken from the data base of
components that saved in program.
After the reaction set has been created, it must be added to the current fluid package in order for HYSYS
to use them.
1. Click on Add to FP button, then click on Add set to the fluid package button.
Once the reaction set is added to the Fluid Package, Click on the Simulation button in the Navigation
Pane, from the new windows click No.
116
10.2.4 Adding a Feed Stream
1. From the Object Palette, click on Column Group, Add the Equilibrium Reactor to the PFD.
2. Name this reactor WGS Reactor and attach Ox_Vap and Steam as feeds. Name the vapor outlet
WGS_Vap and even though the liquid product from this reactor will be zero, we still must name the
stream. Name the liquid product stream as WGS_Liq as shown in Figure 10-18.
Figure 10-18
3. On the Reactions tab, Details page, from Reaction Set select Set-2, Figure 10-19.
117
Figure 10-19
4. Go to the Worksheet tab. On the Composition page, analyze the composition in the WGS_Vap
stream.
Calculate the percentages of the following (compare results from Chapter 9):
CO reduced: _________________________________
In Aspen HYSYS you have the ability to print Datasheets for streams, operations, and workbooks.
From the Variables Group, click Add, From the Variables column Select Master Comp Mole Frac,
118
From the All/Single group select All radio, then click OK, Figure 10-20.
Figure 10-20
The mole fraction of each component will add to the Workbook as shown in Figure 10-21.
Figure 10-21
6. Right-click (object inspect) the Workbook title bar. The Print Datasheet pop-up menu displays.
7. Select Print Datasheet. The Select Datablock view displays, Figure 10-22.
Figure 10-22
119
8. From the list, you can choose to print or preview any of the available datasheets.
Figure 10-23 represents the PFD of the Conversion & Equilibrium Reactors.
Figure 10-23
In the first part of this chapter, we started with a problem to develop a model that represents the water
gas shift reaction. The role of the WGS reaction is to increase the H2 yield and decrease the CO
concentration to cell requirements to prevent the anode being poisoned and the cell efficiency abruptly
drops. The user also learns how to add the equilibrium reactions and reactions sets in HYSYS.
120
Miscellaneous Examples
HYSYS is an interactive process engineering and simulation program. It is powerful program that you
can use to solve all kinds of process related problems. However, since you have to provide various
conditions and choices in order to solve a problem, you cannot use it effectively unless you have good
knowledge about the process and solution procedures.
The toluene production process is started by heating n-heptane from 65 to 800oF in a heater. It is fed to a
catalytic reactor, which operates isothermally and converts 15 mol% of the n-heptane to toluene. Its
effluent is cooled to 65oF and fed to a separator (flash). Assuming that all of the units operated at
atmospheric pressure, determine the species flow rates in every stream. (basis n-heptane = 100 lbmol/h)
Solution
Figure 10-23
5. Double click on reactor. Choose Set-1 as Reaction set and close the window.
6. Now, open worksheet, and type in all the known conditions for the streams. Note that only blue colored
fonts are the values that you specified. If you more information than the degree of freedom allows, it will
give you error messages.
121
Ex.1-b: Modification of Process for the Improvement
Inspection of the calculation results of Example 1-a shows that the cooling duty is comparable to the
heating duty, suggesting that the utility load can be reduced by preheating the feed stream with hot
reactor product. Modify the process by adding a heat exchanger before the heater. This can be
accomplished in the PFD as shown in Figure 10-24:
Figure 10-24
∆p = 0 for both tube side and shell side Heat-exchanger. Choose Weighted Exchanger as Model.
You still need to specify one more condition. Open the Worksheet and specify the temperature of Pre-
Feed stream to 600oF. You may change this temperature to see how it affects the Heat-duty.
You can change the Pre-Feed stream temperature and see how it affects the H-Duty and UA (heat transfer
coefficient x interfacial area). Increasing Pre-Feed temperature can reduce the H-Duty, but it will increase
UA, which means that you need a heat exchanger with more interfacial area (bigger and with more inner
pipes). Obviously, there will be upper limit of Pre-Feed temperature no matter how good your heat
exchanger is. You can see this effect by changing the temperature and recording the change of other
values. This can be done by using Case study function. The process can be described as follows:
a. In the Home Ribbon Open Case Study function. Click Add button, and click Add button, select Pre-
Feed as object, Temperature as Variable and click Add button. Do the same for Heat-Duty as object,
Heat Flow as Variable and E-102 as object, UA as Variable. Close the window.
b. In the Independent variable : Low Bound = 500oF , High Bound = 620oF, and 10oF for Step Size.
c. Click Run. After a few seconds, click Plots, the results are shown in figure 10-25.
122
Figure 10-25
So, the desired temperature of Pre-Feed stream can be evaluated from the figure 10-25 which equal to
582oF.
123
Ex.2: Steam-Methane Reformer (Hydrogen Production)
Steam reformation of methane is often undertaken in conjunction with processes which require large
amounts of hydrogen. Successive reaction stages take advantage of thermodynamics and catalysts to
enhance the production of hydrogen at the expense of the by-product gases carbon monoxide and dioxide.
Finally, remaining carbon oxides are converted back into methane as completely as possible to minimize
CO and CO2 carryover into the downstream process. PFD is shown in Figure 10-26.
Figure 10-26
1. Adding the Reactions: The reactions which take place in this simulation are:
(Equilibrium)
(Equilibrium)
(Equilibrium)
(Conversion)
2. Click on Reaction folder; click Add button, click Add Reaction button, Select Equilibrium and
click Add Reaction button and enter the necessary information which is equivalent to the
following reaction:
4. Click on the Keq tab , click the Keq Vs T table radio and enter the following values:
124
Basis group: Basis = Partial Pressure , Phase = Vapourphase, Basis Units = atm
5. Add the second Equilibrium reaction by selecting the reaction type as Equilibrium:
7. For reaction 2, proceed as above (by enter the stoichiometry for each component) click on Keq tab &
select Keq Vs T table radio and enter the following values for the Equilibrium Constant:
Basis group: Basis = Partial Pressure , Phase = Vapourphase, Basis Units = atm
In the absence of a catalyst and at 430°C (800°F), the rate of reaction number 1 in the Shift Reactor is
negligible, and reaction number 2 becomes the only reaction.
Rename the Set-1 of the reactions to Reformer by Right click on Set-1 in the Navigation pane
125
HYSYS contains a library of some of the most commonly encountered chemical reactions with their
Equilibrium Constants. For the Shift Reactor, you will use the library values for the Equilibrium Constant.
8. Add second Set of reactions by click on Reaction Folder in the Navigation Pane, then click on Add
button and rename this set to Shift
9. Add the third Equilibrium reaction by selecting the reaction type as Equilibrium. On the Library tab,
highlight the reaction with the form CO + H2O = CO2 + H2. Press the Add Library Rxn button. This
adds the reaction and all of the reaction’s data to the simulation. Rename the reaction to Shift1
10. Add Shift Set of reaction to FP By click on Add to FP, then click on Add Set to Fluid Package
14. Move to the Basis tab: Base Component = CO , Co = 100, Phase = Vapourphase
Add Methanation Set of reaction to FP By click on Add to FP, then click on Add Set to Fluid Package
Add a Set Operation: A Set operation is needed to fix the steam rate relative to the methane feed rate.
126
Add a Set operation with the following information:
Add the Steam Reformer: An Equilibrium Reactor will be used to simulate the Steam Reformer.
From the Object Palette, Column Operations, Add the Equilibrium Reactor to the PFD. Make the
following connections:
This will automatically connect the proper reactions to this Reactor and the Reactor will solve.
How much CO and H2 were produced in the reaction; i.e. what is the molar flowrate of these two
compounds in the reactor’s product stream? __________ & __________
127
∆P 24 kPa
∆P 70 kPa
Add a Component Splitter to model the Amine Plant. The purpose of this Splitter is only to remove the
CO2 present in the flow. The connections are shown below:
128
On the Splits page, specify the Fraction to Overhead as 1.0 for Methane, CO, Hydrogen and H2O. The
"Fraction to Overhead" for the CO2 must be 0; this will force all CO2 to the bottom and all other
components to the top.
Add a Heater:
On the Reactions tab, choose the Methanation from the Reaction Set drop down menu.
What is the affects of varying the amount of Steam to the process on the flow of CO2 entering the Amine
Plant and the flow of methane leaving the Methanator Reactor?
Compare steam flow rates of 1.5, 2.0, 2.5, 3.0, 3.5 and 4.0 times that of the natural gas.
1. From the Home Ribbon click on Case Study , click Add button, then click Add button
2. Select Shift Prod as the object, Master Comp Molar Flow as the variable, and CO2 as the variable’s
specific, then click Add.
129
4. In the Independent Variable, Enter values for Low Bound = 1.5 , High Bound = 4 , and
Step Size = 0.25.
5. Click Run button to begin the case study analysis. Press the Plots button to view the calculations as in
figure 10-27.
Figure 10-27
From the Figure above, the optimum value of the Set-1Multiplier = 2.22
130
Ex.3: Synthesis Gas production
In this example the Synthesis Gas production will simulate (synthesis gas is important in Ammonia
synthesis where the molar ratio of H2 to N2 is 3:1). The main role of the synthesis gas plant is to convert
natural gas (primarily methane) into hydrogen. In most synthesis gas plants, four reactors are used.
However, in our simulation five reactors will be used to model this process. This is because the combustor
(single vessel) will be modeled as two reactors in series with two different reaction types. The first reactor
is a Conversion reactor and the second is an Equilibrium reactor. PFD is shown in figure 10-28.
Figure 10-28
Move to the Basis group and enter the information as shown in figure 10-29:
131
Figure 10-29
b) Repeat steps above for another Conversion Reaction. Use the following data:
Add the Reaction Set to the Fluid Package by click Add to FP button then click Add Set to Fluid Package.
a) Add two reactions above to the Combustor Reaction Set by click the arrow beside the Add Reaction
button and select Existing Reaction
Add 3rd reaction to the Combustor Set by click Add Reaction button with the following informations:
Add Combustor Set to the fluid package by click Add to FP button then click Add Set to Fluid Package.
132
3. Add new Reaction Set & rename it to Shift
Add an Equilibrium reaction: the reaction: CO + H2O = CO2 + H2 . Click Library button then click the
Add Library Rxn button. This adds the reaction and all of the reaction’s data to the simulation.
Add Shift Set to the fluid package by click Add to FP button then click Add Set to Fluid Package.
133
2. Add another Set Operation; Set-2
1. Add Conversion Reactor, Named: Reformer, Inlet: Natural Gas and Reformer Steam
Vapor Outlet: Combustor Feed, Liquid Product: Reformer LP, Energy: Reformer Q
3. Select Conversion% radio , Change the Co value for Rxn-1 to 40% and for Rxn-2 to 30%
Name Combustor
6. Add Shift Reactors: The shift Reactor will modeled as Equilibrium Reactor
134
7. Enter another Equilibrium Reactor with the following data:
Name Shifter 1
Name Shifter 2
In order to control the temp of the product leaving the Combustor (2nd Conversion reactor), the flow
rate of the steam to this reactor is controlled. It is desired to have an outlet temp from the 1 st shift reactor
of 930oC.
1. Add the Adjust Operation from the Object palette and enter the data as in figure 10-30:
135
Figure 10-30
10.000
1.0
<Unbounded>
1000
Figure 10-31
136
Figure 10-32
Object = Synthesis Gas, Variable = Master Comp Molar Flow, Variable Specify = Hydrogen
And another time Click Add Import: Object = Synthesis Gas, Variable = Master Comp Molar Flow,
5. Click Spread Sheet tab: Add a ratio formula to an empty cell in the spreadsheet tab,
A second Adjust Operation will be used to control the air flow rate. The air flow rate determines the ratio
of H2 to N2 in the synthesis gas product. We want this value to be at 3.05.
7. On the parameters page, choose a tolerance of 0.001 and step size of 20 kmol/h.
In this case, the two adjust operation might interfere with each other while they are solving. This is
because changing either adjusted variable affects both target variables. To prevent this interference, the
Adjusts can be set to solve simultaneously.
137
8. On the Parameters tab of the ADJ-1 operation, select the Simultaneous Solution checkbox.
10. Start the simultaneous Adjsuts solving by using the Start button on the Adjust or in the
Sim Adj Manager tab.
138
10.3 CSTR (continuous stirred tank reactor)
In this chapter, a flowsheet for the production of propylene glycol is presented. Propylene oxide is
combined with water to produce propylene glycol in a continuously-stirred-tank reactor (CSTR).
The propylene oxide and water feed streams are combined in a mixer. The combined stream is fed to a
reactor, operating at atmospheric pressure, in which propylene glycol is produced.
The next task in defining the Fluid Package is providing the binary interaction parameters.
1. Click the Binary Coeffs tab of the Fluid Package property view.
In the Activity Model Interaction Parameters group, the Aij interaction table appears by default. HYSYS
automatically inserts the coefficients for any component pairs for which library data is available. You can
change any of the values provided by HYSYS if you have data of you own.
In this case, the only unknown coefficients in the table are for the 12C3Oxide/12-C3diol pair. You can
enter these values if you have available data, however, here, you will use one of HYSYS’s built-in
estimation methods instead.
2. Next, you will use the UNIFAC VLE estimation method to estimate the unknown pair.
3. Click on Unknowns Only button. HYSYS provides values for the unknown pair.
The final Activity Model Interaction Parameters table for the Aij coefficients appears in Figure 10-33.
Figure 10-33
4. To view the Bij coefficient table, select the Bij radio button. For this case, all the Bij coefficients will
be left at the default value of zero.
139
10.3.6 Defining the Reaction
The reaction between water and propylene oxide to produce propylene glycol is as follows:
1. In the Reactions group, click on Add button to Add Reaction-Set. The Reaction-Set property view
appears.
2. Click on Add Reaction button. Select Kinetic type & click on Add Reaction button and close the
Reaction property view.
3. Under Active Reactions, Double click on Rxn-1, add the components with its Stoichiometry according
to the reaction above, Figure 10-34.
Figure 10-34
HYSYS provides default values for the Forward Order and Reverse Order based on the reaction
stoichiometry. The kinetic data for this case is based on an excess of water, so the kinetics are first order
in Propylene Oxide only.
3. In the Fwd Order cell for H2O, change the value to 0 to reflect the excess of water, the Rev Order of
the 12-C3diol also set to 0. The stoichiometry tab is now completely defined and appears as shown in
figure 10-35.
140
4. In the Basis group, accept the default value of Molar Concn. Click in the Base Component cell, By
default, HYSYS has chosen the first component listed on the Stoichiometry tab, in this case Propylene
oxide, as the base component.
5. In the Rxn Phase cell, select CombinedLiquid from the drop-down list. The completed Basis tab
appears below.
6. The Arrhenius parameters for the kinetic reaction must be given. In this case, there is no Reverse
Reaction occurring, so you only need to supply the Forward Reaction parameters.
E (activation energy) = 3.24e4 (Btu/lbmole). The status indicator at the bottom of the Kinetic Reaction
property view changes from Not Ready to Ready, indicating that the reaction is completely defined,
Figure 10-35.
Figure 10-35
8. The final task is to make the set available to the Fluid Package, which also makes it available in the
flowsheet. Add the Reaction Set to the Fluid Package by clicking on Add to FP button, then click on Add
set to Fluid Package button.
141
Add another new Material stream with the following values.
Now that the feed streams are known, your next task is to install the necessary unit operations for
producing the glycol.
The first operation is a Mixer, used to combine the two feed streams. Enter the necessary information as
shown in Figure 10-36:
Figure 10-36
Note: In the parameters page of the Mixer there are two options in the Automatic Pressure assignment:
• Equalize All: Selected in the case of one stream has known pressure and the other is unknown.
• Set Outlet to Lower Inlet: Selected in the case of two streams are known pressures.
142
Installing the Reactor
2. Name this reactor CSTR and attach Mix_Out as feed. Name the vapor outlet CSTR Vent and the
liquid product stream as CSTR Product as in Figure 10-37.
Figure 10-37
3. On the Reactions tab, Details page, from the Reaction set, select Set-1. This will automatically
connect the proper reactions to this reactor.
4. The next task is to specify the Vessel Parameters. From Dynamics tab, the reactor has a volume of
280 ft3 and Liq Volume Percent = 85%, Figure 10-38.
143
Figure 10-38
5. Click on the Worksheet tab. At this point, the reactor product streams and the energy stream coolant
are unknown because the reactor has one degree of freedom. At this point, either the outlet stream
temperature or the cooling duty can be specified.
6. Initially the reactor is assumed to be operating at isothermal conditions; therefore, the outlet
temperature is equivalent to the feed temperature, 75oF. In the CSTR Product column, enter 75oF at the
Temperature cell, Figure 10-39.
Figure 10-39
144
7. There is no phase change in the Reactor under isothermal conditions since the flow of the vapor product
stream CSTR Vent is zero. In addition, the required cooling duty has been calculated and is represented by
the Heat Flow of the Coolant stream. The next step is to examine the Reactor conversion as a function of
temperature.
8. Click the Reactions tab, then select the Results page. The conversion appears in the Reactor Results
Summary table, Figure 10-40.
Figure 10-40
9. Under the current conditions, the Actual Percent Conversion (Act.% Cnv.) in the Reactor is 70.88%.
You need to adjust the reactor temperature until the conversion is in the 85-95% range.
Figure 10-41 shows the final process flow diagram (PFD) of the CSTR reactor
Figure 10-41
145
10.4 Plug Flow Reactor
The Plug Flow Reactor (PFR) as shown in figure 10-42, also known as the Tubular Reactor, generally
consists of a bank of cylindrical pipes or tubes. The flow field is modeled as plug flow, implying that the
stream is radially isotropic (without mass or energy gradients). This also implies that axial mixing is negligible.
Figure 10-42
As the reactants flow the length of the reactor, they are continually consumed, hence, there is an axial
variation in concentration. Since reaction rate is a function of concentration, the reaction rate also varies
axially (except for zero-order reactions). To obtain the solution for the PFR (axial profiles of
compositions, temperature, and so forth), the reactor is divided into several subvolumes. Within each
subvolume, the reaction rate is considered to be spatially uniform.
Basis = Partial pressure , Base component = E-Bz, Rxn phase = VapourPhase, Basis Units = bar
In addition to the above reaction the following undesired kinetic reaction is also occurred:
E-Benzene → Bz + Ethylene
Basis = Partial pressure, Base component = E-Bz, Rxn Phase = VapourPhase, Basis Units = bar
The feed enter to the adiabatic PFR (607oC, 138 kPa, 10180 kmol/h & Mass frac E-Bz = 0.33, H2O =
0.67). What is the mass flow rate of the styrene in the product?
146
Solution:
Click on the Reaction tab in the Navigation pane, click Add, and then click Add Reaction
Select Simple rate and Enter the data required as in Figure 10-43:
Figure 10-43
Add another reaction of type kinetics, and enter the data as in figure 10-44:
147
Figure 10-44
Construct the Feed stream, Enter the data of the feed as mentioned in the problem,
PFR reactor
Figure 10-45
Double click on the PFR Reactor, Figure 10-46 will appear, Enter the Feed and Product streams.
148
Figure 10-46
Rating tab: add the data required as follows: Volume = 5 m3 , length = 2 , Void Fraction = 0.445
Figure 10-47
Go to the Reaction tab, Results page and see the conversion of the each reaction.
149
Ex.2: Ammonia Production Plant
Ammonia is produced in a process known as the Haber process, in which nitrogen and hydrogen react in
the presence of an iron catalyst to form ammonia. The hydrogen is formed by reacting natural gas and
steam at high temperatures and the nitrogen is supplied from the air. The PFD is shown in figure 10-48.
Figure 10-48
Basis: Partial Pressure, Base Component: N2 , Rxn Phase: Vapor phase , Basis Unit: atm
Enter Simulation
150
Composition in Mole Flow as in the table
Composition kmol/h
151
Parameter Page: ∆P = 15 kPa tube side
Notes: (1) The process flow diagram can be printed in the PDF file by right click on any space and
select .
(2) The report about the process flow diagram with all conditions and properties …etc can get by select all
unit operations in the (PFD) then right click on any unit operations and select then choose
preview to view the report before the printing.
152
Ex: The hydrogen production is achieved by catalytic reaction:
CH3OH + H2O → CO2 + 3H2
in the PFR and vapor phase. The dimension of reactor: L = 0.12 m, D = 0.04 m.
fluid property: UNQAUC , click on Binary Coeff tab and click on the Unknowns only button.
Figure 10-49
153
Click on the Reaction Rate tab and enter the data as in figure 10-50:
Figure 10-50
154
Enter simulation Environment.
∆P = 30 kPa,
Figure 10-51
155
Figure 10-52
Case Study: Investigate the effect of Reactor length on the hydrogen product. The results will appear as
in Figure 10-53.
Figure 10-53
156
Exercise: It is desired to evaluate the design of a reactor for the production of allyl chloride according
to the main reaction.
The reaction mixture containing 4 mol of propylene per mole of chlorine and get into the reactor at
o
392 F at a rate of 0.85 lbmol/h. The operating pressure can be assumed constant and equal to the feed
with a value of 29.4 psia. Initially it is expected to use a tubular reactor with 2-in. of inside diameter,
jacketed and working with boiling ethylene glycol as a coolant, so that the temperature inside is getting
constant at 392oF. The heat transfer coefficient inside is 5 Btu/h ft2oF.
The speed equations of the two reactions, expressed in lbmol/h ft3, are:
157
Chapter 11
Separation Columns
11.1 Absorber
This chapter introduces the use of HYSYS to model a continuous gas absorption process in a packed
column. The only unit operation contained in the Absorber is the Tray Section, and the only streams are
the overhead vapor and bottom liquid products. There are no available specifications for the Absorber,
which is the base case for all tower configurations. The conditions and composition of the column feed
stream, as well as the operating pressure, define the resulting converged solution. The converged solution
includes the conditions and composition of the vapor and liquid product streams.
CO2 is absorbed into propylene carbonate in a packed column. The inlet gas stream is 20 mol% CO2
and 80 mol% methane. The gas stream flows at a rate of 2 m3/s and the column operates at 60oC and 60.1
atm. The inlet solvent flow is 2000 kmol/h. Use HYSYS to determine the concentration of CO 2 (mole%)
in the exit gas stream, the column height (m) and the column diameter (m). The column equivalent to 10
theoretical trays.
61
158
Add another new Material stream with the following values.
61
7200 m3/h_gas
1. From the Object Palette, Column Operations, Double-click on the Absorber button, which looks like
this
Figure 11-1
3. Click Next, and then enter the following information as shown in Figure 11-2.
159
Figure 11-2
4. Click Next, and then enter the following information as shown in Figure 11-3. Then, click the Done…
button.
Figure 11-3
160
5. By clicking on the done button, HYSYS will bring up a window as shown in Figure 11-4.
Figure 11-4
When the column window as shown in Figure 11-4 pops up, click on the Run button located near the
bottom of the window. The red Unconverged box should turn to green Converged if all the above
procedure was followed. However, the results that are obtained at this point do not represent a true model
for our gas absorption column because the simulation was run using trays, not packing. Now, let’s see
how to replace trays with packing.
1) Packed type
From Fig. 11-4, click on the Internal tab, Figure 11-5 will appear.
161
Figure 11-5
From Internal Type drop list, change the internal to Packed as in Figure 11-5 and Figure 11-6.
From Tray/Packing Type drop list, it can be change the type of the packing, also it can be change the
remains properties of the packing from other properties of the column.
Figure 11-6
To see the detailed results, click on the View button. The results will appear as shown in Figure 11-7.
162
Figure 11-7
The other results can see from Results tab as in Figure 11-8.
Figure 11-8
Write the results required: Section Diameter (m): ___, Section Height (m): ___ , CO2 composition: ___
163
The final PFD will look like figure 11-9. Save your case as Absrber.
Figure 11-9
Note: In practice, efficiency and HETP (Height equivalent to Theoretical plate) values are estimated from
the experience gained with similar processes.
The efficiency in the absorption is low (<0.4) but in the stripping is moderate
If HETP are not available practically, it can be estimated theoretically using the following Eq.
Where: G/ and L/ are molar flowrate of vapor and liquid per unit area (kmol/m2s)
K g/ a
m = Equilibrium constant , = over all mass transfer coefficient kmol/s.m3
HETP = DT
In general, lower values of HETP are achieved with smaller-size random packings, particularly in small-diameter columns.
164
2) Tray type
Double click on the PDF of Absorber column, click on Internal tab, change the Internal Type to Trayed
and from Tray/Packed Type select Sieve as in the figure 11-10.
Figure 11-10
The details of design results cab seen by click on the View button, the window will be as in figure 11-11.
165
Change
tray type
Figure 11-11
Note: There is an error message in bottom of the design window, click on the Message tab in the above
and see how can you correct the errors.
In the first part of this section, we started with a problem to model an absorber that will absorb CO2 into
propylene carbonate in a packed column. In this section the user operated an absorber operation in
HYSYS to model an absorption process.
At the end of this section, the user was asked to use Aspen HYSYS to determine the concentration of CO2
(mole%) in the exit gas stream, the column height (m) and the column diameter (m).
Change the Solvent In flowrate from 2000 kmole/h to 2500 kmol/h. Run the simulation and see how the
column dimension and exit concentration of CO2 have changed.
166
Note: (1) In the case of the tray type column the tray spacing is determined by depending on the
diameter of the column as follows:
(2) The absorption efficiency is favored when the solvent temperature decreases, the pressure is high,
absorbent flow is high and when the molecular weight of the absorbent is low
167
11.2 Distillation
Before initialize the distillation column, it must be know if there is azeotrope between the components
that’s required to be separated.
Before initialize the distillation column, it must be know if there is azeotrope between the components
that’s required to be separated.
Open new case, then Add components: toluene, benzene, water, & ethanol
From the Object Palette, Custom Group: select Equilibrium Unit Operation , Figure 11-12.
Figure 11-12
From the Available group, add Ethanol & Water to Selected group,
168
Figure 11-13
169
Figure 11-14
Figure above shows that there is an azeotrope at concentration about 0.84 by mole for ethanol.
Figure 11-15 shows the XY diagram for Benzene mole fraction with Toluene, it can be seen that there is no
azeotrope between the benzene and toluene
170
Figure 11-15
171
11.2.1 Short Cut Distillation
The Shortcut Column distillation is only an estimate of the Column performance and is restricted
to simple refluxed Columns. For more realistic results the rigorous Column operation should be used.
This operation can provide initial estimations for most rigorous Columns.
Ex: Find the number of trays and reflux ratio of the distillation column required to separate benzene from
feed with flow rate of 850 kg/h at 25oC and 110 kPa contains toluene 60% and benzene 40% by mole. The
purity of benzene in distillate must not be less than 98% and the purity of toluene in the bottom must not
be less than 95%.
Solution:
From the Object Palette, Column Operation; Add Short Cut Distillation and fill the property view
of the column as in Figure 11-25:
Figure 11.25
172
Move to the Parameters page: Fill the parametrs page as in Figure 11-26.
Figure 11.26
The minimum reflux ratio calculated to be 1.78 as shown in figure 11.26 above
Now; the actual reflux ratio = 1.5×minimum reflux ratio = 1.5×1.78 = 2.7 ;
Enter this value in the [External Reflux Ratio] cell to calculate the No. of trays.
Click on performance tab to view the specifications of the column as in Figure 11-27:
173
Figure 11.27
174
11.2.2 Rigorous Separation Columns
Recovery of natural-gas liquids (NGL) from natural gas is quite common in natural gas processing.
Recovery is usually done to:
• Produce transportable gas (free from heavier hydrocarbons which may condense in the pipeline).
• Maximize liquid recovery (when liquid products are more valuable than gas).
HYSYS can model a wide range of different column configurations. In this simulation, an NGL Plant will
be constructed, consisting of three columns:
Figure 11-28
175
- Column Overviews (Figure 11-29 , Figure 11-30 & Figure 11-31)
Figure 11-31
176
- Defining the Simulation Basis
177
- Adding De-Methanizer (Reboiled absorber)
The De-Methanizer is modelled as a reboiled absorber operation, with two feed streams and an energy
stream feed, which represents a side heater on the column.
178
Figure 11-32
6. Supply the following information to the Pressure Estimates page. The values will be 2275 kPa and
2310 kPa, for the Top Stage Pressure and Reboiler Pressure, respectively, Figure 11-33.
179
Figure 11-33
7. Enter the temperature estimates as shown in Figure 11-34. The top stage temperature estimate will be
-87.22oC, and the reboiler temperature estimate will be 26.67oC.
Note: these temperatures are optional, e.i. it can be entered or the program will estimate it.
Figure 11-34
180
8. Click the Next button to continue.
9. For this case, no information is supplied on the last page of the Input Expert, so click the Done
button, Figure 11-35.
Figure 11-35
When you click the Done button, HYSYS will open the Column property view. Access the Monitor page
on the Design tab, Figure 11-36.
Figure 11-36
181
Before you converge the column, make sure that the specifications are as shown in Figure 11-36. You
will have to enter the value for the Ovhd Prod Rate specification. The specified value is 1338 kgmole/h.
Once this value is entered, click on Run button, the column will start running and should converge as in
Figure 11-37.
Figure 11-37
Although the column is converged, it is not always practical to have flow rate specifications. These
specifications can result in columns which cannot be converged or that produce product streams with
undesirable properties if the column feed conditions change.
An alternative approach is to specify either component fractions or component recoveries for the column
product streams.
1. Go to the Specs page on the Design tab of the Column property view as in Figure 11-38.
182
Figure 11-38
2. Click the Add button in the Column Specifications group to create a new specification, Figure 11-39.
Figure 11-39
183
3. Select Column Component Fraction from the list that appears.
Figure 11-40
Note: If the mole fraction required is not specifically 0.96 , but in range (0.955-0.965), then can be
determined this range in the Spec Type tab by give the lower and upper Values, and change the Fixed
Spec to Range in the (Fixed/Ranged Spec) as follows:
The Monitor page of the Column property view shows 0 Degrees of Freedom even though you have just
added another specification. This is due to the fact that the specification was added as an estimate, not as
an active specification.
7. Go to the Monitor page. Deactivate the Ovhd Prod Rate as an active specification and activate the
Comp Fraction specification which you created, Figure 11-41.
184
Figure 11-41
Once the column has converged, you can view the results on the Performance tab, Figure 11-42.
Figure 11-42
185
- Adding a Pump
The pump is used to move the De-Methanizer bottom product to the De-Ethanizer. Install a pump and
enter the following information:
The De-Ethanizer column is modeled as a distillation column, with 16 stages, 14 trays in the column,
plus the reboiler and condenser. It operates at a pressure of 2760 kPa. The objective of this column is to
produce a bottom product that has a ratio of ethane to propane of 0.01.
1. From the Object Palette, double-click on the Distillation Column icon and enter the following
information.
2. Click Done & then click Run button to run the column.
4. Select Column Component Ratio as the specification type and provide the following information:
5. Deactivate the Ovhd Vap Rate specification and activate the C2/C3 specification which you created.
- Adding a Valve
A valve is required to reduce the pressure of the stream DC2 Btm before it enters the final column, the
De-Propanizer.
The De-Propanizer column is represented by a distillation column consisting of 25 stages, 24 trays in the
column plus the reboiler. (Note that a total condenser does not count as a stage). It operates at 1620 kPa.
There are two process objectives for this column. One is to produce an overhead product that contains no
more than 1.50 mole percent of i-C4 and n-C4 and the second is that the concentration of propane in the
bottom product should be less than 2.0 mole percent.
187
2. Run the column.
3. Create two new Column Component Fraction specifications for the column: from the Monitor page,
click Add Spec.. button, then double click on the Column Component Fraction (do this step twice)
1.
2.
5. Activate the i-C4 & n-C4, and C3 specifications which you created.
(1) The efficiency of the tray column can be specified from Parameters tab; Efficiency page
(2) To cool the condensate product to temperature less than condensation temperature, from Design tab,
Subcooling page it can be choose either Degree of Subcool (means how many temperatures are
reducing) or subcool to (means to what temperature are cooled).
(3) It can added another section(s) to the column from Tray Sizing (After complete its procedure) ;
Design tab; Set up page.
(4) The number of flow paths can be specified from Tray Sizing; Design tab; Specs page.
(5) The dimensions of the Condenser and Reboiler can be noted by click on Column Environment button
then double click on the condenser or Reboiler; Rating tab
(6) The column can be modified (e.g. exchange the reboiler or condenser by Heat exchanger) from
parameters tab; Solver page; then select [Modified HYSIM Inside-Out] from the Solving Method Group
(7) To Know the other Inactive Specifications; in the Design tab, Monitor Page; click on Update Inactive
button. From the values appear; it can be active any specification to Run the column.
(8) To Know the Overall Material and Energy balnces; in the Home Ribbon click on the Flowsheet
Summary
(9) To Export the Data of the PFD to Excel file; Click on the Model Summary in the Home Ribbon.
189
Operation of Distillation Column
190
191
192
Ex.1: Propylene Glycol production
Propylene_Glycol is produce from the reaction of proyleneoxide and water. The two feed streams are fed
to the CSTR reactor. Figure 11-43 shows the PFD of the process.
Figure 11-43
- Click on Reaction Folder, Add Reaction Set by click on Add button, then add Reaction by click
on Add Reaction button: choose kinetic, enter the stoichoimetry according the reaction:
- HYSYS provides default values for the Forward Order and Reverse Order based on the reaction
stoichiometry. The kinetic data for this Tutorial is based on an excess of water, so the kinetics are
first order in Propylene Oxide only.
In the Fwd Order cell for H2O, change the value to 0 to reflect the excess of water.
193
Figure 11-44
Close the reaction view and Add Reaction Set to the Fluid Package by click on Add to FP button then
click on Add Set to Fluid Package button
The propylene oxide enters at 75oF , 1.1 atm, 150 lbmol/hr and 100% composition.
The two streams (PropyOxide and water) are mixed and sent to the kinetic reactor.
Add CSTR reactor, fill the design tab (Conection page) as in figure 11-45:
Figure 11-45
194
- Parameters page: the Delta P and the Volume parameters are acceptable at the default values.
- Reactions tab: Attach the Reaction Set that you created in the Basis Environment.
- Dynamic tab: volume = 280 ft3 and Liq volume percent = 85% full.
Initially the Reactor is assumed to be operating at isothermal conditions, therefore the outlet temperature
is equivalent to the feed temperature, 75°F.
Distillation Column: Double-click on the Distillation Column icon in the Object Palette.
Figure 11-46
- Click the Next button then click Next button to advance to the Pressure Profile page.
- Click the Next button to advance to the Optional Estimates page. For this example, no estimates are
required.
- Click the Next button to advance to the five and final page of the Input Expert. This page allows you to
supply values for the default column specifications that HYSYS has created.
195
In general, a Distillation Column has three default specifications. The overhead Vapour Rate and Reflux
Ratio will be used as active specifications, and later you will create a glycol purity specification to exhaust
the third degree of freedom. The third default specification, overhead Liquid Rate, will not be used.
The Flow Basis applies to the Vapour Rate, so leave it at the default of Molar.
Figure 11-47
196
Figure 11-48
Deactivate the Distillate rate specification and Activate the H2O Fraction Specification where the Degree
of Freedom will equal to 0
If the column is not Running; Click on Run button, The column will Converged
197
Ex.2: LNG (Liquefied natural gas) Exchanger (Multi stream heat exchanger)
The LNG (Liquefied Natural Gas) exchanger model solves heat and material balances for multi-stream
heat exchangers and heat-exchanger networks. The solution method can handle a wide variety of specified
and unknown variables.
For the overall exchanger, you can specify various parameters, including heat leak/heat loss, UA or
temperature approaches. Two solution approaches are employed; in the case of a single unknown, the
solution is calculated directly from an energy balance. In the case of multiple unknowns, an iterative
approach is used which attempts to determine the solution which satisfies not only the energy balance, but
any constraints, such as temperature approach or UA.
To install the LNG operation, press F12 and choose LNG or select the LNG Exchanger button in the
Object Palette .
The dynamic LNG exchanger model performs energy and material balances for a rating plate-fin type
heat-exchanger model. Figure 11-49 shows a schematic diagram of LNG exchanger.
Figure 11-49
The process is a Turbo Expander plant. LPG Product is obtained from a Feed Natural Gas stream by
cooling, expanding, separating, and distillation. The remaining gas is then re-compressed for export.
Figure 11-50 shows the PFD of the LPG Process.
198
LPG Product
Figure 11-50
N2 , CO2 , Methane, Ethane, Propane, i-Butane, n-Butane, i-Pentane, n-Pentane, and n-Hexane.
6. Select the Connections page on the design tab and enter the information as in figure 11-51: (Pressure
in kPa)
Note: Any number of Sides may be added simply by clicking the Add Side button.
199
Figure 11-51
Note: If the designated hot pass is actually cold (or vice versa), the operation will still solve properly. The
actual Hot/Cold designation (as determined by LNG) can be found on the Performance tab in the results
page.
7. Switch to the Parameter (SS) page. These values will be left unchanged.
-The Exchanger Design (Weighted) method is used to split the heating curve into intervals. The Step
Type parameter sets whether HYSYS split the curve by temperature, enthalpy, or auto interval to minimize
the errors.
- The Weighted method is required for exchangers with more than two sides.
- By selecting the Dew/Beb pt. checkboxes, HYSYS will points on the heating curve for any phase change
that occur.
9. Go to the Spec (SS) page: Degree of Freedom = 6 = (7 – 1); Add two Specifications:
200
10. Add Cooler as follows:
Name E-100
Inlet 1A
Outlet 2
Parameter: 20 kPa
Pressure drop
Outlet Temp -62oC
Energy Stream Q-100
11. Add Separator: Inlet = 2, Liquid Product = 4, Vapor product = 3
12. Add Expander: Inlet = 3, Outlet = 5, Outlet pressure = 2800 kPA, Adiabatic efficiency = 75%,
Energy Stream = QK-100
13. Add another Separator to the Expander outlet:
Inlet = 5, Vapor product = 7 , Liquid product = 6
14. Add Valve: Inlet = 4, Outlet = 8
15. Add Set operation with the following parameters: Target variable: click select variable:
Object = 8, Variable = pressure
Source: Object = 5
Parameters tab: Multiplier = 1, Offset = 0
201
Type Column Duty Type Colum Vapor Pressure Spec
Name Ovhd Duty Name Btms RVP
Energy Stream Q-102 @ Col1 Stage Reboiler
Spec Value 0.001 kJ/h Type Reid Vapor Pressure
Phase Liquid
Spec Value 200 psia
RVP: Reid vapor Pressure is a volatility measure commonly used in the refining industry. It is the
pressure at which vapor and liquid have a 4:1 volume ratio at 100oF.
Draw 9@Col1
24. Ensure the Ovhd Vap Rate and Btms RVP specs are set to be Active and the Ovhd Duty spec
Estimate.
The other unused spec can be deleted by click on the spec, click view, and then click delete.
23. Click Run if the column is not running. (What is the Duty of the condenser?)
27. Add Set operation: Target: Object = QK-101, Variable = Heat flow
202
Ex.3: Reactive Distillation
With the continuous removal of reaction products, reactive distillation has found acceptance as a means
of improving the technical and economic operation of processes where it is applicable. Advantages of this
process include higher yields, energy savings and reduced capital costs.
In this example, you will study the manufacture of Methyl Acetate, an important oil resin used in the
manufacture of artificial leathers. You will simulate the synthesis of Methyl_Acetate from Methanol and
Acetic Acid in a catalytic distillation column.
Start a new case; components: Methanol, AceticAcid, M-Acetate and Water. The fluid package is Wilson.
6. Add the Rxn Set to the Fluid Package by pressing the Add FP button,
then click on Add Set to Fluid Package
203
Add a Distillation Column with the following connections:
Go to the Monitor page of the Design tab and run the column. Complete the following table with
information from the simulation:
In this case, the reaction occurs on stages 5-10 of the column. Reaction Sets in HYSYS can be attached to
any of the stages in the column, including the Condenser, Reboiler and Side operations.
204
1. Double click on the column PFD: On the Reactions tab, press the New button.
2.
Enter the data as in figure 11-52: (From drop-down list of the Reaction Set select: Set-1)
Figure 11-52
2. Check the Active box on the Column Reaction view. The default column solver in HYSYS "HYSIM
Inside-Out" is not capable of handling reactions in the column. Therefore, HYSYS will change the solver
to Newton-Raphson Inside-Out.
205
3. However, for this column, we need to use the Sparse Continuation Solver. Select this solver on the
Solver page of the Parameters page, figure 11-53.
Figure 11-53
4. Run the column. (The column may run automatically after the solver is changed.)
5. You can view the results on the Performance tab, and the reaction results on the Reactions tab, Results
page.
Does the reaction inside the column have any effect on the compositions of the product streams?
__________
206
Ex.4: Three Phase Distillation
Due to the addition of a second liquid phase in the column, it is very difficult to model this type of
column. For this reason, we must use the Sparse Continuation Solver to solve these systems.
Property Package: NRTL, On the Binary Coefficients tab, press the Unknowns Only button.
Add a Three Phase Distillation Column operation to the simulation by double click on Three phase
Select the Distillation radio button and press the Next button, figure 11-54.
207
Figure 11-54
On the next page, set the “number of stages” = 6, and set the "Two Liquid Phase Check" on the top five
stages (condenser + 4 stages). Figure 11-55, Click next
Figure 11-55
On the third page, enter names for the streams around the condenser. For this column, the Reflux Stream
will be the Lights only; therefore, select the appropriate radio button in this group box, figure 11-56.
208
Figure 11-56
Click Next. On the fourth page, set the vapour rate to 0.0 and the Light Reflux Fraction to 1.0. The
degrees of freedom should read 0, once both of these values are entered, figure 11-57 will appear.
Figure 11-57
Click Next. On the fifth page, attach the Feed stream to stage 1. Also, define the reboiler energy stream
and attach a liquid product stream to the reboiler. From now on, the screens will look familiar to the Input
Expert screens that you have seen previously, figure 11-58.
209
Figure 11-58
Click Next. On the next page, set the condenser pressure at 101 kPa, and the reboiler pressure at
105 kPa.
On the next page, you are asked to supply optional temperature estimates. Enter values of 85, 90, and
95oC for the condenser, top stage and reboiler temperature estimates, respectively.
210
Note
The damping factor serves to reduce the amplitude of oscillations that occur in the solver. Often times
convergence can become cyclic, which can prevent the solver from finding a solution. This is where a
damping factor becomes useful. If you click the Troubleshooting icon on the ribbon under Get Started
and search for ‘damping factor’ you will see the following guidelines.
211
Ex.5: Extractive Distillation for Heptane-Toluene Separation
When the two components in a binary mixture have very close normal boiling points, their relative
volatility is likely to be small if they do not form an azeotrope. For such cases, it may be more efficient to
use extractive distillation with a solvent than normal distillation. In extractive distillation, a less volatile
solvent is used to increase the relative volatilities of the original mixtures, allowing for easier separation.
In this example, phenol is used as the solvent for the separation of n-heptane and toluene.
Figure 11-59
Add Distillation column and fill the 1st page as in figure 11-60:
212
Figure 11-60
Click Next, leave page 2, Next: Enter the pressures of condenser and reboiler = 1 bar
Leave the pages 4, click Done in the page 5
For the Feed stream enter a Vapour Fraction of 0.5, a Pressure of 1 bar, and a Molar Flow of 100
kgmole/h.
Composition: 0.5 for each component
In the Design tab; Monitor page; de-activate the Reflux ratio & Ovhd Vap Rate,
213
Figure 11-61
2. Column component fraction: fill as in figure 11-62
Figure 11-62
The degree of freedom becomes 0, Click Run, the column will converge.
214
Go to the Cond./Reboiler form under the Performance tab. Make a note of both the Condenser and
Reboiler duties. The Condenser Duty is 1.909e7 kJ/h and the Reboiler Duty is 1.909e7 kJ/h.
Create a component list. Add n-Heptane, Toluene, and Phenol to the component list.
Property Package: NRTL, and select RK as the Vapour Model.
Enter the Simulation environment, Add Distillation column
Fill the 1st page as in figure 11-63.
Figure 11-63
In the Worksheet: For the Feed stream enter a Vapour Fraction of 0.5, a Pressure of 1 bar, and a Molar
Flow of 100 kgmole/h. For the Solvent stream enter a Temperature of 181°C, a Pressure of 1 bar, and a
Molar Flow of 60 kgmole/h.
In the Design tab, Monitor page: de-activate the Reflux ratio, In the Ovhd Vap Rate cell enter 50
kmol/hr
The Degrees of Freedom for the column should now be 0. Click Run to begin calculations. The column
should solve.
Add a second distillation column to separate the solvent from the toluene in the Rich-Solvent stream.
On Page 1 enter the information as in figure 11-65 and click Next when complete.
Figure 11-65
Leave the 2nd page, Next: enter the pressure for the condenser and reboiler = 1 bar
Leave the pages 4 empty & click Done in the page 5,
216
In the Design tab, Monitor page; Uncheck the active boxes so that there are no active specifications.
The Degrees of Freedom for the column should now be 0. Click Run to begin calculations. The column
should converge.
Recycle the Lean-Solvent stream back to the first column. Add a Recycle block to the flowsheet from the
Model Palette.
Double click on the recycle block (RCY-1). Select stream Lean-Solvent as the Inlet and stream Solvent as
the Outlet. The recycle block should solve, figure 11- 68.
Figure 11-68
217
Check results. Double click on the first column (T-100) and go to the Cond./Reboiler form under the
Performance tab. Make note of the Condenser and Reboiler Duties.
Double click the second column (T-101) and go to the Cond./Reboiler form under the Performance tab.
Make a note of the Condenser and Reboiler Duties.
The following table will summarize the energy requirements from the case 1 & case 2
For the separation of n-heptane and toluene, extractive distillation has a significant advantage in total
energy requirements. Adding phenol as a solvent increased the relative volatilities of n-heptane and
toluene in the mixture and allowed for a much easier separation. However, extractive distillation required
more equipment in this case. Therefore a further analysis on capital versus operational costs would have to
be performed in order to make a decision as to which design is the better option.
218
Ex.6: Extraction of Acetone-Water with 3-Methylhexane as the Solvent
Water has a high latent heat (heat of vaporization) compared to many other components. For the
separation of a water-acetone mixture (50 wt-% each), it may be more energy efficient to use extraction
instead of direct distillation. In this example, we utilize 3-methylhexane as a solvent to remove water via
liquid-liquid extraction, followed by distillation to remove the solvent from acetone.
Figure 11-70
Click Next, enter the pressure for top and bottom = 1 bar
Next, enter the top temp = 25oC, click Done
Define the Feed stream: T = 25oC, P =1 bar, mass flow = 100kg/h,
Composition: Mass fraction = 0.5 for acetone and water
Click Run for the column.
Check the composition of the Water stream exiting the bottom of the column. You will see that the mole
fraction for water is 1.
Insert a Distillation Column ; On Page 1 enter the information as in figure 11-71 and click Next when
complete.
219
Figure 11-71
Leave page 2, Next; enter Condenser and Reboiler Pressures of 1 bar. Click Next
On Page 4 leave all fields blank. Click Done on the final page to configure the column.
In the Design tab, Monitor page; Enter 1.5 for Reflux Ratio and De-active the Ovhd Vap Rate
Add new specification; Add Column Component Fraction as in figure 11-72.
Enter 0.99999
figure 11-72
The Degrees of Freedom for the column should now be 0. Click the Run button to begin column
calculations.
220
Add a cooler to cool the bottoms stream in order to recycle it back to the mixer
Cooler: Input = Sol-Rec, Outlet = Lean-Sol, Energy = Q-Cool
Lean-Sol temp = 30oC, and its pressure = 1 bar
Recycle the bottoms streams from the second column in order to prevent throwing away acetone product.
Add a Recycle block to the flowsheet : Inlet = Lean-Sol, Outlet = Recycle Sol
Figure 11-73
221
Ex.7: Azeotropic Distillation: Production of Anhydrous Ethanol Using an Entrainer
Ethanol production via fermentation occurs in water, which must later be separated to make
anhydrous ethanol (99.95% ethanol). There is an azeotrope in the ethanol-water system at approximately
95 mol-% ethanol, which is a barrier to separation. Cyclohexane is one of the solvents used for the
production of anhydrous ethanol for food and pharmaceutical usage. It is used as an entrainer: the ternary
mixture forms a ternary azeotrope with a different ethanol concentration, which allows ethanol to enrich in
the other stream. The azeotropic liquid is separated to recover the entrainer and the ethanol that exits the
column in the azeotropic mixture.
Problem Statement
The feed to the separation train is a stream at 100 kgmole/h with 87 mol-% ethanol and 13 mol-%
water. Cyclohexane is added to the column, and > 99.95 mol-% ethanol exits the bottom of the column.
The distillate is then separated in three phase condenser. The cyclohexane-rich stream is recycled directly
to the first column, while the water- and ethanol-rich stream is sent to a second column from which
almost-pure water exits in the bottoms. The distillate of the second column is recycled to the first column.
Open new case: Add Ethanol, Water, and Cyclohexane , property package : PRSV
Enter simulation
Add 5 material streams: Feed, Feed Recycle, Solvent, Make Up, Solvent Recycle
Make Up Stream: T = 25oC, P = 1 bar, mole flow = 0.01 kmol/h, composition = 1 for cyclohexane
Solvent Recycle stream: T = 25oC, P = 1 bar, mole flow = 400 kmol/h,
composition = 0.5 mole fraction for cyclohexane & ethanol
Feed stream: This is the stream that pumps the ethanol-water mixture into the process.
Vapor fraction = 0.3, P =1 bar, mole flow = 100 kmol/h,
Composition: 0.87 for ethanol & 0.13 for water by mole
Feed Recycle stream: This stream will be the ethanol-water mixture that exits the condenser of the second
column. This stream will be fed back to the first column to prevent losses of ethanol.
In the first window that appears select the Distillation radio button. Click Next.
change the Number of Stages to 62, figure 11-76, Make sure that the Condenser is selected to check for
two liquid phases. Click Next when complete.
222
Figure 11-74
Figure 11-75
223
Click Next when complete.
In the fourth window, leave all fields blank and click Next.
The Distillation Column Input Expert window will now appear. Fill the window as in figure 11-76.
Figure 11-76
On Page 3 of the Distillation Column Input Expert enter Condenser and Reboiler Pressures of 1 bar.
Click Next when complete. Leave page 4 empty, Next,
Click Done
Enter a value of 3.5 for the Heavy Reflux Ratio and a value of 1 for the Light Reflux Ratio. First uncheck
the active box for Bot Product Rate and check the active boxes for Light Reflux Ratio and Heavy
Reflux Ratio.
224
Figure 11-77
Click Run.
Before we construct the second column, we will add an Adjust block and a Spreadsheet to find the
correct flowrate for the Make Up stream.
Double click on the spreadsheet (SPRDSHT-1). Go to the Spreadsheet tab. Enter the following text in
cells A1 and A2 as in figure 11-78
Figure 11-78
Right click on cell B1 and select Import Variable. Select the Master Comp Molar Flow of
Cyclohexane in stream ETOH.
225
Right click on cell B2 and select Import Variable, figure 11-79. Select the Molar Flow of the Make Up
stream. Having these two flow rates side by side will easily allow you to check that the amount of solvent
leaving the system is equal to the amount of solvent entering the system.
Figure 11-79
As you can see from the spreadsheet there is more solvent leaving the system than is entering. This will
cause convergence issues when we attempt to close the recycle streams. This is where we will use the
adjust block.
Double click on the adjust block (ADJ-1). Select the Adjusted Variable to be the Molar Flow of the
Make Up stream, select the Target Variable to be the Master Comp Molar Flow (Cyclohexane) of
stream ETOH, and set the Target Value to cell B2 in the spreadsheet as in figure 11-80.
Figure 11-80
226
The adjust block will vary the Make Up stream flowrate until the amount of solvent leaving the system
equals the amount entering the system. Go to the Parameters tab. Change the Tolerance to 0.001
kgmole/h and change the Step Size to 0.01 kgmole/h. Click the Start button to begin calculations. After a
few moments the flowsheet will converge.
Open the spreadsheet and you will see that the solvent leaving the system is now equal to the solvent
entering the system, figure 11-81.
Figure 11-81
Add a Recycle block to close the recycle loop for the solvent.
Figure 11-82
227
We are now ready to add the second column. Add a Distillation Column Sub-Flowsheet from the Model
Palette.
Open the Distillation Column Input Expert window and fill it as 11-83.
Figure 11-83
In the Design tab, Monitor page; de-activate the Feed Rec Rate
Add new specification: Add Column component fraction as follow
228
Click Run to begin calculations.
The last step is to connect the Feed Recycle loop. Add a Recycle block to the flowsheet.
Inlet = Feed Rec, Outlet = Feed Recycle (Wait some minutes for solving)
The final PFD should appear as in figure 11-84.
Figure 11-84
Conclusions
In this example, cyclohexane is used as the entrainer to separate water and ethanol to produce
anhydrous ethanol. By using the proper amount of solvent, we obtain pure ethanol from the bottom of the
first column. The stream from the top of the first column is separated into two streams using a three phase
condenser: One stream is solvent rich and is recycled back to the first column as solvent; the other stream
is well within another distillation region so that we can use the second column to obtain pure water. The
top stream of the second column is recycled back to the first column as feed.
229
Chapter 12
Assay Setup
When setting up a refining case in Aspen HYSYS Petroleum Refining, the first step is to select and
characterize crude assays. This is done within the HYSYS Properties Environment and will determine the
characteristics of the incoming crude feeds before they can be used in the Simulation Environment.
To import the petroleum assay, using the “Petroleum Assays” button, in the Home Ribbon as shown in
Figure 12-1
Figure 12-1
By opening this page, you will get access to the Assay Management ribbon as shown in Figure 12-2.
Figure 12-2
230
Here, you can find the “New Assay” button in the far left. Click on the dropdown arrow to select
between these options: Figure 12-3
Figure 12-3
• Import from File – Import a supported Excel, CSV (comma-separated values file), or AFAM case
with assay data
• Manually Enter – Manually enter data into the Assay Table
By default, clicking on the “New Assay” button will select the “Import from Library” option.
Once the Assay Library has opened, a new window is open as shown below that is required to
choose one of the common assay component lists from the drop down list. From drop down choose
"Assay Components Fahrenheit to 1500F", click OK.
There will be a large number of assays to select from and multiple ways to filter the results. In this
example, we will import an existing assay from the Assay Library called Bachaquero crude. To find the
specific crude to use in this case, type the name “Bachaquero” under “Assay name” and select the correct
crude from the list, which will then be highlighted in yellow. To finalize the selection, click “OK” in the
bottom-right corner, as shown in Figure 12-4.
231
1
Figure 12-4
After the assay has been imported, you can verify stream properties by checking the “Input Assay” and
“Conventional Results” plys in the left-hand Navigation Pane under the Assay name, Fig. 12-5.
Figure 12-5
232
Manually Creating Assays
Since there are a variety of ways that assay information can be entered into a refining model,
Aspen HYSYS Petroleum Refining provides a way for users to manually enter assay data by choosing the
Manually Enter option when creating a new assay. This will open up the New Assay window, giving a
number of custom options to format your assay table, as shown in Figure 12-6.
Figure 12-6
• Multi Cut Properties– This is ideal if you have data for specific assay cuts that you wish to enter, such as
Cut Yields. Here, you can set the number of cuts from your data and enter the boiling point ranges in the
table provided.
• Single-stream Properties– This option is ideal if you want to evaluate an assay from distillation data and
bulk (whole crude) properties to create a single-stream for your model.
• Backblending– This will open a custom form designed specifically for users who have data from product
streams and wish to simulate a backblended feed by entering pure component, distillation, and property
data for specific products.
Characterizing Assays
Under the Assay Management tab, there are also options for different assay characterizations using the
“Characterize”button. The options include matching assay cut properties and whole crude properties, and
molecular characterization, as shown in Figure 12-7.
233
The default option to characterize when clicking the Characterize button will be to match the Whole
Crude properties given in the Input Assay table, which means matching the bulk values that have been
given by the user will take precedence over matching values for specific cuts.
Molecular Characterization
When Molecular Characterization is selected, the user is given a list of hydrocarbon molecules present
in the crude and is shown three editable columns: weight percent (Wt %), and the alpha and beta
coefficients, Figure 12-8. In addition, viscosity coefficients are also available. These three variables can be
modified prior to the molecular characterization step.
Once they are selected, you can click the “Run Evaluation” button to characterize the assay using this
data. After the assay has been characterized using Molecular Characterization, the user can change the
characterization method by unchecking the “Molecular Characterization” box and using the
“Characterize” button.
Figure 12-8
234
The properties that can be plotted for assays are shown under the “Assay management” tab, in the “Plot
Gallery” section, Figure 12-9
Figure 12-9
In this example, we have plotted a distillation curve for the Bachaquero crude, as shown in Figure 12-10.
You can change the graph’s coordinates, type, legend, text, etc. from the ribbon.
Figure 12-10
To compare data from more than one crude assay, the second crude must be added
Again from Home Ribbon click on Petroleum Assay, click Add button, and in the Assay name place:
Shengli, under the Selected Assay; highlight the Shengli and click OK.
In this example, we will compare Bachaquero and Shengli crudes. After Shengli has been added to the
Petroleum Assays list, you can add it to the plot: from Plot Gallery; click on Distillation, and then clicking
“Select Assays” under the Format tab when the Plot is selected, as shown in Figure 12-11.
235
Figure 12-11
Assays can be exported from a HYSYS Petroleum Refining case while in the Properties Environment.
Go into the Petroleum Assays node to access the Assay Summary, select the assays to export, and click
the Export button as shownin Figure 12-12.
Figure 12-12
236
Attaching an Existing Assay
When an assay has been created and characterized in the Properties Environment, it can be immediately
used for modeling in the Simulation Environment. To bring the assay into the simulation, it must be either
attached to a stream or added to a Petroleum Feeder unit operation.
First, add a material stream to the flowsheet using the model palette. Then, open up the stream view by
double-clicking the newly created stream and click on the “Petroleum Assay” option in the Worksheet tab.
Once in the Assay form, select the Attach Existing option which will show a dropdown list where the
previously created assay can be selected as shown in Figure 12-13. Once a characterized assay is attached,
it will define the composition and the properties of the stream so that the user just has to enter the
conditions of the stream for it to be fully defined.
Figure 12-13
If more than one assay will be used or if the assay will be mixed with another stream, they can be added
to a Petroleum Feeder, shown in the palette. In the Feeder, you can select input feed streams and feed
assays (as well as mixing specifications including a fluid package) and Aspen HYSYS will calculate a
product stream that can be used in the flowsheet.
In some cases, creating a new assay on a stream may be necessary. To do so, add a stream to the
simulation and proceed directly to the “Petroleum Assay” option in the stream view. From this form, there
are a number of options to add assays, similar to when creating an assay in the Properties Environment.
The user can select Attach Existing, Import from File, Import from Library, or if they click the “Create”
button, it will open the Manual Entry dialog for the user to begin entering assay data for the stream. Once
the assay has been created on the stream, click the “View Assay” button at the bottom of the stream view
237
at any time to open up the assay table for that specific stream. Figure 12-14 shows the view of the assay
table while in the Simulation Environment.
Figure 12-14
238
Ex: Refinery of Crude Oil
This example models a crude oil processing facility consisting of a pre-fractionation train used to heat
the crude liquids, and an atmospheric crude column to fractionate the crude into its straight run products.
The Main Flowsheet for this process appears in the following figure.
Preheated crude (from a preheat train) is fed to the pre-flash drum, modeled as a Separator, where
vapors are separated from the crude liquids. The liquids are then heated to 650°F in the crude furnace,
modeled as a Heater. The pre-flash vapors bypass the furnace and are re-combined, using a Mixer, with
the hot crude stream. The combined stream is then fed to the atmospheric crude column for separation.
Click OK
239
3- Click on the Oil Manager Folder located in the Navigation Pane, then click on Arrow beside Oil
Manager then click on Input Assay.
2 1
3
4- On Assay tab click Add.. button to create and view a new Assay. From Assay Data Type drop-
down list select TBP.
240
5- From the Light Ends drop-down list, select Input Composition.
6- In the Input Data group, select the Light Ends radio button. Ensure that Liquid Volume% is
selected in the Light Ends Basis drop-down list.
7- Fill the compositions in the following table as follows:
Note: Before entering any of the assay data, you must activate the molecular weight, density, and viscosity
curves by choosing appropriate curve types in the Assay Definition group. Currently, these three curves
are not used.
8- From the Bulk Properties drop-down list, select Used. A new radio button labeled Bulk Props
appears in the Input Data group.
9- From Molecular Wt. Curve drop-down list, select Dependent. A new radio button labeled
Molecular Wt appears in the Input Data group.
10- From the Density Curve and Viscosity Curves drop-down lists, select Independent as the curve
type. For Viscosity, two radio buttons appear as HYSYS allows you to input viscosity assay data at
two temperatures.
Entering Bulk Property Data: Select the Bulk Props radio button, and the bulk property table
appears to the right of the radio buttons. Fill the Molecular weight and Standard Density in the table as
follows:
Note: HYSYS displays a message under the matrix, stating that ‘At least 5 points are required’ before
the assay can be calculated.
2- From the Assay Basis drop-down list, select Liquid Volume. Click the Edit Assay button. The
Assay Input Table property view appears.
3- Fill the table as follows:
242
Click OK button to return to the Assay property view.
Click OK .
Using the same method as for the previous assays, Click Edit Assay, Enter the API gravity curve data
as shown here.
243
Entering Viscosity Data:
1- Select the Viscosity 1 radio button. In the Viscosity Type drop-down list above the assay matrix,
ensure Dynamic is selected.
2- In the Viscosity Curves group, select the Use Both radio buttons. The Temperature field is for
each of the two viscosity curves.
3- Click the Edit Assay button to access the Assay Input Table.
4- Input the Viscosity 1 assay data as shown here. This viscosity curve corresponds to Temperature 1,
100°F.
7- Click the Calculate button at the bottom of the Assay property view. HYSYS calculates the
Assay, and the status message at the bottom of the property view changes to Assay Was
Calculated.
8- Click the Working Curves tab of the Assay property view to view the calculated results. HYSYS
has calculated 50 points for each of the Assay Working Curves.
9- To view the Assay data you input in a graphical format, click the Plots tab. By default, HYSYS
plots the Distillation (TBP) data. You can view any of the other input curves by selecting the
appropriate variable in the Property drop-down list.
244
Cutting the Assay (Creating the Blend) into individual petroleum hypocomponents.
1- Click the Output Blend folder in the Navigation Pane. Click on Add button.
2- Click on Add…> button. There are two results:
• The Assay is transferred to the Oil Flow Information table. (When you have only one Assay, there is
no need to enter a Flow Rate in this table.)
• A Blend (Cut) is automatically calculated based on the current Cut Option.
HYSYS calculated the Blend based on the following default values for the boiling point ranges and
number of cuts per range:
• IBP to 800°F: 25°F per cut, generating [(800-IBP)/25] hypocomponents
• 800 to 1200°F: 50°F per cut, generating 8 hypocomponents
• 1200 to 1400°F: 100°F per cut, generating 2 hypocomponents
The IBP, or initial boiling point, is the starting point for the first temperature range. The IBP is the
normal boiling point (NBP) of the heaviest component in the Light Ends, in this case n-Pentane at 96.9°F.
The first range results in the generation of (800-96.9)/25 = 28 hypocomponents. All the cut ranges
together result in a total of 28+8+2 = 38 hypocomponents.
3- Click the Tables tab to view the calculated properties of these hypocomponents.
4- Return to the Data tab. From the Cut Option Selection drop-down list, select User Ranges.
5- In the Starting Cut Point field, enter 100°F. This is the starting point for the first range.
6- Enter the remaining cut point temperatures and the number of cuts values as shown in the figure below.
245
7- Click the Submit button to calculate the Blend based on the current initial cut point and range values.
8- Click the Tables tab to view the properties of the petroleum hypocomponents.
To view the distribution data, click Table tab, select Oil Distributions from the Table Type drop-down
list.
At the bottom of the Cut Input Information group, the Straight Run radio button is selected, and HYSYS
provides default TBP cut point temperatures for each Straight Run product.
1- Click the Install Oil button. In the Stream Name column Type the name Preheat Crude, press Install.
The stream Preheat Crude was already installed during the oil characterization procedure.
Double click on the stream and enter the date: T = 450oF, P = 75 psia, Std ideal Liquid Volume Flow =
100000 barrel/day.
1- Bottom Steam: T = 375oF, P = 150 psia, Mass Flow = 7500 lb/hr, composition: 100% H2O
2- Diesel Steam: T = 300oF, P = 50 psia , Mass Flow = 3000 lb/hr, composition: 100% H2O
3- AGO Steam: T = 300oF, P = 50 psia , Mass Flow = 2500 lb/hr, composition: 100% H2O
1- Add Separator: Input = Preheat Crude, Vapor Outlet = PreFlashVap, Liquid Outlet = PreFlashLiq
2- Add Furnace (simulated as a heater): Inlet: PreFlashLiq, Outlet: Hot Crude, Energy: Crude Duty
Parameter Page: Delta P: 10 psia , from worksheet , enter the temp of Hot Crude = 650oF
3- Add Mixer: Input: PreFlashVap & Hot Crude, OutPut: Tower Feed
7- Next, Condenser pressure = 19.7psia, Condenser pressure drop = 9 psi, Bottom stage pressure =
32.7psia.
8- Next, Condenser Temp = 100oF, Top stage Temp = 250oF, Bottom Stage Temp = 700oF.
9- Next, vapor rate = 0, Reflux ratio = 1. Click Done
10- Adding Specification Values: Select the Monitor page, In the Column Specifications group, delete
Reflux Rate and Btms Prod Rate Specifications.
11- Click the Add Spec button. Select Column Draw Rate and click Add Spec(s) button.
In the Name cell, type WasteH2O Rate.
No further information is required, Close the property view.
12- Select Connection Page, Click Side Ops, on the Side Stripper page, click Add button, The
completed Side Stripper property view Must appear as shown below. Click Install button, close.
247
13- Use the previous step (step 12) to install the two remaining side strippers DieselSS and AGOSS.
The completed DieselSS and AGOSS side stripper property views appear in the following figures.
248
14- Installing Pump Arounds: select the Pump Arounds page. Click the Add button.
15- In the Return Stage drop-down list, select stage (1_Main TS). In the Draw Stage drop-down list,
select stage (2_Main TS).
16- Click the Install button, and a more detailed Pump Around property view appears.
17- Double click on [PA_1_Rate(Pa)], Change Flow Basis to Std Ideal Vol. Enter Spec Value 50000
barrel/day, close.
50000 barrel/day
18- Double-click in the blank space under the PA_1_Dt(Pa) specification, In the Spec Type drop-
down list, select Duty and enter the spec value (-55e6 Btu/hr).
Close icon to return to Pump Around, Close icon of Pump around.
19- Repeat the previous steps (14-18) to install the two remaining pump around.
249
The summary of previous steps are shown in table below:
20- Select the Connection page on Design tab. (Note: it is not required to export the pump around
duty streams PA_1_Q, PA_2_Q, and PA_3_Q to the Main Flowsheet, so their External Stream
cells remain undefined.)
21- On the Design tab, Specs page: In the Column Specification group, click Add, Select Column
Liquid Flow, click Add Spec(s), Change the name to Overflash, Fill the variables as in figure
below
250
A typical range for the Overflash rate = 3 – 5% (3.5%)
22- Click Add again, Double click on Column Duty , Fill the variables as in figure below
If the column begins to run on its own before you click the Run button, click the Stop button and
continue activating & deactivating specifications.
23- In the Monitor page, double click on Distillate Rate, change flow basis to
24- Activate the Overflash , Kero Reb Duty, and Vap Prod Rate Specifications
25- Deactivate the Reflux Ratio, Waste H2O Rate, and KeroSS BoilUp Ratio Specifications.
251
Ex.2 Refinery of Crude oil 2
Process Description
After passing through a preheat train, 100,000 barrel/day of 29.32o API crude is fed into a pre-
flash separator operating at 450oF and 75 psia. The vapor from this separator bypasses the crude furnace
and is re-mixed with the hot (650oF) pre-flash liquids leaving the furnace. The combined stream is then
fed to the atmospheric crude column.
The column operates with a total condenser, three coupled side strippers, and three pump around circuits.
A naphtha product is produced overhead, a kerosene product is produced from the first side
stripper, a diesel product is produced from the second side stripper, and an atmospheric gas oil (AGO) is
produced from the third side stripper. Both the AGO side stripper and the diesel side stripper are ‘steam
stripped’, while the kerosene side stripper has a reboiler.
Click OK
252
Oil Characterization
1- Click on the Oil Manager Folder located in the Navigation Pane, then click on Arrow beside
Oil Manager then click on Input Assay.
2 1
3
1- On Assay tab click Add button to create and view a new Assay. From Assay Data Type drop-
down list select TBP.
253
2. Change the Bulk Properties setting to Used.
3. Select Bulk Props radio, Complete the Input data for the as shown below:
4. Click the Calculation Defaults tab. In the Extrapolation Methods group, select Lagrange for each
method using the drop-down lists.
Note: HYSYS displays a message under the matrix, stating that ‘At least 5 points are required’ before
the assay can be calculated.
254
Click OK
7. In the Assay Definition group, from Light Ends drop-down list, select Input Composition.
8. In the Input Data group, click the Light Ends radio button.
9. Enter the light ends data as follows.
255
10. Upon completion of characterizing the assay, click the Calculate button. HYSYS will calculate the
Working Curves, which can be viewed on the Working Curve tab.
You can scroll through this table to view all 50 points of the Working Curve.
1. Click the Output Blend Folder in the Navigation Pane (Oil Characterization property view) and click
the Add button. The Blend: Blend-1 view appears.
2. Click the Add…> button. HYSYS will transfer that Assay to the Oil Flow Information table.
As a guideline, each Outlet stream from the crude column should contain a minimum of 5
hypocomponents where the composition is greater than 1.0%. Therefore, a total of 30 components should
fulfil this requirement.
3. From the Cut Option Selection drop-down list, select User Points, then specify the Number of Cuts at
30. HYSYS will calculate the hypocomponents.
1. Click the Install Oil button of the Oil Characterization property view.
2. In the Stream Name cell, type Raw Crude. This is the stream name where you would like to “install”
the oil. Click Install button
3. Enter Simulation Environment
Inlet Stream
Because the composition has been transferred from the Oil Characterization, the stream is automatically
flashed.
Pre-Flash Operations
Install the Separator, Heater, and Mixer and provide the information displayed below:
Separator: Inlet = Raw Crude, Vapor Outlet = PreFlash Vap, Liquid Outlet = PreFlash Liq
Heater: Inlet =PreFlash Liq, Outlet = Hot Crude, Energy = Crude Duty
Parameter: Delta P = 10 psi
Temp of Hot Crude = 650oF
256
Mixer: Inlet = Hot Crude & PreFlash Vap, Outlet = Crude Charge
Before simulating the atmospheric crude tower, the steam feeds and the energy stream (Q-Trim -
representing the side exchanger on stage 28) to the column must be defined.
Install Energy Stream, Name: Q-Trim which does not require any specifications, this will be calculated
by the Column.
Column
For this example, install the 3 Stripper Crude Column custom templates.
1. From the Object Palette, Columns group: Double click on the Blank Column Sub-flowsheet
icon, then click the Read an Existing Column Template button. The Available Column Templates finder
property view appears.
2. In the Files of type drop-down list, select Column Templates (*.col).
3. From the list, select the 3sscrude template file, then click the Open button. The 3sscrude template
installed 40 trays, 29 in the Main Tray section, 3 trays in each of the 3 Side Strippers (1 reboiled and 2
steam stripped), a reboiler, and a condenser.
4. In the Column Property view, connect the Inlet and Outlet streams to the column sub-flowsheet as
shown below (Design tab, Connections page).
257
5. Modify the Draw and Return stages of the Side Strippers and Pump Arounds on the corresponding
page of the Side Ops tab.
6. In the Atmos Tower Column property view, specify the column information below.
258
Specifications
On the Monitor page of the Design tab, make the following changes and input the values into the default
set of specifications supplied with the pre-built 3-Side Stripper Column.
1. Delete the Kero SS BoilUp Ratio and the Residue Rate specs. Open the specification property view
by clicking the View button, then click Delete to delete the specification.
2. Specify the Reflux Ratio spec to have a value of 1. Deactivate the Reflux Ratio
3. Change the following default specifications by selecting the specification in the table and clicking the
View button.
Change the Flow Basis to Std Ideal Volume before entering values.
Change all the Pump Around delta T specifications to a Duty specification.
5. On the Monitor page of the Design tab, add the following new specifications by clicking the Add Spec
button in Column Specifications group.
259
The final specification table will appear as shown below:
6. Once you have provided all of the specifications, Run the column.
260
Chapter 13
Optimizer
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.
Ex: In this example, a simple distillation column to separate Tetrahydrofuran (THF) from Toluene is
simulated. The object of this 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.
1. Open New Case: components: Tetrahydrofuran and Toluene, Fluid Package: Wilson.
2. On the Binary Coeffs tab, press the Unknowns Only button.
3. Enter the Simulation Environment.
4. Add a material stream with the following values:
Name: Feed
Temp = 10oC
Pressure = 140 kPa
Mass Flow = 3700 kg/h
Composition-Mass Fraction: THF = 0.44, Toluene = 0.56
Next
261
7. Create two new Column Component Fraction specifications as shown in figure 13-1:
Figure 13-1
8. These two specifications should be Active. Deactivate the Reflux ratio & Distillate rate
The column should now solve automatically. If it does not, press the Run button to start the solver.
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.
Typically, optimization studies involve an ‘economic model’ composed of Object function or Profit
Function and operating constraints.
262
General a profit function is given by:
1. In the Variable tab Press the Add button to add the first variable, Object = T-100, Variable = Spec
Value, Variable Specification = THF Purity Spec, OK, Figure 13-2.
Figure 13-2
2. Add the second variable, Object = T-100, Variable = Spec Value, Variable Specification = Toluene
Purity Spec, OK
3. Set the Lower Bound at 0.9, and the Upper Bound at 0.99 for each variable, figure 13-3.
Figure 13-3
263
4. Press the Spreadsheet button on the Optimizer view to open the Spreadsheet.
5. Move to the Spreadsheet tab.
To Import any Variable Right click in the cell and select Import variable.
Enter the following constants on the spreadsheet. The comments can be added, if desired.
Note: All common functions must be preceded by a (+ symbol). Special Functions must be preceded by the
@ symbol.
264
Note: the prices of the products increase as the purities increase.
The Objective Function is placed in Cell D8, figure 13-4. The equation is:
+(b3*d4+b5*d6)-b7*d7-(b1*d1+b2*d2)/3600
Figure 13-4
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.
265
Figure 13-5
Parameters tab
The Parameters tab is used for selecting the Optimization Scheme.
• 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.
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. For more information on
the results, record the results:
266
Original Optimizer
Ex: Create the following sample case of multiple heat exchangers (Figure 13-6) to optimize the overall
UA by using the Original Optimizer.
Figure 13-6
267
Outlet = E-100 Out
Stream Specifications
• Flowrate of E-101 Feed = Optimized variable (Initially set to the flow rate of 1670 lbmole/hr)
Mixer: Inlet: E-100 Out, and E-102 Out Outlet: Mixer Out
Optimizing Overall UA
268
The Optimizer determines the optimum Tee flow ratio such that the Overall UA is minimized.
1. Click the Add button in the Variables tab to display the Variable Navigator.
3. Specify the Low and High Bounds (lbmol/hr) as shown in the figure below.
1. Click the SpreadSheet button to display the Optimizer Spreadsheet property view.
2. Click the Connections tab and then click the Add Import button to display the Variable Navigator
property view.
4. Repeat steps#2 and#3 for the Heat Exchangers E-102 and E-103.
+a1+a2+a3. This sums the UAs. In cell A5, enter 0.0. This is used in the constraints.
269
Figure 13-7
You must define the Objective Function and the Constraint Functions. The Objective Function is the
expression being minimized, which in this case is the sum of the Heat Exchanger UAs.
2. Click the Cell drop-down list and select A4. The value of the cell is displayed in the Current Value
field.
Enter constraint functions to ensure the solution is reasonable. Each Heat Exchanger UA must be
greater than zero.
1. Click the Add button three times to add three constraints to the table.
2. In the LHS Cell drop-down list, select the cell A1, A2, and A3 for each of the respective constraints.
3. In the RHS Cell drop-down list, select the cell A5 for each of the constraints, figure 13-8.
270
Figure 13-8
4. Click the Parameters tab. For this example, use the Mixed method, leaving all the parameters at their
defaults.
5. Click the Start button, and then click the Monitor tab to watch the progress of the Optimizer.
An optimum molar flow of 1800 lbmole/hr is obtained for the stream E-101 Feed, corresponding to an
overall UA of about 1.43e5 Btu/F-hr. This compares to the specified value of 1.5e5 Btu/F-hr in the first
part of this example.
271
Chapter 14
The plant takes two different natural gas streams containing carbon dioxide and
(methane through n-butane) and combines and processes them in a simple refrigeration
system. A series of separators and coolers removes the heavier hydrocarbon components
from the natural gas stream, allowing it to meet a pipeline dew point specification. The
heavier liquid component of the gas stream is processed in a depropanizer column, yielding
a liquid product with a specified propane content.
272
- (The mass fraction can be converted to mole fraction using basis button).
- To get the phase diagram and find the relationships between P-T, P-V, V-S…etc. in the worksheet
press attachment button, then utilities, then create and double click on Envelope utility, then
performance
Mixer
Name Mix-100
Inlet Feed1, Feed2
Outlet Mixer out
Separator
Name Inlet sep.
Inlet Mixer out
vapor Out SepVap
Liquid Out SepLiq
Heat Exchanger
Name Gas/Gas
Tube Inlet SepVap
Tube Outlet Cool gas
Shell inlet LTS vap (LTS: low temperature separator)
shell Out Sale gas
- From parameters page, delta P = 10 psi (for shell side and tube side).
- From Rating tab, sizing page: change tube passes per shell to 1, then press OK.
Chiller (cooler)
Name Chiller
Inlet Cool gas
Outlet Cold gas
Energy C3 Duty
- From parameters, delta P = 10 psi.
Note: degree of freedom in the heat exchanger = 2, therefore determine temperature for cold gas as a
default value (0oF) from worksheet.
273
- Return to Gas/Gas heat exchanger by using design tab choose specs page, notice that degree of
freedom = 1. choose add button to display exch. Spec. and change the name to hot side approach, in
the type cell, enter delta temp. In stream (+) cell choose Sep. vap and in stream (-) cell choose sale
Gas and enter 10oF in spec value.
To see the results, open worksheet and notice that temperature of cool Gas = 42.9o F. Also choose
performance tab to view results. Sep Vap
10oF Cool Gas
Hot Side
Sale Gas
LTS Vap
Checking Sale Gas Dew Point:
To ensure that sale Gas is not converted to liquid during the transferring operation, we will check
the dew point temperature using the Balance where the transferring operation line specifications are
15oF & 800 psia while the sale Gas specifications are 50oF & 570 psia therefore we will create a
default line with 15oF & 800 psia.
- In the palette, double click on Balance icon, choose inlet: sale Gas and outlet: sales DP.
- In parameters choose Component Mole flow.
- In worksheet, enter the sales DP vapor fraction as 1.0 and pressure as 800 psia, HYSYS will calculate
temperature for sales DP to be 5.27oF.
This figure can be gained from attachments for the line sales Gas then envelope utilities to see the
phase diagram.
2nd Mixer
Name Mix-101
Inlet LTS liq & Sep. liq
Outlet Tower Feed
274
De-propanizer column (distillation column)
Name De-propanizer
No. of Stages 10
Inlet Tower Feed
Inlet stage 5
Condenser Full Reflux
Top Outlet Ovhd
Bottom Liquid outlet Liquid prod.
Re-boiler Energy Stream Reb Duty
Condenser Energy Stream Cond Duty
Click next & then click Next
Click next
(This page is optional, i.e. it can be left)
275
Deg. of freedom will be = 0
- Click Run. (converged: green color)
Click performance tab and column profile page for more details.
- Adjusted variable: click select var., object: Cold gas, variable: temp & OK
- Target variable: click select var., object: Sale DP, variable: temp & OK
- In the specified target value enter 10oF
- Click parameters tab, in tolerance: 0.1oF & in step size: 5oF.
- Click monitor tab & start button
→ adjust var. = 4.4oF (it was 0oF and now in the dew point it is 10oF)
- See the results from the workbook (save this case “Gas processing”)
Dynamic simulation
When the steady state is finished, we can start the dynamic state but if we wished to start a
dynamic state, we must build a steady state at first as the previous one.
Note:
من األفضل في المبادالت الحراية ان نثبت معدل الجريان ونلغي تثبيت هبوط الظغط (بعد احتسابه في الحالة المستقرة) من حقل.1
الديناميك
) كذلك في المضخات من األفضل ان نثبت البور ونلغي هبوط الضغط (بعد احتسابه في الحالة المستقر.2
في نهاية المسارات النهائية يفضل الغاء معدل الجريان وابقاء (تفعيل مقدار الضغط) من حقل الديناميك او وضع صمامات "هذه الخطوة.3
حيث سيقومMake change button بعد الضغط علىDynamic Assistance البرنامج ممكن ان يفعلها اذا تم نسيانها من خالل
بالتغيير بالنسبة للمواضيع التي مؤشر عليها عالمة صح
- Double click on Feed2 stream in PFD and delete pressure (it will be calculated by dynamic mode).
- Double click on Mix-100: from parameters page select Equalize All and close Mix-100 view
- Insert Sep. valve between Mix-101 & Sepliq stream,
(Double click on Mix-101, Delete Streams: Sepliq, LTS Liq, and TowerFeed)
276
- Added Valve.
Name Sep valve
Inlet Sepliq
Outlet Sepexit
∆p 25psi
277
- Double click on tower feed stream & in temp cell type = 24.73oF.
- Double click on Mix-101: Inlet = LTS Exit, and SepExit, Outlet = Tower Inlet
- double click on de-propanizer column & click column environment (to view details)
- Cut the LiqProd stream from the Reboiler
- After adding the valve click on Run column solver icon in the tools bar.
- It will remove some unit operations that is not needed in the dynamic state,
- Double click on cold gas icon & record T= 4.44oF (this temp. will be controlled in the dynamic state).
- In PFD, delete ADJ-1 , Dew point logical operation (Balance) & sales DP stream.
278
Column sizing
في الحالة المستقرة هبوط الضغط يحدد من قبل المستخدم لكن في الحالة الديناميكية يحسب بأستخدام الحسابات الهيدروليكية الديناميكية
- Double click on the de-propanizer column, From internal type click Add New, Select the Trayed
from Internal Type, and select Valve from Tray/Packing Type.
- Click View button, select results tab to confirm the following parameters:
- Click Rating tab, select sizing page, and enter the parameters as:
Diameter 3 ft.
Tray spacing 24 in
Dynamics Assistant
- Degree of freedom for the flow sheet must be zero, the following sizing parameters must be specified for
these unit operations:
Valves CV valve
Vessels Volume
Cooler / H.E R-valve
- For each valve in PFD: Double click on valve; Rating tab, click Size Valve (Do this for each valve in
PFD)
- Double click on De-Propanizer column click column Environment button, Double click on Re-boil
Valve, Rating tab, click Size Valve, Return to the PFD
The process variable is the actual plant value that the controller
is trying to regulate. The output target object (sometimes referred
to as the manipulated variable), is the parameter that the
controller will change to adjust the process variable to the
set point.
Level control:
Double click on each separator vessel and from dynamic tab; click Add/Configure Level controller
Flow control:
Temperature control:
Stage 9 TC controller:
- Add PID controller,
name: Stage 9 TS; process variable, Flowsheet: De-Propanizer, Object: main Tower,
Variable: stage temp, 9-main Tower, OK
Output target: Flowsheet: De-Propanizer, Object: Reboil valve, Variable: Actuator Desired position
281
- Click parameters tab; Action: direct, Kc =2, Ti = 5min, PV min. = 110oF, PV max.= 260oF,
Mode: auto.
Monitoring in Dynamics
- Double click on each Controller: StripChart tab; from Variable Set select: SP, PV, OP only
- click Create stripchart button, then click Display button
Or it can click on Monitor tab
- Click on Dynamic Ribbon
Then click on Dynamic Mode , yes, then click Make change button
Click on icon
Double Click on any controller and click on Monitor and click Face plate to see the results
In the any graph; right click and select Graph control: to show all variables; In the Axes tab; Axis
Display section; Check the Show All box
1) It must be placed a resistance operation between every pressure node (e.g. valves, pumps, compressor
…etc) to make pressure drop between unit operations to run a dynamic simulation
2) It is recommended to add a valve to all boundary streams, or alternative the pressures in boundary
streams
282
3) In the mixer, it is recommended to set to Equalize all in the parameters page.
4) To make it easier to determine where you need to make pressure and flow specifications, you can
change the color scheme of the flowsheet. In the Display Options group on the Flowsheet/Modify tab,
select the “Dynamic P/F Specs”
5) In order to run dynamic simulations, it is important to have all of your equipment sized properly. There
are two places where you can change the size of vessels: the Rating tab and the Dynamics tab on the
vessel’s property form.
Additionally, you can use the Dynamics Assistant to auto-size your equipment. On the unknown sizes tab,
select volume on the left. Change the residence times for the vessels you want to size and click the Quick
size… button.
283
To review individual changes beforehand, double click any of the items on the list to redirect to the
appropriate tab.
Chapter 15
Economic Evaluation
Open the case of Gas processing (Steady state)
284
1. Enter the buying and selling prices of the feed and product streams. Double click on Feed1 stream and
go to the Cost Parameters form under the Worksheet tab. Select Mass Flow for Flow Basis and enter
0.2 Cost/kg for Cost Factor. By the same way; price of
Feed2 = 0.2 Cost/kg
Sale gas = 0.8 Cost/kg
Liq-Prod = 0.8 Cost/kg
Ovhd = 0.5 Cost/kg
2. To view the total stream costs, go to the Economics tab in the ribbon and select Stream Price.
3. It will estimate costs for utilities. Double click on energy stream C3 Duty. From Utility Type Select
Refrigerant 1.
4. To view the utility summary, click Flowsheet Summary in the Home tab of the ribbon.
Go to the Utility Summary tab. Here you can view the cost of each utility and the total costs of utilities.
The Total Costs of Hot Utilities are $ 5.841/hr, and the Total Costs of Cold Utilities are $8.504/hr.
285
5. The operating profit of this process is equal to:
6. Go to the Economics tab, and select Economics Active. This will enable the Economic Analysis
functionality in Aspen HYSYS.
Note: If the Economic Active box is not checked for any reason, click on Activated Economic Analysis icon and
wait for some minutes
286
When the economic analysis is Activated, click the Map button.
7. When the Map button is clicked, the following window will appear. Press OK to continue.
8. The following window titled Map Preview will allow you to change the mapping for certain unit
operations.
9. Select Heater and click the drop down menu under Equipment Type.
10. A new window will appear, select Heat exchangers, heaters and press OK.
11. From the list appear, select Fixed tube sheet shell and tube exchanger, OK.
Gas/Gas(HEATX) → Heat exchangers, heaters → Heat Exchanger → U-Tube shell and tube exchanger
287
Reboiler_@T-100(HEATER) → Heat exchangers, heaters → Reboiler → Kettle type reboiler with
floating head
13. Next, click on Size icon. The sizing process will complete.
The Economic Evaluation Equipment Summary Grid will open. Go to the Equipment tab.
15. We are ready to evaluate. Click the Evaluate button in the ribbon. The economic engine will perform
the analysis, it may take a few moments.
16. Click View Equipment, and go to the Equipment tab to view any errors that occurred during
evaluation.
These errors will tell you what inputs or changes are required in order to cost the simulation more
realistically.
288
Chapter 16
Solid Operation:
16.1 Cyclone
The Cyclone is used to separate solids from a gas stream and is recommended only for particle sizes
greater than 5 microns. The Cyclone consists of a vertical cylinder with a conical bottom, a rectangular
inlet near the top, and an outlet for solids at the bottom of the cone. It is the centrifugal force developed in
the vortex which moves the particles toward the wall. Particles which reach the wall, slide down the cone,
and so become separated from the gas stream. The solids being separated must be previously specified and
installed as components in the stream attached to this operation.
Figure 16-1
Ex: Stream of air at 25oC and 150 kPa with flow of 150 kmol/h contains 1.5% by weight carbon as
pollutant. It must be removed 95% of the carbon by cyclone. What are the dimensions of cyclone required
for this purpose and what is the number of cyclones required if the maximum diameter of cyclone must
not to be exceeding from 3 m.
Solution:
1- Open new case with above components and with Fluid Package: PRSV.
2- Enter Simulation Environment, Create New Material Stream: Air with the above conditions.
Figure 16-2
5- Parameter Page: Select High Efficiency, Leith/Litch, Particle Efficiency = 95% as shown in figure
16-3
Figure 15-3
Note: for high output and less number of cyclones, select High Output radio
7- Rating tab, Constrain Page: specify the Maximum diameter = 3m and Maximum # Cyclones = 40.
290
It can be seen the Results from Rating tab, Sizing Page
Note: The Polluted liquid can be purified by the same way as in above but using Hydrocyclone instead of
Cyclone.
H.W. Return the above example using the Simple Solid Separator
291
14.2 Rotary Vacuum Filter
The Rotary Vacuum Filter assumes that there is 100% removal of the solid from the solvent stream.
This operation determines the retention of solvent in the particle cake, based on the particle diameter and
sphericity of your defined solid(s). The diameter and sphericity determines the capillary space in the cake
and thus the solvent retention. The solids being separated must be previously specified and installed as
components in the stream attached to this operation. Figure 16-4 shows the Rotary vacuum filter
Figure 16-4
Ex: Stream of slurry at 25oC, 120 kPa and flow rate of 4000 kg/h contains 30% by weight of calcium in
water. It is required to remove the calcium using Rotary vacuum filter. What is the area of the filter if its
diameter is 40 cm?
Solution:
1- Open new case with above components and fluid package: PRSV.
3- From Object Palette, select Rotary Vacuum Filter and enter the information as in the figure
16-5:
292
Figure 16-5
10.00
Figure 16-6
293
Parameter Description
Cycle Time The complete time for a cycle (one complete revolution of the cylinder).
The portion of the cycle between time the cake comes out of the liquid to the time
Dewatering
it is scraped, expressed as a percentage of the overall cycle time.
Submergence The percentage of the overall cycle for which the cake is submerged.
Pressure
Pressure drop across the filter.
Drop
5- Rating tab: Enter the radius of the filter = 0.4 m, Figure 16-7.
Figure 16-7
The constrains of the cake can be enter in the Cake page Design tab
Figure 16-8
H.W.: A rotary drum filter with 30 percent submergence is to be used to filter a concentrated aqueous
slurry of CaCO3 (Mw = 100 , Density = 2800 kg/m3) containing 14.71b of solids per cubic foot of water
(236 kg/m3). The pressure drop is to be 20 in. Hg. If the filter cake contains 50 percent moisture (wet
basis), calculate the filter area required to filter 10 gal/min of slurry when the filter cycle time is 5 min.
The temperature is 20°C.
294
15.3 Baghouse Filter
The Baghouse Filter model is based on empirical equations. It contains an internal curve relating
separation efficiency to particle size. Based on your particle diameter, the reported separation efficiency
for your solids is determined from this curve. The solids being separated must be previously specified and
installed as components in the stream attached to this operation.
Figure 16-9
Ex: Stream of Air at 35oC, 180 kPa and flow rate 60 kmol/h contains 5% by mole sulfur. It is required to
remove the sulfur from air using baghouse filter. What is the efficiency of the sulfur removal? What is
filtration time at which the filter being unable to filtering with 2 kPa pressure drop? What is the relation
between filtration time and pressure drop across the dirty filter?
Solution:
1- Open new case and enter the above components and fluid package: Peng Robinson
2- Enter Simulation Environment and Create New Stream: Air In with the above conditions.
295
Figure 16-10
6- To find the relation between the filtration time and pressure drop across the dirty filter:
2) From Object column, select: Baghouse filter, Variable : Dirty bag pressure drop; click Add
3) Again from Object column: Baghouse filter, Variable : Filtration time; click Add, Close the view.
b. Enter the low bound, high bound and step size as in the figure 16-11 (pressure in kPa)
296
Figure 16-11
e. Click Run then click Plots, the results will appear as in the figure 16-12
Figure 16-12
297
Solid Operations in Aspen Plus
Convective Dryer
Open the Aspen Plus , New Case, select Solids from the right list, Solids with Metric Units,
298
From the Properties: Add the following components in the Components ID Cell; change the Type of the
Al2O3 to Solid as follows:
In the Navigation Pane; click on the Method and select Ideal from the Base Method as follows:
For most cases with only solid components, the “IDEAL” method is likely sufficient
In the “Setup” folder, open the “Solids” form. On the Solid Characterization tab find the section labeled
“Moisture components.” Add Water by highlight the water then click on the Arrow (or double click on
water)
On the PSD Mesh tab of the Solids form, either open the default PSD by select the PSD & click on the
Edit OR select “New...” to generate additional meshes
300
There are four PSD mesh types (equidistant, geometric, logarithmic, and user) which use different
methods to specify the size of the mesh intervals
• User: Specify the size limits for each interval directly in the grid provided.
• Equidistant: Specify the lower and upper limits for the entire distribution, and the number of
intervals, then click Create Mesh. The grid will be filled in with equal-sized intervals.
• Geometric: Specify the lower and upper limits for the entire distribution, then click Create Mesh.
The grid will be filled in with intervals sized so that the volume ratio between consecutive sizes is 2
(since diameter is shown, the ratio of diameters will be the cube root of 2). The upper limit will be
enforced strictly; the number of intervals will be adjusted so that the lower limit appears within the
first interval.
• Logarithmic: Specify the lower and upper limits for the entire distribution, and the number of
intervals, then click Create Mesh. The grid will be filled in with intervals sized so that ratio
between consecutive sizes is constant.
From the PSD Mesh type; select User, No. of intervals = 2, Size units = mm, Fill the table as follows:
301
Build the flowsheet: click on the Main Flowsheet tab
From the Model Palette, select Solid; then select Dryer icon from the list of the Dryer
Click on the Material icon and Connect the appropriate streams with the Dryer as follows:
302
Open the form for the wet solid stream (right click and select Input). You must define this stream using
the specifications on both the Mixed and CI Solid tabs because the stream has both fluid and solid
components. On the CI Solid tab, do not forget to fill out the PSD mesh with the weight fractions for solid
particles in each interval. You can do this by either manually entering values into the matrix, or using a
specified distribution function.
If a stream has only solid components, it is important to only fill out the CI Solid tab within this form
and leave the Mixed tab blank. If a stream is only partly solid, you must fill out both.
303
In the CI Solid tab; enter the data as follows:
Click on the Particle Size Distribution, enter the data required as follows:
Define the dry gas stream: right click on the Gas-In Stream and enter the following data:
304
Mole-frac
Specify that this is a convective dryer under “Dryer type” Next, fill out the rest of the Specifications tab
by defining the dryer’s heat loss and geometry. A typical convective dryer uses the following
assumptions:
• The moisture content and temperature are spatially constant in every particle
305
To describe mass transfer, you can specify either the Sherwood number, mass transfer coefficient, (mass
transfer coefficient)*(surface area), or NTU (number of transfer units). To describe heat transfer you can
either specify the Lewis number or heat transfer coefficient. You also have the option to calculate particle
shrinking.
Navigate to the Drying Curve tab in the Dryer form and fill out the required inputs. You have the option
of generating your normalized drying curve model using either a function or user inputted data.
306
Click Run
Results: In the Navigation Pane; under the Block folder: click Stream Results
you can generate a number of different plots to visually analyze your results. To do so, open the form for
the dryer, navigate to the Plot group on the Home tab on the ribbon, and select the type of plot you wish to
create . The available plots for a Dryer include:
• PSD
• Solid Temperature
• Solid Moisture
• Total Evaporation
• Drying Curve
• Evaporation Rate
307
308
Simple Granulation
309
2
- Enter the Simulation Environment: from the Solids group build the flow sheet below
310
- In Setup>Solids>PSD, define the PSD mesh as Equidistant, 200 intervals, 0-2 mm
Click on Create PSD Mesh
Define the RECYCLE stream on the CI Solid tab (because it is a solid stream)
311
Click on Particle Size Distribution: (Wait some minutes)
- Populate using: A distribution function
-Type of distribution function: RRSB (Rosin-Rammler-Sperling-Bennet)
-Dispersion parameter: 4 (The RRSB distribution is characterized by the parameters n and d63)
-D63: 0.9 mm (Particle diameter (≤ 0.9 mm) which is larger than 63% of solids mass OR = mass
fraction of particles in the population with size ≤ 0.9 mm)
Click Calculate
312
Click on Outlet Flash tab
SCREEN block:
From Navigation Pane; Blocks folder; select SCREEN Or from the flowsheet; double click on the
Deck 1
–Over size stream ID: COARSE Deck 2
–Over size stream ID: PRODUCT
–Selection function: Plitt
–Selection function: Plitt
–Cut size: 1 mm
–Cut size: 0.25 mm
–Separation sharpness: 8
–Separation sharpness: 8
–Offset of Fines: 0
–Offset of Fines: 0
–Entrainment specification: Coarse-split
–Entrainment specification: Coarse-split
–Entrainment fraction: 0
–Entrainment fraction: 0
313
Define CRUSHER block – Specifications
–Outlet PSD calculation method: Select Determine outlet PSD from comminution power and a
distribution function
–Distribution function: Rosin Rammlet Sperling Bennet
–Power specification: 2 kW
–Select parameter(s): RRSB dispersion parameter
–RRSB dispersion parameter: 3
Now the data of the flowsheet is complete; it can be Run the Simulation
Click Run and wait some minutes
315
Screening
A mixture of different particle types in one stream is described by different sub streams (particle types)
– SiO2
– AL2O3
– CaCO3
316
- From the Method in the Navigation Pane; Set property method to ideal
- Go to simulation environment
317
- Go to Solids folder and open the sub streams form
SUBST-A
SUBST-B
SUBST-C
type to CISOLID
– Select PSD as
PSD mesh ID
- Go to the Stream Class tab and open the Stream Class form
SUBST-A
SUBST-B
SUBST-C
318
Define Simulation Mesh
319
Double click on Screen1 & enter the following data Or click on screen1 under the Blocks folder
Double click on Screen2 & enter the following data Or click on screen1 under the Blocks folder
321
- Select “SUBST-C” as sub stream name
322
Results of Screen1
Results of Screen2
323
Spray Drying
particle formulation)
324
Enter the Simulation environment
Note: Dryer from Solids group, GasCyc & Filter (FabFI) from Solid Separators group,
Pump & Blower from Pressure change group, Heater from Exchangers group
Click on the Setup icon: Solids icon; change Stream class to MIXCIPSD
325
Click on the PSD icon;
-Even though no particles enter the process it is necessary to define a particle size mesh, since particles
will be formulated in the spray dryer
–Use of an equidistant mesh with 100 particle size classes from 0 to 250 microns
Streams:
- Feed Definitions
Solution (stream feed) contains 20 wt.-% Ammonium chloride and 80 wt.-% water (defined in the
MIXED sub stream)
326
Drying agent (stream air) is air with ~15 g/kg dry water
- Pump Setup
327
Heater Setup
Specification tab
328
PSD tab
329
Atomization tab
330
- Fabric Filter Setup
331
Efficiency tab
332
- Blower Setup
333
Penicillin Extraction
Open new case, select Blank Simulation , Create
Add the components in Components Name column; rename each in the Component ID as follows
Click on Method folder in the Navigation Pane, from Base Method select NRTL, click NEXT
334
Enter the Simulation Environment and construct the flowsheet as follows:
335
- Enter the data in the Solvent1 as follows
336
Key Components tab:
Streams tab:
337
Pressure tab:
Reactions tab:
338
Click Close, Click RUN
339
Chlorobenzene production
Given: Feed:
100% Cl2
340
Component: Benzene, Chlorine, Chlorobenzene, P-dichlorobenzene, Sodium hydroxide, Sodium chloride,
Water, Hydrogen chloride.
DSTWU2
341
By the same way, enter the conditions of the Cl-Feed:
342
The heaters (HET-Bz & Het-Cl) have the same conditions:
- CSTR: Enter the data as follows: T = 328 K, P = 2.4 bar, Residence time = 0.5 hr
Valid phases: Vapor-liquid
Specification type: Residence time
343
Stream tab: Fill as follows:
Reaction tab: from the Navigation pane click on the Reaction Folder then click new
344
From the new Windows select type: PowerLaw (Kinetic type)
Click New..
345
Fill the windows of reaction as follows
Click Close
Then click New.. again and Fill the data of the 2nd reaction as follows:
Click Close
Kinetics tab:
346
For the 2nd Reaction: Select 2nd reaction from the Menu
Now Add the R-1 to the Reactor as follows: From Navigation Pane; select CSTR
347
Washer Reactor (Stoichiometric Reactor or RStoich)
Reaction tab: Click New.. and Enter the stoichiometric of reaction and Fractional conversion as follows:
Fractional conversion = 1 of component = HCl
348
Short-Cut distillation (DSTWU): Enter the data required as follows
Light key: Bz = 0.99 (Recovery from the feed and not its fraction in the distillate)
Heavy key: Mono = 0.02 (Recovery from the feed and not its fraction in the distillate)
349
For the 2nd Distillation (RadFrac): Enter the data as follows
Configuration tab:
Stream tab
350
Pressure tab:
351
p
Design Specification: Note the flow rate of the mono = 18.4 kmol/hr
We want to be 25 kmol/hr
- From Navigation Pane click on the Flowsheeting options → click on the Design Specs
→ click on the New button
352
Fill the Define tab as follows
353
Vary tab: Fill the data required as follows
What is the mole flow of the Bz required? "see the results in your simulation"
354
References:
1. Mohd. Kamaruddin Abd Hamid, HYSYS®: An Introduction to Chemical Engineering Simulation For
UTM Degree++ Program, 2007
2. Trupti Ambar, Tyagee Chavan , Manali Kavale, S M Walke, Simulation of Process Equipment by
using Hysys, International Journal of Engineering Research and Applications (IJERA) ISSN: 2248-9622
National Conference on Emerging Trends in Engineering & Technology (VNCET-30 Mar’12).
4. Smith J. M., Van Ness H. C., and Abbott M. M., Introduction to Chemical Engineering,
Thermodynamics, Fourth Edition, 1987.
7. Advanced Process Modeling Using Aspen HYSYS, AspenTech Customer Education Training Manual
Work book, Course number EA1000.06.09
9. HYSYS® 2004.2, Tutorials and Applications, Copyright © 1981-2005 by Aspen Technology, Inc.
10. AEA Technology 2000, Process Modeling Using HYSYS with Chemical Industry Focus
11. Jump Start: Acid Gas Cleaning in Aspen HYSYS®, Luisa Herrmann,& Manya Garg, © 2013 Aspen
Technology, Inc. AspenTech®.
12. Jump Start: Rigorous Shell & Tube Exchanger Sizing in Aspen HYSYS® V8, Lauren Sittler, & Tom
Ralston, © 2013 Aspen Technology, Inc. AspenTech®
13. Ammonia Synthesis with Aspen HYSYS® V8.0, Part3 Process Economic Analysis, Copyright © 2012
by Aspen Technolgy.
14. Jump Start: Assay Management in Aspen HYSYS® Petroleum Refining, Luisa Herrmann, Vince Ye,
© 2013 Aspen Technology, Inc. AspenTech®.
15. Process Analysis and Simulation in Chemical Engineering, Iva´n Darı´o Gil Chaves, © Springer 2016.
355
Simulation and Design of Chemical Process Using
Aspen HYSYS
Khalid W. Hameed
356