Knowledge Chemical Engineering Education How can the Strouhal number and fluid time scales be used to evaluate pseudosteady-state conditions in pulsatile or unsteady pilot plant operations?
Author avatar

Tech Team · LABPARK

Updated 2 months ago

How can the Strouhal number and fluid time scales be used to evaluate pseudosteady-state conditions in pulsatile or unsteady pilot plant operations?


Use the Strouhal number as a time-scale test: pseudosteady behavior is reasonable when the imposed pulsation is slow compared with the fluid’s relevant response time. For inertial pilot-plant flows, compare the forcing period with the convective time ( \tau_c=L/U ); for creeping flows, compare it with the viscous diffusion time ( \tau_v=L^2/\nu ). The flow may be treated as pseudosteady when these ratios are sufficiently small, but the appropriate criterion depends on the dominant transport regime.

Pseudosteady modeling is justified when the fluid adjusts much faster than the operating conditions change. Evaluate that adjustment using the convective time scale for inertial flows and the viscous time scale for low-Reynolds-number flows, while checking local peaks and harmonics rather than relying only on average operating conditions.

Identify the Relevant Time Scales

The imposed time scale

For a pulsatile pump, oscillating agitator, or cyclic valve, define the imposed time scale as

[ \tau_i \approx \frac{1}{f} ]

where (f) is the dominant forcing frequency. If the waveform is nonsinusoidal, higher harmonics have shorter time scales and may control the validity of the approximation.

For a periodic operation with period (T), using (\tau_i=T) gives the same physical comparison, with only a constant-factor difference from using (1/f).

The convective time scale

The convective time scale is

[ \tau_c=\frac{L}{U} ]

where (L) is the characteristic flow length and (U) is the characteristic velocity.

It represents the time required for fluid momentum or concentration information to travel through the relevant equipment region, such as a pipe diameter, mixing zone, packed section, or column length.

The viscous response time

For viscous or creeping flow, the relevant response time is the viscous diffusion time:

[ \tau_v=\frac{L^2}{\nu} ]

where (\nu) is the kinematic viscosity.

This is the time required for viscous effects to diffuse across the characteristic length. In a small gap or narrow channel, this time may be short; in a large vessel or highly viscous system, it may be substantial.

Define the Strouhal Number Consistently

The time-scale convention

For evaluating pseudosteady conditions, define the time-scale form as

[ Sr=\frac{\tau_i}{\tau_c} =\frac{U}{fL} ]

Under this convention, a large (Sr) means that the imposed change is slow relative to convection.

Many fluid-dynamics references instead use the conventional Strouhal number

[ St=\frac{fL}{U} =\frac{\tau_c}{\tau_i} ]

Thus,

[ St=\frac{1}{Sr} ]

Both definitions describe the same physics. The critical requirement is to state which convention is being used before interpreting a value as “large” or “small.”

The characteristic length and velocity

The result depends directly on the selected (L) and (U). A pipe diameter may be appropriate for internal pipe flow, while an impeller diameter, blade-passage length, or mixing-zone scale may be more appropriate for an agitated vessel.

For strongly pulsatile systems, use representative peak or cycle-specific velocities in addition to mean velocity. A mean-based value can conceal short intervals in which inertia becomes important.

Evaluate Inertial Flow

Compare forcing and convection

For flows with significant inertia, pseudosteady behavior requires

[ \tau_i \gg \tau_c ]

or, using the time-scale Strouhal convention,

[ \frac{1}{Sr}\ll 1 ]

Equivalently, using the conventional definition,

[ St\ll 1 ]

The fluid then traverses or adjusts across the relevant flow region much faster than the boundary conditions change. At each instant, the velocity field can often be approximated by the steady solution corresponding to the instantaneous flow rate.

Interpret the result physically

If (St) is small, a reciprocating pump’s changing flow rate is slow from the fluid’s perspective. The pressure and velocity fields have time to track the changing boundary condition with relatively little phase lag.

