Determination of Viscoelastic Property in Polyethylene Crystallization Using a Quartz Crystal Resonator

Suh, Yong Kweon; Kim, Byoung Chul; Kim, Young Han

doi:10.3390/s91209544

Open AccessArticle

Determination of Viscoelastic Property in Polyethylene Crystallization Using a Quartz Crystal Resonator

by

Yong Kweon Suh

¹,

Byoung Chul Kim

² and

Young Han Kim

^2,*

¹

Department of Mechanical Engineering, Dong-A University, 840 Hadan-dong, Pusan, 604-714, Korea

²

Department of Chemical Engineering, Dong-A University, 840 Hadan-dong, Pusan, 604-714, Korea

^*

Author to whom correspondence should be addressed.

Sensors 2009, 9(12), 9544-9558; https://doi.org/10.3390/s91209544

Submission received: 26 October 2009 / Revised: 16 November 2009 / Accepted: 26 November 2009 / Published: 30 November 2009

(This article belongs to the Section Chemical Sensors)

Download

Browse Figures

Versions Notes

Abstract

:

A new generalized relationship between the viscoelastic properties of an overlayer placed on the electrode interface of a quartz crystal resonator and its resonant characteristic is developed from the mechanics of the quartz movement. The relationship is used to estimate the viscoelastic properties from the experimentally measured resonant characteristic. It is utilized in the estimation of viscosity and elastic shear modulus of a polyethylene overlayer during its crystallization. The measurements are compared with the viscosity and elastic shear modulus of a polyethylene melt measured using a rheometer. It is found that the development of this study is useful in the determination of viscoelastic property of polymer materials by measuring the resonant frequency and conductance of the polymer overlayer placed on the resonator electrode.

Keywords:

quartz crystal resonator; viscoelasticity; quartz mechanics; polymer crystallization

Graphical Abstract

1. Introduction

The variation of viscoelastic properties of a polymer material during its crystallization is important information in the design of a variety of polymer processing equipments. For instance, the property change has a key role in the design of a mold for an injection molding machine, but the property measurement is difficult while the polymer crystallizes.

A quartz crystal resonator is composed of a thin quartz crystal sandwiched between two metal electrodes that establish an alternating electric field across the crystal, causing vibrational motion of the crystal. The motion is characterized with the resonant frequency and admittance of the resonator, and the characteristic is sensitive to the changes of mass and physical property of an overlayer on its electrode. Because the resonator is so sensitive, it can sense a variety of changes in micro-scale at the electrode interface. In polymerization, the rheological properties of a reactant and product mixture vary continuously as the polymerization proceeds. The resonator has been implemented in the monitoring of a UV photopolymerization by measuring its resonant resistance [1]. The nucleation and crystal formation in a cooling crystallization have also been investigated with the quartz crystal resonator [2,3].

The resonant frequency and admittance of the resonator can be interpreted to the changes of mass and viscoelasticity of an overlayer at its electrode interface. Reed et al. [4] presented a physical description of a viscoelastically loaded resonator for either very thick viscous fluids or very thin elastic overlayers. Buttry and Ward [5] gave an extensive review of the quart crystal resonator including basic principles of an AT-cut resonator and relation of the mechanical system composed of a spring, mass and a damper to the electrical circuit of resistance, inductance and capacitance. Nwankwo and Durning [6] provided an impedance analysis of the resonator for very thick viscoelastic fluids. Kim et al. [7] demonstrated the applicability of a specially designed resonator system to monitor the change of viscoelastic property of thin polymer films subjected to temperature variation. When a PMMA/PVAc blend was coated on the one side of the resonator, the resonant frequency decreased accompanying hysteresis with the increased temperature. Based on the resonant frequency variation and the resonant resistance change, they demonstrated that there are two transient phase transition temperatures between the temperatures of 20 °C and 80 °C. Kunze et al. [8] developed an extended sheet-contact model to describe the change of resonant frequency and the dissipation of very thin viscoelastic solids coated on the electrode. Efimov et al. [9] studied the sensitivity variation of the resonator due to the energy trapping. It was found that energy trapping was insignificant for the small amount of mass loading, but the energy trapping became dominant and an oscillation occurred only in the region of the loading with a large loading. Thermoresponsive viscoelastic property of hydrogel was monitored with the impedance variation of a quartz crystal resonator [10]. A continuum mechanics model was utilized in the analysis of continuous viscoelastic profiles of a liquid film [11]. The frequency shift of viscoelastic overlayer has been interpreted with the small-load approximation [12]. A new set of equations was derived from the complex frequency shift of polymer brushed, and was applied to analyze the dissipation data [13].

