Thermo
Thermo
Thermo
PRO/II 9.3
Thermodynamic Data
Keyword Manual
May 2014
All rights reserved. No part of this documentation shall be reproduced, stored in a retrieval system, or
transmitted by any means, electronic, mechanical, photocopying, recording, or otherwise, without the
prior written permission of Invensys Systems, Inc. No copyright or patent liability is assumed with
respect to the use of the information contained herein. Although every precaution has been taken in the
preparation of this documentation, the publisher and the author assume no responsibility for errors or
omissions. Neither is any liability assumed for damages resulting from the use of the information
contained herein.
The information in this documentation is subject to change without notice and does not represent a
commitment on the part of Invensys Systems, Inc. The software described in this documentation is
furnished under a license or nondisclosure agreement. This software may be used or copied only in
accordance with the terms of these agreements.
For comments or suggestions about the product documentation, send an e-mail message to
ProductDocumentationComments@invensys.com.
All terms mentioned in this documentation that are known to be trademarks or service marks have been
appropriately capitalized. Invensys Systems, Inc. cannot attest to the accuracy of this
information. Use of a term in this documentation should not be regarded as affecting the validity of any
trademark or service mark. Invensys, Invensys logo, PRO/II, and SimSci are trademarks of Invensys
plc, its subsidiaries and affiliates.
Table of Contents
Table of Contents
Chapter 1
Thermodynamic Data Overview
General Information. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-1
Keyword Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-3
Minimum Required User Input . . . . . . . . . . . . . . . . . . . . . . . . . .1-14
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-14
Order of Input . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-16
Notes Statement (optional). . . . . . . . . . . . . . . . . . . . . . . . . .1-17
Thermodynamic Method Sets . . . . . . . . . . . . . . . . . . . . . . . . . . .1-18
Available Thermodynamic Property Generators. . . . . . . . . .1-20
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-23
Predefined Systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-26
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-26
Input Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-27
Multiple Thermodynamic Sets. . . . . . . . . . . . . . . . . . . . . . . . . . .1-30
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-30
The Default Method Set . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-30
Method Set Requirements . . . . . . . . . . . . . . . . . . . . . . . . . . .1-31
Avoiding Inconsistent Methods. . . . . . . . . . . . . . . . . . . . . . .1-31
Multiple Method Sets in Distillation Columns . . . . . . . . . . .1-31
Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-32
Free-water Decant Considerations. . . . . . . . . . . . . . . . . . . . . . . .1-33
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-33
User-Added Water Enthalpy Methods . . . . . . . . . . . . . . . . .1-34
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-35
Vapor-liquid-liquid Equilibrium Considerations . . . . . . . . . . . . .1-38
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-38
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1-43
Chapter 2
Application Guidelines
General Information. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .2-1
Refinery and Gas Processing Applications . . . . . . . . . . . . . . . . . .2-2
Chapter 3
Generalized Correlations
Ideal and Library Methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-1
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-1
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-1
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-2
Grayson-Streed . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-5
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-5
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-5
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-5
Chao-Seader. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-8
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-8
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-8
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-9
Modifications to GS and CS . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-11
ToC-2
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-11
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-11
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-12
Curl-Pitzer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-15
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-15
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-15
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-15
Braun K10 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-18
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-18
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-18
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-18
Johnson-Grayson . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-21
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-21
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-21
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-21
Lee-Kesler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-23
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-23
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-23
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-24
API Liquid Density . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-26
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-26
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-26
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-26
Rackett Liquid Density . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-28
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-28
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-28
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-28
Costald Liquid Density . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-30
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-30
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-30
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-30
Peneloux Density . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-32
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-32
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-32
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-32
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-33
User-supplied K-value Data. . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-34
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3-34
Chapter 4
Equations of State
Soave Modified Redlich-Kwong . . . . . . . . . . . . . . . . . . . . . . . . . . 4-1
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-1
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-1
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-2
Peng-Robinson. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-7
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-7
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-7
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-8
Modified Soave-Redlich-Kwong and Peng-Robinson . . . . . . . . 4-13
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-13
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-13
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-14
UNIWAALS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-22
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-22
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-22
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-23
Filling in Missing Parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-27
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-27
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-27
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-27
Cubic Equation Of State Alpha Formulations . . . . . . . . . . . . . . . 4-28
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-28
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-28
K-value Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-29
Benedict-Webb-Rubin-Starling . . . . . . . . . . . . . . . . . . . . . . . . . . 4-34
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-34
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-34
Associating Hexamer Equation Of State . . . . . . . . . . . . . . . . . . . 4-39
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-39
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-39
Lee-Kesler-Plocker . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-43
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-43
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-43
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-44
ToC-4
Predictive Peng-Robinson 78. . . . . . . . . . . . . . . . . . . . . . . . . . . .4-47
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .4-47
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .4-47
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .4-48
Predictive Soave-Redlich-Kwong (PSRK) . . . . . . . . . . . . . . . . .4-53
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .4-53
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .4-53
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .4-54
Chapter 5
Special Packages
Alcohols . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-1
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-1
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-1
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-2
Glycols . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-14
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-14
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-14
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-15
Sour Water . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-19
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-19
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-19
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-20
GPA Sour Water. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-22
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-22
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-23
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-23
Amines . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-26
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-26
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-27
User-added Subroutines . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-32
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-32
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-32
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-32
CAPE-OPEN . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-34
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-34
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5-34
ToC-6
Filling in Missing Parameters . . . . . . . . . . . . . . . . . . . . . . . . . . .6-52
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-52
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-53
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-53
Henrys Law for Non-condensible Components . . . . . . . . . . . . .6-57
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-57
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-57
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-58
Henrys Law Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-62
Hayden-OConnell Vapor Fugacity . . . . . . . . . . . . . . . . . . . . . . .6-64
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-64
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-64
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-65
Truncated Virial Vapor Fugacity . . . . . . . . . . . . . . . . . . . . . . . . .6-68
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-68
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-68
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-69
IDIMER Vapor Fugacity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-71
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-71
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-71
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-72
Redlich-Kister, Gamma Heat of Mixing . . . . . . . . . . . . . . . . . . .6-75
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-75
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-75
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .6-75
Chapter 7
Solid Solubility Methods
Van't Hoff Solubility . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7-1
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7-1
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7-1
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7-2
User-supplied Solubility . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7-4
Typical Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7-4
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7-4
Input Description. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7-5
Chapter 9
Method-specific Pure Component Data
General Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9-2
Input Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9-2
Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9-11
Index
ToC-8
Chapter 1
Thermodynamic Data Overview
This chapter provides an overview of the Thermodynamic Data
Category. Detailed documentation, along with examples of
common usage of all the thermodynamic features, follows in
Chapters 2 through 9.
General Information
The program offers many thermodynamic methods and options
ranging from the simple to the very complex. However, the thermo-
dynamic data input requirements for most simulations are simple,
requiring only a few input statements. For example, to use the
Soave-modification to the Redlich-Kwong equation of state in a
flowsheet requires only two lines of code in the keyword file:
Example: THERMODYAMIC DATA
METHOD SYSTEM=SRK
and/or
ENTHALPY(L) {BANK = SIMSCI or NONE} or { BANK=PROII_9.3:SIMSCI,
LIBRARY1:USER1, LIBRARY1:USER2, DIPPR, bankid...},
{ALPHA = ACENTRIC or SIMSCI or bankid},
{COMP = SIMSCI, LIBRARY1:USER1, LIBRARY1:bankid...},
{PURE = LIBRARY1:USER1, and/or LIBRARY1:USER2, and/
or LIBRARY1:bankid...}
HMIX = NONE or GAMMA or RK1 or RK2
<optional data statements> ...
Entropy Options, Optional Chapters 4 & 5
ENTROPY(VL) {BANK = SIMSCI or NONE } or { BANK=PROII_9.3:SIMSCI,
LIBRARY1:USER1, LIBRARY1:USER2, DIPPR, bankid...},
{ALPHA = ACENTRIC or SIMSCI or bankid},
{COMP = SIMSCI, LIBRARY1:USER1, LIBRARY1:bankid...},
{PURE = LIBRARY1:USER1, and/or LIBRARY1:USER2, and/
or LIBRARY1:bankid...}
and/or
ENTROPY(V) {BANK = SIMSCI or NONE} or { BANK=PROII_9.3:SIMSCI,
LIBRARY1:USER1, LIBRARY1:USER2, DIPPR, bankid...},
{ALPHA = ACENTRIC or SIMSCI or bankid},
{COMP = SIMSCI, LIBRARY1:USER1, LIBRARY1:bankid...},
{PURE = LIBRARY1:USER1, and/or LIBRARY1:USER2, and/
or LIBRARY1:bankid...}
and/or
( ) = Keyword qualifiers; { } = Optional entries; Values given are defaults;
Underlined keywords are defaults.
NOTES This statement must be placed along with the METHOD statement
of the specific thermodynamic set to which the information
pertains. One NOTES statement is allowed per thermodynamic set.
TEXT This entry specifies the actual notes line(s). The Notes Data can be
specified in multiple lines. Only the first 4096 characters Ar stored.
The rest will be discarded.
Predefined Systems
General Information
The use of predefined systems is a convenient, shortcut method of
specifying thermodynamic sets. A predefined system may be
selected by using the SYSTEM keyword. This option uses a single
entry to declare several property methods at once.
In general, each predefined system includes methods for computing
K-values, vapor and liquid enthalpies, and vapor and liquid densi-
ties. Predefined systems using generalized and equation of state
methods and special package systems also include an entropy prop-
erty method. Predefined liquid activity systems do not include an
entropy method.
As an example, the following METHOD statement uses the SRK
system of thermodynamic generators (see Soave Modified
Redlich-Kwong on page 4-1, for keyword definitions):
METHOD SYSTEM=SRK
From Table 1-6, we can see that this is equivalent to the following
METHOD statement:
METHOD KVALUE=SRK, ENTHALPY=SRK, ENTROPY=SRK, &
DENSITY(V)=SRK, DENSITY(L)=API
or even more explicitly:
METHOD KVALUE=SRK, &
ENTHALPY(V)=SRK, ENTHALPY(L)=SRK, &
ENTROPY(V)=SRK, ENTROPY(L)=SRK, &
Input Overview
Table 1-6 shows the predefined systems available and the property
methods used for each one.
PR VLE1,3 PR PR PR PR PR API
Examples
1.2: SRK is used to model the phase behavior of a
hydrocarbon/ water mixture.
TITLE PROJECT=MANUAL, USER=SIMSCI, PROB=HC-WATER
COMPONENT DATA
LIBID 1,H2O/ 2,C3/ 3,NC4/ 4,NC5/ 5,NC6
THERMODYNAMIC DATA
METHOD SYSTEM=SRK
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=20/20/20/20/20
UNIT OPERATION
FLASH UID=F1
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
SRK
PR
binary
SRKH
interactions
PRH and pure
component
SRKM
alpha values
PRM
SRKP
PRP
SRKK
SRKS
PPR78
BWRS binary interactions
RK1 / RK2 interaction coefficients
COSTALD component critical properties
RACKETT component Rackett properties
HENRY Henrys Law solubility data
binary coefficients for specific
ALCOHOL NRTL
components
binary interactions and pure
GLYCOL SRKM
component alphas
THERMODYNAMIC DATA
METHOD KVALUE(VLLE)=NRTL, ENTHALPY=LIBRARY
THERMODYNAMIC DATA
METHOD KVALUE(VLE)=NRTL, KVALUE(LLE)=NRTL, &
ENTHALPY=LIBRARY
The user may specify any number of the available databanks by
using the BANK option on each property statement or the AZEO-
TROPE option for azeotropic data. BANK=SIMSCI is the default
specification for all properties except for GLYCOL and ALCOHOL
package K-values. The program searches the databanks in the order
they appear on the PROPERTY statement, and automatically
retrieves any data required for the property method (specified on the
METHOD statement).
Any data supplied in the Component Data Category of input or on
data statements in the Thermodynamic Data Category of input are
used instead of the data retrieved from the databanks. The
BANK=NONE entry suppresses all databank operations.
The following thermodynamic data section,
METHOD KVALUE=NRTL, ...
KVALUE BANK=SIMSCI, ALCOHOL
searches the SIMSCI databank, then the ALCOHOL databank for
VLE interaction parameters.
For more information on the use of the FILL option, see Filling in
Missing Parameters on page 4-27 for equations of state and Fill-
ing in Missing Parameters on page 6-52 for liquid activity meth-
ods.
Input Description
The METHOD Statement
METHOD SYSTEM(VLLE)=method, {L1KEY=i, L2KEY=j,}...
Examples
E-1.3: A mixture of furfural, cyclohexane, heptane and nonane at
85 F and 1.5 atmospheres is to be modeled using NRTL
thermodynamics. This system forms two liquid phases. To
avoid arbitrary allocation of key components by the
calculation algorithm, the key components are designated
by the user. The explicit declaration of the dominant
components ensures that the Liquid2 stream contains the
General Information
Usually, there are several thermodynamic methods suitable for any
given application. The user always should try to determine which
methods give the best representation of the whole flowsheet, while
trying to select the simplest, most appropriate thermodynamic
options. The user should bear in mind that the best thermodynamic
method is the one that gives the best agreement with reality. When
laboratory or actual operating data are available, it may be neces-
sary to try several options and compare the results to obtain the best
possible model.
The following guidelines are divided into four basic categories of
applications. These are:
Refinery and Gas Processing
Natural Gas
Petrochemical
Chemical and Environmental
For each application, the various types of unit operations encoun-
tered have recommended thermodynamic methods.
Recommended
Methods Comments
Recommend
ed Methods Comments
GS / GSE / IGS Usually faster than SRK or PR, but less accurate in the
presence of a high concentration of light components.
Substituting LK for CP enthalpies may improve the answers.
SRK / PR Use SRK or PR if light crudes dominate the top of the tower.
Recommended
Methods Comments
Recommended
Methods Comments
SRK / PR / BWRS These methods give good results for most hydrocarbon and
hydrocarbon-water systems.
SRKKD Use SRKKD(VLLE) for high pressure systems containing
hydrocarbons and water. SRKKD uses more CPU time than
SRK. Do not use SRKKD if other polar components, such
as methanol, are present.
SRKM / PRM / Use these methods for systems containing water and other
SRKS polar components, such as methanol. Always use the VLLE
option with these methods for this type of system.
SRKP / PRP Simpler versions of SRKM and PRM. These methods are
not as good as SRKM or PRM and do not significantly
reduce CPU time.
Recommended
Methods Comments
GPSWATER
This method is based on the GPSWAT program for calculating sour
water equilibrium. The GPSWAT method is used for generating K-
values for sour water components (H2O, H2S, CO2, NH3, CO,
CS2, COS, MeSH and EtSH). All other components are calculated
using SRKM. The recommended range of application is:
Temperature (F) 68 < T < 600 F
Pressure P < 2000 psia
Composition xwtNH3 < 0.40
Sour gas partial pressure PCO2+ PH2S < 1200 psia
Amine Systems
The AMINE package is used for generating K-values for aqueous
amine systems and sour gases including H2S, CO2, H2O, MEA,
DEA, DGA, DIPA, and MDEA. All other component K-values are
calculated using SRKM. Water and exactly one amine must be pres-
ent when using the package.
Data for the equilibrium constants are provided for MEA, DEA,
DGA, DIPA, and MDEA. However, the DIPA data are not recom-
mended for use in final designs. For MDEA, the model is modified
to include composition effects.
Petrochemical Applications
Light Hydrocarbon Systems
SRK, PR, or BWRS is recommended for most light-hydrocarbon,
petrochemical applications.
When solubility of water in hydrocarbon and hydrocarbon in water
becomes important, as it would at high pressure, SRKKD(VLLE) is
recommended.
Use the COSTALD method for liquid density. This option is not the
default and must be requested specifically.
Aromatic Systems
For systems containing all aromatic components, use the ideal
method at low pressures as these systems are very close to ideal.
For systems at pressures above 2 atmospheres, use GS, SRK, or PR
for a more accurate result.
The liquid density can be calculated using the ideal, the default API,
or the COSTALD methods. The COSTALD method has data for
many aromatic components, will give good results at higher tem-
peratures, and is better than API if any light components such as
methane are present. For systems at high temperatures with all aro-
matic components, the API method is as good and uses less CPU
time. The ideal method is best at lower temperatures and should not
be used if the temperature is significantly higher than the boiling
point of any one component.
Recommended
Methods Comments
Aromatic/Non-aromatic Systems
Traditionally, these systems are difficult to model accurately.How-
ever, new equation-of-state mixing rules and alpha formulations can
give excellent results when appropriate interaction parameters and
alpha parameters are available. The SIMSCI databank includes a
large number of interaction parameters, but the user should verify
the availability of all critical parameters to ensure getting good
Recommended
Methods Comments
Non-hydrocarbon Systems
These systems typically contain oxygen, halogen, or nitrogen deriv-
atives of hydrocarbon components and tend to be highly non-ideal.
For low pressure systems, use a liquid activity coefficient method.
For single phase systems, the WILSON, NRTL, and UNIQUAC
methods are equally good. The Wilson method is the simplest and
requires the least amount of computer time. Simpler methods, such
Recommended
Methods Comments
HF Systems
A special equation of state, HEXAMER, has been created for sys-
tems containing molecules that hexamerize in the vapor phase. This
method is recommended for HF systems such as the HF alkylation
Recommended
Methods Comments
Recommended
Methods Comments
Environmental Applications
These systems usually involve stripping dilute pollutants out of
water. The normal liquid activity methods, such as NRTL, do not
usually model these very dilute systems with sufficient accuracy. A
better approach is to use Henry's Law data for the components in
water. The Henry's Law data can be obtained from data sources,
such as the U.S. Environmental Protection Agency( EPA), or are
often available in the SIMSCI databank. In order to model the sys-
tems accurately, the user should supply temperature-dependent val-
Solid Applications
Solid-liquid equilibria for most systems can be represented by the
van't Hoff (ideal) solubility method or by using user-supplied solu-
bility data. In general, for those systems where the solute and sol-
vent components are chemically similar and form a near ideal
solution, the van't Hoff method is appropriate. For non-ideal sys-
tems, solubility data should be supplied. For most organic crystalli-
zation systems, which are very near ideal in behavior, the van't Hoff
SLE method provides good results. The VLE behavior can usually
be adequately represented by IDEAL or VANLAAR methods.
Examples
E-2.1: This example models a hydrocarbon-water mixture that
contains less than 10% of water. The SRK method is
chosen.
TITLE PROB=HC-WATER
PRINT INPUT=ALL
COMP DATA
LlBID 1,H2O/2,C2/3,NC4/4,NC5/5,NC6
THERMO DATA
METHOD SYSTEM=SRK
WATER DECANT=ON
STREAM DATA
PROP STREAM=1,TEMP=100,PRES=100, &
COMP=40/100/100/100/100
UNIT OPERATIONS
FLASH
FEED 1
PROD V=2,L=3,W=4
ISOT TEMP=100,PRES=100
END
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=IDEAL
STREAM DATA
. . .
General Information
The IDEAL method calculates K-values, vapor and liquid densities,
and vapor and liquid enthalpies. LIBRARY methods to calculate
liquid entropies were discontinued starting with PRO/II software
version 4.1. These methods should be used with pure component
streams and streams with very similar components; for instance:
o-xylene / m-xylene / p-xylene / ethylbenzenestyrene.
Input Description
The THERMODYNAMIC DATA statement and the METHOD
statement are discussed in full in Chapter 1, Thermodynamic Data
Overview. The keywords relevant to the IDEAL method are dis-
cussed here.
Pi = i Pisat
are used to calculate K-values.
ENTHALPY Selects the method for enthalpy calculations. IDEAL
(VL or methods may be selected for liquid and/or vapor phase
V and/or L) calculations.
Examples
E-3.1: Using the IDEAL methods, model a 50/50 mix of propane
and normal butane at 50 psia and 100 F. By choosing the
IDEAL system, Curl-Pitzer enthalpies are invoked.
TITLE PROB=IDEAL, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,IC4/ 2,NC4
THERMODYNAMIC DATA
METHOD SYSTEM=IDEAL
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
FLASH UID=FL1
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
E-3.3: For the same problem, use SRK for K-values and
enthalpies, and IDEAL densities.
THERMO DATA
METHOD KVALUE=SRK, ENTHALPY=SRK, DENSITY=IDEAL
or
THERMO DATA
METHOD SYSTEM=SRK, DENSITY=IDEAL
E-3.4: For the same problem, use SRK for everything except
liquid enthalpies. Use the IDEAL method instead for
liquid enthalpy calculations.
THERMODYNAMIC DATA
METHOD SYSTEM=SRK, ENTHALPY(L)=IDEAL
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=GS
STREAM DATA
. . .
General Information
The Grayson-Streed method calculates K-values. It is generally use-
ful for light to mid-range refinery hydrocarbons. It is also useful for
vacuum column and coker simulations when API87 is used. Refer
to the PRO/II Reference Manual for additional limitations.
Table 3-2: Attributes of Grayson-Streed Methods
Properties predicted by GS methods
K-values
Required pure component properties1
Molecular weight Liquid molar volume
Critical temperature Acentric factor
Critical pressure Solubility parameter
Input Description
The METHOD Statement
METHOD SYSTEM=GS, ...
or
METHOD KVALUE=GS, .
Examples
E-3.5: Using the GS methods, model a 50/50 mix of propane and
normal butane at 50 psia and 100 F. By choosing the GS
system, Curl-Pitzer enthalpies are invoked.
TITLE PROB=GS, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMO DATA
METHOD SYSTEM=GS
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=CS
STREAM DATA
. . .
General Information
The Chao-Seader method calculates K-values. It is generally useful
for light to mid-range refinery hydrocarbons. SimSci recommends
that the Grayson-Streed method (GS) be used in preference to CS.
Refer to the PRO/II Reference Manual for additional limitations.
Examples
E-3.8: Using the CS methods, model a 50/50 mix of propane and
normal butane at 100 psia and 100 F. By choosing the CS
system, Curl-Pitzer enthalpies are invoked.
TITLE PROB=CS, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
UNIT OPERATION
FLASH UID=F1
FEED 1
PROD V=2, L=3, W=4
TPSPEC TEMP=100
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=IGS
STREAM DATA
. . .
General Information
The Erbar modification to the Chao-Seader and Grayson-Streed
methods offers improved liquid fugacity coefficients for N2, H2S
and CO2. These methods do not support rigorous three-phase calcu-
lations, but they do support free water decant. The Improved Gray-
son-Streed method has better liquid fugacity coefficients for N2,
H2S, CO, CO2, H2O, and O2. This method supports a rigorous
three phase calculation but can also be used with the free water
decant option turned on.
These methods are generally useful for petroleum refinery applica-
tions. Refer to the PRO/II Reference Manual for additional limita-
tions.
Input Description
The METHOD Statement
Improved Grayson-Streed
METHOD SYSTEM(VLE or VLLE)=IGS, ...
or
METHOD KVALUE(VLE and/or LLE or VLLE)=IGS, ...
Examples
E-3.11: Using the IGS methods, model a 50/50 mix of propane
and normal butane at 100 psia and 100 F. By choosing the
IGS system, Curl-Pitzer enthalpies are invoked.
TITLE PROB=IGS, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMO DATA
METHOD SYSTEM=IGS
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=1000, COMP=50/50
UNIT OPERATION
FLASH UID=F1
FEED 1
PROD V=2, L=3
TPSPEC TEMP=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=SRKM, ENTHALPY=CP
STREAM DATA
. . .
General Information
The Curl-Pitzer method calculates enthalpies and entropies. It is
generally useful for refinery hydrocarbons Refer to the PRO/II Ref-
erence Manual for additional limitations.
Input Description
The METHOD Statement
METHOD SYSTEM=GS or CS or IGS or GSE or CSE or BK10, ...
or
METHOD ..., ENTHALPY(VL)=CP,
or
ENTHALPY(V)= CP and/or ENTHALPY(L)=CP,
ENTROPY(VL)=CP, ...
or
ENTROPY(V)=CP and/or ENTROPY(L)=CP, ...
Examples
E-3.14: Using the CS system, model a 50/50 mix of propane and
normal butane at 100 psia and 100 F. By choosing the CS
system, Curl-Pitzer enthalpies and entropies are
invoked.
TITLE PROB=CP, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMO DATA
METHOD SYSTEM=CS
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=BK10
STREAM DATA
. . .
General Information
The Braun K10 method calculates K-values. It is generally useful
for heavy refinery hydrocarbons at low pressures. Refer to the
PRO/II Reference Manual for additional limitations.
Table 3-6: Attributes of the Braun K10 Methods
Properties predicted by BK10 method
K-values
Required pure component properties1
Molecular weight Acentric factor
Critical temperature Normal boiling point
Critical pressure
Suggested application ranges
Pressure - 0 to 100 psia
Temperature - 100 to 1200 F
Composition - C6 & Heavier
Two liquid phase behavior
Free-water decant - Supported
VLLE - Not supported
1
Automatically supplied for library and petroleum components. Must be
supplied by the user for non-library components.
Input Description
The METHOD Statement
METHOD SYSTEM=BK10, ...
or
Examples
E-3.17: Using the BK10 method, model a 50/50 mix of propane
and normal butane at 100 psia and 100 F. By choosing the
SYSTEM=BK10, Johnson-Grayson enthalpies are
invoked.
TITLE PROB=BK10, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMO DATA
METHOD SYSTEM=BK10
STREAM DATA
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=SRKM, ENTHALPY=JG
STREAM DATA
. . .
General Information
The Johnson-Grayson method calculates enthalpies. It is generally
useful for heavy refinery hydrocarbons. When using the Johnson-
Grayson enthalpy method, it is recommended that the Johnson-
Grayson method be used for both liquid and vapor phases.Refer to
the PRO/II Reference Manual for additional limitations.
Input Description
The METHOD Statement
METHOD SYSTEM=BK10, ...
or
METHOD ..., ENTHALPY(VL)=JG, ...
or
ENTHALPY(V)= JG and/or ENTHALPY(L)=JG
Examples
E-3.20: Using the BK10 system, model a 50/50 mix of propane
and normal butane at 100 psia and 100 F. By choosing the
BK10 system, Johnson-Grayson enthalpies are invoked.
TITLE PROB=JG, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMO DATA
METHOD SYSTEM=BK10
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
FLASH UID=D101
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
Lee-Kesler
Typical Usage
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=SRKM, ENTHALPY=LK
STREAM DATA
. . .
General Information
The Lee-Kesler method calculates enthalpies, entropies and densi-
ties. It is generally useful for refinery hydrocarbons. The liquid den-
sity method is not recommended for hydrocarbons heavier than C8.
Refer to the PRO/II Reference Manual for additional limitations.
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=SRKM, DENSITY(L)=API
STREAM DATA
. . .
General Information
The API method calculates liquid densities. It is generally useful for
refinery hydrocarbons. Refer to the PRO/II Reference Manual for
additional limitations.
Input Description
Table 3-9: Attributes of API Methods
Properties predicted by API methods
Liquid densities
Required pure component properties1
Critical temperature
Critical pressure
Normal boiling point
1Automatically
supplied for library and petroleum components. Must be
supplied by the user for non-library components.
Examples
E-3.24: Using the SRK system, model a 50/50 mix of propane and
normal butane at 100 psia and 100 F. By choosing
SYSTEM= SRK, API liquid densities are automatically
invoked.
TITLE PROB=API, PROJ=THERMO
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMO DATA
METHOD SYSTEM=SRK
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
FLASH UID=D121
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
E-3.25: For the same problem, specify SRK K-values, SRK
enthalpies, SRK entropies, SRK vapor densities, and API
liquid densities explicitly.
THERMO DATA
METHOD KVALUE=SRK, ENTHALPY=SRK, &
ENTROPY=SRK, DENSITY(V)=LK, &
DENSITY(L)=API
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=SRKM, DENSITY(L)=RACKETT
STREAM DATA
. . .
General Information
The RACKETT method calculates liquid densities. It is generally
useful for refinery hydrocarbons as well as non-hydrocarbons.
Refer to the PRO/II Reference Manual for additional limitations.
Input Description
The METHOD Statement
METHOD ..., DENSITY(L)=RACKETT, ...
Examples
E-3.26: Using the SRK system and RACKETT liquid densities,
model a 50/50 mix of propane and normal butane at 100
psia and 100 F. The DENSITY(L)=RACKETT overrides
the default API liquid densities.
TITLE PROB=RACKETT, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMO DATA
METHOD SYSTEM=SRK, DENSITY(L)=RACKETT
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
FLASH UID=FL1
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=SRKM, DENSITY(L)=COSTALD
STREAM DATA
. . .
General Information
The COSTALD method calculates liquid densities. It is generally
useful for aromatics and other light refinery hydrocarbons up to
reduced temperatures of 0.95. Refer to the PRO/II Reference Man-
ual for additional limitations.
Input Description
The METHOD Statement
METHOD . ..., DENSITY(L)=COSTALD, ...
Examples
E-3.28: Using the SRK system, and COSTALD liquid densities,
model a 50/50 mix of propane and normal butane at 100
psia and 100 F. The DENSITY(L)=COSTALD overrides
the default (API) liquid densities.
TITLE PROB=COSTALD, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMO DATA
METHOD SYSTEM=SRK, DENSITY(L)=COSTALD
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
FLASH UID=FLSH
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM=SRKM, DENSITY(L)=SRKPENELOUX or
PRPENELOUX , & DENSITY(V)=SRKPENELOUX or
PRPENELOUX
STREAM DATA
. . .
General Information
The SRK/PR Peneloux methods calculate the liquid and Vapor den-
sities. They are generally useful for aromatics and other refinery
hydrocarbons. Refer to the PRO/II online Reference Manual for
additional information.
Input Description
The METHOD Statement
METHOD . ..., DENSITY(L)=SRKPENELOUX or
PRPENELOUX, ...
Examples
E-3.32: Using the SRK system and SRK Peneloux liquid and
Vapor densities, model a 50/50 mix of propane and
normal butane at 100 psia and 100 F. DENSITY(L)=
SRKPENELOUX overrides the default (API) liquid
densities.
TITLE PROB=PENELOUX, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMO DATA
METHOD SYSTEM=SRK, DENSITY(L)= SRKPENELOUX,
DENSITY(V)= SRKPENELOUX
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
FLASH UID=FLSH
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
E-3.33: Using the SRK system and PR Peneloux Liquid and Vapor
densities, model a 50/50 mix of propane and normal
butane at 100 psia and 100 F.
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMODYNAMIC DATA
METHOD KVALUE=DATA
KVALUE(VLE)
KDATA TABU=50, 100, 150/ 1, 0.4, 0.6, 0.7/ &
2, 0.3, 0.4, 10.8/ 3, 0.5, 0.7, 0.9, PREF=50
or
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMODYNAMIC DATA
METHOD KVALUE=DATA
KVALUE(VLE)
KDATA TABU=50, 100, 150/ 1, 0.4, 0.6, 0.7, PREF=50
STREAM DATA
. . .
Input Description
The METHOD Statement
METHOD KVALUE(VLE)=DATA, ...
or
METHOD KVALUE(VLE)=DATA, KVALUE(LLE)=method, ...
or
METHOD KVALUE(LLE)=DATA, SYSTEM(VLLE)=method, ...
or
METHOD KVALUE(LLE)=DATA, KVALUE(VLE)=method, ...
Examples
E-3.30: Using tabular data, predict the bubble point temperature at
50 psia for an equal molar mixture of propane, n-butane
and n-pentane. The K-value data to be entered are as
follows:
TITLE PROB=KDATA BP
...
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMO DATA
METHOD SYSTEM=SRK
STREAM DATA
. . .
General Information
The Soave modified Redlich-Kwong equation of state predicts K-
values, enthalpies, entropies, and vapor densities. It is most often
used in gas and refining processes. Without significant modifica-
tion, it is generally not useful for highly non-ideal systems. Liquid
densities and VLLE behavior can be predicted with the SRK equa-
tion of state, but are not recommended without using an advanced
form (seeModified Soave-Redlich-Kwong and Peng-Robinson on
page 4-13).
Input Description
The METHOD Statement
METHOD SYSTEM=SRK, ...
or
METHOD KVALUE=SRK, ENTHALPY=SRK,
ENTROPY=SRK, DENSITY(V)=SRK,
DENSITY(L)=API,...
Note: DENSITY(L)=SRK is supported but not recommended.
KVALUE(VLLE)=SRK is supported but not recommended.
Note: The SRK and SAxx statements must follow the KVALUE
statement.
BANK This option selects one or more banks from which to retrieve
vapor and/or liquid phase binary interaction data.
PROCESS Selects the SimSci PROCESS databank.
SIMSCI Selects the SimSci standard databank (default).
NONE This option disables all data retrieval from databanks for
interaction parameters.
Examples
E-4.1: Using the SRK method with default interaction parameters
and using acentric factors for alpha formulation, calculate
the temperature of a 50/50 mix of propane and normal
butane at 50 psia and 50% vaporization.
TITLE PROB=SRK, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
E-4.5: For the previous example, also supply parameters for the
Twu-Bluck-Cunningham alpha formulation for
components 1 and 3. See Cubic Equation Of State Alpha
Formulations on page 4-28.
THERMODYNAMIC DATA
METHOD SYSTEM=SRK
KVALUE SRK 1,2,0.01/ 1,3,0.025
SA06 1,0.75,0.93,1.6/ 3,0.61,0.81,2.1
Peng-Robinson
Typical Usage
...
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMODYNAMIC DATA
METHOD SYSTEM=PR
STREAM DATA
. . .
General Information
The Peng-Robinson (PR) equation of state predicts K-values,
enthalpies, entropies, and vapor densities. It is most often used for
gas and refining processes. Without significant modification, it is
generally not useful for highly non-ideal systems. Liquid densities
Input Description
The METHOD Statement
METHOD SYSTEM=PR, ..
or
METHOD KVALUE=PR, ENTHALPY=PR, ENTROPY=PR,
DENSITY(V)=PR, DENSITY(L)=API,...
BANK This option selects one or more banks from which to retrieve vapor
and/or liquid phase binary interaction data.
PROCESS Selects the SimSci PROCESS databank.
SIMSCI Selects the SimSci standard databank (default).
NONE This option disables all data retrieval from databanks for
interaction parameters.
bankid This option selects the user-created databank named
bankid that is created and maintained by the LIBMGR
program.
FILL Selects the method used for estimating values for hydrocarbon/
hydrocarbon binary interaction data missing from the input file and
any selected databank libraries. This option is not valid for all
modified equations of state. See Filling in Missing Parameters on
page 4-27 for further information on the methods given below.
NONE This option disables estimation of any missing binary
interaction data(default).
GOR This option estimates binary interaction parameters
between methane and ethane and heavier hydrocarbons.
Examples
E-4.7: Using the PR method with default interaction parameters and
using acentric factors for alpha formulation, calculate the
temperature of a 50/50 mix of propane and normal butane at 50 psia
and 50% vaporization.
TITLE PROB=PR, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMODYNAMIC DATA
METHOD SYSTEM=PR
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
FLASH UID=F203
FEED 1
E-4.9: Use the SOUR thermo methods but apply PR to vapor density.
THERMODYNAMIC DATA
METHOD SYSTEM=SOUR, DENSITY(V)=PR
E-4.10: Supply PR kij binaries for components 1-2 and 2-3. Note that
unless binary 1-3 is available in the databank, the kij will have a
value of 0.0. Note that the PR statement must follow the KVALUE
statement even though there are no additional entries on the
KVALUE statement.
THERMODYNAMIC DATA
METHOD SYSTEM=PR
KVALUE
PR 1, 2, 0.01/ 2, 3, 0.025
E-4.11: For the previous example, also supply parameters for the Twu-
Bluck-Cunningham alpha formulation for components 1 and 3. See
Cubic Equation Of State Alpha Formulations on page 4-28.
THERMODYNAMIC DATA
METHOD SYSTEM=PR
KVALUE
PR 1, 2, 0.01/ 2, 3, 0.025
PA06 1, 0.75, 0.93, 1.6/ 3, 0.61, 0.81, 2.1
E-4.12: For the previous example, use binary parameters from the
PROCESS databank instead of the default SIMSCI databank for
entropy calculations. Also, supply parameters for the Twu alpha
formulation for entropy calculations for components 1 and 3.
TITLE PROB=PR, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
...
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMO DATA
METHOD SYSTEM=SRKM
STREAM DATA
. . .
General Information
The modified SRK and PR equations of state predict K-values,
enthalpies, entropies, and vapor densities. They are most often used
in gas and refining processes and are generally useful for non-ideal
systems. VLLE behavior can be predicted with the modified SRK
and PR equations of state. Liquid densities can also be predicted but
are not recommended.
Input Description
The METHOD Statement
Huron-Vidal modifications to PR
METHOD SYSTEM(VLE or VLLE)=PRH, ...
or
METHOD KVALUE(VLE and/or LLE or VLLE)=PRH,
ENTHALPY=PRH, ENTROPY=PRH,
DENSITY(V)=PRH, DENSITY(L)=API, ...
Panagiotopoulos-Reid modifications to PR
METHOD SYSTEM(VLE or VLLE)=PRP, ...
or
METHOD KVALUE(VLE and/or LLE or VLLE)=PRP,
ENTHALPY=PRP, ENTROPY=PRP,
DENSITY(V)=PRP, DENSITY(L)=API, ...
SRKM or PRM(K or R) i, j, kija, kjia, kijb, kjib, kijc, kjic, cij, cji /...
or SRKS(K or R) i, j, kija, kjia, kijb, kjib, kijc, kjic, cij, cji /...
or SRKH or PRH(K or KCAL or KJ) i, j, aij, bij, cij, aji, bji, cji, ij, ij/ ...
or SRKP or PRP(K or R) i, j, kija, kjia, kijb, kjib, kijc, kjic/...
or SRKKD(K or R) i, j, kija, kijb, kijc/ ...
Examples
E-4.13: Using the SRKM method with default interaction
parameters, model a 50/40/10 mix of propane, normal
butane, and water at 50 psia and 50% vaporization.
TITLE PROB=SRKM
COMPONENT DATA
LIBID 1,C3/ 2,NC4/ 3,H2O
THERMODYNAMIC DATA
METHOD SYSTEM(VLLE)=SRKM
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=1,50/ 2,40/ 3,10
UNIT OPERATION
FLASH UID=FL2
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
E-4.15: Use the SOUR thermo methods but apply PRM to vapor
density.
THERMODYNAMIC DATA
METHOD SYSTEM=SOUR, DENSITY(V)=PRM...
E-4.16: Supply SRKM kij binaries for components 1-2 and 2-3.
Note that unless binary 1-3 is available in the databank,
the kij will have a value of 0.0. Note that the SRKM
statement must follow the KVALUE statement even
though there are no additional entries on the KVALUE
statement.
THERMODYNAMIC DATA
METHOD SYSTEM=SRKM
KVALUE
SRKM 1, 2, 0.01, -0.02/ 2, 3, 0.025, 0.04
E-4.17: For the previous example, also supply parameters for the
Twu-Bluck-Cunningham alpha formulation for
components 1 and 3. See Cubic Equation Of State Alpha
Formulations on page 4-28.
THERMODYNAMIC DATA
METHOD SYSTEM=SRKM
KVALUE
SRKM 1,2,0.01,-0.02/2,3,0.025,0.04
SA06 1,0.75,0.93,1.6/3,0.61,0.81,2.1...
UNIWAALS
Typical Usage
...
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMODYNAMIC DATA
METHOD SYSTEM=UNIWAAL
STREAM DATA
. . .
General Information
The UNIWAALS equation of state predicts K-values, enthalpies,
entropies, and vapor and liquid densities. It is most often useful for
highly non-ideal systems if group contribution parameters are sup-
plied either from the built-in databanks or by the user. VLLE behav-
ior can also be predicted with the UNIWAALS equation of state.
Input Description
The METHOD Statement
METHOD SYSTEM(VLE or VLLE)=UNIWAAL, ...
or
METHOD KVALUE(VLE and/or LLE or VLLE)=UNIWAAL,
ENTHALPY=UNIWAAL, ENTROPY=UNIWAAL,
DENSITY(V)=UNIWAAL, DENSITY(L)=API, ...
Note: The UFT1, UNIFAC, and VAxx statements must follow the
KVALUE statement. UOM for UNIFAC is K (default) or KCAL/
KMOL or KJ/KMOL.
BANK This option selects one or more banks from which to retrieve
vapor and/or liquid phase binary group contribution data.
SIMSCI Selects the SimSci standard databank (default).
NONE This option disables all data retrieval from databanks
for interaction parameters.
bankid This option selects the user-created databank named
bankid that is created and maintained by the
LIBMGR program.
ALPHA This option allows access to the databank for alpha formulations.
SIMSCI Selects the SimSci alpha form and values supplied from
the SIMSCI databank (default).
ACENTRIC Uses the original Soave acentric form for alpha.
bankid
COMP Selects order of component libraries and databanks when
searching for pure components required for the current Thermo
system. Properties missing from this databanks will be taken from
global component data search order.
PROCESS Selects the SimSci PROCESS databank.
SIMSCI Selects the SimSci standard databank.
bankid This option selects the user-created databank named
bankid that is created and maintained by the LIBMGR
program.
UFT1 This statement allows entry of the temperature dependent Lyngby
and/or modification of the UNIFAC method (UFT1) and non-
UNIFAC temperature dependent UNIFAC group contribution data for
groups m and k for the UNIWAAL equation of state. UFT1 entries
correspond to the following temperature dependent correlation:
Amk = amk + bmk(T-To) + cmk(T ln{To / T} + T - To)
Examples
E-4.20: Using the UNIWAALS method with default group
contribution data, calculate the temperature of a 50/50 mix
E-4.22: For the previous example, also supply parameters for the
Twu-Bluck-Cunningham alpha formulation for
components 1 and 2, propane and butane. See Cubic
Equation Of State Alpha Formulations on page 4-28.
THERMODYNAMIC DATA
METHOD KVALUE=UNIWAAL, ENTHALPY=UNIWAAL, &
ENTROPY=UNIWAAL, DENSITY(V)=UNIWAAL, &
DENSITY(L)=API
KVALUE
VA06 1,0.34,0.85,2.54/ 2,0.19,0.87,2.96
...
COMPONENT DATA
LIBID 1,C1/ 2,C2/ 3,C3/ 4,NC4/ 5,C5
THERMO DATA
METHOD SYSTEM=SRK
KVALUE FILL=GAO $ or GOR or CPHC
STREAM DATA
. . .
General Information
PRO/II software has an extensive facility to backfill missing binary
interaction data between hydrocarbons for cubic equations of state
(i.e., SRK, PR, and their modifications). Table 4-5 shows the meth-
ods developed by SIMSCI that can be used for estimating binary
interaction parameters.
Input Description
The METHOD Statement
METHOD SYSTEM=SRK or PR or SRKM or PRM or SRKS
K-value Data
KVALUE FILL=NONE or GOR or GAO or CPHC
FILL This selects the method used for estimating values for
hydrocarbon/hydrocarbon binary interaction data missing from
the input file and any selected databank libraries. Missing
parameters are regressed using the option selected to fit the
previously specified K-value cubic equation of state method. A
description of the FILL options shown in Table 4-5.
The options are ordered such that each subsequent method in the
table generates more and more binary interaction parameters, so
they provide a spectrum of results that can be selected based on
the goodness of fit to the specific production data being
modeled.
...
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMO DATA
METHOD SYSTEM=SRKM
KVALUE
SA06 2, 0.19, 0.87, 2.96
STREAM DATA
. . .
General Information
The pure component alpha formulations are used for methods based
on the cubic equation of state. The alpha function controls pure
K-value Data
KVALUE ALPHA=SIMSCI or ACENTRIC or bankid, ...
...,
SA01 or PA01 or VA01 i, c1/ ...
or
SA02 or PA02 or VA02 i, c1, c2, c3/ ...
or
SA03 or PA03 or VA03 i, c1, c2/ ...
or
SA04 or PA04 or VA04 i, c1, c2/ ...
or
SA05 or PA05 or VA05 i, c1, c2/ ...
or
SA06 or PA06 or VA06 i, c1, c2, c3/ ...
or
SA07 or PA07 or VA07 i, c1/ ...
or
SA08 or PA08 or VA08 i, c1, c2, c3/ ...
or
SA09 or PA09 or VA09 i, c1, c2, c3/ ...
or
SA10 or PA10 or VA10 i, c1, c2/ ...
or
SA11 or PA11 or VA11 i, c1, c2/ ...
Note: The SAxx, PAxx and VAxx statements must follow the
KVALUE statement.
Table 4-6 gives the equations for the eleven alpha formulations that
are available.
Table 4-6: Alpha Formulations SA01 to SA22 or PA01 to PA11 or VA01 to VA11
xx = Equation Reference
0.5 2
01 = 1 + C1 1 Tr Soave (1972)
C3 2
= C1 + C2 1 Tr
02 Peng-Robinson (1980)
C2
03 = 1 + ( 1 T r ) C 1 + ------- Soave (1979)
Tr
C2
= exp [ C 1 1 T r
04 Boston-Mathias (1980)
2 ( C2 1 ) 2C 2
exp [ C 1 1 T r
05 Twu (1988)
= Tr
C2 C Twu-Bluck-Cunningham
C3 ( C2 1 ) 3
06 (1990) (Recommended
= Tr exp C 1 1 T r
by SimSci)
C3 C 2
08 = T exp C 1 T Alternative for form (06)
r r
09 Mathias-Copeman
0.5 0.5 2 0.5 3 2
= 1 + C1 1 Tr + C2 1 Tr + C3 1 Tr (1983)
0.5 2
10 = 1 + C1 1 Tr + C 2 ( 1 T r ) ( 0.7 T r ) Mathias (1983)
2 Melhem-Saini-Goodwin
0.5
11 = exp C 1 ( 1 T r ) + C 2 1 T r (1989)
Examples
E-4.24: Using the SRK method with default data, model a 50/50
mix of propane and normal butane at 50 psia and 50%
vaporization. Supply parameters for the Twu-Bluck-
Cunningham alpha formulation SA06 for components 1
and 2, propane and butane.
TITLE PROB=ALPHA
COMPONENT DATA
...
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMO DATA
METHOD SYSTEM=BRWS
STREAM DATA
. . .
General Information
The Benedict-Webb-Rubin-Starling equation of state predicts K-
values, enthalpies, entropies, and vapor and liquid densities. It is
most often useful for light to heavy refinery hydrocarbons. Rigor-
ous two liquid phase behavior is not supported with the BWRS
equation of state.
BANK This option selects one or more banks from which to retrieve vapor
and/or liquid phase binary interaction data.
SIMSCI Selects the SimSci standard databank (default).
NONE This option disables all data retrieval from databanks for
interaction parameters.
bankid This selects the user-created databank named bankid that is
created and maintained by the LIBMGR program.
PURE This option selects one or more banks from which pure
component BWRS parameters are retrieved.
bankid This selects the user-created databank named bankid that is
created and maintained by TDM.
COMP Selects order of component libraries and databanks when
searching for pure components required for the current Thermo
system. Properties missing from this databanks will be taken
from global component data search order.
Examples
4.25: Using the BWRS method with default binary interaction
data, model a 50/50 mix of propane and normal butane at
50 PSIA AND 50% vaporization.
TITLE PROB=BWRS
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMODYNAMIC DATA
METHOD SYSTEM=BWRS
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=50/50
UNIT OPERATION
FLASH UID=DRUM
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1,&
RATE, VALUE=0.5
END
4.27: Use the SOUR thermo methods but apply BWRS to vapor
density.
THERMODYNAMIC DATA
METHOD SYSTEM=SOUR, DENSITY(V)=BWRS...
...
COMPONENT DATA
LIBID 1,IBTE/ 2,IC4/ 3,HF
THERMO DATA
METHOD SYSTEM=HEXAMER
STREAM DATA
. . .
General Information
The HEXAMER equation of state1 predicts K-values, enthalpies,
entropies, and vapor densities. It is most often useful for HF alkyla-
tion and refrigerant synthesis. Rigorous two liquid phase behavior
is supported with the HEXAMER equation of state.
BANK This option selects one or more banks from which to retrieve
vapor and/or liquid phase binary group contribution data.
SIMSCI Selects the SimSci standard databank (default).
NONE This option disables all data retrieval from databanks for
interaction parameters.
bankid This option selects the user-created databank named bankid
that is created and maintained by the LIBMGR program.
...
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMO DATA
METHOD SYSTEM=LKP
STREAM DATA
. . .
General Information
The LKP method is derived from a corresponding states approach
combined with the Benedict-Webb-Rubin-Starling (BWRS) equa-
tion of state (see Benedict-Webb-Rubin-Starling on page 4-34).
The LKP method predicts K-values, enthalpies, entropies, and
vapor and liquid densities. It is most often used for light hydrocar-
bons and for reformer systems containing high quantities of hydro-
gen. VLLE behavior can also be predicted with the LKP method.
Table 4-9: Attributes of the LKP Method
Properties predicted by LKP
K-values Vapor densities
Enthalpies Liquid densities
Entropies
BANK This option selects one or more banks from which to retrieve vapor
and/or liquid phase binary interaction data.
SIMSCI Selects the SimSci standard databank (default).
Examples
E-4.29: Using the LKP method, calculate the temperature of a 50/
50 mix of propane and normal butane at 50 psia and 50%
vaporization.
TITLE PROB=SRK, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMODYNAMIC DATA
METHOD SYSTEM=LKP
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
FLASH UID=FLSH
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
E-4.31: Use the SOUR thermo methods but apply LKP to vapor
density.
THERMODYNAMIC DATA
METHOD SYSTEM=SOUR, DENSITY(V)=LKP
...
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMODYNAMIC DATA
METHOD SYSTEM= PPR78
STREAM DATA
. . .
General Information
The Predictive Peng-Robinson 78 (PPR78) equation of state pre-
dicts K-values, enthalpies, entropies, and vapor densities. It is used
for the prediction of phase equilibria of binary mixtures in both the
subcritical and critical regions very accurately. This method has
been applied to mixtures containing alkanes, aromatics, napthenes,
carbon dioxide, hydrogen sulfide, and nitrogen. The present imple-
mentation is limited to PPR78 model - supported groups only. For
unsupported components, the calculations for kij interaction param-
eters will be defaulted to PR EOS. Liquid densities and VLLE
behavior predictions are not recommended without using an
advanced form of the PR equation of state ( See Peng-Robinson
on page 4-7.).
Input Description
The METHOD Statement
METHOD SYSTEM= PPR78, ...
or
METHOD KVALUE=PPR78, ENTHALPY=PPR78,
ENTROPY=PPR78, DENSITY(V)=PPR78,
DENSITY(L)=API,........
BANK This option selects one or more banks from which to retrieve vapor
and/or liquid phase binary interaction data.
SIMSCI Selects the SimSci standard databank (default).
NONE This option disables all data retrieval from databanks for
interaction parameters.
bankid This option selects the user-created databank named
bankid that is created and maintained by the LIBMGR
program.
FILL Selects the method used for estimating values for hydrocarbon/
hydrocarbon binary interaction data missing from the input file and
any selected databank libraries. This option is not valid for all
modified equations of state. See Filling in Missing Parameters on
page 4-27. for further information on the methods given below.
NONE This option disables estimation of any missing binary
interaction data(default).
GOR This option estimates binary interaction parameters
between methane and ethane and heavier hydrocarbons.
Example
E-4.32: Using the PPR78 method with default interaction
parameters and using acentric factors for alpha
formulation, calculate the temperature of a 50/50 mix of
propane and normal butane at 50 psia and 50%
vaporization.
TITLE PROB=PPR78, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,C3/ 2,NC4
THERMODYNAMIC DATA
METHOD SYSTEM=PPR78
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP=50/50
UNIT OPERATION
FLASH UID=F203
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
The table below shows the 15 PPR78 groups defined for reference:
Table 4-11: PPR78 Method Supported Groups
Group No. Group Name IPS UNIFAC
Group
1 CH3 900
2 CH2 901
3 CH 902
4 C 903
5 CH4 2011
6 C 2 H6 None : Must
Convert Ethane
to PPR78 group
6.
7 CHaro 1200
8 Caro 1222
9 Cfused Aromatic 1201
Rings
10 CH2, Cyclic 1000
11 CH , Cyclic 1001
(for CHcyclic)
1002
(for Ccyclic)
12 CO2 2111
13 N2 2199
...
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMODYNAMIC DATA
METHOD SYSTEM= PSRK
STREAM DATA
. . .
General Information
PSRK is a group contribution equation of state based on a combina-
tion of the Soave-Redlich-Kwong equation of state and a mixing
rule whose parameters are determined by the UNIFAC method.
This model can be used for VLE-predictions over a large tempera-
ture and pressure range and is easily extended to mixtures contain-
ing supercritical compounds.
Input Description
The METHOD Statement
METHOD SYSTEM= PSRK, ...
or
METHOD KVALUE=PSRK, ENTHALPY=PSRK,
ENTROPY=PSRK, DENSITY(V)=PSRK,
DENSITY(L)=API,........
Example
E-4.35: Using the PSRK method, calculate the temperature of a 50/
50 mix of propane and normal butane at 50 psia and 50%
vaporization.
TITLE PROB=PPR78, PROJ=THERMO
DESC THERMO MANUAL PROBLEM
PRINT INPUT=ALL
COMPONENT DATA
E-4.37: Use the SOUR thermo methods but apply PSRK to vapor
density.
THERMODYNAMIC DATA
METHOD SYSTEM=SOUR, DENSITY(V)=PSRK
Alcohols
Typical Usage
...
COMPONENT DATA
LIBID 1,ETOH/ 2, MEOH/ 3, H2O
THERMO DATA
METHOD SYSTEM=ALCOHOL
STREAM DATA
. . .
General Information
The ALCOHOL package is used to predict VLE and/or LLE phase
behavior. This method does not support free water decant. This sys-
tem uses a special set of NRTL binary interaction data for systems
containing alcohols, water, and other polar components.
The ALCOHOL package is generally useful for applications
involving alcohols, especially for azeotropic distillation common in
alcohol dehydration plants. Refer to the PRO/II Reference Manual
for information.
When unit is K:
b ij c ij
ij = a ij + ----- + -----2 + d ij InT + e ij T + f ij TInT
T T
+ f ij RTIn ( RT )
AZEOTROPE These statements are used to enter data for binary pairs. This data
and/or will be regressed to the coefficients for the liquid activity method
INFINITE used to generate the K-values. Multiple statements may be used to
and/or enter data with different units, or the data can be entered on a single
statement.
MUTUAL
and/or This statement sets the binary interaction parameters to zero.
IDEAL
Acids
Acetic Acid C2H4O2 HAC, HDAC
Formic Acid CH2O2 HFOR
Ketones
Methyl Ethyl C4H8O MEK
Ketone C3H6O DMK
Acetone
Esters
Ethyl Acetate C4H8O2 EOQC
Methyl Formate C2H4O2 MFOR
Miscellaneous
Acetaldehyde C2H4O ACH
Sulfolane C4H8O2S SULFLN
Light Gases
Hydrogen H2 H2
Nitrogen N2 N2
Oxygen O2 O2
Carbon Dioxide CO2 CO2
Examples
E-5.1: Using the ALCOHOL system, model a 50/25/20/5 mix of
water, benzene, ethanol and cyclohexane at 50 psia and
20% vaporization. Use SRK for vapor phase fugacity
calculations, and supply NRTL binary interaction data.
TITLE PROB=ALCOHOL
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,H2O/2,BENZENE/3,ETOH/4,CHEN
THERMO DATA
METHOD SYSTEM=ALCOHOL, PHI=SRK
KVALUE FILL=UNIFAC, POYNTING=ON
NRTL 1,2,5.256,219.685,-5.645,288.34,0.2/ &
1,3,1.015,536.264,0.4985,-456.0,0.1448/ &
2,3,-2.748,1472.24,-0.449,440.51,0.5355
STREAM DATA
PROP STREAM=1,TEMP=100,PRES=100, &
COMP=50/25/20/5
...
COMPONENT DATA
LIBID 1,H2O/ 2, CH/ 3, TEG
THERMO DATA
METHOD SYSTEM=GLYCOL
STREAM DATA
. . .
General Information
The GLYCOL package is used to predict VLE and/or LLE phase
behavior. This method does not support free water decant. It uses a
special set of SRKM binary interaction data and alpha parameters
for systems containing glycol, water, and other components.
The GLYCOL package is generally useful for applications involv-
ing triethylene glycol, and to a lesser extent, diethylene glycol, and
ethylene glycol. It is useful especially for TEG dehydration plants.
Refer to the PRO/II Reference Manual for information.
SRKM(K or R)i, j, kija, kjia, kijb, kjib, kijc, kjic, cij, cji / ...
Examples
E-5.3: Using the GLYCOL system, model a 50/25/20/5 mix of
water, benzene, ethanol and TEG at 50 psia and 20%
vaporization. Supply SRKM binary interaction data.
TITLE PROB=GLYCOL
COMPONENT DATA
LIBID 1,H2O/2,BENZENE/3,ETOH/4,TEG
THERMO DATA
METHOD SYSTEM=GLYCOL
KVALUE
SRKM 2,4,0.05,0.0,0.0,0.3,0.07,0.0,0.0,0.3
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=50/25/20/5
UNIT OPERATION
FLASH UID=FLSH
FEED 1
PROD V=2, L=3
TPSPEC TEMP=100
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.2
END
E-5.4: For the same problem, use K-values and IDEAL liquid
densities.
THERMO DATA
METHOD KVALUE=GLYCOL, ENTHALPY=SRKM, &
DENSITY(L)=IDEAL, DENSITY(V)=SRKM
...
COMPONENT DATA
LIBID 1,H2O/ 2, NH3/ 3, CO2
THERMO DATA
METHOD SYSTEM=SOUR
STREAM DATA
. . .
General Information
The SOUR package is used to predict VLE and/or LLE phase
behavior. This method does not support free water decant. This sys-
tem uses the API/EPA SWEQ (Sour Water EQuilibrium) method
developed by Grant Wilson to model sour water components NH3,
H2S, CO2, and H2O. SRKM generates phase equilibria for all other
components.
The SOUR package is generally useful for applications involving
sour water containing less than 30% by weight of sour components.
SRKM(K or R) i, j, kija, kjia, kijb, kjib, kijc, kjic, cij, cji /...
Input Description
The METHOD Statement
METHOD SYSTEM(VLE or VLLE)=GPSWATER,...
or
METHOD KVALUE(VLE and/or LLE or VLLE)=GPSWATER,
ENTHALPY(V)=SRKM, ENTHALPY(L)=IDEAL,
ENTROPY=SRKM, DENSITY(V)=SRKM,
BANK This option selects the databank from which to retrieve vapor
and/or liquid phase binary interaction data.
SIMSCI Selects the SimSci standard databank
(default).
NONE This option disables all data retrieval from
databanks for interaction parameters.
bankid This option selects a user-created databank
named bankid that is created and
maintained with the LIBMGR program.
ALPHA This option allows access to the databank for alpha
formulations.
Examples
E-5.7: Using the GPSWATER system, model a 80/5/10/5 mix of
water, H2S, CO2 and NH3 at 50 psia and 20%
vaporization.
TITLE PROB=GPSWATER
COMPONENT DATA
LIBID 1,H2O/2,H2S/3,CO2/4,NH3
THERMO DATA
METHOD SYSTEM=GPSWATER
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100,&
COMP=80/5/10/5
Amines
Typical Usage
...
COMPONENT DATA
LIBID 1, H20/2, H2S/3, MEA/4, C1
THERMO DATA
METHOD SYSTEM=AMINE
STREAM DATA
...
General Information
The AMINE package is used to predict VLE and/or LLE phase
behavior. This method does not support free water decant. This sys-
tem uses the Kent-Eisenberg model for reaction equilibria with
MEA, DEA, or DIPA and an additional residence time correction
for MDEA or DGA. A correction is applied to IDEAL liquid
enthalpies to account for heats of reaction. The DIPA data are not
recommended for final design purposes.
The AMINE package is generally useful for gas sweetening pro-
cesses with a single amine. Only one amine at a time is allowed.
Table 5-7: Attributes of the AMINE Package
Properties predicted by the AMINE method
K-values
Required pure component properties1
Input Description
The METHOD Statement
METHOD SYSTEM(VLE or VLLE)=AMINE,...
or
METHOD KVALUE(VLE and/or LLE or VLLE)=AMINE,
ENTHALPY(V)=SRKM, ENTHALPY(L)=AMINE,
ENTROPY=SRKM, DENSITY(V)=SRKM,
DENSITY(L)=IDEAL, ...
SRKM(K or R) i, j, kija, kjia, kijb, kjib, kijc, kjic, cij, cji /...
BANK This option selects the databank from which to retrieve vapor
and/or liquid phase binary interaction data.
SIMSCI Selects the SimSci standard databank (default).
NONE This option disables all data retrieval from
databanks for interaction parameters.
bankid This option selects a user-created databank
named bankid that is created and maintained
with the LIBMGR program.
ALPHA This option allows access to the databank for alpha formulations.
SIMSCI Selects the SimSci alpha formulation and values
supplied from the SIMSCI databank (default).
ACENTRIC Uses original Soave acentric form for alpha
(default).
bankid This option selects a user-created databank
named bankid that is created and maintained
with the LIBMGR program.
COMP Selects order of component libraries and
databanks when searching for pure components
required for the current Thermo system.
Properties missing from this databanks will be
taken from global component data search order.
PROCESS Selects the SimSci PROCESS databank
SIMSCI Selects the SimSci standard databank.
Examples
E-5.9: Using the AMINE system, model a mix of water, H2S,
CO2, MEA, and NH3 at 50 psia and 20% vaporization.
TITLE PROB=AMINE
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,H2O/2,H2S/3,CO2/4,NH3/5,MEA
THERMO DATA
METHOD SYSTEM=AMINE
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=30/5/10/5/50
UNIT OPERATION
FLASH UID=FL3
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.2
END
E-5.10: For the same problem, use MDEA and supply a residence
time correction factor of 1.0.
COMPONENT DATA
LIBID 1,H2O/2,H2S/3,CO2/4,NH3/5,MDEA
THERMO DATA
COMPONENT DATA
LIBID 1,IC4/ 2, NC4/ 3, H2O
THERMO DATA
METHODKVALUE=U1, ENTHALPY=U2, DENSITY=U3
STREAM DATA
General Information
User-added subroutines can be supplied to calculate equilibrium K-
values and to generate liquid and/or vapor enthalpy data, liquid and/
or vapor entropy data and liquid and/or vapor density data.
Table 5-8: Attributes of User-Added Subroutines
Properties predicted by User-Added Subroutines
K-values Liquid/vapor entropies
Liquid/vapor densities Liquid/vapor enthalpies
Input Description
The METHOD Statement
METHOD
KVALUE(VLE and/or LLE or VLLE)=U1 or U2 or ... U15,
ENTHALPY(VL or V and/or L)=U1 or U2 or ... U15,
ENTROPY(VL or V and/or L)=U1 or U2 or ... U15,
DENSITY(VL or V and/or L)=U1 or U2 or ... U15
Examples
E-5.11: Using a user-supplied liquid enthalpy method U1 and the
GLYCOL package, model a 50/25/20/5 mix of water,
benzene, ethanol and TEG at 50 psia and 20%
vaporization. Supply data for U1 for components 2 and 4.
TITLE PROB=UASENTH
COMPONENT DATA
LIBID 1,H2O/2,BENZENE/3,ETOH/4,TEG
THERMO DATA
METHOD SYSTEM=GLYCOL, ENTHALPY(L)=U1
ENTHALPY(L)
UDATA 2,2.5/4,3.0
STREAM DATA
PROP STREAM=1,TEMP=100,PRES=100, &
COMP=50/25/20/5
UNIT OPERATION
FLASH UID=FL3D
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, RATE, VALUE=0.2
END
. . .
COMPONENT DATA
LIBID 1,PROPANE/2,BENZENE/3,METHANE, &
BANK=PROCESS,SIMSCI
THERMODYNAMIC DATA
METHOD SYSTEM=CO, PID=SIMSCI.THERMOSYSTEM.1, &
PNAME=COC3C6C1SRK, SET=CO01 . . .
. . .
Input Description
The Method Statement
NRTL
Typical Usage
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=NRTL
STREAM DATA
. . .
General Information
The Non-Random Two Liquid (NRTL) liquid activity method is
used to predict VLE and/or LLE phase behavior. This method does
not support free water decant.
The NRTL method is generally useful for non-ideal applications,
especially for partially immiscible systems. Refer to the PRO/II
Reference Manual for additional information.
Input Description
The METHOD Statement
METHOD SYSTEM(VLE or VLLE)=NRTL, PHI=IDEAL,
{HENRY},...
or
METHOD KVALUE(VLE and/or LLE or VLLE)=NRTL,
PHI=IDEAL, {HENRY},..
j ji G ji x j x j G ij k x k kj G kj
In i = ----------------------- + ------------------- ij --------------------------
k G ki x k k G kj x k k G kj x k
j
When unit is K:
b ij c ij
ij = a ij + ----- + -----2 + d ij InT + e ij T + f ij TInT
T T
When unit is KCAL/KMOL or KJ/KMOL:
b ij c ij
ij = a ij + ------- + -----------
- + d ij In ( RT ) + e ij RT
RT R 2 T 2
+ f ij RTIn ( RT )
G ji = exp ( ji ji ) and ij = ij + ij T
AZEOTROPE These statements are used to enter data for binary pairs.
and/or This data will be regressed to the coefficients for the liquid
INFINITE and/ activity method used to generate the K-values. Multiple
or MUTUAL statements may be used to enter data with different units or
the data can be entered on a single statement.
and/or This statement sets the binary interaction parameters to
IDEAL zero.
Examples
E-6.1: Using the NRTL system, model a 50/40/10 mix of DIPE,
IPA and water at 50 psia and 50% vaporization.
TITLE PROB=NRTL
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,H2O/ 2,BENZENE/ 3,ETOH/ 4,CHEN
THERMO DATA
METHOD SYSTEM=NRTL, PHI=SRK
KVALUE FILL=UNIFAC, POYNTING=ON
NRTL 1,2, 5.256, 219.685, -5.645, 288.34, 0.2 / &
1,3, 1.015, 536.264, 0.4985, -456.0, 0.1448 / &
2,3, -2.748, 1472.24, -0.449, 440.51, 0.5355
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, COMP= 50/ 25/ 20/ 5
UNIT OPERATION
FLASH UID=F204
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.2
END
E-6.2: For the same problem, use NRTL K-values, IDEAL
enthalpies and liquid densities and SRK vapor densities.
THERMO DATA
METHOD KVALUE=NRTL, ENTHALPY=IDEAL, &
DENSITY(L)=IDEAL, DENSITY(V)=SRK
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=UNIQUAC
STREAM DATA
. . .
General Information
The UNIQUAC liquid activity method is used to predict both VLE
and VLLE phase behavior. This method does not support free water
decant.
The UNIQUAC method is generally useful for highly non-ideal
applications, especially for partially immiscible systems. Refer to
the PRO/II Reference Manual for additional information.
C R
ln i = ln i + ln i
M
C Z i
ln i = ln -----i + --- q i ln -----i + l i -----
xi 2 i i Xj lj
j=1
M M
j ij
= q i 1 ln
R
ln i j ji -----------------------
M
j=1 j=1
k kj
k=1
U ij (when unit is K)
ij = exp -------
T
U ij = a ij + b ij T
xi qi Z
i = ------------------
- l j = --- ( r j q j ) ( r j 1 )
M 2
xi qi
i=1
xi ri
i = -----------------
M
-
xj rj
j=1
A wi V wi
q i = ------------------9- r i = ------------- , Z = 10
2.5 10 15.17
AZEOTROPE These statements are used to enter data for binary pairs.
and/or This data will be regressed to the coefficients for the liquid
INFINITE and/ activity method used to generate the K-values. Multiple
or MUTUAL statements may be used to enter data with different units or
the data can be entered on a single statement.
and/or This statement sets the binary interaction parameters to
IDEAL zero.
Examples
E-6.3: Using the UNIQUAC system, model a 50/40/10 mix of
DIPE, IPA and water at 50 psia and 50% vaporization.
TITLE PROB=UNIQUAC
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,H2O/ 2,BENZENE/ 3,ETOH/ 4,CHEN
THERMO DATA
METHOD SYSTEM=UNIQUAC
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
UNIFAC
Typical Usage
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=UNIFAC
STREAM DATA
. . .
General Information
The UNIFAC liquid activity method is used to predict both VLE
and VLLE phase behavior. This method does not support free water
decant.
The UNIFAC group contribution method is generally useful for
non-ideal applications at low pressures. It is generally restricted to
components with 10 or fewer different structural groups and sys-
tems containing low molecular weight polymers. Refer to the
PRO/II Reference Manual for additional information.
Input Description
The METHOD Statement
METHOD SYSTEM(VLE or VLLE)=UNIFAC, PHI=IDEAL,
{HENRY},...
or
METHOD KVALUE(VLE and/or LLE or VLLE)=UNIFAC,
PHI=IDEAL, {HENRY},..
where:
i Z i i
ln i = ln ------i + 1 ------ --- q i ln ------ + 1 ------
C
xi xi 2 i i
xj rj xj qj
j=1 j=1
NK NK
i i
ri = k Rk , qi = k Qk
k=1 k=1
11 --- * . No interactions
20 --- @ * V VLE interactions - ( 634)
21 --- V V . L LLE interactions - ( 208)
27 --- . . . . * Both VLE and LLE - ( 91)
30 --- * V @ V . $ Both VLE and T - ( 48)
40 --- V . V V . V @ VLE - LLE and T - ( 103)
43 --- . . V V . V . T TEMP dependent - ( 1)
46 --- V . V . . V . .
50 --- . . $ V . . . . V
52 --- . . $ V . . . . . V
55 --- * * @ V . . . . V V .
56 --- . . V V . V . . . . . .
59 --- . . $ . . . . . V V V . .
60 --- @ V @ V V $ V . V V $ . V V
70 --- @ * @ V . @ . V V . $ V V V @
73 --- V . V . . V . . V . . . . . V V
281 --- . . . . . . . . V . . . .
V . . . . . . . . . . V . V . . . . . .
287 --- . . V . . . . . . . V . .
V V V . . . . . . . . V . V V . . . . V
290 --- . . V . . . . . . . . . .
. . . . . . . . . . . V . V . . . V . .
300 --- V . V . . . . . . . . . .
V . V V . V V V . V . V . V V V V V V V
320 --- V . V . . V . . . . . . .
. . V V . V . V V . . V . V V V . V V .
330 --- . . . . . . . . . . . . .
. . . . . . . . . . . V . . . . . . . .
340 --- . . . . . . . . . . . . .
. . . . . . . . . . . V . . . . . . . .
350 --- . . . . . . . . . . . . .
. . . . . . . . . . . V . . . . . . . .
360 --- . . . . . . . . . . . . .
. . . . . . . . . . . V . . . . . . . .
385 --- L * * . . . . . V . . * .
. . V * . V . . . L . * . * * * V * L .
388 --- . . V . . . . . . . . . .
V . . V . V . V . V . V . V V V V . . .
| | | | | | | | | | | | |
| | | | | | | | | | | | | | | | | | | |
| | | | | | | | | | | | |
| | | | | | | | | | | | | | | | | | | |
1 1 1 1 1 1 1 1 1
1 1 2 2 2 3 4 4 4 5 5 5 5 5 6 7 7 7 8 8 8 8 8 9 1 1 2 2 3 4 4 4 4
0 1 0 1 7 0 0 3 6 0 2 5 6 9 0 0 3 6 1 4 6 7 8 0 0 4 0 2 0 0 1 3 4
Examples
E-6.5: Using the UNIFAC system, model 50/40/10 mix of DIPE,
IPA, and water at 50 psia and 50%vaporization.
TITLE PROB=UNIFAC
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,DIPE/ 2,IPA/ 3,H2O
THERMO DATA
METHOD SYSTEM=UNIFAC
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=50/40/10
UNIT OPERATION
FLASH UID=F2
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
Modifications to UNIFAC
Typical Usage
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=UFT1
STREAM DATA
. . .
General Information
The modified UNIFAC liquid activity methods are used to predict
both VLE and VLLE phase behavior. These methods do not support
free water decant.
The modified UNIFAC group contribution methods are generally
useful for non-ideal applications at low pressures. The Lyngby
modification UFT1, the Dortmund modification UFT2, and UFT3
modified methods can be used if temperature dependent UNIFAC
data is available. The UNFV free volume data method is designed
to model polymer solutions and should not be used for non-polymer
systems. Refer to the PRO/II Reference Manual for additional infor-
mation.
Input Description
The METHOD Statement
METHOD SYSTEM(VLE or VLLE)=UF1 or UFT2 or UFT3 or
UNFV, PHI=IDEAL, {HENRY},...
or
METHOD KVALUE(VLE and/or LLE or VLLE)=UF1 or UFT2 or
UFT3 or UNFV, PHI=IDEAL,{HENRY},..
and/or
T
A mk = a mk + b mk ( T T o ) + c mk T ln -----o + T T o
T
where: T o = 298.25K
23
c i ri xi
ln i = ln -----i + 1 ----- , i = -------------------------
NOC
xi xi
23
rj xj
j=1
UFT2 This statement allows entry of the temperature-dependent
group interaction data for the Dortmund modification of the
UNIFAC method. The liquid activity coefficients are
calculated from the following equations:
Similar to the UNIFAC method, except
2
A mk = a mk + b mk T + c mk T
i Zq i ln i i
ln i = ln -----i + 1 ----- -------- ----------- + 1 ------
c
xi x i 2 i i
34
ri xi
i = -------------------------
NOC
34
rj xj
j=1
UFT3 This selects temperature-dependent group interaction data for
a modification of the UNIFAC method. The liquid activity
coefficients are computed using the equations following:
Similar to the UNIFAC method, except
3
A mk = a mk + b mk T + c mk T
V i
13 ------
- 1
FV V i 1 m
V
ln i = 3C i ln 1 3 -
------------------- C i ---------------------
V m 1 1 ----------1
-
V 1 3
i
NK
Vi
V i = ---------------------- , 1 i
15.17br' i r' i = ------
Mi k Rk
k=1
Vi wi + Vj wj
V m = -------------------------------------------------- , b = 1.28
15.17b ( r' i w i + r' j w j )
where:
Vi = volume per gram of solvent i
Mi = molecular weight of solvent i
Wi = weight fraction of component i
Ci = number of degrees of freedom per molecule of
solvent i (=3.3)
AZEOTROPE These statements are used to enter data for binary pairs.
and/or INFINITE This data will be regressed to the coefficients for the
and/or MUTUAL liquid activity method used to generate the K-values.
Multiple statements may be used to enter data with
different units or the data can be entered on a single
statement.
and/or This statement sets the binary interaction parameters to
IDEAL zero.
Data input using the UNIFAC, UFT1, UFT2, UFT3, UNFV, AZEO-
TROPE, INFINITE, MUTUAL and IDEAL statements is used in
preference to any data retrieved from the databanks or estimated
using FILL options. See Filling in Missing Parameters on
page 52, for further details.
Examples
E-6.7: Using the UFT1 liquid activity method, model a 50/40/10
mix of DIPE, IPA and water at 50 psia and 50%
vaporization.
TITLE PROB=UFT1
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,DIPE/ 2,IPA/ 3,H2O
THERMO DATA
METHOD SYSTEM=UFT1
STREAM
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=50/50
UNIT OPERATION
FLASH UID=FLSH
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=WILSON
STREAM DATA
. . .
General Information
The WILSON liquid activity method is used to predict VLE phase
behavior. This method does not support free water decant.
The WILSON liquid activity method is generally useful for slightly
non-ideal applications. Refer to the PRO/II Reference Manual for
additional information.
Vj
L
a ij (when unit is K)
A ij = -----L- exp ---------
Vi T
Vj
L
a ij (when unit is KCAL/KMOL
A ij = -----L- exp --------- or KJ/KMOL)
Vi RT
AZEOTROPE These statements are used to enter data for binary pairs.
and/or These data will be regressed to the coefficients for the liquid
INFINITE activity method used to generate the K-values. Multiple
and/or statements may be used to enter data with different units, or
the data can be entered on a single statement.
MUTUAL
and/or IDEAL This statement sets the binary interaction parameters to zero.
Examples
E-6.9: Using the WILSON system, model a 50/40/10 mix of
DIPE, IPA and water at 50 psia and 50% vaporization.
TITLE PROB=WILSON
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,DIPE/ 2,IPA/ 3,H2O
THERMO DATA
METHOD SYSTEM=WILSON
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=50/40/10
UNIT OPERATION
FLASH UID=FLSH
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
Van Laar
Typical Usage
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=VANLAAR
STREAM DATA
. . .
Input Description
The METHOD Statement
METHOD SYSTEM(VLE or VLLE)=VANLAAR,
PHI=IDEAL, {HENRY},...
or
METHOD KVALUE(VLE and/or LLE or VLLE)=VANLAAR,
PHI=IDEAL, {HENRY}, ....
where:
Examples
E-6.11: Using the VANLAAR system, model a 50/40/10 mix of
DIPE, IPA and water at 50 psia and 50% vaporization.
TITLE PROB=VANLAAR
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,DIPE/ 2,IPA/ 3,H2O
THERMO DATA
METHOD SYSTEM=VANLAAR
STREAM DATA
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=MARGULES
STREAM DATA
. . .
General Information
The MARGULES liquid activity method is used to predict VLE and
VLLE phase behavior. This method does not support free water
decant.
The MARGULES liquid activity method is generally useful for
slightly non-ideal applications. Refer to the PRO/II Reference Man-
ual for additional information.
2 2
ln i = ( 1 x i ) [ A i + 2 ( B i A i D i )x i + 3D i x i ]
where:
d ij = d ji
AZEOTROPE These statements are used to enter data for binary pairs.
and/or This data will be regressed to the coefficients for the
INFINITE liquid activity method used to generate the K-values.
and/or Multiple statements may be used to enter data with
different units or the data can be entered on a single
MUTUAL
statement.
and/or This statement sets the binary interaction parameters to
IDEAL zero.
Examples
E-6.13: Using the MARGULES system, model a 50/40/10 mix of
DIPE, IPA and water at 50 psia and 50% vaporization.
TITLE PROB=MARGULES
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,DIPE/ 2,IPA/ 3,H2O
THERMO DATA
METHOD SYSTEM=MARGULES
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=50/40/10
Typical Usage
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=REGULAR
STREAM DATA
. . .
General Information
The REGULAR liquid activity method is used to predict VLE and
VLLE phase behavior. It is generally useful for slightly non-ideal
applications. Refer to the PRO/II Reference Manual for additional
information.
L 2
Vi ( i m )
ln i = ----------------------------------
RT
where:
VLi = Liquid molar volume of component i
i = Solubility parameter of component i
and
N
L
xi Vi
i---------------------------
=1
m = N
L
xi Vi
i=1
Examples
E-6.15: Using the REGULAR system, model a 50/40/10 mix of
DIPE, IPA and water at 50 psia and 50% vaporization.
TITLE PROB=REGULAR
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,DIPE/ 2,IPA/ 3,H2O
THERMO DATA
METHOD SYSTEM=REGULAR
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=50/40/10
UNIT OPERATION
FLASH UID=FLSH
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=FLORY
STREAM DATA
. . .
General Information
The FLORY-HUGGINS liquid activity method is used to predict
VLE and VLLE phase behavior. This method does not support free
water decant.
The FLORY-HUGGINS liquid activity method is generally useful
for mixtures of components that differ vastly in size, e.g., polymer
solutions. Refer to the PRO/II Reference Manual for additional
information.
FLORY I, j, cij / ..
.
xi vi
i
where:
vi = the liquid molar volume of component i.
The number of segments is calculated as:
vi
m i = ----------
v min
where:
vmin = the smallest molar volume in the system.
Missing parameters are estimated as:
v min 2
jk = ---------- ( j k )
RT
where:
j and k are solubility parameters of components j and k,
respectively.
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=NRTL
KVALUE FILL=UNIFAC
STREAM DATA
. . .
Input Description
The METHOD Statement
METHOD SYSTEM=NRTL or UNIQUAC or WILSON or
VANLAAR or MARGULES
...
or
K-value Data
KVALUE(VLE)FILL=NONE or UNIFAC or CIUNIFAC or UFT1 or
REGULAR or FLORY,
AZEOTROPE=SIMSCI or NONE or bankid
or
KVALUE(LLE)FILL=NONE or UNIFAC or CIUNIFAC or UFT1 or
REGULAR or FLORY,
AZEOTROPE=SIMSCI or NONE or bankid
or
KVALUE(SLE)FILL=VANTHOFF or ONE or FREE
Examples
E-6.19: Using the VLLE NRTL liquid activity method, and the
Henrys Law option, determine the solubility of O2 in a
methanol/ethanol/benzene/water/IPA system.
NRTL binary interaction parameters are provided for the
ethanol-IPA and benzene-ethanol pairs. Data are provided
for the azeotrope formed between ethanol and water.
Mutual solubility data are provided between water and
benzene, while infinite dilution data are given for the
methanol-water pair.
The pairs ethanol-methanol and methanol-IPA are defined
as ideal pairs in the liquid phase. Supply data for the O2-
IPA pair via a HENDATA statement, and retrieve data for
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, 02
THERMO DATA
METHOD SYSTEM=NRTL, HENRY
STREAM DATA
. . .
General Information
The HENRY option is used to predict gas solubilities, especially for
modeling supercritical components using a liquid activity method.
It is especially useful for environmental applications such as model-
ing trace organics in aqueous streams.
Hydrogen Chloride
Carbon Monoxide
Cartoon Disulfide
Hydrogen Sulfide
Carbon Dioxide
Sulfur Dioxide
Nitrous Oxide
Liquid
Methanethiol
Hydrogen
Acetylene
Propylene
Ammonia
Isobutane
Ethandiol
N-butane
Components
Ethylene
Nitrogen
Methane
Propane
Oxygen
Helium
Ethane
Argon
Water X X X X X X X X X X X X X X X X X X X X X
Acetic Acid X X X X X
Formic Acid X X X X
m-cresol X X
Methanol X X X X X X X X X X
Ethanol X X X X X X X
n-propanol X X X X X X X
isopropanol X X X X X X X
n-butanol X X X X X X
sec-butanol X X X X
tert-butanol X X X X
3-methyl-1butanol X X X X
n-pentanol X X X X
Ethylene Glycol X X
Triethylene Glycol X X X X X X X X X X
Acetaldehyde X X X X
Methyl ether X X X X
Ethyl ether X X X X
Isopropyl ether X X X X
Tetrahydrofuran X
Dibenzofuran X
Methyl formate X X X X
Methyl acetate X X X X
Ethyl acetate X X X X
Hydrogen Chloride
Carbon Monoxide
Cartoon Disulfide
Hydrogen Sulfide
Carbon Dioxide
Sulfur Dioxide
Nitrous Oxide
Liquid
Methanethiol
Hydrogen
Propylene
Acetylene
Ammonia
Ethandiol
Isobutane
N-butane
Components
Nitrogen
Ethylene
Methane
Propane
Oxygen
Helium
Ethane
Argon
Isopentane X X X X
n-pentane X X X X X
2-methylpentane X X X X
n-hexane X X X X X
2-4 X X X X
dimethylpentane
3-methylhexane X X X X
n-heptane X X X X X X
2-4- X X X X
dimethylhexane
n-octane X X
n-nonane X X
n-decane X X X X X X
n-dodecane X
n-tetradecane X
n-hexadecane X X X X X X X X X
Eicosane X
Docosane X
Dotriacontane X
Cyclopentane X X X X
Methylcyclopentane X X X X
Cyclohexane X X X X
Methylcyclohexane X X X X X
1-T-2-DMCP X X X X
1-T-2-C-4MCP X X X X
Bicyclohexyl X X X X X
1-hexene X X X X
Freon 113 X X
Benzene X X X X X
Chlorobenzene X X
Hydrogen Chloride
Carbon Monoxide
Cartoon Disulfide
Hydrogen Sulfide
Carbon Dioxide
Sulfur Dioxide
Nitrous Oxide
Liquid
Methanethiol
Hydrogen
Propylene
Acetylene
Ammonia
Ethandiol
Isobutane
N-butane
Components
Nitrogen
Ethylene
Methane
Propane
Oxygen
Helium
Ethane
Argon
Toluene X X X X X X
m-xylene X X X
1-C1-naphthalene X X X X X
Diphenylmethane X X X X
1,3,5-tri-C1- X
benzene
9,10 X
dihydrophenanthren
e
Tetralin X
2-ethylanthracene X
9-methylanthracene X
Hydrogen fluoride X X X
Acetone X X X X X X
Methyl ethyl ketone X X X X
Methyl n-butyl X
ketone
Methyl propyl X
ketone
n-methyl-2- X X X X X X
pyrrolidone
Quinoline X X
Sulfolane X X X X
Hydrogen chloride X
Ammonia X X
Chloride X
BANK This option selects one or more databanks from which to retrieve
Henrys Law data. Options available are:
SIMSCI Selects the SIMSCI standard databank (default).
NONE This option disables all data retrieval from databanks for
Henrys Law data.
bankid This option selects a user-created databank named bankid
that is created and maintained with the LIBMGR program.
If the data for any component designated as a Henry component is
missing from the supplied databank and is not supplied via a
HENDATA statement, the Henrys law constant is obtained by
extrapolation of the component vapor pressure data.
SOLUTE This statement is used to enter the component ID numbers for the
components to be treated as solutes.
HENDATA This statement allows entry of the Henry coefficients for solute i in
solvent l. Coefficients C1, C2, C3, and C4 apply to Henrys
Law constant Hi (and the K-value) by the following equations
for one or more solutes in a pure solvent:
C2
ln ( H i ) = C 1 + ------ + C 3 ln ( T ) + C 4 ( P )
T
and Hi
K i = -----
P
For a mixture of solvents, the following mixing rules are used for
determining Hi,mix and Ki:
N
ln ( H i, mix ) = xl ln ( Hi, l )
l=1
and Hi
K i = -----
P
Examples
E-6.20: Using the NRTL liquid activity method, and the Henrys
Law option determine the solubility of CO in a methanol/
water system. Data for the CO-methanol pair are supplied
via a HENDATA statement, while data for the other pairs
are retrieved from the SIMSCI databank.
TITLE PROB=HENRY
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1,N2/ 2,CO/ 3,H2O/ 4,MEOH
THERMO DATA
METHOD SYSTEM=NRTL, HENRY
HENRY
SOLUTE 1, 2
HENDATA 2, 4, 152.4, -8000, -20, 0
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=2/3/70/25
UNIT OPERATION
FLASH UID=FLSH
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
...
COMPONENT DATA
LIBID 1, ACETIC/ 2, ACRYLIC/ 3, H2O
THERMO DATA
METHOD SYSTEM=NRTL, PHI=HOCV, &
ENTHALPY(V)=HOCV, DENSITY(V)=HOCV
KVALUE FILL=UNIFAC
STREAM DATA
. . .
General Information
The HOCV method predicts vapor fugacities, vapor enthalpies,
vapor entropies and vapor densities. It is especially useful for sys-
tems where dimers form in the vapor phase, e.g., carboxylic acid
systems. A liquid activity method must be used in conjunction with
the HOCV method.
Examples
E-6.21: Using the HOCV vapor fugacity method with the NRTL
liquid activity method, model a 50/50 mix of acetic and
acrylic acids at 100 psia and 100 F.
TITLE PROB=HOCV
COMPONENT DATA
LIBID 1,ACETIC/ 2,ACRYLIC
THERMO DATA
METHOD SYSTEM=NRTL, PHI=HOCV, &
ENTHALPY(V)=HOCV, ENTROPY(V)=HOCV, &
DENSITY(V)=HOCV
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=50/50
UNIT OPERATION
FLASH UID=FL1
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
E-6.22: For the same problem, use HOCV for vapor fugacities and
SRK for vapor enthalpies, entropies and densities. Specify
API liquid densities.
THERMODYNAMIC DATA
METHOD SYSTEM=NRTL, PHI=HOCV, &
ENTHALPY(V)=SRK, ENTROPY(V)=SRK, &
DENSITY(V)=HOCV, DENSITY(L)=API
...
COMPONENT DATA
LIBID 1, ACETIC/ 2, ACRYLIC/ 3, H2O
THERMO DATA
METHOD SYSTEM=NRTL, PHI=TVIRIAL
STREAM DATA
. . .
General Information
The TVIRIAL method predicts vapor fugacities. It is useful for sys-
tems where dimers form in the vapor phase, e.g., carboxylic acid
systems. A liquid activity method must be used in conjunction with
the TVIRIAL method.
Examples
E-6.23: Using the TVIRIAL vapor fugacity method with the
NRTL liquid activity method, model a 75/25 mix of
ethanol and cyclohexane at 5 atm and 115 C.
TITLE PROB=TVIRIAL
PRINT INPUT=ALL
DIME METRIC, PRES=ATM, TEMP=C
COMPONENT DATA
LIBID 1,ETHANOL/ 2,CYHX
THERMO DATA
METHOD SYSTEM=NRTL, PHI=TVIRIAL
KVALUE POYNTING=ON
NRT6(K) 1, 2, -.1429, 472.71, -.01, 748.33, 0.3681, 2.871E-4
PHI
TVIRIAL 1, 1.0/ 2, 0.0
STREAM DATA
PROP STREAM=1, TEMP=115, PRES=5, &
COMP=75/25
UNIT OPERATION
FLASH UID=FL1
FEED 1
E-6.24: For the same problem, use TVIRIAL for vapor fugacities
and SRK for vapor enthalpies, entropies and densities.
Specify API liquid densities.
THERMODYNAMIC DATA
METHOD SYSTEM=NRTL, PHI=TVIRIAL, &
ENTHALPY(V)=SRK, ENTROPY(V)=SRK, &
DENSITY(V)=SRK, DENSITY(L)=API
...
COMPONENT DATA
LIBID 1, ACETIC/ 2, ACRYLIC/ 3, H2O
THERMO DATA
METHOD SYSTEM=NRTL, PHI=IDIMER, &
ENTHALPY(V)=IDIMER, DENSITY(V)=IDIMER
KVALUE FILL=UNIFAC
STREAM DATA
. . .
General Information
The IDIMER method predicts vapor fugacities, vapor enthalpies,
vapor entropies and vapor densities. It is especially useful for sys-
tems where dimers form in the vapor phase, e.g., carboxylic acid
systems. A liquid activity method must be used in conjunction with
the IDIMER method.
Input Description
The METHOD Statement
METHOD SYSTEM=NRTL or UNIQUAC or UNIFAC or UFT1 or
UFT2 or UFT3 or UNFV or WILSON or VANLAAR or
MARGULES or REGULAR or FLORY, PHI=IDIMER,
ENTHALPY(V)=IDIMER, ENTROPY(V)=IDIMER,
DENSITY(V)=IDIMER, . . .
or
METHOD KVALUE=NRTL or UNIQUAC or UNIFAC or UFT1 or
UFT2 or UFT3 or UNFV or WILSON or VANLAAR or
MARGULES or REGULAR or FLORY, PHI=IDIMER,
ENTHALPY(V)=IDIMER, ENTROPY(V)=IDIMER,
DENSITY(V)=IDIMER, . . .
BANK This option selects one or more databanks from which to retrieve
vapor phase binary interaction data.
SIMSCI Selects the SIMSCI standard databank (default).
NONE This option disables all data retrieval from
databanks for interaction parameters.
bankid This option selects a user-created databank
named bankid that is created and maintained
with the LIBMGR program.
IDIMER This statement supplies interaction data for the IDIMER vapor
fugacity method. This method may be used for components
forming dimers, particularly carboxylic acids. Aii and Bii are the
association parameter for component i, while Aii and Bii are the
solvation parameter for components i and j.
Examples
6.25: Use the IDIMER vapor fugacity method with the NRTL
liquid activity method for a mix of associating acids.
TITLE PROB=IDIMER
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1, HFOR/ 2, ACETIC/ 3, HPRP
THERMO DATA
METHOD SYSTEM=NRTL, PHI=IDIMER, &
ENTHALPY(V)=IDIMER, ENTROPY(V)=IDIMER, &
DENSITY(V)=IDIMER
PHI BANK=SIMSCI
IDIME(K, MMHG, LOG) 1,1,-10.743,3083.0
IDIME(K, MMHG, LOG) 2,2,-10.421,3166.0
IDIME(K, MMHG, LOG) 3,3,-10.843,3316.0
IDIME(K, MMHG, LOG) 1,2,-10.356,3193.0
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=700, &
COMP=0.477/0.094/0.729
UNIT OPERATION
FLASH UID=FL1
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
...
COMPONENT DATA
LIBID 1, IPA/ 2, H2O/ 3, CH
THERMO DATA
METHOD SYSTEM=NRTL
ENTHALPY(L) HMIX=RK1
KVALUE FILL=UNIFAC
STREAM DATA
. . .
General Information
The Redlich-Kister and Gamma heat of mixing methods apply a
correction to IDEAL enthalpy data. A liquid activity method must
be used in conjunction with the Redlich-Kister or Gamma methods.
Input Description
The METHOD Statement
METHOD SYSTEM=NRTL or UNIQUAC or UNIFAC or UFT1 or
UFT2 or UFT3 or UNFV or WILSON or VANLAAR or
MARGULES or REGULAR or FLORY
or
METHOD KVALUE=NRTL or UNIQUAC or UNIFAC or UFT1 or
UFT2 or UFT3 or UNFV or WILSON or VANLAAR or
MARGULES or REGULAR or FLORY, . . . ,
ENTHALPY(L)=IDEAL, . . .
Enthalpy Data
ENTHALPY BANK=SIMSCI or NONE or bankid,
HMIX=NONE or RK1 or RK2 or GAMMA
RK1(K or KCAL or KJ) i, j, aij, bij, cij, dij, eij, fij, gij, hij/ . . .
or
RK2(K or KCAL or KJ) i, j, aij, bij, cij, dij, eij, fij, gij, hij/ . . .
Note: UOM for WILSON is K or KCAL/KMOL or KJ/KMOL.
BANK This option selects one or more databanks from which to retrieve vapor
phase binary interaction data.
SIMSCI Selects the SIMSCI standard databank (default).
NONE This option disables all data retrieval from databanks for
interaction parameters.
bankid This option selects a user-created databank named bankid
that is created and maintained with the LIBMGR program.
HMIX This statement selects the heat of mixing method to be applied to the liquid
enthalpy method.
NONE No mixing enthalpy correction is applied (default).
RK1 or This option selects different formulations of the Redlich-
RK2 Kister correlation.
where
X = xi xj (form 1 of Redlich-Kister)
xi xj (form 2 of Redlich-Kister)
X = --------------
xi + xj
Examples
6.26: Using the NRTL liquid activity method, and correcting for
the heat of mixing using the Redlich-Kister method RK1,
model a 50/50 mix of acetic and acrylic acids at 100 psia
and 100 F.
TITLE PROB=RK1
PRINT INPUT=ALL
COMPONENT DATA
LIBID 1, ACETIC/ 2, ACRYLIC
THERMO DATA
METHOD SYSTEM=NRTL, PHI=HOCV
ENTHALPY(L) HMIX=RK1
STREAM DATA
PROP STREAM=1, TEMP=100, PRES=100, &
COMP=50/50
UNIT OPERATION
FLASH UID=FL1
FEED 1
PROD V=2, L=3
TPSPEC PRES=50
SPEC STREAM=2, RATE, RATIO, STREAM=1, &
RATE, VALUE=0.5
END
6.27: For the same problem, use TVIRIAL for vapor fugacities
and specify API liquid densities.
THERMODYNAMIC DATA
METHOD SYSTEM=NRTL, PHI=TVIRIAL, DENSITY(L)=API
ENTHALPY(L) HMIX=RK1
General Information
The VANTHOFF solubility method is used to calculate solid-liquid
equilibrium K-values for nearly ideal non-electrolyte systems using
the van't Hoff ideal solubility equation.
Input Description
The METHOD Statement
METHOD KVALUE(SLE)=VANTHOFF, ...
where:
Hm= enthalpy change of melting at Tt
7.2: For the same problem, use the VANLAAR method for
VLE calculations.
THERMO DATA
METHOD SYSTEM(VLE)=VANLAAR, &
KVALUE(SLE)=VANTHOFF, TRANSPORT=PURE
General Information
The SOLDATA solubility method is used to calculate solid-liquid
equilibrium K-values using user-supplied solubility data.
KVALUE Selects the method for K-value calculations. Only SLE K-value
calculations are available with the SOLDATA method.
K-value Data
KVALUE(SLE) FILL=VANTHOFF or ONE or FREE
SOLUTE i, j, k,...
SOLDATA(utemp) i, l, c1, c2, c3 { , tmin, tmax}.
FILL The fill option selects the method used to estimate binary interaction
data missing from the selected databank and not given via SOLDATA
statements. The options available are:
VANTHOFF van't Hoff ideal solubility equation. This is the default.
ONE Missing binary solubilities set to 1.0 (i.e., all
unspecified solids completely soluble).
FREE If a solvent k is missing binary solubility data for a
solute, the solvent is ignored in the solubility
calculation. That is, the solvent mole fractions are
normalized to a k-free basis.
SOLUTE Specifies the solutes i, j, k ... in the system.
SOLDATA This statement supplies the molar solubility of a solute
i in a mixture of i and solvent l, as a function of
temperature according to the following equations:
c2
ln x il = c 1 + ----- + c 3 ln T
T
ln ( x i, mix ) = xl ln ( xi, l )
l = 1
li
7.4: For the same problem, use SRK for VLE calculations and
set all missing binary interaction data for SLE calculations
equal to 1.0.
THERMODYNAMIC DATA
METHOD SYSTEM=SRK, KVALUE(SLE)=SOLDATA, &
TRANSPORT=PURE
KVALUE(SLE) FILL=ONE
Transport Properties
Typical Usage
...
COMPONENT DATA
LIBID 1,C3/ 2,IC4/ 3,NC4
THERMO DATA
METHOD SYSTEM=SRK, TRANSPORT=PURE
STREAM DATA
. . .
General Information
The TRANSPORT keyword is used to provide transport properties,
including liquid and vapor viscosities, liquid and vapor thermal
conductivities, and liquid surface tension values. Liquid diffusivi-
ties may be computed by selecting the DIFFUSIVITY(L) keyword.
TRANSPORT
The TRANSPORT keyword selects the method set used to calculate of
transport properties including liquid and vapor viscosities, liquid and vapor
thermal conductivities, and liquid surface tension values. When the
TRANSPORT keyword is absent, no transport calculations occur. Begin-
ning with PRO/II 9.1, when the thermodynamic system is Tacite, the
default transport method is TACITE. When the TRANSPORT keyword is
present, the available options are:
Viscosity Data
VISCOSITY(VL or V and/or L)
This keyword requests the calculation of vapor and liquid viscosities. The
options available with this keyword are:
KVIS This option calculates the viscosity from the values of the
kinematic viscosity and the density. For purposes of this
calculation, the density is computed with the API method.
In order to use this option, a KVIS method must be declared
on the METHOD statement. This option is only available
with the VISCOSITY(L) keyword.
LBC This is the Lohrenz-Bray-Clark prediction method for
calculating liquid viscosities.
U1-U5 This option selects one of the 5 supported user-defined
subroutines that are available to compute transport
properties. This option is available for both vapor and
liquid viscosity calculations.
SIMSCI This option uses SimSci developed liquid viscosity values.
This was known as the Twu method in PROCESS. This
option is available only with the VISCOSITY(L) keyword.
API This option uses liquid viscosities from the API Technical
Data Book. This option is available only with the
VISCOSITY(L) keyword.
BW Bromley-Wilkey correlation (Procedure 11B2.1) suitable
for gaseous mixtures at low pressure. Available only for
vapor viscosity calculations.
TWUB Applies the recommended Twu-Bulls mixing rule in the
API Technical Data Book to the SIMSCI (Twu) method.
This option is available for liquid viscosity.
APIT For liquid viscosities only, this option applies the Twu-
Bulls mixing rule to the viscosity nomograph in the API
Technical Data Book.
Conductivity Data
CONDUCTIVITY (VL or V and/or L)
This keyword requests the calculation of vapor and liquid thermal
conductivities. The options available with this keyword are:
Diffusivity Data
DIFFUSIVITY (L) This specifies a diffusivity method. Diffusivity is not
calculated when this entry is missing. Available method
options are:
WILKE This argument specifies using the Wilke-Chang correlation.
It is the default method used when the DIFFUSIVITY(L)
entry is present without the argument that identifies the
calculation method.
DIFDATA This option requires the user to provide DIFFUSIVITY and
DIFDATA statements after the METHOD statement. The
DIFDATA statement supplies the data to use in the
diffusivity calculations.
DIFDATA(utemp)i, j, c1, c2, c3/ ...
This statement supplies the diffusivity of a solute i in a mixture of
components i and j.
Diffusivity is a measure of the rate at which a solute diffuses through
a given area in a given time period under a concentration gradient.
The dimensions of diffusivity are always in m2/sec. The diffusivity is
given as a function of temperature according to the following
equation:
ln Dij = c1 + c2 / T + c3 ln T
where T is in absolute units (K or R).
Examples
Example 8.1: In this example, of using a METHOD statement, an
HCURVE unit operation is used to generate heating
curves and transport property tables for a crude feed
stream.
TITLE PROBLEM=TRANSPORT
DIMEN ENGLISH, LIQV=BBL
PRINT RATE=M,STREAM=ALL,INPUT=ALL,TBP
COMPONENT DATA
LIBID 1,ETHANE/ 2,PROPANE/ 3,IBUTANE/ 4,BUTANE/ 5,PENTANE
TBPCUTS 115,300,6/ 400,10/ 650,8/ 800,4/ 1500,6
THERMODYNAMIC DATA
METHOD KVAUE(VLE)=BK10, ENTH(V)=JG, ENTH(L)=JG, &
ENTR(V)=CP, ENTR(L)=CP, DENS(V)=IDEA, &
DENS(L)=API, TRANS=PETRO
STREAM DATA
PROP STREAM=1,TEMP=375, PRES=300, PHASE=M, &
RATE(V)=3125, ASSAY=LV
D86 STREAM=1, PRES(MMHG)=760, TEMP=F, &
DATA=8,135/25,210/43,370/67,565/75,665/82,800/92,990
API STREAM=1, AVG=45.37, &
DATA=11.6,80.01/21.6,62.9/41.7,50.6/61.9,38.2/83.8,27.5
MW STREAM=1, AVG=162.9, &
DATA=24.2,99.5/40.5,135/55,184.7/74.8,334.8/100,789
LIGHT STREAM=1, PERCENT(W)=10.4, &
COMP(M)=1,0.1/2,1.4/3,0.65/4,3.15/5,5.1, NORMALIZE
UNIT OPERATIONS
HCURVE UID=HC1, NAME=HEATING CRV
ADIA STREAM=1, TEMP=375, 690, PRES=300, 50, POINTS=20
PROP TRANSPORT
END
Typical Usage
...
COMPONENT DATA
LIBID 1,IC4/ 2,NC4/ 3,NC5
THERMO DATA
METHOD SYSTEM=system, ... propkw(options)=method ...
Optional method-specific property data statements
property(qualifier, basis) {GAMMA=value, REFINDEX=value &
REFVALUE(uprop)=propvalue, NCFILL=ncfill, &
NCBLEND=ncblend}
{ DATA i, datavalue {/ i, datavalue ...} or
INDEX i, indexvalue {/ i, indexvalue ...}}
STREAM DATA
. . .
General Information
The program provides numerous methods to compute special
stream refinery properties. These properties include KVIS which cal-
culates liquid kinematic viscosities, CLOUD and FLPO (CC,OC)
(cloud and flash point temperatures respectively), and SULF (sulfur
content). To include a special property in a simulation, use the
appropriate key word entry on the METHOD statement to designate
a calculation method. Use additional property statements to supply
optional data for the calculations. The keyword associated with
each stream refinery special properties available are given in Table
8-4.
Table 8-4: Keywords for Special Stream Refinery Properties
Keyword Property
KVIS Kinematic viscosity
POUR Pour point temperature
CLOU Cloud point temperature
FRZP Freeze point temperature
CFPP Cold filter plug point
FLPO (CC,OC) Flash point temperature (formerly FLPT) OC
selects open cup, CC is closed cup)
RON Research octane number
MON Motor octane number
IBP Initial boiling point
Note: The D1740 method calculates the Luminometer number (LUMI) from
the Smoke Point (SMOK). Therefore, it requires the SMOK entry to be
provided on the METHOD statement and data supplied. Otherwise, the
LUMI number is reported as missing.
SULF FRAC or PCT or M or WT INDEX or SUM or
LV or LV USINDEX or USFORM
ASUL WT SUM or USFORM
MERC WT SUM or USFORM
NITR TOTA or BASI or WT INDEX or SUM or
NONB USINDEX or USFORM
OXYG WT INDEX or SUM or
USINDEX or USFORM
CCR M or WT INDEX or SUM or
or LV USINDEX or
USFORM
NICK WT INDEX or SUM or
USINDEX or USFORM
SODI WT INDEX or SUM or
USINDEX or USFORM
VANA WT INDEX or SUM or
USINDEX or USFORM
IRON WT SUM or USFORM
ASPH C5 or C7 WT INDEX or SUM or
USINDEX or USFORM
ANEU WT or LV SUM or USFORM
BROM WT or LV SUM or USFORM
MEAB API
CABP API
MOAB API
VABP API
WABP API
NHV API
WAX WT INDEX or SUM or
USINDEX or USFORM
ASH WT INDEX or SUM or
USINDEX or USFORM
NOAC WT INDEX or SUM or
USINDEX or USFORM
For example, when specifying the SIMSCI mixing method for KVIS
using the default liquid volume basis, the entry on the METHOD
statement becomes:
METHOD ..., KVIS=SIMSCI, ...
When specifying the INDEX mixing method for leaded research
octane number (RON) and using the default weight basis, the RON
entry on the METHOD statement becomes:
METHOD ..., RON(L, WT)=INDEX, ...
ln ( KVIS )
Index = -------------------------------------------
ln ( 1000 KVIS )
Index ln ( 1000 )
KVIS = exp --------------------------------------------
( 1 Index )
MAXWELL Converts kinematic viscosity in centiStokes.to an index that is
blended on a liquid volume basis.
Usage Syntax
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM= system, property=propmethod, ...
property1 statement (see Table 8-6)
compProp statement (see Table 9-1)
...
property2 statement (see Table 8-6)
compProp statement (see Table 9-1)
STREAM DATA
...
Typical Usage
...
THERMO DATA
METHOD SYSTEM=SRK
KVALUE . . .
TC 1, 373 . . .
STREAM DATA
. . .
Input Description
Supplying Pure Component Data In A METHOD Set
...
COMPONENT DATA
LIBID 1, IC4/ 2, NC4/ 3, NC5
THERMO DATA
METHOD SYSTEM= system, property=propmethod, ...
property1 data statement (see Table 8-6)
{DATA . . . or INDEX . . . }
compProp1 data statement (see Table 9-1)
...
compPropn data statement (see Table 9-1)
...
property2 statement (see Table 8-6)
{DATA . . . or INDEX . . . }
compProp1 data statement (see Table 9-1)
...
compPropn data statement (see Table 9-1)
...
STREAM DATA
...
Example:
A DENSITY(V)=SRK entry declares the vapor density calculation
method to be the Soave modified Redlich-Kwong equation of state
(SRK).
Table 9-1 shows the only pure component data used by the SRK
generator are critical temperature, critical pressure, and acentric
factor (Tc, Pc, and ). Table 9-2 associates a component data state-
ment keyword with each pure component property listed in Table 9-
1. Here are the associated keywords.
Component Data Symbol Keyword
critical temperature Tc TC
critical pressure Pc PC
acentric factor ACENTRIC
It should be clear that these component properties are the only ones
that apply when computing density with this generator. A partial
keyword listing supplying pure component data for components 1
and 2 could be:
...
THERMODYNAMIC DATA
METHOD . . . DENSITY(V)=SRK
DENSITY(V) . . .
Note the pure component data statements appear after the DENSITY
property statement. The pure component data is used only when
computing vapor density with this method. A different set of pure
component data could be supplied for computing vapor density in
another METHOD set.
Hierarchy
The hierarchy that governs the use of pure component data is:
1. Any method-specific data supplied on the data statements for
each property in each method set. This includes both pure com-
ponent data and property data statements. Only temperature-
independent data may be given in the Thermodynamic Data
Category.
2. Applicable default values supplied on data statements follow-
ing a KVALUE statement in the same method set.
3. Data supplied in the COMPONENT DATA section of keyword
input.
4. Data stored in the component libraries as selected using the
BANK option.
If data are still missing after steps 1-4 are completed, PRO/II issues
an error message and terminates after input processing is complete.
The PROVISION GUI displays the error messages and opens the
keyword input file for correction by the user. .
Table 9-2: Data Card Keywords To Use for Property Headers in Table 9-1
Header Keyword Description Header Keyword Description
Tc TC Critical temperature VISCOSITY Liquid Viscosity
vdW VANDERWAAL van der Waals area and CONDUCT Thermal conductivity
volume
NBP NBP Normal boiling point C Users may supply these data only in the
(temperature) Component Data Category.
H0 ENTHALPY (V) Ideal vapor enthalpy Exponent EXPO Constant for surface
tension calculation
Parachor PCHO Parachor for tacite Peneloux PENE Peneloux for Tacite
surface tension liquid density
Hvap LATENT Latent heat of A The user may supply these data or
vaporization index values either globally in the
Hliq ENTHALPY (L) Saturated liquid Component Data or Stream Data
enthalpies Categories or locally in the Thermo-
VP VP Vapor pressure dynamic Data Category for a given
method sets.
SpGr SPGR Specific gravity
Index- 2
Examples for Liquid Activity Methods , 6-55 GPA Sour Water
General Information for Liquid Activity Attributes of GPSWATER Package , 5-23
Methods , 6-53 Enthalpy, Entropy, and Density Data , 5-24
K-value Data for Cubic EOS , 4-28 Examples , 5-25
K-value Data for Liquid Activity Methods , 6-54 METHOD Statement , 5-23
METHOD Statement for Cubic EOS , 4-27 Method-specific Pure Component Properties , 5-25
METHOD Statement for Liquid Activity Typical Usage , 5-22
Methods , 6-53
Grayson-Streed
Typical Usage for Cubic EOS , 4-27
Attributes of , 3-5
Typical Usage for Liquid Activity Methods , 6-52
METHOD Statement , 3-5
Flory-Huggins Typical Usage , 3-5
Attributes of FLORY Methods , 6-48 Water Handling Options , 3-6
Examples , 6-52
Grayson-Streed Method
K-value Data , 6-49
Examples , 3-6
METHOD Statement , 6-49 Method-specific Pure Component Properties , 3-6
Method-specific Pure Component Properties , 6-51
Typical Usage , 6-48
Free-water Decant Considerations , 1-33 H
Examples , 1-38
General Information , 1-33 Hayden-OConnell Vapor Fugacity
Water Handling Statement , 1-35 Attributes of the HOCV Vapor Fugacity
Method , 6-64
Enthalpy, Entropy and Density Data , 6-66
G Examples , 6-67
K-value Data , 6-65
Generalized Correlations METHOD Statement , 6-65
API Liquid Density , 3-26 Method-specific Pure Component Properties , 6-67
Braun K10 , 3-18 Typical Usage , 6-64
Chao-Seader , 3-8 Vapor Fugacity Data , 6-66
Costald Liquid Density , 3-30 Henrys Law
Curl-Pitzer , 3-15 Attributes of the Henrys Law Option , 6-57
Grayson-Streed , 3-5 Examples , 6-63
Ideal and Library Methods , 3-1 Henrys Law Data , 6-62
Ideal Methods , 3-1 Henrys Law Package Available Pairs , 6-59
Johnson-Grayson , 3-21 METHOD Statement , 6-58
Lee-Kesler , 3-23 Method-specific Pure Component Properties , 6-63
Library Methods , 3-1 Typical Usage , 6-57
Modifications to GS and CS , 3-11 Hierarchy of Pure Component Data , 9-5
Rackett Liquid Density , 3-28
High Pressure Crude Systems , 2-3
Typical Usage , 3-1
Hydrofiners , 2-4
User-supplied K-value Data , 3-34
Glycol Dehydration Systems , 2-6
Glycols
Attributes of the GLYCOL Package , 5-14
I
Components Available for GLYCOL Package , 5- Ideal Method
17
Attributes of , 3-2
Enthalpy, Entropy, and Density Data , 5-16
METHOD Statement , 3-2
Examples , 5-18
Ideal Methods
K-value Data , 5-15
Examples , 3-3
METHOD Statement , 5-15
Method-specific Pure Component Properties , 5-18 IDIMER Vapor Fugacity
Attributes of , 6-71
Typical Usage , 5-14, 5-19
Enthalpy, Entropy and Density Data , 6-73
Index- 4
Attributes of Modified UNIFAC Methods , 6-20 Typical Usage , 6-1
Examples , 6-32
K-value Data , 6-22
METHOD Statement , 6-21 O
Method-specific Pure Component Properties , 6-
26, 6-31 Order of Input , 1-16
Typical Usage , 6-20
UFT1-Lyngby modified UNIFAC , 6-22
UFT2-Dortmund modified UNIFAC , 6-23 P
UFT3-Modified UNIFAC , 6-23
UNFV-Free Volume modification , 6-23 Peng-Robinson , 4-7
Modified SRK and PR Enthalpy, Entropy, and Density Data , 4-10
Enthalpy, Entropy, and Density Data , 4-18 Examples , 4-11
Examples , 4-19 K-value Data , 4-9
Huron-Vidal modifications to PR , 4-15 METHOD Statement , 4-8
Huron-Vidal modifications to SRK , 4-14 Method-specific Pure Component Properties , 4-11
Kabadi-Danner modifications to SRK , 4-15 Typical Usage , 4-7
K-value Data , 4-16 Water Handling Options , 4-11
METHOD Statement , 4-14 Petrochemical Applications
Method-specific Pure Component Properties , 4-19 Alcohol Dehydration Systems , 2-12
Modified Panagiotopoulos-Reid modifications to Aromatic Systems , 2-9
PR , 4-15 Aromatic/Non-aromatic Systems , 2-9
Modified Panagiotopoulos-Reid modifications to HF Systems , 2-12
SRK , 4-14
Light Hydrocarbon Systems , 2-8
Panagiotopoulos-Reid modifications to PR , 4-15
Non-hydrocarbon Systems , 2-10
SimSci modifications to SRK , 4-14
Predefined Systems , 1-26
Typical Usage , 4-13
Property Statements
Multiple Thermodynamic Sets , 1-30
Density Options , 1-7
Avoiding Inconsistent Methods , 1-31
Diffusivity Options , 1-6
Default Method Set , 1-30
Enthalpy Options , 1-8
Examples , 1-32
Entropy Options , 1-8
General Information , 1-30
Henrys Law Options , 1-6
Method Set Requirements , 1-31
Liquid-Liquid Equilibrium Options , 1-6
Multiple Method Sets in Distillation Columns , 1-
31 Solid-Liquid Equilibrium Options , 1-6
Vapor Fugacity Options , 1-6
Vapor-Liquid Equilibrium Options , 1-5
N Pure Component Alpha Formulations , 1-12
Index- 6
Method-specific Pure Component Properties , 6-17 General Information , 7-1
Typical Usage , 6-12 METHOD Statement , 7-2
UNIQUAC Typical Usage , 7-1
Attributes of UNIQUAC Methods , 6-7 Vapor-liquid-liquid Equilibrium Considerations , 1-38
Examples , 6-11 Available Databanks , 1-40
K-value Data , 6-8 Declaring a VLLE set , 1-38
METHOD Statement , 6-8 Examples , 1-44
Method-specific Pure Component Properties , 6-11 K-value Data , 1-44
Typical Usage , 6-7 METHOD Statement , 1-43
UNIWAALS , 4-22 Specifying Key Components , 1-42
Attributes of , 4-22 Specifying Separate VLE and LLE Models , 1-42
Enthalpy, Entropy, and Density Data , 4-25 VLLE Predefined Systems and K-value
Examples , 4-25 Generators , 1-39
K-value Data , 4-23
METHOD Statement , 4-23
Method-specific Pure Component Properties , 4-25 W
Typical Usage , 4-22
Water Handling , 2-2
User-added Subroutine Data , 1-13
User-added Subroutines
Attributes of , 5-32
Enthalpy, Entropy, and Density Data , 5-33
Examples , 5-33, 5-35
K-value Data , 5-33
METHOD Statement , 5-32
Typical Usage , 5-32, 5-34
User-supplied K-value Data , 1-9
Attributes of , 3-32, 3-35
Examples , 3-36
Input Description , 3-35
K-value Data (required) , 3-35
METHOD Statement , 3-35
Typical Usage , 3-34
User-supplied Solubility
Examples , 7-6
General Information , 7-4
K-value Data , 7-5
METHOD Statement , 7-5
Typical Usage , 7-4
V
Van Laar
Attributes of VANLAAR Methods , 6-33
Examples , 6-36
K-value Data , 6-34
METHOD Statement , 6-33
Method-specific Pure Component Properties , 6-36
Typical Usage , 6-32
Vant Hoff Solubility
Examples , 7-3