If (St) approaches unity or becomes larger, local acceleration and delayed momentum transport can no longer be ignored. A fully transient model, or at least an unsteady correction, becomes necessary.

Account for Reynolds number

The Reynolds number is

[ Re=\frac{UL}{\nu} ]

It indicates whether inertial or viscous effects dominate, but it does not by itself determine whether the flow is pseudosteady. A high-Reynolds-number flow can still be pseudosteady if its forcing is sufficiently slow, and a low-Reynolds-number flow can still be unsteady if it is forced rapidly.

The two dimensionless groups answer different questions:

  • (Re) compares inertia with viscosity.
  • (St) or (Sr) compares the imposed time scale with the convective time scale.

Evaluate Creeping or Viscous Flow

Compare forcing and viscous diffusion

For (Re\ll1), the appropriate criterion is

[ \tau_i\gg\tau_v ]

Using the time-scale Strouhal definition,

[ \frac{\tau_v}{\tau_i}

\frac{Re}{Sr} \ll1 ]

This means viscous momentum diffusion is fast compared with the imposed change.

Under these conditions, the transient velocity derivative in the unsteady Stokes equations is small relative to the viscous term. The velocity field responds almost instantaneously to the changing pressure gradient or boundary motion.

Use the local scale

The viscous criterion is sensitive to (L^2). Applying a vessel-scale length to a narrow-gap flow can substantially overestimate the response time, while using a gap width for a full-vessel circulation process can underestimate it.

For pilot plants, evaluate separate regions when necessary:

  • Narrow tubing or membrane channels.
  • Agitator gaps and boundary layers.
  • Porous media or packed beds.
  • Large vessel-scale circulation zones.

The slowest relevant region often controls whether a complete pseudosteady approximation is acceptable.

Apply the Test to Pulsatile Operations

Reciprocating pumps

For a reciprocating pump, calculate the dominant frequency and velocity variation over one cycle. Evaluate (St=fL/U) using both mean and peak velocity where the waveform is strongly asymmetric.

A small (St) supports treating the flow as a sequence of steady states. A larger value indicates that pump acceleration, fluid inertia, line compliance, and pressure-wave propagation may need to be modeled explicitly.

Oscillating agitators

For an oscillating agitator, the characteristic frequency may be the oscillation frequency or a blade-passage frequency. The characteristic length may be the agitator diameter, blade scale, or local mixing-zone dimension, depending on the question being studied.

A vessel may be pseudosteady in its bulk circulation while remaining unsteady near the impeller, baffles, or free surface. A single global Strouhal number should therefore not be assumed to represent every region.

Multi-frequency waveforms

A periodic waveform may contain a fundamental frequency and several harmonics. Since

[ St_n=\frac{n f L}{U} ]

the (n)-th harmonic has a Strouhal number (n) times larger than that of the fundamental.

The highest dynamically significant harmonic should be checked. Otherwise, a waveform may appear pseudosteady based on its average period while its sharper transitions remain genuinely unsteady.

Connect the Criterion to Measurements

Compare predicted and observed phase lag

A pseudosteady model predicts that pressure, flow, and concentration respond nearly in phase with the imposed operating signal.

Measurable phase lag, overshoot, or attenuation indicates that the relevant response time is not negligible. These observations are especially useful for validating the selected length scale and velocity scale.

Check cycle-to-cycle repeatability

A repeatable waveform does not automatically mean the flow is pseudosteady. A system can have a stable periodic solution that remains strongly unsteady.

Use repeatability to establish periodic operation, then use (St), (Re/Sr), pressure data, and local velocity measurements to determine whether each phase can be approximated by a steady state.

Separate fluid dynamics from instrumentation effects

Flowmeters and pressure sensors can introduce their own response delays. The measured phase lag may therefore reflect the instrument, the fluid, or both.

Instrument response should be characterized before using phase comparisons to reject or accept a pseudosteady model.

Understand the Trade-offs

A single threshold is not universal

