Institiute of Space Technology Matlab Oht Pourbaix Diagram of Aluminum
Institiute of Space Technology Matlab Oht Pourbaix Diagram of Aluminum
MatLab OHT
Pourbaix Diagram of Aluminum
Name:
Muhammad Tayyab Ahsan
Department:
Materials science and Engineering
Submitted To:
Sir Arhum
Date of submission:
29 Dec, 2015
Pourbaix diagram:
In chemistry, a Pourbaix diagram, also known as a potential/pH diagram, pE/pH diagram, maps
out possible stable (equilibrium) phases of an aqueous electrochemical system. Predominant ion
boundaries are represented by lines. As such a Pourbaix diagram can be read much like a
standard phase diagram with a different set of axes. Similarly to phase diagrams, they do not allow
for reaction rate or kinetic effects.
Pourbaix diagrams are also known as EH-pH diagrams due to the labeling of the two axes. The
vertical axis is labeled EH for the voltage potential with respect to the standard hydrogen
electrode (SHE) as calculated by the Nernst equation. The "H" stands for hydrogen, although other
standards may be used, and they are for room temperature only.
The horizontal axis is labeled pH for the -log function of the H+ ion activity.
The lines in the Pourbaix diagram show the equilibrium conditions, that is, where the activities are
equal, for the species on each side of that line. ourbaix diagrams are also known as E H-pH diagrams
due to the labeling of the two axes. The vertical axis is labeled EH for thevoltage potential with
respect to the standard hydrogen electrode (SHE) as calculated by the Nernst equation. The "H"
stands for hydrogen, although other standards may be used, and they are for room temperature
only.
The horizontal axis is labeled pH for the -log function of the H+ ion activity.
The lines in the Pourbaix diagram show the equilibrium conditions, that is, where the
activities are equal, for the species on each side of that line. On either side of the line, one
form of the species will instead be said to be predominant.
Uses
Pourbaix diagrams have several uses, for example in corrosion studies, geoscience and in
environmental studies
Code :
%Muhammad Tayyab Ahsan
%120301014
%Pourbaix diagram for Al
E=-3:2;
%defining range of volts
ph=5.07;
%eq(1): pH at which
[Al+3]=[AlO2-1]=1
ph(1:6)=ph;
%equalizing array of ph and
volts to plot graph
plot(ph,E,'r--'),ylim([-2.5 1.5])
%ploting a line at which
concentrations of both are equal
hold on
%to plot more lines on same
graph
Al=1;
%concentration of [Al+3]=1
E3=(-1.66)+((1.971*(10^-2)).*log(Al));
%eq(3): equation for equilibrium
between Al and Al+3
PH=-(E3+1.55)/0.0591;
%Solving equation to fing
maximum point of eq:3 and minimum point of eq:2
PH3=-2:.01:PH;
%defining range for line 3 in
graph
EE3(1:387)=E3;
%equalizing array of ph and
volts according to PH3 values come from calculation
plot(PH3,EE3)
%plotting of line 3 for values
of volts and ph
AlO2=1;
%concentration of [ALO2]=1
E4=E3:0.5:2;
%defining range of volts for
line 4
PH4(1:8)=PH;
%equalizing array of ph and
volts according to the values of volts
plot(PH4,E4)
%plotting of line 4
ph6=14.59+(log(AlO2));
%equation for line 6
pH=PH:0.01:ph6;
%defining range of ph for line 2
E2=-1.55-(0.0591.*pH);
%eq(2): equation for equilibrium
between to solids Al and Al2O3
plot(pH,E2)
%plot for line 2 in the graph
PH6=ph6:16;
%defining range of ph for line 2
E5=-1.26-(0.0789.*PH6)+(0.0198*log(AlO2));
%relation between voltage and ph
i.e. the equilibrium is dependent on both ph and E
plot(PH6,E5)
%poltting line 5 on graph
E6=E5:0.5:2;
%defining range of volts for
line 6
ph6(1:9)=ph6;
%equalizing array of ph for line
6 on graph according to the values of volts calculated
plot(ph6,E6)
%plotting of line 6 on graph
pH_gen=-2:16;
%general range of ph in the
pourbaix diagram
Ea=0.12-(0.0591.*pH_gen);
%eq (2): with initial value of
volts Eo=0.12 from literature
plot(pH_gen,Ea,'g--')
%plotting reference lines
between general value of ph and initial values of volts
Eb=1.30-(0.0591.*pH_gen);
valueof volt Eo=0.30 from literature
plot(pH_gen,Eb,'g--')
xlabel('pH')
ylabel('E (volts) SHE'),
title('POURBAIX DIAGRAM OF ALUMINUM'),
hold off
further proceedings.
Output: