Page 245 - Separation process engineering

Cells(10, 4) = xBbot
Cells(10, 6) = xCbot
' initialize (top stage =1) and start loops
i = 1
yA = xAdist
yB = xBdist
yC = xCdist
' Calculations in enriching section: first equilibrium.
Do While i < feedstage

Eqsum = (yA / alphaAB) + (yB / alphaBB) + (yC / alphaCB)
xA = (yA / alphaAB) / Eqsum
xB = (yB / alphaBB) / Eqsum
xC = (yC / alphaCB) / Eqsum
' Record values on spreadsheet
Cells(i + 12, 1).Value = i
Cells(i + 12, 2).Value = xA
Cells(i + 12, 3).Value = yA
Cells(i + 12, 4).Value = xB
Cells(i + 12, 5).Value = yB
Cells(i + 12, 6).Value = xC
Cells(i + 12, 7).Value = yC
' Top operating line

i = i + 1
yA = LoverV * xA + (1 - LoverV) * xAdist
yB = LoverV * xB + (1 - LoverV) * xBdist
yC = LoverV * xC + (1 - LoverV) * xCdist
Loop
' Calculations in stripping section
Do
' Save values of x from previous stage.
xAold = xA
xBold = xB
xCold = xC
' Equilibrium calculation for x values for current stage number i.
Eqsum = (yA / alphaAB) + (yB / alphaBB) + (yC / alphaCB)
xA = (yA / alphaAB) / Eqsum
xB = (yB / alphaBB) / Eqsum
xC = (yC / alphaCB) / Eqsum

' Record values on spreadsheet.
Cells(i + 12, 1).Value = i
Cells(i + 12, 2).Value = xA
Cells(i + 12, 3).Value = yA
Cells(i + 12, 4).Value = xB
Cells(i + 12, 5).Value = yB
Cells(i + 12, 6).Value = xC
Cells(i + 12, 7).Value = yC
' Test for feed stage too high