05bf56e9-0eb6-4287-9def-30f74045b349 24 Alvaro Diaz Falconi 4 5 false false 0 0 decimal

&[DATE] &[TIME] - &[FILENAME]

pTM 2020-04-17 Thermo Calc Epi

Constants

Rpg 8.314 :

J/K/mol

cal2J 4.184 :

J

Process

Tg 600 273.15 + :

KAs 1.62 1022^* Rpg Tg * (5 -^* :

KSb 4.824 1021^* Rpg Tg * (5 -^* :

Thermo

Ωi 3843 cal2J * :

J/mol

KAs 803.5643

KSb 239.2836

VoIII 3 :

PzAl 1 :

PzSb PzAs VoIII 4 / PzAl * PzSb - :

VoIII 4 / PzAl * 0.75

Partial Pressures

Pas Psb PzSb Pal PzAl 4 PzSb PzAs Pas - (* - 4 PzSb Psb - (* - :

Pas Psb PzSb x PzSb Psb - PzSb PzAs Pas - PzSb + Psb - / :

Solid Mole Fraction

Pas Psb PzSb aAs 1 Pas Psb PzSb x - (e Ωi Rpg Tg * / Pas Psb PzSb x 2^*^* :

Activities

Pas Psb PzSb aSb Pas Psb PzSb x e Ωi Rpg Tg * / 1 Pas Psb PzSb x - (2^*^* :

Pas Psb PzSb KAlAs Pas Psb PzSb aAs Pas Psb PzSb Pal Pas 14 /^* / :

Equilibrium Constants

Pas Psb PzSb KAlSb Pas Psb PzSb aSb Pas Psb PzSb Pal Psb 14 /^* / :

Given

x target Pas Psb 12 mat x 1 el x 2 el 12 mat : Pas Psb PzSb KAlAs KAs - Pas Psb PzSb KAlSb KSb - 21 mat 21 line :

Jx x target x 12 range el vectorize Jacob :

x J Jx eval :

Solve

ε 1012 -^:

PzSb guess f sol guess 0 ε target J al_nleqsolve : sol target normi ε sqrt < sol 00 stack if 21 line :

Check function

guess 1 0.51 21 mat :

0.5 guess f 0021 mat

1.5 guess f 0021 mat

PzSb 0.1102 :

Pas Psb 21 mat PzSb guess f : 0.3965 0.1039 21 mat

Pas Psb 21 mat target 6.2711 1010 -^* 1.8797 1010 -^* 21 mat

Output

PzSb PzAs 0.6398

4 PzSb PzAs PzSb + PzAl / * 3

Review

PzSb PzSb PzAs PzSb + / 0.1469

Pas Psb PzSb x 0.0254

Pas Psb PzSb Pal 0.0015

Pas Psb PzSb aAs 0.976

Pas Psb PzSb aSb 0.2079

Program

N 100 :

ko 01 matrix eval :

k IP M6 k el M6 k 1 - el - M7 k el M7 k 1 - el - / M6 k 1 - el M6 k 2 - el - M7 k 1 - el M7 k 2 - el - / - sign :

k 1 N 1 - (range PzAs VoIII 4 / 1 k N / - (* : PzSb VoIII 4 / k N / * : Pas Psb 21 mat PzSb Pas Psb 21 mat f : M1 k el k : M2 k el PzAs : M3 k el PzSb : M4 k el Pas : M5 k el Psb : M6 k el Pas Psb PzSb Pal : M7 k el PzSb PzAs PzSb + / : M8 k el Pas Psb PzSb x : k 3 ≡ σ k IP : k 3 > (σ k IP ≠ (& ko ko k 2 - stack : σ σ - : 21 line if 121 line for

M7 M4 augment M7 M5 augment M7 ko el M4 ko el o augment M7 ko el M5 ko el o augment 4 1 sys M7 M2 augment M7 M3 augment 2 1 sys M7 M4 augment M7 M5 augment M7 ko el M4 ko el o augment M7 ko el M5 ko el o augment 4 1 sys M7 M2 augment M7 M3 augment 2 1 sys

Inflexion point of the Partial pressure

M7 M8 augment M7 ko el M8 ko el o augment 2 1 sys M7 M6 augment M7 ko el M6 ko el o augment 2 1 sys M7 M8 augment M7 ko el M8 ko el o augment 2 1 sys M7 M6 augment M7 ko el M6 ko el o augment 2 1 sys

Alvaro