Reciprocating Pump Pulsation Simulator

Every reciprocating pump breathes: each cylinder delivers in bursts, and the combined flow ripples around its mean. Pick the throw count (2–9), flip single/double-acting, tune the connecting-rod ratio C = L/R — the simulator draws the instantaneous flow curve q(θ) from slider-crank kinematics and reads the pulsation amplitude against the classic cylinder-count table: duplex ≈ 46%, triplex ≈ 23%, quintuplex ≈ 7%, nonuplex ≈ 3%. Odd throws always beat the neighbouring even ones. A pumpXSolver engineering tool.

q(θ) summed over cylinders mean flow q̄ qmax (above mean) qmin (below mean)

Pulsation readouts — versus the classic table

Suction side — the standing warning

Pump kinematics — pick the cylinders

Cylinders (throws)
Acting

Speed — for the pulse frequency

Speed scales the pulsation frequency, not its amplitude: f = N·n/60 (single-acting), 2N·n/60 (double-acting).

How it works — slider-crank kinematics

A crank and connecting rod drive each plunger. With crank radius \(r\), rod length \(L\), ratio \(C=L/R\) and crank angle \(\theta\), the plunger velocity is

$$v(\theta)=r\omega\left(\sin\theta+\frac{\sin 2\theta}{2C}\right)\qquad X(\theta)=r(1-\cos\theta)+L\left(1-\sqrt{1-(r/L)^2\sin^2\theta}\right)$$

The \(\sin 2\theta\) term is the finite-rod distortion: it sharpens the velocity peak and shifts it ahead of 90° — smaller C, sharper peak, bigger ripple. One cylinder's instantaneous delivery is \(q(\theta)=A\,v(\theta)\) on its discharge stroke; the simulator sums the contributions of \(N\) cylinders — phased at \(2\pi/N\) single-acting, and at \(\pi/N\) double-acting (a duplex DA carries its two cranks 90° apart):

$$q_{SA}(\theta)=\sum_{k=0}^{N-1}\max\!\bigl(0,\;v(\theta-2\pi k/N)\bigr)\qquad q_{DA}(\theta)=\sum_{k=0}^{N-1}\bigl|\,v(\theta-\pi k/N)\,\bigr|$$

The pulsation amplitude is measured against the mean flow \(\bar q\):

$$\text{above}=\frac{q_{max}-\bar q}{\bar q}\qquad \text{below}=\frac{\bar q-q_{min}}{\bar q}\qquad \text{total}=\frac{q_{max}-q_{min}}{\bar q}$$
SymbolMeaningUnits
\(r,\ L,\ C\)crank radius, rod length, ratio \(C=L/R\)m, m, —
\(\omega,\ \theta\)crank angular speed and anglerad/s, °
\(v,\ X\)plunger velocity and displacementm/s, m
\(q,\ \bar q\)instantaneous combined flow and its meanm³/s
\(N\)number of cylinders (throws), phased \(2\pi/N\)
\(A\)plunger area (cancels in the normalised curve)

1 · Cylinder count rules the ripple (C ≈ 6:1)

The classic pulsation table, reproduced by the simulator within a percent — note how the odd counts beat both even neighbours (triplex 23% sits between duplex 46% and quadruplex 33%):

Pump typePlungersAbove meanBelow meanTotalPhase angle
Duplex (double-acting)224%22%46%180°
Triplex36%17%23%120°
Quadruplex411%22%33%90°
Quintuplex52%5%7%72°
Sextuplex65%9%14%60°
Septuplex71%3%4%51.5°
Nonuplex91%2%3%40°

That is why power-pump crankshafts always take an odd number of throws — the duplex double-acting, whose two cranks stand 90° apart and whose rod-side stroke is counted without rod-area correction in the classic table, is the standing exception.

2 · The L/R ratio is the fine adjustment (triplex)

C = L/RAbove meanBelow meanTotal
4 : 18.2%20.0%28.2%
5 : 17.6%17.6%25.2%
6 : 16.9%16.1%23.0%
7 : 16.4%15.2%21.6%

A shorter rod (small C) steepens the velocity peak and worsens the ripple; a longer rod smooths the flow but stretches the frame. The usual industrial compromise is 4:1–6:1 — C also enters volumetric efficiency, dampener sizing, the liquid-separation speed and the acceleration head, so it is never chosen on pulsation alone.

3 · When the ripple is too big — dampeners and acoustic filters

Above 25% total pulsation, fit a discharge dampener (🔴); 10–25% is workable with sturdy piping (🟡); ≤ 10% is smooth (🟢). A bladder dampener precharged to 50–66% of system pressure can attenuate pulsations below ~50 Hz by up to 90%, and a well-designed suction stabiliser makes the suction line behave like only 5–15 pipe diameters of straight pipe. For stubborn cases an acoustic low-pass filter — a volume–choke–volume Helmholtz device — works when its natural frequency sits at or below half the plunger frequency:

$$f=\frac{a\,d}{\pi\sqrt{2}\;D\,L}\;\le\;\frac{f_{\text{plunger}}}{2}$$

a = speed of sound in the liquid, d = choke diameter, D = chamber diameter, L = chamber and choke length. Pipe runs themselves have organ-pipe modes, \(f_N=Na/2L\) (open–open or closed–closed) and \(f_N=Na/4L\) (closed–open): a pulse frequency landing on one of them can multiply vibration twenty-fold — check before blaming the pump.

4 · The suction side pays twice

The kinematics that ripple the delivery also swing the suction pressure. Acceleration head — the pressure spent accelerating the liquid column in the suction line at every stroke — can be the largest single NPSHA component, in bad layouts up to ten times the sum of all other suction losses. Keep the effective NPSHA comfortably above NPSHR with the pulsation margin included, and remember the absolute pressure at the plunger must stay above vapour pressure through the whole cycle, not just on average.

Worked example (anchor)

Triplex single-acting, C = 6:1, 300 rpm. The simulator reproduces the table row:

QuantitySimulatorClassic table
Above mean≈ +6%6%
Below mean≈ −17%17%
Total pulsation≈ 23%23%
Discharge pulses / rev33 (fundamental)
Pulse frequency3 × 300 / 60 = 15 Hz
Verdict🟡 10–25% band — suction stabiliser pays on long suction lines

Slide to quintuplex and the same readouts fall to ≈ 2% / 5% / 7% (🟢); slide to duplex double-acting and they climb to ≈ 24% / 22% / 46% (🔴 — dampener required).

Engineering criteria applied