Argentum Solutions, Inc.

    Sterling guidance on corrosion and materials degradation


 

Potential-pH Diagrams
THERMEXPERT - Potential-pH diagram generator

Intelligent Tools

POLEXPERT - Polarization Scan Artificial Neural Network Expert System

SEQEXPERT - Sequential Immersion Test Artificial Neural Network Expert System

CYLEXPERT - Rotating Cylinder Electrode Intelligent Rotation Rate Calculator

Corrosion Calculator

Corrosion Rate Calculator

Corrosion Economics Estimator

FINCALCULATOR - Corrosion Economic Calculator


TUTORIAL ON CYLEXPERTTM AND THE ROTATING CYLINDER ELECTRODE

David C. Silverman


Table of Contents

Introduction-What is CYLEXPERT?
Overview of the Experimental Technique
Using CYLEXPERT-a step-by-step procedure
Hydrodynamics of a Smooth Cylinder
  1. Boundary Layer
  2. Mass Transfer Correlations
  3. Assessing Mass Transfer Control
         Models Relating Geometries
Effect of Surface Roughness
Wall Shear Stress vs. Mass Transfer
CYLEXPERT (The Intelligent Experimental Design Tool)
List of Symbols



Hydrodynamics of a Smooth Cylinder

1. Boundary Layer


Understanding the physics of the rotating cylinder electrode requires understanding the concept of the boundary layer. When a fluid moves past a stationary, hydraulically smooth  surface, there is an implicit assumption that the fluid is stationary at that surface. That assumption is often called "no slip at the wall". In many flow configurations and especially when turbulent flow conditions prevail, the velocity profile shows an increase from 0 at the solid-fluid interface to the free stream value across a relatively short distance into the fluid. Though in the rotating cylinder the surface is moving, the fluid adjacent to that surface moves at the same velocity as the cylinder so, in effect, the fluid is stationary relative to the surface and fulfills the no slip at the wall criterion. The boundary layer thickness is independent of position on the surface meaning that the velocity profile, momentum transfer, and wall shear stress are independent of position. The rotating cylinder electrode as traditionally constructed cannot examine situations where such uniformity does not exist without some modification.

Mass transfer requires a concentration gradient between the surface and fluid bulk. Such a gradient implies that the concentration changes across some small distance between the surface and the bulk. This region is called the concentration or mass transfer boundary layer. For the large Schmidt Numbers normally encountered in liquids the fully developed mass transfer boundary layer for hydraulically smooth surfaces is much thinner than the fully developed hydrodynamic boundary layer. This relationship is shown in this figure . The hydrodynamic boundary layer is represented as having two sections, a viscous sublayer and a turbulent outer layer and that the fluid is moving relative to a surface. Since the flow is turbulent, the profiles are time-averaged profiles and are not drawn to actual scale. As in the case of the velocity profile, the mass transfer boundary layer thickness is independent of position meaning that the mass transfer rate is independent of position. A full discussion of the meaning of the boundary layer can be found in the literature (H. Schlichting, Boundary-Layer Theory, 7th edition, McGraw-Hill Book Co., New York, 1979).

Top



2. Mass Transfer Correlations


Corrosion being limited by the rate of mass transfer has usually been assumed when the slope of the line calculated from the measured mass transfer rate versus velocity, both in logarithmic format, is approximately 0.7 (Eisenberg, M., Tobias, C. W., and Wilke, C. R., J. Electrochem. Society, 101(6), 306 (1954)). This exponent was derived using the Chilton-Colburn analogy (Chilton, T. H. and Colburn, A. P., Ind. Eng. Chem., 26, 1183 (1934)). The curvature of the friction factor vs. Reynolds number relationship was fit over the range of Reynolods number of 103 to 105 in linear format in log-log coordinates by:

                                                                        (1)


This equation led to the commonly used expression relating the Sherwood number to the Reynolds number of:

                                                                        (2)


But, this relationship is not the only relationship that has been developed for the rotating cylinder electrode. Equation (1) is actually a linearization of a curved function. What may be thought to be a deviation from the 0.7 power as caused by a departure from complete mass transfer control may be nothing more than deviations caused by the limitations of the predictive power of the "linearized" equation. Understanding the predictive limitations is important for making appropriate conclusions about mass transfer control.

Early measurements of drag coefficients on rotating cylinders over the range of Reynolds numbers between about 103 and 106 resulted in the following relationship between the friction factor and Reynolds number:

                                                                        (3)


