Assignment 1 - 2014

Due date(s):	20 January 2014. in class
(PDF)	Assignment questions

<rst> <rst-options: 'toc' = False/> <rst-options: 'reset-figures' = False/>

.. rubric:: Assignment objectives: understanding about feedback systems; dynamic models

.. question:: :grading: 4

Select a system from everyday life that employs feedback in its operation. Describe the system and how feedback is used. A sketch of the system is required.

.. question:: :grading: 8

Select a process that you have studied in a previous course in chemical engineering (e.g. heat exchangers, separation process, reactors) and describe the objectives that should be associated with each of the seven categories of control objectives described in Chapter 2 of the textbook. A sketch of the process is required.

.. question:: :grading: 4

What is the economic cost of an error of 0.5% in the flow rate for a pipeline that carries 100,000 barrels/day of crude oil? (assume the error is not in your favour, i.e. it really is a cost, and not a profit)

.. question:: :grading: 20

Consider the following tank system:

.. image:: ../figures/process-control/CLES/tank-system-CLES.png :width: 500

#. Write out the mass balance for the tank in the form:

.. math::

\text{rate of accumulation of mass} = \text{rate of mass flow in} - \text{rate of mass flow out}

so that you obtain a differential equation relating height :math:`h`, flows :math:`F_i` and :math:`F_0`, and the cross-sectional area, :math:`A`. The flows are considered to be volumetric flows.

#. Let the flows be :math:`F_i = 10\,\text{m}^3\text{.min}^{-1}` and :math:`F_0 = 11.5\,\text{m}^3\text{.min}^{-1}`, the area is :math:`A = 2\,\text{m}^2`, and the current height in the tank is 3 m.

Solve the differential equation by hand to determine an analytical expression of the tank height as a function of time, i.e. :math:`h(t) = ...`

#. Sketch your function over time, for a duration of 3 minutes.

#. What is the tank level at :math:`t=2` minutes?

#. Create a simulation of the system in MATLAB, and integrate the ODE using the MATLAB solver ``ode45(...)``, which you learned about in Chemical Engineering ``3E04``. Print and attach your plot, showing that your simulation matches the solution to part 3 and 4 above.

#. Now we add a simple controller to the system, a *proportional controller*.

.. image:: ../figures/process-control/CLES/tank-system-with-proportional-control-CLES.png :width: 500

A proportional controller modifies the valve position, *in proportion to* the error, or deviation from set point. The set point is :math:`h_s = 2.0\,\text{m}`. When there is no deviation, the controller makes no change to the valve. If there is a high deviation from set point, it makes a greater change. The control algorithm can be represented as:

.. math::

F_0(t) = F_{0,s} + K_c(h - h_s)

where the error term is :math:`h - h_s`, and the controller tuning constant is :math:`K_c = 8\,\text{m}^2\text{.min}^{-1}`. The term :math:`F_{0,s} = 10\,\text{m}^3\text{.min}^{-1}` represents the amount of flow leaving the tank when there is no error (i.e. when the height is equal to the set point).

Add this equation to your simulation from part 5 (either add it as a new equation - *preferred* - or substitute it into the previous equation).

Now re-simulate the height in the tank over time for a period of 3 minutes. What do you notice about your simulation, compared to part 5?

Print out and attach your simulation.

</rst>