Chemical & Process Engineering / Two-Tank Interacting Dynamics
FREE TO EXPLORE
LABORATORY 41 · CHEMICAL & PROCESS ENGINEERING

Two-Tank Interacting Dynamics lab report

A worked two-tank interacting dynamics lab report with sample observation tables, regression plots, calculations, and error analysis.

WORKED EXAMPLE · BENCHMARK OBSERVATIONS

A complete example chemical engineering laboratory report.

Simulated demonstration data.

This report illustrates formal chemical and process engineering data reduction, instrument calibration tares, and unit operations parameter estimation. Replace its observations with your own experimental telemetry when preparing a chemical engineering laboratory report submission.

Student: __________________Student ID: __________________Course: __________________Date: __________________Partner: __________________Instructor: __________________

1. Title and aim

Interacting vs Non-Interacting Two-Tank Liquid Level Dynamics & Coupled State-Space Modeling. Operate a two-tank interacting liquid level system with inter-tank flow coupling, apply step perturbations in primary tank feed flow ΔF_1 from 1.0 to 6.0 L/min, monitor coupled liquid level responses h_1(t) and h_2(t), evaluate second-order overdamped dynamics, and determine the steady-state interacting process gain K_p.

2. Abstract

Six controlled chemical process engineering observations were systematically evaluated across operational settings Feed Perturbation ΔF_1 = 1.0 L/min to Feed Perturbation ΔF_1 = 6.0 L/min for Interacting vs Non-Interacting Two-Tank Liquid Level Dynamics & Coupled State-Space Modeling. Ordinary least-squares regression yielded an empirical Interacting Coupled Process Gain K_p of 4.2500 cm/(L/min) (difference from theoretical reference = 0.0000%). The empirical telemetry conclusively validates the fundamental chemical engineering transport and equilibrium formulation with R² = 1.000000. The investigation demonstrates stable unit operations behavior, hydrodynamic and thermal reproducibility, and strict adherence to classical chemical engineering principles.

3. Apparatus and fixed settings

4. Theory and graph choice

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).

Δh_2 = K_p · ΔF_1 = 4.2500 · ΔF_1

The governing phenomenological transport or equilibrium relationship was linearized into the coordinate form y = mx + b. Here, the independent coordinate x represents Tank 1 feed perturbation ΔF_1 (L/min), and the dependent coordinate y represents Tank 2 level deviation Δh_2 (cm). The gradient m directly determines Interacting Coupled Process Gain K_p, while the intercept b accounts for baseline tare values or reference potentials. Linear transformation enables robust parameter estimation via ordinary least squares without non-linear numerical divergence. Plotting transformed observation pairs displays process discharge coefficients, heat transfer multipliers, mass transfer rates, or kinetic rate constants directly as linear slopes.

5. Method

  1. Fill both tanks with water; drain completely to baseline; calibrate piezoresistive transmitter zeroes to bottom datum.
  2. Open inter-tank coupling valve R_1 completely; verify submerged liquid connection between tanks to enforce interacting mode.
  3. Establish baseline feed flow F_1 = 3.0 L/min; allow levels h_1 and h_2 to stabilize at steady state for 10 minutes.
  4. Apply step increase ΔF_1 = 1.0 L/min (F_1 = 4.0 L/min); record dynamic trajectories h_1(t) and h_2(t) at 2 Hz until new steady state.
  5. Step feed perturbations through ΔF_1 = 2.0, 3.0, 4.0, 5.0, and 6.0 L/min, recording S-shaped second-order step responses for each plateau.
  6. Evaluate steady level deviations Δh_2; plot Δh_2 versus ΔF_1 and perform linear regression to determine interacting process gain K_p.

6. Observations and calculation table

All chemical process observations use strictly controlled 64-bit IEEE 754 precision buffers. Residuals are computed relative to the ordinary least-squares regression line.

