online measurement of diameter and temperature in the melt-spinning and melt blowing process
- 格式:pdf
- 大小:902.48 KB
- 文档页数:9
Online Measurement of Fiber Diameter and Temperature in the Melt-Spinning and Melt-Blowing ProcessesVishnu T.Marla,Robert L.Shambaugh,*and Dimitrios V.PapavassiliouSchool of Chemical,Biological,and Materials Engineering,The Uni V ersity of Oklahoma,100East Boyd Street,SEC T335,Norman,Oklahoma 73019Online measurements of fiber temperature and diameter were made for both the melt-spinning and the melt-blowing processes.The fiber temperature was determined by infrared thermography,and the fiber diameter was determined by high-speed photography.These measurements were then compared with predictions made with mathematical models for melt spinning and melt blowing.There was good agreement between the models and the experimental results,and the agreement was best when heat-transfer correlations developed specifically for fine fibers (cylinders)were used.IntroductionMelt spinning and melt blowing are widely used commercial processes for the manufacture of fibers and nonwovens.For these two processes,experimental measurements of fiber properties along the spinline are useful for understanding the process of fiber formation.Furthermore,these measurements can serve as valuable data for testing against the predictions of mathematical models for the processes.In melt spinning fiber motion is primarily along the spinning direction (1D motion),while in melt blowing the fiber has motion in the spinning direction as well as extensive vibrational motion (so that the fiber has 3D motion).However,the processes are similar in that both have a rapidly attenuating fiber.In the work described herein,experimental data were compared with the predictions of previously developed mathematical models for melt spinning and melt blowing.Some past studies have focused on the online measurement of fiber diameter,fiber temperature,fiber velocity,birefringence,and structural changes along the spinline.Online fiber velocity was measured using laser Doppler velocimetry.1,2Online birefringence was measured with a polarizing microscope 3-5or crossed polarizers.6Online X-ray characterization of the fiber during melt spinning was used to monitor structural changes along the threadline.7Fiber diameter is a key parameter in fiber spinning,and the online measurement of this parameter has been attempted in a number of ways.For example,Lu and Spruiell 3used a polarizing microscope that was mounted sideways such that the running fiber could be seen through the eyepiece of the microscope.By using a calibrated micrometer eyepiece,the fiber diameter could be measured.Bheda and Spruiell 8measured the fiber diameter in melt spinning using a Zimmer diameter monitor based on the principle of electrooptic back illumination.Their setup permitted several measurements at each spinline position,and a statistical analysis was conducted to determine the average diameter.A similar setup for measuring the fiber diameter profile was used by several other research groups.9,10Haberkorn et al.11used a forward laser light scattering (LLS)technique to determine the fiber diameter along the spinline.The fiber diameter was determined from the distance between the interfer-ence maxima on the diffraction pattern that was recorded on a charged couple device (CCD camera).Wallen et al.12used smallangle light scattering to study the fiber diameter and orientation in the melt-blowing process.Bansal and Shambaugh 4,5used high-speed photography for measuring the fiber diameter in both the melt-spinning and melt-blowing processes.Oh et al.13obtained the diameter profile for solid and hollow fibers in the melt-spinning process by making online measurements using a special capturing device that trapped and froze the filament.The captured filaments were cut with a microtome,and the inner and outer fiber diameters were measured with a microscope.Bresee and Qureshi 14used a pulsed laser in conjunction with a rapid frame rate camera to acquire the fiber diameter in the melt-blowing process.Zhao and Wadsworth 15used a 51cm long rigid mechanical arm to collect fiber at different locations along the melt-blowing spinline and analyzed the collected fibers using an optical microscope and image analyzing software.Recently,Golzar et al.16used an infrared camera to measure the fiber diameter in the melt-spinning process.However,due to resolu-tion limitations of the infrared camera,diameter measurements below 25µm were not possible.Moore et al.17used an ensemble laser diffraction (ELD)technique to determine the fiber diameter distribution at a location below a melt-blowing die.One method of measuring the fiber temperature profile is the use of infrared radiation.An infrared technique allows for noncontact measurement of the fiber temperature.For example,Lu and Spruiell 3used a Barnes infrared microscope to measure the temperature of the filament in the melt-spinning process.They used a null-balance technique where the temperature of the filament was compared with that of a small heater at a known temperature.The radiation from the filament was compared with that of the heater,and the filament temperature was assumed to be the same as the heater temperature at the point where the filament radiation matched the radiation from the heater.Apparently,the assumption was made that the emissivity of the heater was the same as that of the filament.Because an IR microscope has a depth of field on the order of a few micrometers,an IR microscope cannot be used to measure fiber temperatures during melt blowing (where the fiber vibrational motion is much greater than a few micrometers).For temperature measurements on fibers during melt blowing,Bresee and Ko 18used a digital IR thermometer with adjustable emissivity.They inserted the thermometer probe into the fiber stream and measured the temperature of the fibers as the fibers collected on the probe surface.However,since this technique requires insertion of a probe in the polymer stream,the technique cannot be applied to melt spinning.Bansal and Shambaugh 4,5*To whom correspondence should be addressed.Tel.:(405)325-6070.Fax:(405)325-5813.E-mail:shambaugh@.Ind.Eng.Chem.Res.2009,48,8736–8744873610.1021/ie900615n CCC:$40.75 2009American Chemical SocietyPublished on Web 08/20/2009and de Rovere and Shambaugh19used a slit response factor or SRF curve to account for the spatial resolution of their IR camera while measuring the temperature offibers in the melt-spinning and melt-blowing processes.The SRF curve corrects for the errors involved in measuring the temperature offine fibers(where,on the view frame of the camera,thefiber diameter has a size that is on the order of a pixel).In the new work described herein,the generic SRF curve provided by the manufacturer of the camera(FLIR)was replaced by a more accurate temperature correction scheme that was specifically tailored to thefine diameterfibers used in our work.Errors in IR measurement can also be attributable to incorrect estimation of object emissivity.Emissivity values are listed in the literature,and values are also provided by manufacturers of IR cameras.However,emissivity has been known to be a function of the thickness in the case of thin polymerfilms,20 the diameter in the case of cylindrical objects,21and the angle of vision.Since the work of Marla et al.21is directly applicable to the type of polymerfilaments produced during melt spinning and melt blowing,the calibration techniques suggested in their paper are an accurate way of measuringfiber temperature in these spinning systems.Both melt spinning and melt blowing involve heat transfer from a polymerfilament to an air stream thatflows across the filament.There are several correlations for heat transfer from relatively large diameter cylinders(wires)to air in parallel flow.22-24The heat-transfer correlation of Kase and Matsuo25 is commonly used in mathematical models for melt spinning and melt blowing;this correlation has a range of applicability that is appropriate for smallfiber diameters.However,use of the Kase and Matsuo correlation might be the cause of overprediction of the temperature(when compared with ex-perimental data).For melt spinning,an example of this overprediction can be found in de Rovere and Shambaugh.19 For melt blowing,see the article by Bansal and Shambaugh.5 Recently,new Nusselt number correlations were developed for heat transfer fromfine cylinders.26These correlations were produced by carefully measuring the change infiber temperature with a calibrated infrared camera.The correlations that were developed from this recent experimental study areandFor comparison,the correlation of Kase and Matsuo25isAs one of the goals of this paper,we evaluated the performance of the newly developed correlations by inserting them into mathematical models for melt spinning and melt blowing and comparing the predicted temperature profiles with online temperature measurements.Details of Experimental SetupMelt Blowing.Figure1shows the experimental setup for measurement of thefiber diameter and temperature in the melt-blowing process.The setup was very similar to that used by Marla and Shambaugh.27A Brabender extruder of19.1 mm(0.75in.)diameter and381mm length was used to melt and pressurize the polymer.The barrel had a20:1L/D ratio and a3:1compression ratio.The polymer exiting from the extruder was then fed to a modified Zenith pump that accurately metered molten polymer through a melt-blowing slot die with a single polymer capillary.The polymer capillary had a diameter of0.420mm.The two air slots were setflush with the nose piece,and the slot widths were0.65mm.Each slot had a length of74.6mm(2.94in.).The die assembly was heated with two250W cartridge heaters.A thermal mass flow meter was used to measure the airflow rate,which was maintained at either100or125slpm(at standard conditions of21°C and atmospheric pressure).The centerline air velocity and temperature below the die were determined using a pitot tube and thermocouple probe,respectively.The details of this equipment can be found in Marla et al.21Thefiber diameter profiles were obtained using high-speedflash photography;see the details in ref5.The camera used was a Nikon N90S equipped with a105mm Nikon macro lens, and the illumination was provided by a Sunpak Auto622flash.The polypropylene(PP)used in the experiments was 88MFR Fina Dypro isotactic polypropylene.This polymer had an M w of165000g/mol and M n of41500g/mol. Polybutylene(PB)resin was also spun.This polybutylene was grade0400manufactured by Basell Polyolefins;the polybutylene MFR was20.The temperature of the polymerfilaments was determined using an IR(infrared)camera.The IR camera used in the study was a ThermaCAM S60manufactured by FLIR Systems.The spectral range for this camera is7.5-13µm(longwave).The IR image displayed had a resolution of640×480pixels,while the actual number of thermal detector sensors or elements was 320×240.The other prominent features of this IR camera can be found on the manufacturer’s Web site and the camera’s user manual(FLIR ThermaCAM S60User Manual).The IR lens on this camera has afield of view(FOV)of24°×18°with a minimum focus distance of0.3m or11.81in.The camera also had a100µm close-up lens for which the minimum focus distance is80mm,while the maximum working distance is 110mm.The close-up lens had the advantage that it could measure the temperature of small targets with greater accuracy. Unless otherwise mentioned,all temperature measurements in the present study were made with the100µm close-up lens. The IR camera had the capability to record60temperature measurements per second,and a real-time plot of theseNu)1.4+0.0141Re for24e Re e150(1) Nu)0.3575Re0.46for150e Re e923(2) Nu)0.42Re0.334(3)Figure1.Experimental equipment used for spinningfibers while measuring fiber diameter and temperature.Ind.Eng.Chem.Res.,Vol.48,No.18,20098737measurements was displayed on the monitor of a PC that was connected to the IR camera by means of an IEEE1394(firewire) connection.Since the100µm close-up lens has a small depth offield and melt-blownfibers exhibit significant motion in a plane perpendicular to the majorfiber axis(the spinning direction),thefiber moved in and out of focus of the IR camera. Whenever thefiber was in focus,the temperature plot displayed a spike in thefiber temperature.This spike was then taken as the rawfiber temperature(i.e.,a temperature in need of correction forfield of view,emisssivity,etc.).Beginning at3 cm below the die,measurements were taken at intervals of1 cm until a position13cm below the die was reached.Beyond 13cm,thefiber temperature was measured at15,17,20,and 25cm below the die.Close to the die,infrared radiation from the die itself can affect the measuredfiber temperature.Thus, no measurements were made closer than3cm to the die.The variation of emissivity withfiber diameter(for both polypro-pylene and polybutylene),the procedure for using the SRF curve,and the operating procedure of the IR camera are described elsewhere.21The base conditions for the experiments, unless otherwise specified,were m p)0.75g/min,T f,die)310°C,T a,die)368°C,and airflow rate)125slpm.Melt Spinning.For melt spinning,the experimental setup was the same as that for melt blowing except that the air supply was turned off and the melt-blowing die was replaced with a single-hole spinneret.The capillary in this spinneret had a diameter of0.407mm(0.016in)and a length of2.97mm (0.1169in).Both polypropylene and polybutylenefibers were spun.The molten polymer exiting the die was drawn by a mechanical take-up roll located1.4m below the spinneret.For experiments with both polypropylene(PP)and polybutylene (PB),the spinneret temperature was kept at200°C,and the polymer throughput was kept at1g/min.For polypropylene, the mechanical take-up roll was used at six spinning speeds ranging from500to1750m/min;for polybutylene the roll speeds were500and1000m/min.Thefiber diameter and temperature were measured in the same way as in melt blowing. Beginning from3cm below the die,the temperature measure-ments were made at1cm intervals up to12cm below the die. Beyond12cm,the measurements were made at15,17.5,and 20cm below the die.Modeling DetailsMelt Spinning.The energy balance over a differential element of size∆y of thefiber at any distance y below the spinneret can be written aswhere m p)mass throughput of the polymer(kg/s),C p,f) specific heat of the polymer(J/kg·K),d y)fiber diameter at any distance y below the spinneret(m),h)convective heat-transfer coefficient(W/m2·K),T a)temperature of the ambient air(°C),and T f)fiber temperature(°C).At very high spinning speeds(in excess of4000m/min),the stress in the polypropylenefiber increases enough to induce crystallization in thefiber.4Equation4neglects the heat gained by the polymer due to crystallization,and for the moderate spinning speeds in our experiments,this is a reasonable assumption.Thefiber diameter profile needed for solving eq4 was obtained using high-speedflash photography.For the conditions of our experiments,the specific heat capacity of both polypropylene and polybutylene do not vary significantly with temperature and can be assumed to be constant.Then,eq4is an ordinary differential equation with the temperature of the fiber,T f,as the only dependent variable.This ODE can be numerically solved for T f using the fourth-order Runge-Kutta method.By using different correlations for calculating the convective heat-transfer coefficient in eq4,the predictedfiber temperature can be compared with the results from online temperature measurements.Melt Blowing.For melt blowing,the model of Marla and Shambaugh28was used to predict the three-dimensional motion of afiber as it was being formed below a melt-blowing die. This model involves the simultaneous solution of the momen-tum,energy,and continuity equations.Predicted parameters includefiber attenuation,vibration frequency,vibration ampli-tude,temperature,and stress.The model was solved using both the Kase and Matsuo heat-transfer correlation and the correlation developed by Marla et al.26Results and DiscussionMelt Spinning.High-speedflash photography was used to measure the onlinefiber diameter.Figure2shows some representa-tive plots of online diameter profiles for polypropylene spun at 500and1750m/min.Thefibers attenuate rapidly close to the spinneret,and thefibers spun at1750m/min attenuate faster than thefibers spun at500m/min.The diameter profiles werefit to fourth-order polynomials so that they could be inserted in eq4as a function of distance y from the spinneret.Table1shows the coefficients of the fourth-order polynomial for each spinning speed for both PP and PB.Figure3shows the results from the online temperature measurements for PP at six different spinning speeds.Thefiber temperatures were determined with the aforementioned IR camera along with the calibration procedure developed in ref 21.Higher spinning speeds result infiner diameters that possess less sensible heat,andfinerfibers cool at a faster rate than thickerfibers.However,in the spinning process thesefinefibers are exposed to the ambient air for less time than the thicker filaments(because,at a constant polymer throughput,thicker filaments have a slower spinning speed).These two effects tend to counterbalance each other such thatfiber temperature shows little dependence on spinning speed.In fact,previous studies4,8 showed that the take-up speed did not significantly affect them p Cp,fd Tfd y)-πhdy(Tf-Ta)(4)Figure2.Fiber diameter as a function of position below the die for themelt spinning of polypropylene.Data are shown for two different windupspeeds.8738Ind.Eng.Chem.Res.,Vol.48,No.18,2009fiber temperature profile.As Figure 3shows,in our present experiments it appears that the spinning speed affects the fiber temperature (the effect is of the order of 10-25°C).For example,at 20cm below the die,the temperature of the fiber spun at 1750m/min is 90°C,while for fiber spun at 500m/min the temperature is 115°C.For the case of PB spinning,Figure 4shows the online temperature measurements at two different spinning speeds.The trends for PB parallel the trends observed for PP.Statistical examination of the data in Figures 3and 4indicates that the standard deviation of the data points is about 10.8°C.Figure 5is a plot of fiber temperature versus position for PP spinning at a windup speed of 500m/min.The diamond symbols are the experimental values of the temperature.These temper-ature values were corrected for fiber diameter (i.e.,corrected for emissivity changes)by using the procedure described by Marla et al.21The solid and dotted lines are predicted temper-ature profiles obtained by solving eq 4.For the solid line,the correlation of Marla et al.26was used in eq 4.For the dotted line,the correlation of Kase and Matsuo 25was used.The temperature profile obtained by using the Kase and Matsuo heat-transfer correlation was high compared to the experimental data.However,the temperature profile obtained by using the cor-relation developed by Marla et pared favorably with experimental data.Figures 6,7,8,9,and 10are similar to Figure 5except that these five figures are for spinning speeds of,respectively,800,1000,1200,1500,and 1750m/min.As was the case for Figure 5,these five figures show that there is a good match between the experimental data and the model predictions made with the correlation of Marla et al.26Also,use of the Kase and Matsuo correlation gives temperature profiles that are higher than the actual temperature measurements.Figures 11and 12are temperature measurements and predictions for the spinning of PB at,respectively,500and 1000m/min.As was the case for PP spinning (Figures 5-10),there is a good match between the experimental data and model predictions made using the correlations 1and 2.As with PP spinning,the use of the Kase and Matsuo correlation gave temperature profiles that were higher than the actual temperature profile.Table 1.Coefficients of the Fourth-Order Polynomial Used to Fit the Diameter Profiles Obtained from High-Speed Flash Photography during the Melt-Spinning Processcoefficients of the polynomial ay 4+by 3+cy 2+dy +epolymer spinning speed (m/min)a b c d e polypropylene500-1.776×10-14-0.143 6.775-112.350878.160800 4.867×10-30.40711.509-145.898938.571000 2.012×10-2 1.03819.707-179.945973.0641200 1.714×10-2-0.91418.146-172.923942.9201500 1.037×10-2-0.60713.674-152.677913.79717509.644×10-3-0.55412.934-150.273892.569polybutylene5000.023-1.07318.297-160.111954.90410000.0070.3728.286-108.846872.800Figure 3.Fiber temperature as a function of position below the die for the melt spinnng of polypropylene.Data are shown for six different windupspeeds.Figure 4.Results for polybutylene melt spinning:fiber temperature as a function of position below the die.Data are shown for two different windupspeeds.Figure 5.Model predictions of fiber temperature for the melt spinning of polypropylene.Two different correlations for Nusselt number were tested in the model.Experimental data are also included for comparison.Ind.Eng.Chem.Res.,Vol.48,No.18,20098739Melt Blowing.To solve the equations that describe melt blowing,air velocity and air temperature are needed as boundary conditions.Figures 13and 14show the centerline air velocityand temperature,respectively,at positions below the melt-blowing die for the conditions shown on the figures.These velocity and temperature measurements were made in the absence of the polymer,i.e.,the polymer flow was stopped during these measurements.Several other researchers 29-32have made measurements of the air velocity and temperature below slot dies in different configurations and developed dimensionless empirical correlations for the velocity and temperature fields.Recently,computational fluid dynamics (CFD)has been used to investigate the flow fields below melt-blowing slot dies.33-35As Figure 13shows,the centerline air velocity decreases along the threadline with the velocity being slightly higher for the 125slpm case compared to the 100slpm case.As Figure 14shows,the centerline air temperature profile shows a trend similar to the air velocity profile.However,the reduction in the air flow rate from 125to 100slpm does not seem to affect the air temperature.As mentioned above,the measurements of air velocity and temperature serve as boundary conditions in mathematical models for melt blowing.27,28,36In the model results shown in this paper,air field measurements were used as boundary conditions.Figure 15shows online,experimental measurements of PP fiber diameter profiles for different polymer throughputs(whileFigure 6.Similar to Figure 5,except for a spinning speed of 800m/min.Figure 7.Similar to Figure 5,except for a spinning speed of 1000m/min.Figure 8.Similar to Figure 5,except for a spinning speed of 1200m/min.Figure 9.Similar to Figure 5,except for a spinning speed of 1500m/min.Figure 10.Similar to Figure 5,except for a spinning speed of 1750m/min.8740Ind.Eng.Chem.Res.,Vol.48,No.18,2009keeping the base conditions constant).High-speed flash pho-tography was used for these measurements.Figure 15shows that,for a higher mass throughput,the fiber diameter is larger at any position below the melt-blowing die.Figure 16shows the diameter predictions from the Marla and Shambaugh model 28for the base operating conditions.The heat-transfer correlations of both Kase and Matsuo 25and Marla et al.26were used in the model.The diameter profile from online measurements (a subset from Figure 15)is also plotted on the figure.For the region of rapid fiber attenuation,the first 5cm below the die,the choice of heat-transfer correlation does not greatly affect the model predictions of fiber diameter profile,and there is good agreement between experimental data and model predictions.Uyttendaele and Shambaugh 37obtained a similar result with an annular melt-blowing die.The choice of heat-transfer correlation is more evident in the prediction of the final fiber diameter.Figure 16lists the final diameters obtained by using the different heat-transfer correla-tions in the model.The final diameter obtained by an off-line experimental measurement of fiber diameter is also listed.As can be seen,use of the Marla et al.correlation gives a better estimate of final fiber diameter.Figure 17shows experimental measurements of fiber diameter at two different air flow rates.The air flow rate does not significantly affect the diameter profile,and this result parallels the findings from the experimental and modeling work of Marla and Shambaugh 28with slot dies.Figure 18shows the experi-mental temperature profiles obtained from online measurements at different polymer throughputs.The fiber temperature profiles show that larger fibers cool more slowly,and increased mass throughput increases the fiber temperature at any spinline position.At any spinline location,there are two effects which counter one another due to an increase in the mass flow rate of the polymer.The first effect involves the thermal inertia (mC p,f )of the fibers.With an increase in the polymer flow rate,the thermal inertia increases due to the increase in diameter (mass).The second effect is the increase of the fiber diameter and,hence,the increase of surface area available for heat transfer.Our experimental data indicate that the first effect is somewhat larger than the second for the heat-transfer process.Similar results were obtained by Bresee and Ko.18For the lowest polymer rate shown in Figure 18(0.75g/min),Figure 19shows predictions of the Marla model for both oftheFigure 11.Polybutylene results:model predictions of fiber temperature for melt spinning.Two different correlations for Nusselt number were tested in the model.Experimental data are also included forcomparison.Figure 12.Similar to Figure 11,except for a spinning speed of 1000m/min.Figure 13.Centerline air velocity below the melt-blowing die.The 60°slot die had a slot length of 7.46cm and an air gap of 0.65mm,and the die had a single polymer orifice.These data,which were taken with a Pitot tube,serve as boundary conditions in the melt-blowingmodel.Figure 14.Centerline air temperatue below the melt-blowing die.These data,which were taken with a fine thermocouple,serve as boundary conditions in the melt-blowing model.Ind.Eng.Chem.Res.,Vol.48,No.18,20098741heat-transfer correlations.As was observed previously,the model predictions made with the Kase and Matsuo correlation were higher than when the correlation of Marla et al.was used.Probably due to the modeling complexities for the melt-blowing process,the predicted temperature is not as close to the actual measured temperature compared to the closeness (of the predicted results to the experimental results)for the process of melt spinning.However,even for melt blowing,the correlation of Marla et al.is definitely more accurate than the correlation of Kase and Matsuo.For example,at 25cm below the die,the model predictions based on the heat-transfer correlation of Marla et al.is 118°C,while the prediction made with the Kase and Matsuo correlation is 190°C.At the same location the experimentally measured fiber temperature is 94°C.The reason the model predictions are high compared to experimental data may be due to the transverse vibratory motion of the fiber.It has been shown that,for heat transfer from wires to air in parallel flow,vibrations cause an increase in the heat-transfer co-efficient.24,38Even in melt spinning,Marla et al.28reported that the fiber exhibits transverse vibrations.However,these vibra-tions were quite low in amplitude compared to the vibrations observed in melt blowing.In an effort to more closely predictfiber temperature,future modeling of melt blowing may include the effects of vibration.The centerline air temperature is also shown in Figure 19.This temperature was measured with a thermocouple in the absence of polymer flow.It can be seen that the fiber remains warmer than the air for the first 25cm below the die.Figures 20and 21are similar to Figure 19,except that Figures 20and 21are for,respectively,polymer rates of 1.00and 1.25g/min.Increasing the polymer rate causes the predicted tem-perature profiles (with either correlation)to be higher.Figure 22,when compared with Figure 19,shows the effect of changing the air flow rate from 100to 125slpm.With either correlation,the predicted profiles are only slightly (a few degrees)lower for the 100slpm flow rate than for the 125slpm flow rate.ConclusionsOnline measurements of fiber diameter and temperature were made in the melt-spinning and melt-blowing processes.In melt spinning,the fibers seemed to cool slightly faster with an increase in spinning speed for both PP and PB.When theheat-Figure 15.Fiber diameter as a function of position below the die for the melt blowing of polypropylene.Data are shown for three different polymerrates.Figure 16.Model predictions of fiber diameter for the melt blowing of polypropylene.Two different correlations for Nusselt number were tested in the model.Experimental data are also included forcomparison.Figure 17.Fiber diameter as a function of position below the die for the melt blowing of polypropylene.Data are shown for two different air flowrates.Figure 18.Fiber temperature as a function of position below the die for the melt blowing of polypropylene.Data are shown for three different polymer rates.8742Ind.Eng.Chem.Res.,Vol.48,No.18,2009。