Assignment 2 - 2014

Due date(s):	27 January 2014, in class
(PDF)	Assignment questions
(PDF)	Assignment solutions

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

.. rubric:: Assignment objectives: gaining an excellent understanding of dynamic systems; using computer models and Laplace transforms to infer and describe the model behaviour over time.

.. question:: :grading: 10

For each transfer function below, what can you say about the time-domain function, :math:`H(t)`, where :math:`0 \leq t \leq \infty`? Do not solve for :math:`H(t)` analytically.

Your answers should include the initial value, final value, and a description of the stability, smoothness, and oscillatory behaviour.

#. :math:`H(s) = \dfrac{6(s+2)}{(s^2 + 9s + 20)(s+4)}` #. :math:`H(s) = \dfrac{10s^2 -3 }{(s^2 - 6s + 10)(s+2)}` #. :math:`H(s) = \dfrac{16s+5}{s^2+9}`

.. question:: :grading: 4

Find the inverse transform of:

.. math::

Y(s) = \dfrac{s+2}{s(s+1)^2}

.. question:: :grading: 20

We want to determine the dynamic response of the liquid level in the open tank shown in the figure.

.. image:: ../figures/process-control/Marlin-slides/Marlin-pump-tank.png :scale: 50 :width: 700px

The flow entering is determined from an upstream process (we cannot control that flow rate). The flow out is determined by the valve percent opening, with a constant speed, centrifugal pump. (The level in the tank in this system does not significantly affect the flow out because of the high pressure supplied by the pump.)

Some data and details are given in the following:

* Initial flow in = :math:`2.0\,\text{m}^3\text{.min}^{-1}` * Initial flow out = :math:`2.0\,\text{m}^3\text{.min}^{-1}` * Total tank volume (when full) = :math:`10.0\,\text{m}^3` * Initial tank level = 50% full * Tank cross-sectional area = :math:`3.0\,\text{m}^2` * Assume that the tank is initially at steady state with the flows in and out equal.

A step change is introduced in the inlet flow; it changes from :math:`2.0\,\text{m}^3\text{.min}^{-1}` to :math:`1.0\,\text{m}^3\text{.min}^{-1}`. No change is made to the valve opening.

#. Without using a computer or analytically derived expression, draw the plot of how you expect the liquid level to change once the step change is made. [There are no grades for this part of the question, so don't modify your graph later on].

#. Determine the actual dynamic response of the liquid level using a computer simulation.

#. Based on your graph and numerical results, what do you conclude about whether the level should be controlled using the feedback principle by adjusting the valve in the output after the pump?

#. The liquid is water at 25°C, and we want to measure the level within ± 1 percent of its true value. What level sensor would you recommend? Give a brief explanation with a citation for information used in your answer. [Hint: use information available at http://pc-education.mcmaster.ca/ and follow links to the Instrumentation section].

#. The flow rate out of the tank must be measured with an accuracy of ± 1 percent from about 0.25 to :math:`2.0\,\text{m}^3\text{.min}^{-1}`. What flow sensor would you recommend? Give a brief explanation with a citation for information used in your answer.

.. question:: :grading: 12

A heat exchanger system has a differential equation below that shows how the temperature changes with a change in flow rate, :math:`q(t)`.

.. math::

2 \dfrac{dT}{dt} = - T + 5 q

The temperature, :math:`T` and :math:`q` are in deviation form already.

#. :math:`q(t)` is changed from 0 to 2.0 at time :math:`t=0`. Sketch the response up to the point where the temperature reaches its new steady state.

#. What is this steady state?

#. How long does it take to reach within 0.1 degree of the new steady state.

#. If instead :math:`q(t)` is changed from 0 to 4.0 at time :math:`t=0`. Sketch the response again, and report the new steady state value of the temperature.

#. Superimpose the plots from part 1 and part 4 of this question and explain why the curves have the same shape.

.. question:: :grading: 15

In this question, you will reinforce the material learned in Chemical Engineering ``3E04``. The system that you will consider is given in Figure 4. It is a CSTR where the reaction :math:`\text{A} \longrightarrow \text{B}` is occurring. You do not have to derive the component material balance, i.e., the differential equation in the figure.

.. image:: ../figures/process-control/Marlin-slides/ode-system.png :scale: 100 :width: 700px

#. First, you must determine the *initial* value for the reactor concentration (which is the same as the concentration leaving the reactor), :math:`C_A`. The system is initially at steady state, with :math:`C_{A,0} = 0.5\,\text{mol.m}^{-3}`.

You may solve this part analytically or numerically. *Hint*: you will solve a quadratic equation.

#. Starting from this steady state, determine the dynamic response of the concentration :math:`C_A` in the reactor if a step increase of :math:`C_{A,0}` is made from :math:`C_{A,0} = 0.5\,\text{mol.m}^{-3}` to :math:`C_{A,0} = 1.0\,\text{mol.m}^{-3}`.

For example, the above occurs when your manager has asked you to increase the concentration of species B that you are producing.

#. How long does it take for :math:`C_{A}(t)` to reach a new steady state?

#. You want to reach the new steady state faster. Explain how you can decrease the time taken. Perform a new simulation with the change(s) implemented, and superimpose the original time-domain plot from part 2 to show how much faster your response is.

*Work ahead*: in the next tutorial and assignment you will create deviation variables :math:`C_A'` for outlet concentration, and :math:`C_{A0}'` for inlet concentration. Create a linearized model, invert it with the Laplace transform, and compare the linearized (approximate) model to the actual model.

</rst>