Lab Report Group 5
Lab Report Group 5
Group No.: 5
Laboratory No.: 1,2,4,5 and 6
All forms of plagiarism, cheating and unauthorized collusion are regarded seriously by the University
and could result in penalties including failure in the Unit and possible exclusion from the University.
Please read the Curtin University Statement on Plagiarism. If you have any doubts, please contact the
Unit Co‐ordinator.
Declaration
Except as indicated, this work being submitted is solely my/our own work.
Name ID Signature
1
Experiment 1 – Leaching of Chalcopyrite
1. Introduction
1.1 Theory of chalcopyrite leaching
The sulphide minerals of copper such as chalcopyrite (CuFeS 2) are not readily
dissolved in dilute acid. Their solubility will increase in the presence of oxidizing
agents such as oxygen, ferric ions, and bacteria. However, acid medium such as sulfuric
acid lack selective properties in the leaching process. Hence, some gangue materials
also dissolved in the liquid medium, necessitating an expensive purification process [1].
Using an alkaline lixiviant seems to have many advantages given the known difficulties
with acidic leaching. Utilization of organic alkaline lixiviant (monosodium glutamate)
to selectively dissolved copper from electric arc furnace dust had been investigated,
with approximately 95% recovery of copper. The use of monosodium glutamate (MSG)
as a reagent in copper alkaline leaching has some advantages such as its wide
availability and biodegradability [2].
Leaching process is considered as a heterogeneous process which involves solid and
liquid phase. The kinetic of leaching process is essential in determining reactor types
and the size of the reactor. Rate of reaction and mean residence time required to
completely dissolve metal in a liquid medium were needed to measure kinetic
properties of leaching processes [3]. The leaching process frequently takes place in
three stages in many hydrometallurgical systems: (a) the mass transfer of reactants and
soluble product from bulk of liquid to the surface of solid; (b) the diffusion of reactant
and soluble product through the pore of solids; and (c) the chemical reaction of reactant
in the solid-liquid surface. The slowest step will be the rate-limiting step for a
particular set of experimental conditions [3].
One of the most reliable models to describe the kinetic of the leaching process is the
shrinking core model. This model states that during leaching, the reaction takes place in
the outer skin and then continuously happens in the deeper zone where an unreacted
core exists. The reacted core will form an ash layer that is assumed to be an inert layer
[4]. The Shrinking Core model covered three mechanisms of rate control during the
leaching process which are diffusion through liquid film control, diffusion through ash
product control, and chemical reaction control. Three types of rate controls covered in
the Shrinking Core model can be summarized in these following equations:
1) When diffusion of reactant A through the liquid boundary layer is rate limiting:
3bkg Cal t
X b= =K f t
ρb R
2) When diffusion of reactant A through the ash (reacted) layer is rate limiting:
2/ 3 6 b D e Cb t
⌊ 1−3(1−X b ) +2(1− X b) ⌋=K d t= 2
ρb R
3) When chemical reaction at the surface of the shrinking core is rate limiting:
2
1/ 3 b K s Cb t
⌊ 1−( 1−X b ) ⌋ =K s t=
ρb R
1.2 Objective
The primary objective of this experiment is to become familiar with chemical
laboratory techniques and with the interpretation of kinetic leaching data.
2. Materials and methods
2.1 Materials
In this experiment, chalcopyrite (CuFeS2) was used as a solid specimen which was
pulverized to obtain -200μm in size. Chalcopyrite was assumed as the only copper
mineral in this specimen. Monosodium glutamate was used as the leaching agent, and
hydrochloric acid (HCl) as pH regulator. Merck brand analytical grade chemicals were
used, and deionized water was used in the preparation of analytical samples.
2.2 Methods
In the preparation of the leaching solution for each test, a 1L volume was achieved by
introducing 600ml of distilled H2O into the 1L volumetric flask. Prior to the addition of
the required quantity of Monosodium Glutamate, 5ml of concentrated HCl was
carefully incorporated. The required amount of monosodium glutamate was dissolved
in distilled water to achieve 0.5 M and 0.25 M reagent concentration. The pH solution
was evaluated within the range of 0.5 to 0.7 by adding concentrated HCl incrementally,
with comprehensive mixing during pH measurements to ensure consistent pH levels
across all solutions.
The solution was then transferred to the reaction flask, placed within a water bath, and
equipped with a thermometer and stirrer. While maintaining stirring, the solution was
gradually brought to the requisite temperature, with particular emphasis on maintaining
a steady temperature reading. In this experiment, the temperature varied between room
temperature and 520C.
C 0.5 M 520C
D 0.25 M 250C
In preparation for data collection, a record sheet was established for time and sample
number documentation, and a stopwatch was initialized. A precisely weighed 1-gram
sample of chalcopyrite was introduced into the reactor, and timing commenced
immediately. Periodically, the stirring rate was adjusted to prevent CuFeS 2 from
floating on the solution's surface. Sample extraction began with the removal of
approximately 6 to 7 ml of the solution using a 10ml pipette, which was then
transferred to a designated sample vial marked as "sample zero."
A total of twelve additional samples, each of the same size, were obtained at specified
elapsed time intervals of 2, 4, 6, 10, 15, 20, 30, 40, 50, 60, 70, and 80 minutes. To
3
separate the solid CuFeS2 from the solution, these samples underwent filtration,
necessitating the use of a clean, dry 10mL pipette, filter paper, and a filter funnel.
3. Results and Discussions (between 3 and 6 pages)
3.1 Fraction of chalcopyrite converted.
The fraction of chalcopyrite converted during specific time was calculated using
equation below:
wt of Cu∈solution
X b= × 100 %
wt of Cu∈initial sample
Initial weight of Cu was calculated by assuming that the sample was pure chalcopyrite.
The weight percentage of Cu in Chalcopyrite was obtained by dividing atomic weight
of Cu and molecular weight of Chalcopyrite (CuFeS 2). The mass fraction of
chalcopyrite converted in this experiment were tabulated in Tables 3 and 4 below.
Wt sample 1 gram
% Cu: 34.626
Wt of Cu 0.34626 gram
346.26 mg
Volume: 1 litre
Table 2 Calculation of Cu weight percentage in the sample
4
Figure 1 Mass conversion of copper with variation of temperature
The results in figure 1 indicate that leaching in room temperature has a higher
conversion rate than 500C. After 300 minutes of leaching time, the amounts of copper
being leached were 0.503% and 0.179% for experiment C and experiment D,
respectively. The possible reason for these results is because experiment C were tested
using more concentrated reagent and higher temperature. Another study stated that the
increase in temperature strongly affects the percentage of copper dissolved in
monosodium glutamate (MSG) reagent. While variation in MSG concentrations had
less significant effects in percent copper extraction [2].
5
3.2 Leaching kinetics
To describe chalcopyrite leaching kinetics, two kinetic models were used to determine
the rate limiting process: diffusion through ash product layer or chemical reaction
process. The plotting results of experimental data to those kinetics models were listed
in Table 5.
Table 5 Correlation coefficients for experimental data fitting using two kinetics models.
Based on the results tabulated in Table 5, diffusion through the product layer was
determined as a rate limiting process. The value of coefficient of correlation (R 2) of the
ash diffusion control model for both experiments was higher than the chemical reaction
control model. It means that the rate of reaction was determined by the rate of diffusion
of reagent through the product layer.
6
Figure 3 Fitting of ash diffusion control model in experiment C and D
The determination of activation energy (Ea) for the leaching process was carried out
using the Arrhenius equation K = Ae (-Ea/RT), where K represents the rate constant through
the unreacted solid layer, A is a surface area constant, T is the leaching temperature (in
Kelvin), and R is the ideal gas constant (8.314 J/K.mol).
The value of K was determined as the slope of the rate limiting step model (ash
diffusion model). If the same step is the limiting mechanism in examined temperature,
then the value of K will be varied following linear regression trend. The value of K for
both different temperatures were estimated using mean square error (MSE) methods.
The results were listed in Table 6.
The value of Ea can be obtained by plotting ln K against (1/T), which yields a straight-
line equation with a slope, as demonstrated in Figure 4. From the slope of the curve in
Figure 4, the activation energy value for the leaching of Chalcopyrite in Monosodium
Glutamate solution was determined to be 60.72kJ.
7
Figure 4 Plotting of Ln K vs 1/T
In the context of high values of fraction converted, it is commonly observed that the
reaction rate exhibits a tendency to plateau, and this phenomenon can be attributed to
several underlying factors. One prominent factor is the depletion of reactants as the
reaction progresses. As the reaction proceeds, the concentration of reactants diminishes,
resulting in a diminishing rate of reaction. At high levels of fraction converted, the
concentration of reactants reaches exceedingly low levels, contributing to the flattening
of the reaction rate curve. Furthermore, equilibrium limitations can significantly impact
reaction kinetics under such conditions. A reaction may approach equilibrium at high
levels of fraction converted, leading to a diminishing reaction rate. These dual factors,
namely reactant depletion and equilibrium limitations, play pivotal roles in elucidating
the observed flattening of reaction rates at elevated values of fraction converted.
8
3.5 Effect of larger particle size in limiting mechanism
As the particle size increases, the diffusion of reactants to the surface of the particle
becomes slower, leading to a decrease in the rate of reaction. Chemical reaction control
is independent of particle size. Based on the shrinking core model equation, the ash
diffusion control mechanism shows the longest time needed to achieve the same
fractional conversion. It equals the square of particle radii. Those the bigger particles
size will lead to longest distance for both reactant and product to diffuse and react with
the unreacted core [4].
3.6 Practical particle size limit
In practice, the optimum particle size for leaching is often determined through a
balance of diffusion limitation, energy consumption, reagent consumption, and solid-
liquid separation. Smaller particles typically have better mass transfer characteristics,
which means that leaching agents can penetrate and react with the material more
rapidly. However, extremely fine particles can suffer from reduced permeability,
leading to slower leaching rates and incomplete extraction. The energy required to
comminute (crush or grind) materials into smaller particles can increase significantly as
particles become finer. Extremely fine particles may require more leaching reagents per
unit of mass due to the increased surface area. After leaching, the solid-liquid
separation process becomes more challenging as particle size decreases. Extremely fine
particles can form sludges or suspensions that are difficult to clarify or filter. It is
essential to conduct thorough laboratory testing and modelling to identify the size range
that maximizes leaching efficiency while considering practical limitations. This may
involve trade-offs between faster kinetics associated with smaller particles and the
operational challenges they present.
3.7 Optimum particle size
One way to determine the optimum size for the leaching process is to conduct
experiments with different particle sizes and observe the leaching efficiency. For
example, in a study on copper recovery from waste printed circuit boards, the optimum
particle size was found to be 150 mesh according to the leaching efficiency [6]. The
smaller the particle size, it has a bigger surface area that leads to an increase in leaching
kinetic. However, the effect of particle size is also correlated with the precedent
processing stage. To get smaller particle size, grinding circuits need a higher energy
consumption per amount of material produced. Hence, there is a trade-off between
leaching recovery and operational cost for grinding circuits.
4. Conclusions
The fraction of chalcopyrite converted during the leaching process was calculated and
presented in Tables 3 and 4. It is evident that leaching at higher temperature (52°C)
exhibited a higher conversion rate compared to a room temperature (25°C). After 300
minutes of leaching, the conversion rates were 0.503% and 0.179% for experiments C
and D, respectively. This variation is attributed to differences in reagent concentration
and temperature. The linear regression analysis results indicate that the rate control of
the reaction is diffusion through the solid product layer, even though the value of
activation energy was 60.72kJ/mol that was probably caused by unreacted minerals.
9
10
EXPERIMENT 2- CURRENT EFFICIENCY FOR ZINC
ELECTRODEPOSITION
1. Introduction.
1.1. Theory:
Electrode reactions are mainly redox reactions forced by a power source (electricity). A
net flow of charge travel into the electrode or out the electrode where reduction and
oxidation respectively take place. A change in the electrochemical potential occur when
the reaction happens, from an initial equilibrium Eeq to a new one Ereaction . this change
can be represented by the term polarisation potential or overpotential and given by
equation 1 below, and the net current that flows is the difference between the anodic i a
and the cathodic i c current as shown in equation 2:
i=i c −i a (2)
There are two types of polarisation potentials based on the cause of occurrence:
a- Activation polarisation: the reaction is chemical controlled and the change is due to
the activation energy for the electrode process.
b- Concentration polarisation: the reaction is transport controlled and the change is
due to a decrease in the concentration of reactant ions in the vicinity of the
electrode.
The relationship between the net current density j ( mA ) and the overpotential can be
2
η=a−b ln j(3)
Where a and b are Tafel constants, unique values for each element redox reaction.
For the theory of the zinc electrodeposition, hydrogen overpotential on the surface of
the cathode can happen. Many factors can affect the hydrogen overpotential such as pH,
current density, and temperature. The current efficiency of electrodeposition has been
shown to be greatly affected by small concentrations of certain impurity ions in the acid
sulfate baths.
The presence of impurities on the cell can increase the hydrogen evolution by:
b- Acting as galvanic couples with zinc, thus accelerating the dissolution of the zinc
electrodeposit with subsequent hydrogen evolution.
11
1.2. Objectives
a- Determining the current efficiency for the electrodeposition of zinc from a pure
zinc sulfate solution at the zinc cathode.
Equipment:
5ml burette, Ammeter, 250ml beaker, Stopwatch, lead anode, Electrical connecting
wire, 7 x zinc cathodes, Switch Voltmeter and Rheostat
Solutions:
Zinc sulfate solution (50g/L Zn and 1M H2SO4)
Cadmium chloride solution (0.1M)
Copper sulfate solution (0.1M)
Stannous chloride solution (0.1M)
2.2. Methods:
Current Efficiency:
A 200ml of the pure zinc sulfate solution was added to 250 ml beaker. Four zinc
cathodes were polished, washed with water and acetone, dried and weighed. The
dimension of each zinc electrode was measured to calculate the area and hence the
current required to give a current density of 400 A /m2.
The reaction was initiated by closing the circuit and inserting the cathode, and the
clock was started. Using the rheostat, the current density was adjusted to 400 A/m 2 and
maintained the current at the set value throughout the experiment. The cell voltage
should be about 3.5V but will vary according to solution resistance, temperature,
distance between electrodes, etc.
After 20 minutes the reaction was stopped. The electrodes were washed with water and
acetone and weighed again when dried.
Effect of impurities:
0.1M Copper and 0.1M Tin impurities are supplied. The same method has been used to
Determine the current efficiency for zinc sulfate solutions containing firstly 2mL then
5mL of one of these impurities using a fresh zinc cathode.
12
3. Results and Discussion.
Z n2 +¿+2 e ¿
Then:
65.38∗0.4∗1200
W= =0.1626
2∗96485
0.1636
≅1
0.1626
The efficiencies for the other five samples are shown in table below:
13
It can be noticed that the presence of impurities generally decreases the current
efficiency of the cell. The hydrogen gas coming out from the electrode was observed
and has increased as impurities presented in the cell. However, Copper didn’t have a
great impact on the efficiency dropping. The suggested reason is that for the efficiency
dropping to happen, the zinc should dissolve in into the solution, and at a certain
concentration the copper deposit would be enough to create a local galvanic cell.
The relation between the concentration of the impurities and the efficiency was also
studied. Figure 1 shows how does the concentration effects the efficiency.
Increasing the concentration of the impurities has effect of the efficiency of the cell, the
direct proportion suggested that increasing the impurities ions interfere with the zinc
ions.
The hydrogen overpotential was calculated for the impurities Cu and Sn using equation
3 and the results are shown in table 2 below:
a b η
4. Conclusion.
As Kerby and Ingraham suggested [7], the results show that the current efficiency is
affected by the concentration of impurities. The more concentrated the impurity ions
14
are, the more drop the efficiency will have. They have classified the impurities to four
main classes, Tin is in class three, which has high hydrogen overvoltage and tend to
alter the crystalline structure of zinc deposit. Copper on the other hand is in class four,
which has low hydrogen potential and small effect on the cell efficiency.
EXPERIMENT 4 – LEACHING OF COPPER IN AN AMMONIACAL
SOLUTION
1. Introduction
1.1 Theory
Ammonia leaching has been widely used in hydrometallurgy for processing of
sulphide and oxide ores, concentrates or by-products, this solvent can be
effectively used for leaching of base metals (Copper, Niquel, Cobalt, Zinc) as well
as for precious metals (Gold, Silver). Additionally. in the case of materials
containing copper-iron and iron minerals (chalcopyrite, bornite, pyrite, pyrrhotite)
it could be beneficial because of its selectivity, which means that allows the
selection of the desired metals to be solubilized while undesired iron can be
precipitated, everything in just one operation. Beside this, it is considered a useful
leaching agent due to several factor such as low grade of toxicity, accessible cost,
and easy regeneration process from alkaline solution [8].
The experiment is going to be addressed with cooper and an ammonia reagent for
the leaching of copper, this base metal (copper) is soluble in an ammoniacal
solution according to the reaction:
−¿¿
2+¿+4 O H ¿
2 Cu+ xN H 3+O 2+2 H 2 O=2 Cu ( NH 3 ) x
2
1.2 Objectives
The aim of this experiment is examining how the reaction is influenced by
variations in the concentration of ammonia and adjustments in stirring speed and
temperature. The ammonium ion concentration will remain constant.
The lab for this group 5 is composed of developing the second analysis based on
effect of stirring speed on leaching rate.
2. Procedure
Effect of Stirring Speed on Leaching Rate Experiment
A copper disc with a diameter of 3.7 cm is going to be used in the experiment, it is
crucial to clean the copper coin by abrading and polishing one or both surfaces, for
this experiment is used a sandpaper and both sides were cleaned. Beside this, the
15
copper coins were placed into two different leaching vessels where the
concentration of ammonia to study was 0.05 M and 0.01 M.
3. Calculation
Solution needed for each concentration of ammonia.
C1V1 = C2V2
- (0.05) x (800mL) = 1(V)
- 40 mL NH3
- (0.01) x (800mL) = 1(V)
16
- 8 mL NH3
Area of the copper coin.
During the experiment was not possible to set the stirring speed to 30 rpm, the
stirring speed used is specified in the table 2.
17
4. Results
To determine the leaching rate in g Cu/cm 2/second for each experiment it is
essential to know the results of the copper concentration that is shows in table 3.
It is possible to notice in the figure 2 the concentration of copper that was leached.
It is notorious to evidence the increment of copper leached when the stirring rate
increased, with 38.2 rpm was a total of 43.4, 83.2 rpm was a total of 74.4 and 119
rpm was a total of 114.7 the leaching rate difference between the starting point and
ending point.
18
Figure 2 Leaching Rate 0.01 Ammonia concentration.
Ti
m
Cu Leachi
Va Sa r e
(m ng
ss mpl p (
g/L Rate *
el e m m
) 10^9
in
)
2 3
0.0 4 1 0.0
5 . 0 3
5 4.7
3
0.0 4 2 0.0
5 . 0 7
5 10.9
3
0.0 4 3 0.0
5 . 0 8
5 12.4
8
0.0 2 4 0.4
5 . 0 3
1 66.7
8
0.0 2 5 0.5
5 . 0 5
1 85.3
0.0 8 6 0.6 99.2
5 2 0 4
19
.
1
1
0.0 7 0.7
1
5 0 4
7 114.7
1
0.0 8 0.8
1
5 0 4
7 130.2
1
0.0 9 0.9
1
5 0 5
7 147.3
Table 4 Experiment 0.05 concentration result
As before, figure 4 shows the concentration of copper that was leached. It is
possible to see the increasing concentration of copper leached when the stirring
rate was settled with to high values such as 34.5 rpm was a total of 12.4, this result
is possibly wrong because it is not going with the trend, theoretically is around
18.6 according to the slop equation where m = 0.62, the 82.1 rpm was a total of
32.5 and 117 rpm was a total of 32.6 the leaching rate difference between the
starting point and ending point. Therefore, it is possible to highlight that apart from
34.5 rpm the trend of the graph was lineal and the stirring as a transport control is
being effective.
5. Conclusion
These linear results graphed suggests that increasing the stirring rate consistently
and proportionally accelerates the rate of transport of the leaching solution, which
consequently affects the amount of substance leached. Therefore, the stirring rate
has a significant impact on the leaching process under these conditions, indicating
transport control as a dominant factor.
20
21
EXPERIMENT 5- SOLVENT EXTRACTION OF COPPER
1. Introduction
1.1 Theory of Solvent Extraction
One of the key options for producing copper from low-grade ore is the solvent
extraction process. Solvent extraction is a technique for selectively removing or
extracting a species from one liquid phase into another. This technique can be used to
separate valuable species (such as copper) from impurities. To minimise excessive
solvent loss, these two solutions must be insoluble and immiscible [9].
The solution from which the extraction is made in hydrometallurgical copper processes
is a dilute acid solution. In copper processing, sulphuric acid is usually used as dilution
[10]. Copper is recovered by mixing a copper sulphate solution with an immiscible
organic solution. In this process, copper ions will be transferred to the organic phase,
thus the organic phase is called "loaded organic". The barren aqueous phase, called
raffinate, is reused as lixiviant. Next, the loaded organic is stripped to separate the
copper ion from the organic phase to the copper solution aqueous phase. The extraction
and stripping process can also be explained by the reaction 1 below [11].
Extraction
[2 RH](org.) + [CuSO4](aq.) [R2Cu](org.) + (2H+SO42-](aq.) (1)
Stripping
1.2 Objective
The main objective of this experiment is to familiarize with the solvent extraction
process on a laboratory scale and the data analysis.
22
2. Materials and methods
2.1 Materials
In this laboratory experiment, three different solutions were used as a feed solution.
The first solution, called solution A, is composed of 1 g/L Cu 2+, 3 g/L Fe 3+, 0.5 M
Na2SO4 adjusted to pH 2.2 with H2SO4. The second solution is solution B, which
consists of 8g/L Cu2+, 0.5M Na2SO4, and natural pH. The last solution is solution C
which contains 20g/l Cu2+, 150g/L H2SO4, 0.5M Na2SO4. An extractant solution is used
with 10% volume of the complexing agent in an organic solvent.
2.2 Methods
There are two main experiments: batch loading experiment and stripping loading
experiment.
2.2.1. Experiments with Batch Loading
30cm3 of aqueous feed solution A was added to five beaker glasses to mix it with
organic extractant with different volumes. 60 cm3 of organic extractant was added to the
first beaker, 30 cm3 to the second, and 15, 7.5, and 3.75 cm3 of organic extractant was
added to the rest of the beaker glass. All mixed solutions were stirred for 15 minutes
with a magnetic stirrer. The aqueous solution from each beaker was sampled after
stirring. The sample was used to determine the copper concentration of the aqueous
phase using Atomic Absorption Spectroscopy (AAS).
All experiments were done at room temperature (28oC). The result of the experiment
was analysed to construct an equilibrium loading isotherm and equilibrium stripping
isotherm. Additionally, McCabe-Thiele Diagrams were also applied to equilibrium
isotherm to determine the number of stages required to achieve a specified copper
separation between the organic and aqueous phases. An equilibrium loading isotherm
for copper is shown in figure 2. Drawing an operational line with slope A/O passing
through the point (xj,o) produces a McCabe-Thiele diagram for extraction, where xj is
the final aqueous solution concentration required. The figure is then completed by
drawing horizontal and vertical lines alternating between the operating and equilibrium
lines, starting with the first horizontal line drawn through (xo, y1). Each triangular step,
1, 2, symbolises a stage of extraction [12].
23
Figure 2. McCabe-Thiele Diagrams for copper extraction (Dyer, 2023)
The pH of raffinate decreases as copper concentrate decreases because the higher the
extracted copper to the organic phase, the higher the H+ ion will produced, which
follows reaction 1. To obtain the copper concentration in the organic phase, the mass
balance method is used by subtracting the copper mass in the feed A from the aqueous
phase. The calculation is shown below:
Determine the copper mass in the feed solution A:
Copper concentration in feed solution A = 1.133 g/L
Volume of solution A = 30 mL
Mass of copper in feed solution A = 1.133 g/L x 0.03 L = 0.0339 gram
24
Determine the copper mass in aqueous solution A:
Volume = 3.75 mL
Copper concentration (after accounting dilution factor) = 100 * 8.27 = 827 mg/L =
0.827 g/L
Copper mass in 3.75 mL solution A = 0.827 g/L* 0.00375 L = 0.02481 gram
Copper mass in that organic phase = 0.0339 gram – 0.02481 gram = 0.00918 gram
Therefore, copper concentrations in organic phase = 0.00918 gram / 0.00375 L = 2.448
g/L
The rest of calculation is done with same method, thus overall result can be shown in
Table 3. From this table, equilibrium loading isotherm can be constructed as illustrated
in Figure 3.
12
10
Cu Concentration (g/L) organic
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Cu Concentration (g/L) aqueous
Results of this experiment show that the higher the volume of organic added to a
solution, the higher the copper concentration extracted from the aqueous solution.
99.5763 % of copper was extracted in 60 mL of organic extractant, while in 3.75 mL,
only 27% of copper was extracted. This result is happening because the larger volume
of organic extractant will contain a higher volume of extracted copper concentrate from
25
the feed solution. Nevertheless, a higher volume of extractant will cost more for the
operation expenditure. Therefore, the optimum volume of extractant is needed to have a
high extraction value and low cost.
3.2 Equilibrium Stripping Isotherm
In the stripping experiment, copper concentration in the aqueous solution C, feed
solution C, feed solution B, and loaded solution B was identified by AAS, shown in
Table 4.
To obtain the copper concentration in the organic phase in the stripping experiment,
mass balance methods are used. The calculation will be shown below:
Determine the copper mass in feed solution B:
Volume = 200 mL
Copper mass = 1,840 mg/L x 0.2 L = 368 mg
26
Volume = 40 mL
Copper concentration = 22,548 mg/L / 1000 = 22.54 g/L
Copper mass = 22.54 g/L x 0.04 L = 0.90192 g
Initial copper mass = 20,900 mg/L x 0.04 L / 1000 = 0.836 g
Copper mass in organic phase = (1.655 g/L x 0.04 L) – (0.90192 – 0.836) = 0.00028 g
Volume organic phase = 40 mL
Copper concentration in organic phase = 0.00028 g / 0.04 L = 0.007 g/L
The rest of calculation is done with same method, thus overall result can be shown in
Table 5. From this table, an equilibrium stripping isotherm can be constructed, as
illustrated in Figure 4.
30
25
Cu Concentration (g/L) aqueous
20
15
10
0
0 0.5 1 1.5 2 2.5
Cu Concentration (g/L) organic
27
3.3. Determination of A/O Ratio Counter Current
To determine the A/O ratio required to reduce 1 g/L Cu 2+ feed to 0.04g/l in three stages,
assuming 0.02g/L Cu2+ in the organic feed (counter-current), McCabe – Thiele diagram
is applied as depicted in Figure 5. Inferred from Figure 6, the aqueous and organic
phases flow in opposite directions. This means that as the aqueous phase (often the feed
solution) flows in one direction, the organic phase flows in the opposite direction. This
mechanism provides maximum contact between the two phases, which enhances mass
transfer and extraction efficiency [13]. However, its design is more complex than the
crossflow system.
The slope of the operating line is determined to obtain three stages. Then, the A/O ratio
is a slope that can be determined by the gradient formula.
y 2− y 1 2.4−0.02
m= m= = 2.38
x 2−x 1 1−0
It is found that the slope or A/O ratio to obtain three stages is 2.38.
3
2.5
1
Cu Concentration (g/L) organic
3 0.5
0
-0.3 0.2 0.7 1.2
Cu Concentration (g/L) aqueous
28
Where A and B are feed solution, C is organic extractant, R is raffinate, and E is loaded
extractant.
2.5
Cu Concentration (g/L) organic
1.5
0.5
1 A/O = 0.5
0
-0.3 0.2 0.7 1.2
Cu Concentration (g/L) aqueous
To determine the A/O ratio required to reduce 1 g/L Cu 2+ feed to 0.04 g/l in three
stages, assuming 0.02g/L Cu2+ in the organic feed with cross flow system, McCabe –
Thiele diagram is utilized as shown in Figure 8. In a crossflow system, the aqueous and
organic phases move in the same direction, parallel to each other, as seen in Figure 9.
Crossflow systems often have lower extraction efficiency because of the shorter contact
period between the phases than counter-current flow systems. This method stands out
for its simplicity and cheaper equipment costs than counter-current flow systems [13].
The slope of the operating line is determined to obtain three stages. Then, the A/O ratio
is a slope that can be determined by the gradient formula.
y 2− y 1 2.2−0
m=
x 2−x 1
m=
0.46−1
= -4.074
29
3
2.5
1.5
1 1 A/O = - 4.074
2
0.5
3
0
-0.3 0.2 0.7 1.2
Cu Concentration (g/L) aqueous
30
EXPERIMENT 6- RESIDENCE TIME AND NON-IDEAL FLOW
1. Introduction
1.1. Theory:
Space time τ is the time required to process one reactor volume. In batch reactors, all
the fluid elements spend the same amount of time in the reactor, and if the density is
constant through the reaction, the space time will be the same as the mean residence
time.
τ =t=t
Same concept can be applied to the plug flow reactor where there is no density charge.
The residence time can be given by equation 1 below:
X Af
d XA
τ =t=C A 0 ∫ (1)
0
−r A
In the back-mix reactors, some of the fluid elements can pass through the reactor
quickly and others can remain in the reactor for long time. In this case the residence
time of the reactor is hard to be determined. However, residence time distribution RTD
can be produced. The mean residence time t can be equated to τ if the density in
constant.
V
τ =t= (2)
v0
There are two main types of stimulus-response or tracer techniques (step and pulse),
they can be used to determine the mixing characteristic of a non-ideal reactor. A
stimulus (tracer) is input to the vessel and the response is measured in the exit stream
over time. Usually, the tracer is a radioactive isotope or an inert chemical.
Inert A
vo
Inert A
vo
31
32
Rate that A exits the reactor:
dA N
=−C A v o =− A v o
dt V -ve indicates leaving reactor
A t
vo
∫ dA =−∫ dt
A NA
o
0 V
v0
Ln N A - Ln N Ao=− t
V
NA
=e
( )
−v0 t
V
N Ao
V
t=
ν0
(3)
∴
NA
=e
( − )
t
t
or
CA
=e
( − )
t
t
N Ao C Ao
The ratio of CA to CAo in the exit stream can then be calculated at any time t.
Note that CAo is the initial concentration in the reactor, not the concentration of the
tracer when it is being injected.
At t=t
CA −1
=e =0.368 C A =0.368 C A 0
CA 0
That means, 36.8% remains in the tank and 63.2% has exited the tank. The exist age
distribution can be represented by a graph using equation 3 as shown in figure 1 below:
33
Figure 3 Relative concentration of A in the reactor as a function of mean time.
It can be noticed from the concentration of A in the CSTR, particles start to leave right
after the injection. This would cause a problem in surface contact processes like
leaching because the particles don’t have enough time to react with the solution.
The Erlang equation gives a representation of the shape of the distribution curve for
multiple CSTRs in series.
CA n
n
n−1 −n t
RTD= = t e r
C A 0 ( n−1 ) ! r
Where:
n: number of reactors in series.
t r: reduced time t /t .
The stimulus-response studies can be used to study the efficiency of the reactors. six
scenarios have been introduced by Levenspiel to diagnose the faulty flow and suggest
causes. Figure 2 shows these suggestions.
34
1.2. Objective:
2.1. Materials:
Equipment
2 continuously stirred tank reactors, Peristaltic pump, 2 x 100 mL measuring cylinders,
2x 1000 mL measuring cylinders, 50 mL injection syringe, 5 x 100 mL volumetric
flasks, 10 x 100mL beakers or sample vials, Water storage tank, pH meter and probe, 1
mL and 10 mL pipettes and filler
Chemicals
100 mL of 5 mol/L acetic acid.
Calibration standards for pH meter
2.2. Methods:
The whole experiment has been conducted using the tap water.
A 100 ml of 5 mol/L acetic acid was used as a tracer, 5 samples of different dilution
factors were prepared, and pH was measured for each to get a calibration curve as
follow:
35
Single Reactor:
Using a single reactor, the volume of the reactor was measured and filled with tap water
by using a measuring cylinder. The Agitation has been started and pump power was
adjusted to 200 mL/min.
The exact volumetric flow rate has been measured to calculate the residence time using
equation 2. At time zero a 20 mL of 5 mol/L acetic acid was injected into the reactor at
the base of the reactor under the impellor.
A 10 ml of samples were collected at 0.1 times the mean residence time for one mean
residence time of the reactor. Then 0.5 times the mean residence time for a further two
mean residence times of the reactor. The pH was measured for each of the samples.
Multiple Reactors:
Two reactors in series were connected. The Agitator has started, and the pump was
adjusted to give a volumetric flow rate of about 400 mL/min.
The exact volumetric flow rate was measured, and the mean residence time was
calculated using the equation 2. At time zero a 20 mL of 5 mol/L acetic acid was
injected into the first reactor at the base of the reactor under the impellor.
Samples of about 10 mL were collected at 0.1 times the mean residence time for one
mean residence time of the reactors. Then 0.5 times the mean residence time for a
further two mean residence times of the reactors. The pH was measured for each of the
samples.
36
3. Results and Discussion.
● Prework
1. Single CSTR plotting using equation below with an assumption t=1.
CA ( −t )
=e t
C Ao
1.2
0.8
CA/CA0
0.6
Theroretical (t*10)
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3 3.5
t (0.1min)
C Ao (n−1)!
0.8
0.7
0.6
0.5
CA/CA0
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5
Time (minutes)
37
3. Plotting of CSTR with single, two, five, and ten reactors in series.
CA n
n
n−1 −n t
RTD= = t e r
C Ao (n−1)! r
1.4
CSTR 5 in Series
1
CSTR 10 in Series
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Mean Residence Time
Figure 5 Plotting of Mean residence time vs CA/CA0 in various CSTR.
Experimental Results
The concentration of tracer (acetic acid) is measured in periodic times. The plotting of the ratio
of tracer concentration at specific time and initial tracer concentration versus mean residence
time is compared with the theoretical line for both single CSTR and Two Reactor in Series. The
results are depicted in Figure 6 and 7, respectively. The experimental results of single reactor
and two reactor is reasonably good flow compared to theoretical results based on the Figure 1
illustration for single reactor and Figure 8 for two reactors. Despite the good flow results, it still
has deviation to theoretical results.
Deviation between experimental line and theoretical line of single CSTR is probably caused by
imperfect mixing that can result in tracer concentration gradient within the reactor [4]. To
provide homogeneous reaction conditions and high conversion rates, single CSTR rely on
efficient reactant mixing. Imperfect mixing can result in concentration gradients, incomplete
reactions, and overall reactor misbehaviour. In some conditions, temperature and pressure
gradients in single reactor can alter reaction rate and selectivity. These gradients can be
caused by insufficient mixing, heat transmission limits, or other circumstances, resulting in the
reactor misbehaving [14]. Another reason for imperfect mixing is the speed of the stirrer
factor. The stirrer speed must be set to the optimum value to achieve the best mixing results.
Suppose the stirrer speed is not in optimum value. In that case, the contact between the acetic
acid and the water will not be in perfect condition, thus reducing the conversion rate and
resulting in deviation results [4] [15] .
38
Deviation in two reactor results is also caused by imperfect conditions and parameters,
resulting in imperfect mixing between acetic acid and water and misbehaving in two reactor
results. The imperfect condition is like the single reactor condition, which is a non-optimum
temperature, pressure, and stirring speed, which leads to deviation results between
experimental and theoretical results. Additionally, In the two-reactors experiment, the
presence of dead zone or stagnant fluid regions in the reactor and pipe is more significant than
in a single reactor [4]. This dead zone can affect the conversion rate performance and lead to
deviation.
1.2
1
Experiment
0.8
Theoritical
CA/CA0
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3 3.5
Mean Residence Time
0.8
0.7
0.6
0.5 Theoretical
Experimental
CA/CA0
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5
Time (min)
39
Figure 8. Misbehaving in plug flow reactor for two reactor experimental results [4].
4. Conclusion.
40
References
[1] J. E. E. O. G. O. B.C. Tanda, “The kinetics of chalcopyrite leaching in alkaline glycine/glycinate solutions,”
Minerals Engineering, vol. 135, pp. 118-128, 2019.
[3] M. E. W. Hong Yong Sohn, Rate processes of extractive metallurgy, New York: Plenum Press, 1979.
[4] O. Levenspiel, “Chemical Reaction Engineering.,” John Wiley & Sons., 1999.
[5] T. Havlik, Hydrometallurgy principles and application.., England: Cambridge International Science
Publishing Limited, 2008.
[6] J. W. X. W. Y. W. Y. &. G. F. Hao, “Optimizing the Leaching Parameters and Studying the Kinetics of Copper
Recovery from Waste Printed Circuit Boards,” ACS Omega, vol. 7, pp. 3689-3699, 2022.
[7] R. a. I. T. KERBY, “Effect of Impurities on the Current Efficiency of Zinc Electrodeposition,” Dept. Energy,
Mines and Resources, Mines Branch,, Ottawa, 1971.
[8] W. J. CHMIELEWSKI T., “ AMMONIA PRESSURE LEACHING FOR LUBIN,” Physicochemical Problems of
Mineral Processing, pp. 5-20, 2009.
[9] V. P. K. &. D. N. A. M. S. Prasad, “Development of SX-EW Process for Copper Recovery,” Mineral Processing
and Extractive Metallurgy Review, vol. 8, pp. 1-4,95-118, 1992.
[10] T. C. Katarzyna Rotuska, “GROWING ROLE OF SOLVENT EXTRACTION IN COPPER ORES PROCESSING,”
Physicochemical Problems of Mineral Processing, vol. 42, pp. 29-36, 2008.
[11] G. A. Kordosky, “Copper solvent extraction: The state of the art,” Journal of Metals, vol. 44, pp. 40-45,
1992.
[12] B. T. B. A. L. D. Hugo Guerra, “Nanofiltration for treatment of Western Australian bore water for mineral
processing operations: A pilot scale study,” Journal of Water Process Engineering, vol. 52, 2023.
[14] F. G. L. v. H. J. v. E. D. S. &. B. J. H. (. .. ,. 3. van der Klis, “From batch to continuous: Au-catalysed oxidation
of d-galacturonic acid in a packed bed plug flow reactor under alkaline conditions,” Reaction Chemistry &
Engineering, vol. 3, no. 4, pp. 540-549, 2018.
[15] R. H. L. &. P. C. A. (. Marinho, “Effect of Stirring Speed on Conversion and Time to Particle Stabilization of
Poly (Vinyl Chloride) Produced by Suspension Polymerization Process at the Beginning of Reaction,” .
Brazilian Journal of Chemical, vol. 35, no. 2, pp. 631-640, 2018.
41