Flow-through and flow-by porous electrodes of nickel foam. II. Diffusion-convective mass transfer between the electrolyte and the foam

The work described here concerns the diffusion-convective mass transfer to flow-through and flow-by porous electrodes of nickel foam. Empirical correlations giving the product of the mass transfer coefficient and the specific surface areaa _e of the material as a function of the pressure drop per unit electrode height and as a function of the grade characterizing the foam are proposed. The performance of various materials are compared in terms of vs the mean linear electrolyte flow velocity.Nomenclature a _e specific surface area (per unit of total volume of electrode) (m^–1) - A, B Ergun law coefficients determined in flow-by configuration - A, B Ergun law coefficients determined in flow-through configurationA, A (Pa m^–3 s²);B, B (Pa m² s^–1) - C _E entering concentration of ferricyanide ions (mole m^–3) - D molecular diffusion coefficient (m² s^–1) - F Faraday number (C mol^–1) - G grade of the foams - I _L limiting current (A) - mean mass transfer coefficient (m s^–1) - n number of stacked foam sheets in the electrode - P/H pressure drop per unit of height (Pa m^–1) - Q _v volumetric electrolyte flow rate (m³ s^–1) - Re Reynolds number - Sc Schmidt number - Sh Sherwood number - T mean tortuosity of the foam pores - mean electrolyte velocity (m s^–1) - V _R electrode volume (m³) - X conversion - dynamic viscosity (kg m^–1 s^–1) - v number of electrons in the electrochemical reaction - v kinematic viscosity (m² s^–1)     

Determination of electrode kinetics by corrosion potential measurements: Zinc corrosion by bromine