In this study a generalized relation between the resonant characteristics of a quartz crystal resonator and the rheological properties of an overlayer applied on the electrode surface are developed from the mechanics of the quartz movement. The elastic shear modulus and viscosity of a polyethylene overlayer are estimated from the relation and the experimentally obtained resonant frequency and conductance of the resonator. The results are compared with the bulk property of polyethylene melt measured with a rheometer.

2. Theoretical Analysis

Consider the thickness-shear motion of a thin circular-disk-shape quartz crystal with thickness h_Q having a pair of concentric electrodes with radius r_e on both sides as shown in Figure 1.

The viscoelastic overlayer attached on the top electrode is assumed to be of axisymmetric shape with radius r_L and thickness h_L. Then the equation of motion for the quartz can be written as:

c_{66} \frac{\partial^{2} u}{\partial y^{2}} + η_{Q} \frac{\partial^{3} u}{\partial t \partial y^{2}} + e_{26} \frac{\partial^{2} ϕ}{\partial y^{2}} = ρ_{Q} \frac{\partial^{2} u}{\partial t^{2}}

(1)

e_{26} \frac{\partial^{2} u}{\partial y^{2}} - ɛ_{22} \frac{\partial^{2} ϕ}{\partial y^{2}} = 0

(2)

where t is time and (r, y) denotes the radial and axial coordinates of the cylindrical coordinate system. Further, u(r, y, t) is the mechanical displacement of the quartz along the x-direction, ø(r, y, t) the electric potential, c₆₆ the elastic shear modulus of the quartz, e₂₆ the piezoelectric constant of the quartz, ε₂₂ the dielectric constant of the quartz, η_Q the viscosity of the quartz, and ρ_Q the volume density of the quartz material.

The equation for the viscoelastic overlayer is similar to Equation (1) but purely mechanical:

μ_{L} \frac{\partial^{2} ν}{\partial y^{2}} + η_{L} \frac{\partial^{3} ν}{\partial t \partial y^{2}} = ρ_{L} \frac{\partial^{2} ν}{\partial t^{2}}

(3)

where v(r, y, t) is the mechanical displacement of the overlayer along the x-direction, μ_L the elastic shear modulus, η_L the viscosity, and ρ_L the volume density.

We have six boundary conditions for the above set of equations. Firstly, no-slip condition must be satisfied at the interface between the quartz and the overlayer (through the intermediate thin electrode of course). Secondly, the shear stress should vanish on the top surface of the overlayer. Next, the potentials on the top and bottom electrode surfaces are specified from a specified AC field. The fifth and sixth boundary conditions can be given by applying the Newton's second law to the mass of the top and bottom electrodes, respectively. Here, the forces acting on the electrodes may include the shear forces from the quartz and/or the overlayer (only for the top electrode).

Solutions to Equations (1), (2) and (3) can be written as:

u (r, y, t) = p (r) \overset{⌢}{u} (y) e^{i ω t}

(4)

ϕ (r, y, t) = {(e_{26} / ɛ_{22}) p (r) \hat{u} (y) + [p (r) E - (2 {\overset{⌢}{ϕ}}_{0} / h_{Q})] y + F (r)} e^{i ω t}

(5)

v (r, y, t) = p (r) \overset{⌢}{v} (y) e^{i ω t}

(6)

where ω and φ̂₀ the angular frequency and the amplitude of the external AC electric potential, respectively. The radial dependence of the displacements is represented by p(r) = exp(-r²/r_e²). Further, we have:

\overset{⌢}{u} = A exp (i k_{Q} y) + B exp (- i k_{Q} y)

(7)

\overset{⌢}{v} = C exp (i k_{L} y) + D exp (- i k_{L} y)

(8)

where k_Q=ω√(ρ_Q /c̃₆₆) denotes the complex wave number and:

{\tilde{c}}_{66} = {\bar{c}}_{66} + i ω η_{Q}

(9a)

{\bar{c}}_{66} = c_{66} + e_{26}^{2} / ɛ_{22}

(9b)

The five unknown constants A, B, C, D, and E and one unknown function F(r) can be determined from the boundary conditions. After some algebra we arrive at the following formula for E:

