From Previous Problem

Download as xlsx, pdf, or txt
Download as xlsx, pdf, or txt
You are on page 1of 12

UNIVERSITY: NNAMDI AZIKIWE UNIVERSITY, AWKA

FACULTY: ENGINEERING
DEPARTMENT: CHEMICAL ENGINEERING
COURSE TITLE: COMPUTER APPLICATION IN CHEMICAL ENGINEERING
COURSE CODE: CHE463
SEMESTER: 1ST SEMESTER 2020/20201 SESSION
LECTURER: ENGR DR UZOH CHIGOZIE FRANCOLLINS
NAME: NZEKWE SAMUEL CHIGOZIE
REG NUMBER: 2017214013

QUESTIONS 2.32:
Find the compressibility factor of ammonia gas at conditions from 50 to 250atm and 400k
using the Redlich-kwong Soave equation of state in Excel. ( Hint: before begining your
spreadsheet, think about how you can organize it so that you can copy formulas from cell to
cell easily.)
Experimental data is available that gives a volume of 0.459l/g.mol at 66.8atm and 446K

From Previous Problem


Parameters - Ammonia Data
Tc 405.5 K T 400
Pc 111.3 atm P 200
R 0.08206 Liter-atm/g mol K
a 4.2527 gmol-1 Tr 0.986436
-1
b 0.0259 gmol α 1.006852
Molar Volume 0.56211270548739 L/gmol Pr 1.796945
Formula Result
A 5/2
0.42747(Pr/Tr ) 0.794818
B 0.08664(Pr/Tr) 0.157828
r AB 0.125445
q 2
B +B-A -0.61208
z -3.672324284E-07 0.317174
ERING

m 50 to 250atm and 400k


before begining your
copy formulas from cell to

l at 66.8atm and 446K

K
atm
UNIVERSITY: NNAMDI AZIKIWE UNIVERSITY, AWKA
FACULTY: ENGINEERING
DEPARTMENT: CHEMICAL ENGINEERING
COURSE TITLE: COMPUTER APPLICATION IN CHEMICAL ENGINEERING
COURSE CODE: CHE463
SEMESTER: 1ST SEMESTER 2020/20201 SESSION
LECTURER: ENGR DR UZOH CHIGOZIE FRANCOLLINS
NAME: NZEKWE SAMUEL CHIGOZIE
REG NUMBER: 2017214013

QUESTIONS 2.42:
Consider the following mixture going into a water-gas shift reactor to make hydrogen for the hydrog
CO, 630; H₂O, 1130; CO2, 189; H₂, 63 kmol/h. The gas is at 1 atm and 500 K. Use Excel (or MATLAB) t
the specific volume using(a) ideal gas law(b) Redlich-Kwong equation of state(c) Redlich-Kwong-Soa
of state.The acentric factors for the RK-Soave method are: CO, 0.049; water, 0.344; CO₂2.0.225; Hyd
-0.22.Where did you get the other data you needed? How do the three answers compare? Is the gas
not? Comment. Then redo the calculations for a pressure of 200 atm and comment on the results.

GIVEN DATA
T(K) P(atm) R(liter-atm/g.mol.k) Total Moles(kmol/hr) Total Moles(mol/s)
500 1 0.08206 2012 558.88
20

FINDING VALUES FOR (CO)

PARAMETERS FORMULA

Tc (k) Critical Temperature


Pc (bar) Critical Pressure
Tr B29/J36
a (RK) 0.42748*((D29^2 * J36^2)/J37)*J42
a(SRK) 0.42748*J43*((D29^2 * J36^2)/J37)
b 0.08664*((D29*J36)/J37)
α (RK) 1 / J38^0.5
α (SRK) (1 + J45*(1 - J38^0.5))^2
ω Pitzer acentric factor
m 0.48508 + 1.55171*J44 - 0.1561*J44^2
Number of moles(n) mass / molecularmass
Mole fraction(Y) J46 / E29
FINDING VALUES FOR (H2O)

PARAMETERS FORMULA

Tc (k) Critical Temperature


Pc (bar) Critical Pressure
Tr B29/J54
a (RK) 0.42748*((D29^2 * J54^2)/J55)*J60
a(SRK) 0.42748*((D29^2 * J54^2)/J55)*J61
b 0.08664*((D29*J54)/J55)
α (RK) 1 / J56^0.5
α (SRK) (1 + J63*(1 - J56^0.5))^2
ω Pitzer acentric factor
m 0.48508 + 1.55171*J62 - 0.1561*J62^2
Number of moles(n) mass / molecularmass
Mole fraction(Y) J64 /E29

FINDING VALUES FOR (CO2)

PARAMETERS FORMULA

Tc (k) Critical Temperature


Pc (bar) Critical Pressure
Tr B29/J72
a (RK) 0.42748*((D29^2 * J72^2)/J73)*J78
a(SRK) 0.42748*((D29^2 * J72^2)/J73)*J79
b 0.08664*((D29*J72)/J73)
α (RK) 1 / J74^0.5
α (SRK) (1 + J81*(1 - J74^0.5))^2
ω Pitzer acentric factor
m 0.48508 + 1.55171*J80 - 0.1561*J80^2
Number of moles(n) mass / molecularmass
Mole fraction(Y) J82 / E29

FINDING VALUES FOR (H2)

