RAYUDU 4a, B, C 5a, B Statement
Uploaded by
Sai VarmaRAYUDU 4a, B, C 5a, B Statement
Uploaded by
Sai VarmaEXPERIMENT NO: 4a D.
SUBBARAYUDU
DATE: 19021A0805
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
𝑑𝑉
= 𝐹0 − 𝐹
𝑑𝑡
𝑑(𝑉𝐶𝐴 )
= 𝐹0 𝐶𝐴0 − 𝐹𝐶𝐴 − 𝑉𝑘𝐶𝐴
𝑑𝑡
𝑑(𝑉𝑇) 𝜆𝑉𝑘𝐶𝐴 𝑈𝐴𝐻
= 𝐹0 𝑇0 − 𝐹𝑇 − − (𝑇 − 𝑇𝐽 )
𝑑𝑡 𝜌𝐶𝑃 𝜌𝐶𝑃
𝑑𝑇𝐽 𝐹𝐽(𝑇𝐽0 −𝑇𝐽) 𝑈𝐴𝐻
= + (𝑇 − 𝑇𝐽 )
𝑑𝑡 𝑉𝐽 𝜌𝐽 𝑉𝐽 𝐶𝐽
𝑑𝐹 𝑑𝑉
𝐹 = 40 − 10(48 − 𝑉); = 10 𝑑𝑡 𝑘
𝑑𝑡 𝐴→ 𝐵
𝑑𝐹𝐽 𝑑𝑇 Figure 4.1 Nonisothermal CSTR system
𝐹𝐽 = 49.9 − 𝐾𝑐 (600 − 𝑇); = 𝐾𝑐
𝑑𝑡 𝑑𝑡
Where, 𝑘 = 𝛼𝑒 (−𝐸⁄𝑅𝑇); 𝐶𝐴0 = 0.55 𝑙𝑏 𝑚𝑜𝑙𝐴⁄𝑓𝑡 3
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 D.SUBBARAYUDU
DATE: 19021A0805
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 D.SUBBARUYUDU
DATE: 19021A0805
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
𝑑𝑉
= 𝐹0 − 𝐹
𝑑𝑡
𝑑(𝑉𝐶𝐴 )
= 𝐹0 𝐶𝐴0 − 𝐹𝐶𝐴 − 𝑉𝑘𝐶𝐴
𝑑𝑡
𝑑(𝑉𝑇) 𝜆𝑉𝑘𝐶𝐴 𝑈𝐴𝐻
= 𝐹0 𝑇0 − 𝐹𝑇 − − (𝑇 − 𝑇𝐽 )
𝑑𝑡 𝜌𝐶𝑃 𝜌𝐶𝑃
𝑑𝑇𝐽 𝐹𝐽(𝑇𝐽0 −𝑇𝐽) 𝑈𝐴𝐻
= + (𝑇 − 𝑇𝐽 )
𝑑𝑡 𝑉𝐽 𝜌𝐽 𝑉𝐽 𝐶𝐽
𝑑𝐹 𝑑𝑉
𝐹 = 40 − 10(48 − 𝑉); = 10 𝑑𝑡 𝑘
𝑑𝑡 𝐴→ 𝐵
𝑑𝐹𝐽 𝑑𝑇 Figure 4.1 Nonisothermal CSTR system
𝐹𝐽 = 49.9 − 𝐾𝑐 (600 − 𝑇); = 𝐾𝑐
𝑑𝑡 𝑑𝑡
Where, 𝑘 = 𝛼𝑒 (−𝐸⁄𝑅𝑇); 𝐶𝐴0 = 0.55 𝑙𝑏 𝑚𝑜𝑙𝐴⁄𝑓𝑡 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 D. SUBBARAYUDU
DATE: 19021A0805
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
𝑑𝑉
= 𝐹0 − 𝐹
𝑑𝑡
𝑑(𝑉𝐶𝐴 )
= 𝐹0 𝐶𝐴0 − 𝐹𝐶𝐴 − 𝑉𝑘𝐶𝐴
𝑑𝑡
𝑑(𝑉𝑇) 𝜆𝑉𝑘𝐶𝐴 𝑈𝐴𝐻
= 𝐹0 𝑇0 − 𝐹𝑇 − − (𝑇 − 𝑇𝐽 )
𝑑𝑡 𝜌𝐶𝑃 𝜌𝐶𝑃
𝑑𝑇𝐽 𝐹𝐽(𝑇𝐽0 −𝑇𝐽) 𝑈𝐴𝐻
= + (𝑇 − 𝑇𝐽 )
𝑑𝑡 𝑉𝐽 𝜌𝐽 𝑉𝐽 𝐶𝐽
𝑑𝐹 𝑑𝑉
𝐹 = 40 − 10(48 − 𝑉); = 10 𝑑𝑡 𝑘
𝑑𝑡 𝐴→ 𝐵
𝑑𝐹𝐽 𝑑𝑇 Figure 4.1 Nonisothermal CSTR system
𝐹𝐽 = 49.9 − 𝐾𝑐 (600 − 𝑇); = 𝐾𝑐
𝑑𝑡 𝑑𝑡
Where, 𝑘 = 𝛼𝑒 (−𝐸⁄𝑅𝑇); 𝐶𝐴0 = 0.55 𝑙𝑏 𝑚𝑜𝑙𝐴⁄𝑓𝑡 3
Table 4.1 Nonisothermal CSTR parameter values
c) Solve the above set of differential equations in MATLAB programming using Runge Kutta 4th
order method.