E = \frac{2 {\overset{⌢}{ϕ}}_{0}}{h_{Q} P_{e}} \frac{- K^{2} [2 (1 - cos ψ_{Q}) + 2 q_{e} ψ_{Q} + q_{L} sin ψ_{Q}]}{(1 - cos ψ_{Q} + q_{e} ψ_{Q} sin ψ_{Q}) (ψ_{Q} cot (ψ_{Q} / 2) - q_{e} ψ_{Q}^{2} - 2 K^{2}) + Q_{L}}

(10)

where:

Q_{L} = q_{L} [ψ_{Q} cos ψ_{Q} - (q_{e} ψ_{Q}^{2} + K^{2}) sin ψ_{Q}]

K^{2} = \frac{e_{26}^{2}}{ɛ_{22} {\bar{c}}_{66}}

q_{e} = \frac{ρ_{e}}{ρ_{Q} h_{Q}}

q_{L} = \frac{γ k_{L} {\tilde{μ}}_{L} tan ψ_{L}}{k_{Q} {\bar{c}}_{66}}

ψ_{Q} = k_{Q} h_{Q}

ψ_{L} = k_{L} h_{L}

{\tilde{μ}}_{L} = μ_{L} + i ω η_{L}

γ = \frac{r_{L}^{2} P (r_{L})}{r_{e}^{2} P (r_{e})}

P (r) = \frac{1}{π r^{2}} \int_{0}^{r} p (r) 2 π rdr

and ρ_e is the areal density of the electrode.

It can be shown that the admittance Y, defined as the ratio of amplitude of the current to that of the voltage applied across the electrodes, is given as:

Y = i ω C_{0} (1 - \frac{h_{Q} P (r_{e})}{2 {\overset{⌢}{ϕ}}_{0}} E)

(11)

where C₀=ε₂₂A_e/h_Q is the static capacitance of quartz and A_e is the area of the electrode. The admittance Y is composed of real part G and imaginary part B called conductance and susceptance, respectively.

The resonant frequency ƒ₀ here is defined as the frequency at which G becomes the maximum. For the case of no overlayer (q_L = 0) and under the assumption of K² ≪ 1 and q_e ≪ 1, the resonance frequency becomes:

f_{0} = f_{00} [1 - (4 K^{2} / π^{2} + 2 q_{e}) (1 - 4 q_{e} + 2 K^{2} q_{e})]

(12)

where:

f_{00} = 1 / (2 h_{Q} \sqrt{ρ_{Q} / {\bar{c}}_{66}})

(13)

It was shown by many investigators that the quartz crystal resonator can be understood in terms of an equivalent electrical circuit. Now, Equation (11) can be written as Y = Y₀+Y_m, where Y₀ = iωC₀ represents the admittance of the static capacitance C₀, and Y_m = –iωC₀h₀P_e/(2φ̂₀) denotes the admittance of motional branch. The impedance of the motional branch Z_m is the inverse of Y_m. Under the assumption that ψ_Q = k_Qh_Q = ψ(1-iξ), where ψ = ωh_Q√(ρ_Q/ c̄₆₆), is very close to π for ƒ ≅ ƒ₀ ≅ ƒ₀₀, it can be shown that the impedance becomes:

Z_{m} = \frac{1}{8 i ω c_{0} K^{2}} {(π^{2} - 8 K^{2}) - [(1 + 4 q_{e}) ψ^{2} + 2 q_{L R} ψ] + 2 i (ξ ψ^{2} + q_{L I} ψ)}

Note that ψ represents the magnitude of ψ_Q, and the ratio ξ is the imaginary part versus the real part of ψ_Q; ξ = ωη_Q /(2c̄₆₆). In this derivation, ξ is also assumed to be very small (ξ ≪ 1). Notice also that q_LR and –q_LI are the real and imaginary parts of q_L, respectively: q_L = q_LR – iq_LI. On the other hand the impedance Z_m is assumed to be composed of resistance R_m, capacitance C_m and inductance L_m:

Z_{m} = R_{m} + i ω L_{m} + 1 / (i ω C_{m})

Comparing these two, we get:

R_{m} = \frac{π ξ + q_{L I}}{8 C_{0} K^{2} f_{0}}

(14a)

L_{m} = \frac{1}{32 C_{0} K^{2} f_{0}^{2}} [(1 - 8 K^{2} / π^{2}) + 2 q_{L R} / π]

(14b)

