Sri 4a, B, C 5a, B Statement
Uploaded by
Sai VarmaSri 4a, B, C 5a, B Statement
Uploaded by
Sai VarmaEXPERIMENT NO: 4a G.A.M.
SRINIVAS
DATE: 19021A0824
4. NON-ISOTHERMAL CSTR SYSTEM
PROBLEM STATEMENT:
Consider the following set of nonlinear differential equations that describe non-isothermal CSTR
system as shown in Figure 4.1 and Table 4.1 gives the parameter values and steady state
conditions
dV
=F 0−F
dt
d (V C A)
=F 0 C A 0−F C A−Vk C A
dt
d (VT ) λVk C A U A H
=F0 T 0−FT − − ( T −T J )
dt ρ CP ρC P
d T J F J (T −T ) U A H
dt
=
VJ
J0
+ J
ρJ V J C J
( T −T J )
AkB
dF dV
F=40−10 ( 48−V ) ; =10 →
dt dt Figure 4.1 Nonisothermal CSTR system
d FJ dT
F J =49.9−K c ( 600−T ) ; =K c
dt dt
3
Where, k =α e(−E / RT ); C A 0=0.55lb molA /ft
Table 4.1 Nonisothermal CSTR parameter values
a) Solve the above set of differential equations in MATLAB programming using ODE45 solver.
EXPERIMENT NO: 5 a,b G.AM. SRINIVAS
DATE: 19021A0824
5. PLUG FLOW REACTOR SYSTEM
Problem Statement:
Plug flow reactor performance:
a. A homogeneous gas reaction A→3R has a reported rate at 215oC
-rA=10-2CA1/2 mol/liter.sec
Find the space time needed for 80% conversion of 50% A- 50% inert feed to a plug flow
reactor shown in Figure 5.1., operating at 215oC and 5 atm, CAo=0.0625 mol/liter.
Figure 5.1: Plug flow reactor
Problem Statement:
Plug flow reactor volume:
b. The homogeneous gas decomposition of phosphine as shown in Figure 5.2,
4PH3(g) P4(g) +6H2
proceeds at 649oC with the first-order rate
-rPH3= (10/hr)CPH3
What size of plug flow reactor operating at 649°C and 460 kPa, can produce 80%conversion
of a feed consisting of 40 mol of pure phosphine per hour? Solve this using Simpsons rule
and Trapezoidal rule in MATLAB Programming.
Figure 5.2: Plug Flow Reactor
EXPERIMENT NO: 4b G.A.M.SRINIVAS
DATE: 19021A0824
4. NON-ISOTHERMAL CSTR SYSTEM
PROBLEM STATEMENT:
Consider the following set of nonlinear differential equations that describe non-isothermal CSTR
system as shown in Figure 4.1 and Table 4.1 gives the parameter values and steady state
conditions
dV
=F 0−F
dt
d (V C A)
=F 0 C A 0−F C A−Vk C A
dt
d (VT ) λVk C A U A H
=F0 T 0−FT − − ( T −T J )
dt ρ CP ρC P
d T J F J (T −T ) U A H
dt
=
VJ
J0
+ J
ρJ V J C J
( T −T J )
AkB
dF dV
F=40−10 ( 48−V ) ; =10 →
dt dt Figure 4.1 Nonisothermal CSTR system
d FJ dT
F J =49.9−K c ( 600−T ) ; =K c
dt dt
Where, k =α e(−E / RT ); C A 0=0.55lb molA /ft 3
Table 4.1 Nonisothermal CSTR parameter values
b) Solve the above set of differential equations in MATLAB programming using Euler’s
method.
EXPERIMENT NO: 4c G.A.M.SRINIVAS
DATE: 19021A0824
4. NON-ISOTHERMAL CSTR SYSTEM
PROBLEM STATEMENT:
Consider the following set of nonlinear differential equations that describe non-isothermal CSTR
system as shown in Figure 4.1 and Table 4.1 gives the parameter values and steady state
conditions
dV
=F 0−F
dt
d (V C A)
=F 0 C A 0−F C A−Vk C A
dt
d (VT ) λVk C A U A H
=F0 T 0−FT − − ( T −T J )
dt ρ CP ρC P
d T J F J (T −T ) U A H
dt
=
VJ
J0
+ J
ρJ V J C J
( T −T J )
AkB
dF dV
F=40−10 ( 48−V ) ; =10 →
dt dt Figure 4.1 Nonisothermal CSTR system
d FJ dT
F J =49.9−K c ( 600−T ) ; =K c
dt dt
3
Where, k =α e(−E / RT ); C A 0=0.55lb molA /ft
Table 4.1 Nonisothermal CSTR parameter values
c) Solve the above set of differential equations in MATLAB programming using Runge Kutta
4th order method.
