Maple Solutions to the Chemical Engineering Problem Set

Reaction Equilibrium for Multiple Gas Phase Reactions
We begin by creating a list of component identities and to calculate the number of components.
> components:=[A,B,C,D,X,Y,Z]; nc:=nops(components);
Next, we create a list of reactions:
> Reactions:=[R1,R2,R3];
components := [ A, B, C, D, X, Y, Z]
nc := 7
Reactions := [ R1, R2, R3 ]
Each reaction may be expressed as a Maple equation:
> Reaction[R1] :=A+B=C+D;
> Reaction[R2] := B+C=X+Y;
> Reaction[R3] :=A+X=Z;
Reaction := A + B = C + D
R1
Reaction := B + C = X + Y
R2
Reaction := A + X = Z
R3
The stoichiometric cofficients can be deduced from the set of reactions as follows:
> for r in Reactions do
for i in components do
nu[i,r]:=coeff(lhs(Reaction[r])-rhs(Reaction[r]),i);
od: od:
print(nu);
table([
( Z, R3 ) = -1
( Y, R2 ) = -1
( X, R1) = 0
( A, R2) = 0
( B, R3) = 0
( X, R3) = 1
( D, R2) = 0
( C, R1 ) = -1
( Z, R2 ) = 0
( Y, R1 ) = 0
( B, R2) = 1
( C, R3 ) = 0
( Y, R3 ) = 0
( X, R2) = -1
( D, R1) = -1
( A, R1) = 1
( A, R3) = 1
Page 20