C_{m} = \frac{8 C_{0} K^{2}}{π^{2} - 8 K^{2}}

(14c)

where ƒ₀ is given as Equation (12). All the three formulas shown above are the same as those derived by Martin et al. [14].

The resonant frequency with the overlayer can be obtained from the three components derived above; ƒ_0L=1/[2π√(L_mC_m)]. Therefore, the shift of resonant frequency Δ ƒ_0L = ƒ_0L – ƒ₀ becomes:

Δ f_{0 L} = - \frac{q_{L R}}{π} f_{0}

(15)

On the other hand, following the R-L-C circuit theory, G_max is given by:

G_{max} = \frac{1}{R_{m}} = \frac{8 C_{0} K^{2} f_{0}}{π ξ + q_{L I}}

(16)

In addition, G_maxΔƒ₁₂ is given as G_maxΔƒ₁₂ =ƒ_0L/(QR_m), where Q = 2πL_m/R_m is the Q-factor, Δƒ₁₂=ƒ₂ – ƒ₁ the band width, and ƒ₁ and ƒ₂ the values of the frequency at which the susceptance becomes the maximum and minimum, respectively. Then, we derive:

G_{max} Δ f_{12} = \frac{1}{2 π L_{m}} = \frac{16 K^{2} ɛ_{22} A_{e} f_{0}^{2}}{π h_{Q}}

(17)

This indicates that G_maxΔƒ₁₂ which is a measure of the static capacitance of the quartz is independent of the properties of the overlayer.

We can see from Equations (15) and (16) that two important parameters that reflect the properties of the material of the viscoelastic overlayer are Δƒ_0L and G_max; the magnitude of the former increases proportionally with the raised q_LR and the latter decreases with the raised q_LI. On the other hand, we have:

q_{L} = M exp (i θ / 2) tan [N exp (- i θ / 2)]

(18)

where:

M = \frac{γ \sqrt{ρ_{L} μ}}{\sqrt{ρ_{Q} {\bar{c}}_{66}}}, N = \frac{π h_{L} \sqrt{ρ_{L} / μ}}{h_{Q} \sqrt{ρ_{Q} / {\bar{c}}_{66}}}

(19a)

and:

μ = \sqrt{μ_{L}^{2} + {(ω_{0} η_{L})}^{2}}, θ = {tan}^{- 1} (ω_{0} η_{L} / μ_{L})

(19b)

in the range of θ being 0 ≤ θ ≤ π/2.

For an effective analysis, we need to introduce two group variables, ∏₁=h_L√(c̄₆₆ /(h_Q√μ)) and ∏₂ =ω₀η_L/μ_L. We note that the asymptotically thin film of the material provides ∏₁ → 0 and the asymptotically thick limit gives ∏₁ → ∞. Similarly, a very hard material means ∏₂ → 0 and the Newtonian liquid ∏₂ → ∞. Further, in the following analysis we assume that γ = O(1) and ρ_L ∼ ρ_Q.

First, for the case of thin film (∏₁ → 0), we get q_LR=q_Lc=π γ ρ_L h_L /(ρ_Q h_Q) and the formula for the frequency shift Δƒ_0L as follows:

Δ f_{0 L} = - γ f_{0} \frac{ρ_{L} h_{L}}{ρ_{Q} h_{Q}}

(20)

which corresponds to the Sauerbrey's formula [15] except for the parameter γ.

For the case of thick (∏₁ → ∞) and elastic (∏₂ → 0) solid, we have:

q_{L R} = M tan N

(21a)

q_{L I} = 0

(21b)

Since q_LI=0, G_max is the largest and independent of the material property. Further, Equation (21a) implies that there can be an infinite number of values of the parameter μ_L that produces a given value of q_LR, because at N = (2n-1) π /2 (n=1,2,3…) q_L becomes infinite. However, the case of thick elastic material may not be encountered so frequently.

For the case of thick layer (∏₁ → ∞) of Newtonian liquid (∏₂ → ∞), we have μ =ω₀ η_L. When the liquid has a low viscosity (η_L → 0), then the frequency shift Δƒ_0L can be given as:

Δ f_{0 L} = - γ f_{0}^{3 / 2} \frac{\sqrt{ρ_{L} η_{L}}}{\sqrt{π ρ_{Q} {\bar{c}}_{66}}}

(22)