Terms such as “much greater than” and “much less than” are engineering criteria, not universal constants. The acceptable limit depends on the required accuracy, the quantity being predicted, and the sensitivity of the equipment to transient effects.

A rough pseudosteady velocity estimate may tolerate a larger unsteady contribution than a pressure-drop, residence-time, mixing, or structural-load calculation.

Vortex-shedding Strouhal numbers are different in purpose

The supplementary vortex-shedding relation

[ St=\frac{fd}{V} ]

is useful for characterizing periodic separation around an obstruction and for interpreting vortex-shedding flowmeters. Values near (0.21) in a broad subcritical regime can help identify shedding behavior and resonance risks.

However, that Strouhal number should not be confused with the forcing-to-convection criterion for pseudosteady modeling. In a pilot plant, vortex shedding can generate an additional frequency that must be checked, but its value does not by itself prove that the overall process is pseudosteady.

Structural resonance can remain important

A flow may be sufficiently slow for a pseudosteady fluid model while still exciting piping or column internals near a structural natural frequency.

Hydrodynamic simplification and structural safety are separate checks. Periodic forcing and vortex-shedding frequencies should be compared with equipment resonance frequencies when vibration is a concern.

Other transient mechanisms may dominate

The time-scale test does not replace checks for compressibility, fluid-structure interaction, pump-line compliance, cavitation, or species and heat-transfer transients.

For liquid systems with long lines or significant pressure-wave effects, acoustic or hydraulic-wave time scales may be more relevant than (L/U) alone.

Making the Right Choice for Your Goal

Start by calculating (Re), (St=fL/U), and the relevant response-time ratio using local and cycle-specific scales.

  • If your primary focus is bulk velocity or pressure prediction: Use a pseudosteady model when (St\ll1) for inertial flow, or (Re/Sr\ll1) for creeping flow, and validate it against phase-lag measurements.
  • If your primary focus is mixing, residence time, or mass transfer: Evaluate the slowest relevant transport time and the highest significant forcing harmonic rather than relying on vessel-average values.
  • If your primary focus is vibration or equipment integrity: Calculate forcing and vortex-shedding frequencies separately and compare them with structural natural frequencies, even when the fluid model appears pseudosteady.
  • If your primary focus is high-fidelity transient prediction: Retain the unsteady terms whenever (St) is not small, (Re/Sr) is not small in viscous flow, or measurements show meaningful phase lag and waveform distortion.

A flow is pseudosteady only when its relevant fluid response time is demonstrably shorter than the time scale of the imposed operation.

Summary Table:

Evaluation factor Key equation or check Practical interpretation
Imposed time scale tau_i = 1/f or T Use the dominant frequency and check significant harmonics.
Convective response tau_c = L/U Relevant for flows with significant inertia.
Viscous response tau_v = L^2/nu Relevant for creeping or viscosity-dominated flow.
Conventional Strouhal number St = fL/U St << 1 supports pseudosteady inertial-flow modeling.
Time-scale Strouhal number Sr = U/(fL) = 1/St Large Sr indicates slow forcing relative to convection.
Creeping-flow criterion Re/Sr << 1 Viscous diffusion is fast relative to the imposed change.
Local validation Check phase lag, attenuation, and overshoot Measured distortion may require transient modeling.
Additional risks Check wave propagation, harmonics, cavitation, and resonance Fluid pseudosteady behavior does not guarantee equipment or structural safety.

Need reliable pilot-plant equipment for dynamic flow, mixing, and transport studies? Contact LABPARK to discuss your application. LABPARK provides Educational and Vocational Unit Operations Pilot Plants for chemical engineering, bioprocess and biotech, and environmental and water treatment applications. Our systems help universities, research institutes, and enterprises evaluate pulsatile operations, validate pseudosteady assumptions, investigate transient effects, and train operators with practical, application-focused equipment. Share your process conditions, target scale, and measurement requirements so our team can help identify an appropriate pilot-plant configuration.


Leave Your Message