Table 1. Six simulated observations; x and y represent transformed process metrics.
No.Tank 1 Flow Perturbation ΔF_1 (L/min)Tank 2 Steady Level Deviation Δh_2 (cm)Tank 1 feed perturbation ΔF_1 (L/min)Tank 2 level deviation Δh_2 (cm)Fitted yResidual
11.04.251.00004.25004.25000.0000
22.08.52.00008.50008.50000.0000
33.012.753.000012.750012.75000.0000
44.017.04.000017.000017.00000.0000
55.021.255.000021.250021.25000.0000
66.025.56.000025.500025.50000.0000

Download the sample observations (CSV) ↓

7. Worked calculation for one observation

For input Feed Perturbation ΔF_1 = 1.0 L/min: Measured Tank 2 Steady Level Deviation Δh_2 = 4.2500 cm. Transformed process coordinates: independent x = 1.0000, dependent y = 4.2500.

8. Graph, gradient and final calculation

1.060.407.971.9514.883.5021.785.0528.696.60Tank 1 feed perturbation ΔF_1 (L/min)Tank 2 level deviation Δh_2 (cm)

Circular markers represent empirical process observations. The solid line represents the ordinary least-squares linear fit.

x̄ = 3.5000; ȳ = 14.8750
Σ(x − x̄)² = 17.5000
Σ(x − x̄)(y − ȳ) = 74.3750
m = Σ(x − x̄)(y − ȳ) / Σ(x − x̄)² = 4.25000000
b = ȳ − mx̄ = 0.00000000
y = 4.2500 · x + 0.0000 ⟹ Target Interacting Coupled Process Gain K_p = 4.2500 cm/(L/min) (theoretical = 4.2500)

9. Uncertainty and comparison

Percentage difference = |estimate − theoretical value| / |theoretical value| × 100 = 0.0000%.

For n = 6 observations, Σresidual² = 0.0000; s² = Σresidual²/(n − 2). Gradient standard error = √[s²/Σ(x − x̄)²] = 0.0000. Intercept standard error = √[s²(1/n + x̄²/Σ(x − x̄)²)] = 0.0000.

Propagating the fit parameters yields a fit standard error of 1.000e-05 cm/(L/min), with an ideal coefficient of determination R² = 1.000000.

In chemical and process engineering experimentation, unit operations verification, and pilot-scale process modeling, empirical uncertainties arise from fluid property variations, slight thermal gradients across heat transfer boundaries, transducer electronic drift (64-bit IEEE 754 precision digital acquisition), and minor flow turbulence eddies. The ordinary least-squares regression across n = 6 standardized test intervals yields unbiased minimum-variance estimators for both slope and intercept. Because the underlying transport phenomena—including Navier-Stokes momentum balances, Fourier conduction, Fickian mass diffusion, and reaction kinetics—adhere strictly to macroscopic conservation laws, residual dispersions remain tightly bounded, yielding an ideal coefficient of determination R² ≥ 0.999999 and confirming the theoretical formulation. Systematic errors (such as sensor zero offsets, thermocouple reference junction drifts, and flowmeter discharge friction) are eliminated using certified calibration standards and differential tare subtraction. Uncertainties in derived transport coefficients, dimensionless numbers, and kinetic constants are evaluated via multivariate error propagation.

10. Result and discussion

Interacting Coupled Process Gain K_p = 4.2500 cm/(L/min)
Theoretical reference value: 4.2500 cm/(L/min)

The empirical laboratory observations provide rigorous quantitative verification of Interacting vs Non-Interacting Two-Tank Liquid Level Dynamics & Coupled State-Space Modeling. In chemical engineering transport phenomena, fluid mechanics, heat transfer, separation processes, and reaction engineering, physical systems are governed by macroscopic mass, momentum, and energy conservation equations. Across all six observation trials, the empirical slope m = 4.2500 confirms the theoretical model with near-zero residual variance. Furthermore, comparison with established chemical engineering literature benchmarks (McCabe, Smith, Harriott; Bird, Stewart, Lightfoot; Fogler; and Levenspiel) underscores vital industrial scale-up principles including continuous steady-state stability, minimal friction losses, optimal reaction conversion, and thermal efficiency. Residual analysis confirms homoscedastic dispersion without systematic non-linear curvature, validating the underlying model assumptions.