which is the same as the one derived by Kanazawa and Gordon [16] except for the parameter γ. In case of the resonator application to the measurement in water solution, the density and viscosity of the solution affect the frequency shift. Muramatsu et al. [17] showed that the experimental results of the shift linearly increase with the raise of square root of solution density and viscosity in water and ethanol-water solution.

For the case of general viscoelastic overlayer, we can investigate the effect of μ_L and η_L on q_L by plotting q_L versus μ_L at different values of η_L. As shown in Figure 2, q_L oscillates with μ_L more pronouncedly at lower values of η_L. A lower value of μ_L produces the smaller amplitude of q_L fluctuation. It again implies that we can expect multiple values of the parameter set (μ_L, η_L) which give rise to a given value of q_L; selecting such parameter set is more highly possible at lower values of η_L and μ_L. From Equations (15) and (16) it is known that the frequency shift is determined from the real part of q_L and the maximum conductance from the imaginary part of q_L. When the experimental measurements of frequency shift and conductance are available, the estimation of μ_L and η_L is possible from the equations as demonstrated in Figure 2. In this case, we must pay attention to the possibility of multiple sets of the properties which has not been reported in the previous studies.

In the practical implementation of a resonator the viscoelastic overlayer is often loaded on the electrode surface. When a viscoelastic overlayer of polymer is applied to the resonator for the measurement of physical property variation along with a change of temperature or polymerization degree, the variation of viscosity and elastic shear modulus of the overlayer can be computed for the frequency shift and maximum conductance of the resonator. In other words, monitoring the frequency and the conductance can explain the polymerization and crystallization processes of a polymer sample. In this case multiplicity in the solutions of the equations may be resolved by a certain strategy. For instance, provided that the resonant property is known at a given parameter set, the viscoelastic property is readily obtained from the equations uniquely: When the experimental condition, such as temperature, varies continuously, the physical properties of overlayer also follow the condition and the determination of the properties from Figure 2 is a simple procedure.

3. Experimental

3.1. Materials

Polyethylene (Sigma-Aldrich Inc., St. Louis, MO, U.S.A.) having a number-average molecular weight of about 7,700 and a melting point of 90 °C was used as received. An AT-cut quartz crystal resonator having a base frequency of 8 MHz (Sunny Electronics Co., Korea) was utilized in this experiment. The electrodes of the resonator were silver finished.

3.2. Experimental Setup

A schematic diagram of a quartz crystal resonator is demonstrated in Figure 3(a). The polyethylene overlayer was placed on the one of electrode surfaces of the resonator, and the resonator was heated and cooled in an oil bath. Because the resonator surfaces can not be in contact with oil for the accurate measurements of its resonant frequency, conductance and susceptance, the resonator was placed in a specially designed module. The cell module illustrated in Figure 3(b) was prepared with two glass holders, two O-rings, two bakelite brackets and four screws. The dimension of the brackets is shown in Figure 3(c). The resonator was placed between two glass holders, and two o-rings keep the oil from wetting the electrode surface. Two brackets tighten the glass holders with four screws.

The detailed dimension of the cell is given in the figure. After the module was assembled, fine particles of polyethylene were obtained by sieving the powder with a sieve of 250 μm. About one third of 0.1 mg of the polyethylene powder was placed on the top electrode of the resonator. The module was put in an oil bath demonstrated in Figure 4. For the better control of resonator temperature, the module was immersed to the level of the upper o-ring in the bath. The bath temperature was adjusted by electric heating, and water cooling was also provided for the cooling cycle in the experiment. The electric heater was controlled with a programmable temperature controller (Hanyoung Electronic Co., Korea, Model NP-200). The resonant frequency, resonant resistance and the temperature of oil bathwere measured using home-made devices, and an A/D converter was employed for signal processing. The digital signals of the resonant frequency, resonant resistance and temperature were provided to a PC for the data analysis.

3.3. Experimental Procedures

The experiment began at a temperature of about 25 °C. After the experimental setup was stabilized for an hour, the bath temperature was raised at a rate of 1 °C/min up to 100 °C, and was lowered at a rate of 1.5 °C/min. The first measurement was conducted at a temperature of 95 °C. At the temperature the setup was steadied for two minutes, and the measurement was conducted for 8 minutes using an impedance analyzer (Agilent Technologies, U.S.A., Model 4192A). The conuctance and susceptance were determined at a resonant frequency between 794 MHz and 804 MHz in a step of 50 Hz. The temperature was lowered by 4 °C down to 55 °C, and the measurement was done in the same manner

