Mathematical Model

Consider the continuous stirred tank reactor (CSTR) shown below, where a simple, liquid-phase,
irreversible, exothermic reaction π΄π΄
ππ
β π΅π΅ takes place. A coolant is flowing in the jacket of the vessel in
order to regulate the temperature of the fluid in the reactor, ππ(π‘π‘), and control in this way the
concentration of component A, πΆπΆπ΄π΄(π‘π‘), in the outlet stream at the desired value. There are two
potential disturbances in the process, the temperature, ππππ
(π‘π‘), and the concentration of component
A πΆπΆπ΄π΄ππ(π‘π‘), in the inlet stream, while the manipulated variable is the inlet temperature of the cooling
stream in the jacket, ππππππ(π‘π‘). The cooling jacket should be considered to be at a uniform temperature
ππππ
(π‘π‘). The sensor used to measure the concentration of component A, πΆπΆπ΄π΄(π‘π‘), is an on-line infrared
spectrometer providing real-time quantification of the concentration.
For the purposes of the assignment, the simplification assumption can be made that, in the expected
operating temperature range, the reaction rate is a lot more affected by changes in the temperature,
ππ(π‘π‘), in the tank in comparison to the concentration of A, πΆπΆπ΄π΄(π‘π‘), so that the reaction rate can be
expressed with a pseudo-zero order expression where:
πππ΄π΄(π‘π‘) = πποΏ½ππ(π‘π‘)οΏ½πΆπΆπ΄π΄(π‘π‘) β πποΏ½ππ(π‘π‘)οΏ½πΆπΆπ΄π΄π΄π΄
Numerical values (assumed constant):
ο§ Overall heat transfer coefficient between the coolant and the fluid in the tank: ππ = 425 W/m2
-K
ο§ Heat transfer area between the coolant and the fluid in the tank: π΄π΄ = 3.5 m2
ο§ Inlet/Outlet flow rate of the fluid in the tank: πΉπΉ = 6.30 Γ 10β4 m3
/s
ο§ The volume of fluid in the tank: ππ = 0.375 m3
ο§ The density of the fluid in the inlet feed stream and in the vessel: ππ = 880 kg/m3
ο§ Specific heat of the fluid in the inlet feed stream and in the vessel: ππππ = 3680 J/kg-K
ο§ Reaction frequency factor: ππππ = 2.87 Γ 105
s-1
ο§ Reaction activation energy: πΈπΈπ΄π΄ = 6.47 Γ 104
J/mol
ο§ The heat of reaction: π₯π₯π¨π¨ππ = β2.79 Γ 104
J/mol
ο§ Density of the coolant: ππππ = 1650 kg/m3
ο§ Specific heat of the coolant: ππππππ = 4190 J/kg-K
ο§ The volume of coolant in the jacket: ππππ = 0.045 m3
F, Cai(t), Ti(t)
F, CA(t), T(t)
M
AC
k
A B ο£§ο£§β
Tci(t)
Tc(t)
CAm(t)
CA(t)
T(t)
ο§ Inlet/Outlet flow rate of coolant in the jacket: πΉπΉππ = 4.15 Γ 10β4 m3
/s
ο§ Steady-state outlet concentration of A: πΆπΆπ΄π΄π΄π΄ = 3300 mol/m3
ο§ Steady-state outlet temperature: πππ΄π΄ = 105 oC
ο§ The transfer function for the final control element: πΊπΊππ
(π π ) = πΎπΎππ, where πΎπΎππ = β10 C/atm
ο§ The transfer function for the online analyzer: πΊπΊππ(π π ) = πΎπΎππ, where πΎπΎππ = 1 (mol/m3
)/(mol/m3
)
a) Starting from the mathematical model of the system, determine all transfer functions that
describe this process in a deviation state in the open-loop. Determine initially any first-order
dependencies between variables and calculate their respective gains and time constants.
Furthermore, combine these first-order transfer functions so that the transfer functions
describing the effect of πΆπΆπ΄π΄ππ
β²
(π‘π‘), ππππ
β²
(π‘π‘) and ππππππ
β²
(π‘π‘) on πΆπΆπ΄π΄
β²
(π‘π‘), ππ
β²
(π‘π‘) and ππππ
β²
(π‘π‘) are derived. Clearly
show all steps taken and clearly indicate the units of all involved variables.
[Marks: 10/22]
b) Based on the obtained transfer functions in a) and using a generic PI controller transfer function
block construct in Simulink the block diagram of the closed-loop process. There are two
alternative but equivalent methods that the open-loop process can be drawn in this block
diagram. Include in the assignment report a screenshot of the Simulink diagram used.
[Marks: 2/22]
c) Consider that initially, the system is at steady-state and the following cases.
ο If your Student ID number is even: At time π‘π‘ = 0 there is a step-change in the inlet
the concentration of component A, πΆπΆπ΄π΄ππ(π‘π‘), equal to 200 mol/m3
.
ο If your Student ID number is odd: At time π‘π‘ = 0 there is a step-change in the temperature, ππππ
(π‘π‘),
equal to -10 oC.
i. Using a P controller, simulate in Simulink the process response, under the above
described step-change, for the following πΎπΎππ values: 0.001, 0.01, 0.1 and 1 atm/(mol/m3
).
Consider an adequate simulation time to reach a new steady state. Plot the responses of
πΆπΆπ΄π΄
β²
(π‘π‘), ππ
β²
(π‘π‘) and ππππ
β²
(π‘π‘) with time (one plot per variable including the responses for all πΎπΎππ
values). Discuss the following: Presence or not and trend, if any, of an offset, presence or
not, and trend, if any, of oscillatory behavior, stability or not of the response.
[Marks: 3/22]
ii. Starting from the general equation describing the closed-loop response of the feedback-controlledΒ system and using the Final Value Theorem, calculate the offset for the various
Kc values considered in part i. Compare the calculated values with those predicted by
Simulink and discuss trends in relation to theory.
[Marks: 2/22]
d) Consider again that the system is at a steady state. There is no change in the disturbances this
time. At time π‘π‘ = 0 there is a step change of 200 mol/m3 in the set-point of the process.
i. Using a PI controller, simulate in Simulink the process response, under the above-described step-change, for a constant πΎπΎππ value of 0.015 atm/(mol/m3
) and for the
following πππΌπΌ values: 1000 and 170 s. Consider an adequate simulation time to reach a
new steady state. Plot the responses of πΆπΆπ΄π΄
β²
(π‘π‘), ππ
β²
(π‘π‘) and ππππ
β²
(π‘π‘) with time (one plot per
variable including the responses for both πππΌπΌ values). Discuss the following: Presence or
not and trend, if any, of an offset, presence or not and trend, if any, of an oscillatory
behavior, stability, or not of the response.
[Marks: 3/22]
ii. Starting from the characteristic equation of this closed-loop system, investigate its
stability for the πππΌπΌ values considered in part i. Compare the findings with those predicted
by Simulink and discuss trends in relation to theory.
[Marks: 2/22]
Submission
1. Submit the report, including plots, calculations, and discussion of the results, as a single word or
pdf file via the Turnitin link available in Aberdeen.
(a) All calculations, equations, etc. should be typed within the file using an appropriate
equation editor. No scanned pictures of hand-written files should be provided.
(b) The filename of the report file should be EX40HC_CA2_XXX.pdf/.doc (XXX to be replaced
2. Not complying with the above specifications can incur penalties or invalidate your submission.
3. Deadline for the submission of all files is 12:00 noon Friday, December 18th, 2020.
4. Penalties for late or non-submission are as follows:
(a) up to one week late, 2 CGS points will be deducted
(b) up to two weeks late, 3 CGS points will be deducted
(c) up to three weeks late, 4 CGS points will be deducted, with a maximum achievable grade
of CGS D3.
(d) up to four weeks late, 5 CGS points will be deducted, with a maximum achievable grade
of CGS D3.
(e) more than four weeks late, NP will be recorded.
5. A completed plagiarism cover sheet should be included in your electronic submission. Advice
about avoiding plagiarism, the Universityβs Definition of Plagiarism, a Checklist for Students, and
instructions for TurnitinUK can be found in the following area of the Student Learning Service
website www.cabin.ac.Uk/sys/online-resources/avoiding-plagiarism. For more information on Mathematical Model check on this: https://en.wikipedia.org/wiki/Mathematical_model

Get 20% Discount on This Paper
Pages (550 words)
Approximate price: -

Try it now!

Get 20% Discount on This Paper

We'll send you the first draft for approval by at
Total price:
\$0.00

How it works?

Fill in the order form and provide all details of your assignment.

Proceed with the payment

Choose the payment system that suits you most.

Our Services

Australia Assessments has gained an international reputation of being the leading website in custom assignment writing services. Once you give us the instructions of your paper through the order form, we will complete the rest.

Essay Writing Services

As we work towards providing the best custom assignment services, our company provides assignment services for any type of academic essay. We will help you develop professionally written essays that are rich in content and free from plagiarism.

Our skilled team of professional writers will ensure that we help you craft a remarkable admission essay for your desired Master's program in your institution of choice. We won't stop there. Once you enter the job market, we will be available to secure you a position at your desired worksite by creating an outstanding portfolio or resume.

Our editorial team is always available for all editing and proofreading services. They check completed papers by our writers and also provide professional opinions to papers completed by our clients.

Technical papers

We harbor professional academic writers with different qualifications in diverse academic fields. As such, we are capable of handling both simple and technical papers. Ensure that you provide us with correct and complete instructions in the order form.