Tutorial 5 - 2013

Due date(s):	All questions due as assignment 3B
(PDF)	Tutorial questions
(PDF)	Assignment solutions

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

=========

• • Objectives of this short assignment**:

We are going to start using software more frequently in the course, because the equations we deal with are too difficult to integrate analytically. So the main purpose of this tutorial is:

* To install and use a numerical integration software package on your computer. You may use any software, but Python, MATLAB and Polymath are recommended candidates. * Ensure you can use the software to successfully integrate one or more differential equations.

.. question::

A can of your favourite beverage, initially at 25°C, is placed in the fridge, where the ambient temperature is 5°C. Plot the temperature profile of the can's contents (assume the contents are well mixed) over time. The temperature of the can's contents, :math:`T`, can be derived and modelled as:

.. math::

\dfrac{dT}{dt} = -K (T - T_f)

where :math:`t` is the time, in minutes, :math:`K` is the heat transfer coefficient that lumps the can's properties into a single variable, and :math:`T_f` is the fridge reference temperature.

#. Plot the temperature profiles for :math:`K = 0.03\,\text{min}^{-1}`, :math:`K = 0.09\,\text{min}^{-1}` and :math:`K = 0.15\,\text{min}^{-1}` on the same plot, on a range to 1 hour.

#. Repeat the same exercise, but use a logarithmic :math:`y`-axis. Comment on the results from this part, as compared to the previous part.

.. question::

This question is a good check to ensure you are using the software correctly: you must be able to get the same answer with the software that you get analytically.

A batch bioreactor is charged with substrate at time :math:`t=0` to start a reaction that will consume the substrate. The reaction mechanism is too complex to model exactly, but previous experience suggests it is roughly a first-order reaction, with apparent rate expression is :math:`-r_S = k_S C_S`, where :math:`k_S = 0.58\,\text{hour}^{-1}`, and :math:`C_S` is the concentration of the substrate being depleted.

#. Show that the rate of change of the concentration of S is given by:

.. math::

\frac{dC_S}{dt} = -k_S C_S(t)

#. The time required to consume the substrate from a starting level of :math:`260\,\text{mol.L}^{-1}` to :math:`5\,\text{mol.L}^{-1}` can be calculated by analytical integration of the above expression. What is the total time required?

#. The tank volume is doubled; how long will the duration be to go from the same starting to final substate concentration?

#. Use software tools and integrate between the same limits and calculate the result numerically. Verify that your result matches the analytical answer.

Show the software source code you wrote to find the solution.

.. question::

From the midterm, the solution to the last question was to write the ODEs:

.. math::

\dfrac{dX}{dW} &= \dfrac{-r_A}{F_{A0}} = \dfrac{(1 \times 10^{-5})(30.07)^2 \dfrac{(1 - X)^2}{(1+\varepsilon X)^2} \left(y\right)^2}{ 0.005\,\text{mol.s}^{-1}} = 1.808\dfrac{(1 - X)^2(y)^2 }{(1+ 0.5X)^2} \\ \dfrac{dy}{dW} &= -\dfrac{\alpha}{2y}(1+\varepsilon X) = -\dfrac{5.17 \times 10^{-4}}{2y}(1+0.5 X) \\ X(W=0) &= 0.0 \\ y(W=0) &= 1.0 \\

#. Integrate the solution from :math:`W=0` zero to :math:`W = 235.6\,\text{kg}` and include the profiles of conversion and :math:`y` (on the same graph) as your answer. #. What is the pressure leaving the reactor? #. What is the conversion leaving the reactor? #. What is the conversion midway along the reactor's length?

*Note*: in the course project and future questions, you will NOT simplify the equations. You will write multiple equations for the knowns and unknowns and let the software do all the work for you.

</rst>