An attractive way of determining the electrode kinetics of very fast dissolution reactions is that of measuring the corrosion potential in flowing solutions. This study analyses a critical aspect of the corrosion potential method, i.e., the effect of nonuniform corrosion distribution, which is very common in flow systems. The analysis is then applied to experimental data for zinc dissolution by dissolved bromine, obtained at a rotating hemispherical electrode (RHE). It is shown that in this case the current distribution effect is minor. However, the results also indicate that the kinetics of this corrosion system are not of the classical Butler-Volmer type. This is explained by the presence of a chemical reaction path in parallel with the electrochemical path. This unconventional corrosion mechanism is verified by a set of experiments in which zones of zinc deposition and dissolution at a RHE are identified in quantitative agreement with model predictions. The practical implications for the design of zinc/bromine batteries are discussed.Notation C _i concentration of species i (mol cm^–3) - D _` diffusivity of species i (cm² s^–1) - F Faraday constant - i _j current density of species j (A cm^–2) - i ₀ ^b exchange current density referenced at bulk concentration (A cm^–2) - J , inverseWa number - N - n number of electrons transferred for every dissolved metal atom - P _m Legendre polynomial of orderm - r ₀ radius of dise, sphere, or hemisphere - s stoichiometric constant - t ₊ transference number of metal ion - V _corr corrosion overpotential (V) Greek letters anodic transfer coefficient of Reaction 21b - _a anodic transfer coefficient of metal dissolution - _c cathodic transfer coefficient of metal dissolution - anodic transfer coefficient of zinc dissolution - velocity derivative at the electrode surface - (x) incomplete Gamma function - , exchange reaction order ofM ⁺ⁿ - , inverseWa number - _a activation overpotential (V) - _c concentration overpotential (V) - polar angle (measured from the pole) (rad) - k solution conductivity (^–1 cm^–1) - kinematic viscosity (cm² s^–1) - ₀ solution potential at the electrode surface (V) - rotation rate (s^–1) - * indicates dimensionless quantities     

Mass transport in a weakly stratified electrochemical cell

A study of natural convection in an electrochemical system with a Rayleigh number of the order 10¹⁰ is presented. Theoretical and experimental results for the unsteady behaviour of the concentration and velocity fields during electrolysis of an aqueous solution of a metal salt are given. The cell geometry is a vertical slot and the reaction kinetics is governed by a Butler-Volmer law. To reduce the effects of stratification, the flush mounted electrodes are located (symmetrically) in the middle parts of the vertical walls. It is demonstrated, both theoretically and experimentally, that a weak stratification develops after a short time, regardless of cell geometry, even in the central part of the cell. This stratification has a strong effect on the velocity field, which rapidly attains boundary layer character. Measured profiles of concentration and vertical velocity at and above the cathode are in good agreement with numerical predictions. For a constant cell voltage, numerical computations show that between the initial transient and the time when stronger stratification reaches the electrode area, the distribution of electric current is approximately steady.List of symbols a _i left hand side of equation system - b _i right hand side of equation system - c concentration (mol m^–3) - c dimensionless concentration - c _i concentration of species i' (mol m^–3) - c₀ initial cell concentration (300 mol m^–3) - c ₀ dimensionless initial cell concentration - c_wall concentration at electrode surface (mol m^–3) - dx increment solution vector in Newton's method - D _i diffusion coefficient of species i (m² s^–1) - D ₁ 0.38 × 10^–9 m² s^–1 - D ₂ 0.82 × 10^–9 m² s^–1 - D effective diffusion coefficient of the electrolyte (0.52 × 10^–9 m² s^–1) - _x unit vector in the vertical direction - _y unit vector in the horizontal direction - F Faraday's constant (96 487 A s mol^–1) - g acceleration of gravity (9.81 m s^–2) - i dummy referring to positive (i = 1) or negative (i = 2) ion - f current density (A m^–2) - f dimensionless current density - i₀ exchange current density (0.01 A m^–2) - J _ij Jacobian of system matrix - L length of electrode (0.03 m) - N _i transport flux density of ion i (mol m^–2 s^–1) - n unit normal vector - p pressure (Nm^–2) - p dimensionless pressure - R gas constant molar (8.31 J K^–1 mol^–1) - R _i residual of equation system - Ra Rayleigh number gL ³ c ₀/D (2.54 × 101¹⁰) - S _c Schmidt number /D (1730) - t time (s) - t dimensionless time - T temperature (293 K) - velocity vector (m s^–1) - dimensionless velocity vector - U characteristic velocity in the vertical direction - V _± potential of anode and cathode, respectively - x spatial coordinate in vertical direction (m) - x dimensionless spatial coordinate in vertical direction - x solution vector for c, and - y spatial coordinate in horizontal direction (m) - y dimensionless spatial coordinate in horizontal direction - z _i charge number of ion i Greek symbols symmetry factor of the electrode kinetics, 0.5 - volume expansion coefficient (1.24 × 10^–4 m³ mol^–1) - _s surface overpotential - constant in equation for the electric potential (–5.46) - _s diffusion layer thickness - scale of diffusion layer thickness - constant relating c/y to the Butler-Volmer law (0.00733) - kinematic viscosity (0.9 × 10^–6 m2 s^–1)     

Behaviour of a tall vertical gas-evolving cell Part II: Distribution of current

Vertical electrolysers with a narrow electrode gap are used to produce gases, for example, chlorine, hydrogen and oxygen. The gas voidage in the solution increases with increasing height in the electrolyser and consequently the current density is expected to decrease with increasing height. Current distribution experiments were carried out in an undivided cell with two electrodes each consisting of 20 equal segments or with a segmented electrode and a one-plate electrode. It was found that for a bubbly flow the current density decreases linearly with increasing height in the cell. The current distribution factor increases with increasing average current density, decreasing volumetric flow rate of liquid and decreasing distance between the anode and the cathode. Moreover, it is concluded that the change in the electrode surface area remaining free of bubbles with increasing height has practically no effect on the current distribution factor.Notation A _e electrode surface area (m²) - A _e,s surface area of an electrode segment (m²) - A _{e, 1–19} total electrode surface area for the segments from 1 to 19 inclusive (m²) - A _e,a anode surface area (m²) - A _e,a,h A _e,a remaining free of bubbles (m²) - A _e,e cathode surface area (m²) - A _e,c,h A _e,c remaining free of bubbles (m²) - a ₁ parameter in Equation 7 (A^–1) - B current distribution factor - B _r B in reverse position of the cell - B _s B in standard position of cell - b _a Tafel slope for the anodic reaction (V) - b _c Tafel slope for the cathodic reaction (V) - d distance (m) - d _ac distance between the anode and the cathode (m) - d _wm distance between the working electrode and an imaginary membrane (m) (d _wm=0.5d _wt=0.5d _ac) - d _wt distance between the working and the counter electrode (m) - F Faraday constant (C mol^–1) - h height from the leading edge of the working electrode corresponding to height in the cell (m) - h _e distance from the bottom to the top of the working electrode (m) - I current (A) - I _s current for a segment (A) - I ₂₀ current for segment pair 20 (A) - I _1–19 total current for the segment pairs from 1 to 19 inclusive (A) - i current density (A m^–2) - i _av average current density of working electrode (A m^–2) - i _b current density at the bottom edge of the working electrode (A m^–2) - i ₀ exchange current density (A m^–2) - i _0,a i ₀ for anode reaction (A m^–2) - i _l current density at the top edge of the working electrode (A m^–2) - n ₁ parameter in Equation 15 - n _s number of a pair of segments of the segmented electrodes from their leading edges - Q _g volumetric rate of gas saturated with water vapour (m³ s^–1) - Q ₁ volumetric rate of liquid (m³ s^–1) - R resistance of solution () - R ₂₀ resistance of solution between the top segments of the working and the counter electrode () - R _p resistance of bubble-free solution () - R _p,20 R _p for segment pair 20 () - r _s reduced specific surface resistivity - r _s,0 r _s ath=0 - r _s,20 r _s for segment pair 20 - r _s, r _s for uniform distribution of bubbles between both the segments of a pair - r _s,,20 r _s, for segment pair 20 - T temperature (K) - U cell voltage (V) - U _r reversible cell voltage (V) - v ₁ linear velocity of liquid (m s^–1) - v _1,0 v ₁ through interelectrode gap at the leading edges of both electrodes (m s^–1) - x distance from the electrode surface (m) - gas volumetric flow ratio - ₂₀ at segment pair 20 - specific surface resistivity ( m²) - _t at top of electrode ( m²) - _p for bubble-free solution ( m²) - _b at bottom of electrode ( m²) - thickness of Nernst bubble layer (m) - ₀ ath=0 (m) - _{0,i 0} ati - voidage - _x,0 atx andh=0 - _0,0 voidage at the leading edge of electrode wherex=0 andh=0 - _0,0 ati _b - _0,0 ati=i _t - _,h voidage in bulk of solution at heighth - _,20 voidage in bubble of solution at the leading edge of segment pair 20 - _lim maximum value of _0,0 - overpotential (V) - _a anodic overpotential (V) - _c cathodic overpotential (V) - _h hyper overpotential (V) - _h,a anodic hyper overpotential (V) - _h,c cathodic hyper overpotential (V) - fraction of electrode surface area covered by of bubbles - _a for anode - _c for cathode - resistivity of solution ( m) - _p resistivity of bubble-free solution ( m)     

Modelling of time-dependent performance criteria in a three-dimensional cell system during batch recirculation copper recovery

A three-dimensional electrode cell with cross-flow of current and electrolyte is modelled for galvanostatic and pseudopotentiostatic operation. The model is based on the electrodeposition of copper from acidified copper sulphate solution onto copper particles, with an initial concentration ensuring a diffusion-controlled process and operating in a batch recycle mode. Plug flow through the cell and perfect mixing of the electrolyte in the reservoir are assumed. Based on the model, the behaviour of reacting ion concentration, current efficiency, cell voltage, specific energy consumption and process time on selected independent variables is analysed for both galvanostatic and pseudopotentiostatic modes of operation. From the results presented it is possible to identify the optimal values of parameters for copper electrowinning.List of symbols a specific surface area (m^–1) - A cross-sectional area (mu²) - a _a Tafel constant for anode overpotential (V) - a _II Tofel constant for hydrogen evolution overpotential (V) - b _a Tafel coefficient for anode overpotential (V decade^–1) - b _H Tafel coefficient for hydrogen evolution overpotential (V decade^–1) - C _e concentration at the electrode surface (m) - C _L cell outlet concentration (m) - C ₀ cell inlet concentration (m) - C ₀ ⁰ initial cell inlet concentration att = 0 (m) - d _p particle diameter (m) - e, e _p current efficiency and pump efficiency, respectively - E specific energy consumption (Wh mol^–1) - E solution phase potential drop through the cathode (V) - F Faraday number (C mol^–1) - h interelectrode distance (m) - i, i _L current density and limiting current density, respectively (A m^–2) - I, I _L current and limiting current, respectively (A) - I _H partial current for hydrogen evolution (A) - k _L mass transfer coefficient (m s^–1) - L bed height (m) - l bed depth (m) - M molecular weight (g mol^–1) - N power per unit of electrode area (W m^–2) - n exponent in Equation 19 - P pressure drop in the cell (N m^–2) - Q electrolyte flow rate (m³ h^–1) - R Universal gas constant (J mol^–1 K^–1) - r _e electrochemical reaction rate (mol m^–2 h^–1) - t _c critical time for operating current to reach instantaneous limiting current (s) - t _p process time to reach specified degree of conversion (s) - T temperature (K) - u electrolyte velocity (m s^–1) - U total cell voltage (V) - U ₀ reversible decomposition potential (V) - U _ohm ohmic voltage drop between anode and threedimensional cathode (V) - V volume of electrolyte (m³) - z number of transferred electrons Greek letters ratio of the operating and limiting currents - _{A, a} anodic activation overpotential (V) - _{c, e} cathodic concentration overpotential (V) - bed voidage - _H void fraction of hydrogen bubbles in cathode - constant (Equation 2) - ₀ electrolyte conductivity (ohm^–1 m^–1) - v electrolyte kinematic viscosity (m² s^–1) - _d diaphragm voltage drop (V) - _H voltage drop due to hydrogen bubble containing electrolyte in cathode (V) - electrolyte density (kg m^–3) - _p particle density (kg M^–3) - reservoir residence time (s)     

Generalized analytical expressions for Tafel slope,reaction order and a.c. impedance for the hydrogen evolution reaction (HER): mechanism of HER on platinum in alkaline media

Following the generally accepted mechanism of the HER involving the initial proton discharge step to form the adsorbed hydrogen intermediate, which is desorbed either chemically or electrochemically, generalized expressions for the Tafel slope, reaction order and the a.c. impedance for the hydrogen evolution reaction are derived using the steady-state approach, taking into account the forward and backward rates of the three constituent paths and the lateral interactions between the chemisorbed intermediates. Limiting relationships for the Tafel slope and the reaction order, previously published, are deduced from these general equations as special cases. These relationships, used to decipher the mechanistic aspects by examining the kinetic data for the HER on platinum in alkaline media, showed that the experimental observations can be consistently rationalized by the discharge-electrochemical desorption mechanism, the rate of the discharge step being retarded on inactive platinum compared to the same on active platinum.Nomenclature C _d double-layer capacity (µF cm^–2) - E _rev reversible electrode potential (V) - F Faraday number (96 487 C mol^–1 ) - R gas constant - T temperature (K) - Y _f Faradaic admittance (^–1 cm^–2) - Y _t Total admittance (^–1 cm^–2) - Z _f Faradaic impedance ( cm²) - i _f total current density (A cm^–2) - i _nf nonfaradaic current density (A cm^–2) - j - k ₀ ¹ rate constant of the steps described in Equations 1 to 3 (mol cm^–2 s^–1 ) - j - qmax saturation charge (µC cm^–2) - Laplace transformed expressions for i, and E - _{1 3} symmetry factors for the Equations 1 and 3 - saturation value of adsorbed intermediates (mol cm^–2) - overpotential - coverage by adsorbed intermediates - angular frequency This paper is dedicated to Professor Brian E. Conway on the occasion of his 65th birthday, and in recognition of his outstanding contribution to electrochemistry.     

Behaviour of a tall vertical gas-evolving cell. Part I: Distribution of void fraction and of ohmic resistance

Electrolysis of a 22 wt % NaOH solution has been carried out in a vertical tall rectangular cell with two segmented electrodes. The ohmic resistance of the solution between a segment pair has been determined as a function of a number of parameters, such as, current density and volumetric rate of liquid flow. It has been found that the ohmic resistance of the solution during the electrolysis increases almost linearly with increasing height in the cell. Moreover, a relation has been presented describing the voidage in the solution as a function of the distance from the electrodes and the height in the cell.Notation A _e electrode surface area (m²) - a _s parameter in Equation 12 (A^–1) - b _s parameter in Equation 12 - d distance (m) - d _ac distance between the anode and the cathode (m) - d _wm distance between the working electrode and an imaginary separator (m) - F Faraday constant (C mol^–1) - h height from the leading edge of the working electrode corresponding to height in the cell (m) - h _e distance from the bottom to the top of the working electrode (m) - h _s height of a segment of working electrode (m) - I current (A) - I ₂₀ current for segment pair 20 (A) - I _1–19 total current for the segment pairs from 1 to 19 inclusive (A) - I _x-19 total current for the segment pairs fromx to 19 inclusive (A) - i current density A m^–2 - N _s total number of gas-evolving pairs - n ₁ constant parameter in Equation 8 - n _a number of electrons involved in the anodic reaction - n _c number of electrons involved in the cathodic reaction - n _s number of a pair of segments of the segmented electrodes from their leading edges - Q _g volumetric rate of gas saturated with water vapour (m³ s^–1) - Q ₁ volumetric rate of liquid (m³ s^–1) - R resistance of solution () - R ₂₀ resistance of solution between the top segments of the working and the counter electrode () - R _p resistance of bubble-free solution () - R _p,20 R _p for segment pair 20 () - r _s reduced specific surface resistivity - r _s,0 r _s ath=0 - r _s,20 r _s for segment pair 20 - r _s, r _s for uniform distribution of bubbles between both the segments of a pair - r _s,,20 r _s, for segment pair 20 - S _b bubble-slip ratio - S _b,20 S _b at segment pair 20 - S _b,h S _b at heighh in the cell - T temperature (K) - V _m volume of 1 mol gas saturated with water vapor (m³ mol^–1) - v ₁ linear velocity of liquid (m s^–1) - v _1,0 v ₁ through interelectrode gap at the leading edges of both electrodes (m s^–1) - W _e width of electrode (m) - X distance from the electrode surface (m) - Z impedance () - Z real part of impedance () - Z imaginary part of impedance () - resistivity of solution ( m) - _p resistivity of bubble-free solution ( m) - gas volumetric flow ratio - ₂₀ at segment pair 20 - _s specific surface resistivity ( m²) - _{s, p s} for bubble-free solution ( m²) - thickness of Nernst bubble layer (m) - ₀ ath=0 (m) - voidage - _x,0 atx andh=0 - _0,0 voidage at the leading edge of electrode wherex=0 andh=0 - _,h voidage in bulk of solution at heighth - ₂₀ voidage in bubble of solution at the leading edge of segment pair 20     

Fundamental studies on a new concept of flue gas desulphurization

A new process for removal of sulphur dioxide from waste gases is proposed consisting of both electrochemical and catalytic sulphur dioxide oxidation. In the catalytic step a part of the sulphur dioxide is oxidized by oxygen on copper producing sulphuric acid and copper sulphate. The other part is oxidized electrochemically on graphite. The cathodic reaction of this electrolysis is used for recovering the copper dissolved in the catalytic step. The basic reactions of this process have been studied experimentally in detail. It has been shown that sulphur dioxide can be electrochemically oxidized on carbon electrodes to sulphuric acid with high current efficiency. The reaction rate of the electrochemical copper deposition is increased by dissolved sulphur dioxide in the electrolyte. The catalytic oxidation of sulphur dioxide on copper has been investigated for different sulphur dioxide concentrations and temperatures. The ratio of the reaction products, sulphuric acid and copper sulphate, varies over a wide range depending on the experimental conditions.Nomenclature SO₂ concentration (gas phase) (vol % SO₂) - SO₂ concentration (electrolyte) (g l^–1) - E potential vs saturated calomel electrode (V) - E _s specific energy consumption (W g^–1 SO₂) - F Faraday constant (A s^–1 mol^–1) - i current density (mA cm^–2) - molecular weight (g mol^–1) - T temperature (° C) - U _c cell voltage (V) - v _e number of electrons being transferred - space-time yield of SO₂-oxidation (g SO₂ h^–1 dm^–3) - _cu space-time yield of Cu-corrosion (g Cu h^–1 dm^–3) - ratio - fractional conversion of SO₂ - current efficiency for SO₂ oxidation     

Joule heating induced local electrodeposition for microelectronic circuits

A fundamental study is performed for local electrodeposition of copper utilizing thermal potential induced by Joule heating. The feasibility of the process for microelectronic applications is assessed by both experiment and mathematical modeling. The results of the investigation show that (i) a copper wire is coated under conditions of a.c. 50 Hz Joule heating in electrolyte containing 1.0 M CuSO₄ and 0.5m H₂SO₄ with relatively high deposition rate of about 0.4 µm min^–1, (ii) the Joule heating current should be kept below the boiling point of the solution to realize uniform deposition, and (iii) results of calculations by the present model based on one-dimensional heat conduction agree well with experimental results.Nomenclature D diameter of wire (m) - D ₀ initial diameter of wire (m) - F Faraday constant (96 487 C mol¹ ) - g acceleration due to gravity (9.807 m s²) - Gr Grashof number - H thickness of electrodeposit (m) - I current (A) - i ₀ exchange current density (Am^–2) - i _n current density normal to electode (Am^–2) - J current density (I/S) (Am^–2) - L length of wire (m) - M molar concentration of electrolyte (mol dm^–3 or M) - m atomic weight (kg mol^–1) - n number of electrons participating - n unit normal vector to boundary - Nu Nusselt number - Pr Prandtl number - q heat per unit volume (W m^–3) - R universal gas constant (8.314 3 J mol^–1 K^–1) - (r, z) cylindrical coordinate (m) - S cross section of wire (m²) - T temperature (K) - T ₀ fixed temperature at both ends of wire (K) - T _y temperature of electrolyte (K) - t time (s) - x longitudinal coordinate over wire (m) Greek symbols heat transfer coefficient (W m^–2 K^–1 - _a,c anodic (a) and cathodic (c) transfer coefficient - thermal expansion coefficient of solution (K^–1) - specific heat (J kg^–1K^–1) - potential (V) - _e electrode potential (V) - thermal conductivity (W m^–1 K^–1 ) - _y ionic conductivity of electrolyte (^–1m^–1) - _e electronic conductivity of electrode (^–1 m^–1) - kinematic viscosity (m²s^–1) - surface overpotential ( _e – ) (V) - time constant (s) - density (kg m^–3) This work was presented at The 7th International Microelectronics Conference, Yokohama, Japan (1992).     

Enhanced mass transfer in electrochemical cells using turbulence promoters

Many electrochemical processes suffer in varying degrees from mass transfer limitations. These limitations may require operation at considerably less than economic optimum current densities. Mass transfer to a surface may be considerably enhanced by insertion of turbulence promoters in the fluid flow path near the affected surface.An instrument was developed to measure local current densities in the hydrodynamically very difficult region near the turbulence promoter. A general method for the relative evaluation of hydrodynamic conditions has been developed. Generalization of the data permits optimization of hydrodynamic cell design using the promoter shapes investigated.

Notation

Symbols A Coefficient for cell power costs, $ m² (As)^–1 - A _c Cell area, m² - a Constant in Equation 4 - B Coefficient for area-proportional costs, $ A (m² s)^–1 - C Coefficient for pumping power costs, $ A (m² s)^–1 - C _b Bulk concentration, kg mol m^–3 - C _bi Inlet bulk concentration, kg mol m^–3 - C _e Energy cost, $ (Ws)^–1 - C _i Interfacial concentration, kg mol m^–3 - ¯C _s Amortized area cost, $ (m² s)^–1 - D Current—density-insensitive costs, $ s^–1 - D _e Equivalent diameter, m - D Diffusion constant, m² s^–1 - e Current efficiency - F _d Cell feed rate, m³ s^–1 - F 96.5×10⁶ A s kg eq^–1 - g Channel width, m - h Channel height, m - i Current density, A m^–2 - i _opt Economic optimum current density, A m^–2 - K Total costs of running cell, $ s^–1 - (K–D)_ideal Total sensitive costs under hydrodynamically ideal conditions, $ s^–1 - k _c Convective mass transfer coefficient, m s^–1 - L Total length of flow path, m - l Promoter spacing, m - N Mass flow rate to surface due to convection, kg mol m² s^–1 - n _e Number of electrons transferred in electrode reaction - P _c Power required by cell, W - P/L Average pressure gradient in channel, N m^–3 - R _av Effective cell resistance, m² - S Open channel cross-section, m² - S ₀ Minimum channel cross-section at promoter, m² - s _i Stoichiometric coefficient of species i - t _i Transport number of species i in solution - ¯t _i Effective tranport number of species at polarized surface - V Average fluid velocity, m s^–1 - x Distance from inception of concentration disturbance, m - ₁ Electrical power conversion efficiency - ₂ Pumping power conversion efficiency - Solution viscosity, kg (m s)^–1 - Solution density, kg m^–3 Dimemionless groups Fanning friction factor - Reynolds number - R h/g Channel aspect ratio - D _e/l Promoter frequency - S/S ₀ Contraction coefficient - Sherwood number - Degree of reaction - Dimensionless total sensitive - Dimensionless current density - Energy cost ratio     

Kinetics of copper electrodeposition in citrate electrolytes

The kinetics of copper electrocrystallization in citrate electrolytes (0.5M CuSO₄, 0.01 to 2M sodium citrate) and citrate ammonia electrolytes (up to pH 10.5) were investigated. The addition of citrate strongly inhibits the copper reduction. For citrate concentrations ranging from 0.6 to 0.8 M, the impedance plots exhibit two separate capacitive features. The low frequency loop has a characteristic frequency which depends mainly on the electrode rotation speed. Its size increases with increasing current density or citrate concentration and decreases with increasing electrode rotation speed. A reaction path is proposed to account for the main features of the reduction kinetics (polarization curves, current dependence of the current efficiency and impedance plots) observed in the range 0.5 to 0.8 M citrate concentrations. This involves the reduction of cupric complex species into a compound that can be either included as a whole into the deposit or decomplexed to produce the metal deposit. The resulting excess free complexing ions at the interface would adsorb and inhibit the reduction of complexed species. With a charge transfer reaction occurring in two steps coupled by the soluble Cu(I) intermediate which is able to diffuse into the solution, this model can also account for the low current efficiencies observed in citrate ammonia electrolytes and their dependencies upon the current density and electrode rotation speed.Nomenclature b, b ₁, b ₁ ^* Tafel coefficients (V^–1) - bulk concentration of complexed species (mol cm^–3) - (si*) concentration of intermediate C^* atx=0 (mol cm^–3) - C concentration of (Cu Cit H)^2– atx=0 (mol cm^–3) - C C variation due to E - C concentration of complexing agent (Cit)^3- at the distancex (mol cm^–3) - C _o concentrationC atx=0 (mol cm^–3) - C _o C _o variation due to E - Cv _s bulk concentrationC (mol cm^–3) - (Cit H), (Cu), (Compl) molecular weights (g) - C _dl double layer capacitance (F cm^–2) - D diffusion coefficient of (Cit)^3- (cm²s^–1) - D ₁ diffusion coefficient of C^* (cm²s^–1) - E electrode potential (V) - f ₁ frequency in Equation 25 (s^–1) - F Faraday's constant (96 500 A smol^–1) - i, i ₁, i ₁ ^* current densities (A cm^–2) - i i variation due to E - Im(Z) imaginary part ofZ - j - k ₁, k ₁ ^* , K₁, K ₁ ^* , K₂, K rate constants (cms^–1) - K rate constant (s^–1) - K ₃ rate constant (cm³ A^–1s^–1) - R _t transfer resistance (cm²) - R _p polarization resistance (cm²) - Re(Z) real part ofZ - t time (s) - x distance from the electrode (cm) - Z _f faradaic impedance (cm²) - Z electrode impedance (cm²) Greek symbols maximal surface concentration of complexing species (molcm^–2) - thickness of Nernst diffusion layer (cm) - , ₁, ₂ current efficiencies - angular frequency (rads^–1) - electrode rotation speed (revmin^–1) - =K ^–1(s) - _d diffusion time constant (s) - electrode coverage by adsorbed complexing species - (in0) electrode coverage due toC _s - variation due to E     

Recovery of pickling effluents by electrochemical oxidation of ferrous to ferric chloride

The characteristics of the effluents from the preparatory pickling step of zinc plating are presented and the various methods of oxidizing ferrous to ferric chloride are briefly considered. An electrochemical oxidation method is proposed to recover these effluents by using an electrochemical cell with three-dimensional electrodes and an anion selective membrane. A near exhausted hydrochloric acid solution was used as catholyte. The experimental data obtained from the proposed cell show a faradic yield of 100% and easy control of the parasitic reactions. The three-dimensional anode was modelled and it is shown that at high values of current only the felt entrance region works efficiently.Nomenclature A membrane surface (cm²) - a specific felt surface (cm^–1) - C concentration difference (mol dm-^–3) - D average diffusion coefficient through the membrane (cm² s^–1) - i _n felt wall flux of species (mol cm^–2 s^–1) - j total current density (A cm^–2) - j ₀ exchange current density (A cm^–2) - j ₁ current density in matrix (A cm^–2) - j ₂ current density in felt solution (A cm^–2) - j _n transfer current density (A cm^–2) - L thickness of felt electrode (cm) - L _m thickness of membrane (cm) - M transport of ferrous and ferric ions through the membrane (mol) - N superficial flux of ion reactant (mol cm^–2 s^–1) - u superficial fluid velocity (cm s^–1) - x distance through felt electrode (cm) - R universal gas constant (8.3143 J mol^–1 K^–1) - T absolute temperature (K) - t time (s) Greek letters _a, _c anodic and cathodic transfer coefficient - local overpotential ( = ₁ – ₂) (V) - conductivity of solution (mS cm^–1) - µ solution viscosity (Pa s) - solution density (g cm^–3) - conductivity of solid matrix (mS cm^–1) - ₁ electrostatic potential in matrix phase (V) - ₂ electrostatic potential in solution (V)     

Impedance of a porous electrode with an axial gradient of concentration

When the impedance is measured on a battery, an inductive impedance is often observed in a high frequency range. This inductance is frequently related to the cell geometry and electrical leads. However, certain authors claimed that this inductance is due to the concentration distribution of reacting species through the pores of battery electrodes. Their argument is based on a paper in which a fundamental error was committed. Hence, the impedance is re-calculated on the basis of the same principle. The model shows that though the diffusion process plays an outstanding role, the overall reaction rate is never completely limited by this process. The faradaic impedance due to the concentration distribution is capacitive. Therefore, the inductive impedance observed on battery systems cannot be, by any means, attributed to the concentration distribution inside the pores. Little frequency distribution is found and the impedance is close to a semi-circle. Therefore depressed impedance diagrams in porous electrodes without forced convection cannot be ascribed to either a Warburg nor a Warburg-de Levie behaviour.Nomenclature A D¦C¦ (mole cm s^–1) - B j+K¦C¦ (mole cm s^–1) - b Tafel coefficient (V^–1) - C(x) Concentration ofS in a pore at depthx (mole cm^–3) - C ₀ Concentration ofS in the solution bulk (mole cm^–3) - C C(x) change under a voltage perturbation (mole cm^–3) - ¦C¦ Amplitude of C (mole cm^–3) - D Diffusion coefficient (cm² s^–1) - E Electrode potential (V) - E Small perturbation inE namely a sine-wave signal (V) - ¦E¦ Amplitude of E(V) - F Faraday constant (96500 A s mol^–1) - F(x) Space separate variable forC - f Frequency in Hz (s^–1) - g(x) KC(x)¦E¦(mole cm s^–1) - I Apparent current density (A cm^–2) - I _st Steady-state current per unit surface of pore aperture (A cm^–2) - j Imaginary unit [(–1)^1/2] - K Pseudo-homogeneous rate constant (s^–1) - K Potential derivative ofK, dK/dE (s^–1 V^–1) - K ^* Heterogeneous reaction rate constant (cm s^–1) - L Pore depth (cm) - n Reaction order - P Reaction product - p Parameter forF(x), see Equation 13 - q Parameter forF(x), see Equation 13 - R _e Electrolyte resistance (ohm cm) - R _p Polarization resistance per unit surface of pore aperture (ohm cm²) - R _t Charge transfer resistance per unit surface of pore aperture (ohm cm²) - S Reacting species - S _a Total surface of pore apertures (cm²) - S ₀ Geometrical surface area - S _p Developed surface area of porous electrode per unit volume (cm² cm^–3) - s Concentration gradient (mole cm^–3 cm^–1) - t Time - U Ohmic drop - x Distance from pore aperture (cm) - Z Faradaic impedance per unit surface of pore aperture (ohm cm²) - Z _x Local impedance per unit pore length (ohm cm³) - z Charge transfer number - Porosity - Thickness of Nernst diffusion layer - Penetration depth of reacting species (cm) - Penetration depth of a.c. signal determined by the potential distribution (cm) - Electrolyte (solution) resistivity (ohm cm) - ₀ Flow of S at the pore aperture (mole cm² s^–1) - Angular freqeuncy of a.c. signal, 2f(s^–1) - Integration constant     

Optimization of an electrochemical hydrogen-chlorine energy storage system

A mathematical model is presented for the optimization of the hydrogen-chlorine energy storage system. Numerical calculations have been made for a 20 MW plant being operated with a cycle of 10 h charge and 10h discharge. Optimal operating parameters, such as electrolyte concentration, cell temperature and current densities, are determined to minimize the investment of capital equipment.Nomenclature A _ex design heat transfer area of heat exchanger (m²) - a _F electrode area (m²) - heat capacity of liquid chlorine (J kg^–1K^–1) - heat capacity of hydrogen gas at constant volume (J kg^–1 K^–1) - c _p,hcl heat capacity of aqueous HCl (J kg^–1 K^–1) - C _$acid cost coefficient of HCl/Cl₂ storage ($ m^–1.4) - C _$ex cost coefficient of heat exchanger ($ m^–1.9) - C _$F cost coefficient of cell stack ($ m^–2) - cost coefficient of H₂ storage ($ m^–1.6) - C _$j cost coefficient of equipmentj ($/unit capacity) - C _$pipe cost coefficient of pipe ($ m^–1) - C _$pump cost coefficient of pump ($ J^–0.98 s^–0.98) - E cell voltage (V) - F Faraday constant (9.65 × 10⁷ C kg-equiv^–1) - F _j design capacity of equipmentj (unit capacity) - G _D design electrolyte flow rate (m³ h^–1) - heat of formation of liquid chlorine (J kg-mol^–1 C1₂) - H _f ⁰ ,_HCl heat of formation of aqueous HCl (J kg-mol^–1HCl) - H _m total mechanical energy losses (J) - I total current flow through cell (A) - i operating current density of cell stack (A m^–2) - L length of pipeline (m) - N number of parallel pipelines - n_HCl change in the amount of HCl (kg-mole) - P pressure of HCl/Cl₂ storage (kPa) - p ₁ H₂ storage pressure at the beginning of charge (kPa) - p ₂ H₂ storage pressure at the end of charge (kPa) - –Q _ex heat removed through the heat exchanger (J) - R universal gas constant (8314 J kg-mol^–1 K^–1) - the solubility of chlorine in aqueous HCl (kg-mole Cl₂ m^–3 solution) - T electrolyte temperature (K) - T ₂ electrolyte temperature at the end of charge (K) - T _max maximum electrolyte temperature (K) - T _min minimum electrolyte temperature (K) - t final time (h) - t _ex the length of time for the heat exchanger operation (h) - Uit _ex overall heat transfer coefficient (J h^–1 m^–2 K^–1) - V _acid volume of HCl/Cl₂ storage (m³) - } volume of H₂ storage (m³) - v design linear velocity of electrolyte (m s^–1) - amount of liquid chloride at timet (kg) - amount of liquid chlorine at timet ₀ (kg) - w _hcl amount of aqueous HCl solution at timet (kg) - W _p design brake power of pump (J s^–1) - X electrolyte concentration of HCl at timet (wt fraction) - X _f electrolyte concentration of HCl at the end of charge (wt fraction) - X _i electrolyte concentration of HCl at the beginning of charge (wt fraction) - X ₀ electrolyte concentration of HCl at timet ₀ (wt fraction) - Y objective function to be minimized ($ kW^–1 h^–1) - _j the scale-up exponent of equipmentj - overall electric-to-electric efficiency (%) - _acid safety factor of HCl/Cl₂ storage - fractional excess of liquid chlorine - _p pump efficiency - average density of HCl solution over the discharge period (kg m^–3)     

Electrochemical machining of refractory materials

The feasibility of the electrochemical machining (ECM) of pure TiC, ZrC, TiB₂ and ZrB₂ has been established. In addition, the ECM behaviour of a cemented TiC/10% Ni composite has been investigated and compared to that of its components, TiC and nickel. ECM was carried out in 2M KNO₃ and in 3 M NaCl at applied voltages of 10–31 V and current densities of 15–115 A cm^–2. Post-ECM surface studies on the TiC/Ni composite showed preferential dissolution of the TiC phase during machining.Nomenclature E ₀ thermodynamic equilibrium potential (V) - F Faraday's constant (96 500 Coul mol^–1) - toolpiece feed rate (cm s^–1 or mm min^–1) - I current (A) - i current density (A cm^–2) - k electrolyte conductivity (^–1 cm^–1) - l interelectrode gap (mm) - mass removal rate (g s^–1 or g min^–1) - M formula weight (g mol^–1) - Q electrolyte flow rate (l min^–1) - t electrolyte temperature (°C) - V applied voltage (V) - V _IR ohmic drop through electrolyte (V) - z apparent valence of dissolution (eq mol^–1) - _i overvoltages (V) - density of refractory materials (g cm^–3)     

The voidage problem in gas-electrolyte dispersions

Assessing the ohmic interelectrode resistance of electrochemical reactors with gas evolution requires data for the gas void fraction of gas-electrolyte dispersions. A voidage equation is derived taking account of the internal liquid flow in stationary electrolytes and at small liquid superficial velocities. The equation is a general form of available voidage equations.Nomenclature C non-dimensional constant, Equation 8 - n exponent, Equation 5 - S cross-sectional area (m²) - v _G gas velocity (m s^–1) - v _L liquid velocity (m s^–1) - v _s rising velocity of a bubble swarm (m s^–1) - v _l terminal rising velocity of a single bubble (m s^–1) - V_G volume flow rate of gas (m³ s^–1) - V_L volume flow rate of liquid (m³ s^–1) - fraction of cross-sectional area - volume (void) fraction of gas - _m geometric maximum of void fraction - maximum of void fraction in infinite gas flow Indices i internal - t total     

Living cationic polymerization of a vinyl ether with a triester pendant

Summary Cationic polymerization of CH₂=CH-O-CH₂CH₂C(COOC₂H₅)₃, a vinyl ether with three pendent esters, initiated by the HI/I₂ system in toluene at –40 °C afforded living polymers with a controlled molecular weight ( = 10³–10⁴) and a narrow molecular weight distribution ( = 1.1–1.2). The number-average molecular weight of the polymers was directly proportional to monomer conversion and the monomer-to-initiator (HI) feed ratio. The polymers obtained with BF₃O(C₂H₅)₂ had a fairly high molecular weight ( 10⁵, 5 × 10⁴) and a broad molecular weight distribution. The triester vinyl ether was similar in reactivity to alkyl counterparts and one order of magnitude more reactive than the corresponding mono- and diester vinyl ethers.     

Computer simulation of high-discharge-rate battery systems

Simulations were carried out for a proposed two-dimensional high-discharge-rate cell under load with an interelectrode gap of the order of 100 m. A finite difference program was written to solve the set of coupled, partial differential equations governing the behaviour of this system. Cell dimensions, cell loads, and kinetic parameters were varied to study the effects on voltage, current and specific energy. Trends in cell performance are noted, and suggestions are made for development of cells to meet specific design criteria. Modelling difficulties are discussed and suggestions are made for improvement.Nomenclature A surface area of unit cell (cm²) - A _k conductivity parameter (cm^{2 –1} mol^–1) - b Tafel slope (V) - c concentration (mol cm^–3) - c ₀ concentration of bulk electrolyte (mol cm^–3) - D diffusivity (cm² s^–1) - D _h lumped diffusion parameter (J s cm^–2 mol^–1) - D _s lumped diffusion coefficient (A cm² mol^–1) - E rest potential of electrode (V) - F Faraday constant (96 500 C mol^–1) - i current density (A cm^–2) - I total current for unit cell (A) - i ₀ exchange current density (A cm^–2) - N flux of charged species (mol cm² s^–1) - R gas constant (8.314 J mol^–1 K^–1) - R _ext resistance external to cell () - t time (s) - T temperature (K) - t ₀ transference number - u mobility (cm² mol J^–1 s^–1) - V volume of an element in the cell (cm³) - V _ext voltage external to cell (V) - z charge on an ion - _c concentration overpotential (V) - _s surface overpotential (V) - conductivity (^–1 cm^–1) - stoichiometric coefficient - electric potential in solution (V)     

Flow-through and flow-by porous electrodes of nickel foam. I. Material characterization

This paper deals with the characterization of three nickel foams for use as materials for flow-through or flow-by porous electrodes. Optical and scanning electron microscope observations were used to examine the pore size distribution. The overall, apparent electrical resistivity of the reticulated skeleton was measured. The BET method and the liquid permeametry method were used to determine the specific surface area, the values of which are compared with those known for other materials.Nomenclature a _e specific surface area (per unit of total volume) (m^–1) - a _s specific surface area (per unit of solid volume) (m^–1) - (a _e)_BET specific surface area determined by the BET method (m^–1) - (a _e)_Ergun specific surface area determined by pressure drop measurements (m^–1) - mean pore diameter (m) - mean pore diameter determined by optical microscopy (m) - mean pore diameter using Ergun equation (m) - e thickness of the skeleton element of the foam (m) - G grade of the foam (number of pores per inch) - P/H pressure drop per unit height of the foam (Pa m^–1) - r electrical resistivity ( m) - R _h hydraulic pore radius (m) - T tortuosity - mean liquid velocity (m s^–1) Greek symbols mean porosity - circularity factor - dynamic viscosity (kg m^–1 s^–1) - liquid density (kg m^–3) - pore diameter size dispersion     

Mathematical modelling of electrode growth

It is known that during electrodeposition or dissolution electrode shape change depends on the local current density (Faraday's law in differential form). Assuming that concentration gradients in the bulk of the solution may be neglected, the current distribution in an electrochemical system can be modelled by a Laplace equation (describing charge transport) with nonlinear boundary conditions caused by activation and concentration overpotentials on the electrodes. To solve this numerical problem, an Euler scheme is used for the integration of Faraday's law with respect to time and the field equation is discretized using the boundary element method (BEM). In this way, and by means of a specially developed electrode growth algorithm, it is possible to simulate electrodeposition or electrode dissolution. In particular, attention is paid to electrode variation in the vicinity of singularities. It is pointed out that the angle of incidence between an electrode and an adjacent insulator becomes right (/2). This is confirmed by several experiments.List of symbols x _i coordinates of a point i belonging to a boundary (m) - t time (s) - h thickness variation at a point belonging to an electrode (m) - M molecular weight (kgmol^–1) - m specific weight (kgm^–3) - z charge of an ion (C) - F Faraday's constant (C mol^–1) - R _a2 impedance of the linearized activation overvoltage on cathode (S2 cm^–2) - efficiency of the reaction - electric conductivity (^–1 m^–1) - U electric potential (V) - rate of mechanical displacement of a point (m s^–1) - V applied potential on an electrode (V) - W Wagner number defined as the ratio of the mean impedance of the reaction and the mean ohmic resistance of the cell given by L/ with L a characteristic length of the cell. - overvoltage (V) - ₁ overvoltage on anode (V) - ₂ overvoltage on cathode (V)