Chemical engineering transport & equilibrium foundations of Two-Tank Interacting Dynamics.
Multi-capacity processes in chemical and process engineering—including multi-stage distillation columns, series heat exchangers, reactor cascades, and multi-compartment separators—exhibit complex dynamic responses resulting from hydraulic and thermodynamic coupling between individual stages. When two liquid storage tanks are connected such that flow between them depends on the difference between both liquid levels, the system is classified as an interacting two-tank process. In an interacting configuration, influent enters Tank 1 (area A_1 = 400.0 cm^2), flowing into Tank 2 (area A_2 = 400.0 cm^2) through an inter-tank laminar restriction of resistance R_1. Because inter-tank flow is driven by the differential hydrostatic head: F_12 = (h_1 - h_2) / R_1, the liquid level in Tank 2 exerts a direct backpressure on Tank 1, coupling both mass balance differential equations: A_1 · (dh_1/dt) = F_1(t) - (h_1 - h_2)/R_1, and A_2 · (dh_2/dt) = (h_1 - h_2)/R_1 - h_2/R_2, where R_2 is the discharge resistance of Tank 2. Transforming these coupled equations into Laplace deviation variables yields the canonical second-order overdamped transfer function: G_2(s) = H_2'(s) / F_1'(s) = K_p / (τ_1 · τ_2 · s^2 + (τ_1 + τ_2 + A_1 · R_2) · s + 1), where individual capacities are τ_1 = A_1 · R_1 and τ_2 = A_2 · R_2. Crucially, the interaction cross-term A_1 · R_2 in the linear damping coefficient increases the effective time constants, making the interacting system more sluggish (higher damping ratio ζ > 1) than an equivalent non-interacting cascade where Tank 1 discharges freely into Tank 2 through an open air gap. However, at steady state (s -> 0), all time derivatives vanish, and mass conservation requires that total effluent from Tank 2 matches primary influent: F_out = h_2 / R_2 = F_1. Consequently, the steady-state process gain K_p = Δh_2 / ΔF_1 equals the final discharge resistance R_2 = 4.2500 cm/(L/min), completely independent of inter-tank resistance R_1. In this laboratory experiment, an instrumented dual-tank benchtop trainer operates under interacting mode across feed flow step changes from 1.0 to 6.0 L/min to evaluate transient S-shaped overdamped response curves and determine K_p = 4.2500 cm/(L/min).
Dynamic State-Space Modeling of Interacting Capacitive Stages
In process systems, multi-capacity units coupled via hydraulic or chemical gradients exhibit interacting dynamics. Consider two vertical cylindrical surge vessels of cross-sectional areas $A_1$ and $A_2$. Feed stream $F_1(t)$ enters Tank 1. Liquid flows from Tank 1 to Tank 2 through an inter-tank capillary pipe of linear hydraulic resistance $R_1$, driven by differential hydrostatic head:
\[ F_{12}(t) = \frac{h_1(t) - h_2(t)}{R_1} \]Tank 2 discharges to ambient through a secondary capillary resistor of resistance $R_2$:
\[ F_{\text{out}}(t) = \frac{h_2(t)}{R_2} \]Applying unsteady-state mass conservation to each tank:
\[ A_1 \, \frac{d h_1}{d t} = F_1(t) - \frac{h_1(t) - h_2(t)}{R_1} \] \[ A_2 \, \frac{d h_2}{d t} = \frac{h_1(t) - h_2(t)}{R_1} - \frac{h_2(t)}{R_2} \]Laplace Transform & Characteristic Polynomial
Defining deviation variables $h_1'(t) = h_1(t) - h_{1,\text{ss}}$, $h_2'(t) = h_2(t) - h_{2,\text{ss}}$, and $F_1'(t) = F_1(t) - F_{1,\text{ss}}$, the Laplace-transformed system equations are:
\[ \left( A_1 s + \frac{1}{R_1} \right) H_1'(s) - \frac{1}{R_1} H_2'(s) = F_1'(s) \] \[ -\frac{1}{R_1} H_1'(s) + \left( A_2 s + \frac{1}{R_1} + \frac{1}{R_2} \right) H_2'(s) = 0 \]Solving for the transfer function relating Tank 2 level to primary feed flow:
\[ G_2(s) = \frac{H_2'(s)}{F_1'(s)} = \frac{R_2}{\tau_1 \tau_2 s^2 + (\tau_1 + \tau_2 + A_1 R_2) s + 1} \]where $\tau_1 = A_1 R_1$ and $\tau_2 = A_2 R_2$ are the uncoupled time constants. The interaction term $A_1 R_2$ appears in the denominator, increasing the linear damping coefficient $\zeta > 1$ and separating the two real poles:
\[ G_2(s) = \frac{K_p}{(\tau_a s + 1)(\tau_b s + 1)} \]where $\tau_a + \tau_b = \tau_1 + \tau_2 + A_1 R_2 > \tau_1 + \tau_2$, demonstrating that interacting coupling slows the dynamic response.
Steady-State Coupled Process Gain
In the limit as $t \to \infty$ ($s \to 0$), all dynamic accumulation terms vanish. Steady-state mass conservation across the entire two-tank envelope requires that net discharge from Tank 2 matches net feed into Tank 1:
\[ F_{\text{out,ss}} = \frac{\Delta h_2}{R_2} = \Delta F_1 \implies \Delta h_2 = R_2 \cdot \Delta F_1 \]Thus, the overall steady-state interacting process gain is governed strictly by the final discharge resistance:
\[ K_p = \frac{\Delta h_2}{\Delta F_1} = R_2 = 4.2500\text{ cm/(L/min)} \]Viva voce questions & analytical theory
What is the physical difference between an interacting and non-interacting two-tank system?
In a non-interacting system, Tank 1 discharges freely into Tank 2 through an air gap, meaning level changes in Tank 2 exert zero backpressure on Tank 1 (one-way information flow); in an interacting system, the tanks are connected below the liquid line, so inter-tank flow depends on (h_1 - h_2), creating two-way dynamic coupling.
Why does an interacting system exhibit more sluggish dynamics than a non-interacting system?
The inter-tank backpressure creates an extra interaction term A_1 · R_2 in the first-order coefficient of the characteristic polynomial: τ_1 τ_2 s^2 + (τ_1 + τ_2 + A_1 R_2) s + 1; this increases the effective damping ratio ζ > 1, shifting the dominant pole closer to the imaginary axis and slowing response time.
Why is the steady-state gain K_p identical for both interacting and non-interacting systems?
At steady state, all liquid entering Tank 1 must eventually discharge from Tank 2 (inflow equals outflow); because outflow from Tank 2 is governed purely by its own discharge resistor F_out = h_2 / R_2, the final level is h_2 = R_2 · F_1, meaning steady-state gain is K_p = R_2 regardless of interaction.
Why does the step response of Tank 2 exhibit an S-shaped inflection point?
Tank 2 is a second-order system; initially, when Tank 1 inflow steps up, Tank 1 level must first rise before inter-tank flow increases, resulting in zero initial slope (dh_2/dt = 0 at t = 0) followed by acceleration and deceleration, producing an S-shaped sigmoidal response.
How can the state-space representation be formulated for the interacting two-tank process?
Defining state vector x = [h_1', h_2']^T and input u = F_1', the state equations are: dx_1/dt = -x_1/(A_1 R_1) + x_2/(A_1 R_1) + u/A_1, and dx_2/dt = x_1/(A_2 R_1) - x_2 · [1/(A_2 R_1) + 1/(A_2 R_2)]; the state matrix A has negative real eigenvalues confirming asymptotic stability.
What happens if the inter-tank valve R_1 is completely closed?
The coupling between tanks is severed: Tank 1 becomes an isolated first-order tank that eventually overflows, while Tank 2 drains to zero level because it receives zero influent.
Experimental precautions and process safety
- Ensure the inter-tank connecting line is completely submerged and purged of air bubbles to prevent airlock disruption of inter-tank flow.
- Verify that both tanks' overflow weirs are connected to the central drainage trough to protect electronics from accidental level surges.
- Calibrate both bottom piezoresistive pressure transmitters simultaneously using identical hydrostatic liquid depths.
- Allow adequate settling time (at least 5 effective time constants, ~8 minutes) for complete steady-state equilibration before recording final levels.
Theoretical assumptions & boundary condition limits
- The transport and thermodynamic models apply within the specified calibrated operating window (1.0 to 6.0 L/min).
- Assumes ideal continuum flow and homogeneous phase mixing without unmodeled parasitic wall losses.
- Continuous sensor readings reflect ordinary least-squares regression across successive observations.
- High-temperature or extreme concentration excursions require non-ideal equation-of-state corrections.