23-Chem-A6 Process Dynamics and Control · May 2016
Question 4 of 8: Two Concentric Tanks — Second-Order Thermal Response
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Reference texts: D. E. Seborg, T. F. Edgar, D. A. Mellichamp & F. J. Doyle III, Process Dynamics and Control (4th ed., Wiley) — Laplace-domain modelling, linearisation, transfer functions, feedback stability, frequency response and IMC design; G. Stephanopoulos, Chemical Process Control: An Introduction to Theory and Practice (Prentice Hall) — Routh test, Nyquist criterion and dead-time systems; D. R. Coughanowr & S. E. LeBlanc, Process Systems Analysis and Control (3rd ed., McGraw-Hill) — first-order thermal/level dynamics, step/ramp/pulse response and block-diagram algebra. Standard control conventions (deviation variables; unity valve/sensor gains unless stated) are used throughout. This sitting is an open-book paper of eight problems of which any five constitute a complete exam; all eight are solved here.
Problem 4: Two Concentric Tanks — Second-Order Thermal Response (20%)
Given. Two well-mixed tanks coupled by a conducting wall:
Quantity
Symbol
Value
Thermal capacitance each tank
$\rho V c_p$
$50\cdot1\cdot2=100\ \mathrm{Btu/{}^\circ F}$
Wall conductance
$UA$
$10\cdot1=10\ \mathrm{Btu/(hr\,{}^\circ F)}$
Outer-flow capacity
$w c_p$
$10\cdot2=20\ \mathrm{Btu/(hr\,{}^\circ F)}$
Step heat input (inner)
$Q$
$0\to500\ \mathrm{Btu/hr}$
Feed & initial temperatures
$T_f,T_{1,0},T_{2,0}$
$100\ ^\circ$F
Find. (a) $T_1(s)$ (inner-tank deviation temperature); (b) $T_1$ at $t=0$ and $t=5$ hr.
Problem 4: the immersed heater adds $Q$ to the inner tank; heat conducts through the wall to the outer tank ($UA(T_1-T_2)$), which is swept by the through-flow $w$. Two capacitances in a feedback path give a second-order response.
Approach. Write an energy balance on each tank in deviation variables (feed constant, so its deviation is zero), Laplace-transform the pair, eliminate the outer-tank temperature to get $T_1(s)$, then invert by partial fractions and evaluate at the two times.
Energy balances (deviation form). With capacitance $100$, wall term $10$, flow term $20$ and $Q$ a step of $500$: $$100\frac{dT_1}{dt}=Q-10(T_1-T_2),\qquad 100\frac{dT_2}{dt}=-20T_2+10(T_1-T_2).$$ (Primes for deviation dropped; $T_f'=0$.)
Transform. $$(100s+10)T_1-10T_2=Q,\qquad -10T_1+(100s+30)T_2=0.$$ From the second equation $T_2=\dfrac{10\,T_1}{100s+30}$.
(a) Eliminate $T_2$. Substituting and clearing, $(100s+10)(100s+30)-100=10000s^2+4000s+200$, so with $Q=500/s$: $$\boxed{T_1(s)=\frac{500\,(100s+30)}{s\,(10000s^{2}+4000s+200)}=\frac{25\,(10s+3)}{s\,(50s^{2}+20s+1)}.}$$
Poles. $50s^2+20s+1=0$ gives $s=-0.0586,\,-0.3414\ \mathrm{hr^{-1}}$ (time constants $17.1$ and $2.93$ hr) — two real, over-damped roots, as expected for a two-capacitance thermal system.
(b) Evaluate. At $t=0$: $75-72.86-2.14=0$, so $\boxed{T_1(0)=100,{}^ rc\text{F}}$ (matches the initial condition). At $t=5$ hr: $75-72.86e^{-0.293}-2.14e^{-1.707}=75-54.36-0.39=20.25$, so $\boxed{T_1(5)=120.3,{}^ rc\text{F}}$. The steady value is $175,{}^ rc$F, confirmed by the overall balance ($Q=wc_p\,\Delta T_2\Rightarrow T_2=125$; $Q=UA\,\Delta T\Rightarrow T_1=175$).
Problem 4(b): the inner tank climbs smoothly (over-damped, no overshoot) from $100,{}^ rc$F toward the steady $175,{}^ rc$F; at $t=5$ hr it has reached $120.3,{}^ rc$F. The slow $17.1$ hr mode dominates the long tail.