4. Application to Polyethylene Crystallization

The formula and analysis presented so far are applied to the experimental study on the measurement of material properties of polyethylene during its crystallization process. As the material and geometrical constant for the quartz, we consider the following parameters many of which are commonly used for an AT-cut crystal:

c_{66} = 2.947 \times 10^{11} dyn / c m^{2}

ɛ_{22} = 40 \times 10^{- 14} F / c m

ρ_{Q} = 2.651 g / c m^{3}

K^{2} = 7.74 \times 10^{- 3}

h_{Q} = 0.0205449 c m

Here, the depth of quartz h_Q was obtained in such a way that the resultant resonant frequency matches with the measured one ƒ₀ = 7,996 [kHz] without any overlayer. Further, we consider the following properties for silver electrodes:

h_{e} = 30 n m

ρ_{e} = 10.5 \times 3 \times 10^{- 5} g / c m^{2}

Then, we have q_e = 0.00594. On the other hand, Equation (13) gives ƒ₀₀ = 8,114 kHz and we obtain the theoretical prediction for ƒ₀, the resonant frequency without the overlayer, from Equation (12) as follows: ƒ₀ = 7,997.7 kHz. On the other hand, the computation with Equation (11) gives ƒ₀ = 7,996.0 kHz. We can see that the error of the theoretical prediction is only 0.02%. The electrode diameter d_e takes different value for each experimental measurement to match the measured G_maxΔƒ₁₂ using Equation (17).

The experiment provided the conductance G and the susceptance B as a function of the frequency ƒ for the case without overlayer at the room temperature and 11 cases with melted polyethylene at different temperatures, and a typical set of data at a temperature of 83 °C is shown in Figure 5(a). First, for the unloaded case, we read G_max from the G curve of Figure 6(a) or from the radius of a circle matching the Nyquist plot as shown in Figure 5(b), and also read Δƒ₁₂ from the B curve in Figure 5(a). From Equation (17) we calculate the electrode diameter; d_e =0.5101cm. Then we find ξ from matching G_max obtained from the experiment with that given by calculating (16). Equation (17) can be used as a guide or initial value for ξ, where R_m =1/G_max is applied. Equation (17) yields ξ = 2.08 × 10^-4, whereas the numerical calculation gives ξ = 2.12 × 10^-4 being only 2% different from the asymptotic calculation.

We next obtain the parameters μ_L and η_L for loaded cases. The measurement was done with the PE particles melted on the electrode surface. The total mass of the particles is 0.033 mg. It is assumed that all the particles contribute to a single particle, its size being determined by summing up all the particles' size; d_L=0.0502 cm. This in turn provides the thickness of the PE particle; h_L=0.0184 cm. Since the melted particles do not change their shape with temperature, d_L and h_L are fixed in the subsequent calculations. It is true that a parameter set (μ_L, η_L) results in one set of parameters q_LR and q_LI and therefore one set of values G_max and Δƒ. However as analyzed in the previous section, the inverse is not true. In fact, infinite number of parameter set (μ_L, η_L) can result in a given parameter set G_max and Δƒ. Figure 6 presents G_max and Δƒ measured from the experiment at 11 different temperatures. We see that G_max decreases monotonously with the temperature decrease, but Δƒ shows fluctuation although in overall it also decreases with the temperature decrease. The sensitivity of the experimental setup is determined by the measurement resolution of the impedance analyzer. The resolution of frequency measurement is 1 Hz, and that of conductance is 1 nS. Therefore, the variation trends of the frequency difference and conductance in Figure 6 are real, while the fluctuation in the frequency is due to the instrumental sensitivity. To find the parameter set (μ_L, η_L) corresponding to each of the 11 data set (G_max, Δƒ), we generate a map on the (μ_L, η_L) space as shown in Figure 7. Intersection of contour lines of G_max and that of Δƒ then gives us the parameters μ_L and η_L. Figure 8 exhibits the parameters obtained in this way at the fundamental mode of the overlayer vibration. We see that as the temperature decreases from the melting point, the shear modulus decreases but the viscosity increases.

