Tr. .mv.cvuar
u.-
rK-nw
m- Ilillillilllllll’l‘lllllmill»
3 1293 01025 0920
This is to certify that the
thesis entitled
SORPTION OF ORGANIC PENETRANTS BY AMORPHOUS POLYAMIDE
presented by
MYUNGHOON LEE
has been accepted towards fulfillment
of the requirements for
MASTER degree inPACKAGING
Major professor
DateélliIL 1. 1994
0-7639 MS U is an Affirmative Action/Equal Opportunity Institution
LIBRARY
Mlchlgan State
Unlverslty
PLACE IN RETURN BOX to romovo this checkout from your record.
TO AVOID FIN IS rotum on or before date duo.
DATE ‘DU E DATE DU E DATE DU E
r ' 45?”:
FEB 0 “
21.059
. “.991
FE! 2
MAGIC 2
W911i.
W“
L‘flll
MSU!.A." .u. .L A - '1 ”.rr
SORPTION OF ORGANIC PENETRANTS BY AMORPHOUS
POLYAMIDE
By
Myunghoon Lee
A THESIS
Submitted to
Michigan State University
in partial fulfillment of the requirements
for the degree of
MASTER OF SCIENCE
School of Packaging
1994
ABSTRACT
SORPTION OF ORGANIC PENETRANTS BY AMORPHOUS
POLYAMIDE
By
Myunghoon Lee
Sorption studies involving the sorption of n-propanol by an amorphous
nylon ( Nylon 61 /6T ) were carried out as a function of sorbate vapor
activity at 23 0 C.
Vapor activity levels from 0.035 to 0.91 were investigated to evaluate
the concentration dependency of sorption mechanism.
Sorption behavior of propanol by Nylon 6I/6T showed distinctive two
mode sorption phenomena as a function of vapor activity . At vapor
activity levels below a=0.11, equilibrium sorption was achievediwithin a
short period of time ( less than 20 hours ) , which can be interpreted as
fellowing a Fickian diffusion model. i
A Langmuir- Flory-Huggins dual mode sorption model can also be
applied at these concentration levels. However , for vapor activities above
a=0.11, the sorption process appeared to be non-Fickian and resulted in a
lack of equilibrium being attained.
ACKNOWLEDGEMENTS
I'd like to express my heartfelt thanks to my major advisor Dr. Jack . R .
Giacin. His devoted and endless supports enabled me to complete this
research. I would also like to thank Dr. Ruben Hernandez and Dr. Eric. A.
Grulke for their precious advice and discussion throughout the course of
research. I also thank to Dr. Heidi and Mr. Bob Hurwitz for helping me to
design and assemble the test system and to provide comfortable
circumstance during the research.
I am deeply indebted to the Korea Institute of Industrial Design and
Packaging , which sent me to study in the School of Packaging. Finally ,
my lovely wife and two daughters deserve to have my deep thanks and love
for their devoted supports during the research.
ii
TABLE OF CONTENTS
LIST OF TABLES
LIST OF FIGURES
INTRODUCTION
LITERATURE REVIEW
Characteristics of Amorphous Polymers
Properties of Amorphous Polyamide
Properties of Penetrant
Sorption Mechanisms
Ideal sorption and diffusion
Non ideal sorption
Modified dual mode sorption
Factors affecting the sorption process
Temperature
Crystallinity
Orientation
Sorption measurement
iii
Page
vi
13
14
14
16
20
22
22
25
27
30
Diffusion and relaxation in glassy polymers
Correlation of sorption with penetrant size
MATERIALS AND METHODS
Polymer film
n-Propanol
Sorption measurement
Density experiment
RESULTS AND DISCUSSION
Data preparation and prelimenary test
Equilibrium sorption
Analysis of sorption at vapor activity a < 0.11
Analysis of sorption at vapor activity a > 0.11
CONCLUSION
APPENDIX
Appendix A Procedure for standard calibration
curve constraction
Appendix B Electrobalance calibration procedure
Appendix C Vapor activity fluctuation
BIBLIOGRAPHY
iv
31
33
35
35
35
36
39
47
49
53
67
68
70
71
76
LIST OF TABLES
Table Title Page
1. Comparison of the properties between Amorphous 12
Nylons and Nylon 66
2. Relationship between rotameter readings and 45
vapor activities
3. Volume fraction vs. vapor activity from Flory- 51
Huggins equation
4. Average vapor activity and standard deviation 52
in each step of vapor activity
Figure
l.
2.
10.
ll.
12.
13.
14.
15.
16.
17.
LIST OF FIGURES
Title
Structure of Nylon 61 / 6T
Structural formulas of various alcohols
Schematic diagram of electrobalance test apparatus
Saturated vapor pressure of n-propanol
Standard calibration graph of n-propanol at 23 0 C
Illustration of sample flushing by nitrogen at 23 0 C
Evaluation of the reproducability of the sorption system
at a = 0.06
Interval sorption of n-propanol at a = O - 0.035
Interval sorption of n-propanol at a = 0 .035- 0.050
Interval sorption of n-propanol at a = 0.050 - 0. 80
Interval sorption of n-propanol at a = 0.080 - 0.110
Interval sorption of n—propanol at a = 0.110 - 0.295
Interval sorption of n-propanol at a = 0.295 - 0.390
Interval sorption of n-propanol at a = 0.390 - 0.490
Interval sorption of n-propanol at a = 0.490 - 0.840
Interval sorption of n-propanol at a = 0.840 - 0.910
Successive interval sorption of n-propanol into Nylon 61/6T
at 23 OC(a=0-0.9l)
vi
Page
10
13
41
42
43
56
56
57
57
58
58
59
59
60
6O
61
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
Successive interval sorption of n- propanol into Nylon 61/6T
at 23 OC(a=0-0.110)
Experimental sorption isotherm for n-propanol in Nylon
6I/6T at the low vapor activities ( a< 0.11 )
Flory - Huggins plot when x = 1.6 from the equation
Ina1=1nV1+ v2 + xv22,v,=1- V3
Plot of Mt/Moo vs. t 1/2 for n-propanol sorption by
Nylon 6I/6T, at 23 0C ( a = 0 - 0.035)
Plot of Mt/Moo vs. t 1/2 for n-propanol sorption by
Nylon 6l/6T, at 23 0C ( a = 0.035 - 0.050)
Plot of Mt/Moo vs. t 1/2 for n-propanol sorption by
Nylon 61/6T, at 23 OC ( a = 0.050 - 0.080)
Plot of Mt/Moo vs. t 1/2 for n-propanol sorption by
Nylon 6I/6T, at 23 0C ( a = 0.080 - 0.110)
Comparison between experimental and calculated
data at vapor activity a =0.035
Deviation from Equation(21) due to non - Fickian
sorption for sorption of n-propanol by Nylon 6I/6T
at 23 0C( a = 0.11 - 0.295): Graphic estimation of
Fickian Moo value
Plot of Mt/Mco vs. t 1/2 for n-propanol sorption by
Nylon 6I/6T, at 23 0C; a = 0.11 - 0.295, Equation
(21) for Moo = 0.275 and t 1/2 = 84 sec.
Plot of Mt/Moo vs. t 1/2 for n-propanol sorption by
Nylon 6I/6T, at 23 0C; a = 0.11 - 0.295 ; points
observed curve,Equation(21) for Moo = 0.275
and t 1/2 = 84 sec.
vii
61
62
62
63
63
65
65
66
66
29.
30.
32.
33.
35.
36.
37.
Distribution of vapor activity at
Distribution of vapor activity at
Distribution of vapor activity at
Distribution of vapor activity at
Distribution of vapor activity at
Distribution of vapor activity at
Distribution of vapor activity at
Distribution of vapor activity at
Distribution of vapor activity at
a = 0.035
a = 0.050
a = 0.080
a = 0.1 10
a = 0.295
a = 0.390
a = 0.490
a = 0.840
a = 0.910
viii
71
71
72
72
73
74
74
75
INTRODUCTION
Over the past decade , the development and use of high barrier
plastic materials has increasingly replaced traditional glass and metal
containers for food and beverage packaging.
The advantages of plastic packaging are numerous and include :
low cost . light weight , a wide range of mechanical properties ,transparency.
flexibility , direct food contact , and general consumeripreference
because of convenience and microwaveability ( Salame , 1986).
However , with the use of plastic packaging there are concomitant
concerns related to product / package interactions. This is a broad based
topic that includes transport of gases and organic vapor of low molecular
weight : (i) from the product through the package , as well as (ii) from
the environment through the package to the product.
The specific mass transport ( permeability ) process may be described as
a function of penetrant concentration , temperature , and time to equilibrium,
as well as the composition of the penetrant / polymer system.
Therefore , polymer structure , free volume , chain stiffness or
segmental mobility , and the availability of specific sites of interaction in the
polymer,plus the physicochemical characteristics of the penetrant (gas, water
vapor , or organic vapor) determine the mode and mechanism of sorption
1
2
and thus . the resultant transport and mechanical properties of the polymer
( Blatz , 1989 ).
The dual mode sorption model has been widely used to describe the
solubility of gases in glassy polymers . and in glassy , polar polymers
as well. This model assumes that the solute molecules in the glassy
polymer consist of non-specifically absorbed and specifically adsorbed
species , which are in dynamic equilibrium in the medium ( Michaels et al.,
1963 )
The solubility of the absorbed species is represented by Henry's law ,
and the solubility of the adsorbed species is described by a Langmuir type
adsorption isotherm. The Langmuir sorption is believed to occur at
specific sites , usually considered to be units of stable free volume present in
the amorphous polymer structure.
Local equilibrium between the absorbed and adsorbed populations is
maintained throughout the polymer matrix. The total amount of solute
sorbed by both mechanisms is :
CH' b
C=CD+CH=KDP+——— (1)
l + bp
where C is the total concentration of sorbed solute in the polymer ,
CD and CH are the solubilities due to absorption ( Henry's Law ) and
Langmuir type adsorption respectively , K1) is the Henry's law
3
constant . b is the hole affinity constant . CH' is the hole saturation
constant and p is the pressure.
Recently , Hernandez et a1. ( 1991 ) have studied the sorption
of water vapor into an amorphous polyamide ( Nylon 61 / 6T ).
The authors have applied a dual mode sorption model , based on
Langmuir and Flory - Huggins equations , to describe the sorption
process. The total amount of water sorbed by the polymer was described
by the summation of a Langmuir type association , and by a solution
component mode given by the Flory-Huggins model , which is given by:
V1 = V1L + V1FH (2)
Where V1 is the total volume fraction of water within the polymer , and the
superscript L and FH refer to Langmuir and Flory - Huggins water sorption
contributions , respectively. These contributions are expressed as :
K a
V1L= (3)
1 + Ba
and
a=exp[1nV1FH+(l~V1FH)+x(1-V1FH)2] (4)
Where a is water activity , x is the Flory - Huggins interaction parameter ,
and K and B are parameters of the Langmuir equation ( Hernandez et a1.,
1991).
The behavior of oxygen solubility within the polymer / water system
4
was related to V114 , according to the following expression.
V0:2 = V* - FV]L (5)
where V02 is the solubility of oxygen in the polymer as a function of a ,
V* is the solubility of oxygen in the dry polymer , and F is a factor that
relates sorption values of oxygen at dry conditions and the fraction volume
of water described by the Langmuir sorption mode , when water activity
equals one. The fact that 80 % of the total oxygen dissolved by the
polymer at dry conditions was displaced by molecules of water associated
with active sites of the polymer matrix , indicated the importance of these
active sites in the mechanism of the solubility of oxygen within the
Nylon 61 / 6T bulk phase. These results also suggested that molecular size
may be an important factor in determining the final equilibrium sorption
values in the three component system of polymer , water and oxygen
( Hernandez et al. , 1991 ).
Based on the results described above , it can be proposed that a similar
model can be applied to the binary system of polymer and organic penetrants
,which has not yet been reported.
A quantitative relationship between molecular weight of sorbates and
active sites concentration in polyethylene was reported by Gedraitite et a1.
( 1989) , who studied the sorption by polyethylene , of selected low
molecular weight compounds dissolved in n - hexane , and in isopropanol.
D
In addition to the limited amount of data describing the effect of
moisture content on the mass transport and solubility of oxygen in
hydrophilic polymers , no references have been found in the literature
which considered the importance of size distribution of sorption sites in
determining the final equilibrium sorption values for the sorption of organic
penetrants by hydrophilic polymers , such as the amorphous polyamide.
There is also a paucity of data describing the effect of organic vapor
concentration on the solubility and diffusion of organic penetrants in
hydrophilic polymers , and the importance of active binding sites associated
with the polymer matrix in determining the final equilibrium sorption value
for the binary system of polymer / organic vapor.
The proposed studies will provide for a better understanding of the
parameters contributing to the diffusion and sorption of organic penetrants
in hydrophilic polymer membranes .
In terms of practical application , the data obtained describing the
relationship between sorption parameters and the molecular weight of
the sorbate molecule may allow calculation of the solubility of a
penetrant in the polymer matrix from knowledge of the penetrants solubility
in a contact phase , and its molecular weight. This could further provide a
means of designing a barrier structure for a specific end use application.
For example , in the packaging of a juice product , when there is
concern for flavor loss due to sorption _. or sorption of ingredients from a
6
pharmaceutical preparation . which could result in loss of product efficacy.
Objectives
1. Determine the equilibrium sorption isotherm for n-propanol in an
amorphous polyamide ( Nylon 61 / 6T ) at water activity 2 0.
2. Provide an empirical framework to study the sorption of organic vapor by
a glassy polymer , and the relationship of Langmuir type association
( active site binding) , to the solubility of organic penetrants in
hydrophilic polymers.
3. Evaluate the appropriateness of the Langmuir-Flory-Huggins dual mode
sorption model to describe the sorption isotherm of n-propanol in the
amorphous polyamide( Nylon 61 / 6T ).
LITERATURE REVIEW
Characteristics of Amorphous Polymers
A vast majority of well-known polymers are characterized by their
structural differences. The conformation of flexible chain molecules in the
melt or glassy state is determined both by the local position and
orientational distribution of chain segments belonging to the same chain
( intramolecular correlations ) or belonging to different chains
(intermolecular correlations ) (Wendorff , 1987 ).
Polymers differ from other low molecular weight organic or inorganic
materials , in that they are composed of long chain molecules , usually
containing a huge number of atoms. This gives rise to a very large number
of internal degrees of freedom for flexible chain molecules ( Blumstein ,
1978 ). The chain is represented by a line in space , the curvature of
which depends upon the persistence length A1. The parameter(Ai) is thus a
measure of the chain stiffness. It corresponds to the distance along the
chain over which correlations in the orientation of successive chain units
extend. The mean square end-to-end distance < h2 > is given by
/2AiL=[1-(l/x)][1-exp(-x)] (6)
where L is the contour length of the chain and x = L /Ai. A worm- like
chain becomes a flexible Gaussian chain for x —> 00 , and a rigid rod- like
7
8
chain for x —> 0 ( Wendorff . I987 )-
Flexible chain molecules have a persistence length of the order of 1 or
2 nm. A flexible chain in a random configuration occupies only a small
portion of the space it pervades. The density of an isolated chain molecule
is typically of the order of l % of that of the condensed state , and it
decreases with increasing chain length ( Wendorff , 1987 ).
Intermolecular forces require , however , that the empty space be filled in
the condensed , amorphous state. This may be achieved in various ways as
described below. i
(l) The chains may collapse , a process which is known to occur in solutions
with poor solvents ( Williams et al. , 1981) ; (2) The chains may aggregate
in bundles in which neighboring chains or chain segments are parallel to
each other as in a nematic liquid crystalline phase ( Pechhold et al. , 1979 ) ;
(3) The chains may remain in their random configuration and many different
chains are able to pervade the space taken up by the reference chain.
This leads to a highly entangled system ( Gennes , 1979 ).
The distribution of the centers of the chains as well as their orientation is
as close to random , as in a real gas. Macroscopic properties of amorphous
polymers are known to depend heavily on the local spatial and orientational
distribution of chain segments , belonging to the same or to a different
chain (Wendorff, 1987).
9
Properties of Amorphous Polyamide
Polyamides containing aromatic residues are normally very high melting,
highly insoluble semi-crystalline materials , which are difficult to process.
In recent years , a number of more easily processed amorphous
polyamides , based on iso— and terephthalic acids have been claimed in the
patent literature. In general , amorphous polyamides can possess a wide
range of good mechanical properties , including tensile and impact
strength , and a high specific modulus , which are desirable for engineering
applications ( Dolden , 1976 ).
The amorphous polyamide ( Nylon 61 / 6T ) which was developed by
the E.I.Du pont De Nemours & Co., Ltd., has a random structure based on
hexamethylenediamine and a 70/ 30 mixture of isophthalic and
terephthalic acids ( Fig. l ), and when dry has a glass transition temperature
(Tg ) of 125 0 C ( Blatz ,1989).
It has good processability by primary extrusion and molding
operations , as well as secondary processing such as thermoforrning.
Film has been produced by both the blown film and flat cast processes.
Coextruded film and sheet can be produced using three to five layers , at
thicknesses from 2 mils to 60 mils.
IO
NHz CH2 - (CH2)4 - CH2 NHz -> Hexamethylene diamine
+
COzH C02”
| + |
\
COzH
|
COzH
lsophthalic acid Terephthalic acid
(70 % ) ( 30 % )
Figure 1. Structure of the Nylon 61 / 6T
Nylon 61 / 6T also has good tensile strength , elongation and tensile
modulus. These properties are similar to those for the other resins used
for quality packaging applications. The tensile strength of amorphous
polyamide is higher than nylon 66 dry , but decreases in a manner
similar to nylon 66 at 50 % RH. The tensile modulus , a measure of
stiffness , shows an atypical trend with increasing humidity. As
expected , the modulus of nylon 66 decreases substantially as the
humidity increases , however the modulus of amorphous polyamide
11
increases slightly. Thus . it maintains its stiffness at high
humidities ( Blatz , 1989 ).
The barrier properties of amorphous polyamide approaches the
levels of polyvinylidene chloride , ethylene vinyl alcohol copolymer
and the extrudable acrylonitrile copolymers. As was previously indicated,
water vapor has an unexpected effect on the permeability of
amorphous polyamide. Whereas the oxygen permeability of nylon 66
film increases on exposure to water vapor , that for the amorphous
polyamide film decreases. Thus , a film of amorphous polyamide has
a significantly better barrier as the humidity increases ( Blatz , 1989 ).
The chemical resistance is usually not as good as the semi-crystalline
resins , such as PET or 6 and 66 Nylon. The resin is unaffected by
dilute acids and bases , but is attacked by acetic acid ( Blatz, 1989 ).
It is also attacked by the common low molecular weight alcohols such as
ethyl alcohol and solutions containing more than 10 % ethyl alcohol.
The amorphous nylon is resistant to both aliphatic and aromatic
hydrocarbons , to chlorinated solvents , ketones , and esters. However ,
methylene chloride will cause swelling ( Blatz, 1989 ).
As with all glassy polymers , as the melt is cooled down and
reaches the glass transition temperature ( Tg ) , the molecular motions
are frozen at that temperature , which results in a void content that is higher
12
than that in a polymer with low Tg . For an amorphous glassy polymer .
the larger the difference in temperature between the Tg and the test
temperature the greater the amount of unrelaxed volume ( sometimes called
the free volume ). Table 1 summarizes selected mechanical and barrier
properties of amorphous nylon ( Nylon 61 / 6T ) and nylon 66.
Items Unit Amorphous Nylon* Nylon 66
Tensile strength K Psi 10.0 8.0
Elongation % 50.0 70.0
Tensile modulus K Psi 330 340
Heat deflection T. @ 66 psi 26] 0 F (128 0 C) 365 0 F (185 0 C)
of packaging resin @ 264 psi 248 0 F (120 0 C) 151 0 F (66 0 C)
Oxygen Dry 3.3 2.5
Permeability (1) 50 % RH 1.4 5.0
CO2 8.0 10.0( 50 % RH)
Permeability (2)
Water vapor 5.0 5.7
Permeability (3)
* SELAR® PA of Du Pont Co.,
(1) CC-mil / 100 SQ. IN. / 24 Hours-Atm. at 30 0 C
(2) CC-mil / 100 SQ. IN. / 24 Hours- Atm. at 80 % RH, 30 0 C
(3) g-mil / 100 SQ. IN. / 24 Hours- Atm. at 90 % RH, 23 0 C
Table 1. Comparison of mechanical and barrier properties of amorphous
nylon* and nylon 66
13
Properties of penetrant : l - Propanol
Alcohols may be viewed as alkyl derivatives of water and are indicated
by the formula R - OH. They are structurally similar to water , but have
one of the hydrogens replaced by an alkyl group.
The reactions of alcohols are those of the -OH functional group , known
as the hydroxyl group .
Figure 2 lists both the IUPAC and common names for a variety of
simple alcohols , which include the 1 - propanol used in this study.
CH3OH CH3CH20H CH3CH2CH20H
methanol ethanol l-propanol
methyl alcohol ethyl alcohol n-propyl alcohol
CH3CHCH3 CH3CHCH20H CH3
| |
OH CH3 CH3COH
l
CH3
2-propanol 2-methyl-l-propanol 2-methyl-2-propanol
isopropyl alcohol isobutyl alcohol tert-butyl alcohol
Figure 2. Structural formulas of alcohols.
14
Sorption Mechanisms
Organic vapors can exhibit concentration dependent mass transport
and sorption processes. The permeant vapor pressure and the type of
vapors that come in contact with the package will also determine the
magnitude of sorption and permeation into and out of polymeric packaging
systems.
The study of the solubility of gases in polymers has a twofold
objective : (1) it permits the establishment of correlations of gas
solubilities , with readily available gas molecular parameters ; and (2) it also
supplies information on the morphology of polymers ( Vieth ,1966 ).
Ideal sorption and diffusion
In ideal gas - polymer systems , both the solubility and diffusion
coefficients are constant at any given temperature , and so the permeability
is also a constant.
Permeation or transport through a polymer film can be described in
terms of its component parts by Equation (7):
P = D S (7)
where S is the solubility coefficient , and D is the diffusion coefficient.
S characterizes the amount of permeant that can be dissolved into the
15
polymer under the given vapor pressure . and D describes the rate at
which the permeant molecules are advancing through the barrier film.
For a simple permeation process , the sorption and desorption steps are
described by the assumption of Henry's law , which relates the
concentration of the penetrant in the polymer , to the vapor pressure in
equilibrium with the polymer. The partial pressure of the penetrant is
further related to the penetrant concentration in the gas phase through the
ideal gas law.
Generally , the application of the ideal gas law is justified since the
concentration of the diffusant in the gas phase is very low. The
diffusion step is described by Fick's first (8) and second laws (9) of
diffusion ( Hopfenberg , 1973 ).
dc
F = - D(c)
(8)
6x
dc 6] 6c\
= —| D(C)- | (9)
dt 6x \ 6x)
where F is the flux or the rate of transfer of penetrant per unit area ,
expressed as mass of diffusant per unit area per time ;c is the
concentration of the penetrant in the film , expressed in the same unit of
mass of diffusant per unit of volume or mass of the polymer.
D is the mutual diffusion coefficient , in ( length )3/ time ; t is time ;
and x is the length in the direction in which transport of the penetrant
16
molecules occurs.
To obtain the flux (F) or the diffusion coefficient (D) from
Equations (8) or (9) , initial and boundary conditions associated with the
experimental method are needed , and the expressions solved to give the
desired values.
It should be recognized that when the diffusion coefficient is calculated
using Equations (8) and (9) , only approximate values will be obtained.
More accurate estimations of this parameter can be carried out by using ,
for example , a nonlinear maximum likelihood sequential method , based
on the Gauss linearization method ( Beck and Arnold, 1977).
The partial pressure of the penetrant is further related to the penetrant
concentration in the gas phase through the ideal gas law. Application of
the ideal gas law is justified since the concentration of the diffusant in the
gas phase is , in general , very low ( Hernandez et al. , 1986 ).
Non-ideal sorption
The sorption of a low molecular weight compound proceeds by two steps.
Firstly the sorbate penetrates into the relatively ordered polymeric
substance , forming a true solution. The concentration of the
homogeneously dissolved compound is related to the concentrations outside
the polymer by Henry's Law. Secondly the dissolved compound is
1 7
reversibly sorbed by certain centers . which are the zones of destruction of
the short range order. The total sorbate concentration in the polymer will
be the sum of the concentrations of the homogeneously dissolved compound.
and of that sorbed by the sorption centers ( Shlyapnikov and Mar'in , 1987 ).
For systems in which the solubility doesn't conform to Henrys law , both
the diffusion and the solubility parameters are concentration dependent.
The solubility coefficient often is essentially constant at low vapor activities
for the more volatile vapors, and only the diffusion coefficient exhibits
significant concentration dependence ( Rogers et a1. , 1960 ). This
concentration dependence can be understood in the context of the free
volume theory ( Cohen and Tumbull , 1959 ; Fujita , 1961; Kwei and Wang,
1972 ; Peterlin , 1975 ) , which gives a relationship between the diffusion
coefficient D and the fractional free volume f :
D=Aexp(-B/f) (10)
where A is a frequency factor and B is a measure of the minimum hole size
for the jump process ( Choy et al. , 1984 ). Fujita ( 1961 ) further assumed
that
f = fo + B Cs (11)
where Cs is the mass fraction of the penetrant and 8 denotes the
effectiveness of the penetrant molecule for increasing the free volume of the
polymer. Substituting Equation ( l 1) into Equation (10) and assuming
BCs << fo , Equation (10) becomes
18
D=Doexp(st) (12)
where
y = 86/ to2 (13)
Peterlin ( 1975 ) has obtained an alternative expression for y in terms of the
fractional free volume of the penetrant , but the physical meaning is
essentially the same.
Equation (13) shows that the concentration coefficient 7 depends on fo
and f3 ; i.e.,y is a function not only of the fractional free volume but also of
the interaction between the penetrant and the polymer. For the two
penetrant-polymer systems , methylene chloride-HDPE and methylene
chloride-LDPE , B does not seem to vary significantly with drawing , so the
change in y is determined largely by the I / fo2 factor , which gives rise to a
substantial increase in y with increasing draw ratio ( Choy et al. , 1984 ).
The idea that the molecules of the penetrant present in a polymer are
of two different types (those mobile , homogeneously dissolved , and those
immobile , sorbed in certain centers ) , was called the dual mode sorption
model. This model has been widely applied to describe nonideal sorption ,
which can explain the irregularities of diffusion and sorption of low
molecular weight compounds ( i.e. , additives ) in glassy polymers ( Paul
and Koros , 1976 ).
The magnitude of the negative enthalpies of sorption reported by
19
Meats ( 1954 ) for neon . nitrogen . oxygen and argon in glassy polyvinyl
acetate and by Barrer et al. ( 1957 ) for organic vapors in ethyl cellulose
were inconsistent with the sorption theories of rubbery systems . and led
Barter et al. to suggest a two-mode , concurrent sorption mechanism for
glassy polymers , namely ordinary dissolution and 'hole' filling. In glassy
polymers , substantial deviations from linearity are often observed in
sorption isotherms at penetrant pressures above 1 atm.
Barrer et al. ( 1957 ) first suggested that two distinct mechanisms may
be operating in the sorption of low molecular mass species by glassy
polymers. They observed isotherms for the sorption of C4 and C5
isomers in glassy ethyl cellulose that curved towards the pressure axis , and
appeared to become asymptotically linear at high penetrant pressures.
These results were interpreted in terms of a "dual mode" sorption model ,
which assumes a combination of Langmuir - type trapping within preexisting
vacancies , plus "true" Henry's law solution. Quantitatively , this may be
written :
C'H P
C = 1(1) p + -——— (14)
l + b p
where C is the total sorbate concentration , K1) is the effective Henry's
law solubility constant , Cu is the Langmuir "capacity factor", b is the
Langmuir site affinity parameter , and p is the gas pressure.
20
Michaels et al. ( 1963 ) also reported the sorption of helium . nitrogen .
oxygen , argon and methane in glassy amorphous , and glassy
semi-crystalline polyethylene terephthalate. For penetrant partial pressures
up to 1.0 Henry's law was followed. However , carbon dioxide at 25 0 C
and 40 0 C , and ethane at 25 0 C, in the same pressure range , deviated from
Henry's law.
The carbon dioxide isotherms prompted Michaels et al. ( 1963 ) to
propose a two - mode sorption model of ordinary dissolution and adsorption
in microvoids ( 'holes' ) , for gas sorption in glassy amorphous polymers.
Michaels and Bixler( I961 ) attempted to adapt the Flory — Huggins
theory of polymer solution to the sorption of simple gases, picturing the
process as first condensation of the gas , followed by mixing of the polymer
and liquid penetrant.
Modified dual mode sorption model
Recently , Hernandez et al. ( 1991 ) also proposed a modified dual
mode sorption model , which described the sorption of water vapor by an
amorphous polyamide at 23 0 C. They modified Eqn. (14) by using the
Flory-Huggins equation to describe non-specific solution , rather than
Henry's Law. This modification allowed the model to fit over the
21
activity range . 0 < a < l . Although other choices for the solution
model are possible , they applied the Flory - Huggins Equation (15)
because of its simplicity.
Ina =an1+ V2 + )(V22 (15)
where x is the interaction parameter. The non-specific sorption term
is nonlinear , the complete model can have an inflection point and should
predict clustering somewhere over the range of activity values
( Hernandez et al., 1990 ). The modified dual mode sorption model is
expressed in terms of volume fractions and solute activity:
V1 = VlFH + V1L (2)
where V 1FH refers to the Flory - Huggins contribution to the solute
volume fraction , and V1L is Langmuir volume fraction contribution.
Since Equation (2) is nonlinear , it is convenient to determine the value
for V 1FH by numerical methods , such as Newton - Raphson technique
and Equation (2) becomes:
K a
V1: FH(a,x)+ (16)
1+ Ba
Factors affecting sorption behavior
Temperature
The temperature dependence of the sorption parameters have been
studied for a number of penetrant / polymer systems. The K1) and b
terms are generally decreasing as a function of T , and appear to be well
correlated at a given temperature with the penetrant Lennard - Jones
energy parameter E/K values ( Pace and Datyner , 1980). The data for
C‘ H are more ambiguous.
Some results suggest that OH is independent , or varies only weakly
with temperature , and others suggest that OH is a strongly decreasing
function of T(Pace, 1980). However, it is reasonable to assume that the
immobile species is contained in microvoids formed by local segmental
motions in the rubbery state which have become ”frozen" in the glass.
This frozen void fraction must be , to a first approximation , at least
independent of T to explain the low thermal expansivity of the glass ,
which is comparable to that of the crystalline state. Hence C'H should
also be essentially independent of T ( Pace and Datyner , 1980 ).
The glass transition temperature ( T8) of any amorphous substance ,
whether polymeric or not , is defined as the point where the thermal
expansion coefficient undergoes a discontinuity. In polymers , there
22
23
expansion coefficient undergoes a discontinuity. In polymers . there
may be more than one discontinuity associated with the thermal expansion
coefficient term.
The largest discontinuity is usually associated with the loss of the
molecular mobility . which permits configurational rearrangements of the
chain backbones ( Ferry , 1970 ).
Below the Tg of a polymer , not enough energy is provided to produce
the Micro-Brownian motion of chain segments 20-50 carbon atoms in
length , and the chains are locked into the position created from the
processing history. Sorption and diffusion in polymers at temperatures
below the polymer's Tg consist of more complex behaviors , and are not
well understood at this time ( Hernandez , 1984).
Above the Tg of a polymer , it seems likely that hole filling is an
important sorption mechanism , which is supported by the fact that the
low-pressure ( Henry's law ) solubility for simple gases varies smoothly
on going through Tg , although the apparent heats of solution are different
and constant in the glassy and rubbery regions.
Meares ( 1957 ) proposed the following idealised concept of the gaseous
sorption and diffusion process in polymers at temperatures above and below
Tg , based on the relative activation energies for diffusion and heats of
solution for polyvinyl acetate in the two temperature regimes , and the zone
24
of activation theory.
Above Tg , gas molecules which dissolve in the polymer must create
their own "holes" by separating the interchain polymer contacts.
The penetrant then diffuses through the polymer matrix along cylindrical
voids created by the synchronized rotation of polymer segments about the
C - C bonds ( Meares , 1957). Below Tg , the polymer consists of regions
of densely packed and arranged chains which have limited freedom for
rotation , separated by less dense regions of disordered chains that form the
'holes' into which the gas sorbs. The gas molecules diffuse between these
'holes' by the slight compressing of localized chains in the dense regions
enabling the gas molecules to pass through. This compression in the glassy
state does not create the long cavities common to the rubbery state,
indicating that the zone of chain activation is much larger in the rubbery
state. This zone of activation was found to be essentially independent of
the size of the penetrating molecule ( Meares , 1954 ).
Above the Tg the polymer is assumed to be in thermodynamic
equilibrium , while below Tg the hole concentration is that which was
frozen in at Tg , provided the glass is formed by slow cooling from the
rubbery region. One therefore expects
C'H=Coexp(—Hh/RT), T>Tg (17a)
C'H:Coexp(’Hh/RTg), T 0.11 ). The rate of sorption in
the initial stage is principally controlled by the Fickian process and
therefore, may be used to estimate a true diffusion coefficient, even when the
total sorption does not follow the Fickian model ( Berens , 1977) .
Furthermore , since the second sorption stage is almost entirely
relaxation - controlled , the kinetic and equilibrium parameters describing
the slow relaxation process can be obtained from long time sorption data.
To determine if the initial portion of the sorption profile , which
exhibited "two stage" sorption behavior , followed Fickian behavior , and to
estimate the diffusion coefficient the following treatment was applied with
the aid of a computer - assisted fitting routine.
Initially , values of M.» were intuitively selected by inspection of a plot
0f the experimental Mt vs. 1 data. This is illustrated for sorption studies
discribed in Figure 26, as shown, the Mac value for a Fickian diffusion
process is assigned , the t 1/ 2 value can be obtained from a plot of Mt / Moo
versus t“2 ( Figure 27 ) and the diffusion coefficient ( Ds ) calculated by
solution of equation (22).
By substituting the Ds value into Equation (21) , calculated values of Mt
are obtained.
The calculated and experimental values of Mt / Moo versus t“2 are then
plotted and the quality of the fit is evaluated by graphical analysis.
If the experimental and calculated values show poor agreement , this
55
procedure is repeated with a new value of Moo until a satisfactory f it
between the calculated and experimental values is obtained ( Figure 26 ).
As shown in Figure 28, an assigned value of Moo = 0.275 ,
gave good agreement ( Figure 28 ) between the experimental and calculated
data for the sorption profile plot at a = 0.295. The estimated diffusion
coefficient for the initial stage of the sorption profile ( i.e, Fickian ) was
1.264 13—10 cm2/sec.
Wt. Gain [m glm 9]
Weight Gainlm g]
56
0.07 _
0.06
0.05
002: —-o-- WLGng
“n.6mn1
I k I l l 1 4
. I
00 50 100 150 200 250
Time(hours)
Figure 7. Evaluation of the reproducability of the sorption
system at a = O - 0.06 ,
001 l
000
000
000
000
000 ‘ ' ‘ 4‘ ‘ '
0 10 20 30
fimaHm)
Figure 8. Interval sorption of n- propanol at a = O - 0.035
Wt. Gain[m 9]]
0.03
0.02 ’
0.01 '
0.00 . 1 1 1 1 1 . . 1
:13 35 37 39
Times(Hrs)
Figure 9. Interval sorption of n- propanol at a =0035 - 0.050
Weight Gain [m g]
0.10
0.05
0.00
0 2 4 6 8 10 12
Times(Hrs)
Figure 10. Interval sorption of n- propanol at a = 0.050- 0.080
Weight Gain [mg]
wt. gainImg]
0.0 A L A I 1 L l l
0.0 5.0 10.0 15.0 20.0
Times(Hrs)
Figurei 1. Interval sorption of n- propanol at a=0.080-0.110
80
Times(Hrs)
Figure12. Interval sorption of n- propanol at a=0.11-0.295
Wt. Gejn[m 9]
WI. Gain[mg]
59
Times(Hrs)
Figure13. Interval sorption of n- propanol at a=0.295-O.390
I l 1 l l L l 4 J
'0 so 100
Times(Hrs)
Figure14. Interval sorption of n- propanol at a=0.39-0.490
Wt. GainIm 9]
WI. Gain[m g]
60
0.0 1 l . l .
50 100 150
Times(Hrs)
Figure15. Interval sorption of n- propanol at a=0.49-0.840
l
0.0 A a 4 1 a a 1 a n A
50 1 00
Times(Hrs)
Figure16. Interval sorption of n- propanol at a=0.84-0.910
61
oooooooooooooo
1.5 - ,.-- 9
.8 .........
I .........
:I'
’3 1 O _ "7 ........... "
E ' f..-
g 6.
=3 5 1: a = 0.035 6: a = 0.390
E o s _ / """"" 2:a=0.050 7:a=0.490
' 3: a=0.080 8: a=0.840
3‘3"; 4:a=0.110 923:0.910
1 2 .---: 5: a = 0.295
0.0 ("""" 1 I 1 1 I 1 1 I L 1 I 1 1 I 1 1 I L 4
0 100 200 300 400 500 600 700
Time(hrs.)
Figure 17. Successive interval sorption of n- propanol into Nylon 61/6T
at 23 0C (a = 0 - 0.91)
P
0.20 -
g a=0.05 - 0.08
I.
‘E’
'8
Q
g 0.10 _ a=0.035 - 0.05
a = 0 - 0.035
0.00 1 1 I 1 1 I
0 50 100
Time(hrs.)
Figure 18. Successive interval sorption of n- propanol into Nylon 6I/6T
at2300(a=0-0.110)
Volume Fraction
Weight Fraction[glgx100)
0.30
0.25
0.20
0.05
0.00
62
1 I 1 J I 4 1 I 1 1 l #4 I J 4
0.02 0.04 0.06 0.08 0.10 0.12
Vapor Activity
Figure 19. Experimental sorption isotherm for propanol in
0.011
0.009
0.007
0.005
0.003
0.001
Nylon 61/6T at low vapor activities (a < 0.11)
I
I ' Y T T V I
4
j A L M J l A j l l L
I I 1
0.06 0.08 0.10 0.12 0.14
Vapor activity
I
0.04
Figure 20. Flory-Huggins plot when I =1.6,from the equation
Ina=InV1+V2+V2
63
8
2
‘2"
1 I 1 I 1 l 1 l 1
0 40 80 120 160 200
t1/2,Sec1/2
Figure 21. Plot of Mthoo vs t1/2 for n- propanol sorption
by Nylon 6|I6T, at 23 00 (a: 0.0 - 0.035)
110 _
100 [-
90 fi-
80 :-
70 :-
6 60 E
5 so 3
5 40 Z-
30 Z-
20 E-
10 r
0 g I 1 L 1 l 1
0 4o 80 120 160
11/2,sec1/2
Figure 22. Plot of Mt/Moo vs t1/2 for n- propanol sorption
by Nylon 6|/6T,at 23 00 (a = 0.035 - 0.050)
MIIMoo
Mtllvloo
110
100
64
I L I m I 1 I
40 80 120 160 200
t1l2,sec1l2
Figure 23. Plot of Mt/Moo vs t1l2 for n- propanol ’sorption
100'
90
80
70
60
50
40
30
20
10
by Nylon 6l/6T,at 23 0C (a = 0.05 - 0.08)
I 1 J
40 80 120 160 200
t1l2,sec1/2
Figure 24. Plot of Mthoo vs t1l2 for n- propanol sorption
by Nylon 6|/6T,at 23 00 (a = 0.08 - 0.11)
65
110
100' ‘1 J
90 ..
80':- 9
7O "
60
50
40
30 A (Mt/Moo)experi x100
28 r - (Mt/Moo)ca|cu. x100
1
O . I 1 I 1 I
0 100 200 300
11/25601/2
'U'
[MtIMoolxIUU
Y'U
Figure 25. Comparison between experimental and calculated
data at vapor activity a = 0.035
0.45
1
0.40 -
0.35
.
0.30
0.25
._. 0.20
CD
3; 0.15
-§ '1 :Moo 0.275 mg
0 0'10 2 :Moo 0.284mg
is; 0.05 3 :Moo 0.298 mg
g 0.0011I11111I14L11I1PIL41.J
0 10 20 30 4O 50 60 70 80
t, Hrs
Figure 26. Deviation from Eqn.(21) due to non-Fickian sorption
for sorption of propanol by Nylon 6l/6T at 230C(a=0.11-
0.295):Graphical Estimation of Fickian Moo value
66
Ivlt I Moo [gig]
OJ
0.0 1 I 1 I 1 I 1 I 1 I 1 I4 [*1 1 I 1
O 20 4O 60 100 120 140 160 180 200
11/2 = 84 sec. t1l2.sec1/2
Figure 27. Plot of Mt/Moo vs. t1l2 for n- propanol sorption by
Nylon 6I/6T ,23 00, a = 0.11 - 0.295,Eqn.(21) for
Moo = 0.275 and 11/2 = 84 sec.
MI I Moo [gig]
0.0 1 41 J
100 200
tI/2,sec1l2
Figure 28. Plot of Mt/Moo vs. t 1/2 for n- propanol sorption by Nylon
6II6T, 23 0C (a = 0.11 - 0.295 ): point observed,
curve, Eqn.(21) for Moo = 0.275 and t 1/2 = 84 sec.
Conclusion
The sorption of propanol by amorphous polyamide, evaluated by an
incremental sorption technique, showed that simple Fickian kinetics were
obeyed and equilibrium was achieved within a short time at very low vapor
activities ( a1 < 0.11 ) . However , at high activity levels ( a1 > 0.11 ) , the
initial,rapid sorption was followed by a slower further uptake of
propanol,which was qualitatively attributed to a relaxation-controlled
swelling process. The initial rapid stage of sorption may be related to the
Fickian diffusion of penetrant into pre-existing and available vacancies or
sites in the glassy polymer. The relaxation processes , which occur over a
long time scale may be related to a structural reordering or redistribution of
free volume elements to provide additional sites of suitable size and
accessability to accommodate more penetrant molecules.
67
Appendix A
Procedure of standard calibration curve construction
In all cases , a standard curve of response vs. penetrant concentration
was constructed from standard solutions of known concentration.
Calibration solutions were prepared by dissolution of known quantities of n-
propanol in dichlorobenzene. The detailed procedure follows:
1. Bake out vials and syringes in oven prior to use to remove any residual
solvent or permeant. Cool to room temperature.
2. Using the gas chromatograph , check the retention time of each solvent
by head space technique to ensure there is no interferring peaks at the
propanol retention times.
3. Prepare dilute solutions for the permeant standards by the following
procedure:
- fill up 100 ml volumetric flask with 10 ul of n-propanol and dilute to
make with dichlorobenzene. This becomes 100 ppm standard solution.
- make up similar standard solution with dichlorobenzene , which
becomes a standard solution of 50 ppm concentration.
- Make 2 or 3 more different concentration solutions by using the same
method
4. From the lowest concentration solution , inject 0.5 111 sample solution
directly into the gas chromatograph and the area response recorded.
68
69
Replicate the run and calculate the average
Plot the gas chromatograph area unit response vs. the number of grams
injected per sample. The slope of this curve equals the calibration
factor. Injected quantity can be calculated as follows
W .—: concentration(v/v) x 0.5 111 x 1 ml/ 1000 111 x propanol density(g/ml)
Setting conditions of gas chromatograph are as follows :
o Oven temp. : 50 0 C 0 Initial time : 5 min.
0 Rate : 3 0 / min. 0 Final time : 0 min.
0 Final temp. : 165 0 C o Inject temp.: 200 0 C
o Detect temp.: 250 0 C o Oven max. : 250 0 C
0 He carrier gas : 33 ml/min. 0 Sample injected : 0.5 ul
0 Head space test : 100 1.11 inject
0 Concentration used : 100 , 50 , 25 , 10 , 5 ppm ( 5 kinds )
Appendix B
Electrobalance Calibration Procedure
Cahn D-200 Electrobalance can be calibrated by the following procedure.
. Suspend the sample container on the hangdown wire
. Press the down arrow key three times and the NET WEIGHT box will
display the weight that the balance is actually seeing.
. Enter "Y" into the REZERO box. Next, press the up arrow key once to
go back to the Balance Calibrate menu.
. Enter into Tare Balance box "Y" and the balance will tare.
. Press the down arrow key to input the calibration weight
. Remove the hangdown tube so that the pan is exposed. Then put the
100 mg calibration weight onto the sample weight pan. Replace the
hangdown tube for stable reading.
. Place the cursur into the Calibrate Balance box and enter "Y" and press
ENTER when the pan with the calibration weight is steady.
. Lower the hangdown tube and remove the calibration weight.
. Suspend the sample and set up run conditions.
70
Yap or Activity
trap or Actin’ty
7 1
Appendix C.
Vapor activity fluctuation
0.10
0.08 -
0.06 -
0.04 "' Bun B D—
002 .. Mean: 0.035
Standard deviation: 0.003
000 I 1 I 1 l - 1
0 10 20 30 40
Time(Hrs)
Figure 29. Distribution of vapor activityat a = 0.035
0.10
0.08
0.06
0.04
0.02
0.00
.
.._B___n—-u—F T
Mean:0.050
" Standard deviation:0.0016
1 l 1 I 1 J 1 I L I 1 I 1
1 2 3 4 5 6 7 8
' Times(Hrs)
Figure 30. Distribution of vapor activity at a = 0.050
i.iap or Activity
trap or Activity
\1
N
0.12 _
0.10 :-
E a E,
0.08 _ 43—43— r
0.06 _-
0.04 :-
I Mean: 0.08
0,02 :. Standard devration: 0.0042
0.00 : J— 1 J l L l 1 L 4
6 8 10 12 14
Times(Hrs)
Figure 31. Distribution of vapor activity at a = 0.080
0.30
0.25 1
0.20 -
0.15 1
0.10 fir—”J; F +
O 05 E Mean: 0.110
' - Standard deviation: 0.003
O’COI- l 1 l 1 l 1 I 1
O 10 15 20 25 3O
Times(Hrs)
Figure 32. Distribution of vapor activity at a = 0.110
Yap or Activity
Yap or Activity
0.50
0.40 -
0.35
0.45 -
h
030E. a 3334 W 5'
Mean: 0.295
0‘25 0 Standard deviation:0.0069
0.2011I11I11I11J11l11I1LI11
0 10 20 30 40 50 60 70 80
Times(Hrs)
Figure 33. Distribution of vapor activity at a = 0.295
0.50
0.45 i-
0.40 - fl 4L
.35
0.35 i' Mean: 0.390
Standard deviation: 0.0084
0.30 LIAAIAIIJIIILIAALAIII;
0 10 20 30 40 so 60 70 80
Times(Hrs)
Figure 34. Distribution of vapor activity at a = 0.390
Yap or Activity
Yap or Activity
0.60
0.58
0.56
0.54
0.52
0.50 in a
0.48
0.46
0.44
0.42
0'40 4 1 1 1 1 L 1 1 1 1 1 I 1 1 1 1
0 50 100 150
Times(Hrs)
U I I I V
ljl
l
E)
' I
mean: 0.490
Standard deviation: 0.015
I
V ‘l I
J
Figure 35. Distribution of vapor activity at a = 0.490
1.00
0.95 '-
O.90 '- a
L In
0.85 . a a
. a a a
0.80 .
: Mean: 0.84
0-75 _” Standard deviation: 0.027
0.70 h 1 l 1 l 1 1 1 l 1 L 1 I 1 l l_ J 1 l 1
0 20 40 60 80 100 120 140 160 180 200
Times(Hrs)
Figure 36. Distribution of vapor activity at a = 0.840
Yap or Activity
1.00
0.95
0.85
0.80
0.90
Mean: 0.91
Standard deviation: 0.0046
l 1 1 l 1 1 l 1 1 l l
1 1 I 1 1 1 1
20 40 ' 60 80 100 i 20 140
Times(Hrs)
Figure 37. Distribution of vapor activity at a = 0.910
BIBLIOGRAPHY
Bagley , F and F. A. Long, 1955. J. Am. Chem. Soc., 77 , 2172.
Barrer, R., J .Barrie and J. Slater, 1957. J. Polymer Sci.,23, 315
Beck,J . V and Arnold. K. J , 1977. "Parameter estimation in engineering
science" , John Wiley and Son , New York , NY.
Berens, A. R. 1975. Angew. Macromol. Chem., 47, 97.
Blatz, P. S. 1989. "Innovations in polyamide technology" AICHE
Meeting, Houston, TX, April 2 - 6.
Blumstein, A., 1978. "Liquid Crystalline Order in Polymers", Academic
Press, New York.
Choy, C. L., W. P. Leung and T. L. Ma, 1984. J. Polym. Sci: Polym.
Physics Edition, Vol. 22, 707-719
Cohen, M. H. and D. Tumbull , 1959. J. Chem. Physics , 31 , 1164.
Dolden, J G., 1976. "Structure property relationships in amorphous
polyamides" Polym. Vol. 17 , pp 875
Gedraitite, G. B., A. P. Martin and V. A. Shlyapnikov, 1989. Eur. Polym.
J ., 25:39.
Gennes,P. G., 1979. "Scaling Concepts in Polymer Physics", Cornell
Univ. Press, Ithaca, NY.
Hartley, G. S. , 1946. Trans. Faraday Soc., 423, 6.
Hernandez, R. J ., J. R. Giacin and A. L.Baner, 1986. J. Plastic Film and
Sheeting, 2 (3) : 187.
Hernandez, R. J ., J. R. Giacin and E. A. Glulke, 1990. J. Membrane
Science
76
77
Hopfenberg,H.B. and V.Stannett, 1973. "The physics of glassy polymers",
R. N. Haward(ed), chapter 9, The diffusion and sorption of gases and vapors
in glassy polymers, Applied science publishers Ltd, London, pp 504—506.
Kwei, T. K. and T. T. Wang , 1972. Macromolecules , 5 , 128.
Meares,P. , 1954. J. Am. Chem. Soc., 76, 3415 - 22.
Meares,P., 1957. Trans. Faraday Soc., 53 , 101 - 6.
Meares,P., 1958. Trans. Faraday Soc., 54 , 40 - 6. l
Michaels. A. S.,W.Vieth and J.Barrie, 1963. J. Appl. Phys.,34(1),1,l3 :
Michaels. A. S. and H. J. Bixler, 1961. J. Polym. Sci., 50, 393
Ohashi , K., R. J. Hernandez , J . R. Giacin and E. A. Gruke , 1991.
" Modeling the sorption of water vapor by a semicrystalline polyamide ",
Thesis , school of packaging , Michigan state univ.
Pace, R. J. and A. Datyner, 1980. "Model of sorption of simple molecules
in polymer", J. Polymer Sci., Vol. 18, 1103-1124
Pechhold, W. and H. P. Grossmann, 1979. Disc. Faraday Soc., 68, 58.
Peterlin , A. 1975. J. macromol. sci. physics , B 11 , 57
Rogers, CE, 1965. " physics and Chemistry of the Organic solid state" ,
(Eds Fox , Labes , ans Weissbergcr) , Interscience , NewYork , Ch 6
Salame, M., 1986. J. Plastic Film and Sheeting, 2 (4) : 321
Vieth, W. R., H. H.Alcalay and A. J. Frabe, 1961. J. Appl. Polymer Sci. 8,
2125.
Wendorff, J. H., 1987. "Studies on the nature of order in the amorphous
polymers", International Symposium.
Williams, C., F. Brochard and H. L. Frisch, 1981. Ann. Rev. Phys. Chem.
32, 433-451.
"‘11111111111111“