AspenPolymers ExsAppsV7 3 Ref
AspenPolymers ExsAppsV7 3 Ref
AspenPolymers ExsAppsV7 3 Ref
Aspen Polymers, Aspen Custom Modeler, Aspen Dynamics, Aspen Plus, Aspen Properties, aspenONE, the
aspen leaf logo and Plantelligence and Enterprise Optimization are trademarks or registered trademarks of Aspen
Technology, Inc., Burlington, MA.
All other brand and product names are trademarks or registered trademarks of their respective companies.
This document is intended as a guide to using AspenTech's software. This documentation contains AspenTech
proprietary and confidential information and may not be disclosed, used, or copied without the prior consent of
AspenTech or as set forth in the applicable license agreement. Users are solely responsible for the proper use of
the software and the application of the results obtained.
Although AspenTech has tested the software and reviewed the documentation, the sole warranty for the software
may be found in the applicable license agreement between AspenTech and the user. ASPENTECH MAKES NO
WARRANTY OR REPRESENTATION, EITHER EXPRESSED OR IMPLIED, WITH RESPECT TO THIS DOCUMENTATION,
ITS QUALITY, PERFORMANCE, MERCHANTABILITY, OR FITNESS FOR A PARTICULAR PURPOSE.
Contents iii
Defining Molecular Structure.................................................................. 32
Specifying Mass Fraction Crystallinity...................................................... 33
Creating Property Sets .......................................................................... 33
Creating Property Tables ....................................................................... 34
Running the Simulation and Examining the Results............................................. 35
Input Summary ................................................................................... 37
References .................................................................................................... 38
iv Contents
Specifying and Characterizing Components.............................................. 84
Specifying Physical Properties ................................................................ 85
Supplying Process Information ............................................................... 86
Creating a Sensitivity Table ................................................................... 87
Running the Simulation and Examining the Results............................................. 89
Input Summary ................................................................................... 89
References .................................................................................................... 91
Contents v
B5 Styrene Butadiene Emulsion Copolymerization Process ................................143
About This Process ....................................................................................... 143
Process Definition......................................................................................... 144
Process Conditions ............................................................................. 144
Physical Property Models and Data ....................................................... 146
Reactors / Kinetics ............................................................................. 146
Process Studies.................................................................................. 147
Selected Simulation Results........................................................................... 151
References .................................................................................................. 153
vi Contents
B10 Methyl Methacrylate Polymerization in Ethyl Acetate ..................................219
About This Process ....................................................................................... 219
Process Definition......................................................................................... 219
Process Conditions ............................................................................. 220
Reactors / Kinetics ............................................................................. 221
Process Studies.................................................................................. 221
Selected Simulation Results........................................................................... 223
References .................................................................................................. 228
Index ..................................................................................................................248
Contents vii
viii Contents
Introducing Aspen Polymers
Note: Some of the simulations described in this book require Fortran files.
These are found in the same location as the simulation files (\xmp or \app
sub-directory of the GUI installation directory). You must compile the
accompanying Fortran before running the files.
Installation Guides
Aspen Engineering Suite Installation Guide
Technical Support
AspenTech customers with a valid license and software maintenance
agreement can register to access the online AspenTech Support Center at:
http://support.aspentech.com
This Web support site allows you to:
Access current product documentation
Search for tech tips, solutions and frequently asked questions (FAQs)
Search for and download application examples
Search for and download service packs and product updates
Submit and track technical issues
Send suggestions
Report product defects
Review lists of known deficiencies and defects
Defining Unit-Sets
This example requires a user defined unit-set for input data and output
results.
To define a unit-set:
1 From the Setup folder in the Menu Tree, double-click the Units-Sets sub-
folder.
2 Click New to create a new set.
3 Enter an ID (for example, SET1) and click OK.
A dialog box appears requesting approval to make SET1 the global unit
set, click Yes.
4 Confirm that SI is entered in the Copy from field. If it is not selected,
select SI from the pulldown list. This means the simulation will use SI
units as the basis for your new unit set.
5 Set the following options:
o Mass Flow = kg/hr
o Mole Flow = kmol/hr
o Temperature = C
o Pressure = atm
Specifying
ing Other Simulation Option
ions
To enter other sim
simulation options for this example:
1 From the Setu
Setup folder in the Menu Tree, click the Simulat
lation Options
form.
2 Click the Limi
imits tab.
3 Enter 1000 in the Simulation time limit in CPU seconds
ds field.
4 Click the Syst
System tab, and click the Print Fortran traceba
ebacks when a
Fortran erro
rror occurs option.
Specifyi
ifying Components
You use Compon
ponents forms to select chemical components for
fo your
simulation and sp
specify component types (for example, convent
entional, solid,
assay, blend, poly
olymer, segment, oligomer, and pseudocompononent).
To select compon
onents:
1 From the Data menu, click Components.
2 On the Select
lection sheet, enter the components as shown:
4 From the Menu Tree, double-click the FEEDB sub-folder. Click Input.
5 On the Specifications sheet, enter:
o Temperature =70 C
o Pressure =2 atm
o Total Flow (Mass) = 24000 kg/hr
o Composition = Mass Frac
AIBN = 0.004
STY = 0.1
ACN = 0.4
XYLENE = 0.496
Specifying Kinetics
Kinetic inputs are specified in the Reactions folder. This example uses the
free-radical kinetics model.
To specify the free-radical polymerization inputs:
1 From the Menu Tree, double-click the Reactions folder.
2 From the Reactions folder, click the Reactions sub-folder.
3 Click New to create a new reaction.
4 In the Enter ID field, type R1, in the Select type field, select FREE-
RAD. Click OK.
5 On the Species sheet, enter information as shown:
Modifying Reactions
To review and edit the reactions:
1 Click the Reactions tab.
An Auto Generation dialog box appears.
2 Click Yes to turn off reaction generation.
3 On the Reactions sheet, review the automatically generated reactions
and delete any unnecessary equations by selecting them in the table and
clicking the Delete button . For this example, you can delete all
reactions for:
o Chat-agent
o Chat-sol
o Term-dis
4 Click the Rate Constants tab and enter reaction rate constants as shown
for each reaction:
6 Click Close.
Note: To enter data in the Sentence ID field you can also use the right
mouse button to click in the field, then click Create from the menu.
4 Click OK to accept the ID 1. Enter the information for the gel effect
correlation as shown:
Repeat the above procedures for Specifying Kinetics, Modifying Reactions, and
Specifying Gel Effect to specify the kinetics for R2. Refer to the Input
Summary REACTIONS paragraph for R2 on page 20 for the appropriate
comment, reactions, rate constants, and gel-effect correlations.
Improving Convergence
Polymer reaction kinetics present very difficult convergence problems, so the
standard convergence options are frequently insufficient. To resolve this
problem, we use the RCSTR block Convergence form. To do this:
1 From Blocks folder, B1 sub-folder in the Menu Tree, click the
Convergence form.
2 Click the Parameters tab and change the Solver to Newton.
3 Click Newton Parameters and change the Stabilization Strategy to
Line-Search. Click Close.
Note: You can also open a control panel from the Aspen Plus toolbar, by
clicking the button, then the button.
As the run proceeds, status messages appear in the Control Panel. When
the calculations are complete, the message Results Available appears in the
status bar at the right corner of the Aspen Plus main window.
You select Results from the pulldown list between the and buttons
and use the button to navigate to the next form with results.
1 After reviewing the Material sheet, and selecting Results from the
pulldown list, click the button.
The Poly. Curves sheet appears. You can review the polymer structural
property distribution results for the streams.
2 Click the button.
The Stream FEEDA Results form appears.
3 Use the button to review the remainder of the results.
Input Summary
The input summary for this example is shown here:
SYS-OPTIONS TRACE=YES
RUN-CONTROL MAX-TIME=1000.0
DESCRIPTION "
This example describes how to put together a polymer
simulation model.
"
COMPONENTS
AIBN C8H12N4 /
EB C8H10-4 /
STY C8H8 /
ACN C3H3N /
XYLENE C8H10-3 /
STYSEG C8H8-R /
ACNSEG C3H3N-R /
SAN SAN
FLOWSHEET
BLOCK B1 IN=FEEDA OUT=GRADEA
BLOCK B2 IN=FEEDB OUT=GRADEB
BLOCK M IN=GRADEA GRADEB OUT=PRODUCT
PROPERTIES POLYNRTL
PROP-DATA DATA1
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=atm TEMPERATURE=C PDROP='N/sqm'
PROP-LIST MW
PVAL AIBN 164.0
PVAL SAN 104.150
POLYMERS
SEGMENTS STYSEG REPEAT / ACNSEG REPEAT
POLYMERS SAN
ATTRIBUTES SAN SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM &
FMOM SMOM LDPN LZMOM LFMOM LSFLOW LSFRAC LEFLOW &
LEFRAC LPFRAC
DISTRIBUTION SAN CHAIN-SIZE
STREAM FEEDA
SUBSTREAM MIXED TEMP=70.0 PRES=2.0 MASS-FLOW=24000.0
STREAM FEEDB
SUBSTREAM MIXED TEMP=70.0 PRES=2.0 MASS-FLOW=24000.0
MASS-FRAC AIBN .0040 / STY .10 / ACN .40 / XYLENE &
.4960
BLOCK M MIXER
PARAM PRES=1.0 NPHASE=1 PHASE=L
BLOCK B1 RCSTR
PARAM VOL=10. TEMP=70.0 PRES=2.0 NPHASE=1 PHASE=L &
SCALING=COMPONENTS ALGORITHM=MIXED
CONVERGENCE SOLVER=NEWTON STAB-STRAT=LINE-SEARCH
BLOCK-OPTION TERM-LEVEL=7
REACTIONS RXN-IDS=R1
BLOCK B2 RCSTR
PARAM VOL=10. TEMP=70.0 PRES=2.0 NPHASE=1 PHASE=L &
SCALING=COMPONENTS ALGORITHM=MIXED
CONVERGENCE SOLVER=NEWTON STAB-STRAT=LINE-SEARCH
BLOCK-OPTION TERM-LEVEL=7
REACTIONS RXN-IDS=R2
EO-CONV-OPTI
CONV-OPTIONS
PARAM CHECKSEQ=NO
REACTIONS R1 FREE-RAD
DESCRIPTION "Grade A SAN Polymerization Kinetics"
PARAM
SPECIES INITIATOR=AIBN MONOMER=STY ACN CHAINTAG=EB &
SOLVENT=XYLENE POLYMER=SAN
MON-RSEG STY STYSEG / ACN ACNSEG
INIT-DEC AIBN .00003710 0.0 0.0
CHAIN-INI STY 4820.0 0.0 0.0
CHAIN-INI ACN 225.0 0.0 0.0
PROPAGATION STY STY 4820. 0.0 0.0
PROPAGATION STY ACN 10277. 0.0 0.0
PROPAGATION ACN STY 7165.6 0.0 0.0
PROPAGATION ACN ACN 225. 0.0 0.0
CHAT-MON STY STY 0.289 0.0 0.0
CHAT-MON STY ACN 0.289 0.0 0.0
CHAT-MON ACN STY .0060 0.0 0.0
CHAT-MON ACN ACN .0060 0.0 0.0
TERM-COMB STY STY 13900000. 0.0 0.0
TERM-COMB STY ACN 358000000. 0.0 0.0
TERM-COMB ACN STY 358000000. 0.0 0.0
TERM-COMB ACN ACN 10200000. 0.0 0.0
GEL-EFFECT TERMINATION 2 MAX-PARAMS=10 GE-PARAMS=1.0 0.0 &
2.570 -.005050 9.560 -.01760 -3.030 .007850 0.0 2.0
REACTIONS R2 FREE-RAD
Plotting Distributions
Follow these instructions to:
Display and plot the distribution data for a polymerization reactor
Display a distribution table for a stream or an entire flowsheet
To display and plot the distribution data for a polymerization reactor:
1 From the Blocks folder B1 sub-folder in the Menu Tree, click the Results
form.
2 Click the Distributions tab.
3 From the Aspen Plus menu bar, click Plot, then Plot Wizard.
The Aspen Plus Plot Wizard appears.
4 Click Next, then the icon for Chain Size Distr. plot type.
5 Click Next. If you want to review further plot options, click Next, or
proceed directly to the next step.
6 Click Finish.
The plot is displayed.
Note: You can modify the plot by using the right mouse button to click the
objects in the plot window. From the menu, click Edit, Properties, or Modify
to make necessary changes.
To display a distribution table for a stream or for the entire flowsheet:
1 From the Menu Tree, double-click the Results Summary folder.
2 Click the Streams form.
3 Click the Poly. Curves tab. You may need to use the scroll arrows to
locate the tab.
Note: To activate this feature when the plot window is active, from the Aspen
Plus menu bar, click Edit, then Live Plot.
9 Click Finish.
The plot appears.
You are now ready to make another simulation and see the plot change.
To change the specifications for B1:
1 Keeping the plot window on your screen, from the Blocks folder B1 sub-
folder in the Menu Tree, click the Setup form.
2 On the Specifications sheet, change temperature to 5C and reactor
volume to 70 cum.
3 To rerun the simulation, click the button.
Watch the distribution curve and average properties change in the plot
window.
The repeat segment considered in this example is not available in the Aspen
Polymers databanks. Therefore, the properties of the segment are estimated
using Van Krevelen group contribution methods, by providing the groups
constituting this segment.
Property tables and property sets are used to tabulate liquid heat capacity,
heat capacity and density of the polymer.
PPO-A OLIGOMER
PPO-40 OLIGOMER POLY(PHENYLENE-OXIDE) PPO
PPO-REP SEGMENT POLY(PHENYLENE-OXIDE) PPO
* O
Note: In Aspen Polymers, the Liquid identifier is also used to refer to all
types of solid polymer in the mixed substream. The property models switch
equations (for liquid, glassy, crystalline, or amorphous polymer) depending on
the temperature and the predicted transition temperatures.
9 Click the Tabulate tab, and click PS-1, then click the button to move
it to Selected Prop-Sets.
Your simulation is now complete.
2750
2500
MASS-HEAT-CA J/kg-K
2000 2250
LIQUID CP PPO-A
1750
LIQUID CP PPO-40
1500
Thermal Conductivity
0.18 0.19 0.2 0.21 0.22 0.23 0.24 0.25 0.26 0.27 0.28 0.29 0.3
Semi-Crystalline
THERMAL-COND Watt/m-K
LIQUID K PPO-A
LIQUID K PPO-40
Liquid
Amorphous Solid
Glassy Solid
1100
LIQUID RHO PPO-A
1050
LIQUID RHO PPO-40
1000
MASS-DENSITY kg/cum
900 950850
Input Summary
The input language summary for this example is shown here:
TGS
DYNAPLUS
DPLUS RESULTS=ON
DESCRIPTION "
Estimate properties of 2,6-DM PPO using Van Krevelen
group contribution methods.
"
COMPONENTS
PPO-A PPO /
PPO-40 PPO /
PPO-REP
PROPERTIES POLYNRTL
STRUCTURES
VANKREV PPO-REP 118 1 / 149 1
PROP-DATA PURE-1
IN-UNITS SI PRESSURE=atm TEMPERATURE=C PDROP='N/sqm'
POLYMERS
SEGMENTS PPO-REP REPEAT
OLIGOMERS PPO-A PPO-REP 100. / PPO-40 PPO-REP 100.
References
Van Krevelen, D. W. (1990). Properties of Polymers, 3rd Ed. Amsterdam:
Elsevier.
This example demonstrates how to use the data regression (DRS) capabilities
to fit the mixture parameters of an equation of state (EOS) model to binary
vapor-liquid equilibrium (VLE) data.
The steps covered include:
Defining the Simulation
Creating a New Run
Specifying Setup and Global Options
Specifying and Characterizing Components
Specifying Physical Property Method
Entering Experimental Data
Specifying a Regression Case
Specifying Physical Property Parameters
Running the Simulation and Examining the Results
Correlative models used to describe thermodynamic properties of mixtures
often contain binary interaction parameters. These parameters account for
mixture non-idealities, and are necessary for accurate representation of the
mixture behavior. For each constituent pair of a multicomponent mixture,
these parameters are obtained by regressing some form of binary
experimental information.
In this example, binary VLE data of ethylene-polyethylene mixture is
regressed to obtain the two binary interaction parameters of the Sanchez-
Lacombe EOS, k ij and ij . You can find the details of this model in Volume 2
of the Aspen Polymers User Guide. Ethylene-polyethylene binary mixture is
encountered in polyolefin production, and at high pressures the
thermodynamic behavior of this mixture can be described by an EOS such as
the Sanchez-Lacombe model.
Defining Unit-Sets
To define a unit-set:
1 From the Setup folder in the Menu Tree, click the Units-Sets sub-folder.
2 Click New.
3 To accept US-1 as the unit set ID, click OK.
A dialog box appears requesting approval to make US-1 the global unit set
4 Click No.
5 On the Standard sheet, from the Copy from pulldown list, click Eng.
6 Set the following options:
o Temperature = K
o Pressure = bar
7 Click the Transport tab
8 From the Density pulldown list, click kg/cum.
9 Repeat these steps to create unit set US-3 (step 3) in which:
o Copy from=SI (step 5)
o Pressure=bar (step 6)
o Delta P=bar (step 6)
10 Click .
A Component Specifications Selection sheet appears.
PE Oligomer POLY(ETHYLENE) PE
ETHYLENE Conventional ETHYLENE C2H4
C2H4-R Segment ETHYLENE-R C2H4-R
2 From the Components folder in the Menu Tree, click the Polymers sub-
folder.
3 On the Segments sheet, from the Type pulldown list, click REPEAT.
4 Click the Oligomers tab, select C2H4-R as the segment and enter 1132
for PE in the oligomer row.
5 Click .
The Properties Specifications Global sheet appears.
STD-DEV 0 1% .001 0
DATA 399.15 455.8 .0018 0.9982
DATA 399.15 790.30 .0037 0.9963
DATA 399.15 1135 .0055 0.9945
DATA 399.15 1479 .0075 0.9925
DATA 399.15 1824 .0107 0.9893
DATA 399.15 2168 .0136 0.9864
DATA 399.15 2513 .0158 0.9842
DATA 399.15 2857 .0175 0.9825
DATA 399.15 3202 .0198 0.9802
DATA 399.15 3546 .0221 0.9779
DATA 399.15 3891 .0242 0.9758
DATA 399.15 4235 .0255 0.9745
DATA 399.15 4580 .0285 0.9715
DATA 399.15 4924 .0305 0.9695
DATA 399.15 5269 .0330 0.967
DATA 399.15 5613 .0359 0.9641
After the first four entries, your data sheet should look like this:
STD-DEV 0 1% .001 0
DATA 413.15 455.8 .0015 0.9985
DATA 413.15 790.30 .0034 0.9966
DATA 413.15 1135 .0048 0.9952
DATA 413.15 1479 .0068 0.9932
DATA 413.15 1824 .0087 0.9913
DATA 413.15 2168 .0112 0.9888
DATA 413.15 2513 .0131 0.9869
DATA 413.15 2857 .0151 0.9849
DATA 413.15 3202 .0166 0.9834
DATA 413.15 3546 .0189 0.9811
DATA 413.15 3891 .0208 0.9792
DATA 413.15 4235 .0235 0.9765
DATA 413.15 4580 .0250 0.975
DATA 413.15 4924 .0277 0.9723
DATA 413.15 5269 .0296 0.9704
DATA 413.15 5613 .0328 0.9672
DATA 428.15 455.8 .0013 0.9987
DATA 428.15 790.30 .0029 0.9971
DATA 428.15 1135 .0039 0.9961
DATA 428.15 1479 .0055 0.9945
DATA 428.15 1824 .0074 0.9926
DATA 428.15 2168 .0090 0.991
DATA 428.15 2513 .0105 0.9895
DATA 428.15 2857 .012 0.988
DATA 428.15 3202 .0146 0.9854
DATA 428.15 3546 .0164 0.9836
DATA 428.15 3891 .0178 0.9822
IN-UNITS SI
DESCRIPTION "
The objective of this example is to demonstrate how to
use the data regression capabilities to fit the mixture
parameters of an EOS model to binary VLE data. "
COMPONENTS
ETHYLENE C2H4 ETHYLENE /
PE PE PE /
C2H4-R C2H4-R C2H4-R
PROPERTIES POLYSL
PROP-DATA MW
IN-UNITS SI PRESSURE=bar PDROP=bar
PROP-LIST MW
PVAL PE 31756.8563
PROP-DATA PCES-1
IN-UNITS ENG DENSITY='kg/cum' PRESSURE=bar TEMPERATURE=K &
PDROP=psi
PROP-LIST SLPSTR / SLRSTR / SLTSTR
PVAL C2H4-R 4250 / 887 / 673
PVAL ETHYLENE 3339 / 660 / 291
PROP-DATA SLETIJ-1
IN-UNITS SI
PROP-LIST SLETIJ
BPVAL PE ETHYLENE -.4040175320
BPVAL ETHYLENE PE -.4040175320
PROP-DATA SLKIJ-1
IN-UNITS SI
PROP-LIST SLKIJ
BPVAL PE ETHYLENE -.0900000000
BPVAL ETHYLENE PE -.0900000000
PARAMETERS
BIPARAMETER 1 SLKIJ PE ETHYLENE 1 1.00000000E-02 &
-9.00000000E-02 1.10000000E-01 1.00000000E+00
BIPARAMETER 2 SLETIJ PE ETHYLENE 1 -9.00000000E-02 &
-9.90000000E-01 8.10000000E-01 1.00000000E+00
CASE R-1
DATA-GROUP C2PE399
IN-UNITS SI PRESSURE=kPa PDROP='N/sqm'
SYSTEM-DEF TPXY ETHYLENE PE COMPOSITION=MASS-FRAC
PHASE-EQ VL ETHYLENE
DATA 1 399.15 455.8 .0018 1 /
2 399.15 790.30 .0037 1 /
3 399.15 1135 .0055 1 /
4 399.15 1479 .0075 1 /
5 399.15 1824 .0107 1 /
6 399.15 2168 .0136 1 /
7 399.15 2513 .0158 1 /
8 399.15 2857 .0175 1 /
9 399.15 3202 .0198 1 /
10 399.15 3546 .0221 1 /
11 399.15 3891 .0242 1 /
12 399.15 4235 .0255 1 /
13 399.15 4580 .0285 1 /
14 399.15 4924 .0305 1 /
15 399.15 5269 .0330 1 /
16 399.15 5613 .0359 1
STD-DEV 1 0 -1 .001 -1
DATA-GROUP C2PE413
IN-UNITS SI PRESSURE=kPa PDROP='N/sqm'
SYSTEM-DEF TPXY ETHYLENE PE COMPOSITION=MASS-FRAC
PHASE-EQ VL ETHYLENE
DATA 1 413.15 455.8 .0015 1 /
2 413.15 790.30 .0034 1 /
3 413.15 1135 .0048 1 /
4 413.15 1479 .0068 1 /
5 413.15 1824 .0087 1 /
6 413.15 2168 .0112 1 /
7 413.15 2513 .0131 1 /
8 413.15 2857 .0151 1 /
9 413.15 3202 .0166 1 /
10 413.15 3546 .0189 1 /
11 413.15 3891 .0208 1 /
12 413.15 4235 .0235 1 /
13 413.15 4580 .0250 1 /
14 413.15 4924 .0277 1 /
15 413.15 5269 .0296 1 /
16 413.15 5613 .0328 1
STD-DEV 1 0 -1 .001 -1
DATA-GROUP C2PE428
IN-UNITS SI PRESSURE=kPa PDROP='N/sqm'
SYSTEM-DEF TPXY ETHYLENE PE COMPOSITION=MASS-FRAC
PHASE-EQ VL ETHYLENE
DATA 1 428.15 455.8 .0013 1 /
2 428.15 790.30 .0029 1 /
3 428.15 1135 .0039 1 /
4 428.15 1479 .0055 1 /
5 428.15 1824 .0074 1 /
POLYMERS
SEGMENTS C2H4-R REPEAT
OLIGOMERS PE C2H4-R 1132
References
Hao, W., Elbro, H. S., & Alessi, P. (1992). Polymer Solution Data Collection.
Chemistry Data Series, XIV, Part 1, DECHEMA.
The feed stream consists of 100 kg/hr of vinyl acetate, 80 kg/hr of methanol,
0.0258 kg/hr of initiator at one atmosphere and 64.7C. The initial charge is
determined by multiplying the flowrates by the cycle time needed to prepare
the charge.
Defining Unit-Sets
To define a unit-set:
1 From the Setup folder in the Menu Tree, click the Units-Sets sub-folder.
2 Click New.
3 To accept US-1 as the unit set ID, click OK.
A dialog box appears requesting approval to make US-1 the global unit set
4 Click Yes.
5 On the Standard sheet, from the Copy from pulldown list, click SI.
6 Set the following options:
Mass Flow = kg/hr
Mole Flow = kmol/hr
Pressure = atm
2 Click .
3 On the Segments sheet, from the Type pulldown list, click REPEAT to
define VAC-SEG as a repeat segment.
4 Click the Polymers tab, from the Built-in attribute group pulldown list,
click Free-radical selection.
5 To modify the attribute group, click Edit.
6 On the Polymer Attributes dialog box, click to clear the attributes that
do not appear below.
Defining a Prop-Set
To enter a Prop-Set:
1 From the Properties folder in the Menu Tree, click the Prop-Sets sub-
folder.
2 In the Object Manager, click New.
3 Enter DENSITY as the new name of the property set. Click OK.
4 On the Properties tab, click in the Physical properties cell to select
RHO, click the Units cell to select kg/cum. Repeat this step to select
RHOMX, with units of kg/cum.
5 Click the Qualifiers tab from the Phase pulldown list, click and Liquid.
6 From the pulldown lists in the Component row, click to select VAC,
POLY, and MEOH.
Specifying Gel-Effect
To specify gel effect:
1 Click the Options tab, and click the Gel Effect option.
2 Click the Gel Effect tab.
3 From the No. pulldown list, click New.
Note: To enter data in the No. field you can also use the right mouse button
to click in the field, then click Create from the menu.
4 Enter the information for the gel effect correlation as shown:
1 0.0 1.554E-4
2 0.5 1.551E-4
3 1.0 1.546E-4
4 1.5 1.542E-4
5 2.0 1.538E-4
6 2.5 1.532E-4
7 3.0 1.526E-4
8 3.5 1.522E-4
9 4.0 1.518E-4
10 4.5 1.513E-4
11 5.0 1.509E-4
6 Click the Initial Conditions tab, enter:
o Temperature = 65C
o Pressure = 1 atm
The Reactor type and Block name fields already display RBATCH and
BATCH respectively.
To create the second data set:
1 From the Menu Tree Model Analysis Tools folder Data Fit sub-folder,
click the Data Set sub-folder.
2 In the Object manager, click New.
3 In the Enter ID field, enter D-PROPAG, and in the Select Type field,
click to select PROFILE-DATA. Click OK.
4 On the Define sheet, enter:
o Variable name = DUTY
o Variable = DUTY
o Model and block name = BATCH
RBATCH is already selected as the Block name.
5 Click the Data tab, enter a Std-Dev of 5.0.
Enter the time and heat release data shown:
Data Time (hr) Instantaneous Heat Release
Point (Watt)
2 0.5 -3750.0
3 1.0 -3700.0
4 1.5 -3300.0
5 2.0 -3250.0
6 2.5 -3100.0
7 3.0 -3000.0
8 3.5 -2800.0
9 4.0 -2650.0
DESCRIPTION "
This example illustrates how to fit kinetic rate constant
parameters in Aspen Polymers. "
COMPONENTS
MPO C4H6O2-1 /
VAC C4H6O2-1 /
MEOH CH4O /
POLY PVAC /
VAC-SEG C4H6O2-R-3
FLOWSHEET
BLOCK BATCH IN=FEED OUT=PRODUCT
PROPERTIES POLYNRTL
PROP-DATA DATA1
IN-UNITS SI
PROP-LIST MW
POLYMERS
SEGMENTS VAC-SEG REPEAT
POLYMERS POLY
ATTRIBUTES POLY SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM &
FMOM SMOM LSFLOW LEFLOW
STREAM FEED
SUBSTREAM MIXED TEMP=64.70 PRES=1.0
MASS-FLOW MPO .028 / VAC 100 / MEOH 80
EO-CONV-OPTI
SEQUENCE S-1 BATCH INITIA BATCH (RETURN INITIA) PROPAGA BATCH &
(RETURN PROPAGA)
REGRESSION INITIA
DATA D-INITIA
VARY REACT-VAR REACTION=REAC-1 VARIABLE=IDPRE-EXP &
REGRESSION PROPAGA
DATA D-PROPAG
VARY REACT-VAR REACTION=REAC-1 VARIABLE=PRPRE-EXP &
SENTENCE=PROPAGATION ID1=VAC ID2=VAC
LIMITS 1E3 5E4
ALGORITHM AFCTOL=1.0 INIRIN=0
PROFILE-DATA D-INITIA
PARAM BLOCK=BATCH TEMP=65.0 PRES=1 UNITS=hr
DEFINE ZZTEMP BLOCK-VAR BLOCK=BATCH SENTENCE=PARAM &
VARIABLE=TEMP
DEFINE ZZPRES BLOCK-VAR BLOCK=BATCH SENTENCE=PARAM &
VARIABLE=PRES
VECTOR-DEF MFINIT BLOCK-VEC BLOCK=BATCH SENTENCE= &
REGR-C-PROF VARIABLE=MASSFRAC-L ID2=MPO
PROFILE-DATA D-PROPAG
A5 Fractionating Oligomers 69
Creating a New Run
To start Aspen Polymers:
1 Start Aspen Plus from the Start Menu or by double clicking the Aspen Plus
icon on your desktop.
The Aspen Plus Startup dialog box appears.
2 On the Aspen Plus Startup dialog box, click the Template option. Click
OK.
The New dialog box appears. You use this dialog box to specify the
simulation template and Run Type for the new run. Aspen Plus uses the
Simulation Template you choose to automatically set various defaults
appropriate to your application.
3 For this example, click Polymers with English Units for your template.
The default Run type, Flowsheet, is appropriate for this example.
4 Click OK.
The Aspen Plus main window is now active.
70 A5 Fractionating Oligomers
2 Move the mouse to the Process Flowsheet window, then click to place
the block (B1).
To place streams:
1 Click the Material STREAMS icon on the left corner of the Model Library
and move the mouse to the Process Flowsheet window.
Red and blue arrows appear on the unit operation blocks. These arrows
indicate the location of the required (red) and optional (blue) stream
connection ports.
2 Place the + cursor on the red feed arrow and click to make the
connection. Move the mouse away from the block and click again to place
stream 1.
3 Repeat step 2 to create the other streams.
4 After placing all streams, turn off the insert mode by clicking the Select
Mode button in the upper left corner of the Model Library.
Renaming Streams
To rename streams:
1 Click the stream to be renamed, such as 5, and click the right mouse
button to display the menu.
2 Click Rename Stream.
3 Enter the Stream ID as 1P and click OK.
You are now ready to enter input data for your simulation.
A5 Fractionating Oligomers 71
This example illustrates the use of Aspen Polymers for modeling
polymer/oligomer fractionation process.
Defining Unit-Sets
To define a unit-set:
1 From the Setup folder in the Menu Tree, click the Units-Sets sub-folder.
2 Click New.
3 In the Enter ID field, enter SET1, then click OK.
A dialog box appears requesting approval to make SET1 the global unit set
4 Click No.
5 On the Standard sheet, from the Copy from pulldown list, click Eng.
6 Set the following options:
o Pressure = atm
o Temperature = K
72 A5 Fractionating Oligomers
Note: Styrene segment is needed to define the structure of oligomers.
To characterize components:
1 From the Menu Tree Components folder, double-click the Polymers sub-
folder.
2 In the Polymers sub-folder, click Characterization.
3 On the Segments sheet, in the Type field for segment STY-SEG, click
REPEAT.
4 Click the Oligomers tab.
You need to provide the number of segments for each oligomer. Enter
these as shown below using even numbers from 2-20.
A5 Fractionating Oligomers 73
Defining Pure Component Parameters
You enter three pure component properties in this example. To enter the first
pure component constant:
1 From the Menu Tree Properties folder, double-click the Parameters sub-
folder.
2 From the Parameters sub-folder, click the Pure Component sub-folder.
3 In the Object manager, click New.
4 Click to select T-dependent Correlation and expand Liquid vapor
pressure to select PLXANT-1. Click OK.
5 On the Input sheet, enter data shown in the Input Summary PROP-DATA
PLXANT-1 paragraph on page 80. A partial example is shown here:
74 A5 Fractionating Oligomers
5 From the unit field pulldown list on menu bar, click SET1.
A5 Fractionating Oligomers 75
Defining Properties to be Calculated
To select the properties to be calculated:
1 From the Menu Tree Properties folder, click the Prop-Sets sub-folder.
2 In the Object manager, click New.
3 To accept PS-1 as the new name, Click OK.
4 On the Properties sheet, in the Physical properties field, click GAMMA.
5 Click the Qualifiers tab, from the Phase pulldown list, click Liquid.
76 A5 Fractionating Oligomers
Specifying Block Conditions
To specify the Flash3 block:
1 From the Menu Tree, double-click the Blocks folder.
2 From the Blocks folder, doubl-click the B1 sub-folder.
3 From the B1 sub-folder, click Input.
4 On the Specifications sheet, enter:
o Temperature = 101.48 F
o Pressure = 14.7 psi
5 Click the Key Components tab, from the Key component in the 2nd
liquid phase pulldown list, click SOLVENT2.
A5 Fractionating Oligomers 77
As the run proceeds, status messages appear in the Control Panel. When
the calculations are complete, the message Results Available appears in the
status bar at the right corner of the Aspen Plus main window.
To examine the results:
1 From the Menu Tree, double-click the Results Summary folder.
The Summary sheet appears.
2 From the pulldown list between the and buttons, click Results.
3 Click to navigate to the next form with results.
After flashing, the quasi-ternary mixture separates into two liquid phases, a
solvent rich phase and an antisolvent rich phase. Using the Flory-Huggins
model, the equilibria compositions in solvent and antisolvent phases are
computed. The Flory-Huggins predictions for liquid-liquid equilibria of
polystyrene-Benzene-Ethanol system are shown here:
78 A5 Fractionating Oligomers
The solid line is the mass distribution of polystyrene oligomer for the feed
stream 1P, dashed lines are for the product streams.
When ethanol is added to the polystyrene/benzene mixture, the majority of
high molecular weight PS remains in the solvent rich stream while most of the
low molecular weight ones dissolve into the antisolvent stream.
Input Summary
The input language summary for this example is shown here:
IN-UNITS ENG
DESCRIPTION "
This example illustrates the use of Aspen Polymers
for modeling polymer/oligomer fractionation process. "
COMPONENTS
SOLVENT1 C6H6 /
SOLVENT2 C2H6O-2 /
STY-SEG C8H8-R /
OL2
OL3
OL4
OL5
OL6
OL7
OL8
OL9
OL10
OL11
FLOWSHEET
BLOCK B1 IN=1 1P OUT=2 3 4
PROPERTIES POLYFH
PROP-DATA U-3
IN-UNITS ENG PRESSURE=atm TEMPERATURE=K PDROP=psi
PROP-LIST TC / PC / ZC
PVAL OL2 500.0 / 20.0 / .20
PVAL OL3 500.0 / 20.0 / .20
PVAL OL4 500.0 / 20.0 / .20
PVAL OL5 600.0 / 20.0 / .20
PVAL OL6 700.0 / 20.0 / .20
PVAL OL7 800.0 / 20.0 / .20
PVAL OL8 850.0 / 20.0 / .20
A5 Fractionating Oligomers 79
PVAL OL9 900.0 / 20.0 / .20
PVAL OL10 1000.0 / 20.0 / .20
PVAL OL11 1100.0 / 20.0 / .20
PROP-DATA PLXANT-1
IN-UNITS ENG
PROP-LIST PLXANT
PVAL OL2 -10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PVAL OL3 -10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PVAL OL4 -10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PVAL OL5 -10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PVAL OL6 -10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PVAL OL7 -10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PVAL OL8 -10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PVAL OL9 -10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PVAL OL10 -10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PVAL OL11 -10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PROP-DATA FHCHI-1
IN-UNITS ENG
PROP-LIST FHCHI
BPVAL SOLVENT1 SOLVENT2 1.7390 0.0
BPVAL SOLVENT1 OL2 .2210 0.0
BPVAL SOLVENT1 OL3 .2210 0.0
BPVAL SOLVENT1 OL4 .2210 0.0
BPVAL SOLVENT1 OL5 .2210 0.0
BPVAL SOLVENT1 OL6 .2210 0.0
BPVAL SOLVENT1 OL7 .2210 0.0
BPVAL SOLVENT1 OL8 .2210 0.0
BPVAL SOLVENT1 OL9 .2210 0.0
BPVAL SOLVENT1 OL10 .2210 0.0
BPVAL SOLVENT1 OL11 .2210 0.0
BPVAL SOLVENT2 OL2 1.4510 0.0
BPVAL SOLVENT2 OL3 1.4510 0.0
BPVAL SOLVENT2 OL4 1.4510 0.0
BPVAL SOLVENT2 OL5 1.4510 0.0
BPVAL SOLVENT2 OL6 1.4510 0.0
BPVAL SOLVENT2 OL7 1.4510 0.0
BPVAL SOLVENT2 OL8 1.4510 0.0
BPVAL SOLVENT2 OL9 1.4510 0.0
BPVAL SOLVENT2 OL10 1.4510 0.0
BPVAL SOLVENT2 OL11 1.4510 0.0
POLYMERS
SEGMENTS STY-SEG REPEAT
OLIGOMERS OL2 STY-SEG 2 / OL3 STY-SEG 4 / OL4 STY-SEG &
6 / OL5 STY-SEG 8 / OL6 STY-SEG 10 / OL7 STY-SEG &
12 / OL8 STY-SEG 14 / OL9 STY-SEG 16 / OL10 &
STY-SEG 18 / OL11 STY-SEG 20
STREAM 1
SUBSTREAM MIXED TEMP=101.480003 PRES=14.6959488 &
MASS-FLOW=2.10
MASS-FLOW SOLVENT2 2.10
80 A5 Fractionating Oligomers
STREAM 1P
SUBSTREAM MIXED TEMP=101.480 PRES=14.6959488 MASS-FLOW=1.0
MASS-FRAC SOLVENT1 .708990 / OL2 .020 / OL3 .060 / OL4 &
.080 / OL5 .060 / OL6 .040 / OL7 .020 / OL8 .010 &
/ OL9 .0010 / OL10 .00010 / OL11 .000010
BLOCK B1 FLASH3
PARAM TEMP=101.480003 PRES=14.6959488 L2-COMP=SOLVENT2
References
Wu, A. H., & Prausnitz, J. M. (1990). J. Appl. Pol. Sci., 39, pp. 629-637.
A5 Fractionating Oligomers 81
A6 Calculating End-Use
Properties
Defining a Prop-Set
To be able to list desired polymer properties, you need to create property sets
and property tables, to do this:
To enter a Prop-Set:
1 From the Properties folder in the Menu Tree, click the Prop-Sets sub-
folder.
2 In the Object Manager, click New.
3 In the Enter ID field, enter IV as the new name of the property set. Click
OK.
4 On the Properties sheet, click in the Physical properties cell to select
IV.
5 Repeat steps 1-4 to create Prop-Sets named MI-K, MI-S and ZVIS (step
3). In step 4, select the corresponding user property name from the
Physical properties cell.
6 Click .
A Required Properties Input Complete dialog box appears.
11 Click the Tabulate tab, in Column numbers 1-4, enter IV, MI-K, MI-S,
and ZVIS respectively.
12 Click the Table Format button, enter:
Alternately, you can copy the files USERFORT.DLL and USERFORT.DLOPT from
the Aspen Polymers application folder (app) to your working directory. Then,
from the Aspen Plus menu bar, chose Run | Settings | Engine Files and
enter USERFORT.DLOPT in the Linker Options box before running the
simulation.
Follow these instructions to:
Run the simulation
Examine the results
To run the simulation:
1 From the Aspen Plus toolbar, click .
The Control Panel appears.
2 To run the simulation, click .
As the run proceeds, status messages appear in the Control Panel. When
the calculations are complete, the message Results Available appears in the
status bar at the right corner of the Aspen Plus main window.
To examine the results:
1 From the Menu Tree, double-click the Results Summary folder.
The Summary sheet appears.
2 From the pulldown list between the and buttons, click Results.
3 Click to navigate to the next form with results.
Input Summary
The input language summary for this example is shown here:
IN-UNITS MET
COMPONENTS
HDPE PE /
C2H4-R C2H4-R
FLOWSHEET
BLOCK B1 IN=1 OUT=2
PROPERTIES POLYFH
POLYMERS
SEGMENTS C2H4-R REPEAT
POLYMERS HDPE
ATTRIBUTES HDPE SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM &
FMOM SMOM
USER-PROPERT IV SUBROUTINE=USRPRP
PROP-SET IV IV SUBSTREAM=MIXED
STREAM 1
SUBSTREAM MIXED TEMP=420. <K> PRES=1.
MASS-FLOW HDPE 1000.
COMP-ATTR HDPE SFRAC ( 1. )
COMP-ATTR HDPE DPN ( 1000. )
COMP-ATTR HDPE PDI ( 3.61 )
BLOCK B1 HEATER
PARAM TEMP=420. <K> PRES=1.
EO-CONV-OPTI
SENSITIVITY END-USE
DEFINE IV STREAM-PROP STREAM=2 PROPERTY=IV
DEFINE MIK STREAM-PROP STREAM=2 PROPERTY=MI-K
DEFINE MIS STREAM-PROP STREAM=2 PROPERTY=MI-S
DEFINE ZVIS STREAM-PROP STREAM=2 PROPERTY=ZVIS
TABULATE 1 "IV" COL-LABEL="IV"
TABULATE 2 "MI-K" COL-LABEL="MI-K"
References
Karol, F. J., Brown, G. L., & Davison, J. M. (1973). Chromocene-Based
Catalysts for Ethylene Polymerization: Kinetic Parameters. J. of Polymer
Science: Polymer Chemistry Edition, 11, pp. 413-424.
Sinclair, K. B. (1983). Characteristics of Linear LPPE and Description of UCC
Gas Phase Process. Process Economics Report. Menlo Park, CA: SRI
International.
Process Definition
Styrene is polymerized in a reactor train consisting of two CSTRs followed by
a plug flow reactor. All of the reactors are considered liquid filled, and are
therefore modeled without taking into account vapor-liquid equilibrium.
Unreacted monomer is flashed in a devolatilizer to be modeled as an ideal
flash unit.
Process Conditions
The first CSTR operates at 120C, 1 atm, and the second operates at 160C, 1
atm. Both have a volume of 20 m3. The plug flow reactor operates at 1 atm,
with a temperature range of 160-200C from the inlet to the outlet.
The process conditions are:
Components Name Databank Description
Feeds
Temperature (C) 25
Pressure (atm) 1
Mass Flow (kg/hr) 7000
CSTR-1 120 1 20 m3
CSTR-2 160 1 20 m3
PLUG 160-200 1 80 m length by
0.40 m diam
DV-H1 220 1
FLASH-1 220 1
FLASH-2 220 1
Reactors / Kinetics
The kinetics of bulk and solution polymerization of styrene have been studied
extensively. It has been reported (Albright, 1985) that chain transfer to the
monomer and the Diels-Alder dimers primarily controls the molecular weight
of polystyrene. Furthermore, based on chemical evidence and kinetic
modeling, it has been reported that termination of the growing chains occurs
principally by combination over disproportionation.
The reactions included from the built-in free-radical kinetics are:
Description Reaction
Initiator decomposition I R
Thermal initiation M C P1
Process Studies
The model is used to study the effect of feed flow rate on styrene conversion,
polymer properties, and recycle styrene flowrate. In order to determine the
effect of feed flowrate on styrene conversion, polymer properties, and recycle
flowrate, a sensitivity study is carried out with feed mass flow as the varied
parameter.
Input Summary
An input language summary for this example is shown here:
; Free radical bulk polymerization of styrene by thermal and chemical
; initiation
; QSSA case
;
DYNAPLUS
DPLUS RESULTS=ON
TITLE &
'Bulk Polymerization of Styrene by Thermal & Chemical Initiation'
SYS-OPTIONS TRACE=YES
RUN-CONTROL MAX-TIME=2000.0
COMPONENTS
TBP C8H18O2 /
CINI C8H8 /
STY C8H8 /
PS PS-1 /
STY-SEG C8H8-R /
EB C8H10-4 /
DDM C12H26S /
H2O H2O
FLOWSHEET
BLOCK CSTR-1 IN=FEED OUT=R1-P
BLOCK CSTR-2 IN=R1-P OUT=R2-P
BLOCK PLUG IN=R2-P OUT=R3-P
BLOCK FLASH-1 IN=R3-PH2 OUT=REC-STY R3-OLI
BLOCK FLASH-2 IN=R3-OLI STRIP-AG OUT=VAPOR PRODUCT
BLOCK DV-H1 IN=R3-P OUT=R3-PH2
PROPERTIES POLYNRTL
PROP-DATA MW
IN-UNITS MET ENTHALPY-FLO='kcal/hr' PRESSURE=bar &
TEMPERATURE=C DELTA-T=C MOLE-ENTHALP='kcal/kmol' &
MASS-ENTHALP='kcal/kg' HEAT=kcal PDROP=bar &
INVERSE-PRES='1/bar'
PROP-LIST MW
PVAL TBP 216.320
POLYMERS
SEGMENTS STY-SEG REPEAT
POLYMERS PS
ATTRIBUTES PS SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM FMOM &
SMOM LDPN LZMOM LFMOM LSFLOW LSFRAC LEFLOW LEFRAC &
LPFRAC
DISTRIBUTION PS
STREAM FEED
SUBSTREAM MIXED TEMP=25.00000 PRES=1.013250 MASS-FLOW=7000. &
NPHASE=1 PHASE=L
MASS-FRAC TBP 0.0003 / STY .980 / EB 0.019 / DDM &
0.0007
STREAM STRIP-AG
SUBSTREAM MIXED TEMP=25.00000 PRES=2.026500 &
MASS-FLOW=144.0000 NPHASE=1 PHASE=L
MASS-FLOW H2O 144.0000
EO-CONV-OPTI
SENSITIVITY S1
IN-UNITS SI
DEFINE R3MWN COMP-ATTR-VAR STREAM=R3-P SUBSTREAM=MIXED &
COMPONENT=PS ATTRIBUTE=MWN ELEMENT=1
DEFINE R3MWW COMP-ATTR-VAR STREAM=R3-P SUBSTREAM=MIXED &
COMPONENT=PS ATTRIBUTE=MWW ELEMENT=1
DEFINE R3PDI COMP-ATTR-VAR STREAM=R3-P SUBSTREAM=MIXED &
COMPONENT=PS ATTRIBUTE=PDI ELEMENT=1
DEFINE R2MWN COMP-ATTR-VAR STREAM=R2-P SUBSTREAM=MIXED &
COMPONENT=PS ATTRIBUTE=MWN ELEMENT=1
DEFINE R2MWW COMP-ATTR-VAR STREAM=R2-P SUBSTREAM=MIXED &
COMPONENT=PS ATTRIBUTE=MWW ELEMENT=1
DEFINE R2PDI COMP-ATTR-VAR STREAM=R2-P SUBSTREAM=MIXED &
COMPONENT=PS ATTRIBUTE=PDI ELEMENT=1
DEFINE R1MWN COMP-ATTR-VAR STREAM=R1-P SUBSTREAM=MIXED &
COMPONENT=PS ATTRIBUTE=MWN ELEMENT=1
DEFINE R1MWW COMP-ATTR-VAR STREAM=R1-P SUBSTREAM=MIXED &
COMPONENT=PS ATTRIBUTE=MWW ELEMENT=1
DEFINE R1PDI COMP-ATTR-VAR STREAM=R1-P SUBSTREAM=MIXED &
COMPONENT=PS ATTRIBUTE=PDI ELEMENT=1
DEFINE R3STY MASS-FRAC STREAM=R3-P SUBSTREAM=MIXED &
COMPONENT=STY
DEFINE R3PS MASS-FRAC STREAM=R3-P SUBSTREAM=MIXED &
COMPONENT=PS
DEFINE R2STY MASS-FRAC STREAM=R2-P SUBSTREAM=MIXED &
COMPONENT=STY
DEFINE R2PS MASS-FRAC STREAM=R2-P SUBSTREAM=MIXED &
COMPONENT=PS
DEFINE R1STY MASS-FRAC STREAM=R1-P SUBSTREAM=MIXED &
COMPONENT=STY
DEFINE R1PS MASS-FRAC STREAM=R1-P SUBSTREAM=MIXED &
CONV-OPTIONS
PARAM CHECKSEQ=NO
SEQUENCE S-1 CSTR-1 CSTR-2 PLUG S1 CSTR-1 CSTR-2 PLUG DV-H1 &
FLASH-1 FLASH-2 (RETURN S1)
REACTIONS R1 FREE-RAD
IN-UNITS SI
DESCRIPTION "Free-Radical Kinetic Scheme"
PARAM QSSA=YES
SPECIES INITIATOR=TBP COINITIATOR=CINI MONOMER=STY &
CHAINTAG=EB DDM POLYMER=PS
MON-RSEG STY STY-SEG
INIT-DEC TBP 1.6220E+11 1.1530E+08 0.0 EFFIC=.80 NRADS=2
INIT-SP STY CINI 438000.0 1.1480E+08 0.0
CHAIN-INI STY 1.0510E+07 2.9570E+07 0.0
PROPAGATION STY STY 1.0510E+07 2.9570E+07 0.0
CHAT-MON STY STY 3310000.0 5.3020E+07 0.0
CHAT-AGENT STY EB 1051.0 2.9590E+07 0.0
CHAT-AGENT STY DDM 1051.0 2.9590E+07 0.0
TERM-COMB STY STY 1.2550E+09 7017000.0 0.0
INIT-SP-EFF STY COEFFA=0.0 COEFFB=3.0 COEFFC=0.0
;
1
0.9
Styrene Conversion
0.8
0.7
0.6 CSTR-1
0.5 CSTR-2
0.4 PLUG
0.3
0.2
0.1
0
3000 4000 5000 6000 7000 8000 9000
Feed Flow Rate, kg/hr
The effect of feed flow rate on the number average molecular weight is shown
here:
Number Average MW
180000
CSTR-1
CSTR-2
160000
PLUG
140000
120000
3000 4000 5000 6000 7000 8000 9000
The effect of flow rate on the weight average molecular weight is shown here:
350000
Weight Average MW
300000 CSTR-1
CSTR-2
PLUG
250000
200000
3000 4000 5000 6000 7000 8000 9000
The effect of feed flow rate on the polydispersity index is shown here:
Polydispersity Index
1.9
CSTR-1
CSTR-2
1.8
PLUG
1.7
1.6
3000 4000 5000 6000 7000 8000 9000
Feed Flow Rate, kg/hr
The effect of feed flow rate on the recycle flow rate is shown here:
1500
Recycle Flow Rate, kg/hr
1000
500
0
3000 4000 5000 6000 7000 8000 9000
Feed Flow Rate, kg/hr
SUBSTREAM: MIXED
PHASE: LIQUID LIQUID LIQUID LIQUID LIQUID
COMPONENTS: KG/HR
TBP 2.1000 0.0 1.2516 6.5145-02 0.0
CINI 0.0 0.0 0.0 0.0 0.0
STY 6860.0000 147.1938 5147.1336 2515.9679 485.4496
PS 0.0 5645.3447 1713.7193 4346.0858 5645.3446
EB 133.0000 12.7095 132.9956 132.9818 47.5461
DDM 4.9000 3.9310 4.8998 4.8993 4.4821
H2O 0.0 12.4952 0.0 0.0 0.0
COMPONENTS: MASS FRAC
TBP 3.0000-04 0.0 1.7881-04 9.3065-06 0.0
CINI 0.0 0.0 0.0 0.0 0.0
STY 0.9800 2.5284-02 0.7353 0.3594 7.8516-02
PS 0.0 0.9697 0.2448 0.6209 0.9131
EB 1.9000-02 2.1831-03 1.8999-02 1.8997-02 7.6900-03
DDM 7.0000-04 6.7523-04 6.9998-04 6.9990-04 7.2493-04
H2O 0.0 2.1463-03 0.0 0.0 0.0
TOTAL FLOW:
KMOL/HR 67.1522 56.4492 67.1565 67.1624 59.3342
KG/HR 7000.0000 5821.6741 7000.0000 7000.0000 6182.8225
L/MIN 129.9401 105.8946 137.8228 133.2545 114.6276
STATE VARIABLES:
TEMP C 25.0000 220.0000 120.0000 160.0000 220.0000
PRES BAR 1.0133 1.0133 1.0133 1.0133 1.0133
VFRAC 0.0 0.0 0.0 0.0 0.0
LFRAC 1.0000 1.0000 1.0000 1.0000 1.0000
SFRAC 0.0 0.0 0.0 0.0 0.0
ENTHALPY:
KCAL/KMOL 2.4055+04 2.1462+04 2.5663+04 2.3270+04 2.3116+04
KCAL/KG 230.7652 208.1063 246.2074 223.2675 221.8392
KCAL/HR 1.6154+06 1.2115+06 1.7235+06 1.5629+06 1.3716+06
ENTROPY:
CAL/MOL-K -79.6026 -72.0269 -70.9953 -72.6021 -71.3946
CAL/GM-K -0.7636 -0.6984 -0.6811 -0.6966 -0.6851
DENSITY:
MOL/CC 8.6132-03 8.8845-03 8.1211-03 8.4003-03 8.6271-03
GM/CC 0.8978 0.9163 0.8465 0.8755 0.8990
AVG MW 104.2408 103.1312 104.2342 104.2250 104.2034
COMPONENT ATTRIBUTES:
PS SFRAC
STY-SEG MISSING 1.0000 1.0000 1.0000 1.0000
SFLOW
STY-SEG MISSING 54.1828 16.4459 41.7088 54.1828
DPN
DPN MISSING 1212.4496 1821.7072 1320.3382 1212.4496
DPW
DPW MISSING 2309.3927 3177.1815 2469.6733 2309.3927
PDI
PDI MISSING 1.9047 1.7441 1.8705 1.9047
MWN
MWN MISSING 1.2628+05 1.8973+05 1.3752+05 1.2628+05
MWW
References
Albright, L.F. (1985). Processes for Major Addition-Type Plastics and Their
Monomers. Krieger Publishing Co.
Brandrup, J., & Immergut, E. H. (1989). Polymer Handbook, 3rd Ed. Wiley
Intersciences, 57-82.
Friis, N., & Hamielec, A. E. (1976). Gel-Effect in Emulsion Polymerization of
Vinyl Monomers. ACS Symp. Ser., 24.
Gaur, U., & Wunderlich, B. (1982). J. Phys. Chem. Ref. Data, 11, (2), p. 313.
Hui, A. W., & Hamielec, A. E. (1972). Thermal Polymerization of Styrene at
High Conversions and Temperatures. An Experimental Study. J. of Applied
Polymer Science, 16, pp. 9-769.
Van Krevelen, D. W. (1990). Properties of Polymers, Their Correlation With
Chemical Structure; Their Numerical Estimation and Prediction from Additive
Group Contributions. New York: Elsevier.
Van Krevelen, D. W. (1990). Properties of Polymers, 3rd Ed., Amsterdam:
Elsevier.
The polystyrene with styrene monomer distillation model illustrates the use of
Aspen Polymers to model and predict the extent of polymer formation during
styrene distillation. The model examines the optimum level of inhibitor to
minimize polymer formation.
Process Definition
The process consists of three columns in series:
Process Conditions
The process conditions are:
Component Name Databank Description
Temperature 80 C
Pressure 4.5 Bar
Mass Flow 18000 Kg/hr
Mass fraction of 0.63
styrene
Feed - stage 2
Thermal initiation 3M R*
Chain initiation R* + M P1
Propagation Pn + M Pn+1
Chain transfer to monomer Pn + M Dn + P1
Chain thansfer to EB Pn + EB Dn + P1
Termination by combination Pn + Pm Dn+m
Termination by inhibition Pn + X Dn
Inhibitor
Common inhibitors used in this process are nitrobenzene, dinitro-o-bencene,
dinitro-m-bencene, dinitro-p-bencene, 4,6 dinitro-o-cresol, etc. In this
example 4,6 dinitro-o-cresol is used as inhibitor. (Mw = 198.135)
Process Studies
Input Summary
The input language summary for this example is shown here:
SIM-OPTIONS
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
VOLUME-FLOW='cum/hr' PRESSURE=bar TEMPERATURE=C DELTA-T=C
&
PDROP-PER-HT='mbar/m' PDROP=bar
SIM-OPTIONS FLASH-TOL=1.0000E-06
COMPONENTS
STY C8H8 /
PS PS-1 /
STY-SEG C8H8-R /
EB C8H10-4 /
BEN C6H6 /
TOL C7H8 /
INHIBIT C8H8 /
CINI C8H8
FLOWSHEET
BLOCK BT-COL IN=BT-FEED BT-V1 OUT=BT-TOP BT-L1
BLOCK BT-REB IN=BT-L1 OUT=BT-V1 BT-BOT
BLOCK EB-COL IN=EB-V1 EB-FEED OUT=EB-TOP EB-L1
BLOCK EB-REB IN=EB-L1 OUT=EB-V1 EB-BOT
BLOCK PUMP-1 IN=BT-BOT OUT=EB-FEED
BLOCK PUMP-2 IN=EB-BOT OUT=STY-FEED
BLOCK STY-COL IN=STY-FEED STY-V1 OUT=STY-TOP STY-L1
BLOCK STY-REB IN=STY-L1 OUT=STY-V1 RESIDUE
PROPERTIES POLYNRTL
PROP-DATA MW
IN-UNITS MET VOLUME-FLOW='cum/hr' ENTHALPY-FLO='MMkcal/hr' &
HEAT-TRANS-C='kcal/hr-sqm-K' PRESSURE=bar TEMPERATURE=C &
VOLUME=cum DELTA-T=C HEAD=meter MOLE-DENSITY='kmol/cum' &
MASS-DENSITY='kg/cum' MOLE-ENTHALP='kcal/mol' &
MASS-ENTHALP='kcal/kg' HEAT=MMkcal MOLE-CONC='mol/l' &
PDROP=bar
PROP-LIST MW
PVAL INHIBIT 198.1350
PROP-DATA PLXANT-1
IN-UNITS SI PRESSURE=bar TEMPERATURE=C PDROP='N/sqm'
PROP-LIST PLXANT
PVAL INHIBIT -30.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PROP-DATA NRTL-1
IN-UNITS ENG
PROP-LIST NRTL
BPVAL STY EB -.93850 48.950099610 .30 0.0 0.0 0.0 &
135.80600290 206.60000230
BPVAL EB STY 1.31730 -66.017519470 .30 0.0 0.0 0.0 &
POLYMERS
SEGMENTS STY-SEG REPEAT
POLYMERS PS
ATTRIBUTES PS SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM FMOM &
SMOM LDPN LZMOM LFMOM LSFLOW LSFRAC LEFLOW LEFRAC &
LPFRAC
STREAM BT-FEED
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
VOLUME-FLOW='cum/hr' PRESSURE=bar TEMPERATURE=C &
DELTA-T=C PDROP-PER-HT='mbar/m' PDROP=bar
SUBSTREAM MIXED TEMP=80.0 PRES=4.50 MASS-FLOW=18000.0
MASS-FRAC STY .630 / EB .3180 / BEN .020 / TOL .0320
STREAM BT-L1
SUBSTREAM MIXED TEMP=107.746962 PRES=.40 &
MASS-FLOW=25896.6115
MASS-FRAC STY .616547893 / EB .365240172 / BEN &
.008181312 / TOL .010030622
STREAM BT-V1
IN-UNITS MET MASS-FLOW='kg/sec' MOLE-FLOW='kmol/sec' &
VOLUME-FLOW='cum/hr' ENTHALPY-FLO='MMkcal/hr' &
HEAT-TRANS-C='kcal/hr-sqm-K' PRESSURE=bar TEMPERATURE=C &
TIME=sec VOLUME=cum DELTA-T=C HEAD=meter &
MOLE-DENSITY='kmol/cum' MASS-DENSITY='kg/cum' &
MOLE-ENTHALP='kcal/mol' MASS-ENTHALP='kcal/kg' &
HEAT=MMkcal MOLE-CONC='mol/l' PDROP=bar
SUBSTREAM MIXED TEMP=109 PRES=0.4 MASS-FLOW=2.15351197
MASS-FRAC STY 0.560602271 / EB 0.390454005 / BEN &
0.00198335909 / TOL 0.0469603649
STREAM EB-L1
SUBSTREAM MIXED TEMP=112.427233 PRES=.430 &
MASS-FLOW=31872.8612
MASS-FRAC STY .619641479 / EB .369950341 / BEN &
1.0888E-10 / TOL .000287885
STREAM EB-V1
STREAM STY-L1
SUBSTREAM MIXED TEMP=83.8200436 PRES=.140 &
MASS-FLOW=51412.4723
MASS-FRAC STY .985046601 / EB .001068554 / BEN &
1.0782E-29 / TOL 1.9846E-14
STREAM STY-V1
SUBSTREAM MIXED PRES=0.092 VFRAC=1 MOLE-FLOW=341.040394
MASS-FRAC STY 0.994096437 / EB 0.00590356315
EO-CONV-OPTI
SENSITIVITY S-1
PARAM BASE-CASE=FIRST
DEFINE INHI MASS-FRAC STREAM=BT-FEED SUBSTREAM=MIXED &
COMPONENT=INHIBIT
DEFINE R1PS MASS-FLOW STREAM=BT-BOT SUBSTREAM=MIXED &
COMPONENT=PS
DEFINE R2PS MASS-FLOW STREAM=EB-BOT SUBSTREAM=MIXED &
COMPONENT=PS
DEFINE R3PS MASS-FLOW STREAM=RESIDUE SUBSTREAM=MIXED &
COMPONENT=PS
DEFINE STY MASS-FLOW STREAM=BT-FEED SUBSTREAM=MIXED &
COMPONENT=STY
DEFINE FEEDI MASS-FLOW STREAM=BT-FEED SUBSTREAM=MIXED &
COMPONENT=INHIBIT
C
F PS1 = R1PS / STY
F PS2 = (R2PS - R1PS )/ STY
F PS3 = (R3PS - R2PS )/ STY
F PSTOT =R3PS/STY
F PS0 = 731 / STY
C
F PPMIN = INHI * 1.0E6
C
F P = PS0 - PSTOT
C
F PROFIT = P - 100 * FEEDI/STY
TABULATE 1 "PPMIN" COL-LABEL="INHIBIT" "PPM"
TABULATE 2 "PS1" COL-LABEL="PS1" "POLYMER" "BT-COL"
TABULATE 3 "PS2" COL-LABEL="PS2" "POLYMER" "EB-COL"
TABULATE 4 "PS3" COL-LABEL="PS3" "POLYMER" "STY-COL"
TABULATE 5 "PSTOT" COL-LABEL="PSTOT" "TOTAL" "POLYMER"
TABULATE 6 "PROFIT" COL-LABEL="PROFIT"
VARY MASS-FLOW STREAM=BT-FEED SUBSTREAM=MIXED &
COMPONENT=INHIBIT
RANGE LIST=0.0 0.2 0.4 0.5 0.55 .60 0.65 0.7 0.75 1. &
2.0 3.0
CONV-OPTIONS
PARAM TEAR-METHOD=BROYDEN CHECKSEQ=YES
PROPERTY-REP PARAMS
The effect of the inhibitor feed mass fraction on polymer production is shown
here:
0.01 0.02 0.03 0.04 0.05 0.06 0.07
The effect of the inhibitor feed mass fraction on process profitability is shown
here:
0 10 20 30 40 50 60 70 80 90 100 110
Concentration of Inhibitor in Feed, PPM
SUBSTREAM: MIXED
PHASE: LIQUID LIQUID LIQUID LIQUID LIQUID
COMPONENTS: KMOL/HR
STY 108.8798 4.1382-07 105.5295 88.5771 12.7131
PS 0.0 0.0 3.3508 0.0 7.5900
EB 53.9149 9.0635-05 53.9148 53.9138 9.7391-04
BEN 4.6087 4.5123 9.6354-02 9.6354-02 0.0
TOL 6.2513 2.5789 3.6725 3.6725 0.0
INHIBIT 0.0 0.0 0.0 0.0 0.0
CINI 0.0 0.0 0.0 0.0 0.0
COMPONENTS: KG/HR
STY 1.1340+04 4.3100-05 1.0991+04 9225.4371 1324.0935
PS 0.0 0.0 348.9875 0.0 790.5152
EB 5724.0000 9.6225-03 5723.9916 5723.8877 0.1034
BEN 360.0000 352.4725 7.5265 7.5265 0.0
TOL 576.0000 237.6205 338.3835 338.3835 0.0
INHIBIT 0.0 0.0 0.0 0.0 0.0
CINI 0.0 0.0 0.0 0.0 0.0
COMPONENTS: MASS FRAC
STY 0.6300 7.3038-08 0.6313 0.6032 0.6261
PS 0.0 0.0 2.0045-02 0.0 0.3738
EB 0.3180 1.6306-05 0.3288 0.3742 4.8894-05
BEN 2.0000-02 0.5973 4.3231-04 4.9208-04 0.0
TOL 3.2000-02 0.4027 1.9436-02 2.2123-02 0.0
INHIBIT 0.0 0.0 0.0 0.0 0.0
CINI 0.0 0.0 0.0 0.0 0.0
RESIDUE STY-TOP
---------------
SUBSTREAM: MIXED
PHASE: LIQUID VAPOR
COMPONENTS: KMOL/HR
STY 0.6196 12.0586
PS 7.6249 0.0
EB 2.2485-06 9.7166-04
BEN 0.0 0.0
TOL 0.0 0.0
INHIBIT 0.0 0.0
CINI 0.0 0.0
COMPONENTS: KG/HR
STY 64.5362 1255.9241
PS 794.1484 0.0
EB 2.3872-04 0.1032
BEN 0.0 0.0
TOL 0.0 0.0
INHIBIT 0.0 0.0
CINI 0.0 0.0
COMPONENTS: MASS FRAC
STY 7.5157-02 0.9999
PS 0.9248 0.0
EB 2.7801-07 8.2131-05
BEN 0.0 0.0
TOL 0.0 0.0
INHIBIT 0.0 0.0
CINI 0.0 0.0
TOTAL FLOW:
KMOL/HR 8.2446 12.0596
KG/HR 858.6848 1256.0273
CUM/HR 0.8909 6647.3689
STATE VARIABLES:
TEMP C 120.0000 59.1993
PRES BAR 9.2000-02 5.0000-02
VFRAC 0.0 1.0000
LFRAC 1.0000 0.0
SFRAC 0.0 0.0
ENTHALPY:
KCAL/MOL 18.2651 36.2480
KCAL/KG 175.3702 348.0308
MMKCAL/HR 0.1506 0.4371
ENTROPY:
CAL/MOL-K -82.6169 -43.9790
CAL/GM-K -0.7932 -0.4223
DENSITY:
KMOL/CUM 9.2543 1.8142-03
KG/CUM 963.8521 0.1890
AVG MW 104.1515 104.1517
References
Scheeline, H. W., & Pons, J. I. (1977). Styrene, Supplement B (Process
Economics Program Report No. 33B). Menlo Park: Stanford Research
Institute.
Process Definition
The first stage of the EPS suspension polymerization process is modeled using
the batch reactor model, RBatch. The recipe is defined in terms of
concentrations at 25C in the monomer phase. Therefore, the charge stream,
which is specified in overall concentrations, is fed to a three-phase block at
25C.
The charge of individual components is adjusted until the equilibrium
concentrations in the monomer phase match the required recipe and the
dispersed phase volume fraction equals 0.4. All streams from the Flash3 block
are then passed to a heater that heats the charge stream to the initial
reaction temperature.
Note: In the delivered model, the calculator block is hidden. If the calculator
block is revealed, a Fortran compiler is required to run the example.
The process conditions are:
Components Name Databank Description
Recipe
Temperature (C) 90
Pressure (KPa) 275 (n-pentane=0), 825 (n-pentane>0)
Volume Flow (m3/hr) 3.2E-3
Mass parts of initiator .27
Mass parts of polystyrene 0.0
Mass parts of coinitiator 0.0
Mass parts of monomer 100, 92.5, 85
Mass parts of water 170.1
Mass parts of n-pentane 0.0, 7.5, 15
Kinetics FREE-RAD model (options specify reactions in L1 phase only)
Reactors / Kinetics
The kinetics of suspension EPS are similar to those of the bulk polystyrene
(See B1 Polystyrene Bulk Polymerization by Thermal Initiation application
example).
The reactions included from the free-radical built-in kinetics are:
Description Reaction
Thermal initiation 3M 2 R
Initiator decomposition I 2R
Chain initiation R M P1
Propagation Pn M Pn 1
Chain transfer to monomer Pn M Dn P1
Termination by combination Pn Pm Dn m
The units for the rate constants for the thermal initiation reaction are
m6/kmol2/s. For initiator decomposition the units are /s. The units for the rate
constants for the other reactions are m3/kmol2/s.
Thermal initiation is defined using the special initiation reaction and by setting
third-order thermal initiation with respect to monomer. A coinitiator, which is
required for the special initiation reaction, is included in the list of
components, but its feed rate is set to zero so that it will not influence the
rate for the thermal initiation reaction.
As mentioned previously, the gel and glass effects are modeled using a user
gel effect correlation:
Description Correlation Number
The correlations contained within the user Fortran routine, usrgel.f, calculate
the gel factor for each of these reactions using the free volume theory
described by Villalobos (1993). The gel factor is multiplied with the rate
constant calculated from the pre-exponental and activation energy to give the
diffusion limited rate constant.
Process Studies
The model is used to study the effect of initial n-pentane concentration on
conversion and molecular weight profiles.
Input Summary
An input language summary for this example is shown here:
TITLE 'SUSPENSION EXPANDED POLYTSTYRENE BATCH PROCESS'
IN-UNITS SI
DIAGNOSTICS
HISTORY SIM-LEVEL=4
SYS-OPTIONS TRACE=YES
RUN-CONTROL MAX-TIME=2000.0
;
; THE COMPONENT SYSTEM DEFINITION FOLLOWS
;
DESCRIPTION "
EXPANDED POLYSTYRENE SUSPENSION POLYMERIZATION BATCH
PROCESS. THIS EXAMPLE USES THE KINETIC MODEL AND RATE
CONSTANTS DESCRIBED IN VILLALOBOS,M.A.,HAMIELEC,A.E,AND
WOOD, P.E.,BULK AND SUSPENSION POLYMER IZATION OF STYRENE
IN THE PRESENCE OF N-PENTANE. AN EVALUATION OF MONO
FUNCTIONAL AND BIFUNCTIONAL INITIATION,
J.APPL.POLYM.SCI,50,327-343 (19 93)
"
COMPONENTS
FLOWSHEET
BLOCK BATCH IN=4 OUT=PRODUCT
BLOCK FLASH3 IN=CHARGE OUT=1 2 3
BLOCK HEATER IN=1 2 3 OUT=4
PROPERTIES POLYNRTL
PROP-DATA DATA1
IN-UNITS SI
PROP-LIST MW
PVAL BPO 242.230
PROP-DATA PLXANT-1
IN-UNITS SI
PROP-LIST PLXANT
PVAL PS -10.0
PROP-DATA NRTL-1
IN-UNITS MET VOLUME-FLOW='cum/hr' ENTHALPY-FLO='MMkcal/hr' &
HEAT-TRANS-C='kcal/hr-sqm-K' PRESSURE=bar TEMPERATURE=K &
VOLUME=cum DELTA-T=C HEAD=meter MOLE-DENSITY='kmol/cum' &
MASS-DENSITY='kg/cum' MOLE-ENTHALP='kcal/mol' &
MASS-ENTHALP='kcal/kg' HEAT=MMkcal MOLE-CONC='mol/l' &
PDROP=bar
PROP-LIST NRTL
BPVAL BPO WATER -176.71520 10542.57030 .20 0.0 25.52320 &
0.0 279.150 338.150
BPVAL WATER BPO 150.57740 -5675.06590 .20 0.0 -21.68180 &
0.0 279.150 338.150
BPVAL CINI WATER -176.71520 10542.57030 .20 0.0 25.52320 &
0.0 279.150 338.150
BPVAL WATER CINI 150.57740 -5675.06590 .20 0.0 -21.68180 &
0.0 279.150 338.150
BPVAL STY WATER -176.71520 10542.57030 .20 0.0 25.52320 &
0.0 279.150 338.150
BPVAL WATER STY 150.57740 -5675.06590 .20 0.0 -21.68180 &
0.0 279.150 338.150
BPVAL STY-SEG WATER 0.0 500.0 .30 0.0 0.0 0.0 0.0 &
1000.0
BPVAL WATER STY-SEG 0.0 500.0 .30 0.0 0.0 0.0 0.0 &
1000.0
BPVAL WATER NPENTANE 12.38660 -791.79130 .20 0.0 0.0 0.0 &
273.150 303.150
BPVAL NPENTANE WATER -10.68920 5051.72750 .20 0.0 0.0 0.0 &
273.150 303.150
POLYMERS
SEGMENTS STY-SEG REPEAT
POLYMERS PS
STREAM CHARGE
IN-UNITS MET VOLUME-FLOW='cum/hr' ENTHALPY-FLO='MMkcal/hr' &
HEAT-TRANS-C='kcal/hr-sqm-K' PRESSURE=bar TEMPERATURE=C &
VOLUME=cum DELTA-T=C HEAD=meter MOLE-DENSITY='kmol/cum' &
MASS-DENSITY='kg/cum' MOLE-ENTHALP='kcal/mol' &
MASS-ENTHALP='kcal/kg' HEAT=MMkcal MOLE-CONC='mol/l' &
PDROP=bar
SUBSTREAM MIXED TEMP=25.0 PRES=8.250 VOLUME-FLOW=.00320 &
FREE-WATER=NO NPHASE=3 PHASE=V
MASS-FRAC BPO .2780 / STY 100.0 / WATER 170.10
EO-CONV-OPTI
CALCULATOR PHI
IN-UNITS MET VOLUME-FLOW='cum/hr' ENTHALPY-FLO='MMkcal/hr' &
CONV-OPTIONS
PARAM CHECKSEQ=NO
REACTIONS R1 FREE-RAD
IN-UNITS MET
DESCRIPTION "EXAMPLE FREE-RADICAL INPUT"
PARAM QSSA=YES
SPECIES INITIATOR=BPO COINITIATOR=CINI MONOMER=STY &
POLYMER=PS
MON-RSEG STY STY-SEG
INIT-DEC BPO 3.8160E+12 27233.0 0.0 EFFIC=.60 NRADS=2
INIT-SP STY CINI 219000.0 27440.0 0.0
CHAIN-INI STY 1.0213E+07 7068.0 0.0
PROPAGATION STY STY 1.0213E+07 7068.0 0.0
CHAT-MON STY STY 1.0213E+07 13450.0 0.0
TERM-COMB STY STY 1.2583E+09 1677.0 0.0
INIT-SP-EFF STY COEFFA=0.0 COEFFB=3.0 COEFFC=0.0
GEL-EFFECT TERMINATION 4 MAX-PARAMS=3 GE-PARAMS=.50 .3480 &
1.750
GEL-EFFECT PROPAGATION 5 MAX-PARAMS=2 GE-PARAMS=1.0 .04650
GEL-EFFECT CHAT-MON 3 MAX-PARAMS=1 GE-PARAMS=1.750
SUBROUTINE GEL-EFFECT=USRGEL
OPTIONS PHASE=L1;
0.9
0.8
0.7
Conversion
0.6
0.5
0.4
0.3
0.2
0.1
0 1 2 3 4 5 6 7 8
Time(hours)
300
275
(Thousands)
250
225
200
175
150
Mn,Mw
125
100
75
50
25
Conversion
One would expect the addition of a solvent to reduce the molecular weight.
However, the n-pentance reduces this effect by lowering the chain transfer to
monomer rate. Additionally, because it is a stronger electron donor than
styrene, n-pentane can obstruct the hydrogen abstraction necessary for chain
transfer to solvent to occur. Therefore, n-pentane addition leads to a
molecular weight build-up in the system.
Process Definition
The copolymerization of styrene and ethyl acrylate is carried out in a batch
reactor. The batch reactor is charged with a pre-mixed stream of styrene,
ethyl acrylate, and the initiator azo-bis-isobutyronitrile. The reactor is
operated at a constant temperature and pressure. The flowsheet is shown
here:
Temperature (C) 60 60
Pressure (atm) 1 1
Mass flow kg/h 1000
Mole flow of AIBN 0.05 MR/L in reaction mixture
Mole fraction of styrene 0.762
Mole fraction of ethyl acrylate 0.238
Kinetics FREE-RAD model
Operating Conditions
BATCH reactor
Temperature (C) 60
Pressure (atm) 1
Reaction time (hr) 25
Cycle time (hr) 1
Reactors / Kinetics
The reaction set used in this process and the initial rate constant parameters
are:
Description k0 Ea V EFFIC
(m3/kmol) (J/kmol) (m3/kmol)
Gel effect is applied to the initiator efficiency and termination reaction using
correlation No. 2.
Reaction Correlation Parameters
INIT-EFF 2 -17.40
0.05528
17.8240
-05090
0.0
0.0
0.0
0.0
0.0
2.0
TERMINATION 2 1.0
0.0
2.570
-.005050
9.560
-.01760
-3.030
.007850
0.0
2.0
r2 k p22 / k p21
Where:
STY=1
EA=2
Four termination reactions are considered in this scheme. Cross-
termination reactions involving STY and EA segments were set to have the
same rate constant. This reduces specification of the rate constants to
three. If the homo-termination rate constants are known, the cross-
termination rate constants can be estimated using the following
equations:
k12 k11 k 22
k 21 k12
In this example homo-termination rate constants were regressed using
the molecular weight data.
In the kinetic scheme termination by combination was selected. It is
possible that termination by disproportionation can also occur. One needs
to know apriori the termination mechanism. As a general rule if
polydispersity (PDI) is 1.5 termination by combination is controlling and if
PDI is two termination is controlled by the disproportionation reaction.
At high conversion termination is diffusion controlled. This is modeled
using the gel-effect option. Therefore use of low conversion data is
recommended for regression of the kinetic parameters without using gel-
effect. High conversion data is then used to regress the gel-effect
parameters.
Monomer conversion data was used to regress the initiator efficiency,
propagation and the termination rate constants. Molecular weight data
INIT-DEC kd Initiator
PROPAGATIO kp Monomer X X X
N
CHAT-MON ktrm X X
CHAT-SOL ktrs X X
CHAT-AG ktra X X
TERM-DIS ktd Monomer X X
TERM-COM ktc Monomer X X
INHIBITION kx Monomer X X
CHAT-POL ktrp X X
SC-BRANCH kscb X
Process Studies
Input Summary
The input language summary for this example is shown here:DYNAPLUS
DPLUS RESULTS=ON
RUN-CONTROL MAX-TIME=1200000.0
DESCRIPTION "
Batch reactor data for styrene-ethyl acrylate
polymerization is used to regress the propagation
COMPONENTS
AIBN C8H8 /
STY C8H8 /
EA C5H8O2 /
STY-SEG C8H8-R /
EA-SEG C5H8O2-R-2 /
POLYMER PS-1
FLOWSHEET
BLOCK BATCH IN=FEED OUT=POLYMER
BLOCK MIX IN=MONOMER INIT OUT=FEED
PROPERTIES POLYNRTL
PROP-DATA MW
IN-UNITS ENG
PROP-LIST MW
PVAL AIBN 164.2120
POLYMERS
SEGMENTS STY-SEG REPEAT / EA-SEG REPEAT
POLYMERS POLYMER
ATTRIBUTES POLYMER SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM &
FMOM SMOM
PROP-SET INITCONC
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=atm TEMPERATURE=C PDROP='N/sqm'
PROPNAME-LIS MOLECONC UNITS='mol/l' SUBSTREAM=MIXED &
COMPS=AIBN
PROP-SET POLY-MF
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=atm TEMPERATURE=C PDROP='N/sqm'
PROPNAME-LIS MASSFRAC SUBSTREAM=MIXED COMPS=POLYMER PHASE=L
STREAM INIT
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=atm TEMPERATURE=C VOLUME=l &
MASS-DENSITY='gm/cc' PDROP='N/sqm'
SUBSTREAM MIXED TEMP=50.0 PRES=1.0
MOLE-FLOW AIBN .05710370
STREAM MONOMER
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=atm TEMPERATURE=C VOLUME=l &
MASS-DENSITY='gm/cc' PDROP='N/sqm'
SUBSTREAM MIXED TEMP=60.0 PRES=1.0 MASS-FLOW=1000.0
MOLE-FRAC STY .7620 / EA .2380
EO-CONV-OPTI
CALCULATOR F-1
DEFINE F PARAMETER 1
F F = 5
EXECUTE FIRST
CALCULATOR K-PROP
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=atm TEMPERATURE=C PDROP='N/sqm'
DEFINE PROP11 REACT-VAR REACTION=STY-EA VARIABLE=PRPRE-EXP &
SENTENCE=PROPAGATION ID1=STY ID2=STY
DEFINE PROP12 REACT-VAR REACTION=STY-EA VARIABLE=PRPRE-EXP &
SENTENCE=PROPAGATION ID1=STY ID2=EA
DEFINE PROP21 REACT-VAR REACTION=STY-EA VARIABLE=PRPRE-EXP &
SENTENCE=PROPAGATION ID1=EA ID2=STY
DEFINE PROP22 REACT-VAR REACTION=STY-EA VARIABLE=PRPRE-EXP &
SENTENCE=PROPAGATION ID1=EA ID2=EA
F R1 = 0.717
F R2 = 0.128
F PROP12 = PROP11/R1
F PROP21 = PROP22/R2
READ-VARS PROP11 PROP22
WRITE-VARS PROP12 PROP21
CALCULATOR K-TERM
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=atm TEMPERATURE=C PDROP='N/sqm'
DEFINE TC11 REACT-VAR REACTION=STY-EA VARIABLE=TCPRE-EXP &
SENTENCE=TERM-COMB ID1=STY ID2=STY
DEFINE TC12 REACT-VAR REACTION=STY-EA VARIABLE=TCPRE-EXP &
SENTENCE=TERM-COMB ID1=STY ID2=EA
DEFINE TC21 REACT-VAR REACTION=STY-EA VARIABLE=TCPRE-EXP &
SENTENCE=TERM-COMB ID1=EA ID2=STY
DEFINE TC22 REACT-VAR REACTION=STY-EA VARIABLE=TCPRE-EXP &
SENTENCE=TERM-COMB ID1=EA ID2=EA
DEFINE F PARAMETER 1
F TC12 = F*(TC11*TC22)**0.5
F TC21 = F*(TC11*TC22)**0.5
READ-VARS TC11 TC22
WRITE-VARS TC12 TC21
CONV-OPTIONS
PARAM CHECKSEQ=NO
REGRESSION K-PROP
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=atm TEMPERATURE=C PDROP='N/sqm'
DATA WPOLY
VARY REACT-VAR REACTION=STY-EA VARIABLE=PRPRE-EXP &
SENTENCE=PROPAGATION ID1=STY ID2=STY
LIMITS 5E4 1E7
ALGORITHM MXPASS=3000
DATA-SET MWN
IN-UNITS ENG
DEFINE TIME BLOCK-VAR BLOCK=BATCH VARIABLE=VALUE &
SENTENCE=STOP ID1=1
DEFINE MWN COMP-ATTR-VAR STREAM=POLYMER SUBSTREAM=MIXED &
COMPONENT=POLYMER ATTRIBUTE=MWN ELEMENT=1
INPUT TIME
RESULT MWN
USE STD-DEV .010 5000.0 / DATA 2.0 96148.0 / DATA 4.0 &
96939.0 / DATA 6.0 98547.0 / DATA 8.0 102418.0 / &
DATA 10.0 101941.0 / DATA 12.0 104354.0 / DATA 14.0 &
105066.0 / DATA 15.50 112415.0 / DATA 17.0 112159.0 / &
DATA 18.50 112001.0 / DATA 20.0 114564.0 / DATA &
21.50 135065.0 / DATA 23.50 133690.0 / DATA 27.0 &
120899.0
DATA-SET WPOL-SFR
IN-UNITS ENG
DEFINE TIME BLOCK-VAR BLOCK=BATCH VARIABLE=VALUE &
SENTENCE=STOP ID1=1
DEFINE WPOLY MASS-FRAC STREAM=POLYMER SUBSTREAM=MIXED &
COMPONENT=POLYMER
DATA-SET WPOLY
IN-UNITS ENG
DEFINE TIME BLOCK-VAR BLOCK=BATCH VARIABLE=VALUE &
SENTENCE=STOP ID1=1
DEFINE WPOLY MASS-FRAC STREAM=POLYMER SUBSTREAM=MIXED &
COMPONENT=POLYMER
INPUT TIME
RESULT WPOLY
USE STD-DEV .010 .010 / DATA 2.0 .05680 / DATA 4.0 &
.10790 / DATA 6.0 .16090 / DATA 8.0 .21060 / DATA &
10.0 .26160 / DATA 12.0 .37750 / DATA 14.0 .3660 / &
DATA 15.50 .39360 / DATA 17.0 .45970 / DATA 18.50 &
.54680 / DATA 20.0 .61440 / DATA 21.50 .68660 / &
DATA 23.50 .77460
;
Styrene homo- Pre-exponential 2.154 x 107 19169.6 2.1168 x 106 2.1919 x 106
propagation factor
EA homo-propagation Pre-exponential 3.0 x 106 6.65 x 106 5 x 105 1 x 107
factor
Styrene homo- Pre-exponential 1.26 x 109 --- --- ---
termination factor
EA homo-termination Pre-exponential 2.6 x 109 --- --- ---
factor
Process Definition
The process flowsheet consists of the batch reactor with an initial batch
charge and a continuous feed for the addition of monomers and other
ingredients:
This model provides the base case, which can be used to study various
process variables: effect of initiator and emulsifier levels, temperature,
Smith-Ewart kinetics, etc.
Process Conditions
The process conditions are:
Components Name Databank Description
Feeds
Operating
Conditions
B1 Pressure (atm) 10
FLASH Option NO
Temperature
profile
Time (hr) Temperature
(C)
0.0 65.0
0.5 65.0
0.5 70.0
1.0 70.0
1.0 75.0
Note: Prior to release 12, all initiators were required to be in the aqueous
phase for emulsion polymerizations. However, NRTL parameters are now
required to define the initiator APS in the aqueous phase.
The NRTL binary interaction parameters are:
NRTL
H2O STY 3.6260 1513.50 .360 0.0 0.0 0.0 0.0 1000.0
STY H2O -4.4360 2869.70 .360 0.0 0.0 0.0 0.0 1000.0
H2O BD 3.5890 862.820 .30 0.0 0.0 0.0 0.0 1000.0
BD H2O -.8520 702.170 .30 0.0 0.0 0.0 0.0 1000.0
APS STY 0 5000 .01 0.0 0.0 0.0 0.0 1000.0
STY APS 0 5000 .01 0.0 0.0 0.0 0.0 1000.0
APS BD 0 5000 .01 0.0 0.0 0.0 0.0 1000.0
BD APS 0 5000 .01 0.0 0.0 0.0 0.0 1000.0
Reactors / Kinetics
There are considerable data available on the kinetics of emulsion
copolymerization of styrene and butadiene. It is assumed that the primary
chain transfer reaction is similar to that of a monomer.
The rate constants for the kinetic scheme are obtained from Broadhead and
are summarized here (Broadhead, 1984; Ponnuswamy & Hamielec, 1997):
Reaction Phase Comp 1 Comp 2 k0 Ea
(J/kmol)
Other parameters affecting the kinetics are the emulsifier parameters and the
monomer partitioning. These are listed here:
CMC kmol / m3 Area m3 / kmole
Process Studies
The model was used to examine the following process parameters as a
function of time through user profiles:
Glass transition temperature of the polymer (K)
Average number of radicals per particle ---
Percent of soap coverage of particles %
Volume of the monomer phase (m3)
Concentration of monomer 1 in monomer phase (Kmol/ m3)
Concentration of monomer 2 in monomer phase (Kmol/ m3)
Input Summary
The input language summary for this example is shown here:
TITLE 'EMULSION COPOLYMERIZATION OF STYRENE AND BUTADIENE'
SYS-OPTIONS TRACE=YES
DESCRIPTION "
THIS EXAMPLE ILLUSTRATES THE USE OF ASPEN POLYMERS TO MODEL
THE COPOLYMERIZATION OF STYRENE AND BUTADIENE IN A BATCH
REACTOR.
"
COMPONENTS
H2O H2O /
STY C8H8 /
BD C4H6-4 /
POLYMER POLYMER /
STY-SEG C8H8-R /
BD-SEG C4H6-R-1 /
APS H2O /
SLS H2O /
TDM C8H8
FLOWSHEET
BLOCK REACTOR IN=CHARGE CONTFEED OUT=PRODUCT
PROPERTIES POLYNRTL
PROP-DATA TGVK
IN-UNITS SI MASS-FLOW='kg/hr' PRESSURE=bar TEMPERATURE=C &
TIME=hr PDROP='N/sqm'
PROP-LIST TGVK
PVAL STY-SEG 100.0
PVAL BD-SEG -54.0
PROP-DATA NRTL-1
IN-UNITS SI
PROP-LIST NRTL
BPVAL H2O STY 3.6260 1513.50 .360 0.0 0.0 0.0 0.0 &
1000.0
BPVAL STY H2O -4.4360 2869.70 .360 0.0 0.0 0.0 0.0 &
1000.0
BPVAL H2O BD 3.5890 862.820 .30 0.0 0.0 0.0 0.0 1000.0
BPVAL BD H2O -.8520 702.170 .30 0.0 0.0 0.0 0.0 1000.0
BPVAL APS STY 0 5000 .01 0.0 0.0 0.0 0.0 1000.000
BPVAL STY APS 0.0 5000 .01 0.0 0.0 0.0 0.0 1000.000
BPVAL APS BD 0 5000 .01 0.0 0.0 0.0 0.0 1000.000
BPVAL BD APS 0.0 5000 .01 0.0 0.0 0.0 0.0 1000.000
POLYMERS
SEGMENTS STY-SEG REPEAT / BD-SEG REPEAT
POLYMERS POLYMER
ATTRIBUTES POLYMER SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM &
FMOM SMOM PSDZMOM PSDFMOM DIAV
DISTRIBUTION POLYMER CHAIN-SIZE CLD NPOINTS=100 FUNCLOG=NO &
UPPER=100000.0
STREAM CHARGE
SUBSTREAM MIXED TEMP=20.0 PRES=10.0
MASS-FLOW H2O 1000.0 / STY 300.0 / BD 300.0 / POLYMER &
0.0 / APS 7.0 / SLS 30.0 / TDM 30.0
STREAM CONTFEED
SUBSTREAM MIXED TEMP=20.0 PRES=10.0
MASS-FLOW H2O 1000.0 / STY 700.0 / BD 700.0 / POLYMER &
0.0 / APS 3.0 / SLS 100.0
EO-CONV-OPTI
CONV-OPTIONS
PARAM CHECKSEQ=NO
0.14
1.4
Number of Particles (PSD 0th moment)
0.12
Average Particle Diameter, micron
1.2
0.1
1
0.02
0.2
0
0 1 2 3 4 5 6 7 8 9 10
Time, hr
Styrene, Kg
1000
Butadiene, Kg
Polymer, Kg
Mass Inventory, Kg
500 750
250
0 1 2 3 4 5 6 7 8 9 10
Time, hr
1.6 1.2
Phase Volume, m3
0.8
Monomer
Aqueous
Polymer
0.4
0 1 2 3 4 5 6 7 8 9 10
Time, hr
As shown, the monomer droplets are completely depleted after 2.5 hours. As
expected, the volume of the swollen polymer phase decreases after the
depletion of the monomer droplets.
The number average and weight average degree of polymerization and the
polydispersity is shown here:
6000
DPN
2000 3000 4000 5000
DPW
3.5
Degree of Polymerization
PDI
Polydispersity Index
2.5 3
1000
2
0
0 1 2 3 4 5 6 7 8 9 10
Time, hr
Styrene
Butadiene
0 1 2 3 4 5 6 7 8 9 10
Time, hr
Distrib Moments
DPN = 1155.97 968.44
DPW = 3926.32 3242.01
0.0001 0.00015
Area = 1.0021
5e-5
References
Blackley, D. C. (1983). Synthetic Rubbers: Their Chemistry and Technology.
New York: Applied Science Publishers.
Broadhead, T. O. (1984). Dynamic Modeling of the Emulsion Copolymerization
of Styrene/Butadiene. M. E. Thesis, McMaster Univ.
Ponnuswamy, S. R., & Hamielec, A. E. (1997). Emulsion Polymerization
Theory and Practice. Lecture Notes for Intensive Short Course on Polymer
Reaction Engineering, Burlington, Ontario Canada.
The styrene butadiene rubber (SBR) ionic solution process model illustrates
the use of Aspen Polymers to model an anionic polymerization of styrene and
butadiene. The model describes the production of tapered block copolymer
using a batch reactor.
Process Definition
This application demonstrates production of tapered block SBR polymerization
in a continuous reactor. A schematic of a batch reactor using the ionic kinetics
is shown here:
Even though there is more reactive styrene present in the feed, butadiene is
first consumed during the polymerization. After butadiene is consumed, the
reaction rate jumps as soon as styrene homopolymerization is initiated. This
behavior is governed by cross propagation reactions as discussed earlier. The
key outputs predicted by the model are:
Total amount of polymer formed 100 kg/hr
Mole fraction of styrene in final polymer 0.341
DPn of final polymer 1281
PDI of polymer out of CSTR 1.06
Fraction of growing polymer that is growing 0.146
Fraction of associated polymer 0.853
DPn of associated polymer 1382
Kinetics
A comprehensive kinetic scheme for ionic copolymerization of any number of
monomers has been built into Aspen Polymers. The block, tapered-block, and
random SBR can be modeled using this scheme. The built in kinetics can also
handle the presence of multiple propagating species (for example: free-ions,
ion pairs, and dormant esters).
The polymerization kinetic scheme includes a subset of the reactions available
in the Aspen Polymers ionic kinetic model:
Reaction Type Site ID Reacting Description of Reaction
Species
Process Studies
Input Summary
An input language summary for this example is shown here:
TITLE 'SBR (tapered-block) Solution Polymerization'
OUT-UNITS MET
DESCRIPTION "
This example demonstrates tapered block SBR made
by charging the ingredients in a batch reactor
using the Ionic kinetic reaction scheme.
"
COMPONENTS
BULI-6 C24H38O4-D1 /
BULI-1 C4H5N-1 /
STYRENE C8H8 /
STY-SEG C8H8-R /
BUTDIENE C4H6-4 /
BUT-SEG C4H6-R-1 /
HEXANE C6H10-2 /
SBR SBR
FLOWSHEET
BLOCK BATCH IN=FEED OUT=PROD
PROPERTIES POLYNRTL
PROP-DATA MW
IN-UNITS MET
PROP-LIST MW
PVAL BULI-6 384
PVAL BULI-1 64
POLYMERS
PARAM NSITES=1 SITCOM=INITIATOR
SEGMENTS STY-SEG REPEAT / BUT-SEG REPEAT
INITIATOR BULI-1
POLYMERS SBR
ATTRIBUTES SBR SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM &
FMOM SMOM LEFLOW LSFLOW LSFRAC LDPN LSMOM LPFRAC &
SSFRAC SSFLOW SDPN SDPW SPDI SMWN SMWW SZMOM SSMOM &
SPFRAC LSSFLOW LSSFRAC LSEFLOW LSEFRAC LSDPN LSSMOM &
LSPFRAC ASEFLOW ASSFLOW ASSMOM ASDPN ADPN ASPDI APDI &
ASFLOW ASPFRAC APFRAC DSEFLOW DSSFLOW DSSMOM
ATTRIBUTES BULI-1 P0FLOW PT0FLOW CIONFLOW
STREAM FEED
IN-UNITS MET
SUBSTREAM MIXED TEMP=50. <C> PRES=1.
MASS-FLOW BULI-6 0.13 / BULI-1 0. / STYRENE 50. / &
BUTDIENE 50. / HEXANE 1000. / SBR 0.
EO-CONV-OPTI
100
45
40
80
60
MASS kg
MASS kg
25
20
40
15
10
20
5
0
The figure shows the cumulative segment fraction of monomer in the polymer
as a function of time. As shown, during the first phase of the polymerization,
only about 10% of styrene is incorporate in the polymer. The final molar
composition of styrene in the polymer is 34%.
The degree of polymerization for live (LDPN), composite (DPN), and
associated (ADPN) is shown here:
1 6 11 16 21 26 31 36 41 46
NPOINT
The HDPE high temperature solution process model illustrates the use of
Aspen Polymers to model a Ziegler-Natta solution polymerization of ethylene.
The model is used to study the effect of feed flow rate on conversion and
polymer properties, and to study the effect of hydrogen concentration on melt
index and polydispersity index.
Process Definition
In this example, the solution polymerization of ethylene is carried out at
160C using n-hexane as solvent. Hydrogen is used as a chain transfer agent
to control molecular weight. The flowsheet consists of two reactors and a
flash unit as shown:
The first reactor produces a high molecular weight polymer while the second
produces a low molecular weight polymer. Solvent, unreacted monomer, and
hydrogen are removed from the product in a flash tank.
An intermediate feed stream going to the second reactor is used to set the
concentration of hydrogen in that reactor to be several times higher than in
the first reactor. A total of four site types are used for the catalyst. Two sites
are considered to be active in the first reactor and all four sites are active in
the second reactor.
Process Conditions
The process conditions are:
Components Name Databank Description
Reactors / Kinetics
The set of reactions included from the built-in kinetics for this model are:
Description Reaction
Propagation P n,j + M i P n + i i
k
Process Studies
In order to determine the effect of feed flowrate on conversion and polymer
properties, and the effect of hydrogen concentration on melt index and
polydispersity index, a sensitivity study is carried out. Melt index is calculated
in a user supplied subroutine (Xie et al., 1994) and is chosen as one of the
sampled parameters in a SENSITIVITY block. Other sampled parameters
include the polymer properties, the component fractions used to determine
conversion, and the flowrates.
DYNAPLUS
DPLUS RESULTS=ON
DIAGNOSTICS
TERMINAL SIM-LEVEL=4 CONV-LEVEL=4 COST-LEVEL=4 PROP-LEVEL=4 &
ECON-LEVEL=4 STREAM-LEVEL=4 SYS-LEVEL=4
SYS-OPTIONS TRACE=YES
RUN-CONTROL MAX-TIME=10000.0
COMPONENTS
TICL4 TICL4 /
TEA C6H15AL /
C2H4 C2H4 /
H2 H2 /
HEXANE C6H14-1 /
HDPE HDPE /
C2H4-R C2H4-R
FLOWSHEET
BLOCK CSTR-1 IN=R1FEED OUT=R1PROD
BLOCK CSTR-2 IN=R1PROD FEED2 OUT=R2PROD
BLOCK FLASH IN=R2PROD OUT=RECYCLE POLYMER
;
; THE PROPERTY CALCULATION METHOD (OPTION SET) SPECIFICATION FOLLOWS
;
PROP-DATA PCSAFT
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=atm TEMPERATURE=C PDROP='N/sqm'
PROP-LIST PCSFTM
PVAL TICL4 20
PVAL TEA 20
POLYMERS
PARAM NSITES=4
SEGMENTS C2H4-R REPEAT
CATALYST TICL4 .00010
POLYMERS HDPE
ATTRIBUTES HDPE SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM &
FMOM SMOM LDPN LZMOM LFMOM LSFLOW LSFRAC LEFLOW &
LEFRAC LPFRAC SSFRAC SSFLOW SDPN SDPW SPDI SMWN SMWW &
SZMOM SFMOM SSMOM SPFRAC LSDPN LSZMOM LSFMOM LSSFLOW &
LSSFRAC LSEFLOW LSEFRAC LSPFRAC
ATTRIBUTES TICL4 CPSFLOW CPSFRAC CVSFLOW CVSFRAC CDSFLOW &
CDSFRAC
DISTRIBUTION HDPE CHAIN-SIZE CLD NPOINTS=100 FUNCLOG=YES &
UPPER=100000.
USER-PROPERT IV SUBROUTINE=USRPRP
PROP-SET IV IV SUBSTREAM=MIXED
STREAM FEED2
SUBSTREAM MIXED TEMP=70.0 PRES=200.0 MASS-FLOW=10000.0 &
NPHASE=1 PHASE=L
MASS-FRAC TICL4 .00010 / TEA .000150 / C2H4 .20 / H2 &
.00020 / HEXANE 0.79955
COMP-ATTR TICL4 CPSFLOW ( 0.0 )
COMP-ATTR TICL4 CPSFRAC ( 1.0 )
COMP-ATTR TICL4 CVSFLOW ( 0.0 0. 0. 0. )
COMP-ATTR TICL4 CVSFRAC ( 0.0 0.0 0.0 0.0 )
COMP-ATTR TICL4 CDSFLOW ( 0.0 )
COMP-ATTR TICL4 CDSFRAC ( 0.0 )
STREAM R1FEED
EO-CONV-OPTI
SENSITIVITY RTIME
PARAM BASE-CASE=LAST
DEFINE RHDPE2 MASS-FLOW STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=HDPE
DEFINE RMW2 COMP-ATTR-VAR STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=MWW ELEMENT=1
DEFINE RMN2 COMP-ATTR-VAR STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=MWN ELEMENT=1
DEFINE RPDI2 COMP-ATTR-VAR STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=PDI ELEMENT=1
DEFINE RE22 MASS-FLOW STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=C2H4
DEFINE RHDPE1 MASS-FLOW STREAM=R1PROD SUBSTREAM=MIXED &
COMPONENT=HDPE
DEFINE RMW1 COMP-ATTR-VAR STREAM=R1PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=MWW ELEMENT=1
DEFINE RMN1 COMP-ATTR-VAR STREAM=R1PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=MWN ELEMENT=1
DEFINE RPDI1 COMP-ATTR-VAR STREAM=R1PROD SUBSTREAM=MIXED &
SENSITIVITY RTIME2
PARAM BASE-CASE=LAST
DEFINE RHDPE2 MASS-FLOW STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=HDPE
DEFINE RMW2 COMP-ATTR-VAR STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=MWW ELEMENT=1
DEFINE RMN2 COMP-ATTR-VAR STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=MWN ELEMENT=1
DEFINE RPDI2 COMP-ATTR-VAR STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=PDI ELEMENT=1
DEFINE RE22 MASS-FLOW STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=C2H4
DEFINE RHDPE1 MASS-FLOW STREAM=R1PROD SUBSTREAM=MIXED &
COMPONENT=HDPE
DEFINE RMW1 COMP-ATTR-VAR STREAM=R1PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=MWW ELEMENT=1
DEFINE RMN1 COMP-ATTR-VAR STREAM=R1PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=MWN ELEMENT=1
DEFINE RPDI1 COMP-ATTR-VAR STREAM=R1PROD SUBSTREAM=MIXED &
COMPONENT=HDPE ATTRIBUTE=PDI ELEMENT=1
DEFINE RE21 MASS-FLOW STREAM=R1PROD SUBSTREAM=MIXED &
COMPONENT=C2H4
DEFINE RMI2 STREAM-PROP STREAM=R2PROD PROPERTY=MI-KAR
DEFINE RMI1 STREAM-PROP STREAM=R1PROD PROPERTY=MI-KAR
DEFINE R2H2 MASS-FRAC STREAM=R2PROD SUBSTREAM=MIXED &
COMPONENT=H2
DEFINE R1H2 MASS-FRAC STREAM=R1PROD SUBSTREAM=MIXED &
COMPONENT=H2
F CONV2 = RHDPE2/(RHDPE2+RE22) * 100
F CONV1 = RHDPE1/(RHDPE1+RE21) * 100
CONV-OPTIONS
PARAM CHECKSEQ=NO
SEQUENCE S-1 CSTR-1 CSTR-2 RTIME CSTR-1 CSTR-2 (RETURN RTIME) &
RTIME2 CSTR-2 FLASH (RETURN RTIME2)
80
70
CSTR1
Ethylene Conversion
60 CSTR2
50
40
30
20
10
0
0 500000 1000000 1500000
Feed Flow Rate to CSTR-1, kg/hr
Properties of the product exiting from CSTR-2 are influenced by the flowrates
of both feeds. Increasing the feed to CSTR-1 decreases the residence time in
CSTR-1 and CSTR-2, as shown, this leads to a decrease in conversion in both
reactors.
Although conversion decreases with increases in CSTR-1 feed, the next two
figures show an increase in the CSTR-2 polymer Mn and Mw. This is due to the
concentration of hydrogen in CSTR-2. As CSTR-1 feed flow increases,
hydrogen concentration in CSTR-2 decreases. As a result, chain transfer rate
decreases, allowing longer polymer chain length.
35000 CSTR1
CSTR2
Mn
30000
25000
20000
0 500000 1000000 1500000
Feed Flow Rate to CSTR-1, kg/hr
120000
115000
110000
105000
Mw
100000
95000 CSTR1
CSTR2
90000
85000
80000
0 500000 1000000 1500000
Feed Flow Rate to CSTR-1, kg/hr
The effect of feed flow rate on the hydrogen concentration (PPM) is shown
here:
H2 Concentration (ppm)
50
45 CSTR2
40
35
30
25
20
15
10
5
0
0 500000 1000000 1500000
Feed Flow Rate to CSTR-1, kg/hr
The effect of feed flow rate on the polydispersity index is shown here:
4.2
CSTR1
Polydispersity Index
3.8 CSTR2
3.6
3.4
3.2
3
0 500000 1000000 1500000
Feed Flow Rate to CSTR-1, kg/hr
The effect of feed flow rate on the melt index is shown here:
2.5
2.25
CSTR1
2
CSTR2
Melt Index 1.75
1.5
1.25
0.75
0.5
0.25
0 500000 1000000 1500000
Feed Flow Rate to CSTR-1, kg/hr
The results for the sensitivity studies to determine the effects of hydrogen
flowrate on polymer properties are shown in the three figures that follow. In
this case, hydrogen flowrate in the intermediate feed to CSTR-2 is varied. As
hydrogen flowrate increases, average molecular weight of HDPE in CSTR-2
decreases as shown in the first figure.
The effect of hydrogen concentration on Mn and Mw is shown here:
140000
120000
100000
Molecular Weight
80000
60000
Mn
40000 Mw
20000
0
0 50 100 150 200 250 300
Hydrogen Conc. in CSTR-2, ppm
The effect of hydrogen concentration on melt index and PDI is shown here:
The molecular weight distribution for polymer produced at the different sites
is shown here:
0.5
0.4 Composite
Site_1
Weight Fraction
0.3 Site_2
Site_3
0.2 Site_4
0.1
0.0
1 10 100 1000 10000 100000
Polymer Chain Length
As shown, although the polymer at each site follows the Flory distribution, the
composite polymer MWD is quite broad. In conclusion, the studies show the
importance of choosing the right operating conditions in order to optimize
reactor performance and product properties.
SUBSTREAM: MIXED
PHASE: LIQUID LIQUID LIQUID LIQUID LIQUID
COMPONENTS: KG/HR
TICL4 1.0000 7.0000 6.0000 6.0000 7.0000
RECYCLE
-------
STREAM ID RECYCLE
FROM : FLASH
TO : ----
SUBSTREAM: MIXED
PHASE: VAPOR
COMPONENTS: KG/HR
TICL4 3.0502-09
TEA 4.5753-09
C2H4 2146.7657
H2 2.1806
HEXANE 5.9925+04
HDPE 0.0
COMPONENTS: MASS FRAC
TICL4 4.9139-14
TEA 7.3707-14
C2H4 3.4584-02
H2 3.5129-05
HEXANE 0.9654
HDPE 0.0
TOTAL FLOW:
KMOL/HR 772.9750
KG/HR 6.2074+04
CUM/SEC 0.6172
STATE VARIABLES:
TEMP C 160.0000
PRES ATM 10.0000
VFRAC 1.0000
LFRAC 0.0
SFRAC 0.0
ENTHALPY:
References
MacAuley, K. B., MacGregor, J. F., & Hamielec, A. E. (1990). AIChE J., 36, p.
837.
Xie, T., McAuley, K. B., Hsu, J. C. C., & Bacon, D. W. (1994). Gas phase
ethylene polymerization: Production processes, polymer properties and
reactor modeling. Ind. Eng. Chem. Res., 33, pp. 449-479.
The LDPE high-pressure process model illustrates the use of Aspen Polymers
to model a high pressure free-radical process for low-density polyethylene.
The model considers the tubular jacketed reactor along with the post-
polymerization devolatization. Among the parameters studied in the model
are reactor temperature profile, monomer conversion, polymer molecular
weight and degree of branching.
Process Definition
In this example, polymerization of ethylene is carried out in a jacketed
tubular reactor divided into four sections, with two initiator injection points.
High and low pressure separators are also used. The flowsheet for the process
is shown here:
Reactors / Kinetics
The built-in free-radical reactions used in this example are:
Description k0 Ea V
(m3/kmol) (J/kmol) (m3/kmol)
Input Summary
The input language summary for this example is shown here:
DYNAPLUS
DPLUS RESULTS=ON
TITLE &
'High-Pressure LDPE tubular reactor flowsheet with four
sections'
SYS-OPTIONS TRACE=YES
RUN-CONTROL MAX-TIME=5000.0
DESCRIPTION "
This flowsheet models a high-pressure LDPE tubular
reactor flowsheet with four sections, two initiator
and monomer injection points, and high and low
pressure separators.
ASPCOMP USRHPE
"
FLOWSHEET
BLOCK PFR1 IN=FEED1 CIN1 OUT=OUT1 COUT1
BLOCK PFR2 IN=OUT1 CIN2 OUT=OUT2 COUT2
BLOCK FDMIX2 IN=OUT2 INIFD2 E2FD2B OUT=FEED2
BLOCK PFR3 IN=FEED2 CIN3 OUT=OUT3 COUT3
BLOCK PFR4 IN=OUT3 CIN4 OUT=OUT4 COUT4
BLOCK HPS IN=OUT4 OUT=HP-RECY HPSB
BLOCK LPS IN=HPSB OUT=LP-RECY LPSB
BLOCK FDMIX1 IN=E2FD1B INIFD1 OUT=FEED1
BLOCK FSPLIT1 IN=E2FD OUT=E2FD2A E2FD1A
BLOCK EFD1 IN=E2FD1A OUT=E2FD1B
BLOCK EFD2 IN=E2FD2A OUT=E2FD2B
PROPERTIES POLYSL
PROPERTIES STEAM-TA
PROP-DATA POLYSL
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=bar PDROP='N/sqm'
PROP-LIST SLTSTR / SLPSTR / SLRSTR
PVAL E2 333.0 / 2400.0 / 631.0
PVAL LDPE 667.70 / 3500.0 / 894.0
PVAL INI1 530 / 3040 / 1120
PVAL INI2 530 / 3040 / 1120
PVAL WATER 623.0 / 2.6871390E+04 / 1105.0
PROP-DATA TGTM
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=bar TEMPERATURE=C PDROP='N/sqm'
PROP-LIST TGVK / TMVK / CRITMW
PVAL LDPE -36.0 / 141.60 / 3500.0
PROP-DATA VLSTD
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=bar TEMPERATURE=C PDROP='N/sqm'
PROP-LIST VLSTD
PVAL E2-SEG 0.0
PROP-DATA CPIG-1
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=bar TEMPERATURE=C PDROP='N/sqm'
PROP-LIST CPIG
PVAL E2-SEG 35342.0 70.20
PVAL E2 42291.0 47.3550
PROP-DATA MULMH-1
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=bar TEMPERATURE=C PDROP='N/sqm'
PROP-DATA SLETIJ-1
IN-UNITS SI
PROP-LIST SLETIJ
BPVAL LDPE E2 -.15
BPVAL E2 LDPE -.15
PROP-DATA SLKIJ-1
IN-UNITS SI
PROP-LIST SLKIJ
BPVAL LDPE E2 .04
BPVAL E2 LDPE .04
POLYMERS
SEGMENTS E2-SEG REPEAT
POLYMERS LDPE
ATTRIBUTES LDPE DPN DPW PDI MWN MWW ZMOM FMOM SMOM &
SFLOW SFRAC FLCB FSCB LCB SCB LSFLOW LSFRAC LEFLOW &
LEFRAC
STREAM CIN1
SUBSTREAM MIXED TEMP=160.0 PRES=100.0 NPHASE=1 PHASE=L
MASS-FLOW WATER 160000.0
STREAM CIN2
SUBSTREAM MIXED TEMP=160.0 PRES=100.0 NPHASE=1 PHASE=L
MASS-FLOW WATER 160000.0
STREAM CIN3
SUBSTREAM MIXED TEMP=160.0 PRES=100.0 NPHASE=1 PHASE=L
MASS-FLOW WATER 160000.0
STREAM CIN4
SUBSTREAM MIXED TEMP=160.0 PRES=100.0 NPHASE=1 PHASE=L
MASS-FLOW WATER 160000.0
STREAM E2FD
SUBSTREAM MIXED TEMP=100. PRES=2020.
MASS-FLOW E2 65000.
STREAM E2FD1B
SUBSTREAM MIXED TEMP=140. PRES=2010.
MASS-FLOW E2 40000.
STREAM E2FD2B
SUBSTREAM MIXED TEMP=50. PRES=2000.
MASS-FLOW E2 25000.
STREAM FEED1
SUBSTREAM MIXED TEMP=140. PRES=2000.0 MASS-FLOW=63000.0 &
NPHASE=1 PHASE=L
MASS-FRAC INI1 2.2E-005 / E2 0.999976 / INI2 1.1E-005
STREAM INIFD1
SUBSTREAM MIXED TEMP=0. PRES=2010.
STREAM INIFD2
SUBSTREAM MIXED TEMP=0. PRES=2010. NPHASE=1 PHASE=L
MASS-FLOW INI1 12. / INI2 3.
EO-CONV-OPTI
CONV-OPTIONS
PARAM CHECKSEQ=NO
PROPERTY-REP NOPCES
REACTIONS R1 FREE-RAD
DESCRIPTION "Free-Radical Kinetic Scheme"
PARAM QSSA=YES
SPECIES INITIATOR=INI1 INI2 MONOMER=E2 POLYMER=LDPE
MON-RSEG E2 E2-SEG
INIT-DEC INI1 2.5000E+14 1.2600E+08 0.0 EFFIC=0.4 NRADS=2
INIT-DEC INI2 5.93E+018 1.9500E+08 0.0 EFFIC=0.4 NRADS=2
CHAIN-INI E2 2.5000E+08 3.5300E+07 0.0
PROPAGATION E2 E2 2.5000E+08 3.5300E+07 -.02130
CHAT-MON E2 E2 1250000.0 4.5400E+07 0.0
CHAT-POL E2 E2 1240000.0 3.0400E+07 .00160
B-SCISSION E2 6.0700E+07 4.5300E+07 0.0
TERM-DIS E2 E2 2.5000E+09 4190000.0 .0010
TERM-COMB E2 E2 2.5000E+09 4190000.0 .0010
SC-BRANCH E2 E2 1.3000E+09 4.1600E+07 0.0
Degree of Polymerization
3500
DPN
3000
DPW
2500
2000
1500
1000
500
0
0 50 100 150 200 250
Length, m
Initiator decomposition rate versus reactor length for PFR1 is shown here:
2.00E-05
INI1
Initiator Mole Fraction
1.50E-05 INI2
1.00E-05
5.00E-06
0.00E+00
0 50 100 150 200 250
Length, m
Coolant and reactant temperature profiles for PFR1 are shown here:
300
Temperature (C)
250
Process
200 Coolant
150
100
0 50 100 150 200 250
Length, m
The remaining ethylene mole fraction profile in PFR1 and PFR3 is shown here:
PFR1
0.95
Ethylene Mole Fraction
PFR3
0.9
0.85
0.8
0.75
0 50 100 150 200 250
Length, m
300
PFR1
Temperature (C)
PFR2
250
PFR3
PFR4
200
150
100
0 50 100 150 200 250
Length, m
SUBSTREAM: MIXED
PHASE: LIQUID LIQUID LIQUID LIQUID LIQUID
COMPONENTS: KMOL/HR
INI1 0.0 0.0 0.0 0.0 0.0
E2 0.0 0.0 0.0 0.0 0.0
LDPE 0.0 0.0 0.0 0.0 0.0
INI2 0.0 0.0 0.0 0.0 0.0
WATER 8881.3496 8881.3496 8881.3496 8881.3496 8881.3496
COMPONENTS: KG/HR
INI1 0.0 0.0 0.0 0.0 0.0
E2 0.0 0.0 0.0 0.0 0.0
LDPE 0.0 0.0 0.0 0.0 0.0
INI2 0.0 0.0 0.0 0.0 0.0
WATER 1.6000+05 1.6000+05 1.6000+05 1.6000+05 1.6000+05
COMPONENTS: MASS FRAC
INI1 0.0 0.0 0.0 0.0 0.0
E2 0.0 0.0 0.0 0.0 0.0
LDPE 0.0 0.0 0.0 0.0 0.0
INI2 0.0 0.0 0.0 0.0 0.0
WATER 1.0000 1.0000 1.0000 1.0000 1.0000
SUBSTREAM: MIXED
PHASE: LIQUID LIQUID LIQUID VAPOR VAPOR
COMPONENTS: KMOL/HR
INI1 0.0 0.0 0.0 0.0 0.0
E2 0.0 0.0 0.0 2316.9800 1425.8338
LDPE 0.0 0.0 0.0 0.0 0.0
INI2 0.0 0.0 0.0 0.0 0.0
WATER 8881.3496 8881.3496 8881.3496 0.0 0.0
COMPONENTS: KG/HR
INI1 0.0 0.0 0.0 0.0 0.0
E2 0.0 0.0 0.0 6.5000+04 4.0000+04
LDPE 0.0 0.0 0.0 0.0 0.0
INI2 0.0 0.0 0.0 0.0 0.0
WATER 1.6000+05 1.6000+05 1.6000+05 0.0 0.0
COMPONENTS: MASS FRAC
INI1 0.0 0.0 0.0 0.0 0.0
E2 0.0 0.0 0.0 1.0000 1.0000
LDPE 0.0 0.0 0.0 0.0 0.0
INI2 0.0 0.0 0.0 0.0 0.0
WATER 1.0000 1.0000 1.0000 0.0 0.0
TOTAL FLOW:
KMOL/HR 8881.3496 8881.3496 8881.3496 2316.9800 1425.8338
KG/HR 1.6000+05 1.6000+05 1.6000+05 6.5000+04 4.0000+04
CUM/SEC 4.9225-02 4.9225-02 4.9225-02 3.2683-02 2.0113-02
STATE VARIABLES:
TEMP C 170.0000 170.0000 170.0000 100.0000 100.0000
PRES BAR 96.0000 96.0000 96.0000 2020.0000 2020.0000
VFRAC 0.0 0.0 0.0 1.0000 1.0000
SUBSTREAM: MIXED
PHASE: LIQUID VAPOR LIQUID LIQUID LIQUID
COMPONENTS: KMOL/HR
INI1 0.0 0.0 0.0 2.4770-02 4.9539-02
E2 1425.8338 891.1461 891.1461 1425.8338 2147.5521
LDPE 0.0 0.0 0.0 0.0 169.6952
INI2 0.0 0.0 0.0 1.0258-02 2.0516-02
WATER 0.0 0.0 0.0 0.0 0.0
COMPONENTS: KG/HR
INI1 0.0 0.0 0.0 6.0000 12.0000
E2 4.0000+04 2.5000+04 2.5000+04 4.0000+04 6.0247+04
LDPE 0.0 0.0 0.0 0.0 4760.5891
INI2 0.0 0.0 0.0 1.5000 3.0000
WATER 0.0 0.0 0.0 0.0 0.0
COMPONENTS: MASS FRAC
INI1 0.0 0.0 0.0 1.4997-04 1.8455-04
E2 1.0000 1.0000 1.0000 0.9998 0.9266
LDPE 0.0 0.0 0.0 0.0 7.3214-02
INI2 0.0 0.0 0.0 3.7493-05 4.6138-05
WATER 0.0 0.0 0.0 0.0 0.0
TOTAL FLOW:
KMOL/HR 1425.8338 891.1461 891.1461 1425.8689 2317.3174
KG/HR 4.0000+04 2.5000+04 2.5000+04 4.0007+04 6.5023+04
CUM/SEC 2.0772-02 1.2570-02 1.2129-02 2.0803-02 3.3515-02
STATE VARIABLES:
TEMP C 140.0000 100.0000 50.0000 140.6072 145.4816
PRES BAR 2010.0000 2020.0000 2000.0000 2000.0000 1900.0000
VFRAC 0.0 1.0000 0.0 0.0 0.0
LFRAC 1.0000 0.0 1.0000 1.0000 1.0000
SFRAC 0.0 0.0 0.0 0.0 0.0
ENTHALPY:
J/KMOL 5.5884+07 5.3717+07 5.1043+07 5.5875+07 4.8595+07
J/KG 1.9920+06 1.9148+06 1.8195+06 1.9914+06 1.7319+06
WATT 2.2134+07 1.3297+07 1.2635+07 2.2131+07 3.1281+07
ENTROPY:
J/KMOL-K -1.1342+05 -1.1907+05 -1.2647+05 -1.1331+05 -1.1659+05
SUBSTREAM: MIXED
PHASE: VAPOR LIQUID LIQUID LIQUID VAPOR
COMPONENTS: KMOL/HR
INI1 0.0 0.0 2.4770-02 4.9539-02 0.0
E2 1829.6532 51.0738 0.0 0.0 50.9159
LDPE 0.0 437.0550 0.0 0.0 0.0
INI2 0.0 0.0 1.0258-02 2.0516-02 0.0
SUBSTREAM: MIXED
PHASE: LIQUID LIQUID LIQUID LIQUID LIQUID
COMPONENTS: KMOL/HR
INI1 0.0 0.0 0.0 0.0 0.0
E2 0.1579 1256.4059 1256.4059 1880.7270 1880.7270
LDPE 437.0550 169.6952 169.6952 437.0550 437.0550
INI2 0.0 0.0 0.0 0.0 0.0
WATER 0.0 0.0 0.0 0.0 0.0
COMPONENTS: KG/HR
INI1 0.0 0.0 0.0 0.0 0.0
E2 4.4300 3.5247+04 3.5247+04 5.2761+04 5.2761+04
LDPE 1.2261+04 4760.5891 4760.5891 1.2261+04 1.2261+04
INI2 0.0 0.0 0.0 0.0 0.0
WATER 0.0 0.0 0.0 0.0 0.0
COMPONENTS: MASS FRAC
INI1 0.0 0.0 0.0 0.0 0.0
E2 3.6118-04 0.8810 0.8810 0.8114 0.8114
LDPE 0.9996 0.1190 0.1190 0.1886 0.1886
INI2 0.0 0.0 0.0 0.0 0.0
WATER 0.0 0.0 0.0 0.0 0.0
TOTAL FLOW:
KMOL/HR 437.2129 1426.1012 1426.1012 2317.7820 2317.7820
KG/HR 1.2265+04 4.0008+04 4.0008+04 6.5023+04 6.5023+04
CUM/SEC 4.8136-03 2.1985-02 2.1512-02 3.6219-02 3.5501-02
STATE VARIABLES:
TEMP C 264.7096 243.8088 202.7850 278.5977 237.8896
PRES BAR 1.0000 1900.0000 1800.0000 1800.0000 1700.0000
VFRAC 0.0 0.0 0.0 0.0 0.0
LFRAC 1.0000 1.0000 1.0000 1.0000 1.0000
SFRAC 0.0 0.0 0.0 0.0 0.0
ENTHALPY:
J/KMOL -4.1755+07 4.9898+07 4.7084+07 4.4866+07 4.2011+07
J/KG -1.4884+06 1.7787+06 1.6783+06 1.5993+06 1.4975+06
WATT -5.0710+06 1.9767+07 1.8652+07 2.8886+07 2.7048+07
ENTROPY:
References
Brandup, J., & Immergut, E. H. (1989). Polymer Handbook, 3rd Ed. NewYork:
John Wiley.
Carr, N. L., Parent, J. D., & Peck, R. E. (1955). Viscosity of Gases and Gas
Mixtures at High Pressures. Chem. Eng. Progr. Symposium Ser. No. 16, 51,
91.
Danner, R. P., & High, M. S. (1982). Handbook of Polymer Solution
Thermodynamics. DIPPR.
Process Definition
In this example, the top portion of the VK column (zone 1) will be modeled as
two stirred tank reactors. The bottom portion of the column (zone 2 and
below) will be modeled as a plug flow reactor. The distillation column will be
included in the simulations as a multistage separator. A flash tank will be
inserted after each stirred tank in zone 1 to account for VLE calculations. The
plug flow reactor will be considered liquid filled.
Process Conditions
The process conditions are:
Components Name Databank Description
Feeds
Temperature 260
(C)
Pressure (atm) 1
Flowrate 40
(kg/hr)
Caprolactam 99.0% by weight
Water 1.0% by weight
Kinetics STEP-GROWTH model
Reaction Kinetics
The main reactions used in the step-growth kinetics are:
Description Reaction
1 Ring-opening of
caprolactam NH(CH2)5 CO + H2O H2N (CH2)5 COOH
2 Polycondensation
NH2 + HOOC NHCO + H2O
3 Polyaddition of
NH2 + CO (CH2)5 - NH NHCO - (CH2)5 - NH2
caprolactam
4 Ring-opeining of cyclic
dimer NH(CH2)5CONH(CH2)5CO + H2O H2N(CH2)5 CONH(CH2)5 COOH
5 Polyaddition of cyclic
dimer NH2 + NH(CH 2)5CONH(CH2)5CO NHCO(CH2)5NHCO(CH2 )5 NH2
The model generates the reaction network based on these functional groups.
To keep the model consistent with Tai and Tagawa, we have cleared the
reaction generation option for re-arrangement reactions on the Options
sheet, as shown here:
Input Summary
The input language summary for this example is shown here:
TITLE 'CAPROLACTAM POLYMERIZATION'
DIAGNOSTICS
HISTORY SIM-LEVEL=3 PROP-LEVEL=3
RUN-CONTROL MAX-TIME=15000.0
COMPONENTS
H2O H2O /
CL C6H11NO /
ACA C6H11NO /
CD C6H11NO /
NYLON NYLON6 /
T-NH2 C6H12NO-E-1 /
T-COOH C6H12NO2-E-1 /
B-ACA C6H11NO-R-1
FLOWSHEET
BLOCK STAGE1 IN=FEED RECYCLE R2-VAP OUT=R1-OLIGO R1VAP
BLOCK DISTIL IN=R1VAP OUT=COND RECYCLE
BLOCK STAGE2 IN=R1-OLIGO OUT=R2-VAP R2-OLIGO
BLOCK PLUG IN=R2-OLIGO OUT=POLYMER
PROPERTIES POLYNRTL
PROP-DATA MW
IN-UNITS SI
PROP-LIST MW
PVAL ACA 131.174760
PROP-DATA PLXANT-1
IN-UNITS SI PRESSURE=atm PDROP='N/sqm'
PROP-LIST PLXANT
PVAL ACA -40.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PVAL CD -40.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1000.0
PROP-DATA NRTL-1
IN-UNITS ENG
PROP-LIST NRTL
BPVAL H2O CL -4.01840 1634.2356470 .30 0.0 0.0 0.0 &
105.80000320 508.99999990
BPVAL CL H2O -.79760 3537.9406520 .30 0.0 0.0 0.0 &
105.80000320 508.99999990
BPVAL H2O ACA -4.01840 1634.2356470 .30 0.0 0.0 0.0 &
105.80000320 508.99999990
BPVAL ACA H2O -.79760 3537.9406520 .30 0.0 0.0 0.0 &
105.80000320 508.99999990
BPVAL H2O CD -4.01840 1634.2356470 .30 0.0 0.0 0.0 &
105.80000320 508.99999990
BPVAL CD H2O -.79760 3537.9406520 .30 0.0 0.0 0.0 &
105.80000320 508.99999990
POLYMERS
SEGMENTS T-NH2 END / T-COOH END / B-ACA REPEAT
POLYMERS NYLON
ATTRIBUTES NYLON SFLOW ZMOM FMOM DPN MWN EFRAC SFRAC
STREAM FEED
SUBSTREAM MIXED TEMP=260.0 PRES=1.0 MASS-FLOW=3600. &
NPHASE=1 PHASE=L
MASS-FRAC H2O .010 / CL .990
STREAM RECYCLE
SUBSTREAM MIXED TEMP=116.523631 PRES=1 MASS-FLOW=145.425133
MASS-FRAC H2O 0.217256374 / CL 0.782743626
EO-CONV-OPTI
SENSITIVITY FLOW
PARAM BASE-CASE=LAST
DEFINE R1CL MASS-FLOW STREAM=R1-OLIGO SUBSTREAM=MIXED &
COMPONENT=CL
DEFINE R1NYL MASS-FLOW STREAM=R1-OLIGO SUBSTREAM=MIXED &
COMPONENT=NYLON
DEFINE R1MN COMP-ATTR-VAR STREAM=R1-OLIGO SUBSTREAM=MIXED &
COMPONENT=NYLON ATTRIBUTE=MWN ELEMENT=1
DEFINE R2CL MASS-FLOW STREAM=R2-OLIGO SUBSTREAM=MIXED &
COMPONENT=CL
DEFINE R2NYL MASS-FLOW STREAM=R2-OLIGO SUBSTREAM=MIXED &
COMPONENT=NYLON
DEFINE R2MN COMP-ATTR-VAR STREAM=R2-OLIGO SUBSTREAM=MIXED &
COMPONENT=NYLON ATTRIBUTE=MWN ELEMENT=1
DEFINE R3MN COMP-ATTR-VAR STREAM=POLYMER SUBSTREAM=MIXED &
COMPONENT=NYLON ATTRIBUTE=MWN ELEMENT=1
DEFINE R3CL MASS-FLOW STREAM=POLYMER SUBSTREAM=MIXED &
COMPONENT=CL
DEFINE R3NYL MASS-FLOW STREAM=POLYMER SUBSTREAM=MIXED &
COMPONENT=NYLON
DEFINE R3ACA MASS-FLOW STREAM=POLYMER SUBSTREAM=MIXED &
COMPONENT=ACA
DEFINE R3CD MASS-FLOW STREAM=POLYMER SUBSTREAM=MIXED &
COMPONENT=CD
DEFINE R3F STREAM-VAR STREAM=FEED SUBSTREAM=MIXED &
VARIABLE=MASS-FLOW
DEFINE R1COND MASS-FLOW STREAM=COND SUBSTREAM=MIXED &
COMPONENT=H2O
DEFINE R3FW MASS-FRAC STREAM=R2-OLIGO SUBSTREAM=MIXED &
COMPONENT=H2O
DEFINE R2COND MASS-FLOW STREAM=R2-VAP SUBSTREAM=MIXED &
COMPONENT=H2O
DEFINE R2FW MASS-FRAC STREAM=R1-OLIGO SUBSTREAM=MIXED &
COMPONENT=H2O
DEFINE FWF MASS-FRAC STREAM=FEED SUBSTREAM=MIXED &
CONV-OPTIONS
PARAM TRACE=.00010 CHECKSEQ=YES
WEGSTEIN MAXIT=50
Stage 1
0.82
Stage 2
Product
As flow rate increases, the reactor residence time decreases. The average
molecular weight of the polymer at each point in the reactor decreases with
the residence time.
22000
Number Average Molecular Weight (Mn)
Stage 1
Stage 2
20000
Product
16000 18000
Extraction value (EV) is the sum of the extractables from the polymer. These
include unreacted caprolactam, aminocaproic acid, cyclic dimer and water.
Extraction value is an indicator of the cost of recycling unreacted monomers.
Since conversion decreases with an increase in feed flow rate, EV increases.
Therefore, one has to find the optimum feed flow rate that increases the
production rate while minimizing the percentage of extractables.
The effect of feed flow rate on extraction value is shown here:
SUBSTREAM: MIXED
PHASE: LIQUID LIQUID LIQUID LIQUID VAPOR
COMPONENTS: KMOL/HR
H2O 1.4138 1.9983 0.4199 0.4074 3.1676
CL 0.0 31.4954 2.5820 4.3448 1.0059
ACA 0.0 0.0 9.3725-04 1.3062-03 0.0
CD 0.0 0.0 0.1336 9.1208-02 0.0
NYLON 0.0 0.0 28.6715 26.9949 0.0
COMPONENTS: KG/HR
H2O 25.4704 36.0000 7.5647 7.3389 57.0643
CL 0.0 3564.0000 292.1742 491.6540 113.8274
ACA 0.0 0.0 0.1229 0.1713 0.0
CD 0.0 0.0 30.2269 20.6420 0.0
NYLON 0.0 0.0 3244.4503 3054.7330 0.0
COMPONENTS: MASS FRAC
H2O 1.0000 1.0000-02 2.1163-03 2.0531-03 0.3339
CL 0.0 0.9900 8.1738-02 0.1375 0.6661
ACA 0.0 0.0 3.4394-05 4.7932-05 0.0
CD 0.0 0.0 8.4562-03 5.7747-03 0.0
NYLON 0.0 0.0 0.9077 0.8546 0.0
TOTAL FLOW:
KMOL/HR 1.4138 33.4937 31.8078 31.8396 4.1735
KG/HR 25.4704 3600.0000 3574.5391 3574.5392 170.8917
L/HR 26.4708 4128.3629 4076.1662 4082.4945 1.8156+05
Process Definition
The batch reactor is modeled using RBatch, and the devolitizer using Flash2:
Feed Streams
Temperature (C) 0
Pressure (atm) 1
Methylmethacrylate 100 (25 wt%)
(kg/h)
Ethyl Acetate (kg/hr) 300 (75 wt%)
Initiator (kg/hr) 3.757
The feed is sent to RBATCH as a batch feed
Kinetics FREE-RAD model
Temperature (C) 50 60
Pressure 1 atm 100 mmHg
Process Studies
The PMMA process model is used to examine process variables such as
conversion number and weight average molecular weight as a function of
time, at three different reactor temperatures. The model predictions are
compared against literature data (Ellis, 1990).
The model is also used to study the effects of initiator concentration and
polymerization temperature on the following: reaction time to reach 99%
monomer conversion, peak heat load, number and weight average molecular
weight.
A sensitivity study is also performed to show the percent solids as a function
of outlet temperature.
Input Summary
The input language summary for this example is shown here:
DYNAPLUS
DPLUS RESULTS=ON
IN-UNITS MET
DESCRIPTION "
Polymerization of methyl methacrylate(MMA) in a batch.
reactor "
COMPONENTS
MMA C5H8O2-D3 /
EA C4H8O2-3 /
AIBN C5H8O2-D3 /
PMMA PMMA /
MMA-R C5H8O2-R-1
FLOWSHEET
BLOCK FLASH IN=PROD OUT=FV FL
BLOCK RBATCH IN=FEED OUT=PROD
PROPERTIES POLYNRTL
PROP-DATA MW
IN-UNITS MET
PROP-LIST MW
PVAL AIBN 164.210
POLYMERS
SEGMENTS MMA-R REPEAT
POLYMERS PMMA
ATTRIBUTES PMMA SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM &
FMOM SMOM LDPN LZMOM LFMOM LSFLOW LSFRAC LEFLOW &
LEFRAC LPFRAC
DISTRIBUTION PMMA CHAIN-SIZE CLD NPOINTS=100 UPPER=10000.
STREAM FEED
SUBSTREAM MIXED TEMP=273.0 PRES=1.0
MASS-FLOW MMA 100.0 / EA 300.0 / AIBN 3.7570
EO-CONV-OPTI
CALCULATOR F-1
DEFINE FTEMP STREAM-VAR STREAM=FEED SUBSTREAM=MIXED &
VARIABLE=TEMP
DEFINE RTEMP BLOCK-VAR BLOCK=RBATCH VARIABLE=TEMP &
SENTENCE=PARAM
F FTEMP = RTEMP
EXECUTE BEFORE BLOCK RBATCH
CONV-OPTIONS
PARAM CHECKSEQ=NO
100
T = 323K
80 T = 333K
PMMA Mass (Kg)
T = 338K
60 T = 343K
40
20
0
0 1 2 3 4 5
Time (Hr)
70000
60000
T = 323K
50000
T = 333K
PMMA MWN
T = 338K
40000
T = 343K
30000
20000
10000
0
0 1 2 3 4 5
Time (Hr)
50
45 AIBN Charge = 3
Time for 99% Conversion
40
AIBN Charge = 3.757
35
AIBN Charge = 6
30
AIBN Charge = 10
25
20
15
10
5
0
320 325 330 335 340 345 350 355
Temperature (K)
25000
AIBN Charge = 3
20000 AIBN Charge = 3.757
MWN at 99% Conversion
AIBN Charge = 6
15000 AIBN Charge = 10
10000
5000
0
320 325 330 335 340 345 350 355
Temperature (K)
4.4
4
AIBN Charge = 6
3.8 AIBN Charge = 10
3.6
3.4
3.2
2.8
320 325 330 335 340 345 350 355
Temperature (K)
-500
Max. Duty (cal/s)
-1000
For the case of a reactor with specified temperature, the maximum heat duty
occurs immediately after reactions start. In the simulation results presented
above, the heat duty reported at 0.5 hr is taken to be the maximum. This is
the second profile point computed when the time interval is set to 0.5 hour.
Technically it is not the true peak heat duty. Smaller time intervals between
profile points may be used to obtain more accurate peak heat duty.
The results of the sensitivity study on the flash block are shown here:
0.96
0.94
PMMA Mass Fraction
0.92
0.9
0.88
0.86
0.84
320 325 330 335 340 345 350 355
Flash Temperature (K)
As the flash temperature varies, the polymer content in the liquid output
stream changes accordingly. No binary parameters are used in the polymer
NRTL model. If liquid-vapor equilibria data is available, the binary parameters
may be fitted to allow more accurate predictions.
Process Definition
An Aspen Polymers model is developed to simulate a PP gas-phase UNIPOL
flowsheet. The flowsheet includes the fluidized bed reactor, the gas
recycle/cooling loop, discharge, and purge units:
Process Conditions
The process conditions are:
Components
Temperature 69 C
Pressure drop 200 KPa
Total volume 90 m3
Polymer Phase volume 60 m3
The reactor feed also includes a large recycle stream with partially condensed
propylene and other volatile components.
An Aspen Plus CSTR reactor model is used to represent the fluid-bed reactor.
The reactor is considered to have two phases; a vapor phase and a polymer
phase.
Reactors / Kinetics
The Ziegler-Natta kinetic scheme used in the model includes spontaneous site
activation, chain initiation, propagation, chain transfer to hydrogen, monomer
and co-catalyst, and spontaneous site deactivation. Four site types are used
in this example to represent the broad molecular weight distributions that are
typically observed for Ziegler-Natta polymers. The actual number of site types
necessary to model a given catalyst-polymerization system is determined by
deconvolution of the polymer molecular weight distribution curve obtained
from a GPC analysis.
The rate parameters used in the model are:
Reaction Site ID Comp. ID Comp. ID Pre- Activation Order Ref.
1 2 exp Energy Temp.
(kcal/mol) (C)
DEACT-SPON 3 0.001 1 1 69
DEACT-SPON 4 0.001 1 1 69
ATACT-PROP 1 C3H6 C3H6 13.2 7.2 1 69
ATACT-PROP 2 C3H6 C3H6 31.8 7.2 1 69
ATACT-PROP 3 C3H6 C3H6 8 7.2 1 69
ATACT-PROP 4 C3H6 C3H6 18.5 7.2 1 69
The units for the frequency factor and activation energy are in SI units.
Process Studies
Input Summary
The input language summary for this example is shown here:
TITLE 'Ziegler-Natta Gas-Phase Polymerization of Propylene'
SIM-OPTIONS
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
VOLUME-FLOW='cum/hr' PRESSURE=kPa TEMPERATURE=C TIME=hr &
PDROP=kPa
SIM-OPTIONS FLASH-TOL=1E-006
SYS-OPTIONS TRACE=YES
RUN-CONTROL MAX-TIME=10000.0
DESCRIPTION "
Ziegler-Natta gas-Phase polymerization of propylene
"
COMPONENTS
TICL4 TICL4 /
TEA C6H15AL /
C3H6 C3H6-2 /
C3H8 C3H8 /
PP PP /
H2 H2 /
N2 N2 /
FLOWSHEET
BLOCK REACTOR IN=RFEED1 OUT=VAP-A POWDER1
BLOCK FMIX1 IN=CAT COCAT H2FEED C3FEED N2FEED CYCGAGB &
OUT=RFEED1
BLOCK EXC1 IN=H2O-IN VAP-B OUT=H2O-OUT CYCGAGB
BLOCK COMP1 IN=VAP-A OUT=VAP-B
BLOCK STRIP1 IN=POWDER1 OUT=GAS1 POWDER2
BLOCK STRIP2 IN=POWDER2 STRIP-N2 OUT=GAS2 PP-PROD
PROPERTIES POLYPCSF
PROP-DATA PCSAFT
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
TEMPERATURE=C TIME=hr DELTA-T=C
PROP-LIST PCSFTM
PVAL TICL4 20
PVAL TEA 20
PROP-LIST PCSFTV
PVAL C3H6-R 4.1473
POLYMERS
PARAM NSITES=4
SEGMENTS C3H6-R REPEAT
CATALYST TICL4 0.00045
POLYMERS PP
ATTRIBUTES PP SFRAC SFLOW DPN DPW PDI MWN MWW ZMOM FMOM &
SMOM LDPN LZMOM LFMOM LSFLOW LSFRAC LEFLOW LEFRAC &
LPFRAC SSFRAC SSFLOW SDPN SDPW SPDI SMWN SMWW SZMOM &
SFMOM SSMOM SPFRAC LSDPN LSZMOM LSFMOM LSSFLOW LSSFRAC &
LSEFLOW LSEFRAC LSPFRAC ATFLOW ATFRAC SATFLOW SATFRAC
ATTRIBUTES TICL4 CPSFLOW CPSFRAC CVSFLOW CVSFRAC CDSFLOW &
CDSFRAC
DISTRIBUTION PP
STREAM C3FEED
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
VOLUME-FLOW='cum/hr' PRESSURE=kPa TEMPERATURE=C TIME=hr &
PDROP=kPa
SUBSTREAM MIXED TEMP=30.0 PRES=3000.0 MASS-FLOW=45000.0
MASS-FRAC C3H6 .9980 / C3H8 .0020
STREAM CAT
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
VOLUME-FLOW='cum/hr' PRESSURE=kPa TEMPERATURE=C TIME=hr &
PDROP=kPa
SUBSTREAM MIXED TEMP=30.0 PRES=3000.0 MASS-FLOW=300.0 &
NPHASE=1 PHASE=L
MASS-FRAC TICL4 .010 / C3H6 .9880 / C3H8 .0020
COMP-ATTR TICL4 CPSFLOW ( 0.0 )
COMP-ATTR TICL4 CPSFRAC ( 1.0 )
COMP-ATTR TICL4 CVSFLOW ( 0.0 0.0 0.0 0.0 )
COMP-ATTR TICL4 CVSFRAC ( 0.0 0.0 0.0 0.0 )
COMP-ATTR TICL4 CDSFLOW ( 0.0 )
COMP-ATTR TICL4 CDSFRAC ( 0.0 )
STREAM CYCGAGB
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
VOLUME-FLOW='cum/hr' PRESSURE=kPa TEMPERATURE=C TIME=hr &
PDROP=kPa
SUBSTREAM MIXED TEMP=59.0 PRES=3000.0
MASS-FLOW C3H6 800300.0 / C3H8 15930.0 / H2 3998.70 / &
N2 5098.0
STREAM H2FEED
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
VOLUME-FLOW='cum/hr' PRESSURE=kPa TEMPERATURE=C TIME=hr &
PDROP=kPa
SUBSTREAM MIXED TEMP=30.0 PRES=3000.0 MASS-FLOW=4.0
MASS-FRAC H2 1.0
STREAM H2O-IN
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=kPa TEMPERATURE=C TIME=hr DELTA-T=C PDROP=kPa
SUBSTREAM MIXED TEMP=35.0 PRES=3000.0
MASS-FLOW H2O 847500.0
STREAM N2FEED
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
VOLUME-FLOW='cum/hr' PRESSURE=kPa TEMPERATURE=C TIME=hr &
PDROP=kPa
SUBSTREAM MIXED TEMP=30.0 PRES=3000.0 MASS-FLOW=100.0
MASS-FRAC N2 1.0
STREAM RFEED1
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=kPa TEMPERATURE=C TIME=hr DELTA-T=C PDROP=kPa
SUBSTREAM MIXED TEMP=65.0 PRES=3000.0
MASS-FLOW TICL4 3.0 / TEA 10.0 / C3H6 817500.0 / C3H8 &
16100.0 / H2 4000.0 / N2 5100.0
COMP-ATTR TICL4 CPSFLOW ( 0.0 )
COMP-ATTR TICL4 CPSFRAC ( 1.0 )
COMP-ATTR TICL4 CVSFLOW ( 0.0 0.0 0.0 0.0 )
COMP-ATTR TICL4 CVSFRAC ( 0.0 0.0 0.0 0.0 )
COMP-ATTR TICL4 CDSFLOW ( 0.0 )
COMP-ATTR TICL4 CDSFRAC ( 0.0 )
STREAM STRIP-N2
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=kPa TEMPERATURE=C TIME=hr DELTA-T=C PDROP=kPa
SUBSTREAM MIXED TEMP=30.0 PRES=300.0
MASS-FLOW N2 500.0
EO-CONV-OPTI
CALCULATOR SETFEED
IN-UNITS SI MASS-FLOW='kg/hr' MOLE-FLOW='kmol/hr' &
PRESSURE=kPa TEMPERATURE=C TIME=hr DELTA-T=C PDROP=kPa
DEFINE RC3E MASS-FLOW STREAM=CYCGAGB SUBSTREAM=MIXED &
COMPONENT=C3H6
DEFINE RC3 MASS-FLOW STREAM=CYCGAGB SUBSTREAM=MIXED &
COMPONENT=C3H8
DEFINE RH2 MASS-FLOW STREAM=CYCGAGB SUBSTREAM=MIXED &
COMPONENT=H2
DEFINE RN2 MASS-FLOW STREAM=CYCGAGB SUBSTREAM=MIXED &
COMPONENT=N2
DEFINE FC3E MASS-FLOW STREAM=C3FEED SUBSTREAM=MIXED &
COMPONENT=C3H6
DEFINE FC3 MASS-FLOW STREAM=C3FEED SUBSTREAM=MIXED &
COMPONENT=C3H8
CONV-OPTIONS
PARAM CHECKSEQ=NO
SUBSTREAM: MIXED
PHASE: LIQUID LIQUID LIQUID MIXED VAPOR
COMPONENTS: KMOL/HR
TICL4 0.0 1.5815-02 0.0 0.0 0.0
TEA 0.0 0.0 8.7591-02 0.0 0.0
C3H6 481.8696 7.0436 0.0 1.8938+04 93.4567
C3H8 2.0007 1.3607-02 0.0 363.1075 1.8182
PP 0.0 0.0 0.0 0.0 0.0
H2 0.0 0.0 0.0 1983.2900 0.2875
N2 0.0 0.0 0.0 181.9875 6.7409-02
H2O 0.0 0.0 0.0 0.0 0.0
COMPONENTS: KG/HR
SUBSTREAM: MIXED
PHASE: VAPOR VAPOR LIQUID LIQUID VAPOR
COMPONENTS: KMOL/HR
TICL4 0.0 0.0 0.0 0.0 0.0
TEA 0.0 0.0 0.0 0.0 0.0
C3H6 9.3084 0.0 0.0 0.0 0.0
C3H8 0.1833 0.0 0.0 0.0 0.0
PP 0.0 0.0 0.0 0.0 0.0
H2 4.1346-04 0.9551 0.0 0.0 0.0
N2 17.8059 0.0 0.0 0.0 6.7702-02
H2O 0.0 0.0 4.7043+04 4.7043+04 0.0
COMPONENTS: KG/HR
TICL4 0.0 0.0 0.0 0.0 0.0
TEA 0.0 0.0 0.0 0.0 0.0
C3H6 391.7052 0.0 0.0 0.0 0.0
C3H8 8.0813 0.0 0.0 0.0 0.0
PP 0.0 0.0 0.0 0.0 0.0
H2 8.3348-04 1.9253 0.0 0.0 0.0
N2 498.8065 0.0 0.0 0.0 1.8966
H2O 0.0 0.0 8.4750+05 8.4750+05 0.0
COMPONENTS: MASS FRAC
TICL4 0.0 0.0 0.0 0.0 0.0
TEA 0.0 0.0 0.0 0.0 0.0
C3H6 0.4359 0.0 0.0 0.0 0.0
C3H8 8.9932-03 0.0 0.0 0.0 0.0
PP 0.0 0.0 0.0 0.0 0.0
H2 9.2754-07 1.0000 0.0 0.0 0.0
N2 0.5551 0.0 0.0 0.0 1.0000
H2O 0.0 0.0 1.0000 1.0000 0.0
TOTAL FLOW:
KMOL/HR 27.2981 0.9551 4.7043+04 4.7043+04 6.7702-02
KG/HR 898.5939 1.9253 8.4750+05 8.4750+05 1.8966
CUM/SEC 0.2096 2.2702-04 0.2391 0.2412 1.5692-05
STATE VARIABLES:
TEMP C 60.0000 30.0000 35.0000 44.0000 30.0000
PRES N/SQM 1.0000+05 3.0000+06 3.0000+06 3.0000+06 3.0000+06
VFRAC 1.0000 1.0000 0.0 0.0 1.0000
LFRAC 0.0 0.0 1.0000 1.0000 0.0
SFRAC 0.0 0.0 0.0 0.0 0.0
ENTHALPY:
J/KMOL 7.4863+06 1.8496+05 -2.8683+08 -2.8602+08 -4.6019+04
J/KG 2.2742+05 9.1750+04 -1.5921+07 -1.5877+07 -1642.7572
WATT 5.6767+04 49.0676 -3.7481+09 -3.7376+09 -0.8654
ENTROPY:
J/KMOL-K -3.9913+04 -2.7710+04 -1.6580+05 -1.6323+05 -2.8251+04
J/KG-K -1212.4972 -1.3746+04 -9203.5570 -9060.3732 -1008.4905
DENSITY:
KMOL/CUM 3.6182-02 1.1686 54.6449 54.1828 1.1985
KG/CUM 1.1910 2.3557 984.4427 976.1187 33.5736
VAP-A VAP-B
SUBSTREAM: MIXED
PHASE: VAPOR VAPOR
COMPONENTS: KMOL/HR
TICL4 0.0 0.0
TEA 0.0 0.0
C3H6 1.8938+04 1.8938+04
C3H8 363.1075 363.1075
PP 0.0 0.0
H2 1983.2900 1983.2900
N2 181.9875 181.9875
H2O 0.0 0.0
COMPONENTS: KG/HR
TICL4 0.0 0.0
TEA 0.0 0.0
C3H6 7.9692+05 7.9692+05
C3H8 1.6012+04 1.6012+04
PP 0.0 0.0
H2 3998.0747 3998.0747
N2 5098.1034 5098.1034
H2O 0.0 0.0
COMPONENTS: MASS FRAC
TICL4 0.0 0.0
TEA 0.0 0.0
C3H6 0.9695 0.9695
C3H8 1.9478-02 1.9478-02
PP 0.0 0.0
H2 4.8637-03 4.8637-03
N2 6.2018-03 6.2018-03
H2O 0.0 0.0
TOTAL FLOW:
KMOL/HR 2.1466+04 2.1466+04
KG/HR 8.2203+05 8.2203+05
CUM/SEC 4.6162 4.2580
STATE VARIABLES:
TEMP C 69.0000 74.8156
PRES N/SQM 2.8000+06 3.0500+06
VFRAC 1.0000 1.0000
LFRAC 0.0 0.0
SFRAC 0.0 0.0
ENTHALPY:
J/KMOL 1.6384+07 1.6641+07
J/KG 4.2784+05 4.3456+05
WATT 9.7693+07 9.9228+07
ENTROPY:
J/KMOL-K -1.4945+05 -1.4924+05
J/KG-K -3902.5890 -3897.1742
DENSITY:
KMOL/CUM 1.2917 1.4004
KG/CUM 49.4656 53.6268
AVG MW 38.2939 38.2939
248 Index
entering pure component input summary 90
parameters 57 EPS See expanded polystyrene
predicting properties for 2939 equations of state
specifying 9, 31 fitting parameters 4051
convergence Sanchez-Lacombe 40
for blocks 61, 65 ethanol 70
improving 17 ethyl acetate polymerization 221
crystallinity 34 31
custom property parameters 11 ethyl acrylate copolymerization
13044
D ethyl benzene 95, 106
ethylene
data conversion 174
experimental 44 HDPE polymerization 16583
fitting 5269 LDPE polymerization 184202
regression 4051 ethylene-polyethylene 40
sets 62 examples
data fit calculating end-use properties
folder 62 8392
specifying regression cases 64 expanded polystyrene
data regression suspension polymerization
entering experimental data 44 12029
example for property parameters fitting kinetic parameters 5269
4051 fractionating oligomers 7082
input summary 49 HDPE polymerization 16583
results 66 LDPE polymerization 184202
specifying 62 methyl methacrylate
specifying regression cases 46 polymerization 22131
density nylon 6 caprolactam
predicting 29 polymerization 20221
descriptions for simulations 9 polypropylene gas phase
distillation of polystyrene 10619 polymerization 23250
distributions polystyrene bulk polymerization
for flowsheets 22 93105
for reactors 22 polystyrene with styrene
for streams 22 monomer distillation 10619
linking with Excel 27 predicting physical properties
live plots 23 2939
plots for multiple simulations 25 regressing property parameters
plotting 22 4051
tables 22 styrene butadiene emulsion
tracking 11 copolymerization 14456
di-tert-butyl peroxide 95 styrene butadiene ionic
dodecyl mercaptan 95 polymerization 15764
DRS See data regression styrene ethyl acrylate
copolymerization 13044
E Excel
linking 27
emulsifiers 145
expanded polystyrene
emulsion polymerization of styrene
conversion effects 128
butadiene 14456
example flowsheet 121
end-use properties
example input summary 124
calculating 8392
example kinetics 123
Index 249
example process conditions 122 G
example process studies 124
gas-phase polymerization 23250
example results 128
gel effect specification 15
free-radical polymerization 120
global options
number average molecular
defining unit-sets 8
weight 128
overriding values 17
polymerization 12029
report options 9
experimental data 44
specifying 8
group contribution method 29
F
feed flow rate H
effect on ethylene conversion
HDPE
174
conversion results 174
effect on extraction value 218
example flowsheet 166
effect on hydrogen concentration
example input summary 168
175
example kinetics 167
effect on melt index 176
example model 16583
effect on Mn 174, 217
example process conditions 166
effect on Mw 175
example process studies 168
effect on number average
example results 174
molecular weight 101
hydrogen concentration results
effect on polydispersity 102, 176
175, 177
effect on recycle rate 102
melt index results 176
effect on styrene conversion 100
Mn results 174
effect on weight average
molecular weight distribution 178
molecular weight 101
Mw results 175
feed streams specification 12
polydispersity results 176
fitting
stream summary 178
input summary example 67
heat capacity
flory-Huggins model 70
predicting 29
flow rate See feed flow rate
high-density polyethylene See
flowsheets
HDPE
creating 6
hydrogen 166
distribution data 22
hydrogen concentration 175, 177
placing blocks 6
placing streams 7
renaming blocks 7 I
renaming streams 7 inhibitors 109
Fortran subroutines 83, 90 initiation
fractionation results 66
example 7082 thermal 93105
input summary 80 initiators
free-radical polymerization ammonium persulfate 145
of ethylene 184202 decomposition vs reactor length
of expanded polystyrene 12029 193
of methyl methacrylate 22131 methyl peroxide 52
of polystyrene 93105 input summary
of styrene butadiene 14456 for data regression example 49
of styrene ethyl acrylate 13044 for end-use properties example
90
for expanded polystyrene
polymerization 124
250 Index
for fitting kinetic parameters reactant temperature profiles
example 67 193
for HDPE polymerization 168 reactor temperature profiles 194
for LDPE polymerization 188 stream summary 195
for methyl methacrylate live plots 23
polymerization 224 low-density polyethylene See LDPE
for nylon 6 caprolactam
polymerization 212 M
for oligomer fractionation 80
for polypropylene polymerization mass fraction crystallinity 34
237 melt index 83
for polystyrene bulk methanol 52
polymerization 96 methyl methacrylate
for polystyrene with styrene example flowsheet 223
monomer distillation 109 example input summary 224
for predicting physical properties example kinetics 224
example 38 example process conditions 223
for styrene ethyl acrylate example process studies 224
copolymerization 136 example results 226
for styrene-butadiene emulsion flow vs time 226
copolymerization 150 free-radical polymerization 222
for styrene-butadiene ionic heat duty 229
polymerization 160 MWn vs temperature 228
intrinsic viscosity 83 MWw vs temperature 229
ionic polymerization number-average molecular
of styrene butadiene 15764 weight vs time 227
polydispersity 229
polymerization in ethyl acetate
K
22131
kinetics reaction time for monomer
fitting parameters 5269 conversion 227
modifying reactions 13 methyl peroxide 52
polymerization 58 mixers 18
specifying 13 mixtures
specifying gel effect 15 ethylene-polyethylene 41
molecular structure 33
L molecular weight distribution
HDPE example 178
LDPE
coolant temperature profiles 193
N
DPn and DPw vs reactor length
192 n-hexane 166
ethylene mole fraction 194 n-pentane 121
example flowsheet 185 nylon 6
example input summary 188 caprolactam conversion 217
example kinetics 187 caprolactam polymerization 202
example model 184202 21
example process conditions 186 example flowsheet 204
example process studies 188 example input summary 212
example results 192 example kinetics 206
free-radical polymerization 184 example process conditions 205
initiator decomposition vs reactor example process studies 212
length 193 example results 217
Index 251
example stream summary 218 of polystyrene by thermal
extraction value results 218 initiation 93105
Mn results 217 of styrene butadiene 14456,
step-growth polymerization 203 15764
of styrene ethyl acrylate 13044
O solution 16583, 22131
suspension 12029, 23250
OLE linking 27 polymers
oligomers 2,6 dimethyl phenol 30
characterizing 10 characterizing 10
fractionation example 7082 end-use properties 8392
ethylene-polyethylene 40
P fractionation example 7082
polystyrene 70
parameters
predicting formation during
binary interaction 40
styrene distillation 10619
entering pure component 57
specifying mass fraction
fitting kinetic 5269
drystallinity 34
fitting mixture 4051
vinyl acetate 52
Flory-Huggins 70
polymethyl methacrylate See
for polypropylene example 236
methyl methacrylate
regression for styrene ethyl
polyolefins 41
acrylate example 133
polypropylene
specifying 47
example flowsheet 233
specifying binary interaction 76
example input summary 237
specifying custom 11
example kinetics 236
physical properties
example parameters 236
predicting parameters 2939
example process conditions 234
regressing parameters 4051
example results 242
plots
example stream summary 242
distributions 22
gas phase polymerization 232
live plots 23
50
manually 25
process types 232
multiple simulations 25
polystyrene See also expanded
poly(phenylene oxide) 30
polystyrene
polydispersity 70
bulk polymerization 93105
polyethylene
example flowsheet 94
HDPE polymerization 16583
example input summary 96
intrinsic viscosity 83
example kinetics 96
LDPE polymerization 184202
example process conditions 94
melt index 83
example process studies 96
zero-shear viscosity 84
example results 100
Ziegler-Natta polymerization 165
example stream summary 103
polymerization
fractionation 70
emulsion 14456
fractionation results 79
gas-phase 23250
free-radical polymerization 93
ionic 15764
number average molecular
of expanded polystyrene 12029
weight 101
of HDPE 16583
polydispersity 70, 102
of LDPE 184202
recycle rate 102
of methyl methacrylate 22131
styrene monomer distillation
of nylon 6 caprolactam 20221
10619
of polypropylene 23250
252 Index
weight average molecular weight cases 46, 64
101 for styrene ethyl acrylate
polystyrene distillation example 133
example flowsheet 106 reports
example input summary 109 defining options 9
example kinetics 109 displaying a prop-set 58
example process conditions 107 results
example results 114 chain length distribution (SBR)
example stream summary 116 155
inhibitor 109 for expanded polystyrene
polymer production results 115 polymerization 128
profitability results 116 for HDPE polymerization 174
sensitivity study 115 for initiation 66
process models See simulations for LDPE polymerization 192
propagation results 66 for methyl methacrylate
properties polymerization 226
creating property sets 34 for nylon 6 caprolactam
creating property tables 35 polymerization 217
end-use 8392 for polymer heat capacity 36
estimating 30 for polypropylene polymerization
predicting 29 242
specifying calculations 77 for polystyrene bulk
specifying custom parameters 11 polymerization 100
specifying parameters 47, 76 for polystyrene with styrene
property distillation 114
methods 11, 43 for propagation 66
property analysis 30 for styrene ethyl acrylate
property method copolymerization 140
POLYSL 43 for styrene-butadiene emulsion
property parameters copolymerization 153
pure components 47 for styrene-butadiene ionic
property sets 34 polymerization 162
property tables 35 fractionation of polystyrene 79
prop-sets plotting 22
defining 58 reviewing 18
Fortran subroutines 83 saving 28
pseudo-component approach 70 run types
pure components data regression 41
constant 75 property analysis 31
parameters 57 runs
predicting properties for 2939 new 5
starting 18
R
S
rate constants
example input screen 14 Sanchez-Lacombe 40
for styrene-butadiene emulsion SBR See styrene butadiene
example 149 segments
reactions characterizing 10
modifying 13 defining molecular structure 33
rate constants 14 sensitivity tables 88
reactor distribution data 22 setup options
regression adding descriptions 9
Index 253
adding titles 8 solution polymerization 16583,
defining unit-sets 8 22131
report options 9 solvents
specifying 8 benzene 70, 106
specifying time limits 9 methanol 52
simulations n-hexane 166
calculating end-use properties toluene 106
8392 species characterization 10
creating flowsheets 6 sstyrene distillation
creating models 528 example stream summary 116
creating new runs 5 step-growth polymerization of
defining report options 9 caprolactam 20221
defining unit-sets 8 stream summary
entering descriptions for 9 for HDPE polymerization 178
entering time limits 9 for LDPE polymerization 195
entering titles for 8 for nylon 6 caprolactam
examples 592 polymerization 218
expanded polystyrene for polypropylene polymerization
suspension polymerization 242
12029 for polystyrene bulk
fitting kinetic parameters 5269 polymerization 103
fractionating oligomers 7082 for polystyrene with styrene
HDPE polymerization 16583 distillation 116
LDPE polymerization 184202 streams
methyl methacrylate distribution data 22
polymerization 22131 entering conditions 61
nylon 6 caprolactam in flowsheets 7
polymerization 20221 renaming 7
plotting multiple 25 report options 58
polypropylene gas phase specifying 12
polymerization 23250 styrene
polystyrene bulk polymerization bulk polymerization 93105
93105 butadiene emulsion
polystyrene with styrene copolymerization 14456
monomer distillation 10619 butadiene ionic polymerization
predicting physical properties 15764
2939 conversion results 100
regressing property parameters ethyl acrylate copolymerization
4051 13044
results 18 example flowsheet 94
running 18, 36 example input summary 96
saving 28 example kinetics 96
styrene butadiene emulsion example process conditions 94
copolymerization 14456 example process studies 96
styrene butadiene ionic example results 100
polymerization 15764 example stream summary 103
styrene ethyl acrylate free-radical polymerization 93
copolymerization 13044 monomer distillation of
templates 6 polystyrene 10619
site-based species characterization number average molecular
10 weight 101
sodium lauryl sulfate 145 polydispersity 102
recycle rate 102
254 Index
suspension polymerization 120 T
weight average molecular weight
templates
101
for simulations 6
styrene distillation
tertiary dodecyl mercaptane 145
example flowsheet 106
thermal conductivity
example input summary 109
predicting 29
example kinetics 109
time limits 9
example process conditions 107
titles for simulations 8
example results 114
toluene 106
inhibitor 109
polymer production results 115
profitability results 116 U
sensitivity study 115 UNIPOL flowsheet 233
styrene ethyl acrylate unit operation blocks specification
copolymerization example 130 16
44 unit-sets
example flowsheet 131 defining 8
example input summary 136
example kinetics 132
V
example process conditions 131
example results 140 Van Krevelen group contribution 29
molecular weight 142 vapor-liquid equilibrium 40
monomer conversion 141 vinyl acetate 52
parameter regression 133 viscosity
styrene-butadiene intrinsic 84
chain length distribution 155 zero-shear 84
copolymer composition 163 VK column 203
degree of polymerization 164
example flowsheet 146, 158 Z
example input summary 150,
160 zero-shear viscosity 83
example kinetics 148, 159 Ziegler-Natta catalysts 165
example process conditions 147, Ziegler-Natta polymerization
159 of ethylene 16583
example process studies 150 of polypropylene 23250
example rate constants 149
example results 153, 162
monomer mass 154
particle growth 153
particle nucleation 153
phase volumes 154
polymer composition 155
polymer mass 154
polymer properties 154
rate of consumption 162
styrene-butadiene rubber See
styrene-butadiene
subroutines
Fortran 83
user 90
suspension polymerization 12029,
23250
Index 255