11. Precautions and limitations

The transport and thermodynamic models apply within the specified calibrated operating window (1.0 to 6.0 L/min). Evaluated using an instrumented dual acrylic column process rig (A_1 = A_2 = 400.0 cm^2 = 0.0400 m^2) with inter-tank laminar capillary valve and final laminar discharge tube. Independent piezoresistive pressure level transmitters (±0.05 cm) on both vessels.

12. Conclusion

This laboratory investigation quantitatively verified the mathematical formulation and physical behavior of Interacting vs Non-Interacting Two-Tank Liquid Level Dynamics & Coupled State-Space Modeling. Through rigorous tracking across six observation points, we confirmed an empirical Interacting Coupled Process Gain K_p of 4.2500 cm/(L/min) matching the theoretical reference 4.2500 cm/(L/min) within 0.0000% error. The investigation demonstrates why this chemical engineering method serves as an indispensable cornerstone in industrial process design, unit operations scaling, and plant optimization.

In-Depth Technical Supplement & Historical Unit Operations Derivation

1. Objective & Practical Significance

This laboratory investigation determines the dynamic coupling, transfer function poles, and steady-state process gain of a two-tank interacting liquid level system ($A_1 = A_2 = 400.0\text{ cm}^2 = 0.0400\text{ m}^2$) subjected to step perturbations in Tank 1 influent flow rate $\Delta F_1 = 1.0$ to $6.0\text{ L/min}$. By monitoring piezoresistive level responses $h_1(t)$ and $h_2(t)$, we evaluate the overdamped second-order sigmoidal response, verify the linear gain model $\Delta h_2 = K_p \cdot \Delta F_1$, and determine the coupled process gain $K_p = 4.2500\text{ cm/(L/min)}$.

2. Governing Transport Equations & Phenomenological Relations

The interacting transfer function is: $G_2(s) = \frac{H_2'(s)}{F_1'(s)} = \frac{K_p}{\tau_1 \tau_2 s^2 + (\tau_1 + \tau_2 + A_1 R_2) s + 1}$. The steady-state level deviation model is:

\[ \Delta h_2 = K_p \cdot \Delta F_1 = 4.2500 \cdot \Delta F_1 \]

where $K_p = R_2 = 4.2500\text{ cm/(L/min)} = 255.0\text{ s/m}^2$.

3. Experimental Apparatus & Instrumentation Specification

The apparatus consists of two matched acrylic cylinders ($22.56\text{ cm}$ ID, $A = 400.0\text{ cm}^2$, height $80\text{ cm}$), an inter-tank laminar capillary valve ($R_1 = 2.500\text{ cm/(L/min)}$), a Tank 2 discharge resistor ($R_2 = 4.250\text{ cm/(L/min)}$), two base piezoresistive pressure level transmitters ($\pm 0.02\text{ cm}$), a servo-metered gear pump, and a 16-bit DAQ interface sampling at $10\text{ Hz}$.

4. Standard Operating Procedure & Measurement Protocols

Transmitters were hydrostatically zeroed. The inter-tank valve was opened fully to establish interacting mode. Baseline flow $F_{1,\text{initial}} = 3.0\text{ L/min}$ was established until levels $h_{1,\text{ss}} = 20.25\text{ cm}$ and $h_{2,\text{ss}} = 12.75\text{ cm}$ were steady. Step increases in feed flow $\Delta F_1 = 1.0, 2.0, 3.0, 4.0, 5.0,$ and $6.0\text{ L/min}$ were applied. Transients $h_1(t)$ and $h_2(t)$ were logged until steady state, recording final level deviations $\Delta h_2$.

5. Safety Analysis, Hazard Mitigation & Personal Protection

Both tanks were equipped with continuous overflow lines draining to the central reservoir. Electrical wiring for sensors was conduit-sealed and elevated. Protective eyewear, slip-resistant laboratory footwear, and waterproof lab coats were worn.

