When Operating Chemical Engineering Pilot Plants, How Is Draining Time Calculated For Unsteady-State Flow?

Draining time for a pilot-plant tank is not a steady-state problem. The liquid level falls, the discharge velocity slows, and these changes are coupled. To calculate the time required, you combine an instantaneous material balance with the Bernoulli equation to relate velocity to height, then integrate the resulting differential equation over the changing liquid level. This gives you the exact draining time—or the level at any moment—without needing a constant-flow assumption.

The core principle is that the tank’s volume decrease rate must equal the instantaneous outflow through the pipe. Bernoulli links that outflow velocity directly to the current liquid height, turning a dynamic mass balance into an integrable differential equation. Once integrated, this model accurately predicts the draining time and can be verified against actual pilot-plant measurements.

How the Draining Time Equation is Built

The Dynamic Material Balance

Over an infinitesimal time step $d\theta$, the liquid height drops by $dh$. The volume lost from the tank is $-A,dh$, where $A$ is the tank’s cross-sectional area (assuming a vertical-walled vessel).

This lost volume must exit through the discharge pipe. With pipe cross-sectional area $a$ and instantaneous velocity $u$, the outflow volume in the same interval is $a,u,d\theta$. Equating the two gives the material balance relationship:

$$-A,dh = a,u,d\theta$$

This equation is not immediately usable because $u$ is a function of the changing level $h$.

Linking Velocity to Height with Bernoulli

Apply the mechanical energy balance (Bernoulli) between the liquid surface and the pipe exit. For a simple bottom discharge with negligible inlet velocity and a common exit pressure, you get:

$$u = C_d \sqrt{2g(h - h_f)}$$

Here, $h$ is the instantaneous height above the pipe exit, $g$ is gravity, and $h_f$ accounts for friction and minor losses in the pipe. The discharge coefficient $C_d$ bundles all non-idealities—typically determined empirically or from standard loss correlations.

If the pipe friction is small relative to the static head (common in short pilot-plant drains), $h_f$ can be neglected or absorbed into $C_d$.

The Integration Step

Substituting the velocity expression into the material balance separates the variables:

$$-A,dh = a,C_d\sqrt{2g(h - h_f)},d\theta$$

Rearrange to isolate $h$ and $\theta$:

$$d\theta = -\frac{A}{a,C_d\sqrt{2g}},\frac{dh}{\sqrt{h - h_f}}$$

Integrate the left side from $0$ to the total draining time $\theta_f$, and the right side from the initial height $h_0$ to the final height $h_f$ (usually $0$ or the pipe exit level adjusted for $h_f$). The result for a constant-area tank and negligible backpressure:

$$\theta_f = \frac{2A}{a,C_d\sqrt{2g}}\left(\sqrt{h_0 - h_f} - \sqrt{h_f - h_f}\right)$$

For a tank draining completely to the exit level, $h_f = 0$, this simplifies to $\theta_f = \frac{2A\sqrt{h_0}}{a,C_d\sqrt{2g}}$.

This integrated equation directly gives the total draining time—the number you need for scheduling a pilot-plant batch or designing a drain-down procedure.

Applying the Method in a Pilot Plant Setting

Accounting for Friction and Minor Losses

Pilot-plant drain lines often include elbows, valves, and a length of pipe that induce significant frictional losses. These directly reduce the effective driving head. You can either:

Include a friction term $h_f$ as a function of velocity (making the integration slightly more complex), or
Absorb all losses into an experimentally determined $C_d$ by measuring the actual flow rate at a known head.

For accuracy, measure $C_d$ from a calibration run; this captures both pipe friction and entrance/exit effects without requiring detailed hydraulic modeling.

Experimental Verification

You can validate your calculation by tracking the water level over time on the tank’s graduated scale. Compare the theoretical curve (from the integrated equation) with the actual level drop. Discrepancies often point to incorrect assumptions about $A$ (e.g., a domed bottom), a varying $C_d$, or vortex formation that draws air and reduces effective outflow.

Understanding the Trade-offs and Practical Pitfalls

This unsteady-state method gives precision, but it carries assumptions that can mislead if ignored.

Constant cross-sectional area: If the tank has a dished bottom or internal structures, $A$ varies with height. The integration must account for that geometry; otherwise, time estimates will be off.
Friction coefficient variability: A fixed $C_d$ works for a smooth drain, but if the pipe’s roughness changes (e.g., due to scaling) or the flow regime shifts, the coefficient drifts.
Vortex and air entrainment: When the liquid level drops near the drain opening, a vortex can form, allowing air to enter the pipe and drastically reducing flow rate. The Bernoulli-based model does not predict this—real draining often takes longer.
Backpressure and siphon effects: If the discharge pipe empties into a closed vessel or a siphon forms, the driving head changes. The simple open-to-atmosphere assumption must be adjusted.

These trade-offs mean the calculated time should be treated as a best-case hydraulic answer. For safety-critical operations, always add a margin of 15–30 % unless your $C_d$ was determined specifically under similar conditions.

Making the Right Choice for Your Goal

Your draining time calculation strategy depends on your exact need and available data.

If your primary focus is operator training or educational demonstration: Use the full differential equation and integrate it. Compare the theoretical level-time curve with real-time readings from the graduated vessel to build intuition about unsteady-state behavior.
If your primary focus is a rapid field estimate for a standard tank: Use the simplified integrated formula with a known discharge coefficient (e.g., 0.6 for a sharp-edged orifice) and a constant $A$. Accept an uncertainty of ±20 %.
If your primary focus is designing a reliable drain-down sequence: Start with the integrated Bernoulli model, then incorporate a measured $C_d$ from a small-scale test. Add a safety factor to account for vortex effects, air locking, or unexpected pipe blockage during real operation.
If your primary focus is optimizing a piping layout to minimize drain time: Run the differential model iteratively while varying pipe diameter $a$, length, or fittings to see the impact on $C_d$ and total time—this pinpoints the most effective upgrade.

By anchoring your calculation in the dynamic material balance and a realistic discharge coefficient, you turn a seemingly complex unsteady-state problem into a manageable, defensible engineering result.

Summary Table:

Parameter	Symbol	Description & Role
Tank Area	$A$	Cross-sectional area of the vessel
Pipe Area	$a$	Cross-sectional area of the discharge pipe
Discharge Coefficient	$C_d$	Empirical value accounting for friction & losses
Liquid Height	$h$	Instantaneous head driving the flow
Total Draining Time	$\theta_f$	Calculated duration for level to drop from $h_0$ to $h_f$

Enhance Your Engineering Programs with LABPARK

Putting fluid dynamics theory into practice requires reliable, real-world systems. LABPARK provides state-of-the-art Educational and Vocational Unit Operations Pilot Plants in chemical engineering, bioprocess & biotech, and environmental & water treatment. Specially designed for universities, research institutes, and enterprises, our systems offer hands-on validation of unsteady-state flow, heat transfer, and mass balance principles.

Ready to elevate your training and research labs? Contact LABPARK today to find the perfect pilot plant for your institution!