(Theodorsen, T. and Regier, A., "Experiments on Drag of Revolving Disks, Cylinders,and Streamline Rods at High Speeds" Nat. Advisory Comm. Aeronaut., p. 367, Report No. 793, U.S. Government Printing Office, Washington, D.C., 1945.) This equation itself was derived as a best fit to friction factor data plotted as vs. . The relationship between friction factor and Reynolds number implied by this equation and the data are actually curved on log-log coordinates. From the Chilton-Colburn analogy the Sherwood number or mass transfer coefficient can be related to the friction factor by:

                                                                        (4)


where the exponent c-1 has the value of -0.644 to -0.69 depending on author and study and St is the Stanton number. Equation (4) is substituted into equation (3) to obtain the Sherwood number versus Reynolds number relationship from the friction factor. The resulting equation would require an iterative calculation of the Sherwood number because the Sherwood number (mass transfer coefficient) would appear on both sides of the equation. Equation (2) is one approximation of this relationship.

The desire of a linear approximation has led to a number of Sherwood number vs. Reynolds number relationships being developed for the hydraulically smooth rotating cylinder. The following table is a compilation of a number of the relationships. CYLEXPERT uses this type of relationship for the rotating cylinder electrode portion of the calculation.

Experimental Conditions
Ferricyanide-ferrocyanide on nickel electrode
103≤Re≤105
Ferricyanide-ferrocyanide on platinum electrode
103≤Re≤2x104
Theoretical calculation compared to data from
ferricyanide-ferrocyanide on nickel electrode
3x104≤Re≤1.3x106
Theoretical calculation with constant assumed
to be 0.0791
Cathodic deposition of copper from copper sulfate
104≤Re≤5x105
Copper-copper sulfate deposition on copper electrode
103≤Re≤1.9x104
Compilation of 290 mass transfer data points
102≤Re≤105
Copper-copper sulfate deposition on copper electrode
103≤Re≤1.9x104
Oxygen reduction in sodium chloride on Monel
and sodium sulfate on steel
2x103≤Re≤1.2x105
Derived from curve-fit of Theodorsen and Regier over range
2x102≤Re≤4x106


This figure. shows the various relationships in (Sherwood number)/(Schmidt number)c versus Reynolds number format. The exponent on the Schmidt number has been assumed to be constant across the relationships even though some variation in that exponent exists in the correlations. Changing the exponent on the Schmidt number would change the coefficient slightly but would not affect the conclusions (see the next to the last correlation in the table). The heavy dashed line in the figure was derived (Silverman, D. C., Corrosion, 59(3), 207 (2003)1    (328k)) by iteratively solving equation (3) followed by using equation (4) to estimate the quantity Sh/Scc. The area denoting mass transfer control in the diagram found in Method B for assessing mass transfer control in CYLEXPERT is the area surrounding all of these curves. As discussed in Assessing Mass Transfer Control or in Using CYLEXPERT-a step by step procedure the actual points measured in any system should be placed on a plot such as this one to more fully assess the degree of mass transfer control.

Though the Sherwood number vs. Reynolds number relationships are clustered fairly closely together, the relative difference between the largest and smallest estimates at a Reynolds number of 100 is about a factor of 3 and the difference at a Reynolds number of 106 is about a factor of 2. While equation (3) does provide a reasonable representation of the relationship and a reasonable estimate of the mass transfer coefficient, the estimate is just that, an estimate. As shown in the figure some of the equations shown in the table (Silverman, D. C., Corrosion, 59(3), 207 (2003)1    (328k)) could actually provide a slightly better prediction at higher Reynolds numbers but the difference still lies within the extremes. The implication of this discussion is that if equation (3) or any of the equations in the above table are used and measured data points do not lay on it but are still within the extremes of the variability, then mass transfer control is likely a dominant contributor to the mechanism.

Top



3. Assessing Mass Transfer Control

To determine if mass transfer controls the corrosion rate, one might measure the actual corrosion rate as function of rotation rate either electrochemically or by mass loss. Then, one of two methods could be used to characterize the mass transfer influence on the corrosion process. In all cases, if the data points or slope are similar to those expected from the literature-derived relationship for the rotating cylinder, mass transfer probably controls the corrosion rate. These methods are presented in CYLEXPERT and are also found in the section on Using CYLEXPERT-a step by step procedure.

Method #1

