Chemical engineering transport & equilibrium foundations of Gravity-Drained Tank Dynamics.
Liquid level control in surge tanks, flash drums, reflux accumulators, and chemical storage reservoirs is a ubiquitous unit operation in process engineering. Maintaining stable tank levels prevents pump cavitation from run-dry conditions while avoiding hazardous overfills. A vertical, uniform cross-sectional surge tank (cross-sectional area A = 0.050 m^2) drained by gravity through a laminar restriction tube is the canonical physical realization of a first-order dynamic process. Applying the unsteady-state mass conservation principle to incompressible liquid (constant density ρ) gives: ρ · A · (dh/dt) = ρ · F_in(t) - ρ · F_out(t), which simplifies to: A · (dh/dt) = F_in(t) - F_out(t). When discharge occurs through a capillary tube or packed porous plug under laminar flow (Re < 2000), the Hagen–Poiseuille equation dictates that effluent flow rate is directly proportional to hydrostatic pressure head: F_out = h / R, where R is the linear hydraulic discharge resistance (with dimensions of head per volumetric flow rate, cm/(L/min)). Substituting this linear outflow relation yields the first-order ordinary differential equation: A · R · (dh/dt) + h(t) = R · F_in(t). Defining the process time constant τ_p = A · R, this standard first-order transfer function in Laplace deviation variables is: G(s) = H'(s) / F_in'(s) = K_p / (τ_p · s + 1), where steady-state process gain K_p equals hydraulic resistance R. In response to a step increase in inflow ΔF_in, the level response is an exponential relaxation: h'(t) = R · ΔF_in · [1 - exp(-t / τ_p)]. At long times (t -> inf), the new steady-state liquid level is h_ss = R · F_in. In this laboratory investigation, an instrumented transparent acrylic surge tank equipped with an ultrasonic level transmitter and variable-speed magnetic gear pump is operated across inflow setpoints from 2.0 to 12.0 L/min to evaluate transient step responses, calculate τ_p, and determine the linear hydraulic resistance R = 3.8500 cm/(L/min).
Dynamic Mass Conservation & Lumped-Parameter Tank Modeling
The dynamic behavior of single-phase liquid storage vessels is modeled by macroscopic mass conservation. For a constant-density incompressible liquid ($\rho = \text{const}$) inside a vertical cylinder with uniform cross-sectional area $A$:
\[ \frac{d(\rho A h)}{d t} = \rho \, F_{\text{in}}(t) - \rho \, F_{\text{out}}(t) \implies A \, \frac{d h}{d t} = F_{\text{in}}(t) - F_{\text{out}}(t) \]In laminar capillary or porous plug discharge, the effluent volumetric flow rate $F_{\text{out}}$ follows the Hagen–Poiseuille equation, scaling linearly with hydrostatic driving head $h(t)$:
\[ F_{\text{out}}(t) = \frac{h(t)}{R} \]where $R$ is the linear hydraulic discharge resistance ($R = \frac{128 \mu L}{\pi \rho g d^4}$). Substituting into the mass balance yields the first-order linear ordinary differential equation:
\[ A \, \frac{d h}{d t} + \frac{h(t)}{R} = F_{\text{in}}(t) \implies A R \, \frac{d h}{d t} + h(t) = R \, F_{\text{in}}(t) \]Defining the process time constant $\tau_p = A \, R$ and steady-state process gain $K_p = R$:
\[ \tau_p \, \frac{d h}{d t} + h(t) = K_p \, F_{\text{in}}(t) \]Transfer Function & Analytical Step Response
Transforming into Laplace deviation variables ($h'(t) = h(t) - h_{\text{initial}}$, $F_{\text{in}}'(t) = F_{\text{in}}(t) - F_{\text{in,initial}}$):
\[ (\tau_p s + 1) H'(s) = K_p F_{\text{in}}'(s) \implies G(s) = \frac{H'(s)}{F_{\text{in}}'(s)} = \frac{K_p}{\tau_p s + 1} = \frac{R}{A R s + 1} \]For an instantaneous step increase of magnitude $\Delta F_{\text{in}}$ at $t = 0$ ($F_{\text{in}}'(s) = \Delta F_{\text{in}} / s$):
\[ H'(s) = \frac{R \, \Delta F_{\text{in}}}{s (\tau_p s + 1)} \implies h'(t) = R \, \Delta F_{\text{in}} \left( 1 - e^{-t / \tau_p} \right) \]At $t = \tau_p$, $h'(\tau_p) = R \, \Delta F_{\text{in}} (1 - e^{-1}) = 0.6321 \, \Delta h_{\text{ss}}$. At $t = 3 \tau_p$, the response reaches $95.0\%$, and at $t = 5 \tau_p$, $99.3\%$.
Steady-State Hydraulic Resistance Scaling
At steady state ($d h / d t = 0$), accumulated mass is zero, requiring $F_{\text{out}} = F_{\text{in}}$. Therefore:
\[ h_{\text{ss}} = R \cdot F_{\text{in}} = 3.8500 \cdot F_{\text{in}} \]Plotting the steady-state liquid level $h_{\text{ss}}$ versus influent flow rate $F_{\text{in}}$ yields a straight line through the origin with slope equal to the hydraulic resistance $R = 3.8500\text{ cm/(L/min)}$.
Viva voce questions & analytical theory
Why is the discharge flow rate through an orifice non-linear while a capillary tube is linear?
Flow through a sharp-edged orifice is turbulent and inertia-dominated, governed by Torricelli's law F_out = C_d A_o sqrt(2 g h) (proportional to h^0.5); in contrast, flow through a long, narrow capillary tube is viscous and laminar, governed by the Hagen–Poiseuille law F_out = (π ρ g d^4 / (128 μ L)) · h (linearly proportional to h).
What is the physical interpretation of the process time constant τ_p in a liquid level system?
The time constant τ_p = A · R represents the characteristic capacitive-resistive relaxation time of the vessel; physically, it is the time required for liquid level to complete 63.2% of its ultimate transition following an instantaneous step change in feed flow.
How does increasing tank cross-sectional area A affect system dynamics and steady state?
Increasing cross-sectional area A increases the liquid capacitance and time constant τ_p = A · R, slowing down dynamic response and filtering high-frequency flow disturbances; however, because steady state is governed solely by resistance (h_ss = R · F_in), tank area has zero effect on the final steady-state level.
How are deviation variables defined and why are they used in process dynamics?
A deviation variable is the difference between an instantaneous process variable and its initial steady-state value: h'(t) = h(t) - h_ss; using deviation variables eliminates constant terms, allowing differential equations to transform directly into linear Laplace transfer functions with zero initial conditions.
What is self-regulation in a process control context?
A self-regulating process naturally converges to a new steady state following an unmeasured disturbance without requiring active feedback control; a gravity-drained tank is self-regulating because higher level increases hydrostatic discharge pressure until outflow naturally re-balances inflow.
How does liquid temperature and viscosity affect the discharge resistance R?
In laminar capillary discharge, R = 128 μ L / (π ρ g d^4); because water viscosity μ decreases significantly with rising temperature (~2.5% per °C), warmer water reduces hydraulic resistance R, resulting in lower steady-state liquid levels for the same volumetric inflow.
Experimental precautions and process safety
- Ensure the discharge capillary manifold is purged of all air bubbles prior to testing, as trapped air pockets artificially elevate discharge resistance.
- Never allow the liquid level to exceed 85 cm to prevent overflowing into the overflow weir and flooding the laboratory floor.
- Verify that ultrasonic level sensor transducer face is clean and mounted perpendicularly to the liquid surface to avoid echo distortion.
- Maintain water sump temperature within ±0.5 °C to avoid viscosity-induced shifts in laminar hydraulic resistance.
Theoretical assumptions & boundary condition limits
- The transport and thermodynamic models apply within the specified calibrated operating window (2.0 to 12.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.