| 1 | package Reactor
|
|---|
| 2 | // Package for study of nonlinear reactor adapted from Seborg et al.
|
|---|
| 3 | // author: Bernt Lie
|
|---|
| 4 | // University College of Southeast Norway
|
|---|
| 5 | // June 16, 2016
|
|---|
| 6 | //
|
|---|
| 7 | model ModReactor
|
|---|
| 8 | // Reactor model in Ex 2.5 of Seborg et al. (2011),
|
|---|
| 9 | // Process Dynamics and Control
|
|---|
| 10 | //
|
|---|
| 11 | // author: Bernt Lie
|
|---|
| 12 | // University College of Southeast Norway
|
|---|
| 13 | // June 16, 2016
|
|---|
| 14 | //
|
|---|
| 15 | // Parameters
|
|---|
| 16 | parameter Real V = 100 "Reactor vessel volume; L";
|
|---|
| 17 | parameter Real k0 = exp(8750/350) "Preexponential factor, reaction; 1/min";
|
|---|
| 18 | parameter Real EdR = 8750 "Activation 'temperature'; K";
|
|---|
| 19 | parameter Real dHr = -5e4 "Specific reaction enthalpy; J/mol";
|
|---|
| 20 | parameter Real rho = 1e3 "Density, reactor fluid; g/L";
|
|---|
| 21 | parameter Real cp = 0.239 "Specific heat capacity, reactor fluid; J/(g.K)";
|
|---|
| 22 | parameter Real UA = 5e4 "Heat transfer parameter; J/(min.K)";
|
|---|
| 23 | // Initial state parameter
|
|---|
| 24 | parameter Real cA0 = 0.5 "Initial cA; mol/L";
|
|---|
| 25 | parameter Real nA0 = cA0*V "Initial 'real' state nA; mol";
|
|---|
| 26 | parameter Real T0 = 350 "Initial T; K";
|
|---|
| 27 | // Declaring variables
|
|---|
| 28 | // -- states
|
|---|
| 29 | Real nA(start = nA0, fixed = true); // initial # moles of A in vessel; mol
|
|---|
| 30 | Real T(start = T0, fixed = true); // initial temperature in vessel; K
|
|---|
| 31 | // -- auxiliary variables
|
|---|
| 32 | Real ndAi "Influent of A; mol/min";
|
|---|
| 33 | Real ndAe "Effluent of A; mol/min";
|
|---|
| 34 | Real ndAg "Generation rate of A; mol/min";
|
|---|
| 35 | Real r "Rate of reaction; mol/(L.min)";
|
|---|
| 36 | Real Qd "Heat flow from cooling fluid to reactor vessel; J/min";
|
|---|
| 37 | // Inputs
|
|---|
| 38 | input Real Vd "Volumetric flow through reactor vessel; L/min";
|
|---|
| 39 | input Real cAi "Reactor vessel influent concentration of A; mol/min";
|
|---|
| 40 | input Real Ti "Reactor vessel influent temperature; K";
|
|---|
| 41 | input Real Tci "Cooling circuit influent temperature; K";
|
|---|
| 42 | // Outputs
|
|---|
| 43 | output Real cA "Reactor concentration; mol/L";
|
|---|
| 44 | equation
|
|---|
| 45 | nA = V*cA;
|
|---|
| 46 | ndAi = Vd*cAi;
|
|---|
| 47 | ndAe = Vd*cA;
|
|---|
| 48 | ndAg = -V*r;
|
|---|
| 49 | r = k0*exp(-EdR/T)*cA;
|
|---|
| 50 | Qd = UA*(Tci-T);
|
|---|
| 51 | //
|
|---|
| 52 | der(nA) = ndAi - ndAe + ndAg;
|
|---|
| 53 | V*rho*cp*der(T) = rho*cp*Vd*(Ti-T) + (-dHr)*r*V + Qd;
|
|---|
| 54 | end ModReactor;
|
|---|
| 55 | // End package
|
|---|
| 56 | end Reactor;
|
|---|