The measurements of elastic shear modulus and viscosity shown in Figure 5 are compared with those of bulk polyethylene in melt state [18], and the results are listed in Table 1. The top three rows are of the instrumental measurements in the low rate of shear, the estimated values of this study is followed, and the bottom row is the estimation for the high frequency viscoelasticity from the low frequency measurements. Because the melt properties are of a large difference of frequency, a direct comparison between the measurement and the estimation of this study is difficult. However, the two rows show comparable values of viscoelasticity. The tendencies of the increase of elastic shear modulus and the decrease of viscosity with the frequency elevation of sample movement show that the measurements of this study are same to the melt properties. In addition, the effects of temperature variation on the elastic shear modulus and viscosity are the same in the studies. While the instrumental measurement of rheological property is limited to polymer melt, the proposed measurement technique of this study can be implemented to the wide range of applications for various materials.

5. Conclusions

A generalized relationship between the viscoelastic property of a polymer overlayer placed on the electrode of a quartz crystal resonator and its resonant charateristic is developed from the mechanics of the quartz movement and applied to the determination of the viscoelastic properties during the polymer crystallization. The developed procedure is utilized in polyethylene crystallization, and the estimated viscosity and elastic shear modulus from the measurements of resonant frequency and conductance of the resonator are compared with those of a polyethylene melt determined instrumentally. The comparison indicates that the estimation is comparable to that of the melt, implying the possibility of wide application of the proposed technique in the field of polymer processing.

Acknowledgments

Financial support from the National Research Foundation of Korea through the National Research Laboratory Program (2005-1091) and the Basic Science Research Program (2009-72473) is gratefully acknowledged.

References

Choi, H.S.; Kim, Y.H.; Shin, S.M. In-Line monitoring of uv photopolymerization with a quartz crystal resonator. J. Polym. Sci. Chem. 2006, 44, 2428–2439. [Google Scholar]
Kim, B.C.; Kim, Y.H.; Fukui, K. Crystallization monitoring in supersaturated solution with a quartz crystal sensor. Anal. Chim. Acta 2003, 491, 71–80. [Google Scholar]
Maeda, K.; Kim, B.C.; Kim, Y.H.; Fukui, K. Nucleation monitoring in cooling crystallization with a quartz crystal oscillator. Anal. Chim. Acta 2006, 558, 337–344. [Google Scholar]
Reed, C.E.; Kanazawa, K.K.; Kaufman, J.H. Physical description of a viscoelastically loaded AT-cut quartz resonator. J. Appl. Phys. 1990, 68, 1993–2001. [Google Scholar]
Buttry, D.A.; Ward, M.D. Measurement of interfacial processes at electrode surfaces with the electrochemical quartz crystal microbalance. Chem. Rev. 1992, 92, 1355–1379. [Google Scholar]
Nwankwo, E.; Durning, C.J. Fluid property investigation by impedance characterization of quartz crystal resonators: Part II: Parasitic effects, viscoelastic fluids. Sens. Actuat. 1999, 72, 195–202. [Google Scholar]
Kim, J.-M.; Chang, S.-M.; Muramatsu, H. Monitoring changes in the viscoelastic properties of thin polymer films by the quartz crystal resonator. Polymer 1999, 40, 3291–3299. [Google Scholar]
Kunze, M.; Shull, K.R.; Johannsmann, D. Quartz crystal microbalance studies of the contact between soft, viscoelastic solids. Langmuir 2006, 22, 169–173. [Google Scholar]
Efimov, I.; Hillman, A.R.; Schultze, J.W. Sensitivity variation of the electrochemical quartz crystal microbalance in response to energy trapping. Electrochim. Acta 2006, 51, 2572–2577. [Google Scholar]
Nakano, Y.; Sato, I.; Seida, Y.; Nakano, Y. Viscoelastic behavior of thermoresponsive polymer hydrogel with organic adsorbate observed using quartz crystal microbalance. Chem. Lett. 2007, 36, 1204–1205. [Google Scholar]
Zhu, D.M.; Fang, J.J.; Wu, B.; Du, X.B. Viscoelastic response and profile of adsorbed molecules probed by quartz crystal microbalance. Phys. Rev. E 2008, 77, 031605. [Google Scholar]
Johannsmann, D. Viscoelastic, mechanical, and dielectric measurements on complex samples with the quartz crystal microbalance. Phys. Chem. Chem. Phys. 2008, 10, 4516–4534. [Google Scholar]
Zhang, Y.; Du, B.; Chen, X.; Ma, H.W. Convergence of dissipation and impedance analysis of quartz crystal microbalance studies. Anal. Chem. 2009, 81, 642–648. [Google Scholar]
Martin, S.J.; Granstaff, V.E.; Frye, G.C. Characterization of a quartz crystal microbalance with simultaneous mass and liquid loading. Anal. Chem. 1991, 63, 2272–2281. [Google Scholar]
Ward, M.D.; Delawski, E.J. Radial mass sensitivity of the quartz crystal microbalance in liquid media. Anal. Chem. 1991, 63, 886–890. [Google Scholar]
Kanazawa, K.K.; Gordon, J.G. II. Frequency of a quartz microbalance in contact with liqui. Anal. Chem. 1985, 57, 1770–1771. [Google Scholar]
Muramatsu, H.; Tamiya, E.; Karube, I. Computation of equivalent circuit parameters of quartz crystals in contact with liquids and study of liquid properties. Anal. Chem. 1988, 60, 2142–2146. [Google Scholar]
Xanthos, M.; Tan, V.; Ponnusamy, A. Measurement of melt viscoelastic properties of polyethylenes and their blends-a comparison of experimental techniques. Polymer Eng. Sci. 1997, 37, 1102–1112. [Google Scholar]