This method is more qualitative and may be less accurate.

  1. Estimate the "corrosion rate" at several rotation rates, three rotation rates being the minimum. Mass loss, corrosion current, or polarization resistance could be used to estimate the corrosion rate. The hydrodynamics must be turbulent at all times. That corrosion is at steady state is implicitly assumed. Generating this information over several days can ensure that corrosion is at steady state.

  2. Plot the logarithm of the corrosion rate estimate versus the logarithm of the rotation rate as measured in rpm or in cm/s. This relationship should be fairly linear when plotted in this format.

  3. Perform a linear regression of the data points. Place the regression line on the same plot as the data points. Examine the slope and the agreement of the regression with the measured data. Significant curvature in the measured data relative to the regression may suggest additional processes are occurring or the measurement is in error. If agreement between the regression and the measured values is poor ensure that the experimental set-up and measurement technique and analysis are correct. If they are, other processes not measured by the rotating cylinder electrode could be occurring. Proceed on if measured and calculated values in the regression are in reasonable agreement.

    • If the slope of the regression lies between about 0.65 and 0.75, the corrosion rate is likely to be significantly influenced or even controlled by mass transfer.
    • If the slope of the regression is much greater than 0.75, artifacts such as surface roughness or particulate erosion could be affecting the results. Examine the electrode under magnification for surface damage or surface roughening.
    • If the slope is much less than about 0.65 but is still greater than zero, mass transfer might be having some influence but other factors may be detracting from fluid motion being the only influence on the corrosion rate.
    • If the slope is close to zero, mass transfer is probably not having an effect.


Method #2

This method is more quantitative and could be more accurate but requires more knowledge of the physical properties of the environment.

  1. Estimate the "corrosion rate" at several rotation rates, three rotation rates being the minimum. Mass loss, corrosion current, or polarization resistance could be used for the measurement. For this calculation, the rate must be converted to moles/s to be consistent with the units in the equations below. The fluid motion must be turbulent at all rotation rates and the equipment must be able to handle the rotation rate. Generating this information over several days can ensure that corrosion is at steady state.

  2. Estimate the mass transfer coefficient from the corrosion rate. This estimate requires use of the following equation:

    where
    kcyl is the mass transfer coefficient, cm/s.
    |Cbulk-Csurf| is absolute value of the difference in concentration of the rate limiting species between its concentration in the bulk fluid and its concentration at the surface. One of the values is usually zero when the reaction is controlled by mass transfer. The value that is not zero usually has to be measured. These concentrations are usually expressed as mole/cm3.
    Acyl is the active area of the rotating cylinder, cm2.
    Corrosion Rate is the corrosion rate in the same units as the concentration, moles/cm3. The corrosion current is converted to these units by using the Faraday constant and the number of electrons transferred for each molecule that reacts. If the rate is in mass loss per unit time, the equivalent weight of the alloy is used for the conversion to moles per unit time.

  3. Calculate the Sherwood Number, Shcyl by

    The Sherwood number is dimensionless. Ensure that the units for each parameter are consistent and that the Sherwood number as calculated is dimensionless.

  4. Calculate the Schmidt Number, Sc, by

    The diffusion coefficient can sometimes be found in physical property tables. Otherwise, a reasonable value to assume is often in the range of 10-5 cm2/s. Ensure that the units for each parameter are consistent and that the Schmidt number as calculated is dimensionless.

  5. Calculate the Reynolds Number by

    Ensure that the units for each parameter are consistent and that the Reynolds number as calculated is dimensionless.

  6. Divide the Sherwood number by the Schmidt number raised to the appropriate power for the rotating cylinder electrode correlation you are using. The most widely accepted value of that power is 0.356. Plot the logarithm of that calculation versus the logarithm of the Reynolds number.

  7. For complete mass transfer control on a hydraulically smooth cylinder, the results should fall within the area labeled as such on the following plot:

    Click for larger image

    If the results fall well below the region but still have a slope, then there is probably some mass transfer influence but additional resistances to corrosion probably exist. If the results are significantly above the region, then surface roughness may be playing a role. Examples of this approach can be found in <Silverman, D. C., Corrosion, 40(5), 220 (1984)1    (467k) and Silverman, D. C. and Zerr, M. E., Corrosion, 42(11), 633(1986) 1    (456k))






Previous Page: Using CYLEXPERT-a step-by-step procedure

Next Page: Models Relating Geometries

Return to Table of Contents



1 © NACE International publication and year shown in citation above. All rights reserved. Displayed with permission from NACE International, Houston, TX (http://www.nace.org). Published in Corrosion, in the month and year shown in the citation above.






David C. Silverman, Ph.D. - Primary Consultant
E-Mail:     dcsilverman@argentumsolutions.com
Phone:     314-576-3586
Fax:         314-754-9825
Address:   The Argentum House
                14314 Strawbridge Ct.
                Chesterfield, MO 63017