SummaryIn a copolymerization process two (or more) monomer species M1 and M2 chemically react to form long-chained polymer molecules. During this process the ratio of the rates of reaction of the monomers M1 and M2 usually does not equal the ratio of the monomer concentrations in the solution. In other words, the composition of the copolymer chains produced in a very small time period usually differs from the monomer feed composition. For that reason the mole fractions of the momomers in the feed, X and (1–X), change as the copolymerization reaction proceeds. Consequently, the mole fractions of the momomers found in the polymer chains that are being produced, Y and (1–Y), change in course of time too. Considering a reactive site at the end of the chain of a growing bipolymer, it could be a radical, or ionic species, or a coordination complex. Moreover, the active site can be either a remainder of monomer M1, or a remainder of M2, depending on which one of both monomers has been the latest to couple with the growing chain. Using a rewritten version of the so-called "copolymer equation", which is applicable to many polymer systems, the program calculates the copolymer composition Y from a given X and reactivity ratios r1 and r2. To track the composition drift, it performs reiterative calculations in which the computed instantaneous copolymer composition Y is used to compute the next feed composition X. The composition drift is visualized graphically by plotting the copolymer composition versus the feed composition. As the program runs, the compositions passed through are being represented as a lengthening curve, starting from the given X and ending as soon as the compositions no longer change. To set an example: the program lets you track the composition drift coming up for discussion in problem 7 of the preliminary examination example Polymeren I. Here the initial value of the mole fraction in the feed X equals 0,50 and the reactivity ratios r1 and r2 are 0,0 and 1,0 respectively (see also the solution help at page 3 of the document). You can download Polymeren I from here. |
Program Listing |
Elucidation and Help |
PROGRAM:COPOLYM | Title |
GridOn | > 2nd, FORMAT > GridOn, ENTER |
AxesOn:LabelOn | |
FnOff | > VARS > Y-VARS > On/Off... > FnOff |
Plotsoff | > 2nd, STAT PLOT > PlotsOff |
0 → Xmin | "→" Symbolizes the STO> key > VARS > Window... > Xmin |
94/62 → Xmax | |
.2 → Xscl : 0 → Ymin | |
1→Ymax : .2→Yscl | |
1 → Xres | |
5/62→K | |
1.5 → R | Example values r1=1.5, r2=0.1, X=0.8 |
.1 → S | |
.8 → F | Initial mole fraction X (monomer M1) |
ClrDraw | > 2nd, DRAW > ClrDraw |
Lbl PS | > PRGM > CTL > Lbl (Label) PS: Polymer System |
ClrHome | > PRGM > I/O > ClrHome |
Disp "Enter R1, R2, X",R,S,F | |
Input R | Reactivity ratio r1 |
Input S | Reactivity ratio r2 |
Input F | Initial feed composition X |
max(.001,R) → R | > MATH > NUM > max( |
max(.001,S) → S | |
Line (0,1,1,1) | > 2nd, DRAW > Line( |
Line (1,1,1,0) | |
Line (0,0,1,1) | |
F → X | |
X+.0001 → U | |
(RX²+X(1–X)) / (RX²+2X(1–X)+ S(1–X)²) + .0001 → V |
|
Lbl CD | CD: Composition Drift |
Repeat abs(V–Y)<K | > PRGM > CTL > Repeat > MATH > NUM > abs( > 2nd, TEST > "<" |
U+.7(X–U)→X | |
(RX²+X(1–X))/(RX²+2X(1–X)+S(1–X)²) → Y | |
End | > PRGM > CTL > End |
Line (U,V,X,Y) | |
abs(Y–X) → J | |
If J<.0005 | > PRGM > CTL > If |
Goto CC | > PRGM > CTL > Goto |
(X–Y) / J → H | |
Y → V : X → U | |
X+Hmin(K,.9J) → X | > MATH > NUM > min( |
min(1,X) → X | |
max(0,X) → X | |
Goto CD | |
Lbl CC | CC: Composition remains Constant |
Circle(X,Y,1/62) | > 2nd, DRAW > Circle( |
Pause | > PRGM > CTL > Pause |
Goto PS |
While running the program, a composition drift curve is drawn. When the program execution is interrupted, you may press the ENTER key to continue and add another curve to the graph; or press the ON key to stop the program. |
Hein ten Horn
Top | Programs Texas Instruments | Home