6. Raw Experimental Observations & Systematic Data Records

Simulated demonstration data. The following coupled level measurements were recorded at $T = 20.0^\circ\text{C}$ ($A_1 = A_2 = 400.0\text{ cm}^2$, $F_{1,\text{initial}} = 3.0\text{ L/min}$):

RunPerturbation $\Delta F_1$ (L/min)Tank 1 Deviation $\Delta h_1$ (cm)Tank 2 Deviation $\Delta h_2$ (cm)Total Feed $F_1$ (L/min)Inter-Tank Flow $F_{12}$ (L/min)Apparent Gain $K_p$
11.06.754.254.004.004.250
22.013.508.505.005.004.250
33.020.2512.756.006.004.250
44.027.0017.007.007.004.250
55.033.7521.258.008.004.250
66.040.5025.509.009.004.250

7. Quantitative Mathematical Analysis & Regression Models

Linear regression of Tank 2 level deviation $\Delta h_2$ versus feed flow perturbation $\Delta F_1$ yields:

\[ \Delta h_2 = 4.25000 \cdot \Delta F_1 + 0.00000 \]

with coefficient of determination $R^2 = 1.000000$ and zero intercept. The experimental slope $K_p = 4.25000\text{ cm/(L/min)}$ matches calibration target ($4.2500$) with $0.00\%$ deviation. Evaluating Tank 1 steady-state gain:

\[ K_{p,1} = \frac{\Delta h_1}{\Delta F_1} = R_1 + R_2 = 2.5000 + 4.2500 = 6.7500\text{ cm/(L/min)} \]

which matches the observed Tank 1 slope ($6.7500\text{ cm/(L/min)}$) exactly.

8. Rigorous Error Analysis & Uncertainty Propagation

Piezoresistive pressure transmitter precision is $\pm 0.02\text{ cm}$, and feed pump delivery error is $\pm 0.015\text{ L/min}$. Propagating into process gain yields an experimental uncertainty in $K_p$ of $\pm 0.026\text{ cm/(L/min)}$ ($0.61\%$), reporting $K_p = 4.250 \pm 0.026\text{ cm/(L/min)}$.

9. Engineering Discussion & Phenomenological Interpretation

The measured Tank 2 step response exhibited clear S-shaped inflection with zero initial derivative ($dh_2/dt|_{t=0} = 0$), confirming second-order overdamped dynamics. The effective time constants extracted from the transfer function were $ au_a = 210.8\text{ s}$ and $ au_b = 40.2\text{ s}$, yielding a total settling time of $\sim 900\text{ s}$. The identical match between $K_p$ and $R_2$ proves that inter-tank resistance $R_1$ affects only transient response speed without altering overall steady-state gain.

10. Industrial Scale-Up & Unit Operation Applications

Interacting multi-capacity dynamics govern distillation column tray hydraulics, continuous counter-current liquid-liquid extraction columns, anaerobic digester cascades, and automated boiler drum level control networks.

11. Self-Assessment Conceptual Questions & Model Answers

Q1: Why does the Tank 1 level deviation depend on both $R_1$ and $R_2$ ($\Delta h_1 = (R_1 + R_2) \Delta F_1$)?
A1: At steady state, flow through resistor 1 is $(h_1 - h_2)/R_1 = \Delta F_1$, meaning $h_1 - h_2 = R_1 \Delta F_1$; substituting $h_2 = R_2 \Delta F_1$ yields $h_1 = (R_1 + R_2) \Delta F_1$.

Q2: How can the interaction factor be minimized if a fast dynamic response is required?
A2: By decoupling the vessels: elevating Tank 1 above Tank 2 so effluent discharges freely into the open atmosphere above Tank 2, eliminating the backpressure cross-term $A_1 R_2$.

12. Concluding Summary & Verified Engineering Metric

Dynamic step response analysis of the interacting two-tank liquid level system validated coupled second-order process theory. The experimental process gain $K_p = 4.2500\text{ cm/(L/min)}$ confirmed steady-state hydraulic resistance and state-space conservation modeling.