CHE 317: COMPUTER PACKAGES FOR
CHEMICAL ENGINEERS
(EXCEL SECTION)
LECTURE IV
Instructors:
Engr. (Dr) Efeovbokhan, V.E and Mr. A. T. Ogunbiyi
Solver method : simultaneous linear
equations
Determine the solution to the following system of
equations using Solver.
Solver method : simultaneous linear
equations
Change equations to the form: f(x)=0
F1(x):
3X1 + 2X2- X3 4 = 0
F2(x):
2X1 - X2 + X3 3 = 0
F3(x):
X1 + X2 - 2X3 + 3 = 0
Finding the sum of the squares of each:
Y= F12 +
F22 +
F32
Minimizing the sum of the squares
If F1 = 0, F2 = 0 and F3 = 0, then Y will also equal zero.
We will use Excels Solver feature to determine the
values of X1, X2 and X3 that make the target function, Y
zero.
Solver method : simultaneous linear
equations
Initial guess for X1
Enter this formula
in B6. Remember
to begin with =
Y = F12 + F22 + F32
F1(X):
F2(X):
F2X):
3X1+2X2-X3 = 4
2X1-X2+X3 = 3
X1+X2-2X3 = -3
Solving Simultaneous Nonlinear Equations
in Excel Using Solver
Class Exercise
Determine the solution to the following two
simultaneous nonlinear equations using Solver.
Restrict the range of the independent variables to
nonnegative values.
Molar Volume and Compressibility Factor from Van der
Waals Equation
- - - - - - - -- 1
Rearrange eqn.1 in the form f(x)=0
Calculate the constants for the given gas and apply goalseek.
*It is very important to arrange your presentation properly
*NEAT PRESENTATION OF YOUR WORK IS VERY ESSENTIAL
Equation of State (EOS) Property Calculation
INPUT DATA
Chemical Name
METHANE
Critical Temperature
( Tc )
190.60
Deg K
Critical Pressure
( Pc )
46.00
Bar
Accentric Factor
()
0.0080
Molecular Weight
16.043
Peng Robinson
(PR)
Equation of State
( EOS )
Pressure
(P)
1.00
Temperature
(T)
100.00
RESULT
Bar
Deg K
Molar Volume and Compressibility Factor from Van der
Waals Equation
Problem 1
Value
P = Pressure in atm=Pr*Pc
V = molar volume in liters/g-mol
T= temperature in K
R= gas constant atm-liter/g-molK
Tc= citical temp for ammonia
Pc=critcal pressure for ammonia
56
0.574892698
450
0.08206
405.5
111.3
Preduced
Z=compressibilty factor=PV/RT
4.043126685
0.871827961
a
b
4.196945872
0.037371215
Value of RT
Value of P+a/V^2
Value of V-b
LHS
f(V)=0
36.927
68.69871519
0.537521483
36.92703527
3.52743E-05
Preduced
56
0.503 0.575084 0.871819
Molar Volume and Compressibility Factor
from Redlich Kwong Equation
Rearrange equation in the form f(x)=0
LHS RHS=0
SOLVER CAN NOW BE APPLIED
YOU MAY ALSO SOLVE for the RHS=P, since P is a constant
Molar Volume and Compressibility Factor from
REDLICH KWONG Equation
Value
P = Pressure in atm=Pr*Pc
V = molar volume in liters/g-mol
T= temperature in K
R= gas constant atm-liter/g-molK
Tc= citical temp for ammonia
Pc=critcal pressure for ammonia
56
0.68531
450
0.08206
405.5
111.3
Preduced
Z=compressibilty factor=PV/RT
4.04313
1.03928
a
b
0.00362
0.0259
Value of RT/(V-b)
Value of a/(V+b)VT.5
RHS
LHS
f(V)=0
56
0.00035
55.9996
56
3.9E-07