PARAMETERS FORMULA

Tc (k) Critical Temperature


Pc (bar) Critical Pressure
Tr B29/J90
a (RK) 0.42748*((D29^2 * J90^2)/J91)*J96
a(SRK) 0.42748*((D29^2 * J90^2)/J91)*J97
b 0.08664*((D29*J90)/J91)
α (RK) 1/J92^0.5
α (SRK) (1 + J99*(1 - J92^0.5))^2
ω Pitzer acentric factor
m 0.48508 + 1.55171*J98 - 0.1561*J98^2
Number of moles(n) mass / molecularmass
Mole fraction(Y) J100 / E29

SUMMATION OF a and b VALUES

SUMMATION FORMULA
∑a (RK) (J47*J39^0.5 + J65*J57^0.5 + J83*J75^0.5 + J101*J93^0.5)^2
∑a (SRK) (J47*J40^0.5 + J65*J58^0.5 + J83*J76^0.5 + J101*J94^0.5)^2

∑b J47*J41 + J65*J59 + J83*J77 + J101*J95

A) Calculation of specific volume using the ideal


gas law

PARAMETER FORMULA

AT 1 atm V D29*B29/C29
AT 200 atm v D29*B29/C30

B) Calculation of specific volume using REDLICH KWONG EQUATION OF STATE


(RK)
PARAMETER
FORMULA

F(v)

AT 1 atm RESULTS 2.90285396442602E-07


AT 200 atm 2.92489830117093E-06

C) Calculation of specific volume using REDLICH


KWONG SOAVE EQUATION OF STATE (SRK)
C) Calculation of specific volume using REDLICH
KWONG SOAVE EQUATION OF STATE (SRK)

PARAMETER FORMULA

F(v)

AT 1 atm RESULTS 1.77377983902327E-05


AT 200 atm 1.88832752115453E-06

DISCUSSION
I obtained the other data needed for the calculation from the ChemSep software, versi
The three answers obtained appear to be fairly close together, as they can be easily
approximated to 41 l/g mol. Therefore, it can be concluded that the gas behaves ideally
on which value tends to approach the ideal gas value, the Redlich-Kwong-Soave EOS ap
to be more accurate than the Redlich-Kwong EOS.
Unlike at 1 atm, the specific volume values here are not quite close together. Thus, as th
pressure is increased, the gas begins to deviate from ideal conditions. Consequently, th
Redlich-Kwong EOS appears to show greater deviation from the ideal gas value than the
Redlich-Kwong-Soave EOS.   Also, the specific volume decreases as the pressure is incre
from 1 atm to 200 atm.              
ydrogen for the hydrogen economy.
se Excel (or MATLAB) to compute
(c) Redlich-Kwong-Soave equation
.344; CO₂2.0.225; Hydrogen,
ers compare? Is the gas ideal or
ment on the results.

MOLE CONVERSIONS
1mol/s = 3.6Kmol/hr

RESULTS

132.7
34.532
3.7678975131877
0.7562221585157
0.3294288803167
0.0273211607691
0.5151698748957
0.224420605967
0.048
0.5592024256
630
0.3131212723658
RESULTS

647
220.5
0.7727975270479
6.2164999097474
6.8647193939542
0.0208615053279
1.1375412080448
1.2561572115593
0.343
0.9989515211
1130
0.5616302186879

RESULTS

304.21
73.8
1.6436014595181
2.8156014307636
2.1256037580091
0.0293067109223
0.7800128204075
0.5888611094731
0.224
0.8248305664
189
0.0939363817097

RESULTS

33.19
13
15.064778547755
0.0628446920603
0.0832982633145
0.0181515558535
0.2576431640855
0.3414962731073
-0.215
0.1442466275
63
0.0313121272366

RESULTS

3.378514495038

3.2299831083625

0.0235926186323

LA RESULTS

C29 41.03
C30 2.0515

40.97123
1.992287
RESULTS

40.97485
1.996075

emSep software, version 8.


they can be easily
he gas behaves ideally. Based
h-Kwong-Soave EOS appears

se together. Thus, as the


ons. Consequently, the
deal gas value than the
s the pressure is increased
UNIVERSITY: NNAMDI AZIKIWE UNIVERSITY, AWKA
FACULTY: ENGINEERING
DEPARTMENT: CHEMICAL ENGINEERING
COURSE TITLE: COMPUTER APPLICATION IN CHEMICAL ENGINEERING
COURSE CODE: CHE463
SEMESTER: 1ST SEMESTER 2020/20201 SESSION
LECTURER: ENGR DR UZOH CHIGOZIE FRANCOLLINS
NAME: NZEKWE SAMUEL CHIGOZIE
REG NUMBER: 2017214013

QUESTIONS 2.51:
Consider a mixture of 25% ammonia, the rest nitrogen and hydrogen in a 1:3 ratio. The gas
is at 270 atm and 550 K. Use Excel (or MATLAB) to compute the specific volume using(a)
ideal gas law(b) Redlich-Kwong equation of state(c) Redlich-Kwong-Soave equation of
state.Where did you get the data you needed? How do the three answers compare? Is the
gas ideal or not? Comment on the reasons the RKS equation might be expected to be better.
3 ratio. The gas
me using(a)
quation of
ompare? Is the
ted to be better.

You might also like