The Kármán Spacing Ratio of the Vortex Street
This appendix proves Equation (35.54) of Phenomenon 35.61 in Fluid Dynamics: among all staggered double rows of point vortices, exactly one ratio of row separation to streamwise spacing is not destroyed by its own induced motion, namely
[Karman:1911]. What the calculation settles and what it does not is stated in the chapter at Remark 35.62, and the reader should have that remark in view throughout: nothing here bears on the Strouhal number Equation (35.53), and the stability obtained is neutral, not asymptotic. The closing Remark A61.10 returns to the point with the calculation in hand.
Throughout, the fluid is ideal, unbounded and two-dimensional; the vortices are point vortices in the sense of Section 35.7.1, each of circulation \(\pm\Gamma\) in \(\mathrm{m}^{2}/\mathrm{s}\), moving with the local velocity induced by all the others. The complex variable is \(z=x+\ii y\), the complex potential \(w\) is that of Equation (35.57), and a vortex of circulation \(\Gamma\) counterclockwise at \(z_{0}\) has \(w=-\left(\ii\Gamma/2\pi\right)\log\left(z-z_{0}\right)\), so that
the second equation being the statement that each vortex is carried by the flow of the others.
One summation formula does all the work
Every lattice sum below is a special case of a single identity, which is proved from the Fourier series of an elementary function and needs nothing else.
Let \(\zeta\in\C\setminus\Z\) and \(0<\varphi<2\pi\). Then
and consequently, differentiating with respect to \(\zeta\),
whose limit as \(\varphi\rightarrow0^{+}\) is \(\Sigma\left(0,\zeta\right)=\pi^{2}/\sin^{2}\pi\zeta\). Rests on Theorem 12.6 and Equation (12.3).
Derives Lemma A61.1. Let \(F(\varphi):=\left(\pi/\sin\pi\zeta\right) \ee^{\ii\left(\pi-\varphi\right)\zeta}\) on \((0,2\pi)\) and compute its Fourier coefficients:
the elementary integral of \(\ee^{-\ii\left(\zeta+n\right)\varphi}\). Since \(n\) is an integer, \(\ee^{-2\pi\ii\left(\zeta+n\right)} =\ee^{-2\pi\ii\zeta}\), and \(1-\ee^{-2\pi\ii\zeta} =\ee^{-\ii\pi\zeta}\left(\ee^{\ii\pi\zeta}-\ee^{-\ii\pi\zeta}\right) =2\ii\,\ee^{-\ii\pi\zeta}\sin\pi\zeta\) by Equation (12.3). The two exponentials and the two sines cancel, leaving \(c_{n}=1/(n+\zeta)\), which is Equation (A61.3).
The series Equation (A61.4) converges absolutely and uniformly in \(\zeta\) on compact subsets of \(\C\setminus\Z\), so it may be obtained from Equation (A61.3) by differentiating term by term, using \(\dd\left(n+\zeta\right)^{-1}/\dd\zeta=-\left(n+\zeta\right)^{-2}\) and the product rule on the right-hand side. For the limit, set \(\varphi=0\) in Equation (A61.4): the bracket becomes \(\pi\left(\cos\pi\zeta-\ii\sin\pi\zeta\right)/\sin^{2}\pi\zeta =\pi\,\ee^{-\ii\pi\zeta}/\sin^{2}\pi\zeta\), which cancels the prefactor \(\ee^{\ii\pi\zeta}\) exactly.
∎The complex potential of vortices of circulation \(\Gamma\) placed at \(z=na\) for every \(n\in\Z\) is
and it is the only such potential whose velocity is bounded as \(\abs{\operatorname{Im}z}\rightarrow\infty\). The induced velocity tends to \(\mp\Gamma/2a\) in the \(x\) direction as \(\operatorname{Im}z\rightarrow\pm\infty\). Rests on Equation (A61.2) and Theorem 12.18.
Derives Corollary A61.2. The function \(\cot\left(\pi z/a\right)\) is meromorphic with simple poles exactly at \(z=na\), where \(\sin\left(\pi z/a\right)\) has its simple zeros, and near \(z=na\) its principal part is \(a/\pi\left(z-na\right)\). Hence \(\dd w/\dd z\) in Equation (A61.5) behaves near \(z=na\) as \(-\left(\ii\Gamma/2a\right)\cdot a/\pi\left(z-na\right) =-\ii\Gamma/2\pi\left(z-na\right)\), which by Equation (A61.2) is precisely a vortex of circulation \(\Gamma\) there, and it is holomorphic elsewhere. As \(\operatorname{Im}z\rightarrow\pm\infty\), \(\cot\left(\pi z/a\right)\rightarrow\mp\ii\), so \(\dd w/\dd z\rightarrow\mp\Gamma/2a\); since \(\dd w/\dd z=v_{x}-\ii v_{y}\), the row drives the fluid backwards above it and forwards below it, at the stated speed. Any other candidate differs from Equation (A61.5) by a function that is entire (the singularities coincide) and bounded (both velocities are), hence constant by Liouville's theorem (Theorem 12.18).
∎The street and its rigid translation
Fix \(a>0\) and \(h>0\), both in \(\mathrm{m}\). The staggered vortex street consists of the rows
for all \(n\in\Z\). Write \(s:=\pi h/a\), \(c:=\cosh s\) and \(t:=\tanh s\) throughout. Rests on Corollary A61.2.
The configuration Equation (A61.6) is an exact solution of Equation (A61.2): every vortex moves with the same velocity
parallel to the rows, so the street is carried downstream without change of shape. Rests on Definition A61.3 and Corollary A61.2.
Derives Proposition A61.4. Take the vortex \(A_{0}=0\). Its own row induces nothing there: the contributions of \(A_{n}\) and \(A_{-n}\) in Equation (A61.2) are equal and opposite. Row \(B\), of circulation \(-\Gamma\) and origin \(a/2-\ii h\), contributes by Equation (A61.5)
using \(\cot\left(\theta-\pi/2\right)=-\tan\theta\) and \(\tan\left(\ii s\right)=\ii\tanh s\). The result is real and positive, so \(v_{x}=\Gamma t/2a\) and \(v_{y}=0\). Repeating at \(B_{0}=a/2-\ii h\) under row \(A\) gives \(-\left(\ii\Gamma/2a\right)\cot\left(\pi/2-\ii s\right) =-\left(\ii\Gamma/2a\right)\tan\left(\ii s\right) =\Gamma t/2a\), the same velocity; and again the vortex's own row contributes nothing. Every vortex therefore moves with Equation (A61.7).
∎Equation (A61.7) is worth reading before the stability calculation, because it is directly visible in the photographs: the street moves more slowly than the stream that made it, since \(U_{\mathrm{s}}\) is the velocity relative to the fluid at infinity and it is directed against the free stream in the frame of the cylinder. In the frame translating with the street the base state is an equilibrium, and since the base velocity is a constant it drops out of the linearized equations altogether; the analysis below may be read in either frame.
Linearization
Displace every vortex,
with \(\abs{\alpha_{n}},\abs{\beta_{n}}\ll a\). Expanding Equation (A61.2) to first order with \(\left(D+\delta\right)^{-1}=D^{-1}-\delta D^{-2}+O(\delta^{2})\), and using Proposition A61.4 to cancel the zeroth-order terms against the rigid motion, gives for each \(k\)
the signs following from the circulations \(\pm\Gamma\) carried by the inducing vortices.
Equations (A61.9) and (A61.10) are linear over \(\R\) but not over \(\C\): the left-hand sides carry the complex conjugates of the unknowns while the right-hand sides carry the unknowns themselves. The single-mode ansatz \(\alpha_{n}=\alpha\,\ee^{\ii n\varphi}\) is therefore inconsistent, since \(\bar{\alpha}_{n}=\bar{\alpha}\,\ee^{-\ii n\varphi}\) belongs to the mode \(-\varphi\): matching the \(k\) dependence of the two sides would force \(\varphi\equiv0\). The wavenumbers \(\varphi\) and \(-\varphi\) close on each other and on nothing else, so the smallest closed system is four-dimensional. Overlooking this is the one way the calculation can be got wrong while looking right, and it changes the answer: the two-dimensional system obtained by ignoring it makes the street unstable at every spacing, Equation (A61.1) included.
Accordingly set, for a fixed \(\varphi\in(0,2\pi)\),
Three lattice sums appear. With \(m=n-k\) and \(\zeta=\tfrac{1}{2}-\ii h/a\), so that \(\pi\zeta=\pi/2-\ii s\) and therefore \(\sin\pi\zeta=\cosh s=c\) and \(\cos\pi\zeta=\ii\sinh s\):
where the closed forms come from Lemma A61.1 and
For Equation (A61.12) use \(\sum_{m\ge1}\cos m\varphi/m^{2} =\pi^{2}/6-\pi\varphi/2+\varphi^{2}/4\) on \([0,2\pi]\); for Equation (A61.14) substitute \(\sin\pi\zeta=c\), \(\cos\pi\zeta=\ii\sinh s\) and \(\ee^{\ii\psi\zeta}=\ee^{\ii\psi/2}\ee^{\psi s/\pi}\) in Equation (A61.4). The companion sum with \(\varphi\) replaced by \(-\varphi\), i.e. \(\psi\) by \(-\psi\), is
Two further sums reduce to these: replacing \(n\) by \(-n\) shows that \(\sum_{n}\left(\left(n-\tfrac{1}{2}\right)a+\ii h\right)^{-2}=S_{0}\) and \(\sum_{n}\ee^{\ii n\varphi} \left(\left(n-\tfrac{1}{2}\right)a+\ii h\right)^{-2}=\tilde{S}_{1}\).
Substituting Equation (A61.11) into Equations (A61.9) and (A61.10) and matching the coefficients of \(\ee^{\ii k\varphi}\) and \(\ee^{-\ii k\varphi}\) separately — noting \(P(-\varphi)=P(\varphi)\) and \(S_{1}(-\varphi)=\tilde{S}_{1}(\varphi)\) — yields, with \(K:=\Gamma/2\pi\) and
the closed four-dimensional system
in the variables \(p:=\alpha\), \(q:=\beta\), \(r:=\bar{\gamma}\), \(u:=\bar{\delta}\), with
The growth rates
Solutions of Equation (A61.18) proportional to \(\ee^{\sigma t}\) have \(\sigma^{2}\) an eigenvalue of \(\mathsf{B}\mathsf{C}\). Three real quantities suffice to write it down:
the last being real because the two phases in Equations (A61.14) and (A61.16) cancel and \(\left(\ii\pi\right)^{2}=-\pi^{2}\); note \(XY=\Pi^{2}\).
The eigenvalues \(\lambda=\sigma^{2}\) of \(\mathsf{B}\mathsf{C}\) satisfy
Derives Lemma A61.6. Multiplying the matrices of Equation (A61.19),
whose trace is \(-K^{2}\left(-2Q^{2}+Y+X\right)\), i.e. the coefficient displayed. Its determinant is the product of the two determinants. Now \(\overline{S_{1}}=\ee^{\ii\varphi}S_{1}\) — replace \(n\) by \(-n-1\) in the series defining \(\overline{S_{1}}\) and compare — and likewise \(\overline{\tilde{S}_{1}}=\ee^{-\ii\varphi}\tilde{S}_{1}\), so the determinant of \(\mathsf{B}\) is \(-K^{2}\left(-Q^{2} +\overline{\tilde{S}_{1}}\,\overline{S_{1}}\right) =K^{2}\left(Q^{2}-\Pi\right)\), and that of \(\mathsf{C}\) is the same. The characteristic polynomial of a \(2\times2\) matrix is \(\lambda^{2}-\left(\text{trace}\right)\lambda +\left(\text{determinant}\right)\), which is Equation (A61.21).
∎The mode \(\varphi\) neither grows nor decays — all four roots \(\sigma\) are purely imaginary — if and only if
and it grows otherwise. Rests on Lemma A61.6.
Derives Proposition A61.7. Since \(\sigma^{2}=\lambda\), the four exponents are purely imaginary exactly when both roots of Equation (A61.21) are real and non-positive; if either root is complex, or real and positive, some \(\sigma\) has positive real part. Write \(\mathcal{T}=2Q^{2}-X-Y\) and \(\mathcal{D}=\left(Q^{2}-\Pi\right)^{2}\ge0\) for the two coefficients divided by \(K^{2}\) and \(K^{4}\). The discriminant factorizes,
because \(\mathcal{T}\pm2\left(Q^{2}-\Pi\right)\) are the two brackets. Put \(\pi_{+}:=\sqrt{X}\), \(\pi_{-}:=\sqrt{Y}\), so \(\Pi=\pm\pi_{+}\pi_{-}\) by \(XY=\Pi^{2}\), the sign being that of \(-MM'\). If \(MM'>0\) then \(\Pi=-\pi_{+}\pi_{-}\) and the factorization reads \(-\left(\pi_{+}+\pi_{-}\right)^{2} \left[4Q^{2}-\left(\pi_{+}-\pi_{-}\right)^{2}\right]\), non-negative exactly when \(2\abs{Q}\le\abs{\pi_{+}-\pi_{-}}\); if \(MM'<0\) then \(\Pi=+\pi_{+}\pi_{-}\) and it reads \(-\left(\pi_{+}-\pi_{-}\right)^{2} \left[4Q^{2}-\left(\pi_{+}+\pi_{-}\right)^{2}\right]\), non-negative exactly when \(2\abs{Q}\le\pi_{+}+\pi_{-}\). Both cases are covered by \(2\abs{Q}\le\abs{EM-M'/E}\pi/a^{2}\), since \(EM\) and \(M'/E\) have the same sign as \(M\) and \(M'\) respectively; and the closed form on the right of Equation (A61.22) follows from Equation (A61.15),
Finally the second requirement, \(\mathcal{T}\le0\), is implied by the first: in either case \(4Q^{2}\le\left(\pi_{+}\pm\pi_{-}\right)^{2} \le2\left(X+Y\right)\), so \(2Q^{2}\le X+Y\).
∎Only one spacing survives
The staggered street of Definition A61.3 is neutrally stable to every infinitesimal disturbance Equation (A61.11) if and only if
which is Equation (35.54). For every other ratio some mode grows exponentially. Rests on Proposition A61.7 and Definition A61.3.
Derives Theorem A61.8. Necessity. Take \(\varphi=\pi\), that is \(\psi=0\): neighbouring vortices of a row are then displaced in opposite senses, the disturbance of longest wavelength that is not a rigid translation. Then \(u=0\), \(E=1\) and \(M=M'=\pi t/c\), so the right-hand side of Equation (A61.22) vanishes identically and the criterion demands \(Q=0\). By Equation (A61.17) with \(\psi=0\) this is \(\pi^{2}/2=\pi^{2}/c^{2}\), i.e.\ Equation (A61.23). Solving, \(\sinh s=\sqrt{c^{2}-1}=1\) and \(s=\log\left(c+\sinh s\right)=\log\left(1+\sqrt{2}\right)\), whence Equation (A61.1).
Sufficiency. Let \(c=\sqrt{2}\), so \(t=1/\sqrt{2}\), \(\sinh s=1\) and \(s=\log\left(1+\sqrt{2}\right)\). Then \(Q=-\psi^{2}/2a^{2}\) by Equation (A61.17), and Equation (A61.22) becomes, after multiplying by \(a^{2}\) and using \(c=\sqrt{2}\),
Both sides are even in \(\psi\), so take \(0\le\psi\le\pi\), where the bracket is negative and Equation (A61.24) is the assertion \(G(u)\ge0\) for \(0\le u\le s\), with \(\psi=\pi u/s\) and
obtained by dividing Equation (A61.24) by \(\pi^{2}\) and substituting \(\sqrt{2}=\cosh s\), \(1=\sinh s\). Now \(G(0)=0\) and \(G(s)=\cosh^{2}s-\sinh s-1=2-1-1=0\); moreover
Since \(2\cosh s/s>1\), every \(u\)-dependent term of \(G''\) is strictly increasing on \([0,s]\), so \(G''\) is; and \(G''(0)=-2/s^{2}<0\) while \(G''(s)=2\sqrt{2}/s+1-2/s^{2}>0\) numerically (\(s=0.88137\) gives \(3.209+1-2.575\)). Hence \(G''\) changes sign exactly once, so \(G'\) decreases and then increases. As \(G'(0)=\cosh s/s-1>0\) and \(G'(s)=\cosh^{2}s/s+\cosh s\sinh s-\cosh s-2/s=0\) — the two terms \(\cosh^{2}s/s=2/s\) and \(-2/s\) cancel, as do \(\cosh s\sinh s=\sqrt{2}\) and \(-\cosh s=-\sqrt{2}\) — the function \(G'\) is positive on \((0,u^{*})\) and negative on \((u^{*},s)\) for a single \(u^{*}\). Therefore \(G\) rises from \(G(0)=0\) and falls back to \(G(s)=0\) without crossing zero in between, so \(G\ge0\) on \([0,s]\), which is Equation (A61.24). Equality holds only at \(\psi=0\) and \(\psi=\pi\), the mode of Equation (A61.23) and the rigid translation.
Instability elsewhere. If \(c\neq\sqrt{2}\) then \(Q\neq0\) at \(\psi=0\), while \(X=Y=-\Pi=\pi^{4}t^{2}/a^{4}c^{2}\) there, so the discriminant computed in Proposition A61.7 equals \(\left[4\Pi\right]\left[4Q^{2}\right]=16\Pi Q^{2}<0\): the roots \(\lambda\) are a complex conjugate pair, their square roots are not purely imaginary, and the mode \(\varphi=\pi\) grows.
∎If the second row is placed at \(B_{n}=na-\ii h\), directly beneath the first, the mode \(\varphi=\pi\) grows for every \(h>0\). Rests on Theorem A61.8.
Derives Proposition A61.9. The only change is in the offset, so Lemma A61.1 is applied at \(\zeta=-\ii h/a\) instead of \(\tfrac{1}{2}-\ii h/a\), giving \(\sin\pi\zeta=-\ii\sinh s\) and \(\cos\pi\zeta=\cosh s\); hence \(S_{0}=-\pi^{2}/a^{2}\sinh^{2}s\), and \(S_{1}\), \(\tilde{S}_{1}\) come out real. At \(\psi=0\) they are equal, with \(\sqrt{X}=\pi^{2}\cosh s/a^{2}\sinh^{2}s\) and \(\Pi=+X\), so the second case of Proposition A61.7 applies and neutrality would require \(\abs{Q}\le\sqrt{X}\). Here \(Q=\left(\pi^{2}/a^{2}\right)\left(\tfrac{1}{2} +1/\sinh^{2}s\right)\), so the requirement reads \(\tfrac{1}{2}\sinh^{2}s+1\le\cosh s\), i.e.\ \(\left(\cosh s-1\right)^{2}\le0\) after using \(\sinh^{2}s=\cosh^{2}s-1\). That holds only at \(s=0\), where the two rows coincide. For every real separation the unstaggered arrangement therefore grows, and no spacing rescues it — which is why only the staggered pattern is ever photographed.
∎Four limitations are built into the statement proved above and none of them is repaired by any refinement of the algebra.
First, the stability is neutral. Theorem A61.8 says that at the critical ratio no linear mode grows; it does not say that any decays, and it says nothing at second order, where the street is in fact unstable. The observed consequence is that the photographed spacing drifts slowly downstream instead of locking to Equation (A61.1), and the agreement with the photographs — close to a few percent — is better than a neutral result has any right to expect.
Second, the vortices are points in an unbounded ideal fluid. Real vortices have cores of finite size, diffuse by Equation (35.52), and decay; the calculation has no viscosity in it anywhere.
Third, and most important, the cylinder never appears. The configuration Equation (A61.6) is postulated, not derived: nothing here explains why two boundary layers separating from a bluff body should roll up into two staggered rows, nor with what strength \(\Gamma\), nor at what rate. In particular nothing here predicts the Strouhal number, and Remark 35.62 says so at length: Equation (35.53) is a measurement, disciplined by the similarity argument of Phenomenon 35.49 into the form \(\mathrm{St}=\mathrm{St}(\mathrm{Re})\) and no further.
Fourth, the analysis is of a rigid infinite pattern already formed. The selection it performs is therefore of the same kind as the selection of a preferred wavelength in Section 35.8.1: it says which configurations can persist, not which one nature will build.
The Kármán Spacing Ratio of the Vortex Street discharges Equation (35.54) of Phenomenon 35.61 in Fluid Dynamics, and only that equation: Equation (35.53), the Strouhal number that shares the phenomenon with it, is reported there as measurement and is not derived here or anywhere else in this treatise. The two halves of that phenomenon are of different epistemic kinds and Remark 35.62 is where the difference is set out; the reader who has followed the algebra above should return to it, because the calculation just performed makes concrete how much idealization purchased the one number that could be derived. The measured record for the cylinder wake is collected in Experiment: Fluid Flow and Turbulence.