Wright 93
Wright 93
Wright 93
Smelting Cells
by
A. Wright B.Sc.
June, 1993
ABSTRACT
I would like to express my sincere thanks to the following people who have
enabled me to carry out this research:
The staff of the Department of Chemical and Process Engineering for their
help and advice.
George Guelfo, David Salt, Tony Sargent and Mike Wilson of Anglesey
Aluminium for their support and practical assistance in developing the
model and their patience when showing me the ropes of aluminium
production.
The technical staff of Anglesey Aluminium also for their advice and
assistance.
Finally, I would like to thank my wife, Suzanne for her moral support and
encouragement during the overly long course of this research.
To my wife, Suzanne and our son, Toseland ...
Contents
Page
1.0 History of the Hall-Heroult Process 1
1.1 Main Costs Of Production 1
1.1.1 Fundamental relationships 2
1.2 Types of Cell 3
1.3 Operation of a Prebake Cell 6
1.3.1 Data Sampling 6
1.3.2 Alumina feed 7
1.3.3 Voltage control 7
1.3.4 Anode effect termination 8
1.3.5 Metal tapping 9
1.3.6 Addition of raw materials other than alumina 10
1.3.7 Anode change 11
1.3.8 Sampling 12
II
Page
4.2.5 Stochastic disturbances 77
4.3 Alumina Dissolution 78
4.4 Electrolyte Mass Balance 81
4.4.1 Bath volatilisation 82
4.4.2 Gaseous emissions 82
4.4.3 Particulate emissions by volatilisation 84
4.4.4 Aluminium dissolution and back reaction 88
4.4.5 Current efficiency during anode effect 91
4.4.6 Current efficiency losses at the cathode 93
4.4.7 Overall cu rrent efficiency 94
4.4.8 Mass balances 95
4.4.9 Material in suspension and as sludge 95
4.4.10 A13+ ions in solution 97
4.4.11 0 2 - ions in solution 98
4.4.12 Aluminium metal in electrolyte 99
4.4.13 Aluminium in the metal pad 99
4.4.14 Sodium fluoride 100
4.4.15 Calcium fluoride 101
4.4.16 Alumina in solution 102
4.4.17 Aluminium fluoride 102
4.4.18 Overall mass balance on electrolyte 103
4.5 Electrolyte Conductivity 103
4.6 Electrolysis Charge Balance: Reaction Modelling 105
4.6.1 Determination of Tafel constants and mass
transfer coefficient 111
4.6.2 Reactor modelling 114
4.6.3 Anode effect 117
4.7 Cell Voltage, Energy Consumption and Cell
Resistance 120
4.8 Heat Balance 121
4.9 Electrolysis Heat Balance 122
III
Page
IV
Page
v
Page
5.23.2 Point feed cells: adaptive feed strategy 168
5.23.3 Search routine 169
5.23.4 Overfeed/Underfeed strategy (O/U strategy) 169
5.23.5 Starve routine 170
5.24 ACD Control Strategies 171
5.24.1 Fixed ACD 171
5.24.2 ACD optimisation 171
5.25 Summary 172
VI
Page
Appendices 203
References 220
Nomenclature 229
VII
List of Figures
Figure Following
number page
1.1 Cost of Production against Time 1
1.2 Schematic Diagram of the Hall Heroult Cell 4
1.3 Schematic Diagram of the Soderberg Cell 4
1.4 Anode Effect Termination Sequences 9
3.1 Schematic of the Current Efficiency Model 38
3.2 Cell Current Efficiency vs Alumina Concentration 44
3.3 Aluminium Concentration across CO 2 Boundary
Layer 49
3.4 Heat Flow and Temperature Gradient across Frozen
Sidewalls 51
3.5 Variation of Freeze Thickness with Parameter f 53
3.6 Local Nusselt Numbers as a Function of Local
Rayleigh Number 57
3.7 Fanning Factors for Determination of Effective
Anode Area 59
3.8 Schematic Diagram of the Model for Resistance
across the Bath 60
3.9 Models for Electrolyte Conductivity 62
3.10 Results of Fellners Simple Model for Binary Systems 63
3.11 Results of Fellners Simple Model of Tertiary System 63
4.1 Conceptual Model of the Aluminium Reduction Cell 75
4.2 Electrolyte Mass Balance Processes 81
4.3 Random Variation on Mass of Alumina Fed to a Cell 97
4.4 Schematic Representation of Reactions at an
Electrode 105
4.5 Schematic diagram of Tafel Derivation i+ »i- 106
VIII
Following
page
4.6 Overvoltages from Zuca 111
4.7 Overpotential vs Current Density 113
4.8 Overpotential vs Alumina Concentration 114
4.9 Overpotential by Mixed Model vs Alumina
Concentration 114
4.10 Current Density-Potential Curve for a Mixture of Two
Compounds Oxidising at Different Potentials 118
4.11 Schematic Diagram of Heat Flow in an Aluminium
Cell 121
4.12 Heat balance across the cathode 125
4.13 Heat balance across an anode 129
4.14 Cell Sidewall Profiles in Old and New Cells 138
4.15 Current Drawn by a Cold Anode 143
5.1 CElTROl distributed control system 156
5.2 Voltage Control Block Diagram 160
5.3 Anode Effect Termination Sequences 165
5.4 Cell Current Efficiency vs. Alumina Concentration 166
IX
Following
page
6.5 Simulated Ratio vs Temperature for 151 kA Centre
Break Cell 180
6.6 Simulated Bath and Pad Depths for 151 kA Centre
Break Cell 181
6.7 Alumina Concentration over First 24 Hours of the
'Standard' Run 182
6.8 Measured Variation of Alumina Content in Reduction
Cell Bath 182
6.9 Measured Variation of Alumina Content in Reduction
Cell Bath 182
6.10 Reduction Cell Operation with Two-Hour Break/Feed
Cycle 182
6.11 Simulation of Alumina Concentration, Freeze
Thickness and Cell Voltage 182
6.12 Simulated Effect of Anode Change upon Alumina
Concentration 183
6.13 Comparison between Simulated and Actual Voltage
Trace During Anode Effect in Cell 2A 18 183
6.14 Voltage Trace of Anode Effect Predicted by Model 184
x
Following
page
6.23 Typical Temperature-Time Curves for Two Different
Feed Modes 186
6.24 Simulated Temperature Profile over a 67 Minute
Feed Cycle 186
6.25 Average Temperature in an Industrial Cell After a
0.1V Increase in Operating Voltage 186
6.26 Temperature in an Industrial Cell During Anode
Effect 189
XI
List of Tables
Table Following
number page
1.1 Major Cost Centres in Aluminium Production 1
3.1 Viewpoints on the Rate Determining Step for
Alumina Dissolution 37
3.2 Current Efficiency vs Alumina Concentration:
Literature Data 45
3.3 Commonly Used Dimensionless Groups for
Determination of Heat Transfer Coefficients 54
4.1 Experimental Tafel Coefficients for Anodes of
Varying Porosity 111
XII
1.0 History of the Hall-Heroult Process
The electrolytic process for the industrial production of aluminium was
patented in 1886 by P.L.T Heroult in France and, independently, by C.M.
Hall in the U.S.A. The process involves the electrolysis of a molten
solution of alumina in cryolite and utilises the reaction
Figure 1.1 shows how the cost of production and the price commanded for
aluminium has changed over the 55 years prior to 1985,[1] It is clear that
overall there has been a trend of decreasing real profitability since 1930.
This is not due to a similar trend in capital costs but rather to the increase
in world production especially due to the building of large smelting plants
in developing nations to utilise cheap electricity and labour.
1
1890 1910 1930 1950 1970
Alumina 35
Electrical Energy 25
Labour 10
Electrolyte Materials 3
Cathode Replacement 2
It can be seen that the cost of the alumina fed to the cell constitutes the
largest percentage of the cost of producing aluminium. In the case of
Anglesey Aluminium, the alumina is purchased upon the open market
and, although alumina quality affects the operation of the cell, control over
alumina costs does not rest directly with the smelter.
1. 1. 1 Fundamental relationships
There are three fundamental relationships that dictate the way in which an
aluminium cell should be operated. These are
2
increasing current and increasing current efficiency gives increased
productivity
I i + CE i ~ productivity i
CE i + V .t ~ energy .t
cd i ~ energy i
The economic climate will dictate which of these relationships are the
most important. At times when aluminium can command a high price it is
desirable to increase productivity at the expense of energy efficiency,
while at times when the price of aluminium is low it is most important to
reduce production costs by increasing energy efficiency. Any control
strategy must be able to accommodate both these situations without the
need for major changes to the hardware.
3
The aluminium cell or IIpotll consists of a steel box lined with a refractory
brick. A carbon cathode block lies on top of this and forms the floor of the
cell. The cell is partially filled with the molten electrolyte and carbon
anodes are immersed in the electrolyte. A layer of molten aluminium, the
ll
IIpad forms at the surface of the cathode block protecting the carbon from
the corrosive action of the electrolyte and taking over the role of the
cathode.
The sidewalls of the cell are designed to allow a crust or IIledge ll of frozen
cryolite to form around the cell. This contains the molten electrolyte and is
very important in protecting the cell lining from corrosion and in
maintaining the thermal stability of the cell.
4
Fume Collection Anode Bus Bar
Removeable Covers
/
Carbon Anode
Graphite Rammed
Sidewall Carbon
Current
Collector Bar
Thermal Insulation
Thermal Insulation
5
1.3 Operation of a Prebake Cell
Measured data
Measurable data includes such information as electrolyte composition,
cell voltage and electrolyte temperature. These can be measured with
different degrees of accuracy and frequency. For example, cell voltage
can be measured almost continuously and with good accuracy. At
Anglesey Aluminium, the control computer, GELTROl, samples
individual cell voltages at a 100 ms interval to 4 significant figures
accuracy. Bath composition, on the other hand, is measured by mass
spectroscopy methods which are time consuming and cannot be
performed continuously. Although the accuracy of this method is good
for individual samples, it is not necessarily representative of the whole
cell.
Observed data
Appearance of the cell. eg flame colour, crust cover, movement of the
electrolyte. This type of data is subjective and depends upon the
experience of the operator for its usefulness.
6
1.3.2 Alumina feed
Alumina is fed to a cell by breaking the surface crust and allowing the
crust along with a mass of alumina to fall into the electrolyte. The alumina
on the surface is then replenished and a new crust allowed to form.
In early cells, the crust was manually broken using metal bars and alumina
was shovelled onto the surface. The process was first mechanised by the
development of large vehicular crust breakers that were driven along the
sides of the cells breaking the crust by means of a cogged wheel. Alumina
was replenished from a following truck.
7
1.3.4 Anode effect termination
Anode effect occurs when the alumina concentration in the electrolyte falls
below a critical level. During anode effect, the composition of the bath
results in dewetting of the anodes and the formation of large bubbles on
the anode surface. The bubbles are predominately carbon tetrafluoride
and present a physical barrier to the passage of current into the
electrolyte. Anode effect cannot, therefore, be terminated by the simple
addition of alumina to the melt thus increasing the alumina concentration
above the critical level, but requires the breaking of the gas bubble layer
for alumina addition to be effective. This may be achieved by a number of
methods.
raise anode
lower anode
break duct end
break tap end
feed duct end
feed tap end
8
A typical sequence for a half break cell is shown in figure 1.4. Actions may
be repeated a number of times by means of a loop construct or bypassed
on the condition of the anode effect being terminated by means of a jump.
These sequences are very effective with a 95% kill efficiency on point feed
cells. If the sequence fails in terminating anode effect, however, the
traditional method of using a green pole is still used.
9
0 1 2 3 4 5 6 7 8 9
0 BRAK LOWR 7.0 W006 JUMP 15 RAIS 5.0 DUMP BRAK
1 LOWR 5.0 W006 LOPO 1 6 -RAIS 3.0 WOOS RAIS 3.0
2 W006 BRAK WOOS DUMP DUMP DUMP DUMP
3
4 WOOS BRAK W010
5 DUMP DUMP DUMP DUMP
0 1 2 3 4 5 6 7 98
0 DUMP BRAK LOWR 2.0 DUMP BRAK JUMP 13 LOWR 2.0
1 DUMP BRAK W001 LOPO 18 JUMP 24 RAIS 3.0 DUMP
2 BRAK LOWR 3.0 W001 LOP1 217 BRAK DUMP -RAIS 2.0
3 BRAK DUMP WOO2 RAIS 2.0 W002 RAIS 0.0 -W04S BRAK
4 DUMP DUMP
5
Nomenclature
1. Sodium oxide in the alumina feed reacts with aluminium fluoride in the
melt to form sodium fluoride.
3. Both sodium fluoride and aluminium fluoride are volatilised off the bath
surface. This becomes more severe at higher temperatures.
4. The melting and freezing of the ledge causes changes in the ratio. This
is because the electrolyte freezes with a ratio closer to that of pure
cryolite at about 1.40. In an electrolyte with a ratio less than that of
cryolite this effectively means that sodium fluoride freezes preferentially
to aluminium fluoride with the result that freezing causes the ratio to
decrease and vice versa.
This tends to result in a ratio greater than that desired and it is necessary
to add aluminium fluoride to the bath. This is simply done by adding
bagged aluminium fluoride to the cell at a fixed interval. A typical example
would be 25 kg aluminium fluoride per day. The bag of aluminium fluoride
is simply thrown into the centre trench of the cell and broken into the
electrolyte by the action of the crust breakers.
The overall mass of electrolyte in the cell may also need replenishing. This
can be due to volatilisation of the electrolyte in a hot or poorly covered cell
or due to the splashing of bath out of the cell by anode movements. The
need to add electrolyte is not common in mature cells but is necessary in
cells shortly after start up. This is done by siphoning electrolyte from a
10
mature cell using the same method and crucible as for metal tapping and
blowing it into the receiving cell using compressed air.
A spent anode is pulled from the cell by a small Icrane ' . If necessary, the
hole is then enlarged to take the new anode by drilling out the frozen
sidewall. The new cold anode is then placed in the electrolyte and bolted
to the bus bar. The positioning of the anode on the bus bar is adjusted to
give a uniform inter electrode gap within the cell. As the metal pad has a
humped profile due to the electromagnetic forces in the cell, the depth of
the anode relative to the bus bar is determined by its position within the
cell.
A new anode is much colder than the electrolyte into which it is being
placed, about 200 0
e compared to 970 o
e. A layer of frozen electrolyte
immediately forms around the anode effectively insulating it both
electrically and thermally. This frozen layer melts as the anode warms
over the next 30 hours but in order to increase the speed of melting it is
common practise to set a new anode low so that the frozen layer can
contact the metal pad. This prevents the freeze from growing and heat is
then conducted more speedily through the freeze and into the anode. On
the day following the anode change, the new anode is raised to the
correct position.
11
1.3.8 Sampling
Liquid levels are determined on a daily basis by inserting an iron rod into
the bath and metal. When the rod is removed from the cell, the aluminium
and electrolyte freeze, the metal forming a thinner layer. The depth of
each section may then be read off against a rule on the side of the cell
and the values used to decide the amount of metal that should be tapped.
12
expensive. If a thermocouple fails, it cannot be replaced, and some means
of discarding the results from failed thermocouples is therefore necessary.
The severity indices for individual cells are stored and can be used In
The potline survey can be carried out by the line supervisor and cell
operators and includes an assessment of the flame coloration, pot voltage
and stability, carbon dusting in the bath and liquid levels. A severity index
of 0.5 to 3 is assigned to a cell having one or more abnormal attributes
and the required corrective action may then be taken.
References
13
3 Manaktala, S., "Process Control Techniques for Reduction Cell
Operation", AIME Light Metals, 1971, pp165
14
2.0 Process Modelling
The usual method for evaluating the effect of changes in operating
parameters upon a system is by actual physical experimentation. This
type of experimentation can be very successful but its range is limited by
cost and safety. An alternative approach is to construct a mathematical
model of the process and experiment with the model instead of reality.
This has a number of advantages:
2) Economics
In industry, the cost of raw materials, labour etc can restrict extensive
experimentation. A suitably detailed process model however, is able to
quantitatively describe the behaviour of a system under a variety of
different operating conditions. Simulation studies can provide large
amounts of data at little cost and in a fraction of the time required for
experimentation. The cost of simulation of an experiment that results in
failure is no more than for one that is successful.
15
4) Safety
5) Computer control
The advent of cheap powerful computers in the past few years has
permitted vastly more detailed and realistic analysis to be carried out for a
reasonable degree of effort and cost. More complex and fundamental
physical principles have been employed to formulate mathematical models
of processes. This has resulted in many processes being comprehensively
better monitored and controlled and has allowed control systems to be
designed without undue disturbance to the process.
6) Operator training
Process modelling makes possible the training of process operators using
simulation rather than real plant. Operators can gain first hand experience
in dealing with conditions that may occur only infrequently In process
operation as well as the implementation of safety procedures.
16
2) Model description
4) Model evaluation
5) Application
If the model is able to accurately describe the system, it can then be used
to simulate process behaviour under a wide range of conditions. Any
assumptions made in development of the model must be taken into
consideration, however, as the model may not be valid if these do not
apply.
In practise, real problems are often ill defined and difficult to specify. Lack
of experimental data often makes it difficult to construct a model that
provides a balanced image of reality. Too simple a model may fail
because it does not adequately describe the process, while too complex a
model will contain a great number of unknown parameters and variables
and may not be mathematically tractable. The complexity of model
constructed therefore depends upon the accuracy required and the
application for which the model is intended.
17
2.2 Types of Model
18
The complexity of distributed parameter models may be reduced by
reducing the model dimensionally. In the tubular reactor example, if flow
through the reactor is turbulent then the system may be assumed well
mixed in the radial direction. Only one set of differential equations,
representing conditions along the length of the reactor, need then be
solved. What is lost in neglecting the radial distribution may be gained in
the simplicity of the resulting model.
For example, static models have been developed for the cathode lining of
an aluminium cell. These models are complex and can accurately predict
the temperature profiles for a particular set of conditions. The model is
used to design cathode linings based upon average conditions expected
in operating the cell.
Dynamic models on the other hand, often sacrifice high levels of detail
accorded by multi dimensional distributed models for the sake of both
mathematical tractability and computation time. An obvious example of
this approach is that used in classical control system analysis where linear
ordinary differential equations are used throughout. These equations are
often 'fitted ' to process responses to known disturbances. Although this
approach may provide sufficient detail for controller tuning, it provides little
data for other engineering design.
19
2.3 Choice of Model
References
20
3.0 Mathematical Models of Aluminium Reduction Cells
Process models may be divided into two main groups, static and dynamic
models. A static model represents a steady state system, or a system at a
given instant in time, for a defined average of states.
For example, many static models exist for the heat balance in a cell. In
fact the heat flows will vary with changes in temperature, instantaneous
cell voltage, electrolyte composition and depth, cell design and ambient
temperature. These models assume constant values for all such
parameters.
Whilst this simple approach is adequate for the investigation of the long
term effects of changes within some operational parameters it cannot be
used to describe the dynamic effect of a cell to frequent changes such as
alumina feed and anode effect.
21
both the electrolyte and metal pad and peripheral heat flows. The pot is
represented by a cross section of unit thickness and divided into a number
of elements. In the particular instance published, 1683 intersections were
necessary for the required accuracy. Temperatures and potentials are
determined at each intersect by iterative calculation from an arbitrary initial
solution.
22
nodal models like that presented by Sulmont and Hudault but accurately
describes the working of the cell. The accuracy of this type of model
depends upon the accuracy of the individual assumptions and
expressions used within it.
23
3. Dissolution from a stagnant bottom layer of loose alumina or sintered
alumina (Gerlach[7])
Experiments have also been done on the effect of the alumina feeding
method. Maeda et a1l 10] did bench scale experiments using industrial
composition bath at 970 ce. Alumina was poured into the melt through a
funnel and the dissolution time measured visually. No stirring was done
and a large number of alumina samples were collected world wide.
The rate determining step for alumina dissolution in industrial cells has not
been clarified. According to the Thonstad[8] it may be
• mass transfer
• heat transfer
• chemical reaction
(1 )
1
where r = rate of dissolution of alumina g cm- 2 min-
ko = reaction rate constant cm min- 1
E = activation energy kcal mole- 1
R = universal gas constant kcal mole- 1 K-1
To = reference temperature oK
T = reaction temperature oK
24
For concentrations greater than 3-6 wt% alumina they seem to agree with
a first order model. The accuracy of measurement by both visual and
chemical methods is not very great and it may be that at low
concentrations the driving force for dissolution is such that a first order
mechanism is involved. Other workers[9] have found first order correlations
at all concentrations.
(2)
Much work has also been done on the dissolution of alumina from
sludge[14,15,16,17,18]. All workers agree that the rate of dissolution from
sludge is slower than the initial rate for dispersed alumina particles in
suspension. This would be expected as the specific area of particles in
suspension will be many times greater than that of sludge but may also be
due to sub cooling and the formation of agglomerates. The rate of
dissolution is also reduced by the presence of aluminium between the
sludge and bath forming a physical barrier. Thonstad[14] determined that
sludge dissolved by a first order reaction when exposed to bath with 4-5 %
25
conditions on the dissolution of alumina by using a voltammetric technique
for the in-situ measurement of dissolved alumina. Alumina was fed to
molten cryolite in a crucible by a variety of different methods and the
alumina concentration measured at five second intervals. The
temperatures of the bath and the crucible were also monitored. In many
cases a two stage dissolution pattern was observed, indicative of fast
initial dissolution followed by the formation, melting and dissolution of
agglomerates. The time for remelting of the agglomerates seen as a delay
in the dissolution of alumina can be calculated from
(3)
(4)
The time for remelt is highly dependent upon the mass of the agglomerate
and its position in the cell. Bagshaw has calculated times of 4-6 seconds
for individual particles exposed to turbulent bath with high heat transfer
coefficient and 10°C superheat. Under the same conditions a much larger
aggregate of 1.2 g alumina would take 230 seconds to melt. However, in a
stagnant zone with a consequently lower heat transfer coefficient this time
would be doubled. The circulation of bath would also be required to break
up the agglomerate to prevent further freezing.
26
These experiments show the validity of a mUlti-stage model. The freezing
effect may be simplified by assuming a lower specific area of alumina in
suspension than the powder to account for the presence of agglomerates.
This approach would reduce the accuracy of a multi-nodal model but will
be adequate for a well mixed model in which an average heat transfer
coefficient will be assumed throughout the bath zone. This negates the
necessity of determining the sizes and flow patterns of agglomerates.
Blake[3] and Entner et aH20] assume three states of existence. These are
as dissolved alumina, dispersed alumina or suspension and bottom
sludge. Entner assumes that alumina is first dispersed into suspension
from which it dissolves by a first order model dependent upon the mass of
alumina in suspension, the temperature of the electrolyte and the time
since anode effect. The latter factor is to account for the enhanced
dissolution of alumina due to the perturbation of the bath during anode
effect. Alumina in suspension also sediments out to form a bottom sludge
from which it dissolves into solution. The rate at which it dissolves from
sludge is unclear but appears to have the same form as dissolution from
suspension.
Blake also defines the same three states as Entner but a proportion of the
alumina introduced is added to each of these states. Mass transfer is
27
assumed to occur sequentially from sludge to suspension and from
suspension to solution. Alumina sludge is transferred to suspension
according to the relationship
(5)
Blake1s model was primarily designed for a side break cell in which
feeding takes place mainly from the sludge. In centre break cells where
the amount of sludge should not be so great and alumina crusts may form
at the metal/bath interface, the rate at which alumina can dissolve directly
from sludge may be greater than the dispersion of alumina in suspension
by the movement of electrolyte over the sludge. This would account for
the difference between the models of Blake and Entner.
Fluoride is lost from the cell during normal operation as well as by the
production of fluorocarbons during anode effect. The fluoride is often
classified as gaseous and particulate. Gaseous fluorides are those that
continue to be gases at ambient temperature, namely HF, SiF4 and
fluorocarbons. Entrained and volatilised electrolyte becomes particulate
on cooling. The exhaust gases are captured by the cells hood and passed
through the baghouse where a proportion of the fluoride is trapped and
recycled to the cells with the feed. Haupin[21J developed a mathematical
model for the prediction of fluoride evolution from Hall-Heroult cells using
thermodynamics, kinetics and the physical properties of the electrolyte.
28
3.4. 1 Model of W.E. Haupin
The water is assumed to react with AIF3 in the bath by the reaction
[1 ]
since the equilibrium constant for H2 0 reacting with AIF3 is 103 times
greater that for H2 0 reacting with Na3AIF6 and 108 times greater than that
for H2 0 reacting with NaF.
(6)
29
where PHF = partial pressure of HF gas
L\G = Gibbs free energy of reaction
R = universal gas constant
T = temperature of electrolyte
PH 20 = partial pressure of H2 0 in pot gas
aAIF3 = activity of AIF3 in electrolyte
aAI20 3 = activity of AI 2 0 3 in electrolyte
F = 2310000
G %CE
_ 955000 oR o[
P
( 8444JJ
exp 7.553- T
b 0
%H 20 in AI 20 3 + O/oH in anodesJ12 0
[ 37.44 21 .. 5
c ]-0.462
exp( 0.44199 - 3.1733R b + O. 78127R~) C:~~03 (7)
0
[ AI 20 3
30
between the model and an extrapolation of a relationship derived by linear
regression.
It has been demonstrated that the vapour over the electrolyte is NaAIF4
with a small but significant amount of dimerisation[22]
[2]
Kd = -P2
d
= exp (21414
- +15.6 J (8)
Pm T
(9)
where
2
A = 7101.6 + 3069.7· Ratio + 635.77 . Ratio + 51.22· (%UF)
764.5·(0/oAI 2 0 3 ) (0 C F )
-24 638·Ratio·(0/oUF) + )] + 13.2· Yo a 2
. [1 + 3.2029·(0/oAI 2 0 3
3
B = 7.0184 + 0.6844· Ratio + 0.08464 . Ratio + 0.01085· (%UF)
1.1385· (0/oAI 2 0 3 ) (0/ C F )
_ 0 005489· Ratio· (%UF) + )] + 0.0068· /0 a 2
. [1 + 3.2029· (0/oAI 2 0 3
31
Assuming that the pot gas is saturated with NaAIF4 and Na2AI2 Fa at their
calculated vapour pressures Haupin arrives at the following equation for
the evolution of particulate fluoride
F - 204000 (Pm + Pd )
Vp - %CE· P kg F I tonne AI (10)
b
The calculated value of particulate fluoride by this method was always less
than that measured. The difference between the calculated and measured
values was assumed to be entrained electrolyte and was correlated
against bath composition and temperature.
32
Being derived by difference this correlation contains all the errors of
measurement and the predictive error of the other models. Additionally it
predicts negative entrainments for electrolytes with bath ratios above 1.4
at temperatures greater than 990 ae. This demonstrates the danger of
extrapolating an empirical correlation.
The total fluoride evolution is given by the sum of the gaseous and
particulate evolution
(12)
This is the total fluoride evolved from the cell into the exhaust gases.
These gases are passed to scrubbers where much of the fluoride is
recaptured and recycled. The total fluoride irrevocably lost from the cell is
therefore given by
(13)
33
3) Electronic conductivity of the electrolyte and physical short circuiting
reduces the actual current used for electrolysis.
5) Metal being removed from the product before the amount formed can
be assessed. Aluminium could be reoxidised in the electrolyte either
by dissolution from the metal pad or as a dispersion of droplets.
Obviously the effect of some of these mechanisms will be small and may
be ignored in the modelling of the process. The purity of the metal
produced is typically greater than 99.8% indicating that process (1) has
only a small effect on current efficiency.
34
the rate of dissolution of aluminium carbide in the metal and electrolyte
with which it comes in contact. Although it has a significant solubility in the
electrolyte, its solubility in metal is very low. Under normal conditions the
cell sidewalls are protected by a ledge of frozen electrolyte and carbide on
the cathode surface is protected from the bath by the metal pad. The
effect upon current efficiency is therefore limited by the rate of dissolution
into the metal and is very small. If the sidewalls become exposed, as after
an extended anode effect, then aluminium carbide will form and dissolve
causing erosion of the sidewall and loss of current efficiency until the
frozen ledge is once more established. These effects are of great
significance during the very early life of a cell or in a sick pot, but can be
considered negligible during normal operation.
The effects of all the above processes are fairly small and it is universally
accepted that the irreversible removal of metal by process (5) is the
primary loss mechanism. It occurs through the dissolution of one or both
of the electrode products and transportation to a reaction zone.
+ 50 (14)
35
byproduct of the oxidation reaction. The metal reoxidation reaction IS
usually depicted
36
1. Dissolution-diffusion step. The dissolution of metal at the metal-
electrolyte interface and the diffusion of the dissolved species through
a stagnant cathodic layer. As the cathodic layer is usually very thin
these two steps may be considered together.
Qiu also lists various viewpoints presented in the literature for the rate
determining step. This is duplicated in table 3.1.
T = bath temperature °C
XA1F3 = excess AIF3 in the electrolyte wt%
37
authors year experimental method rate determining
step
laboratory cells
The model is based upon the assumption that the reoxidation reaction
occurs between dissolved CO 2 and aluminium dissolved or suspended in
the bath. The two products are brought together by diffusion and mass
transport in the bath. The model is illustrated by figure 3.1. The anode to
cathode gap is divided into five regions
38
o
oOcfooo
00 000 0
0o~ o0
000
0
0
o 00000 0 0
Q..,: 0 0
~UOOOOOO
00 0 00
o gOo 0 0
0000000 0 0
o 0
On-O
0"V0000 0 0
00
0
'0 0000 0 OJ
o ~ °O'b
00
~O0~0?:''d0000Cb~00'cX)00CfOo~62-'
0 0
o 0 "0 0 '<3 0 0
0 :l, ';;0 'I!
000.'l.J'% 00 00
o~goo
0
~ifo?~oo °00~8000000'tJ00
_'I>oooo,!;/ 0 0 0
0~n~,ectoQ
0 0 v-o ..'cf: • 2 fO 0 0 00 dO80
00 000
0 0 0
0 0 00
o Reaction 0 c ...... cGradient 0 0 0
0
00
0
0 - !O 0 o 00 0 0 0
0
0ic'-J
0 0
0000
00000
00
0
00000
0 00 0 0
0
0 0
0
Metal Transport e 0 0 0 0 0 0 0
by Convection
Boundary Layer,
Metal transport
by Diffusion p / / / / / ~/ / / / / '/ / / / /
/ / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / ,
/ / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / ,
/ / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / ,
/ / / / / / / / / / / / / Aluminium Pad · / / / / / / / / / / / / / / ,
/ / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / ,
,
/ / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / /
/ / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / ,
= °
m (17)
dC
- -0 . - (18)
m dz
(19)
+ [0.9. Om . U· p. de +.u 2 ]
39
where Jl = viscosity of the electrolyte
P = density of the electrolyte
U = average velocity of electrolyte reference to metal
d = interelectrode distance
de = critical dimension for Reynolds number
Z = distance normal to interface
Cf = friction factor
Cs = factor for skewness of velocity profile due to unequal
surface roughness
In the steady state, if the reaction between the two species is assumed to
be instantaneous then the rate of metal loss is equal to the rate at which
metal is transported to the reaction zone. The current efficiency is then
calculated as
(20)
40
An advantage of the model is that it is not restricted to the assumption of
any particular rate determining step. However, for known conditions in
certain Alcoa cells, the model without dissolution rate limits applied
predicted current efficiencies that were too low. Robl concluded that in
many instances it is likely that the dissolution rates for metal and/or gas
are the limiting factors.
Evans[38] assumed that the rate controlling step for the reoxidation
reaction is transport through a surface tension induced boundary layer at
the aluminium electrolyte interface. Using a relationship derived by
Levich[39] they correlated the rate of transport through the interface to the
turbulent kinetic energy. The kinetic energy was obtained from their k-e
model used to calculate fluid flow. This lead to the relationship
O.5
CE (22)
41
approximations for many of the unknown parameters reduces these
benefits.
Like the previous models Liliebuen[401 assumed that the reaction between
aluminium and carbon dioxide was solely responsible for current efficiency
losses. In addition he makes the following assumptions:
2) simple film theory is valid and transportation of the heat of reaction and
the reaction products is disregarded
The dissolution rates for aluminium and carbon dioxide taking into account
the reaction can be written
(24)
42
As in the previous models the steady state current efficiency is given by
r.CO 2 -
-
3 r
2" AI (26)
Attempts to calculate the mass transfer coefficient for gas bubbles, kC02 ,
by correlations proved unsuccessful due to lack of experimental data for
CO 2 diffusivity and bubble diameters. Lillebuen therefore assumes that
kco = kAI based upon the work of Aarebrot[42] with the ratio of gas bubble
2
area to aluminium area being estimated as 0.1.
k ·2£
where Sh = Sherwood number = _AI__
DA'
p·v·2f
Re = Reynolds number - ~--
11
11
Sc = Schmidt number
£ = interelectrode gap
p = density of electrolyte
oAI = diffusivity of aluminium in the electrolyte
43
v = interfacial velocity between electrolyte and metal pad
J.l = viscosity of electrolyte
0.03027 - 0.0006985· wt
(30)
1- 0.08989· wt
44
95
Temperature 965°C
Weight Ratio 1.35
if94 ~ \
Calcium Fluoride 5%
-
>
()
-
ACO 4cm
c: 93
.-
CD
.-
()
=
w
C 92
...CD...
::::J
0
91
90
1 2 3 4 5 6 7 8
Alumina Concentration (wt%)
Although this model shows reasonable results for the variation of current
efficiency with alumina concentration it does not include a sufficiently
accurate hydrodynamic model to show the effect of interelectrode gap.
Current efficiency has been measured for differing anode to cathode
distances in industrial cells by Rolseth et aH55] and at Anglesey Aluminium.
Both have found that the current efficiency varies linearly with ACO at
values above 3-4 cm but drops off rapidly below this. Lillebuen's model for
the dissolution of aluminium shows the rate of dissolution and hence the
current efficiency to be virtually independent of the anode cathode
distance. It may be that there are two mechanisms responsible for current
efficiency dependent upon the distance between the electrodes. At low
45
Reference Type of Cell Correlation Range
Gjerstad[46] Laboratory Cell Linearly increasing 0.5 - 12%
Belyaev[47] Laboratory Cell Minimum at 5% 1-11%
G rjotheim[48] Laboratory Cell Minimum at 4.5% 0-11%
Fi rsanova[49] Laboratory Cell Minimum at g% 4 - 10.5%
Bersimenko[5O] Laboratory Cell Linearly increasing 3 - 12%
Szeker[51] Laboratory Cell Minimum at 4.5% 2 - 12%
Abramov[52] Laboratory Cell Linearly increasing 5 - 12%
Schmitt[53] Side-break Cell Increasing 2 - 7%
Thonstad[54] Side-break Cell Increasing 2 - 7%
Poole[55] Side-break Cell Increasing 2 - 7%
Lillebuen[4O] Side-break Cell Increasing or const 0.5 - 6%
46
The experiments resulted in the following correlation
The volumetric flowrate of anode gases may be determined from the rate
of electrolysis at the anode. The number of moles of gas formed
1 10
(33)
2 Zo ·F
The electrolysis of oxide and fluoride are the only reactions at the anode
so
I
- (35)
4·F
Assuming that the gases behave ideally then according to the universal
gas law
R·T
- n .-- (36)
9 P
47
where v9 =volumetric flowrate m 3 S-1
R = universal gas constant J mol- 1 K-1
T = temperature of gas oK
P = Pressure of gas N m- 2
The velocity of the anode gases is the volumetric flowrate for unit surface
area of anode
- -.- R·T
-.-1 i·R·T
- (37)
4·F P AA 4·F·P
The speed of the stirrer has no direct counterpart in the cell but the
velocity at its tip may be equated to the maximum velocities found in the
metal pad.
Inserting the relationships and the values for the physical properties of
water into the equation gives
48
3.5.7 Model of Haupin
The rate of dissolution of metal through the metal boundary layer is given
The rate of convective transport of metal from the metal to carbon dioxide
boundary layers is given by
(40)
(41 )
49
Dissolved CO2
ANODE
Dissolved AI
C=o
L...
c / concentration
CATHODE
Deo2 _
--=- 1. 703
__ - 0.7630
_ _ _In _Reeo
_--=-2 Re o.s Sco.s (42)
Deo 2 1- 0.6272 In Re . cO 2.
eo2
50
3.6 Models for the Estimation of Frozen Bath Shape
Aluminium reduction cells are designed such that a layer of frozen bath
forms around the sidewalls of the cell. This protects the material of the
sidewalls from the corrosive effects of the electrolyte and the metal pad.
The frozen layer also helps to regulate the temperature of the cell by
freezing and melting. The freeze profile also affects the flow patterns
within the cell and the current distribution.
Models for the precise estimation of freeze shapes have been developed
by a number of workers[57,58,59]. These involve the division of the cell into
a number of triangular elements for which the temperature distribution can
be calculated by a finite difference method. The accuracy of the simulation
depends upon the number of lattice points used. This type of model is
therefore unsuitable for dynamic simulation due to the time required for
these calculations.
A simpler model for heat transfer through the side walls assumes
unidirectional heat flow through the walls. This model has been used by a
number of workers for the calculation of steady state freeze
thicknesses[60,61] and for the dynamic behaviour of freezes[2,5,62]. The heat
flow and temperature gradients are shown schematically in figure 3.4.
The fundamental equations for the flow of heat through the side walls are
qm - Am . hm. (t m- t 1 ) (43)
qs - As . hs . (t s - ta ) (45)
51
electrolyte
lj
/
freeze
~
carbon
sidewall
•
firebrick
t
steel
casing
" /
Aluminium
pad
Tb
Figure 3.4 Heat Flow and Temperature Gradient across Frozen Sidewalls
and suffix
m represents the melt ( bath or metal)
n represents a section of the wall or ledge with n=1 being
the freeze and n=m being the outer shell
s represents the outer shell
a represents the cell surroundings
The heat transfer at the inner and outer surfaces is assumed to be purely
by convection.
qm - Am ·h m·(t m-t 1 )
Most models have assumed that the temperature at the bath ledge
interface is the liquidus temperature of the bath. This is a reasonable
assumption if the melt can be assumed to be well mixed and the diffusion
of components away from the boundary layer fast.
52
bath. This indicates that there is diffusion of individual components from
the boundary layer to the melt.
aC j _ 0.. a2 c j u. aC j
at I ax 2
ax (47)
(a) The mass of the components moving between the freeze and the
boundary layer is calculated from the amount of the electrolyte which
freezes or melts.
The mass of the components moving between the bulk electrolyte and the
boundary layer is then calculated
(48)
53
20
- E
(,)
t/)
t/)
C1l
c:
18
~
.~
.c:
I-
C1l
N
C1l 16
...
C1l
11. f = 0.00001
14
o 1 2 3 4 5 6
1000
/ Bath temperature
--
o
...C1l
:::::J
~
980
940
o 1 2 3 4 5 6
Time I hr
It can be seen from figure 3.5 that as f is reduced the freeze becomes
more stable to changes in temperature. If it is to be assumed that f is very
little for the metal zone then the freeze profile in this zone will be stable.
Steady state models have found the ledge in the metal zone to be at least
as responsive to temperature as that in the bath zone if not more so. The
assumption that the diffusivity is less in the metal zone is also not borne
out by the measurements at Anglesey. These found the same ledge
composition in both zones.
Heat transfer coefficients will vary throughout the cell dependent upon
local conditions. An approximate heat transfer coefficient may be obtained
either by calculation using average conditions or by averaging the values
obtained for expected extremes of conditions.
Bennet and Meyers[63] give the following correlation for the heat transfer
coefficient of fluids, other than metals, flowing over a plate at Re < 5 x 105 .
where the physical dimension for evaluating the Nusselt and Reynolds
numbers is the interface dimension parallel to the fluid flow.
For aluminium flowing through a pipe Filimonov et aH64] gives the following
correlation
54
Number Symbol Equation
h·L
Nusselt Nu -
k
L3. p3 .g.~T
Grashof Gr 112
Cp'l1
Prandtl Pr
k
Rayleigh Ra Gr·Pr
L·v·p
Reynolds Re
11
The equivalent diameter for fluid flowing in a shallow stream is used in the
correlation for the metal. This is given by
4·H m ·R m
L = (51 )
Applying typical properties for the metal and electrolyte at 980°C Haupin
obtained values of hb = 370 W m- 2 K-1
and hm = 1500 W m- 2 K-1. The
result for hb compares well with values calculated from measured heat
fluxes but the calculated heat flux in the metal zone was 30-40% high.
This discrepancy was attributed to the additional thermal resistance of the
film of molten electrolyte always present at the metal/ledge interface.
The effective heat transfer coefficient is then given by
1
(52)
55
and anode-ledge distance on the heat transfer coefficient between the
bath and ledge. This resulted in the relationship
GF .AIJk2
hb(gas) - k1· ( (53)
ALD
This gives heat transfer coefficients in the region of 200-500 W m- 2 K-1 but
does not include the rotational effect of the bath due to electromagnetic
forces. By applying a correction factor to account for horizontal flow the
equation becomes
(54)
56
natural convection and natural convection plus superimposed velocity of
approximately 0.15 m S-1 is shown in figure 3.6. Whilst it probably
represents the most reliable data published, the accuracies for the heat
transfer coefficients is only +20%.
where a and m vary with orientation and size. Typically a varies between
0.2 and 1.3 and m is usually between 0 and 0.33.
For forced convection with laminar flow the correlation is usually of the
form
(56)
For the combined model with fluid velocities up to 0.15 m S-1 the range in
heat transfer coefficient is
This data shows that the heat transfer coefficient can vary by a factor of
two for various locations on a cell IS sidewall. Although these values are
higher than other estimates published in literature Welch claims them to
have been substantiated by independent unpublished data.
57
1000
, ...
- •• - •• Regression Line for Run 1
::l
Z 100 Regression Line for Run 4
Natural plus superimposed
forced convection
, .
~
..-
.. ... .
ci
-
z
a;
In
In
- - - - Flat Plate Correlation for Free
Convection
__ ----_e- -----_. - -.
• _.-D---a--- ---
..
....... . .
::l
Z
cau _.'
__ - - . -
0----- . . []- -- .0 _ - _ 0-
o
..J 10 •- ---
Natural convection
1
1.00E+04 1.00E+05 1.00E+06 1.00E+07 1.00E+08
Local Rayleigh No. Ra
The voltage drop across the cell may be broadly split into external and
internal voltages. The external voltages are those external to the
electrolyte and metal pad and include the ohmic drops across the buswork
of the cell, the anodes and the cathode. The internal voltage drops include
the ohmic drops across the electrolyte and the metal pad and the voltage
required for the electrochemical reactions in the bath.
(57)
The electrical resistance for the anode assembly, cathode assembly and
external buswork can all be simply measured on an operating cell.
Although the ohmic drop across the external buswork does not affect the
heat balance on the cell, it must be included to determine the power
efficiency of the cell.
The general equation for the resistance for a homogenous material with
uniform temperature is
58
1(.£
R - (58)
A
The resistance of the bath and metal pad may be approximated by the
above relationship when an effective area for current flow is used. The
effective area is an approximation to the true current distribution of the
cell. Accurate models for the calculation of current distribution and voltage
profiles are too complex for use in dynamic simulation. In models where
the temperature and composition of the electrolyte is assumed to be
uniform throughout, the current may be assumed to flow uniformly through
an effective area based upon the area of the electrode. Haupin[67]
correlated the effective anode cross sectional area with anode
environment by determining the isopotential planes under and along the
sides of each anode in a laboratory cell. This was done using an
aluminium coated tungsten reference electrode and a micrometer
positioning device. The effective area was correlated to anode enviroment
by fanning factors, defined as an additional length or width added to the
actual dimensions of the anode. The effective area of an anode is then
(59)
Haupin assumes that the effective area for the bath is equal to the total
effective area of the bank of n anodes and is given by
59
12
This value is used to calculate the current density at the anode and the
electrical resistance of the electrolyte.
(61 )
d-8 (62)
AB .1(
(63)
- 1( . ( 1- 1. 26fJ . A A
Haupin quotes this equation valid up to fc = 0.65. This model for the
electrolyte resistance is shown schematically in figure 3.8.
60
Rs - Re + Rba
Re - d - 0 o - ta
+ As K (1_£)1.5
As K
Rba - til
K (1 - 1.26 fc) AA
Anode
tA
000
d
Electrolyte
Aluminium Pad
1
Figure 3.8 Schematic Diagram of the Model for Resistance across the Bath
The accuracy of this model is dependent upon the relationships used for
fc' £, ta and 8. In calculations at the Trondheim course the following
relationships for fc and £ are used and ta and 8 are assumed constant.
100
0/0 of anode surface (64)
61
to the number density of the effective electrical charges in the melt. In
cryolite melts with large complex anions the current is carried mainly by
cations which may be assumed to move by the influence of an electric
field relative to the framework formed by the anions.
1\:mix (66)
The numerical value of 93.5 Q-1 cm 2 mol- 1 for the proportionality factor
was determined experimentally from the conductivity of pure cryolite by
extrapolation to its value at 1273 K. This factor is assumed constant when
the anion framework is not significantly affected by the additions.
The two other models assume that the Charge-carrying species do not
mutually influence each other. The conductivity is then the sum of the
contributions of the pure melts placed in two conductivity cells connected
in series or parallel. A schematic diagram of these two models is shown in
figure 3.9. For two melts with conductivities 1\:1 and 1\:2 and volumes V 1 and
V 2 in conductivity cells with cross sectional area A the resistivities R1 and
R2 are given
(67)
In the series model, the conductivity cells are connected in series and the
total resistivity, Rtotal is the sum of R1 and R2
+ (68)
where "-mix,s denotes the conductivity of the system for the two cells
1;""
62
Series model Parallel model
~x.I ·v~I
£.
K:mix,s (69)
Fellner states that this model can be shown to give a good description of
systems that exhibit nearly ideal thermodynamic behaviour.
The parallel model, following the procedure for the series model, gives the
following expression for the electrical conductivity "mix,p
The calculated data using all three of the models show good agreement
with experimentally measured values for the electrical conductivity at
various bath compositions used industrially. As shown in figure 3.10 the
simple model shows very good agreement for binary cryolite melts at
1273 K with the exception of the Na3AIF6 - LiF system containing more
than 5 wtO/o LiF. The experimental data for this system falls between the
predicted values of the series and parallel models approximating to the
former at low wtO/o LiF and to the latter at high wtO/o LiF.
Data calculated from the simple model for the tertiary system Na3AIF6 -
AI 2 0 3 - LiF is shown in figure 3.11. In this system AI 2 0 3 decreases the
conductivity while LiF increases it. This then provides a good test for the
suitability of the simple model to predict conductivity data. The data shows
good agreement with measurements by Matiasovsky et aH72 1.
63
3.4
3.2
3
,....
I:(,) 2.8
a:c..
2.6 MgF2 ----
2.4
2.2 AI203
2
0 5 10 15 20 25 30 35 40 45 50
Na3A1Fs + Additions, mol%
3 x o rrolO/o AI203 u
--
0
0.....;--
,....
Ie
2.8 F
x 10 rrol% AI203
0 -
(,) 0
a "II
2.6 u
~
0
x = 20 rrol% AI203
2.4 0
2.2
2 I I
o 10 20 30 40
NCia AIFs + x AI 2 0 3 mol % LiF
~G¢'
(71 )
z·F
R· T ( rra~ )
+ --In (72)
z·F
There is also some excess voltage above the reversible range required for
the overall reaction. This is the reaction overpotential and may be
considered to consist of an anodic part and a cathodic part.
[4]
(73)
64
Haupin[74] assumes the activities of all reactants and products, except
alumina, to be unity. This is valid if these components may be assumed to
exist in their standard states, that is that the aluminium and carbon are
essentially pure condensed phases and that nearly pure carbon dioxide
contacts the anode at approximately 1 atm pressure. Alumina in solution
will have unit activity only at saturation and Haupin uses the following
relationship for alumina activity
(74)
65
Under different circumstances the concentration and reaction
overpotentials have greatly different values and the reaction can become
concentration or surface limited
with coefficient a ranging from 0.19 to 0.21 and f3 ranging from 0.23 to
0.25. Stirring lowered the overpotential by 40-50% which Haupin[30] takes
to indicate a concentration overvoltage. Thonstad and Rolseth found the
cathodic overpotential in industrial cells to be between the values for
unstirred and vigorously stirred laboratory cells at about 0.1 V. Alcoa data
for industrial cells[77] resulted in the following empirical relationship valid
for current densities above 0.1 A cm- 2 .
66
References
10 Maeda, H., Matsui, S., and Era, A., "Measurement of Dissolution Rate
of Alumina in Cryolite Melt", AIME Light Metals, 1985 pp763-780
67
11 Asbj0rnsen, O. A., Andersen, J. A., "Kinetics and Transport Processes
in the Dissolution of Aluminium Oxide in Cryolite Melts", AIME Light
Metals 1977, pp137-152
15 Hove, S.J. and Kvande, H., "Centre-Break Alumina Feeding and Sludge
Control of Prebaked Cel/s", AIME Light Metals, 1982, pp513
68
21 Haupin, W. E., "Mathematical Model of Fluoride Evolution from Hall-
Heroult Cells", The 4th Int. Course on Process Metallurgy of Aluminium,
Trondheim, 1985
24 Newitt, Dombrowski and Knellman, Trans. Inst. Chem. Engr., vol 32,
1954,p244
31 Qiu, Z. and Fan, L., "The Rate Determining Step of Metal Loss in
Cryolite-Alumina Melts", AIME Light Metals, 1984, pp789
69
32 Poole, R.T., Etheridge, C., "Aluminium Reduction Cell Variables and
Operations in Relation to Current Efficiency", AIME Light Metals, 1975,
pp163-183
35 Bratland, D., Jose del Campo, J. and Cho, K., "Current Efficiency
Measurements in Laboratory Aluminium Cells. Influence of Bath
Volume, Anode Size and Anode Quality", AIME Light Metals, 1981,
pp281-308
36 Berge, B., Grjotheim, K., Krohn, C., Naeumann, R. and T0rklep, K.,
"The Influence of Operating Parameters on the Current Efficiency in the
Aluminium Reduction Cells", AIME Light Metals, vol 1, 1976, pp23
38 Evans, J.W., Zundelevich, Y. and Sharma, D., Met. Trans, vol 12B,
1981
40 Lillebuen, B., Ytterdahl, S. A., Huglen, R., and Paulsen, K. A., "Current
Efficiency and Back Reaction in Aluminium Electrolysis", Electrochimica
Acta, vol 25, No, 2, February 1980, p131
70
42 Aarebrot, E. et ai, Metall. vol 32, no 41, 1978
46 Belyaev et ai, Izv, Vyssh. Ucheb. Zav., Tsvetn. Met., vol 4, no 3, 1961,
p67
50 Szeker, Acta Tech. Acad. Sci. Hung., vol 10, 1955, p91
71
55 Rolseth, S. and Thonstad, J., "On the Mechanism of the Reoxidation
Reaction in Aluminum Electrolysis", AIME Light Metals, 1981, pp289
57 Arita, Y., Urata., N. and Ikeuchi, H., "Estimation of Frozen Bath Shape
in an Aluminum Reduction Cell by Computer Simulation", AIME Light
Metals, 1978, pp59
62 Paulsen, K.A., Huglen, R., Andersen, J.A., Lillebuen, B., and Ytterdahl,
S.A., "Variation of Side Lining Temperature, Anode Position and
CurrenWoltage Load in Aluminium Reduction Cells", AIME Light
Metals, 1980, pp325-339
72
65 Solheim, A. and Thonstad, J., "Heat Transfer Coefficients between Bath
and Side Ledge in Aluminum Cells. Model Experiments", AIME Light
Metals, 1983, pp425
69 Sides, P.J. and Tobias, C.W., J. Electrochem. Soc., vol 129, 1982,
p2715
73
75 Dewing, E.W., "Liquidus Curves for Aluminum Cell Electrolyte v.
Representation by Regression Equations", J. Electrochem. Soc., 117,
1970, pp780-781
74
4.0 Development of a Dynamic Model
A conceptual model of the aluminium reduction process is shown in figure
4.1. The interconnections between the functional blocks represent the flow
of mass, energy and information. Each block represents a set of equations
describing the relationships between the inputs and outputs of that block.
These relationships take the form of differential and algebraic expressions
derived from analysis of published work, plant data and mathematical
modelling techniques.
For a stirred tank reactor with a number of streams entering and leaving
the vessel this can be written
~dCj - fC j dV (1)
V.L..J- -
j=1 dt j=1 dt
The mass balance for an electrochemical cell differs primarily in that there
are both chemical and electrochemical reactions to be considered. The
reaction term in equation (1) therefore becomes
75
LI.~--------~~~~=::":'--------~
Bath temperature
Freeze melt rate Freeze
r-1 Sludge
Formation
Anode
change
AI product
Breaker
action 1: ) -\ Dissolution I Net
-( 1: 1 -I
Electrolyte
mass
balances
I
r---1 Heat
losses I
.....f - - -
Dissolution
Rate AI 2 0 3
Stochastic
disturbance AI 2 0 3 concentration
D1 concentration estimated
Anode beam
position 1: )
t -
Icond~ctivity
Electrolyte ~
rop 1:
Cell R =V / I
volts Energy = 12 R
Pot heat
balance Bath
temperature
Stochastic
disturbance
D2
Electrolysis Electrolysis
charge heat
balance balance
Line current
In the aluminium smelter, the feed streams into the cell are the feed rates
of alumina, cryolite and aluminium fluoride into the electrolyte. These
flowrates are not continuous and the cell can be treated as a semi batch
reactor where the flowrates, Ok, and the feed concentrations, Cfk,j' vary
discontinuously with time. The main stream out of the reactor, 0 1, is the
periodic tapping of aluminium from the metal pad.
In addition to the heat and mass balances for the reactor, the model will
also include many algebraic equations describing the relationships
between parameters in the model, for example, the dependence of density
on the composition and temperature of the electrolyte. The remainder of
this chapter will describe in detail the modelling of each block in the
conceptual model and the equations used therein. This represents the
physico-chemical model of the aluminium cell. The model of the control
system to be applied in the operation of the cell will be considered in the
following chapter.
76
4.2.2 Anode beam position
The anode beam position is controlled in the same way as the crust
breaker by pneumatic pumps. Input is as a number of seconds activating
the raise or lower pumps. The anode beam position governs the anode to
cathode distance (ACO) and hence the voltage drop across the
electrolyte. It also affects the mass balance as it is one of the parameters
that affects the reoxidation of aluminium.
All cells in series draw the same current. Anode effects in other cells in a
line can cause current fluctuations in individual cells. Current fluctuations
are generally small and in this model the current will be assumed to be
constant throughout the simulation. As well as determining the power
input to the cell by the ohmic drop across the cell, the current determines
the rate of the electrochemical reactions.
77
4.3 Alumina Dissolution
It has been shown in the previous chapter that the dissolution of alumina
is dependent upon the method of feeding. This is due to the formation of
agglomerates and sludge in the electrolyte. Laboratory experiments have
differentiated between two extremes of dissolution, that from finely
dispersed particles of alumina and that from alumina agglomerates and
sludge. In an industrial cell dissolution will be by a combination of effects.
This is reflected in the models used by Blake[1] and Entnen2] in which
alumina is assumed to exist as dissolved ions, a suspension of small
particles and sludge. The models differ in that Blake assumes a sequential
mechanism from sludge to suspension to solution whilst Entner assumes
a parallel mechanism from sludge to solution and suspension to solution.
Blake assumes a proportion of alumina to enter each state on addition to
the electrolyte whilst Entner assumes a settling of suspended alumina to
sludge.
The sequential model for dissolution must assume the stripping of alumina
from sludge into suspension by the movement of electrolyte across the
surface of the sludge. This is more applicable to sludge at the side walls
than to large agglomerates floating on the metal pad and is therefore a
more realistic model for side break cells in which large amounts of sludge
are present at the sidewalls than in half break and point feed cells.
A model for the settling of alumina in the cell will produce an alumina
concentration distribution but will increase the number of calculations
required for solution. The accuracy of the model will depend upon the
number of zones used and the assumptions for the settling of alumina.
This in turn will depend upon the geometry and hydrodynamics of
individual cells. If a uniform distribution of alumina is assumed then this
model approximates to that used by Blake.
78
detects the presence of the oxygen atom by precipitation of alumina and
so empirical relationships do not depend upon the form that dissolved
alumina takes.
(3)
79
T = reaction temperature oK
*
c = saturation concentration of alumina in melt g cm-3
c = concentration of alumina in melt g cm-3
Values for the parameters in the above equation have been determined by
a number of sources and have been summarised by Asbj0rnsen[51. The
data from Kachanovskaya will be used in the simulation, the values of the
parameters being ko = 0.00124 cm S-1, To = 1383 oK and
E = 18.153 kcal mol- 1 . For simplicity, alumina dissolving from both sludge
and suspension will be assumed to obey the same rate expressions with
the same values for kinetic parameters. The dissolution rates will be
different due to differing areas of contact between the electrolyte and the
alumina as sludge and in suspension
where
(6)
The contact area for alumina in suspension will be less than the specific
area of the alumina as measured in a typical analysis due to the presence
of alumina/bath agglomerates. Sludge will tend to form piles on the bottom
of the cell and the contact area will be further reduced. The contact areas
will be based upon the specific area for the particular alumina used but will
be fitted to produce the desired dissolution curve. This will allow the
simulation of the effects of different types of alumina and different cell
geometries allowing the formation of different sized agglomerates.
Alumina fed to the bath contains an amount of Na20 which reacts with
AIF3 in the electrolyte
80
Although the proportion is small, typically 0.5°10 Na20 by weight, over a
period of time it has a noticeable effect on the weight ratio NaF/AIF3 in the
electrolyte. For simplicity, all alumina in suspension and sludge will be
assumed to contain Na20 in the same proportion as the alumina feed. The
dissolution rate of Na20 into solution will be assumed to be proportional to
that of alumina and given by
In many cases the alumina is also blended with material captured in the
l
exhaust scrubbers. This Icatch consists predominately of NaAIF4 from
particulate fluorides and AIF3 from the chemisorbtion of gaseous fluorides
onto A1 20 3. It may therefore be considered to be a mixture of NaF and
AIF3 and the dissolution rates will be determined in the same way as that
for Na20.
Again the actual amounts of each compound in the feed are small but
significant with the total fluoride content of a blend being typically 0.5-1.0
wt°lo.
81
Bath
Volatilisation
Alumina Mass ~
Dissolution Balances
Back
Reaction
[2]
82
F
G
= 2310000
%CE
_ 955000. R .[
P b
( 8444)]
exp 7.553 - T .
b
0/oH 20 in AI 20 3 + %H in anOdeS]~ .
[ 37.44 21.5
(10)
%CE
.- moles AI S-1
100 3·F
%CE 27
--.-.
100 3· F 1x 1 0 6
83
dMF(Q)
dt
%H 2 0 in AI 2 0 3 + %H in anOdeS]~ .
[ 37.44 21..5
dMO(g) 1 dMF(g) Wo
- _e . _
gO s-1 (13)
dt 2 dt WF
dMO(g)
where dt = rate of addition of oxide to bath by gaseous
evolution of fluoride
Wo = atomic weight of oxygen = 16 g mol- 1
WF = atomic weight of fluorine = 19 g mol- 1
(14)
84
(15)
and
-_.& (18)
4·F Pb
dM vo'
dt
g S-1 (19)
85
The vapour pressure of molten cryolite with additives has been extensively
studied and there is good agreement among recent investigators. The
data has been combined by Haupin to produce the following formula
which will be used in this work
kPa (20)
where
The total vapour pressure above the electrolyte is also the sum of the
partial pressures of the monomer and dimer.
Pv = Pm + Pd (21 )
The equilibrium constant for dimerization has been correlated from mass
spectrographic measurements[7] and is given by Haupin as
= exp (
-21414
T +15.6
J (22)
86
_ [-1 + (1+4·K d.PV)O.5]
Pm (23)
2·K d
introducing equations (21) and (23) into (19) and rearranging gives
dM VO'
dt (24)
dMAI~~, dMFvo' W A1
- g AI S-1 (27)
dt dt WNaAIF4
dMO G
where = rate of addition of 0 2 - ions to electrolyte by the
dt
evolution of HF gas g S-1
dMNaFvO' f
-----'-~ = rate 0 volatilisation of NaF g S-1
dt
3+
dMAl vo,
= rate of volatilisation of A13+ ions g S-1
dt
The pot gases from all cells are collected by an exhaust system and
passed into a baghouse where particulate fluorides are filtered and
gaseous fluorides are chemisorbed onto alumina. The product, termed
fines, is added to the alumina fed to the cells.
87
4.4.4 Aluminium dissolution and back reaction
All the models proposed for the rate of the back reaction of aluminium
require the estimation of a number of parameters. In general, these
parameters will differ from cell to cell and the ability to tune the model to
the real cell will depend upon the number of parameters that must be
estimated and the complexity of the model. The complexity of the model
also determines the speed of calculation in dynamic simulation.
The models of Robl[8] and Evans[9] may be capable of estimating the rate
of reaction in different sections of the cell but the difficulty in estimating the
parameters required reduces the accuracy under dynamic conditions. The
complexity of these models also makes them unsuitable for dynamic
simulation.
The assumption that the bulk of the electrolyte is well mixed leads to a
uniform distribution of all species. Under these conditions the mixed rate
models for the back reaction of aluminium proposed by Haupin and Robl
reduce to the dissolution rate controlled model of Liliebuen[10]. This model
has shown reasonable agreement with industrial cells when used for the
calculation of steady state current efficiencies and shall be used in this
work.
[3]
(28)
88
A = interfacial area of species and electrolyte
k = dissolution rate constant
C = concentration of species
C* = saturation concentration of species
suffix AI refers to aluminium
CO 2 refers to carbon dioxide
The models presented for the dissolution of aluminium based upon the
hydrodynamics of the system (Robl, Evans) are complex and require the
estimation of a large number of unknown quantities. The advantage of
these models is offset by the lack of any data relating parameters in the
equation to operating conditions.
89
- 2.147x10- ( T
11 i. T )1.42
.(1+313.3V 2.5 ).
(
~::
J
O 33
.
·D~~7 (30)
The dissolution rate constant for carbon dioxide has been proposed by
Haupin as
Haupin does not, however, state how this relationship was determined
and so no comment may be made upon its validity. There is as yet no
data relating bubble diameters to operating conditions or the geometry of
the cell. As the current efficiency is determined from the dissolution rate of
carbon dioxide these relationships must be known in order to simulate the
effect of the parameters. Until such data is available Haupin's relationship
has limited use in dynamic simulation.
90
Lillebuen proposed that the dissolution rate constants for aluminium and
carbon dioxide were equal based upon the assumption that the interfacial
area for carbon dioxide dissolution is ten times that for aluminium and the
ratio of kc02 'Ac02 to kAI'AAI is 0.1. Vasiliadis, however, found an average
value of this ratio to be 0.06 and hence kAI / kco = 0.6.
2
The interfacial area for carbon dioxide dissolution has since been related
to the concentration of alumina in the melt but unfortunately to no other
parameters such as temperature. Using Haupin's relationship however
it can be seen that the ratio of areas will vary with alumina concentration
and anode current as the aluminium pad area remains approximately
constant. The ratio of mass transfer coefficients will therefore depend
upon the average operating conditions. In this work, due to the lack of
carbon dioxide bubble data, the dissolution rate constant for carbon
dioxide will be determined
(33)
91
this reaction is no longer possible. If this is the only mechanism for current
efficiency losses then during anode effect an efficiency of 100% would be
obtained. General operating experience on central point feeding cells has
indicated high current efficiency at low alumina concentrations ie close to
anode effect although this could be due solely to the reduction in bubble
surface area at low alumina concentrations as shown by Haupin1s
relationship and not the presence of an alternative reaction. The current
efficiency during anode effect in industrial cells cannot be determined by
gas analysis and the Pearson and Waddington equation because the gas
contains species other than CO and C0 2 J11]
(34)
(35)
92
Due to the lack of data available for the dissolution of CF 4 into molten
electrolyte the dissolution rate constant and the saturation concentration
of CF4 will be taken to be the same as the values calculated for CO
2
dissolution. As fluorocarbons cause dewetting of the anode surface and
coalesce to form a single large bubble, the interfacial area between
electrolyte and CF 4 gas will be approximately equal to the surface area of
the anode covered by the gas. Defining a fractional coverage of the anode
surface by CF 4 as <I> gives the interfacial area for CF4 dissolution as
(36)
In the build up to anode effect both oxygen and fluoride ions are oxidised
and during this period aluminium will be consumed by both the above
reactions. The net rate of consumption of dissolved aluminium is then
In- + ne ~ J (38)
at an electrode is
93
(39)
- CEAlo-- (40)
ZAl oF
The overall current efficiency is defined as the ratio of the actual rate of
production of aluminium to the theoretical rate and differs from that
determined by the Pearson Waddington equation in that it includes effects
other than chemical reaction of aluminium at the anodeo Here, the current
efficiency losses are assumed to be due to back reaction and cathodic
losses. The net rate of aluminium production is therefore
(41 )
CE =
CE _ rback (42)
- AI r.
t
I (43)
-
ZAl oF 3F
94
4.4.8 Mass balances
- R.J (44)
The mass balances for individual species in the aluminium reduction cell
are developed in the following sections.
95
than 2.0 wt% total, but have a marked effect upon the operation of the cell
and therefore must be included. For simplicity the rates of dissolution of
each compound are considered to be proportional to the rate of
dissolution of alumina. It is therefore assumed that all sludge and
suspension has the same chemical composition as the feed and the mass
balances will be performed as on a homogenous material.
dMsuspn (45)
dt
dMs,udge
- s f' [dM
dt
feed
- dM
dt
ac
- rdiss(sludge) ]
(46)
dt
Alumina additions by the crust breakers and anode change are not
continuous and it will be assumed that any alumina introduced to the
electrolyte will instantly enter the sludge and suspension states. The
addition rates may therefore be represented by a series of impulses with
varying period and height. The period represents the frequency of
operation of the crust breakers or anode change and will be set by the
control system. The height represents the amount added to the cell and is
subject to a random proportionality factor to simulate the uncertainty
96
between the actual mass added and the intended mass. The mean and
variance of the amount of alumina added to the cell may be measured by
sampling a working cell. The random variable may then be generated from
a distribution function based upon these variables. An example of the form
the rates of addition might take is shown in figure 4.3.
dMAI 3+ bag
-_ [ 2 . ( rdiss(sludge) + rdisS(Suspn) ) .
+ r AIF3 dlss
+ dM
dt ]. WAI
dt
dMAI~~1 dMAI~+ dMAI!~ _ CE. .W + ~ (47)
- d AI F AI back
dt dt t ZAI'
97
-25 -20 -15 -10 -5 o 5 10 15 20 25
Percentage Overfeed
dM0 2 -
= [3. (rdiSS(SIUdge) + rdiSS(SUSPn)) + rNa20diss ]. Wo
dt
(48)
dM0 2f -
where ------.:.- = rate of transfer of 0 2- ions into frozen
dt
sidewalls, +ve for freezing
dM0 2 -
_------=a:.::--c = rate of transfer of 0 2- ions into frozen
dt
layer on cold anode, +ve for freezing
98
4.4. 12 Aluminium metal in electrolyte
Aluminium dissolves from the metal into the electrolyte where it can react
with dissolved CO 2 or CF 4 , The dissolved aluminium is assumed to be non
ionic. The mass balance is then
Mathematically this is
dMAI
- rAldiss - rback (49)
dt
Mathematically this is
99
dMAl tap
where dt = rate of metal tapped from the cell
[5]
At the operating temperature and composition of the cell this reaction will
be fast and the rate determining step will be assumed to be the dissolution
of Na20. Two moles of NaF will therefore be formed for every mole of
Na20 dissolved.
100
dMNaF dMNaFf dMNaF dMNaF
dt
= ~ ,
NaFdlss
-
dt -
dt ac
- dt vol
+ ~Na 2Od'ISS ' 2 , WNaF (51 )
dMNaFf
where dt = rate of transfer of NaF into frozen
dMNaFac
dt = rate of transfer of NaF into frozen layer
dMNaFVO'
------=.=- = rate of volatilisation of NaF
dt
As with sodium ions, the calcium ions in the electrolyte are assumed to be
associated only with calcium fluoride. In practise the alumina fed to the
cell contains a small percentage of CaD which reacts with AIF3 in the melt
to form CaF2. The amount of CaD in the feed is typically a tenth of the
amount of Na20 and as such may be considered negligible. Similarly, the
mass of calcium co-deposited at the cathode is a fraction of the mass of
sodium and will also be assumed negligible. In this model the calcium
fluoride content of the molten electrolyte will therefore only be affected by
the rates of freezing at the sidewalls and on cold anodes.
dMCaF
where - - - - ' -2f' = rate of transfer of CaF2 into frozen
dt
sidewalls, +ve for freezing
101
dMCaF2ac
dt ' = rate of transfer of CaF2 into frozen layer
on cold anode
Since this work assumes that the only oxide containing species dissolved
in the melt is alumina, (other compounds, eg Na20, react to form alumina)
it follows that the mass of alumina is proportional to the mass of 0 2- ions.
From the chemical formula of alumina, A1 20 3, it can be seen that each
mole of oxide is equivalent to one third moles of alumina. The mass
balance on alumina in solution is therefore
1 dM0 2 - WAI203
- -. . (53)
3 dt Wo
W AI
2
0
3
= molecular weight of alumina = 102 g mol- 1
This work assumes that cryolite is a mixture of NaF and AIF3 in the mole
ratio 3: 1. Therefore, the only species in the electrolyte that contain the A!3+
ion are AIF3 and AI 20 3 with mole ratios 1 mole A!3+ per mole AIF3 and 2
102
moles A13+ per mole AI 2 0 3 respectively. The A13+ ions associated with
alumina may be determined from the oxide content of the melt and any
remaining A13+ ions must therefore be associated with AIF . The
3
aluminium fluoride mass balance is
dMAIF3
dt (54)
The rate of change of the total mass of the bath is the summation of the
rates of change of mass of the individual species that make up the
electrolyte. That is
103
1 [2.0156 - 0.0207· %A1 2 0 3 - 0.005· %caF2 ]
p
= 1.06· exp + 0.4349. Rb _ 2068.4 (56)
T
f
R = (57)
A·p
The area for current flow through the electrolyte will be affected by the
design of the electrochemical cell. A good approximation for the design of
cell used at Anglesey Aluminium has been found to be the assumption of
a two inch shadow around the block of anodes. The effective anode area
is then
(58)
This relationship assumes that all the anodes are drawing the same
current. This is not the case however when an anode has been removed
from the bath or a cold anode has been introduced to the bath and is
insulated by a layer of frozen electrolyte. The heat balance across the
anodes and the algorithms developed in a later section for the modelling
of an anode change, calculate a fraction of the surface area of each
anode that is effectively available for current flow. In this work, the
effective anode area will be further modified by subtracting the area of
each anode that is effectively insulated. The calculation of the effective
anode area then becomes
104
where fj = fraction of anode able to conduct current
This modification maintains the two inch shadow around each anode.
I+~
A + ne ( ) B (60)
f-I
by the equation
zF zF
-a-ll (1-a)-1l
- 1 ' e RT _ e RT
0 (61 )
105
~~~~~~~~~~
As .t------ Bs
reaction via
intermediates
mass mass
transport transport
Ab Bb
Figure 4.4 Schematic Representation of Reactions at an Electrode
The Butler Volmer equation can be rewritten for the individual rates
zF
-a-1J
- ioe RT
(62)
and
(63)
Usually the conditions at the electrode are well removed from equilibrium
and i+ or _ » Lor +. The overall rate then becomes
(64)
Note
The Tafel relationship for an electrochemical reaction is analogous
to the Arrhenius rate expression for a first order chemical reaction
where the rate of reaction is given
E
r - A·e RT
106
'-
f\
f\
I
W
C!:l. '-c::+
I ·2
.....
Q)
~
~ .~
Q)
W
II
........
..
-.-ctS -
Q
~
~
...... 't--
0
c: E
Q)
...... ~
0)
0 .~
a. -0
.~
.....
ct!
E
Q)
..c::
(.)
C/)
LO
~
@
:J
.~
LL
I
•
'AI!SUaa IUaJJ n8
The rate of reaction of A is given by Faraday's law as
z·F (65)
(66)
At steady state the rate of transport of A to the surface will be equal to its
rate of reaction at the electrode
dNA
-- (67)
dt
(68)
107
From (65) and (66) the following simultaneous equations can be obtained ,
(69)
iA = z F k [Als e-j3E
(70)
- [AJ IA iA
----'-'-e -fiE
b z.F.k L (71 )
z·F·k
and
1 1 (72)
+
z·F·k·e -fiE
z·F·k L
Although the model has been derived for the simple reaction
A+ze ~ B [6]
it applies to any reaction that exhibits Tafel type behaviour, and can be
shown to be valid for many more complex reactions[ 121.
0 2- ~ 0 + 2e [7]
108
Most workers have fitted the experimental results to Tafel type behaviour
although some have included a concentration overvoltage effect to
accommodate the rise in voltage at low concentrations of alumina. Some
of the relationships derived will be considered here and the parameters to
be used in this work will be determined.
where
R·T I
11 SA ---·In-
a·z·F io
Io - Y + 8·c b,A1 0 2 3
_ R.T.lni
z·F ie- i
_ [5.5+0.018(T -1323)].A-o.1,[c~:!'203 -0.4]
_ £. (C~·5 - Jl) for a given temperature
cb,A1
20 3
=concentration of alumina in the electrolyte
A = anode area
I Ie (74)
11 = A·ln-. + B·ln--
ie -i
Io
109
Rewriting equation (74) in terms of the overpotential , 1] -
-
E - Etio'1', gives
.
I =
1 1 (75)
+
z·F·k L
rearranging gives
17 - 1 1 [ z·F·k·i
-·In-· L
]
b a z . F . kL . c b - i (76)
1 [i', i ]
17 - b1 ·In a' icc - i
1 .,
- -.In~ + 2. lni 1 1
+ -·In--
b a b b ic- i
- a + fi·lni + fi.ln-.1_.
Ie -I
Comparison of equations (74) and (77) shows that although the two
models are similar, they differ in several important respects and it is not
possible to match the coefficients in the models and so determine the
parameters required for this work.
The differences between the Trondheim model and the theoretical model
are probably due to the derivation of the Trondheim model. The
Trondheim model is highly empirical and has been derived over a number
of years by amalgamating several empirical relationships as research has
made them available. Whilst this approach can lead to good results over
the specific range of parameters investigated, it may no longer be possible
to relate the derived relationship directly to the underlying theory, as in this
case, bringing the applicability of the relationship under a wider range of
conditions into question. For this reason the theoretical relationship will be
used in this work with the required Tafel coefficients and mass transfer
coefficients being determined from the available empirical relationships for
polarisation voltage.
110
4.6. 1 Determination of Tafel constants and mass transfer coefficient
Tafel coefficients for a wide variety of carbon anodes with varying porosity
have been determined by Zuca et aH13] and figure 4.6 shows the variation
of overvoltage with current density for different carbon types at 965 ce.
The Tafel plots determined by Zuca are also shown on the graph. The
Tafel coefficients for the equation 11 = a + b log i and the range of current
density over which Tafel behaviour was assumed are shown in table 4.1.
The cryolite composition used was similar to that used in industrial cells
with 5°1o alumina, 4°1o calcium fluoride and a cryolite ratio of 1.35.
The experimental data in figure 4.6 clearly shows the effect of mass
transfer at higher current densities with the rise in potential as the reaction
becomes mass transfer limited. It is evident from the graph that the range
over which the Tafel coefficients are calculated includes some effect due
to mass transfer and so the coefficients cannot be used directly in
equation 72.
Assuming that anode effect is initiated when the oxide reaction becomes
mass transfer limited then at this point the partial current density for the
reaction will be given by
I = z·F·k L ·c b (78)
(79)
111
Anode Porosity Tafel intercept Tafel slope Tafel range
0/0 V V/decade A cm- 2
0.7
>
co 0.6
...c
.-
~ 0.5
~o
I o vitreous carbon
a.
a..
Q) • baked carbon
> 0.4
0 o porosity 8.90/0
0.3 x porosity 12.82%
0.1
0
0.01 0.1 1 10
Current Density Alcrrt2
Having obtained a value for the mass transfer coefficient, the required
Tafel coefficients, a and /3, can be obtained by solving equation (77) at
two different current densities, with the overpotential being calculated from
the equations derived by Zuca.
For a baked carbon anode with porosity of 14.73% , Zuca derives the Tafel
equation
over the range i=0.05 to i=0.5 A cm- 2 . The anodic overpotential at these
limits is then
11 - a + /3 i .J
.In(-. (81 )
Ie -I
Substituting the values for the composition at which the experiment was
carried out and the calculated value of kL into the equation gives ic =
2.3163 A cm-2 . Substituting this value and the values obtained for 11 at 0.5
and 0.05 A cm-2 into equation (80) gives the two simultaneous equations
112
Solving these gives a = 0.653 and f3 = 0.080 and
11 - 0.653 + 0.080.ln( I
752.7,c b -i
J (82)
Equations (77) and (82) are plotted in figure 4.7 for an electrolyte
composition of 5 wtOlo alumina, 4 wt°/o calcium fluoride, cryolite weight ratio
of 1.35 and at 965°C. The Tafel plot derived by Zuca for an anode
porosity of 14.73°/0, equation (80), is also shown.
As would be expected, the equations (80) and (82) give very close results
for the range of current density over which Tafel behaviour was assumed.
As the current density increases, the rate of reaction becomes mass
transfer controlled and the potential rises more rapidly.
There is a very great similarity between the shapes of equations (74) and
(82) although the equation combining 1JSA and 1JCA' (74), produces an
overpotential some 0.15 V less than that measured by Zuca. This could be
due to the reference used in calculating the overpotential.
R.T.(1.375-0.125.CR)ln( ic J (83)
1.5·F 0.257
The combined results of equations (74) and (83) are shown in figure 4.7
with the results from equation (82). The expression for the cathodic
overpotential is only valid at cathodic current densities above 0.1 A cm- 2 .
Inclusion of the cathodic overpotential causes the two relationships to
converge at higher current densities but results in a greater divergence in
the Tafel region of the graphs. Information on the reference potential for
the derivation of equation (83) was not available and it is possible that this
could conflict with the reference used by Zuca or those used in the
113
Wt% Alumina 5%
Temperature 975°C
1.2
Ratio 1.35
Calcium Fluoride 4% 1
->---
Q)
o Tafel only
X Tafel + mass transfer
t:, Haupin, sa + ca 0.8
-...o
C)
co o Haupin, sa + ca + cc
>
c
...
o
CO
0.4
-en
~
~X
-o
CO 0.2
a..
0.1 10
-0.2
In this work, equation (82) will be used as the model for the anodic
overpotential of the oxide reaction. The effect of current density upon the
polarisation voltage at different alumina concentrations is plotted in figure
4.9. This graph clearly shows that the onset of anode effect occurs at
higher alumina concentrations as the current density is increased as is
experienced in practice.
m m
- LR j=1
j (chemical) + LR
j=1
j (electrochemical)
(84)
For a normal chemical reaction the kinetic rate expression would take the
form
_ k. [ ]a. [ ]b (85)
R j (chemical)
114
1.2
1.1
Current Density 0.7 A/cm2
1 Temperature 965°C
-
~
Q)
0.9 o
Ratio
Calcium Fluoride
Haupin, sa+ca
1.35
4%
C) 6. Haupin, sa+ca+cc
~ 0.8 o Tafel + mass transfer
0
>
c 0.7
-
0
ca
~
~
ca
0.6
~ 0.5
0.4
0.3
0.2
1 2 3 4 5 6 7 8 9 10
Alumina Concentration (wt°jo)
1.1
Temperature 965°C
1 Ratio 1.35
-
~ 0.9
(1) o
Calcium Fluoride
o id = 0.5 A/cm2
id = 0.7 A/cm2
4%
C)
X id = 0.9 A/cm2
~ 0.8
0 I:::. id = 1.1 A/cm2
>
c 0.7
0
+:
co
en 0.6
";:
co
~~------X.-- ___- 5<-
~ 0.5
0.4
0.3
0.2
0 1 2 3 4 5 6 7 8 9 10
Alumina Concentration (wt%)
I·J
R j (electrochemical) - (86)
z·F
The total current is the sum of the partial currents for all electrochemical
reactions in the reactor and a current balance must be done on the
electrode su rface
(87)
m
_ ~
LJ a J. . x J.. p~at
J
(88)
j=1
115
Xj = mole fraction of species j in the gas
PJ~at = sat urate d vapour pressure of species j
= A.exp(- B ) (89)
C+T
In the electrochemical reactor the partial current density for each species
at an electrode may be related to activation (electrode potential) and mass
transport by the polarisation curve. The equation for this curve has been
derived in the previous section and the relationship for a reaction
exhibiting Tafel behaviour is described by equation (72).
116
The solution of these simultaneous equations can become complex for
multiple reactions but can be simplified for the aluminium reduction cell by
making the following assumptions:
(91 )
and can be easily calculated. The partial current for the anode effect
reaction, the oxidation of fluoride ions, is then the remainder of the
total current through the cell
(92)
117
As has been stated previously, the anode effect in an aluminium reduction
cell results from a change in the electrochemical reactions occurring at the
anode. During normal operation, electrolysis depletes the concentration of
oxygen ions in the electrolyte until the concentration reaches a level
where the reaction becomes mass transfer limited and the reaction is
unable to proceed at the rate imposed by the current. As the cells are
operated at a constant current alternative ions will be oxidised at the
anode to maintain the overall rate of reaction. In the cryolite melt it is
fluorine ions that are discharged.
The fluorine discharged reacts with the carbon anode to form various
carbon fluorides, ie CF4 , C2FS etc. Unlike the carbon dioxide formed under
normal operating conditions, fluorine is very effective at preventing the
wetting of the anode surface by the electrolyte. The fluorocarbon bubbles
'stick' to the anode surface and coalesce into an insulating gas film. This
film spreads rapidly over the anode surface resulting in a high ohmic
resistance across the anode-electrolyte interface.
The fluorocarbon gas film is tenacious, and whilst it covers the surface,
normal anode reactions cannot be resumed, hence the requirement of
mechanical agitation of the anodes in anode effect quenching routines.
118
..5!2
~
c::
~
-
0
0..
c::
~
~
'-
-
a
Ql
0')
·s
.(J)
"5
~
~
I c::
a:l ::J
0
0
~
()
W 0
-.-ttS- ~
'+-
0
~
«
+-'
c:
Q)
+-'
0
-
::J
.~
~
a.. Ql
"-
.E
~::J
«
I
-c::
()
.Ql
OJ <C .t:::
E:-
_--I
• ...:::::- I
I
I
E:-
___ I
• ...:::::- I
I
I
-Q)
0
0..I
.0
'-
I I
I
I (J)
I
c::
Q)
-
a
c::
~
"-
::J
()
a
.,.-.
~
~
:J
~
The cell voltage during anode effect in the cells at Anglesey Aluminium is
around 30 V, far in excess of the decomposition voltage associated with
the discharge of fluorine. This is due to the insulating effect of the
fluorocarbon gas film. As the gas film spreads across the anode, so less
of the anode area is able to conduct the current into the electrolyte and
the effective current density increases. An increase in current density
increases the rate and polarization voltage for the fluoride reaction and so
the film continue to grow until a point is reached where the ohmic
resistance of film is less than the resistance due to the reaction. The
current then arcs across the gas film and it is assumed that under these
conditions a conductive gas plasma is formed. Electrolysis will then take
place at the gas/electrolyte interface.
l
--·CEF·l/J moles S-1 (93)
4·z F ·F
The surface coverage of the anode is the ratio of the bubble area to the
effective anode area. Assuming that the bubble is of uniform thickness
and that the gas will obey the ideal gas law, the rate of change of anode
coverage during anode effect becomes
119
The partial current density for each reaction at the anode during anode
effect then becomes
(95)
1-8 (96)
(97)
The cell voltage can be used to calculate the energy consumption of the
cell by Ohm's law
The cell voltage is not simply the sum of the electrode potentials but must
include the ohmic drop across electrolyte and the cell structure. The total
voltage drop across the cell is therefore
120
The polarisation voltages and voltage drop across the electrolyte have
been considered in detail above. The voltage drop across the cell
structure is due to the resistance of the materials used in the cell
construction and includes terms for the anodes and cathode, anode and
cathode bus, anode riser and cathode flex. The resistance of each section
was provided by Anglesey Aluminium and the data that will be used in this
work is tabulated in table 4.2.
The voltage measured in practice is the cell voltage and includes all the
above components. The cell resistance is calculated from this and can be
used to infer conditions such as the alumina concentration in the
electrolyte.
Reell (100)
In spite of thermal disturbances to the cell, the heat balance of the cell
over a long period follows equation (101) for the steady state treatment.
(101 )
where Q res is the average power input to the cell, Q R the net heat of
reaction, Q H the average heat required for the heating and dissolution of
raw materials and Qout the average total heat dispersed from the cell.
However, for the simulation of the dynamic behaviour of a cell, the heat
121
Component Resistance Q
Freeze
Electrolyte Tb
A Aluminium Pad Tm
Qc
The output from this block is the net heat of reaction due to electrolysis.
The net heat of reaction is given
n
- ~ r,I ·~HRT·
L..J . .I
(102)
i=1
In this work four reactions have been defined as the main reactions
occurring in the aluminium reduction cell. The primary reaction during
normal operation is given
[8]
122
2 4
+ -3 . rco 2dlSS
. . L1H 2 - -3 . rCF diss . L1H 3
4 (103)
aA + bB ~ cC + dD
(104)
(106)
123
The heats of reactions for [9] and [10] are similarly calculated.
In this work heat losses from the cell are defined to be the heat transferred
from the molten electrolyte and aluminium pad to the surrounding. In
figure 4.11 the losses are shown as the heat transferred from the
electrolyte to the anodes, OA' from the electrolyte to the frozen sidewalls,
Obf, from the electrolyte to the crust, 0bc, from the exposed bath to the
cavity beneath the shields, Oe' from electrolyte through the centre trench,
Oct, from the metal pad to the cathode block, Oc and from the molten pad
to the frozen sidewalls, 0mf' In this work it has been assumed that the
resistance to heat transfer between metal pad and electrolyte IS
insignificant with the result that heat transfer between these zones IS
Heat may be lost from the cell by convection and, in the case of exposed
bath, by radiation. Although the heat transfer coefficient for convection
may be calculated from the Sieder-Tate equation for variations in physical
properties and temperature this method is more suited for local heat
transfer coefficients. As physical properties change little in the range of
operating conditions and the velocity of the electrolyte and metal pad are
largely affected by the anode configuration and cell design, a
representative coefficient will be taken for each section and fixed for the
duration of the simulation.
124
ledge. This film circulates with the bulk of the electrolyte and will be
assumed to have the same composition as the bulk electrolyte. The heat
transferred to the ledge is then
(108)
(109)
The heat transferred to the cathode block is by convection from the bulk of
the metal pad to the surface of the carbon cathode. The heat transferred
from the metal pad is then
(110)
125
Zone Cs
Qcoll
Cathode Block
Alumina Zone C2
Firebrick
Qbott
_ 12 kc x C2 + K'cXCS (T _ T )
A A Cs C2
C C
126
(113)
For a centre break cell, heat transfer in the centre trench is alternately
from exposed bath after breaker action and by convection to an insulating
layer of alumina after a feed has been dumped onto the surface. As the
alumina has good insulating properties and does not remain on the
surface long enough to form a complete agglomeration with the electrolyte
by capillary action there is little formation of a frozen layer as there is with
the crust. The heat lost through the centre trench will then be by
convection and radiation for exposed bath and predominately by
convection and conduction for alumina cover. For simplicity, the heat
transfer through the alumina cover will be assumed to be at steady state
and represented by an overall heat transfer coefficient including the heat
transfer coefficients at both surfaces and the thermal conductivity of the
alumina layer. Therefore, for exposed bath
(115)
(116)
127
Uct = overall heat transfer coefficient in trench
for alumina cover
The overall heat transfer coefficient is given for steady state flow of heat
through a composite wall as
1
U (117)
The overall heat transfer coefficient across the alumina cover in the centre
trench, where the area for heat transfer is assumed to be equal across
each surface, is then
1 1
(118)
where suffix
bt refers to heat transfer from bath to alumina
ac refers to alumina cover
tc refers to heat transfer from trench to cavity
Heat is also lost from exposed bath whilst the anode is removed from the
cell when a spent anode is being replaced. The operation of anode
removal lasts for several minutes and during this period the shield above
the anode must also be removed. Heat transfer is therefore to the air
surrounding the cell which may be at a significantly lower temperature
128
than the cavity beneath the shields. The heat transferred from the
electrolyte will therefore be
Heat losses to the anodes constitutes one of the major losses from the
electrolyte. Heat is transferred from the electrolyte to the anode by
convection and conducted away from anode/bath interface to the upper
surface where it is dispersed. The replacement of a spent carbon anode in
the cell with a new one at a temperature several hundreds of degrees
Celsius below that of the cell constitutes a massive thermal disturbance.
The initial effect is the formation of a frozen layer of electrolyte upon the
anode surface which effectively insulates the anode and reduces the rate
of heat loss from the bath. This layer initially grows, the latent heat of
fusion being given up in heating the anode, until either the anode has
reached a temperature such that the heat conducted away from the frozen
surface is less than that transferred from the bath and the layer melts, or
the frozen layer contacts the aluminium pad. In the latter case, the very
low liquidus temperature of aluminium prevents further freezing and heat
is conducted through the frozen layer and into the anode. The speed at
which the anode is heated is thus increased and anodes are frequently
'set low' to take advantage of this effect. As the anode warms and starts to
draw current, Joule heat is generated within the anode further warming it.
The complex nature of this problem requires that a heat balance must be
performed on each of the anodes and for this purpose each anode will be
divided into a number of unequal zones as shown in figure 4.13. As with
the cathode heat balance, the dynamics at the anode surfaces are
considered to be most important and the zones in these regions are
consequently the narrowest. Heat transfer is assumed to be uni-directional
through the anode with no heat losses from the sides.
There are three methods of heat transfer from the electrolyte to an anode
that must be considered. These are as follows:
129
Anode zone 2
Anode zone 1
(120)
(121 )
and the heat conducted through the frozen layer into the anode surface is
(122)
Any difference between the heats into and out of the boundary layer
results in freezing or melting of the layer. As with the heat balance upon
the crust the rate of freezing can be determined by a performing a heat
balance upon the boundary layer
130
(123)
dX af
dt (124)
dM af
where dt = rate of freezing of frozen layer at anode surface kg S-1
surface m S-1
= QAS -
~.A A .("[b - TAs ) (125)
Xaf
(126)
131
- Kcov . xcov . (T - T )
A A2 A3 - Ustub · Astub . (TA2 - TcaJ (127)
A
The total rate of heat transferred from the electrolyte to the block of
anodes is simply the sum of the individual rates. For n anodes this is
(129)
The total heat lost from the cell is the sum of the individual components:
(130)
132
The overall heat balance on the pot is:
(131 )
133
Qfeed = heat required to raise feed materials to the
temperature of the electrolyte
The heat due to the voltage drop across the electrolyte must be calculated
from the cell voltage excluding the voltage drop across the cell structure.
(133)
This is the only Joule heat that enters the electrolyte. The resistance of
the cell structure causes heat to be generated within the structure itself,
for example in the cathode block, and does not directly affect the molten
electrolyte.
The only other raw material added to the electrolyte in any quantity is
aluminium fluoride. As this is added infrequently and in relatively small
amounts compared to the alumina feed, 25 kg AIF3 per day as to 90 kg
AI 2 0 3 every 70 minutes, the relative disturbance to the heat balance is
very small. In this work, the aluminium fluoride is assumed to dissolve
instantly and the heat of dissolution is assumed negligible.
134
4. 11.2 Specific heat of feed materials
The feed materials added to the cell are much colder than the molten
electrolyte at about 200 0 e and heat is required to raise them to the
temperature of the melt, about 970 o e. The amount of heat required to
raise a substance from temperature T 1 to T2 is given
(135)
1
H927 - 103.24 kJ mol-
H 1027 - 116.08 kJ mol- 1
H200 - 16.41 kJ mol- 1
and
The heat loss due to addition of aluminium fluoride is much less significant
and for simplicity, the specific heat capacity of this component will be
assumed constant. The heat required to raise AIF3 additions to the
operating temperature, based upon data from Aylward and Findlay[17] is
135
4.12 Frozen Sidewalls
Models for the precise estimation of the frozen sidewall profile have been
presented[18,19,20] but require the division of the cell into a number of
triangular sections. A mass and heat balance must be performed across
each section and the time required for these calculations makes these
models unsuitable for dynamic simulation.
The heat transferred from the molten bath to the freeze and from the cell
outer wall to the surroundings is assumed to be purely by convection.
Heat transferred through the freeze and composite wall is by conduction.
qb -- hb ·A b .(Tb - T, ) (139)
The primary area of interest is the heat flow into and out of the boundary
layer at the surface of the frozen sidewall as this controls the melting and
freezing of the frozen sidewalls. The resistance to heat transfer of the
136
freeze, composite sidewall and outer surface may be approximated to an
overall heat transfer coefficient U, by
(145)
where the heat transfer coefficient is based upon the internal area of the
freeze. The equations for heat transfer across the boundary layer then
become
In qb = hb Ab ( Tb - Td (146)
Out qt = U At ( TI - T b ) (147)
(148)
dm (149)
dt
(150)
137
The cell wall in the metal pad region consists of a sloping carbon section
and two layers of firebrick. As the cell ages the profile of the carbon
section changes as shown in figure 4.14 with the slope in the metal pad
area becoming less steep. In this model, the wall thickness in the metal
pad area will be assumed to be of a uniform thickness as is more the case
in an older cell. The assumptions made in general in this model do not
favour the unpredictable behaviour of a cell that has been newly brought
on line and the assumption of a freeze profile similar to that of an
established cell is therefore justified.
The diffusion of species between the frozen ledge and the melt during
freezing has been modelled by Hashimoto[21] and requires the estimation
of a diffusion coefficient for each species and the velocity of solidification.
A further constant, related to the diffusivity of the melt, must also be fitted
and estimates by Hashimoto range between 0.1 and 0.00001 for
commercial cells. The complexity of this model and the difficulty in
obtaining meaningful values for these parameters make it unattractive.
Chemical analysis of the composition of the frozen layer in spent cells at
Anglesey suggest that the composition of the freeze remains
approximately constant throughout the life of the cell, with weight ratios of
1.36 to 1.41 being obtained irrespective of the operating ratio in the bath.
This allows a much simpler model to be developed for the migration of
species between molten bath and frozen sidewall.
dM A , At dM
-- - . - -f (151 )
dt 100 dt
dM c, C dM
-- - f. - -f (152)
dt 100 dt
dM A'F3,' (100-A f -C f ) dM f
- .-- (153)
dt 100·(R f -1) dt
dM NaF,I dM f dM A, dM c, dM A'F3,1
-- (154)
dt dt dt dt dt
138
Frozen Sidewall
/ / / / Electrolyte / / / / /
/ / / / / / / / /
Aluminium Pad
:I Cathode
The heat transfer coefficients between the molten bath and the frozen
ledge and between the aluminium pad and the frozen ledge will vary
throughout the cell dependent upon the geometry and local conditions.
The object of this work is not to calculate a detailed freeze profile and the
hydrodynamic and temperature models already presented here would not
be sufficient for the task. Here, average heat transfer coefficients will be
assumed, with constant values for the bath and metal pad sections. This
will allow an average thickness of freeze and the total heat lost through
the sidewalls to be calculated in each section.
Both Arita and Haupin[22] have presented average values for the heat
transfer coefficients in bath and metal regions and these are shown in
table 4.3. The effective value for the metal zone presented by Haupin is
calculated by assuming a constant film of molten electrolyte to exist
between the ledge and the metal pad. Using these values as a starting
point, the values used in simulation may be tuned for individual cells to
give a stable ledge profile.
139
source bath metal effective
The heat lost from the electrolyte during an anode change and whilst the
anode is warming may be divided into three stages.
1) Heat lost from exposed bath whilst the anode is removed from the
cell
2) Heat transferred through a frozen layer of bath to the cold anode
3) Direct warming of the anode once the frozen layer has melted
The heat losses at this stage will last for the time it takes to replace an
anode in the cell. At Anglesey, this takes approximately four minutes.
2) When a cold anode is lowered into the cell, the heat lost due to direct
radiation and convection to the surroundings is immediately stopped. Heat
is now lost by convection to the surface of the cold anode. This causes a
layer of frozen electrolyte to form at the anode surface and the heat losses
will be treated in the same way as for the frozen sidewalls.
(156)
140
dX f
dt (157)
When the thickness of the frozen layer becomes greater than the anode to
cathode distance, the frozen surface will come into contact with the metal
pad. The heat transfer coefficient in this section is much greater than that
for the electrolyte and the heat transferred to the freeze will become
(158)
141
and resistance throughout the anode and freeze and this is beyond the
scope of this work. The model may be simplified by assuming a profile for
the surface coverage of the anode by an insulating layer. The resistance
of the anode can be modelled as two resistances in series, R1 relating to
the resistance across the anode and R2 relating to the resistance of the
frozen layer and the lower conductivity of a cold anode. The total
resistance of the anode is then
Ranode = R1 + R2 (159)
Assuming that the resistivity of the frozen layer is much greater than that
of the molten electrolyte in the same region gives
- PI
(160)
f·A
The fraction of the anode surface area unable to conduct current can be
related to time by an exponential function
f = a·e bt (161 )
where the pre exponential and exponential factors, a and b, are chosen to
match the current distribution profile for a cold anode obtained on working
cells.
1 n 1
= L-
i=1 Ri
(162)
I . Rblock (163)
- tot R.
I
142
The total area for current flow through the cell is the sum of the effective
areas through each anode
n
A Anodes - ~ 1,I ·A.I
L..J
i=1
(164)
Itot
AAnodes (165)
143
4 9000
3.5 8000
- ~ 7000 ;(
-"g
-...
E
(.)
CI)
>-
3 ~
6000
CI)
as
.!! 2.5
cCI) .c
N
5000 g'
...
0
..... 2 ...0
.c
.....0 4000 i
rn 0
~ 1.5 ~
c 3000 C
~
(.) ......
CI)
.c 1
t- 2000 8
0.5 1000
t , ! ! ! I I I I I I I I I I ,
0
I I
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40
Time since anode change (hours)
-12.5))-1~~
3
-4.3x10- -(%AI2 0 3 g cm- 3 (167)
4.14.3 Density
144
4. 14.5 Viscosity of the electrolyte
In lib ~ C~~} In liN" + x{ dy- x)exp[ -by {dy- x)']ex+ :~}n Ii.
+ Y( d x - Y)exp[ -b x(d x - y)2 Jexp[ -cxx]ln,uy
where
100X A1F
X= 3
In,uy = 0.01502 + 47.22t by = 8.515x10-4
XNaF + XA1F 3
4.15 Summary
References
145
3 Johnston, A. R., "Alumina Crusting and Dissolution in Molten
Electrolyte", AIME Light Metals, 1981, pp373
9 Evans, J.W., Zundelevich, Y. and Sharma, D., Met. Trans, vol 12B,
1981
146
13 Zuca, S., Herdlicka, C. and Terzi, M.,"On Porosity-Overvoltage
Correlation for Carbon Anodes in Cryolite-Alumina Melts",
Electrochimica Acta, vol 25, February 1980, pp211-216
18 Arita, Y., Urata., N. and Ikeuchi, H., "Estimation of Frozen Bath Shape
in an Aluminum Reduction Cell by Computer Simulation", AIME Light
Metals, 1978, pp59
147
23 Kvande, H. and RC:Hvik, H., "The Influence of Bath Density in Aluminium
Electrolysis", AIME Light Metals, 1985, pp671-678
148
5.0 The Control System
The overall control aim of an aluminium reduction cell is to maintain good
current and energy efficiencies. The system is, however, highly complex
and control involves the manipulation of several process variables at
regular and irregular intervals.
There are three fundamental problems which the process presents to all
process control systems. These are:
These limitations have resulted in the reliance of the control system upon
considerable manual intervention based upon the judgement of
experienced cell operators. The use of digital computers and better
mechanisation has, however, improved controllability of several but not all
of the process variables.
The main tasks that must be performed in order to control the cell are as
follows:-
149
The way in which these tasks are implemented constitutes the control
strategy.
The mass flow of materials into and out of the cell consists of the flow of
solid raw materials to the cell and the flow of reaction products from the
cell. Although the production of aluminium in the reduction zone is
continuous the addition of the major raw material, alumina, is batchwise
and it is this task that must take priority in the control scheme. If the level
of alumina in the electrolyte rises to a high level, a sludge is deposited at
the cathode causing poor operation of the cell. Ultimately, overfeeding
may result in the density of the electrolyte exceeding that of the liquid
aluminium and the two layers will invert with disastrous consequences.
Alternatively, if the alumina concentration is allowed to fall below a critical
level, anode effect occurs with the production of fluorinated carbon
compounds at the anode accompanied by a large increase in cell voltage.
The removal of aluminium product from the cell is also a batch operation
with material being siphoned from the metal pad at regular intervals. The
depth of the metal pad has an important effect upon the heat balance of
the cell and upon the current efficiency by affecting the pad stability. More
importantly though, the amount of aluminium removed gives an indication
of the rate of production of aluminium in the cell and is used to calculate
the overall current efficiency. This is used as a measure of the
performance of the cell over a long period of time and is used in tuning
other control variables.
The amount of metal removed from the cell is determined by the depth of
the aluminium pad within the cell. Both the depth of pad and the amount of
150
metal tapped are generally measured by a crude dip stick method. The
molten aluminium is siphoned off into a crucible. The amount tapped is
measured by level in the crucible and is calculated to leave a desired level
of metal pad depth in the cell. The recent introduction of electronic scales
to weigh the crucible has improved the accuracy of the tapping operation
and hence the calculated current efficiency but the control of metal pad
depth is still very crude.
151
The inter electrode gap may be altered by moving the anode bus bar
relative to the cathode. This is done by pneumatic motors that may be
activated either manually or automatically. The motors act upon the whole
anode block, not individual anodes. Any anodes that are badly positioned
relative to the block must therefore be reset manually.
For a cell operating at a constant current the cell voltage determines the
energy input and consequently temperature. Electrolyte conductivity is,
however, affected by temperature and so to some extent, temperature
affects cell voltage. Another effect of temperature is to alter the freeze
profile and consequently cell voltage. These two parameters therefore
show considerable interaction.
A low voltage drop across a cell is desirable in order to reduce the energy
required to produce the aluminium. A reduced ACD may, however, have
an adverse affect upon electrolysis by increasing the rate of back
reduction and reducing the current efficiency. At times when the market
152
price of aluminium is high it is preferable to maximise the production rate
over the energy required for production and vice versa.
The frozen sidewalls contain the molten electrolyte and aluminium and
protect the cell walls from their corrosive effects. They also regulate the
temperature of the cell by melting or freezing. A good freeze profile,
however, is necessary for good current distribution and efficiency. If the
freeze extends too far under the metal pad then the area for current flow is
reduced. If the freeze becomes too narrow then the stability of the metal
pad is reduced and ultimately the sidewalls become corroded.
The design of the cell structure and the materials used in its construction
very much determine the shape of the frozen ledges, but the profile is also
affected by the temperature and composition of the electrolyte. Although
the freeze profiles cannot be measured easily they may be controlled by
maintaining the electrolyte temperature and composition about the design
conditions. The correctness of the operating temperature is, however, only
reflected in the long term current and energy efficiency of the cell and it is
this information which is used to alter the voltage set point.
153
of control relies greatly upon the judgement and commitment of the
operator and for a reasonable level of control it would be necessary to
have an operator assigned to only a small number of cells.
154
must concentrate upon achieving a long term average voltage by
maintaining an approximately constant ACO. Movement of the anode bar
is then only necessary when a setpoint change is required or during
tapping when the metal pad level is falling rapidly.
Central to the manual control strategy is the IIpotline surveyll. This survey
is conducted by the potline general foreman and involves the inspection of
individual pots on the line on a shift to shift basis. The specific attributes
inspected are:
flame colour
pot voltage and stability
carbon dusting in the bath
liquid levels.
All abnormalities are recorded and a severity factor ranging from 0.5-3.0 is
assigned to the cell based upon the overall condition of the cell.
Assignment of severity values is prescribed by the standard practice.
On completion of surveys upon all cells within the line, the individual cell
severity values and abnormal characteristics are tallied and summarised.
A complete record of the potline operating condition is thus available for
applying compensatory corrective action on individual cells. Review of
abnormal attributes and their previous history provides information for
formulating cause removal corrective action.
In the late 1950's automation of some of the more tedious and laborious
hand operations on reduction cells began to take place. This initially
consisted of the mechanisation of crust breakers on prebake cells and
voltage control by means of mechanised anode movement controlled by
pneumatic motors activated by operators from the floor. This was soon
followed by analog electro-mechanical devices which broke crust, fed
alumina and controlled cell voltage on some predetermined schedule set
by potroom operating supervision. Automatic anode effect termination was
155
added as an additional task and the development of point feed cells
allowed smaller amounts of alumina to be fed at more frequent intervals
on a semi continuous basis.
In the mid 1960's digital computers were first incorporated into pot line
operation and process control. Although the environment is hostile to
these devices due mainly to the strong magnetic fields and lack of
cleanliness in the pot rooms this was overcome by the use of distributed
systems.
156
POTROOM
SUPERVISORY
COMPUTER
MODCOMP
ffl w ~
Each GELTROl unit contains all the input/output functions and logic
necessary to control a reduction cell. The hardware allows for cell voltage
input and operator inputs through toggle switches, push buttons and
selector switches. These switch inputs are used for altering control
parameters at the cell, acknowledging alarm conditions and viewing
setpoints and some historical data. Outputs to the cell allow anode bus
movement and break and feed operations at the tap and duct ends of the
cell.
157
equipment in the rectifier and is sampled by buffer software every 100
milliseconds. These values are digitally filtered and broadcast to each
GELTROl every 500 milliseconds where it is used to calculate the cell
resistance.
158
and historical operating data for the line and the facility is available to
make control setpoint changes on individual cells. When not in use the
terminals display any exceptional conditions that exist on the line allowing
operators to know at a glance which cells require attention.
159
different operating conditions e.g. anode effect, high alumina
concentration and sick pot conditions. A specific program must be written
to implement each of these control strategies. The current strategy applied
for half break cells at Anglesey aluminium will be described in detail here.
The voltage is sampled by the controller every 100 ms and the cell
resistance calculated from this and the line current. This value is known as
the Irawl resistance and is subject to high frequency variation. The raw
resistance is filtered through 1 and 3 minute digital filters to produce a
smooth resistance. Voltage corrections are made proportional to the error
of this value from a set point by movement of the anode beam. This beam
is moved by pneumatic motors that are operated by solenoids activated by
GELTROL. The voltage noise, leakage within the air pipes activating the
pneumatic motors and the inertia of the anodes reduces the accuracy of
the control such that small errors would cause unnecessary movement of
the anodes. To prevent this, a dead band around the resistance set point
is defined in which no anode movement is allowed to take place.
160
Cell voltag eV
Voltage + L Beam Anode beam
set po int Controller position
~-
movement
V
Under normal operation the smooth resistance is compared with the set
point every 3 minutes. If the absolute error is greater than the defined
dead band, typically 1.2 !J,Q, then the solenoid activation time is calculated
(1 )
161
5. 15.4 Anode effect
At the onset of anode effect, the cell voltage rises due to the increase in
polarisation potential and cannot be rectified by controlling ACO. Any
attempt at control by the normal AVC routine would cause the anodes to
be driven into the cathode. In order to prevent this, AVC is terminated if
the cell voltage exceeds an upper limit. If the voltage continues to rise and
exceeds 9V the cell is deemed to be in anode effect and an anode effect
termination routine is activated. This routine incorporates both alumina
feed and voltage control and is described in section 5.17.
The tapping of liquid metal from the cell lowers the pad depth so
increasing the ACO and hence cell voltage. The rate of voltage increase is
too great to be controlled by the normal AVC routine based upon a voltage
filtered over 3 minutes and use of this routine would likely result in the cell
voltage exceeding the upper limit of the controllable band.
The removal of metal from the cell is not directly detected by the controller
and any routine for controlling the voltage during tapping must be
activated by the operator. Use of normal AVC routines applied directly to
the unfiltered resistance at 18 second intervals has been utilised with
some success but often during the metal tapping operation the operator
controls the voltage manually.
162
alumina concentration, control centres upon the control of anode effect
frequency.
CE 1 I
- -.-. ·W ·t kg
100 2 z AI . F AI 2 0 3 f (2)
The cell may be under fed by either increasing the feed time or decreasing
the batch weight. Although this would theoretically give a regular and
predictable anode effect frequency, the error in the batch fed and
variations in current efficiency do not allow this in practice.
163
Theoretically, this routine could be used to always detect the time at which
a cell should be fed. Unfortunately, abnormal conditions in the cell may
give rise to changes in cell voltage that can prematurely trigger feeding
and cause serious overfeeding problems. For this reason, the strategy is
subject to a number of restrictions that are designed to prevent this from
happening.
Whilst the controller is searching for rate of change of voltage due to the
decreasing alumina concentration there must be no movement of the
anodes to mask these changes. AVe is therefore temporarily suspended
whilst search is in progress.
Noise control and the metal tapping operation could also trigger feeding
and for this reason Search is temporarily suspended should these
conditions occur. As soon as the terminating condition is removed Search
is able to restart.
The starve routine is included in this strategy to ensure that the alumina
concentration and sludge content of the cell is unable to build up. This
routine forces a scheduled anode effect by suspending all alumina feeding
until anode effect has occurred and is only activated if the cell has not has
experienced anode effect within the previous 9 hours.
The combined effects of the demand feed strategy and starve routines is
to limit anode effect frequency whilst ensuring that the average alumina
concentration remains at a reasonably low level. Low levels of alumina in
centre fed prebaked cells, have been found by some researchers to be
associated with good current efficiency [1] and occasional anode effects
may be considered a better risk than the less controllable effects of too
much sludge at the cathode.
164
5.17 Anode Effect Termination
The anode effect termination, AET, routine is activated if the cell voltage
exceeds 9V. A typical AET routine is shown in figure 5.3. The routine
consists of alternately raising and lowering the anodes in order to disperse
the bubble layer, and the addition of alumina to increase the alumina
concentration. These operations may be repeated a number of times, but
if the failure location is reached then AET is deemed to have failed and
manual assistance is required. The success of the AET routine results in
an addition to the cell of double the standard feed weight of alumina.
Point feeder cells are able to feed small quantities of alumina on a semi-
continuous basis and are therefore, in theory, able to maintain the alumina
concentration at a more constant level than a half break cell. The alumina
control strategy for a point feed cell is basically the same as that for a half
break cell with Search and Starve routines enhancing a fixed feed rate.
The short feed times, however, require that these routines are subject to
different rules to those in the half break routines.
The voltage control routine is the same as that used in a half break cell.
165
0 1 2 3 4 5 6 7 8 9
0 BRAK LOWR 7.0 W006 JUMP 15 RAIS 5.0 DUMP BRAK
1 LOWR 5.0 W006 LOPO 1 6 -RAIS 3.0 WOOS RAIS 3.0
2 W006 BRAK WOOS DUMP DUMP DUMP DUMP
3
4 WOOS BRAK W010
5 DUMP DUMP DUMP DUMP
0 1 2 3 4 5 6 7 8 9
0 DUMP BRAK LOWR 2.0 DUMP BRAK JUMP 13 LOWR 2.0
1 DUMP BRAK W001 LOPO 18 JUMP 24 RAIS 3.0 DUMP
2 BRAK LOWR 3.0 W001 LOP1 217 BRAK DUMP -RAIS 2.0
3 BRAK DUMP W002 RAIS 2.0 W002 RAIS 0.0 -W04S BRAK
4 DUMP DUMP
5
Nomenclature
The only modifications to the Starve routine for point feed cells are the
conditions required for activation. Starve will only occur if there has been
no anode effect within the past 48 hours and the time is within the
Search/Starve window.
The basic form of Anode Effect Termination, AET, in a point feed cell is
the same as that in a half breaker. A typical sequence for a point feed cell
is shown in figure 5.3.
The above strategy is that currently used at Anglesey and has been
developed over many years. Clearly there are alternatives, both in use
and under development. The alternative strategies reflect not just the
different types of cell but also different control philosophies. The basic
routines for feeding alumina and AVC are, however, much the same
throughout the industry.
166
95
Temperature 965°C
Weight Ratio 1.35
~-
- t \ 94
Calcium Fluoride
ACD
5%
4cm
~
CJ
c: 93
.!!!
CJ
w
==
1: 92
......::::J
Q)
0
91
90
1 2 3 4 5 6 7 8
Alumina Concentration (wr%)
The alumina feed strategies may be broadly split into two groups, open
loop and closed loop strategies.
Open loop strategies are the simplest and most easily implemented as the
feed of alumina is completely independent of the alumina concentration.
This is achieved by activating alumina feed with a simple timer device.
The large batch weights of half break cells often give an average
concentration of around 4%, that is close to the minimum of the current
efficiency curve. The smaller additions of point feed cells reduces the
variation in alumina concentration but still cannot operate at the optimum
using this strategy.
167
centre break cells where the build up of sludge in the centre of the metal
pad has a detrimental effect upon cell operation.
The search routine is activated 40 minutes prior to the base feed time but
feeding may not occur earlier than 20 minutes prior to base. If 25 minutes
after the base feed time Search has not been successful then feeding
occurs anyway. AVC is suspended during Search and for 5 minutes
following a feed.
168
alumina and the search and starve routines are used to provide a
reference point.
A number of operations disturb the cell and following these, the search
routine is activated in order to find the appropriate alumina concentration
for the O/U strategy. These Searches are not subject to the Search/Starve
window restriction and will start 60 minutes after tapping or anode setting.
On successful termination of Search, the cell receives a feed cycle at both
ends. It will also trigger a change to overfeed and will suspend Ave for 30
minutes. This enables the voltage to stabilise and return to the level prior
to the end of Search.
169
-
!
~
c
~
-'iii -------+----
o I
I
I
I
I
I
I I
.-:
i \ Target Alumina
:
I
: Concentration
Anode
effect Alumina Concentration
Figure 5.5 Target Voltage for Adaptive Feed Strategy for Point Feed Cells
is used as the limits for under and overfeeding. A typical voltage trace for
the strategy is shown schematically in figure 5.6.
When cell voltage passes outside the dead band, AVe raises or lowers
the bus to bring the voltage back inside the dead band. These movements
are used to trigger over and underfeeding. Four consecutive raises
changes the feed from over to underfeeding. Two consecutive lowers
changes feed from under to overfeeding. In order to prevent frequent
changes on a noisy cell, no change is made if the four raises occur within
9 minutes of a lower or the two lowers occur within 9 minutes of a raise.
The cell is over or underfed by changing the feed time relative to a base
value. Typical setting for the amounts of over and underfeed are:
The standard AVe routine is designed to return the voltage to the centre
of the dead band. In this strategy, however, it is required to bring the
voltage back to only just inside these limits so that over and underfeeding
will control the voltage within the dead band. This is achieved by halving
the value of Kp in equation (1) used in the standard AVe routine thus
halving the bus movement called for by AVe.
170
I
I
I
I
I
I
I
I
I
....... J ... . . . . .. . . . . . . . . . L..-_____ . . . . . . . . . . . . . . . . . .
I
I
I
I
I
I
CD
I
I
E
. . . . . . . .. ....;.._~ ... J ..... .
i=
I
I
I
I
I
"C
Q)
~
Q)
"C
c:
=>
overfeed/underfeed routines. This routine is activated if there has been no
anode effect within the previous 48 hours. When the feed routines are
working well, all anode effects will be due to Starve giving an average of
0.5 anode effects per pot day. Test cells have been run successfully with
one anode effect every 10 days but this requires consistent alumina feed
and good bath level control.
Ouring the tapping operation and after anode effect terminations, the ACO
is returned to its optimum value by returning the cell to a target voltage.
This strategy assumes that this is the point of optimum ACO and controls
voltage just above the point of instability. Periodically, the anodes are
lowered until noise becomes excessive and the voltage setpoint adjusted.
171
Voltage control is then by standard AVC routines. This method of control
is not currently used at Anglesey Aluminium.
5.25 Summary
This chapter has discussed some control strategies than are used in
industry and that can be investigated in simulation. The next chapter will
discuss the programming of the aluminium cell model and the control
systems and will compare the results of the simulation to experimental
data in order to validate the model.
References
3 Belyaev et ai, Izv, Vyssh. Ucheb. Zav., Tsvetn. Met., vol 4, no 3, 1961,
p67
5 Firsanova et ai, Izv. Vyssh. Ucheb. Zav., Tsvetn. Met., vol 5, no 3, 1962,
p52
6 Szeker, Acta Tech. Acad. Sci. Hung., vol 10, 1955, p91
172
6.0 Dynamic Simulation of an Aluminium Reduction Cell
One of the main requirements of the simulation is that it should be able to
be used both on-site at Anglesey Aluminium and at Newcastle University.
With the advent of cheap powerful microcomputers this presents no
problem and most of this work was carried out on a Compaq 32 bit IBM
compatible personal computer. A wide range of programming languages
are available for these machines and the choice of language is based
upon the applications for which the simUlation is intended.
The simulation is required to predict the dynamic behaviour of the cell over
long periods of time. In order that this takes as short a run time as
possible a compiled, rather than interpretative, language must be used.
The simulation of different control systems and types of cell requires that
parts of the model be altered periodically. It is therefore desirable that the
language used allows the program to be specified as a number of
modules eg process model, control system, user interface etc. If changes
are made to any module, only that module need be re-compiled not the
whole program.
173
As well as providing multi tasking capabilities, Modula-2 has a number of
other advantages over PASCAL. These are as follows:
6.2 Initialisation
The complexity of the overall process model does not allow for all of the
variables to be initialised by steady state analysis. Initial conditions were
therefore obtained by running the simulation with the approximate starting
conditions and allowing the program to converge to a steady state. The
process variables at this point were obtained as a memory dump from the
program and provide the 'steady state conditions of a 'standard ' run. This
I
Once this start point has been achieved, memory dumps of the process
variables represent the system at any point in time. The conditions at the
174
end of one simulation can therefore be used as the starting point of
another simulation. In effect, the simulation can be continued without
interruption.
Separate memory dumps are produced for the control and model
variables. In addition, the control variables may be initialised from
complete or partial memory dumps. The complete version contains all
historical data for the control system and allows the simulation to be
continued from the point at which it was stopped. Use of the partial dump,
however, allows for changes to be made in control parameters and
simulates the resetting of aGELTROl device.
175
The simulation results are also stored in ASCII format files on the
computer disk. This produces vast quantities of data that must be
processed. The best means of interpreting the data is graphically but the
commercial spreadsheet programs available at the time of this work were
unable to deal successfully with the amount of data produced. It was
therefore necessary to produce a graphics package to do this. The
program allows the user to produce graphs of combinations of parameters
and includes facilities to zoom inion areas of interests. Hard copy is
I
provided on a Hewlett Packard plotter and full details of the package are
given in appendix C.
176
of the control algorithms must be produced to allow comparison of the
simulation with measurements from the cell.
Initially the full control program for the operation of half break cells ,
including the starve and search routines, were programmed as a Modula-2
module exactly as they were implemented in Modcomp FORTRAN. This
implementation of GELTROl was used as the 'standard ' version for
verifying and tuning the simulation.
The complete set of initialisation variables used in the standard run are
listed in appendix A and the results of a 170 hour simulation shown in
figures 6.1 to 6.6. The simulation is for a mark IV half break cell operating
at a target voltage of 4.39V with Anglesey Aluminium's search and starve
routines activated. Alumina is fed at a rate of 90 kg on a base feed time of
67 minutes. The current through the cell is constant at 151 kA and the cell
is configured for eighteen 53" anodes. A comparison of the average cell
performance for the first half of 1986 with the simulated performance is
shown in table 6.1.
The voltage trace for the simulation is shown in figures 6.1 and 6.2. This
shows an anode effect frequency of 1.55 AEs per day with anode effects
occurring at 7, 16, 35, 53.5, 64.5, 83, 101.5, 121, 139.5, 158 and 169
hours. Those occurring at 7, 35, 53.5, 83, 101.5, 121, 139.5 and 158
hours are due to the starve routine and successful searches have
occurred at 28.3, 47.5, 77.2, 95.6 and 151.5 hours. The average voltage
over the simulated period is 4.41 V the deviation from target voltage being
due to anode effects. This voltage is that measured by GELTROl and is
not directly comparable with the gross operating voltage which includes
the voltage drop across the buswork between individual cells. This
accounts for approximately 0.03V at 151 kA and gives the simulated gross
177
Parameter Actual Simulation
I
I
(j)
tJ
~~
r-
o
> ~
....jl ""
I
''''111 " -J
~
~~
r'u.. L
-.....vvI " ~~
... "
~
_
~
u.
-- ~ ~~ ~~u~
3.00 , , , , I j
_ .- Target
...
o
> .... . .A .l l
3 .0 )
o 4 8 12 16 20 24 28 32 36 40 44 48
Time ( hrs )
6 . 00)
~
u
o
> ./. A I .A L
3 .00
50 54 58 62 66 70 74 78 82 86 90 94 98
Time (hrs)
6 . 00)
.
u
o
> ~ . . I ! . -"
3.0 0
100 10 . 1 09 11 2 116 120 124
T1me ( hrs)
128 132 136 140 , .. 149
Figure 6.2 Simulated Voltage Profile for 151 kA Centre Break Cell
operating potential as 4.44V. This compares very favourably with the
average value for the first half of 1986 of 4.48V[11.
Metal tapping was set to a 24 hour tap rate with a target pad of 14 cm.
The first tap occurred at 22 hours and the total mass of metal tapped over
the 170 hour period was 7980 kg. The amount of metal is used by the
industry to calculate the overall current efficiency over a period of time.
The calculation is based upon the assumption that the metal tapped is
equal to the net amount of metal produced over the period giving
CE I
Mtap - -. ·W·t (1 )
100 Z.F
AI
AI
rearranging and inserting values for W AI, ZAI and F gives the overall
current efficiency
M
CE = 1.072x106.~ (2)
I· t
The current efficiency by this method for the seven day period simulated in
the 'standard' run is then 93.7 % •
This value for current efficiency is somewhat higher than that for
commercial cells at Anglesey Aluminium which operate at an overall
efficiency of about 90%. The discrepancy between the two values is due
to the short time period of the simulation. In the industry, the current
efficiency is calculated as quarterly and yearly averages. This removes
discrepancies caused by the method of tapping.
Although the aim of tapping is to return the pad depth to a target level the
amount tapped is not directly proportional to the deviation from this set
point. The schedule used for determining the metal tapped is shown in
table 6.2
178
Deviation from target Metal tapped
em kg
0 1140
1 1320
>1 1500
179
..
w
u
50 LI____ ~ ____ ~ ______ ~ ____ ~L_ ____ ~ ____ ~~ __ ~ __ _L~ ______ ~ _ __L_ _ _ _ _ _ ~ ____ ~ ______ ~ ____ ~L_ ________ ~ __ ~ __ ~ ________ ~ ______ ~
r-l __ -
6.00I
~~ ~r--~
-- --..- ...-- l 1 r---'-
E
u
o
u
<
0 .0 J
o 17 34 51 68 85 102 119 136 153 170
T1me (hrs)
8 . 00
...o
N
<
. .
:t
o oo LI______~______~______~____~______~____________________~______~____________________~______~--------------------~------~--------------------~
o 17 34 51 68 85 102 119 136 153 170
T1me (hrs)
Figure 6.3 Simulated Current Efficiency, Concentration and ACO Profiles for 151 kA Centre Break Cell
The value for the simulation is therefore 14.6 DC kWh I kg which is close
to the actual value for the mark IV cell of 14.8 DC kWh I kg.
The weight ratio of NaF/AIF3 is shown in figure 6.5. The graph shows a
sharp rise in ratio at anode effect corresponding to the evolution of
gaseous fluorides by the anode effect reaction and also an increase in the
rate of volatilisation due to increased temperature. The addition of bagged
aluminium fluoride to the cell is characterised by a sharp fall in ratio. Over
the period of simulation, the ratio is stable and fluctuates about an
average of 1.23 0/0.
The average values for temperature, 971°C, and Ratio, 1.23, obtained by
simulation both compare favourably with the historical data for the year
1986, 969°C and 1.22 Ratio, although the simulated temperature is 2 °C
greater than the actual value. This is due to the slightly higher value of
Ratio which raises the liquidus temperature hence requiring a higher
temperature to maintain the heat lost through the sidewalls of the cell.
180
0
. ~- . --
=-=--
-.: . ~ r--..
r-I
r0
L.[)
r-I
\0
r0
r-I
...:::J
Q)
-...ctI
Q)
(j)
r-I
a.
E r-I
~
en
:::J N
"0
:::J 0
0- r-I
:.:J
L.() v
ro (1)
E
'r
ro I-
\0
o
o o
o avln+t:?vladwaJ..
o
Ii
1. 25
o
..
.
at:
:E
1.20 ~1 ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~
1000
u
..
L
....
::>
..
L
.....
Co
E
950 LI____ ~ ____ ~ ______ ~ ____L __ _ _ _ ~_ _ _ _~_ _ _ _~_ _ _ _~_ _ _ _~_ _ _ _ _ _~_ _ _ _L __ __ _~_ _ _ _~_ _ _ _ _ L _ _ _ _~_ _ _ _~_ _ _ _~~_ _ _ _~_ _ _ _L __ _ __J
Figure 6.5 Simulated Ratio vsTemperature for 151 kA Centre Break Cell
commercial cells but is calculated from the voltage drop across the
electrolyte by the relationship
(4)
The relationship between the gross operating voltage and the voltage drop
across the electrolyte is given by the equation
(5)
The calculations for the actual and simulated values differ only in the value
used for the overpotential. In Anglesey Aluminium's calculation, the
overpotential is calculated from a Tafel expression for the reaction
[1 ]
based solely upon the current density at the anodes. The effect of alumina
concentration is disregarded and the relationship appears to have been
determined for a concentration of 5 wt% alumina.
181
30.00 rl------------------------------------------------------------------------------------------------------------------------------------,
E
u
VI
...a.
~
.
c
Bath
10.00 1L-_____ L_ _ _ _ _ _ ~ _ _ _ _ _ L _ _ _ _ _ _~ _ _ _ _ _ L _ _ _ _ _ _~ _ _ _ _ _ L _ _ _ _ _ _ ~ _ _ _ _ _ L _ _ _ _ _ _~ _ _ _ _ _ L _ _ _ _ _ _ ~ _ _ _ _ _ L _ _ _ _ _ _ ~ _ _ _ _ _ L _ _ _ _ _ _ ~ _ _ _ _ _ L _ _ _ _ _ _ ~ _ _ _ _ _ L _ _ _ _ ~
10.00
.
u
..,'"
-'
.
..,
.
"-
5 OO LI____ ~ ____--'______ ~ ____ ~ _____'_____ ~ ____ ~ ____ ~ ______'______ ~ ____ ~ ____ ~ _____'_____ ~ ____ ~ __ --~------'------~----~-----'
Figure 6.6 Simulated Bath and Pad Depths for 151 kA Centre Break Cell
the very low anode-cathode distances due to anode movement during
anode effect termination.
The 'standard' run shows good agreement between the simulation and
industrial data and can therefore be used as a reference for the validation
of the more detailed dynamic behaviour of aluminium reduction cells. The
lack of comprehensive data for individual cell requires that the model be
tuned to average conditions at the Anglesey plant. In validating the effect
of disturbances to operating conditions it is the resulting trend in the cell
parameters that is of importance rather than the absolute values attained.
Simulation results will be compared with both experimental data and the
results obtained by other models.
The results of measurements for a 150 kA centre break cell with a two
hour feed cycle[4] are shown in figure 6.10 and concentration curves from
the simulation by Ek and Fladmark[5] in figure 6.11. These all have similar
shapes over individual feed cycles.
182
6.00
M
o
N
r-
«
~
+J
3
l .ool~~~~~~~~~~~~~-L-L-L~~~~_ _L-L-~~~~~-L-L~~~
'"'-:J
o
.r:
o •
- <I)
E
<J CO
<J
~6
.-'"
'">- '-
ro <I)
c '-'
<! <I)
<J: -ro
"l:
ro
$Y~O
u
E
.-Ec
<I) :J
.r:
u
-<!
0 <J N
o
o
N
~ Co
III
c: ...,CIJ
<l: CIJ
::E
~
III III
U
.-Ec:
E
CIJ :J
-
CO
. "----CL OJ ..c
u
o <J
~
<l:
<J .
-
-.D
~eaJ8 ~ ~6J-'
,-<P <J
'.~~-~~~
o
N
Dutil Temperature, o C Current Efficiency, %
lD lD lD lD (0 lD
~ U1 U1 (j) (j) -...J 00 00 00 to
U1 0 U1 o U1 () -...J 'Y' to 0
'\
00
.~ r>
?
o
o \
:-rJ
.~
cO"
c:: '"'
(t) ;:)
~
)
o
!J)
-.&.
o \
a I-"
N \•
JJ o .~
CD
g.
C"')
......
::>
1:1
:::: ,>nn \ Ji. [>
Effect
0" N /." i ::=--" <-
:::l o
o
.\
~ •
.
[>U---::-O
~
:::::::
A
o · '; \
/~ s:E:
CD
ti3
o
'" f>l> _cr
~ ~:'> ~~O/"
......
0" O'l Anode Effect
:::l
~"
o
c ( J>--U
S '\ . y-a-CYcY
o~
00
o
I
o .'- :;;= t>~ /\node F.f fect
~
c::
....,
OJ
I-"
o
o
o
I• "
\
• y~-1:::U
(t)
III
I-"
Iv \
~ \
;:)
;:)
•
'c> ~
~
>
CD
CD :.:
\.
Q..
~
N
o
I' '. "--
\.~
/L>~·D.-'
C"')
<t~
o
r'~ D»I> o.~.U Anode Effect
A
o
o
.(\ ~.
~ ? 0"
.<"0 Metnl
. ....~ L.J
(j)
;:)
Cell Voltage
'\
.~. ~ [> J ~.
~
o ~ ... \.Y' AI umi na,%
~
m
I~. 5 r ------,----,---.------__
4.4
4.00n-----------7G---------~-----------~
12 18
t (hr)
The cell voltar.e.
10.0r----------r--------~-----------
7.0
6.0~-----J---- _ _ _ L _ _ _ _ _ _~
o 12 18
t (hr)
Thic~~ess o~ unpcr freeze.
5.0r-----------~r-----------~--------------
p (::)
ox
3.0
2.0~ __________~~----------~------------~
o 6 12 lR
. t (hr)
Consentration 0: Al 0 in
2 3
t~e bath
183
5 . 00 1
Anode
Change
1 . 00LI______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______L __ _ _ _ ~_ _ _ _ _ _~_ _ _ _~_ _ _ _ _ _L __ _ _ _~
.s::
\
.,
c-
..,
Ol
..
.s::
.,
II)
w
Anode
Change
E ~
15 0 0
85 0 85 .2 85 . 4 85 . 6 85 . 8 86 . 0 86 . 2 86 . 4 86 . 6 86 . 8 8 7. 0
Time (hr,)
4.02 4 .12 4.22 4.32 4.42 4.52 4.62 TMR FED EVENTS
0: 0 : . . . . : . . . . : . . . . : . . . . : . . . . : .. .. : . . . . X. Z .. : . . . . : . . . . : . . . . :
'" ........... 25 68
0: 3 x Z 28 68
0: 6 TX SR: 31 68
0: 9 T XZ 34 68
0:12 T Z 37 68
0:15 T ZR 40 68
0:18 T ZR 43 68
0:21 T Z 46 68 L
0:24 T X S R 49 68
0:27 T SXR 52 68 L
o : 30 : . ... : .... : .. . . : . ... : .... : .... : .... T .. . . : .. . . : SRX. : .. .. : .... : .... : .... 55 68
0:33 T S R X 2 68 A
0: 36 X S R T 5 68 A
0: 39 X T R S 3 68
0:42 X Z T 6 68
0:45 X SR T 9 68 R
0:48 X S R T 12 68 R
0:51 SR T X 15 68 R
0:54 Z T X 18 68
0:57 Z: T X 21 68
1: 0 : . . . . : . ... : . . . . : . . . . : .. .. : ... SR . .. . T. X .. : . . . . : . . . . : . . . . : . . . . : . . . . : . . . . 24 68
1: 3 Z: T X 27 68
1: 6 Z: T X 30 68
1: 9 SR T X 33 68
1:12 SR T X 36 68
NOTATION
The voltage change during anode effect is shown in more detail in figure
6.14. The experimental data used here is from a 10 second scan for cell
1A3. The simulation again shows good agreement with the industrial data.
184
10
- - - Simulation
9 - - - - - - Target \()Itage
8 o measured
-->
Q)
7
C)
as
==0 6
>
5
JIIIIIIIIIIIIIIIIIIIIIIIII~_- - - _L - -- - -- - -- - -- - -- - -- - -- - -- - -- - -- - -- - -- - --
3
29 31 33 35 37 39 41 43 45
25 27
Time (minutes)
: tt '{I
.J
CO
\ ('tJ
q::
(()
(D
"
"a"-
N
N
' 41
...,t-..
~
'-
...
: t€ '{I
-
.,--
0
0..
'+....
r 0
I- .... 0')
a .:I ~
a ·S
(T) '", Q
.;;
'4
'4. ... : ~l '{I Q
~
., (
~
N
-
I\J
~
'(
I)
~ 0')
,\ .c::
\., ... t::
• f'- ., : .1 '{I .g
Q)
(D ()
~
(T)
<
~
: . '{I Q)
I- ()
a ~
-c::
'-
ll. It: r
~ ~ ctI
oc
a ~~
,. ... .~
ll.1 CJ)
'tI
-"
Q)
.J i~
:r '
(
16 ---'r "K'i!1
ex:
" a... 1- l{)
.,--
J
•, uI I
•o ·
0~ '
eO
-
\0, n::
Ul
(J)
: tt ~I
@
::J
.0')
t
LL
.;:
..., t:
(J)
...
w : t€ 'i!l
a:: .
"'«
<.
a::
I
I • r tl '(I
0'12 v'02 8 ·61 2 "61 9 "8t 0"81 po ·L 1 8"91 2"9t 9 "S l o 'S t
S~HO O~Jlt4 3JNY1SIS3~
E
u
.
r:
0-
o•
.
r:
.·
r --------------------------------------------------------------------------------------------------------------------------------------------------------------------------,
10.00 '
•
"•"
.J
·
".
S.OO L'____________ ~ ______ ~ __________________ ~~ ____ ~ ______ ~ ____ ~ ____________ ~ ______ ~ ____ ~ ______ ~ ________________________________ ~~ ____ ~ ______ ~ ____ ~ _____________'
o 4 8 12 16 20 24 28 32 36 40 44 48
Time (hrs)
98S 'r----------------------------------------------------------------,---------------------------------------y---------------------------------------------------------------,
·
L
~
·•"
L
·
(
S' IL____ ~ ____ ~ ______ ~ ____ ~ __________ ~ ____ ~ ____ ~ __________ ~ ____________ ~ __________ ~ ____ ~ __ --------~----~----~~----------~----------~ _____'
o 4 8 12 16 20 24 28 32 36 40 44 48
Time (hrs)
Figure 6.16 Change in Temperature and Bath Depth after a Simulated Decrease in the Target Voltage of O. 1V
25.00 1 .1 L
.
.c
0.
o
o
..
.c
·
"
15.00 tl______
o ~ ________________________
4 8
~ ______
12
=_----~----::----~----~~----~----_;:_--------__~~----~----:;----------~:_----~----~;_----~----::~--~----~
16 20 24 28 32 36 40 44 48
Time (hrs)
10.00 rr- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ,
D•
"•
.J
·
"
Q.
5. oo,L~_~~_~~_~~____=__::__~________::__~_::_~_:;_~___;;_~___:~~____::;__~____::~~_:,
o 4 8 12 16 20 24 28 32 36 40 44 48
T1me (hrs)
990
·..,
L
L
·•
0.
·
E
960L
I ________________________ ~ __________________ ~ ______ ~ __________________ ~ ____ ~ ______________ ----------~----~~----~------------------~----~----------~
o A 8 12 16 20 24 28 32 36 40 44 48
T1me (h,.. .. )
Figure 6.17 Change in Temperature and 8ath Depth after a Simulated Increase in the Target Voltage of O.2V
In the simulation of the step change in target voltage to 4.49 V the bath
depth initially drops as the anodes are raised in the bath then increases as
the ledges and crust melt into the electrolyte due to the overall rise in
temperature associated with a voltage increase.
The effect of anode movement upon bath depth is shown in figure 6.18
which shows the variation in depth due to the large anode movements
during anode effect termination. As expected, lowering the anodes is
associated with a rise in bath depth due to displacement of the electrolyte
by the anode blocks. Returning the anodes to their position prior to the
anode effect retu rns the bath level to its previous depth.
185
28.00
E
u
.,
.r=
c.
<II
o
.,
.r=
'"
co
1 8 . 001~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____ ~ ______ ~ ____- L_ _ _ _ _ _ ~ _ _ _ _~
.,'"
'0
>
t ~
3.00 1L-____L __ _ _ _ ~
v
_ _ _ _~_ _ _ _~_ _ _ _~_ _ _ _~_ _ _ _ _ _~_ _ _ _L __ _ _ _~_ _ _ _~_ _ _ __ L_ _ _ _~_ _ _ _~_ _ _ _~_ _ _ _ _ _~_ _ _ _L __ _ _ _~_ _ _ _~_ _ _ __ L_ _ _ __J
Figure 6.18 Simulated Effect of Anode Movement During Anode Effect Termination
1. 25
o
.
+"'
'"
+"'
os
1 . 20 L'____L __ __ L_ _ _ _L __ _ ~ _ _~~_ _~_ _~_ _ _ _~_ _~_ _ _ _~_ _~_ _ _ _~_ _~_ _ _ _~_ _~_ _ _ _L __ __ L_ _ _ _~_ _~_ _~_ _ _ _~_ _~_ _ _ _~_ _~_ _ _ _~_ _~_ _ _ _~_ __ L_ _ _ _L __ _~
0.0 1.6 3.2 4 .8 6.4 8.0 9.6 11. 2 12.8 14.4 16 . 0 17 .6 19. 2 20.8 22.4 24. 0
Time ( hr s)
990
'"'-
.."
+"'
'-
'a.E"
'"
I-
---"-~~~~~~----'----;;-~~~~~-'---~~::-~~--'---'------'-------'-----'-----'--
i
960:-
o 0 1.6 3.2 4.8 6.4 8.0 9.6 11.2 12 . 8 14 . 4 16 . 0 17 .6 19.2 20.8 22.4 24.0
Time (hrs)
0
+:;
ctS
a:
1.25 r
~ / \ ~ - -
_ - - -
/ , --
- - 510:=0.0028 Ihr
,,0.0041 (O C
-
~
ctS
en
1.2
1.15
1.1
6' 1000
-
0
(]) 980
-
"-
::J
ctS
"-
(])
a. 960
--
E
~ 940
1000
....-..
~
~ 900
---.....
~
::::l
c. 800
c:::
.....
rn
ANODE
Q) 700 CHANGE
vsp'*"-· .. • • • • • • • •
VSP VSP AE ~.2 Voit8
I
ITt AE .1 .VoIt8
TAP AE
+0.1 .VoIt8
TAP
600
0 10 20 30 40 50 60
Time (hours)
- -- -- --
1.25 f ~-----
-
-- -- --
0
:;::;
~
1.2 r slope = 0.0033 Ihr
II:
1.15
-- " 0.00861° C
1.1
0 10 20 30 40 50 60
-
()
0
1000
....
Q) 980
--....
:::J
(Ij
Q)
960
I ~ ~~ "'-/ slope = 0.38° Clhr
a.
E
~
940
0 10 20 30 40 50 60
-
I
1000
I It I 1270 kWh
~
----
.::&.
:::J
900
a.
c 800 ~ vsp
+0.1 Volts
--
(Ij
Q)
700
Tap,AE
----...----
..
Anode Change
-
.
Tap
- --VSp
.
+0.1 Volts
vsp
~.2 v~ .. TAP,AE
I
600
0 10 20 30 40 50 60
Time (hours)
1.25
0
:;::; ~ = 0.0053
~
II:
1.2 r ~ ... ...
...... slope Ihr
" 0.0088 I" C
1.15
1.1
0 10 20 30 40 50 60
-
()
0
1000
0) 980
'-
:::J
......
'-
C13
I c- ..G---...c.
~~---V slope = 0.6
0
C/hr
O)
a. 960
E
~
940
a 10 20 30 40 50 60
-
I
1000
~
-
~
......
:::J
a.
900
800
c
......
VSP 1\ ;\ ~2
vSP
Volt
In the simulation the ratio has increased by about 0.03 for a temperature
rise of 6 °C giving the temperature dependence of ratio as 0.005 °C-1. The
Anglesey data shows some agreement with this, with the temperature
dependence of ratio in figures 6.20 to 6.22 ranging between 0.005 °C-1
and 0.0075 °C-1.
6.6.11 Temperature
Typical temperature time curves for two different feed frequencies are
shown in figure 6.23[61. The temperature is measured at the tap end and
due to imperfect mixing in the cell the effect of tap end feed is greater than
that for duct end feeding. In order to compare this to the simulation which
assumes the electrolyte to be thermally well mixed, the tap and duct ends
of the cell will be assumed to be independent of each other. The graph
can then be compared with simulation of about a 70 minute feed time.
186
9
I
0
0
I
W 980
0::
:;)
~
~ \FEED DUCT
0:: END
~ •
w
.-
975
\
~
4
~
I
I \.V /
~I\"
960
I
0 20 40 60 80 I~ 120 140 160 180 200
TIME - MINUTES
Figure 6.23 Typical Temperature- Time Curves for Two Different Feed Modesf6)
o
N
~ ~
~
~
(.)
'" '" ~
Q)
~
~
:::J
~I'-..
~ ~
(()
~
~
~
~
0
Q)
~
0
"'
L-
.<:
~
"'
L-
.<:
e
a..
QI QI
E
E
.... >= @
:::J
......
~
Q)
CD CD Q..
0 0
E
~
~
Q)
......
'" ~
0
~
0
:::J
.§
CI)
~
~
C\J
~ 0 (()
0
@
:::J
.~
l(
N
N
0
,;
% "
Z 0 LIJ
~~
-' w Wu~
-' 0 OWl
0 LIJ LIJ LIJ
::f~ ~
LIJ W OI.4..W LIJ
~ z 0 Z~O 0 0 I-
wI- ~ ~ ex ~ ~W~ ~ ~ ~ "-
o(/)
W W W W W W LIJ Cl.:) W W
~
"-
0
~
~
"-
o
~
l-
~ e "- ~
~5
0
~ ::1
w
a::
:J
~
ex
a:: 97
w
a..
~
w
~
~
t·37
ex
CD
CELL VOL. TS INCREASED - 0.1
BA TH TEMP. - 976°C A BATH TEMP.- 981°C AVG.
12 13 \4 15 16 17 18 19 20 21 22 23 24 2 3 8 9
HOUR OF DAY
Figure 6.25 Average Temperature in an Industrial Cell After a O.tV Increase in Operating Voltage/61
a 0.2 V increase in target voltage upon the temperature and ratio of the
electrolyte are shown in figure 6.19. It is difficult to judge the rate of
increase of temperature from this graph due to the fluctuation throughout
a feed cycle but the overall increase in average temperature is about 6 0C
from 968°C to 974 °C. The accompanying increase in ratio occurs over
the first 13 hours of simulation which suggests that this is the length of
time that the cell takes to reach its new average temperature. Based upon
this, the temperature rise is 3 °C / 0.1 V at a rate of 0.46 °C hr1.
The temperature data for Anglesey's cells shown in figures 6.20 to 6.22
was measured at 2 hourly intervals and it is therefore difficult to determine
an accurate rate of increase of temperature due to changes in cell voltage.
The overall rate of temperature rise for an increase of 0.2 V is in the range
0.38-0.67 °C hr1 based upon the time taken for the cell to achieve a new
steady temperature. The overall temperature rise, however, is about
7.5 °C / 0.1 V.
The differences between the two sets of industrial data reflect the different
dynamic behaviour of different designs of reduction cell. There is
agreement, however, in the overall temperature increase with increase in
cell voltage. The temperature rise in simulation is somewhat less, being
approximately half that found in commercial cells at 3 °C / 0.1 V.
(6)
The temperature increase in the cell is then determined by dividing the net
power input by the thermal mass of the cell
dT ~QI.net (7)
-
dt M bath . CPbath + Mpad . CP pad
187
where Mbath = mass of electrolyte
CPbath = specific heat of electrolyte
Mpad = mass of aluminium pad
CP pad = specific heat capacity of aluminium pad
The main parameters affecting the temperature increase are therefore the
volumes of aluminium and electrolyte in the cell and the composition of the
electrolyte. These directly affect the thermal mass of the cell and
contribute to the heat transfer and heat loss terms in the determination of
the increase in net power input.
(8)
(9)
The area for heat transfer is directly related to the depth of the electrolyte
and aluminium pad and hence the volume of these two sections. Both the
heat transfer coefficient and liquidus temperature of the electrolyte are
dependent upon the composition of the electrolyte. As these parameters
are not available for the industrial cells it is impossible to determine the
accuracy of the temperature rise predicted although it is to be expected
that the Anglesey Aluminium cells are operated at a composition close to
that used in the model. It is therefore likely that the heat transfer
coefficients used in the model require some tuning with the simulated
results being compared to additional industrial measurements where the
above parameters have been recorded.
188
A contributing factor to this discrepancy is also the increase in current
efficiency in the simulation causing a reduction in the heat input to the cell.
In practise, an increase in temperature is usually accompanied by a
reduction in current efficiency caused by an increase in the rate of the
back reaction of aluminium to alumina and an associated increase in heat
input. In the simulation the increase in the interelectrode gap to acheive
the new target voltage appears to give a greater increase in current
efficiency than the decrease accorded by the rise in temperature. This
reduction in the rate of the back reaction reduces the heat output of the
chemical reaction and hence the net heat input. It may be that the model
for the dissolution of aluminium into the electrolyte has too great a
dependency upon the interelectrode gap and not enough upon the
temperature and this should be investigated further.
189
t-
o
W
Ll..
Ll..
W
W
1030 o
o
z
c:t
1020
;010'
100
o
w
W
Ll..
980
t
97Gl -
35 30 25 20 15 10 5 10 15 20 20 30 3"5 40
MIN BEFORE AE AE MIN AFTER AE
u
...
.....
.,
::I
..
c.
E
>-
.,
..
.c
<D
960~
1 ____ ~ __ ~ ____ ~ __ ~ __ ~~ __ ~ __ ~ ____ ~ __ ~ ____ ~ __ ~ ____ ~ __ ~ ____ ~ __ ~ ____ ~ __ ~ __ ~~ __ ~ __ ~ ____ ~ __ ~ ____ ~ __ ~ ____ ~ __ ~ ____ ~ __ ~ __ ~~ __ ~
o 4 8 12 16 20 24 28 32 36 40 44 48 52 56 60
Time (min)
30.00
[G
...'"
o
>
O. OOLI__~____~__~__~__~__~____~__~__~__~__~__~~__~__~__~__~__~____~__~__~__~--~--~~--~--~--~--~--~----~~
o 4 8 12 16 20 24 28 32 36 40 44 48 52 56 60
T1me (m1n)
Since this work was completed at Newcastle University, the simulation has
been in use at Anglesey Aluminium to investigate the cell operation. An
early study considered the effect of changing the metal tap interval upon
the temperature and current efficiency of the cell. The results of the
simulation study are shown in figures 6.28 to 6.30.
It was unfortunate that the tap at around 50 hours coincided with the cell
being in search mode. The automatic tapping routines are unable to
successfully restore the cell voltage to its target value under these
conditions and this leads to an artificially long period of high voltage
following the tap. The average cell conditions are distorted by this, but
comparison of the tapping regimes is still valid.
190
100
en 90
fa
<D
C\I
N
OJ 80
>-
(,)
cQ)
'0
!E 70
UJ
C
~
~
() 60
50
o 10 20 30 40 50 60 70 80 90 100
30
~
en
fa
E
~ 25
oC\I
o/j
E
~
lli
,....
20
15
10
o 10 20 30 40 50 60 70 80 90 100
990
985
en
>
co
<D 980
o
"-
OJ
975
~
~
"§ 970
Q)
Cl.
E 965
~
£
iii
CD
960
955
950
70 80 90 100
o 10 20 30 40 50 60
01 90
>
OJ
CD
C\J
(\j
Ol 80
>-
()
c
Q)
'0
:E 70
w
C
~
::::l
0 60
50
0 10 20 30 40 50 60 70 80 90 100
30
01
>
OJ
E
()
-=t 25
ci
C\J
0<:1
E Bath depth
()
-=t
o.ri 20
,....
.c
Cil
CD
u Metal pad
c 15
OJ
~
Q)
Tap Tap Tap
~
10
0 10 20 30 40 50 60 70 80 90 100
990
01
iO 980
CD
ci
f'..
Ol
~ 970
::::l
~
Q)
Cl..
E
~
.c 960
Cil
CD
950
o 10 20 30 40 50 60 70 80 90 100
90
dl
~
C'?
N
en 80
()
c
Q)
T5
~ 70
w
C
~
~
0
60
50
o 10 20 30 40 50 60 70 80 90 100
30
dl
~
E 25
0
v
ci
N
o/l
§
N
r--: 20
.J::.
11i Metal pad
CD
"0
c 15
C1l
Tap Tap
~
::2
10
0 10 20 30 40 50 60 70 80 90 100
990
01 980
>
C1l
<D
0)
<D
en
~ 970
~
~
Q)
a.
E
~
.J::.
11i 960
CD
950
70 80 90 100
o 10 20 30 40 50 60
6.8 Conclusions
It can be concluded from the results of this work that the simulation
successfully simulates both the behaviour of the cell and the control
system, although further tuning of the heat transfer parameters would
seem to be indicated. It would also be advisable to further investigate the
dependence of the current efficiency upon the cell temperature as this
does not seem to be sensitive enough in the current model. The modular
nature of the program makes the changing of individual algorithms
straightforward and facilities have been included to allow the integration of
additional user menus an easy process.
191
al6c3 al6c4 al6c5 al6c6 al6c7
Table 6.3 Average Parameters From Five 170 Hour Simulations of a Half
Break Cell
References
192
7.0 Recommendations For Further Work
An area of work worthy of investigation with the model is real time control
of the aluminium cells using the model to predict parameters in the cell
that cannot be measured or are not available on a continuous basis, for
example, the alumina concentration. By connecting a real time version of
the simulation in parallel with the cell as shown in figure 7.1 the simulation
could output the required control signals based upon continuous results of
simulation. This would allow the optimal control strategy to be developed
continuously for any particular cell.
The model would continually correct its predicted parameters based upon
any available data from the cell, eg voltage, tap time, length of anode
change, periodic temperatures and electrolyte composition.
The first stage in this work was to develop a real time version of the
simulation that was controlled by a GELTROl panel. This was then
connected in parallel with a working cell with the aim being the
comparison of available experimental data against the predicted
parameters. This open loop system would allow the tuning of the model to
the particular cell aid in the validation of the model.
1. In order for the model to be used successfully for real time control, the
simulation must run several times faster than real time. In order to co-
ordinate the model and the logging of experimental data for comparison a
real time kernel must be developed. The kernel has a number of functions
as shown in figure 7.2.
193
Voltage
Fe edtime
Cell Model
Feeds
Operator
Interface
Screen I Keyboard
Datalogging Ioo-o------t
IControl
Disk I Printer
Figure 7.2 Use of Model in Real Time Control
The real time kernel must schedule the events depending upon a pre-set
timing and priority. In the simulation the order priority of the tasks is
2) Model For good agreement between the model and reality the
process must be simulated at a frequent interval. Additionally, it will be
necessary to run the model a number of times adjusting various
parameters until the results agree with the logged data. The
simulation of the process must therefore have a high priority and must
also be fast. It may not be necessary however to solve the equations
faster than the base interval for the control signals. The noise in the
measured data, eg Voltage, prohibits the use of raw data logged at
the same interval as the control signals and a filtered signal should be
used. It would be sufficient to compare the 'results' at intervals of say
30s, so long as all the control actions are known.
A very simple real time kernel can consist of a timing loop that calls each
of the tasks in turn , where the time taken for each task is known. If the
lower priority tasks are known to take longer than the interval required for
194
the higher priority tasks then the higher priority tasks may also be called
by the lower priority tasks. For example, if there are three tasks
DATALOGGING, MODEL and STORAGE with descending priorities the~
the kernel could be written
Call STORAGE
Timer loop
Call DATALOGGING
Call MODEL
Call DATALOGGING
MODEL:
Call DATALOGGING
Call DATALOGGING
DATALOGGING: ...
where the ellipsis represents the required code for each task.
This approach was originally tried in order to validate the response of the
model to the control actions of the CELTROL computer. The program was
written in PASCAL and connected to the CELTROL panel via a TeS
Signal processor.
195
The CElTROl panel outputs digital signals for the control of the raising
and lowering of the crust breakers, the operation of the feed dump and the
raising and lowering of the anodes. It inputs an analog signal of the cell
voltage in the range 0 to 30 V. The CElTROl 10 signals are represented
schematically in figure 7.3.
The simulation was connected via the TCS signal processor to the
CElTROl panel. Communication between then signal processor and the
computer was via the RS232 serial port. This is shown schematically in
figure 7.4.
The simulation outputs the simulated cell voltage to the CElTROl panel
which returned the control signals for the simulation. Initial trials with an
isolated CElTROl panel were successful with the simulation responding
correctly to the signals from the panel to raise/lower anodes and feed the
cell. The next stage was to connect the simulation in parallel with a
working cell so that the simulation would receive the same control signals
as the cell in order that the results from the simulation might be compared
with the data from the cell. This configuration is shown in figure 7.5.
In this case the CElTROl computer operated upon the actual cell
voltage. It was also necessary for the simulation to use the actual line
current that was being supplied to the cell. This was obtained as an
analog signal through the TCS signal processor.
A more sophisticated real time kernel can be written using languages that
allow concurrency, that is allow multiple tasks to be run simultaneously.
196
R'
alse A..nodes
Lower An odes
...
Line Current
... Crust Breaker (Tap End)
...
CELTROL Crust Breaker (Duct End)
...
...
Dump (1ap End)
Ce II Voltage ......
Dump (D uct End)
......
mrs
CELTROL
The flow diagram for the scheduler is shown in figure 7.6. The scheduler
should be assigned to the timer interrupt DOS interrupt 1CH which is
called 1024 times every second. The timers should be set to the number
of timer ticks required between each call to each process and the timer
tests must be performed in order of priority. The lowest priority task may
be run continuously as the main process to be interrupted. In the
aluminium simulation, the operator interface and data storage may be run
as the main tasks with the high priority tasks of data logging and the
model being handled by the scheduler.
197
Return to
decrement timers current process
no
no
.....-
Return to
current process
is
DOS in
critical
\ . ------
section?
no
Reset timer P
Transfer to
process P
In the model a much greater variety of variables are available for control,
most of which it is not possible to measure continuously on the process.
An example is the control of alumina feed based upon the concentration of
alumina in the cell. At present it is not possible to continuously measure
the alumina concentration in the melt and feed back control methods to
the alumina feed time are based upon the rate of rise of cell voltage as it
approaches anode effect. If the alumina concentration was known,
however, this could be directly fed back to the alumina feed algorithm to
control the alumina concentration at much more precise levels.
This use of the model is shown in figure 7.7. The simulated alumina
concentration is fed back to the control system to determine the amount of
alumina to be fed. The control signal is then passed to both the model and
the process. All available measurable parameters from the process, in this
example, the cell voltage and temperature, are then compared with the
predicted values from the model. Discrepancies will arise between these
parameters due to uncertainties in the control action. In this example,
there is a variation on the amount of alumina fed to the cell. This
discrepancy is used to alter effects of the control action in the model until
the difference between the measured and predicted variables is within the
required tolerance. The simulation is then continued and this procedure is
repeated at the next integration interval.
198
, I'\.
predicted
h• [A1 20 31
. £ U
-, l: controller l: model correction to
mass of
T, V
AI20 3 fed
comparator/
corrector
~
T, V
process
(1 )
199
~
·s
E
~
«
1:)
c:
~
~ .0
r--"'I
M '-
(J)
c:
0C\I Q)
a......
-
« c:
........ W ~
10....:
:J
... ()
-CO c::
0) 0
:§ '-............
--
+-' ~
c: .§ c::
Q)
.....J
Q) (.)
+-' c:
Q)
c:
0
0
a.. ~
......
()
Q)
co
.9-
.c:
(J)
c:
'- 0
......
~
Q)
E
._1 cc:
-..;::-- co
1
1
t-..:
~
:J
.0)
Li:
lUaJJn8
- 'Al!SUaa
I
The alumina concentration meter comprises a probe, that inserts into the
molten electrolyte, and a voltage-current programmer that measures the
voltage drop across the probe for a programmed series of imposed current
steps. The probe is shown in figure 7.9 and consists of graphite anode
and cathode separated with a boron nitride insulator encased in a steel
sheath.
200
+
-
Graphite Anode
and Cathode
Steel Tube
Boron Nitride
Insulator
The rod was inserted back in the anode, and the gap between the rod and
the anode was packed with crushed cryolite. The anode was then loaded
into the cell normally and the connections made to the potentiostat device.
Although installation appeared to have been successful, the trial did not
yield any satisfactory results and the anode was removed early due to
poor performance. It was impossible to determine whether the anode
failure was due to the insertion of the probe or simply due to the anode
itself.
201
CELTROL
Potential
Control
Device
/ / / /
/ / / /
/ / / /
/ / / / / / / / / / / / /
Molten Electrolyte ,.
/ / / / / / / / / / / / / /
Aluminium Pad
References
2 Johnston, T.J. and Richards, N.E., "Laboratory and Field Studies with an
Improved Direct Reading Alumina Concentration Meter", AIME Light
Metals, vol 1, 1974, pp94
202
Appendix A Initialisation of the Simulation Program
At the end of every simulation, the program dumps the model and control
parameters to disk in ASCII format. These 'end ' files can then be used as
the initialisation files for subsequent simulations.
(*This is the end file of a 170 hr run with model Sc. It will be used as
(*an initialisation file for further runs on that model.
**OATA**
(* Model Version
SIXC
(* Cell type
HALF
(* current A
1.5100E+005
(* random number seed absolute variance on feed ± 20%
3 2.0000E-001
(* mass of feed g/dump
2.3000E+004
(* mass of bagged AI F3 added g
2.5000E+004
(* pad depth bath depth acd cm
203
18 8.0000E+001 1.3460E+002 4.0000E+001 2.2S00E+001
S.OOOOE+OOO
1.S600E+000
204
4.0000E-003 6.9000E-003
(* initial masses
(* alumina in suspn, sludge added by anode change aluminium in solution
6.5510E-003 1.6442E+001 2.2000E+004 1.2924E+003
(* dissolution rate parameters kO tempO
1.2450E-003 1.3830E+003
(* specific areas for alumina in sludge and suspension
7.5000E+001 1.5000E+002
(* sludge factor
5.0000E-001
(* Tafel coefficients for oxide reaction ao bo and anode effect af bf aev
4.7250E-004 6.9500E-005 2.5500E+OOO 3.0000E+001 5.0000E+OOO
(* Anode effect bubble parameters thickness,leakage,conductivity
6.1000E-001 1 .OOOOE+OOO 2.0300E-002
(* Aluminium dissolution parameters kl padvelO cm/s kc02diss/kaldiss
3.9030E-003 1.3500E+001 1.0000E+OOO
(* cathodic current efficiency
1.0000E+OOO
(* beam speed, (up/down)
2.6000E-001 4.0000E-001
(* tap rate g/s
5.0000E+003
(* composition of ledge ratio, wt% alumina, wt% CaF2
1.4000E+OOO 4.0000E+OOO 4.0000E+OOO
(* density a1203, ledge, cathode, alumina
4.0000E+OOO 2.5000E+OOO 1.6000E+OOO 1.0000E+OOO
205
(* htc cathode side and bottom
6.1400E-005 7.0000E-005
206
9.6571 E+002 9.5426E+002 9.4280E+002 3.9225E+002
(* anode #7 f, xf
1.8224E-003 4.6010E+000
(* Tf, T1, T2, TTop
5.6439E+002 5.4695E+002 5.2978E+002 2.9026E+002
(* anode #8 f, xf
1.0000E+OOO O.OOOOE+OOO
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* anode #9 f, xf
1.0000E+000 O.OOOOE+OOO
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* anode #10 f, xf
1.0000E+OOO O.OOOOE+OOO
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* anode #11 f, xf
1.0000E+000 O.OOOOE+OOO
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* anode #12 f, xf
1.0000E+000 O.OOOOE+OOO
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* anode #13 f, xf
1.0000E+000 O.OOOOE+OOO
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* anode #14 f, xf
1.0000E+000 O.OOOOE+OOO
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* anode #15 f, x
1.0000E+000 O.OOOOE+OOO f
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* anode #16 f, xf
207
1.0000E+000 O.OOOOE+OOO
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* anode #17 f, xf
1.0000E+000 O.OOOOE+OOO
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* anode #18 f, xf
1.0000E+000 O.OOOOE+OOO
(* Tf, T1, T2, TTop
9.6574E+002 9.5433E+002 9.4292E+002 3.9230E+002
(* output parameters
(* plot frequency tenth minutes
10
(* tap file name
al6c3.tap
(* data file name
al6c3.dat
(* Data file output
ON
(* Tap file output
ON
(* Printer output
OFF
(* Autotap switch
ON
(* tap frequency hrs
24
(* time to next tap hrs
22
(* AIF3 addition frequency hrs
24
(* time to next AIF3 addition hrs
22
208
(* Time for anode out of bath mins
4
(* Simulation Time hours
1.7000E+002
(* feedtimer
25
(* target call em
19
(* pad target em
14
Although the initialisation files allow all the parameters in the model to be
initialised, it is a cumbersome method for dealing with those variables that
most commonly require changing, for example feed times or target
voltage. It also makes the program difficult to operate by anyone without a
very detailed knowledge of the initialisation files. In order to facilitate the
use of the software, a front end menu system was developed to allow the
main operating parameters to be easily changed.
The menu system is accessed by function keys from the top level menu
shown below.
Fl Output Options
F2 Tapping Schedule
F3 Chemical Additions
F4 Control Variables
FS Model Variables
F6 Run
F7 Quit
209
The version of the control algorithms and the cell model linked into the
program is listed at the top of the screen for reference. The functions keys
allow access to the following menus.
The output options menu controls the data that is stored by the program,
the names of the data files used for storage and the length of the
simulation. The frequency at which data is stored is controlled by the Print
Data option which allows the data to be stored as frequently as every
second. This can be used to increase the resolution of simulated events
such as anode effect or to reduce the amount of data stored for longer
runs.
The simulation will also produce a voltage trace with a three minute
interval in the same format as that produced by Anglesey Aluminium
control and datalogging system. This can be directed to a printer or to file.
Output Options
210
A.2.2 Tapping schedule
The parameters pertaining to the tapping of aluminium from the cell can
be specified from this menu.
Target Pad 14 cm
Automatic Tapping ON
Tap Frequency 24 hrs
Target Call 19 cm
The metal tap rate specifies the rate at which the molten aluminium is
removed from the cell into the crucible. The remaining parameters pertain
to the automatic tapping routine. When the routine is enabled the tapping
sequence will be initiated at the specified tap frequency. The amount
removed is determined by the deviation of the pad depth from the target
pad and is based around the target call parameter. This is the amount of
aluminium that should be removed from the cell, expressed as the
required depth in the crucible, if the pad depth is at the target.
211
A.2.3 Chemical additions
The chemical additions menu allows the user to control the amount and
frequency of aluminium fluoride added to the cell to maintain the
composition of the electrolyte. It is intended that this menu could be
extended to include the addition of other compounds, fo r example lithium
fluoride should this be required.
F2 Cell Type
212
A.2.S Current and voltage
This menu allows the line current and initial target voltage to be specified
for the simulation.
Cell Type
Duct End 4
Tap End 4
The basic alumina feed control parameters are also specified from this
menu. The base feed time is the nominal time between feeds but this may
be altered by the control system, for example during search or starve
cycles.
213
module. Should features be added to the model that require user input it is
then a simple matter to include additional menu items the menu system
without losing the general menu structure. Only the model module need
then be altered.
Fl Bath Composition
F2 Cell Geometry
In version SIXC of the aluminium reduction cell model, there are two menu
items. These allow the following parameters to be changed.
The bath composition menu allows the initial composition of the electrolyte
and the bath temperature to be specified. The alumina is specified in three
states, as a weight percentage of dissolved alumina, as a mass of alumina
held as agglomerates in suspension and as a mass of alumina in sludge
pockets around the edge of the cell.
214
A.2.9 Cell geometry
The cell geometry menu allows the initial shape of the frozen sidewalls
containing the electrolyte to be specified. As well as requiring the average
thicknesses of the ledge in the regions of the bath and the metal pad the
depths of the electrolyte and the aluminium pad must be specified.
The final parameters on this menu are the anode cathode distance and
the operating volts at the start of the simulation. These two parameters
are interdependent and the user can specify which should be calculated
by the program.
215
Appendix B Output from the Simulation Program
In addition to continuously recording the conditions with in the simulated
cell in data files on the computer disc, the simulation outputs many of the
parameters to the computer screen. This allows the user to follow the
progress of the simulation and to make changes to some of the
parameters should this be required. The format of the output screen is
shown below.
Day
Cel l Conditions a t 0 : 2 : 38
DAY BATH PAD TAPPED
em em kg
OPERATING VOLTS 4 . 39 V
PAD DEPTH 12 . 4 em
PAD LEDGE 4 . 6 em
BATH LEDGE 8 . 6 em
4 . 8 em
ACD
WT% AL 2 03
3.95 %
216
would be desirable to monitor continuously, for example the alumina
content of the electrolyte, are also displayed.
In the top left hand corner of the screen, the status of the various control
routines is continuously updated. Reading from left to right these are
The ability to change operating conditions during the simulation allows the
user to investigate the behaviour of the cell to disturbances in the system.
For example, use of the manual tapping option would allow the effect of
over or under tapping to be investigated.
217
Appendix C. Graphical Presentation of Results
A necessary part of the simulation of a process is the ability to present the
data in an intelligible format. For dynamic simulation this generally
requires the process variables to be displayed graphically against time.
218
Once the required graphs have been produced and investigated on
screen, hard copy may be obtained on a Hewlett Packard plotter.
219
References
Aarebrot, E. et ai, Metall. vol 32, no 41, 1978
Arita, Y., Urata., N. and Ikeuchi, H., "Estimation of Frozen Bath Shape in
an Aluminum Reduction Cell by Computer Simulation", AIME Light Metals,
1978,pp59
Bagshaw, A. N., Kuschel, G., Taylor, M. P., Tricklebank, S. B., and Welch,
B. J., "Effect of Operating Conditions on the Dissolution of Primary and
Secondary (Reacted) Alumina Powders in Electrolytes", AIME Light
Metals, 1985, pp649-660
Belyaev et ai, Izv, Vyssh. Ucheb. Zav., Tsvetn. Met., vol 4, no 3, 1961,
p67
Berge, B., Grjotheim, K., Krohn, C., Naeumann, R. and T0rklep, K., "The
Influence of Operating Parameters on the Current Efficiency in the
Aluminium Reduction Cells", AIME Light Metals, vol 1, 1976, pp23
220
Blake, S. R., '~ Model for the Variation of Alumina Concentration in
Industrial Hall-Heroult Cells", MSc Thesis, Department of Chemical
Engineering, University of Newcastle upon Tyne, June 1983
Bratland, D., Grjotheim, K., Krohn, C. and Motzfeldt, K., Journal of Metals,
vol 19, 1967, pp13
Bratland, D., Jose del Campo, J. and Cho, K., "Current Efficiency
Measurements in Laboratory Aluminium Cells. Influence of Bath Volume I
Anode Size and Anode Quality", AIME Light Metals, 1981, pp281-308
Calandra, A.J., Castellano, C.E., Ferro, C.M. and Cabo, 0., "Experimental
and Theorectical Analysis of the Anode Effect in Industrial Cells", AIME
Light Metals, 1982, pp345-358
221
Evans, J.W., Zundelevich, Y. and Sharma, D., Met. Trans, vol 128, 1981
Firsanova et ai, Izx. Vyssh. Ucheb. Zav., Tsvetn. Met., vol 5, no 3, 1962,
p52
Frovlova, E.B., Dobrokhotov, V.B. and Tsyplakov, A.M., Trudy V.A.M.I, 89,
vol 36, 1971
Gerlach, J., Hennig, U. and Kern, K., "The Dissolution of Aluminium Oxide
in Cryolite Melts", AIME Light Metals, vol 1, 1974, pp49-61
222
Haupin, W.E., "Overvoltages Determined by Nonlinear Regression
Analysis of Cell Volts vs. Current Density", AIME Light Metals, 1980,
pp183-190
Hertzberg, T., T0rklep, K. and 0ye, H.A., "Viscosity of Molten NaF- AIF3-
AI2 0 3 -CaF2 Mixtures", AIME Light Metals, 1980, pp159-170
Hove, S.J. and Kvande, H., "Centre-Break Alumina Feeding and Sludge
Control of Prebaked Cells", AIME Light Metals, 1982, pp513
Johnston, T.J. and Richards, N.E., "Laboratory and Field Studies with an
Improved Direct Reading Alumina Concentration Meter", AIME Light
Metals, vol 1, 1974, pp94
223
Keller, R., "Alumina Dissolution and Sludge Formation", AIME Light
Metals, 1984, pp513
Lillebuen, B., Ytterdahl, S. A., Huglen, R., and Paulsen, K. A., "Current
Efficiency and Back Reaction in Aluminium Electrolysis", Electrochimica
Acta, vol 25, No, 2, February 1980, p131
Maeda, H., Matsui, S., and Era, A., "Measurement of Dissolution Rate of
Alumina in Cryolite Melt", AIME Light Metals, 1985 pp763-780
Matiasovsky, K., Danek, V. and Malinovsky, M., "Effect of LiF and Li~/F6
on the Electrical Conductivity of Cryolite-Alumina Melts", J. Electrochem.
Soc., 116, 1969, pp1381-1383
224
Newitt, Dombrowski and Knellman, Trans. Inst. Chem. Engr., vol 32, 1954
p244 '
Paulsen, K.A., Huglen, R., Andersen, J.A., Lillebuen, B., and Ytterdahl ,
S.A., "Variation of Side Lining Temperature, Anode Position and
CurrenWoltage Load in Aluminium Reduction Cells", AIME Light Metals,
1980, pp325-339
Qiu, Z. and Fan, L., "The Rate Determining Step of Metal Loss in Cryolite-
Alumina Melts", AIME Light Metals, 1984, pp789
Saksvikronnig, T., Valsvik, G., and Hove, S. J., "Anode Effect Quenching
by Anode Tilting", AIME Light Metals, 1982, p553-558
225
Schmitt, H., "Extractive Metallurgy of Aluminium", vol2 Intersc' P b N
' I. U., ew
York, 1963, p169
Sides, P.J. and Tobias, C.W., J. Electrochem. Soc., vol 129, 1982, p2715
Sidorov, L.N., Kolosov, E.N. and Shol'ts, V.B., Russian Journal of Physical
Chemistry, 42, pp1382-1386, 1968
Szeker, Acta Tech. Acad. Sci. Hung., vol 10, 1955, p91
226
:honst~d,. J. and ROlset~, S., "Equilibrium between Bath and Side Ledge
In AlumInIum Cel/s. BasIc Principles", AIME Light Metals, 1983, pp415
Winkhaus, G., "On the Dissolution Rate of Alumina in Cryolite Melts", 99th
AIME Annual Meeting, Denver,1970, Paper A70-25
227
Zuca, S., Herdlicka, C. and Terzi, M., "On Porosity-Overvoltage Correlation
for Carbon Anodes in Cryolite-Alumina Melts", Electrochimica Acta, vol 25,
February 1980, pp211-216
228
Nomenclature
aJ activity of species j in the electrolyte
A area
ACO anode to cathode distance
Cj concentration of species j in the electrolyte
c~ saturation concentration of species j
J
229
total current
I·J current associated with reaction j
current density
l'
J
partial current density associated with reaction j
K electrical conductivity
k rate constant
k·J thermal conductivity of component j
Kd equilibrium constant for dimerisation of NaAIF4 vapour
kL mass transfer coefficient
Kp constant of proportionality
kt turbulent kinetic energy
L physical dimension for description of system
M·J mass of component j
n number of anodes
N stirrer speed
Vf velocity of anode gases
Nu Nusselt number
Pj partial pressure of species j
Pb barometric pressure
PE scrubber efficiency for particulates
Pr Prantl number
Pv total vapour pressure
q heat flow
Ok volumetric flowrate of stream k
R·I rate of reaction i
R·J electrical resistance of component j
R universal gas constant
Ra Rayleigh number
rAI 20 3 radius of AI 20 3 agglomerate
rate of convective transport of aluminium into the
rAlconv
electrolyte
aluminium dissolution rate
rAldiss
Rb weight ratio NaF/AIF3
Re Reynolds number
230
SCt turbulent Schmidt number
Sf fraction of feed becoming sludge
Sh Sherwood number
T temperature
To reference temperature
ta thickness of bubble layer at anode
td dead time
tf fixed feed time
tremelt time taken for electrolyte frozen onto a surface to melt
ts activation time for solenoids
U overall heat transfer coefficient
Ut turbulent viscosity
V volume
W·J molecular weight of species j
wt weight percent alumina in electrolyte
x thickness of a layer or zone
Xj mole fraction of species j
xa anode width
XA1F3 AIF3 in electrolyte expressed as a weight percentage
excess of that in pure cryolite
anode length
distance normal to the interface
Charge number
Tafel constant
Tafel constant
bubble diameter
gas fraction in free bubble layer
emissivity in centre trench
11 overpotential
f inter electrode gap
A'J latent heat of fusion of species j
Il viscosity
v interfacial velocity
e surface coverage by bubbles at the anode
p density
Stefan Boltzmann constant
molar volume of species j
231
suffixes
A Anode
ac alumina cover
af frozen layer at anode surface
AI Aluminium
AI 20 3 Alumina
A13+ Aluminium ions
AIF3 Aluminium fluoride
As anode surface
b bath
ba bath
bf bath / frozen ledge interface
bt bath / alumina interface
C cathode
cav cavity beneath shields
CO2 carbon dioxide bubbles
conv convective
cov alumina cover on top of anode
ct centre trench
d dimer Na2AI2 Fa
diss dissolution
e electrolyte
ex exposed bath
ext external
F fluoride ions
f freeze
g gas
h horizontal component
H2O water
HF Hydrogen fluoride gas
hole hole in the crust
I liquid
m monomer NaAIF4
mf metal pad / frozen ledge interface
NaAIF4 Sodium Aluminium Fluoride
0 oxide reaction or ions
sludge alumina as sludge at the cathode
stub anode stub
232
suspn alumina in suspension in the electrolyte
tc centre trench / cavity interface
v vertical component
vol volatilisation
W water
sat saturation
233