Figure 1. Sketch of a quartz crystal resonator with electrodes on both sides and a viscoelastic overlayer attached on the external surface of an electrode.

Figure 2. Sketch of a quartz crystal resonator with electrodes on both sides and a viscoelastic overlayer attached on the external surface of an electrode.

Figure 3. Schematic diagrams of quartz crystal resonator (a), resonator cell module (b) and holding bracket (c).

Figure 4. A schematic diagram of experimental setup.

Figure 5. Comparison between calculated (lines) and measured (symbols) conductance (G) and susceptance (B) at T = 83 [°C] (a) in the ƒ – G and ƒ – B planes; (b) in the G – B plane (Nyquist plot).

Figure 6. Experimental results of Δƒ and G_max obtained at 11 different temperatures.

Figure 7. Typical map for determining μ_L and η_L of a viscoelastic overlayer. Solid lines are contour lines for Δƒ = −740[Hz] and dashed lines for G_max = 0.0013[S] at 83 °C.

Figure 8. Estimated values of μ_L and η_L of a viscoelastic overlayer at each temperature with fundamental mode of oscillation. The temperature starts from 95 °C decreasing down to 55 °C with 4 °C increment.

Table 1. Experimental and estimated values of viscosity and elastic shear modulus.

**Table 1.** Experimental and estimated values of viscosity and elastic shear modulus.
Frequency (Hz)	Viscosity (dyn s/cm²)	Shear modulus (dyn/cm²)
Bulk [18]
0.1	2 × 10⁵	1.5 × 10⁴
1	1.1 × 10⁵	1 × 10⁵
4	6 × 10⁴	1.5 × 10⁵
Overlayer
8 × 10⁶	185	7.54 × 10¹⁰
From bulk data
8 × 10⁶	346	6.47 × 10⁹

© 2009 by the authors; licensee Molecular Diversity Preservation International, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).

Share and Cite

MDPI and ACS Style

Suh, Y.K.; Kim, B.C.; Kim, Y.H. Determination of Viscoelastic Property in Polyethylene Crystallization Using a Quartz Crystal Resonator. Sensors 2009, 9, 9544-9558. https://doi.org/10.3390/s91209544

AMA Style

Suh YK, Kim BC, Kim YH. Determination of Viscoelastic Property in Polyethylene Crystallization Using a Quartz Crystal Resonator. Sensors. 2009; 9(12):9544-9558. https://doi.org/10.3390/s91209544

Chicago/Turabian Style

Suh, Yong Kweon, Byoung Chul Kim, and Young Han Kim. 2009. "Determination of Viscoelastic Property in Polyethylene Crystallization Using a Quartz Crystal Resonator" Sensors 9, no. 12: 9544-9558. https://doi.org/10.3390/s91209544

APA Style

Suh, Y. K., Kim, B. C., & Kim, Y. H. (2009). Determination of Viscoelastic Property in Polyethylene Crystallization Using a Quartz Crystal Resonator. Sensors, 9(12), 9544-9558. https://doi.org/10.3390/s91209544

Article Menu

Determination of Viscoelastic Property in Polyethylene Crystallization Using a Quartz Crystal Resonator

Abstract

1. Introduction

2. Theoretical Analysis

3. Experimental

3.1. Materials

3.2. Experimental Setup

3.3. Experimental Procedures

4. Application to Polyethylene Crystallization

5. Conclusions

Acknowledgments

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.