courseshub.world · Lives in Mathematics and Physics

Theodore von Kármán

Kármán Tódor (szőllőskislaki Kármán Tódor)

A course in two strands: the life of the Budapest prodigy who carried Göttingen aerodynamics to Aachen and Pasadena, founded the Jet Propulsion Laboratory and advised the US Air Force and NATO, and the physics of buckling, vortex streets, lattice vibrations, turbulence and high-speed flight that bears his name.

Born 11 May 1881, Budapest Died 6 May 1963, Aachen Vortex street Göttingen, 1911 Level advanced undergraduate → graduate
Kármán's model of the wake, computed live: two staggered rows of point vortices seen in the frame moving with them. Only the spacing \(h/l=\tfrac1\pi\operatorname{arsinh}1\approx0.281\) is free of linear instability.
Lesson 00

How to use this course

Theodore von Kármán was trained as an engineer and thought like an applied mathematician. His gift was to find the one simplified model of a messy problem that could be solved and that still captured the physics: a pair of vortex rows for the wake of a cylinder, a chain of springs for a crystal, a similarity hypothesis for turbulence, a mass balance for a rocket motor. He was also a builder of institutions on three continents. He directed the aeronautical institute at Aachen, then the Guggenheim laboratory at Caltech, co-founded the Jet Propulsion Laboratory and the Aerojet company, chaired the US Air Force's scientific advisers, and founded NATO's aeronautical research group.

The course follows both strands. Part I covers his formation in Budapest and Göttingen and his years at Aachen. Part II presents his main contributions to the physics of solids and fluids. Part III treats his American career, his advisory work and his legacy. The page on his teacher Ludwig Prandtl is a natural companion; several results here extend those of Prandtl's course.

Prerequisites for Part II
  • Vector calculus, ordinary differential equations and complex variables.
  • The Navier–Stokes equations, Bernoulli's equation and the idea of a boundary layer at the level of a first fluids course.
  • For Lesson 04, beam theory and linear elasticity.
  • For Lesson 06, the harmonic oscillator and elementary statistical mechanics (the Planck distribution).
  • For Lesson 09, the thermodynamics of an ideal gas and the speed of sound.

Each technical lesson opens with a paragraph for the general reader, then gives precise statements, derivations and exercises with folded solutions. The lessons run roughly in the order of his work: structures (1908–1932), the vortex street (1911–12), lattice dynamics (1912), laminar boundary layers (1921), turbulence (1921–1948), compressible and supersonic flow (1932–1956), and rockets and aeroelasticity (1940–41).

Reference

Timeline

1881Born 11 May in Budapest, son of Mór Kármán, professor of education and founder of the Minta Gymnasium, and Helene Konn.
1902Diploma in mechanical engineering from the Royal Joseph Polytechnic in Budapest. A year of service as an artillery cadet follows.
1903–06Assistant in hydraulics at the Polytechnic; consulting engineer for the Ganz works.
1906Fellowship of the Hungarian Academy of Sciences takes him to Göttingen and Prandtl.
1908Doctorate at Göttingen on the buckling strength of columns. In March, a visit to Paris and an early powered flight turn him toward aeronautics.
1909Habilitation at Göttingen; Privatdozent until 1912.
1910Inelastic buckling paper; encyclopedia article containing the Föppl–von Kármán plate equations.
1911–12Theory of the vortex street. Lattice dynamics with Max Born (1912). One term at the mining academy of Selmecbánya.
1913February: professor and director of the Aeronautical Institute at the Technische Hochschule Aachen.
1914–18Serves in the Austro-Hungarian army; at Fischamend near Vienna he helps design a tethered observation helicopter.
1919Heads the department of technical higher education in the Hungarian ministry of education during the revolutionary period; returns to Aachen after the fall of the Soviet Republic.
1921Über laminare und turbulente Reibung: the momentum integral, the rotating disk, the turbulent flat plate.
1922Organizes the first international congress of applied mechanics at Innsbruck, the forerunner of IUTAM.
1926First visit to Caltech, to advise on a wind tunnel.
1930Full-time director of the Guggenheim Aeronautical Laboratory at Caltech (GALCIT). Similarity theory of turbulence presented at Stockholm.
1932Kármán–Moore theory of supersonic slender bodies.
1936Supports Frank Malina's student rocket group at GALCIT.
1938Kármán–Howarth equation. Elected to the National Academy of Sciences.
1939Josiah Willard Gibbs Lecture of the American Mathematical Society.
1940–41Member of the board investigating the Tacoma Narrows Bridge collapse. Kármán–Tsien compressibility rule (1941). 12 August 1941: first American jet-assisted take-off.
1942March: co-founds the Aerojet Engineering Corporation, with Kármán as president.
1944The Jet Propulsion Laboratory takes its name, with Kármán as director. Chairs the Army Air Forces Scientific Advisory Group.
1946Foreign Member of the Royal Society; US Medal for Merit.
1947Similarity law of transonic flow.
1949Resigns his Caltech directorships; professor emeritus.
1952AGARD, NATO's Advisory Group for Aeronautical Research and Development, is created with Kármán as chairman.
1956Training centre for experimental aerodynamics near Brussels, later the von Karman Institute. Paper on aerodynamic heating from which the "Kármán line" was later named.
1957First Ludwig-Prandtl-Ring of the German aeronautical society.
1963February: first National Medal of Science, from President Kennedy. Dies 6 May in Aachen, aged 81.
Part I · Formation
Lesson 01

Budapest, 1881–1906

Theodore von Kármán was born Kármán Tódor in Budapest on 11 May 1881, into a Jewish family. His father, Mór Kármán, was a professor of philosophy and education at the University of Budapest and a reformer of Hungarian secondary schooling; he was ennobled in 1907, which gave the family the prefix "von" in German. His mother, Helene Konn, came from a family of scholars. Theodore was one of several children; his sister Josephine later kept house for him in Aachen and Pasadena.

At about six he could multiply large numbers in his head. His father, worried that he would become a performing prodigy, steered him away from mathematics for several years and had him taught at home in geography, history and poetry. At nine he entered the Minta Gymnasium, the model school his father had founded on his own educational principles, where pupils were led to derive rules from examples rather than memorize them. In his final year he won the Eötvös Prize, awarded by competition to the best student in mathematics and science in Hungary.

Engineering in Budapest

His father pushed him toward a practical profession, and he studied mechanical engineering at the Royal Joseph Polytechnic in Budapest, graduating in 1902. His diploma work analysed the motion of a heavy rod resting on a rounded end on a horizontal plane. After a year as an artillery cadet he became an assistant in hydraulics at the Polytechnic and consulted for the Ganz engineering works. In 1906 he published a paper on the buckling of long slender columns, and the same year a two-year fellowship from the Hungarian Academy of Sciences took him to Göttingen.

Lesson 02

Göttingen, 1906–1912

Göttingen in 1906 was the centre of world mathematics. Felix Klein was building applied science around it, and David Hilbert and Hermann Minkowski lectured there. Ludwig Prandtl had arrived two years earlier with the boundary-layer idea. Kármán joined Prandtl's group. His doctoral work, completed in 1908, combined theory and experiment on the buckling of columns stressed beyond the elastic limit, using a large hydraulic press (Lesson 04). The thesis was published in 1910 as Untersuchungen über Knickfestigkeit.

In March 1908 he made a short trip to Paris and watched one of the early powered flights near the city. The experience redirected his interests toward aeronautics. He habilitated at Göttingen in 1909 and taught there as a Privatdozent, lecturing on mechanics and thermodynamics.

Two famous papers

Kármán's best-known Göttingen story concerns a doctoral student of Prandtl, Karl Hiemenz, who was trying to measure the pressure distribution on a cylinder in a water channel. The flow behind the cylinder kept oscillating, and Hiemenz blamed his apparatus. Kármán suspected the oscillation was real and asked why the wake should prefer a staggered pattern of vortices. He calculated that symmetric vortex rows are always unstable and that staggered rows are stable only at one spacing ratio. The papers of 1911 and 1912 (the second with Henri Rubach) made the "Kármán vortex street" a standard object of fluid mechanics (Lesson 05).

The second paper took him into physics. With Max Born, then also a Privatdozent at Göttingen, he wrote in 1912 a theory of the vibrations of a crystal lattice that introduced the periodic boundary conditions now used throughout solid-state physics (Lesson 06). In the same period he wrote, for Klein's Encyklopädie der mathematischen Wissenschaften, a survey of strength problems in machine design containing the nonlinear plate equations now named after August Föppl and himself.

In 1912 he accepted a chair of applied mechanics at the mining academy of Selmecbánya (now Banská Štiavnica) in Hungary. He found it poorly equipped and returned to Göttingen after one term.

Lesson 03

Aachen, war and revolution, 1913–1929

In February 1913 the Technische Hochschule in Aachen appointed him professor of aeronautics and mechanics and director of its Aeronautical Institute. He built up equipment, including a wind tunnel, and a research group. The First World War interrupted the work. From 1914 he served in the Austro-Hungarian army. At the military aviation centre at Fischamend near Vienna he worked on aircraft and propellers, and with István Petróczy and Vilém Žurovec he took part in building a tethered helicopter intended to replace observation balloons; prototypes were built and flown in 1917–18.

At the end of the war he was in Budapest. During the revolutionary months of 1918–19 he headed the department of technical higher education in the ministry of education and worked on university reform. When the Hungarian Soviet Republic fell in August 1919 he left Hungary, and by the end of the year he was back in Aachen with his mother and his sister Josephine.

The Aachen school

The 1920s were the years in which Kármán became an international figure. Aachen students built competitive gliders for the Rhön contests, and Wolfgang Klemperer, who designed them, took his doctorate with Kármán. In 1921 Kármán published his paper on laminar and turbulent friction (Lessons 07 and 08). In September 1922 he organized at Innsbruck the first international congress of applied mechanics, which led to the International Union of Theoretical and Applied Mechanics. German aerodynamics now had two centres, Göttingen and Aachen, and the friendly rivalry between Prandtl and his former student, especially over turbulence, sharpened both.

In 1926 Robert Millikan invited him to Caltech to advise on the new Guggenheim wind tunnel. Over the next years he divided his time between Aachen and Pasadena. Rising antisemitism in Germany made the decision easier, and in 1930 he moved to Caltech full time, with his mother and sister.

Part II · The physics
Lesson 04

Buckling, plates and shells

For the general reader. A column pressed from both ends stays straight until a critical load, then suddenly bows sideways. Euler computed that load in the eighteenth century for perfectly elastic columns. Real columns made of steel often fail after the metal has started to yield, where Euler's formula overestimates their strength. Kármán's thesis measured and explained this inelastic buckling. He went on to write the equations that govern thin plates with large deflections, and later to explain why a buckled thin sheet in an aircraft can still carry load.

Goals
  • Derive the Euler load and the reduced-modulus load of an inelastic column.
  • State the Föppl–von Kármán equations for thin plates.
  • Use Kármán's effective-width formula for buckled plates.

Elastic and inelastic columns

A pinned column of length \(L\), Young's modulus \(E\) and second moment of area \(I\), under axial load \(P\), has lateral deflection \(w(x)\) with \(EIw''+Pw=0\), \(w(0)=w(L)=0\). A non-zero solution \(w=\sin(\pi x/L)\) exists first at the Euler load \(P_E=\pi^2EI/L^2\).

Above the proportional limit the stress–strain curve flattens, with tangent modulus \(E_t=d\sigma/d\varepsilon<E\). Friedrich Engesser proposed in 1889 to replace \(E\) by \(E_t\). He then noted that during bending one side of the column unloads elastically, with modulus \(E\), while the other side loads with modulus \(E_t\). Kármán worked out this double-modulus theory independently in his thesis and tested it on carefully made steel columns.

Reduced-modulus load (Engesser 1895, Kármán 1908–10)

Assume the column stays straight until it buckles at constant load, so that on bending the concave side loads with modulus \(E_t\) and the convex side unloads with modulus \(E\). Then the critical load is \(P_r=\pi^2E_rI/L^2\). For a rectangular section,

\[ E_r=\frac{4EE_t}{\big(\sqrt E+\sqrt{E_t}\big)^2}, \qquad E_t<E_r<E. \]

Sketch. In the buckled state the neutral axis shifts toward the unloading side, to the place where the resultant of the extra stresses vanishes. For a rectangle of depth \(d\), let the unloading part have depth \(d_1\) and the loading part \(d_2\), \(d_1+d_2=d\). Force balance gives \(Ed_1^2=E_td_2^2\), so \(d_1/d_2=\sqrt{E_t/E}\). The bending moment is then \(M=-E_rI\,w''\) with \(E_r=4(Ed_1^3+E_td_2^3)/d^3\), which simplifies to the formula above (Exercise 4.1). In 1947 F. R. Shanley showed that a real column starts to bend at the lower tangent-modulus load \(P_t=\pi^2E_tI/L^2\), with increasing load, and that its maximum load lies between \(P_t\) and \(P_r\). Kármán accepted the refinement.

The Föppl–von Kármán equations (1910)

A thin elastic plate of thickness \(h\) whose deflection \(w\) is comparable to \(h\) stretches in its own plane as it bends. Let \(\Phi\) be the Airy stress function of the in-plane (membrane) stresses, \(\sigma_{xx}=\Phi_{yy}\), \(\sigma_{yy}=\Phi_{xx}\), \(\sigma_{xy}=-\Phi_{xy}\), and let \(D=Eh^3/[12(1-\nu^2)]\) be the bending stiffness.

Föppl–von Kármán equations

Under a transverse load \(q(x,y)\) per unit area,

\[ \begin{aligned} D\,\nabla^4w-h\big(\Phi_{yy}w_{xx}+\Phi_{xx}w_{yy}-2\Phi_{xy}w_{xy}\big)&=q,\\ \nabla^4\Phi+E\big(w_{xx}w_{yy}-w_{xy}^2\big)&=0. \end{aligned} \]

The first equation balances bending against the transverse component of the membrane stresses. The second is the compatibility condition: the Gaussian curvature \(w_{xx}w_{yy}-w_{xy}^2\) of the deflected plate forces in-plane strain.

Linearizing about a flat plate under prescribed in-plane compression gives the classical buckling problem. For a simply supported plate of width \(b\), compressed along its length, the critical stress is

\[ \sigma_{cr}=k\,\frac{\pi^2D}{b^2h}=k\,\frac{\pi^2E}{12(1-\nu^2)}\Big(\frac hb\Big)^2,\qquad k\ge4, \]

with \(k=4\) for long plates. The nonlinear equations explain what happens next. Unlike a column, a buckled plate keeps carrying extra load, because the edges, held straight by the stiffeners, take up the stress.

Effective width (von Kármán, Sechler and Donnell, 1932)

A long, simply supported plate of width \(b\), loaded in compression beyond buckling until the stress at its supported edges reaches \(\sigma_e\), carries approximately the same total load as a strip of width

\[ b_e=b\sqrt{\frac{\sigma_{cr}}{\sigma_e}} \]

stressed uniformly to \(\sigma_e\). Equivalently, \(b_e\) is the width of a plate whose own buckling stress equals \(\sigma_e\).

The second form is Kármán's argument: since \(\sigma_{cr}\propto 1/b^2\), the width at which a plate buckles at \(\sigma_e\) is \(b\sqrt{\sigma_{cr}/\sigma_e}\). The formula became the basis of stressed-skin aircraft design.

Exercise 4.1Complete the reduced-modulus calculation for a rectangle of width \(c\) and depth \(d\). Show that \(E_r=4EE_t/(\sqrt E+\sqrt{E_t})^2\), and evaluate \(E_r/E\) when \(E_t=E/4\).

Show solution
With curvature \(\kappa\), the extra strain at distance \(s\) from the neutral axis is \(\kappa s\). Zero net extra force: \(E\kappa c\,d_1^2/2=E_t\kappa c\,d_2^2/2\), so \(\sqrt E\,d_1=\sqrt{E_t}\,d_2\) and \(d_1=d\sqrt{E_t}/(\sqrt E+\sqrt{E_t})\), \(d_2=d\sqrt E/(\sqrt E+\sqrt{E_t})\). The moment is \(M=\kappa c(Ed_1^3+E_td_2^3)/3\). Substituting, \(Ed_1^3+E_td_2^3=d^3EE_t(\sqrt{E_t}+\sqrt E)/(\sqrt E+\sqrt{E_t})^3=d^3EE_t/(\sqrt E+\sqrt{E_t})^2\). With \(I=cd^3/12\), \(M=E_rI\kappa\) where \(E_r=4EE_t/(\sqrt E+\sqrt{E_t})^2\). For \(E_t=E/4\): \(E_r=E^2/(1.5\sqrt E)^2=E/2.25\approx0.444E\), compared with \(E_t=0.25E\).

Exercise 4.2A long aluminium panel (\(E=70\) GPa, \(\nu=0.33\)) of width \(b=150\) mm and thickness \(h=1.5\) mm is simply supported on its long edges. Find \(\sigma_{cr}\). If the stiffeners reach \(\sigma_e=200\) MPa, what is the effective width?

Show solution
\(\sigma_{cr}=4\pi^2\cdot70\times10^9/(12\cdot0.891)\cdot(0.01)^2\approx25.8\) MPa. Then \(b_e=150\sqrt{25.8/200}\approx54\) mm. The panel carries about a third of the load of an unbuckled panel at 200 MPa, but far more than its buckling load alone would suggest.

Exercise 4.3Show that if \(w\) depends on \(x\) only (cylindrical bending), the compatibility equation forces \(\nabla^4\Phi=0\). Why is this the physically expected result?

Show solution
Then \(w_{yy}=w_{xy}=0\), so the Gaussian curvature \(w_{xx}w_{yy}-w_{xy}^2\) vanishes and \(\nabla^4\Phi=0\). A surface bent into a cylinder is developable: it can be made from a flat sheet without stretching, so bending alone produces no membrane strain. Doubly curved deflections, like a dome, cannot avoid stretching. This is why plates are stiffer than beams once deflections are large.
Lesson 05

The vortex street

For the general reader. Wind blowing past a wire or a chimney sheds swirls alternately from each side. They form two rows behind the body, with each swirl in one row opposite a gap in the other. The shedding makes telegraph wires sing and can shake structures. Kármán asked why the swirls choose this staggered pattern. His answer was a stability calculation: of all arrangements of two rows, only the staggered one with a particular ratio of width to spacing can persist.

Goals
  • Compute the velocity of a double row of point vortices.
  • State Kármán's stability result and the ratio \(h/l\approx0.281\).
  • Use the Strouhal number to estimate shedding frequencies.

Rows of vortices

A point vortex of circulation \(\Gamma\) at the origin has complex potential \(W=\frac{\Gamma}{2\pi i}\log z\). An infinite row of such vortices at \(z=nl\), \(n\in\mathbb Z\), has \(W=\frac{\Gamma}{2\pi i}\log\sin\frac{\pi z}{l}\), so the velocity it induces at a point \(z\) is

\[ u-iv=\frac{dW}{dz}=\frac{\Gamma}{2il}\cot\frac{\pi z}{l}. \]

A single row induces no velocity on its own members, by symmetry. Now add a second row of vortices of circulation \(-\Gamma\), displaced by a vertical distance \(h\). In the symmetric arrangement each vortex faces one of the other row; in the staggered arrangement the second row is also shifted by \(l/2\). In both cases the whole pattern translates rigidly, with speed

\[ V_{\text{sym}}=\frac{\Gamma}{2l}\coth\frac{\pi h}{l},\qquad V_{\text{stag}}=\frac{\Gamma}{2l}\tanh\frac{\pi h}{l}, \]

relative to the fluid at infinity. Mallock (1907) and Henri Bénard (1908) had observed such wakes; Kármán supplied the dynamics.

Theorem (Kármán, 1911–12)

For two infinite rows of point vortices of opposite circulation in an ideal fluid, the symmetric arrangement is linearly unstable for every \(h/l\). The staggered arrangement is linearly unstable unless

\[ \cosh\frac{\pi h}{l}=\sqrt2,\qquad\text{that is}\qquad \frac hl=\frac{\operatorname{arsinh}1}{\pi}=\frac{\ln(1+\sqrt2)}{\pi}\approx0.2806, \]

at which value it is neutrally stable to all small periodic disturbances. The street then moves at \(V=\Gamma/(2\sqrt2\,l)\) relative to the fluid at infinity.

Method. Displace the \(n\)th vortex of each row by small amounts proportional to \(e^{in\varphi}\), \(0\le\varphi\le2\pi\), and linearize the induced velocities. Each \(\varphi\) gives a \(4\times4\) linear system. Its eigenvalues come in pairs \(\pm\lambda\), with \(\lambda^2\) equal to \(\Gamma^2/l^4\) times an explicit function of \(\varphi\) and \(h/l\). Neutral stability requires \(\lambda^2\le0\) for every \(\varphi\). For the staggered street this happens only at \(\cosh(\pi h/l)=\sqrt2\). The result is subtle: a finer analysis shows that even this street is unstable to some nonlinear and finite-core perturbations. Its success lies elsewhere. The measured spacing ratio in real wakes is close to 0.28, and the model gave a first theory of the drag of bluff bodies from the momentum carried away by the vortices.

The Strouhal number

In 1878 Vincenc Strouhal found that the pitch of a wire singing in the wind is proportional to the wind speed divided by the wire diameter. In dimensionless form, the shedding frequency \(f\) behind a body of width \(D\) in a stream \(U\) satisfies

\[ \mathrm{St}=\frac{fD}{U}\approx0.2 \]

for a circular cylinder over a wide range of Reynolds numbers. Regular shedding begins at \(\mathrm{Re}=UD/\nu\approx47\), when the steady wake loses stability; above \(\mathrm{Re}\approx190\) the wake becomes three-dimensional, and near \(\mathrm{Re}\sim10^5\) shedding becomes irregular.

Exercise 5.1Using \(u-iv=\frac{\Gamma}{2il}\cot\frac{\pi z}{l}\), derive \(V_{\text{stag}}=\frac{\Gamma}{2l}\tanh\frac{\pi h}{l}\) and \(V_{\text{sym}}=\frac{\Gamma}{2l}\coth\frac{\pi h}{l}\).

Show solution
Put the row of strength \(-\Gamma\) at \(z=nl\) and look at a vortex of the other row. For the staggered street it sits at \(z_0=l/2+ih\). Then \(\cot(\pi z_0/l)=\cot(\pi/2+i\pi h/l)=-\tan(i\pi h/l)=-i\tanh(\pi h/l)\). The induced velocity is \(u-iv=\frac{-\Gamma}{2il}\cdot(-i\tanh\frac{\pi h}l)=\frac{\Gamma}{2l}\tanh\frac{\pi h}{l}\), real, so the motion is parallel to the rows. For the symmetric street, \(z_0=ih\) and \(\cot(i\pi h/l)=-i\coth(\pi h/l)\), giving \(\frac{\Gamma}{2l}\coth\frac{\pi h}{l}\). The same calculation for the other row gives the same velocity, so the pattern translates without deforming.

Exercise 5.2Verify that \(\cosh(\pi h/l)=\sqrt2\) gives \(h/l\approx0.2806\) and \(V=\Gamma/(2\sqrt2\,l)\approx0.354\,\Gamma/l\).

Show solution
\(\cosh x=\sqrt2\) gives \(\sinh x=1\), \(x=\ln(1+\sqrt2)\approx0.8814\), so \(h/l=0.8814/\pi\approx0.2806\). Then \(\tanh x=1/\sqrt2\) and \(V=\Gamma/(2\sqrt2\,l)\approx0.354\,\Gamma/l\). The hero figure uses this spacing and adds a uniform stream of speed \(V\), so that the vortices are at rest.

Exercise 5.3A telephone wire 3 mm thick is in a 10 m/s wind. Estimate the frequency of the tone it sings, and the Reynolds number (\(\nu=1.5\times10^{-5}\) m²/s).

Show solution
\(f\approx0.2\times10/0.003\approx670\) Hz, a clearly audible tone. \(\mathrm{Re}=10\times0.003/1.5\times10^{-5}=2000\), well inside the shedding range. The sound is the oscillating lift force on the wire, which alternates at the shedding frequency.
Lesson 06

Lattice vibrations with Born

For the general reader. A crystal is a regular array of atoms held together by forces that act like springs. Heat makes the atoms vibrate, and the vibrations travel through the crystal as waves. In 1912 Max Born and Kármán computed these waves for a real lattice and used them to explain how the heat capacity of a solid falls toward zero at low temperature. Their trick of wrapping a finite crystal around on itself, so that it has no surface, is still used in almost every calculation in solid-state physics.

Goals
  • State the Born–von Kármán periodic boundary condition and count normal modes.
  • Derive the dispersion relation of a monatomic chain.
  • Compute the vibrational density of states and the low-temperature heat capacity.

The problem in 1912

Einstein had explained in 1907 that the heat capacity of solids falls at low temperature because atomic vibrations are quantized. His model gave every atom the same frequency, and it predicted an exponential fall that experiment did not show. In 1912 Peter Debye treated the crystal as an elastic continuum with a cut-off frequency. Born and Kármán, in the same year, treated the actual lattice of atoms.

The monatomic chain

Consider \(N\) atoms of mass \(m\) on a line with spacing \(a\), nearest neighbours joined by springs of stiffness \(K\). With displacements \(u_n\),

\[ m\ddot u_n=K(u_{n+1}-2u_n+u_{n-1}). \]

A finite chain has ends, and the end atoms behave differently from the rest. Born and Kármán removed the ends by declaring \(u_{n+N}=u_n\): the chain is closed into a ring, or equivalently the crystal is repeated periodically.

Born–von Kármán boundary condition and mode counting

Impose \(u_{n+N}=u_n\). Then the normal modes are plane waves \(u_n\propto e^{i(kna-\omega t)}\) with

\[ k=\frac{2\pi j}{Na},\quad j=0,1,\dots,N-1\ (\text{equivalently }-\tfrac\pi a<k\le\tfrac\pi a),\qquad \omega(k)=2\sqrt{\frac Km}\,\Big|\sin\frac{ka}{2}\Big|. \]

There are exactly \(N\) modes. In three dimensions, for a crystal of \(N_1N_2N_3\) cells with \(s\) atoms per cell, the allowed wavevectors are \(\mathbf k=\sum_i (j_i/N_i)\mathbf b_i\) in one Brillouin zone, and there are \(3s\) branches \(\omega_\alpha(\mathbf k)\), giving \(3sN_1N_2N_3\) modes in all.

For long waves, \(\omega\approx a\sqrt{K/m}\,|k|\): ordinary sound with speed \(c=a\sqrt{K/m}\). Near the zone boundary the waves become dispersive and their group velocity falls to zero. Each mode is a quantum oscillator, so the thermal energy is

\[ E(T)=\sum_{\text{modes}}\frac{\hbar\omega}{e^{\hbar\omega/k_BT}-1}=\int_0^{\infty}g(\omega)\frac{\hbar\omega}{e^{\hbar\omega/k_BT}-1}\,d\omega, \]

with \(g(\omega)\) the density of states. Debye's model replaces \(g\) by its long-wave form up to a cut-off; the Born–von Kármán lattice gives the true \(g\), with sharp features (van Hove singularities) where the group velocity vanishes. At low temperature only long waves are excited, and both theories give a heat capacity \(C\propto T^3\) in three dimensions.

Exercise 6.1Derive \(\omega(k)=2\sqrt{K/m}\,|\sin(ka/2)|\) and show that \(k\) and \(k+2\pi/a\) describe the same motion of the atoms.

Show solution
Substituting \(u_n=Ae^{i(kna-\omega t)}\): \(-m\omega^2=K(e^{ika}-2+e^{-ika})=-2K(1-\cos ka)=-4K\sin^2(ka/2)\). So \(\omega^2=(4K/m)\sin^2(ka/2)\). Replacing \(k\) by \(k+2\pi/a\) multiplies \(u_n\) by \(e^{2\pi in}=1\): the displacements of the atoms, the only physical variables, are unchanged. This is why \(k\) can be restricted to one Brillouin zone.

Exercise 6.2Show that for large \(N\) the chain has density of states \(g(\omega)=\dfrac{2N}{\pi}\dfrac{1}{\sqrt{\omega_m^2-\omega^2}}\) on \(0<\omega<\omega_m=2\sqrt{K/m}\), and check that it integrates to \(N\).

Show solution
The allowed \(k\) are spaced by \(2\pi/(Na)\), so there are \(Na/(2\pi)\) per unit \(k\); both \(\pm k\) give the same \(\omega\). Thus \(g\,d\omega=2\cdot\frac{Na}{2\pi}\,dk\). From \(\omega=\omega_m\sin(ka/2)\), \(d\omega/dk=\frac a2\omega_m\cos(ka/2)=\frac a2\sqrt{\omega_m^2-\omega^2}\). Hence \(g=\frac{Na}{\pi}\cdot\frac{2}{a\sqrt{\omega_m^2-\omega^2}}=\frac{2N}{\pi\sqrt{\omega_m^2-\omega^2}}\). Since \(\int_0^{\omega_m}d\omega/\sqrt{\omega_m^2-\omega^2}=\pi/2\), the total is \(N\). The divergence at \(\omega_m\) is a van Hove singularity, which a continuum model misses.

Exercise 6.3Show that at low temperature the heat capacity of the chain is \(C\approx\dfrac{2\pi N k_B^2T}{3\hbar\omega_m}\), linear in \(T\). What power of \(T\) would you expect in two and three dimensions?

Show solution
For \(k_BT\ll\hbar\omega_m\) only small \(\omega\) matter, where \(g\approx2N/(\pi\omega_m)\). Then \(E\approx\frac{2N}{\pi\omega_m}\int_0^\infty\frac{\hbar\omega\,d\omega}{e^{\hbar\omega/k_BT}-1}=\frac{2N}{\pi\omega_m}\frac{(k_BT)^2}{\hbar}\frac{\pi^2}{6}=\frac{N\pi(k_BT)^2}{3\hbar\omega_m}\). Differentiating, \(C=\frac{2\pi Nk_B^2T}{3\hbar\omega_m}\). In \(d\) dimensions the low-frequency density of states grows like \(\omega^{d-1}\), so \(E\propto T^{d+1}\) and \(C\propto T^d\): \(T^2\) and \(T^3\).
Lesson 07

Boundary layers and the rotating disk

For the general reader. Prandtl's boundary-layer equations of 1904 were simpler than the full equations of viscous flow, but still hard to solve. In 1921 Kármán found a shortcut: instead of satisfying the equations at every point of the layer, satisfy only the balance of momentum for the layer as a whole. With a reasonable guess for the shape of the velocity profile, the problem reduces to one ordinary differential equation. Engineers used the method for forty years. In the same paper he solved exactly the flow next to a spinning disk, which pumps fluid like a fan.

Goals
  • Derive the momentum integral with a pressure gradient.
  • Apply the Kármán–Pohlhausen method and predict separation.
  • Reduce the rotating-disk flow to ordinary differential equations and compute the torque.

The momentum integral

Integrating Prandtl's equation \(uu_x+vu_y=U_eU_e'+\nu u_{yy}\) across the layer, using continuity to eliminate \(v\), gives an exact relation between integral thicknesses and the wall stress \(\tau_w\).

Kármán momentum integral (1921)

For steady two-dimensional incompressible boundary layers, with \(\delta^*=\int_0^\infty(1-u/U_e)\,dy\) and \(\theta=\int_0^\infty\frac u{U_e}(1-\frac u{U_e})\,dy\),

\[ \frac{d}{dx}\big(U_e^2\theta\big)+U_e\,\delta^*\frac{dU_e}{dx}=\frac{\tau_w}{\rho}. \]

It holds for laminar layers and, with \(u\) the mean velocity and \(\tau_w\) the mean wall stress, for turbulent ones.

Kármán's idea was to insert an assumed profile with one free parameter. In the same 1921 volume of the Zeitschrift für angewandte Mathematik und Mechanik, Ernst Pohlhausen used a quartic in \(\eta=y/\delta\),

\[ \frac{u}{U_e}=\big(2\eta-2\eta^3+\eta^4\big)+\frac{\Lambda}{6}\big(\eta-3\eta^2+3\eta^3-\eta^4\big),\qquad \Lambda=\frac{\delta^2}{\nu}\frac{dU_e}{dx}, \]

chosen to satisfy the wall compatibility condition \(\nu u_{yy}(0)=-U_eU_e'\) and smooth matching at \(y=\delta\). The wall shear is proportional to \(2+\Lambda/6\), so the method predicts separation at \(\Lambda=-12\). The Kármán–Pohlhausen method gives the flat-plate drag within a few per cent and locates separation roughly; it was the standard engineering tool until computers made marching the full equations routine.

The rotating disk

An infinite disk at \(z=0\) rotates with angular velocity \(\Omega\) in fluid at rest far away. Kármán sought solutions in cylindrical coordinates \((r,\phi,z)\) of the form

\[ v_r=r\Omega F(\zeta),\quad v_\phi=r\Omega G(\zeta),\quad v_z=\sqrt{\nu\Omega}\,H(\zeta),\quad p=p(z),\qquad \zeta=z\sqrt{\Omega/\nu}. \]
Kármán's rotating-disk flow (1921)

The Navier–Stokes equations reduce exactly to

\[ 2F+H'=0,\qquad F''=F^2-G^2+HF',\qquad G''=2FG+HG', \]

with \(F(0)=H(0)=0\), \(G(0)=1\), \(F,G\to0\) as \(\zeta\to\infty\). The solution (computed accurately by W. G. Cochran in 1934) has \(F'(0)=0.510\), \(G'(0)=-0.616\), \(H(\infty)=-0.884\).

The layer has constant thickness of order \(\sqrt{\nu/\Omega}\). Friction drags fluid round with the disk; centrifugal force throws it outward; and fluid is drawn axially toward the disk with speed \(0.884\sqrt{\nu\Omega}\) to replace it. It is an exact solution of the full Navier–Stokes equations with all three velocity components non-zero, and it models rotating machinery, viscometers and the boundary layer under a spinning fluid.

Exercise 7.1Use the cubic profile \(u/U=\tfrac32\eta-\tfrac12\eta^3\) (\(\eta=y/\delta\le1\)) and the momentum integral with \(U_e=U\) constant to find \(\delta(x)\) and \(c_f\). Compare with the Blasius values \(4.91\) and \(0.664\).

Show solution
\(\theta=\delta\int_0^1(\tfrac32\eta-\tfrac12\eta^3)(1-\tfrac32\eta+\tfrac12\eta^3)\,d\eta=\tfrac{39}{280}\delta\), and \(\tau_w=\tfrac32\mu U/\delta\). Then \(\tfrac{39}{280}\delta\,\delta'=\tfrac32\nu/U\), so \(\delta^2=\tfrac{840}{39}\nu x/U\) and \(\delta=4.64\,x/\sqrt{\mathrm{Re}_x}\). Finally \(c_f=\tau_w/(\tfrac12\rho U^2)=3\nu/(U\delta)=0.646/\sqrt{\mathrm{Re}_x}\), within 3% of Blasius.

Exercise 7.2Check that the Pohlhausen quartic has \(u/U_e=1\) and \(\partial u/\partial y=0\) at \(\eta=1\), and that the wall shear vanishes at \(\Lambda=-12\).

Show solution
At \(\eta=1\): \(2-2+1=1\) and \(1-3+3-1=0\), so \(u=U_e\). Derivatives: \(\frac{d}{d\eta}(2\eta-2\eta^3+\eta^4)=2-6\eta^2+4\eta^3\), zero at 1; \(\frac{d}{d\eta}(\eta-3\eta^2+3\eta^3-\eta^4)=1-6\eta+9\eta^2-4\eta^3\), zero at 1. At \(\eta=0\) the slope is \(2+\Lambda/6\), which vanishes at \(\Lambda=-12\).

Exercise 7.3Show that the torque needed to turn a disk of radius \(R\) (both faces wet) is \(M\approx0.616\,\pi\rho R^4\sqrt{\nu\Omega^3}\), neglecting edge effects. Evaluate it for \(R=0.1\) m, \(\Omega=100\) rad/s in water (\(\rho=1000\) kg/m³, \(\nu=10^{-6}\) m²/s).

Show solution
The azimuthal wall stress is \(\tau=\mu\,\partial v_\phi/\partial z|_0=\mu r\Omega G'(0)\sqrt{\Omega/\nu}\). One face: \(\int_0^R|\tau|\,r\,2\pi r\,dr=2\pi\mu\Omega\sqrt{\Omega/\nu}\,(0.616)R^4/4=\tfrac\pi2(0.616)\rho R^4\sqrt{\nu\Omega^3}\). Two faces double it. Numerically \(\sqrt{\nu\Omega^3}=\sqrt{10^{-6}\cdot10^6}=1\), so \(M\approx0.616\pi\times1000\times10^{-4}\approx0.19\) N·m.
Lesson 08

Turbulence

For the general reader. Turbulence was the problem Kármán returned to most often. He attacked it in three ways. First, with empirical laws and the momentum integral, he computed the drag of turbulent boundary layers for engineers. Second, he argued that turbulence near a wall looks the same at every scale, which leads to the logarithmic velocity profile and the constant that bears his name. Third, with Leslie Howarth, he wrote the first exact equation for the statistics of turbulence, the starting point of the modern theory.

Goals
  • Derive the 1921 turbulent flat-plate law from the one-seventh power profile.
  • Derive the logarithmic law from Kármán's similarity hypothesis.
  • State the Kármán–Howarth equation and extract the dissipation rate.

The turbulent flat plate (1921)

Using Heinrich Blasius's empirical pipe-friction law, Prandtl and Kármán independently inferred that turbulent velocity profiles follow roughly \(u/U=(y/\delta)^{1/7}\), with wall stress \(\tau_w=0.0225\,\rho U^2(\nu/U\delta)^{1/4}\). Inserting these in the momentum integral gives, for a plate turbulent from its leading edge,

\[ \delta\approx0.37\,x\,\mathrm{Re}_x^{-1/5},\qquad c_f\approx0.0577\,\mathrm{Re}_x^{-1/5},\qquad C_D\approx0.072\,\mathrm{Re}_L^{-1/5}, \]

valid for \(\mathrm{Re}\) up to about \(10^7\). Fitted to plate measurements, the drag constant becomes 0.074.

The similarity hypothesis (1930)

Prandtl's mixing length \(\ell\) (see his course) had to be specified by hand. Kármán proposed that, away from the wall, turbulence is locally similar: the mean flow near each point is described by \(\bar u'\) and \(\bar u''\) alone, and the only length they define is \(\bar u'/\bar u''\). Hence

\[ \ell=\kappa\left|\frac{\bar u'}{\bar u''}\right|,\qquad \tau_t=\rho\kappa^2\frac{(\bar u')^4}{(\bar u'')^2}. \]
Logarithmic law and friction laws (Kármán, Prandtl, 1930–32)

In the near-wall region where \(\tau_t\approx\tau_w\), with \(u_\tau=\sqrt{\tau_w/\rho}\), the similarity hypothesis gives

\[ \frac{\bar u}{u_\tau}=\frac1\kappa\ln\frac{yu_\tau}{\nu}+B,\qquad \kappa\approx0.41,\ B\approx5. \]

Integrated over a pipe or a plate, the log law yields friction laws of logarithmic form, valid at all large Reynolds numbers, such as the Kármán–Schoenherr law for a smooth flat plate, \(1/\sqrt{C_D}=4.13\log_{10}(\mathrm{Re}_LC_D)\).

Derivation. Setting \(\rho\kappa^2(\bar u')^4/(\bar u'')^2=\tau_w\) gives \((\bar u')^2=-u_\tau\bar u''/\kappa\) for a profile with \(\bar u''<0\). Put \(\bar u'=p(y)\): \(p'=-\kappa p^2/u_\tau\), so \(1/p=\kappa y/u_\tau+c\). With \(c=0\) (the profile measured from the wall), \(\bar u'=u_\tau/(\kappa y)\) and the logarithm follows. Kármán's similarity length coincides with Prandtl's \(\ell=\kappa y\) for this profile, which is why the two men reached the same law at the same time and argued about priority. The number \(\kappa\) is the Kármán constant.

Isotropic turbulence (1938)

G. I. Taylor had introduced in 1935 the statistical theory of homogeneous isotropic turbulence. Kármán and Howarth showed that isotropy reduces the two-point velocity correlation to one scalar function. Let \(u'^2=\overline{u_1^2}\), and define the longitudinal correlations

\[ \overline{u_1(\mathbf x)\,u_1(\mathbf x+r\mathbf e_1)}=u'^2f(r,t),\qquad \overline{u_1^2(\mathbf x)\,u_1(\mathbf x+r\mathbf e_1)}=u'^3k(r,t). \]
Kármán–Howarth equation (1938)

For homogeneous isotropic turbulence in an incompressible fluid,

\[ \frac{\partial}{\partial t}\big(u'^2f\big)-u'^3\Big(\frac{\partial k}{\partial r}+\frac{4k}{r}\Big)=2\nu u'^2\Big(\frac{\partial^2f}{\partial r^2}+\frac4r\frac{\partial f}{\partial r}\Big). \]

The equation is exact but not closed: it links the double correlation to the unknown triple correlation \(k\), the closure problem in its simplest form. Kármán and Howarth looked for self-preserving solutions in which \(f\) keeps its shape as the turbulence decays. Andrey Kolmogorov derived his four-fifths law (1941) from the same equation, rewritten for structure functions. In 1948 Kármán proposed an interpolation for the energy spectrum, \(E(k)\propto k^4/(1+k^2/k_e^2)^{17/6}\), which behaves like \(k^4\) for large eddies and like Kolmogorov's \(k^{-5/3}\) in the inertial range; it is still used to generate synthetic turbulence and in atmospheric gust models.

Exercise 8.1Derive \(\delta\approx0.37\,x\,\mathrm{Re}_x^{-1/5}\) from the one-seventh power profile and the wall-stress law above.

Show solution
For \(u/U=\eta^{1/7}\), \(\theta=\delta\int_0^1\eta^{1/7}(1-\eta^{1/7})d\eta=\delta(\tfrac78-\tfrac79)=\tfrac7{72}\delta\). The momentum integral \(U^2\theta'=\tau_w/\rho\) gives \(\tfrac7{72}\delta'=0.0225(\nu/U)^{1/4}\delta^{-1/4}\). Integrate from \(\delta(0)=0\): \(\tfrac45\delta^{5/4}=\tfrac{72}{7}(0.0225)(\nu/U)^{1/4}x\), so \(\delta^{5/4}=0.289(\nu/U)^{1/4}x\) and \(\delta=0.289^{4/5}(\nu/U)^{1/5}x^{4/5}=0.37\,x\,\mathrm{Re}_x^{-1/5}\).

Exercise 8.2Solve the Kármán–Schoenherr law for a ship hull with \(\mathrm{Re}_L=10^7\) by fixed-point iteration, and compare with the one-seventh-power formula \(0.074\,\mathrm{Re}_L^{-1/5}\).

Show solution
Iterate \(C_D\leftarrow\big(0.242/\log_{10}(\mathrm{Re}_LC_D)\big)^2\), since \(1/4.13=0.242\). Starting from 0.003: \(\log_{10}(3\times10^4)=4.48\), \(C_D=(0.242/4.48)^2=0.00292\); next \(\log_{10}(2.92\times10^4)=4.47\), \(C_D=0.00293\). The power law gives \(0.074\times10^{-1.4}=0.00295\). The two agree here; at \(\mathrm{Re}_L=10^9\) the power law is badly off, while the logarithmic law, built on the log profile, remains accurate.

Exercise 8.3Near \(r=0\), \(f=1-r^2/(2\lambda^2)+O(r^4)\) (Taylor microscale \(\lambda\)) and \(k=O(r^3)\). Show from the Kármán–Howarth equation that the dissipation rate of kinetic energy per unit mass is \(\varepsilon=15\nu u'^2/\lambda^2\).

Show solution
At \(r=0\), \(f=1\), so the left side's first term is \(d(u'^2)/dt\). Since \(k\sim r^3\), \(\partial_rk+4k/r\to0\). On the right, \(f''(0)=-1/\lambda^2\) and \(f'/r\to-1/\lambda^2\), so the bracket is \(-5/\lambda^2\). Thus \(du'^2/dt=-10\nu u'^2/\lambda^2\). The kinetic energy per unit mass is \(\tfrac32u'^2\), so \(\varepsilon=-\tfrac32\,du'^2/dt=15\nu u'^2/\lambda^2\).
Lesson 09

High-speed aerodynamics

For the general reader. As aircraft approached the speed of sound in the 1930s, designers needed to know how compressibility changes the forces on a wing. Kármán and his student Qian Xuesen (Hsue-shen Tsien) produced a correction formula that worked better than Prandtl's simpler rule. Kármán also computed the wave drag of slender supersonic bodies such as shells and rockets, and found the scaling law that organizes transonic flight, the awkward regime just below and above Mach 1.

Goals
  • Compare the Prandtl–Glauert and Kármán–Tsien compressibility rules.
  • Use linearized supersonic theory and Kármán's slender-body drag formula.
  • Derive the transonic similarity parameter.

Subsonic compressibility (1939–41)

The Prandtl–Glauert rule \(C_p=C_{p,0}/\sqrt{1-M_\infty^2}\) comes from linearizing the potential equation. It underestimates suction peaks at higher subsonic Mach numbers. Tsien (1939) and Kármán (1941) used the hodograph method: in the plane of velocity components the equations are linear, and replacing the true adiabatic pressure–density curve by its tangent at free-stream conditions (formally \(\gamma=-1\)) makes them solvable in closed form.

Kármán–Tsien rule

For a thin aerofoil at subsonic free-stream Mach number \(M_\infty\), with \(\beta=\sqrt{1-M_\infty^2}\) and \(C_{p,0}\) the incompressible pressure coefficient at the same point,

\[ C_p=\frac{C_{p,0}}{\beta+\dfrac{M_\infty^2}{1+\beta}\,\dfrac{C_{p,0}}{2}}. \]

For small \(C_{p,0}\) it reduces to Prandtl–Glauert; for suction (\(C_{p,0}<0\)) it predicts a larger correction, in better agreement with experiment up to the critical Mach number.

Supersonic slender bodies (1932)

Above Mach 1 the linearized equation \((M_\infty^2-1)\phi_{xx}=\phi_{yy}+\phi_{zz}\) is a wave equation, and disturbances are confined to Mach cones. Kármán and Norton Moore represented an axisymmetric body by a line of sources along its axis, the supersonic analogue of Rankine's method, and computed pressure and drag on projectiles.

Kármán's slender-body wave drag

For a pointed slender body of length \(\ell\) with cross-sectional area \(S(x)\), closed or ending in a cylinder (\(S'(0)=S'(\ell)=0\)), the linearized supersonic wave drag is

\[ D=-\frac{\rho U^2}{4\pi}\int_0^\ell\!\!\int_0^\ell S''(x)\,S''(\xi)\,\ln|x-\xi|\,dx\,d\xi, \]

independent of Mach number to this order.

The formula depends only on the area distribution, so two bodies with the same \(S(x)\) have the same wave drag. This is the root of the area rule found by Richard Whitcomb in 1952, and the shapes of least drag for given volume (Sears–Haack bodies) or given base area (the von Kármán ogive) come from minimizing it.

Transonic similarity (1947)

Near \(M_\infty=1\) the linear term \((1-M_\infty^2)\phi_{xx}\) is small, and the leading nonlinearity must be kept:

\[ \Big[1-M_\infty^2-(\gamma+1)\frac{\phi_x}{U}\Big]\phi_{xx}+\phi_{yy}=0, \]

with \(\phi_y=U\tau g'(x)\) on a body of thickness ratio \(\tau\) and shape \(g\) (here \(M_\infty\approx1\) in the coefficient of the nonlinear term).

Kármán's transonic similarity law

For affinely related thin aerofoils, \(\dfrac{(\gamma+1)^{1/3}}{\tau^{2/3}}\,C_p\) is a function of \(x/c\) and of the single parameter

\[ K=\frac{1-M_\infty^2}{\big[(\gamma+1)\tau\big]^{2/3}}. \]

Sketch. Put \(y=Y/a\) and \(\phi=Ub\,\Phi(x,Y)\). The equation becomes \((1-M_\infty^2)b\Phi_{xx}-(\gamma+1)b^2\Phi_x\Phi_{xx}+a^2b\Phi_{YY}=0\) and the boundary condition \(ab\,\Phi_Y=\tau g'\). Choosing \(ab=\tau\) and \((\gamma+1)b=a^2\) removes \(\tau\) and \(\gamma\) from all but one coefficient: \(a=[(\gamma+1)\tau]^{1/3}\), \(b=\tau^{2/3}/(\gamma+1)^{1/3}\), and the remaining parameter is \((1-M_\infty^2)/a^2=K\). Then \(C_p=-2\phi_x/U=-2b\Phi_x\).

Exercise 9.1At a point on an aerofoil \(C_{p,0}=-0.5\). Find \(C_p\) at \(M_\infty=0.7\) by Prandtl–Glauert and by Kármán–Tsien.

Show solution
\(\beta=\sqrt{0.51}=0.714\). Prandtl–Glauert: \(-0.5/0.714=-0.700\). Kármán–Tsien: denominator \(0.714+\frac{0.49}{1.714}(-0.25)=0.714-0.0715=0.643\), so \(C_p=-0.778\). The larger suction brings the local flow closer to sonic, so the critical Mach number is lower than the linear rule suggests.

Exercise 9.2In linearized supersonic flow the pressure coefficient on a surface inclined at small angle \(\theta\) to the stream is \(C_p=\pm2\theta/\sqrt{M_\infty^2-1}\) (Ackeret, 1925). Find \(C_L\) and the wave-drag coefficient of a flat plate at angle of attack \(\alpha\). Why is this drag present even with zero thickness?

Show solution
The lower surface has \(C_p=+2\alpha/\sqrt{M^2-1}\), the upper \(-2\alpha/\sqrt{M^2-1}\). The net normal force coefficient is \(4\alpha/\sqrt{M^2-1}\); resolved, \(C_L=4\alpha/\sqrt{M^2-1}\) and \(C_D=C_L\alpha=4\alpha^2/\sqrt{M^2-1}\). Making lift requires turning the flow, and in supersonic flow the turning radiates Mach waves that carry momentum to infinity: wave drag due to lift.

Exercise 9.3Kármán asked at what height an aircraft must fly so fast that it is effectively in orbit. Model the atmosphere as \(\rho=\rho_0e^{-z/H}\) with \(\rho_0=1.2\) kg/m³, \(H=7\) km. An aircraft with wing loading \(W/S=3000\) N/m² and \(C_L=1\) flies at orbital speed \(v_o=\sqrt{gR}\) (\(R=6.37\times10^6\) m). At what altitude does the lift at that speed just equal the weight?

Show solution
Lift equals weight when \(\tfrac12\rho v_o^2C_L=W/S\), so \(\rho=2(W/S)/(C_LgR)=6000/(9.81\times6.37\times10^6)\approx9.6\times10^{-5}\) kg/m³. Then \(z=H\ln(\rho_0/\rho)=7\ln(12500)\approx66\) km. Above this height a vehicle at orbital speed needs no wings. The answer depends on the vehicle, which is why the line is a convention; the round figure of 100 km was adopted by the FAI. Kármán's own 1956 estimate, about 84 km, came from limits set by aerodynamic heating.
Lesson 10

Rockets and structures in the wind

For the general reader. In the late 1930s rockets were the hobby of enthusiasts, and many engineers thought them unreliable. Solid-fuel rockets in particular tended either to fizzle or to explode. Kármán and his student Frank Malina showed with a short calculation that a solid motor can burn steadily, and the first American rocket-assisted take-off followed in 1941. In the same years Kármán was called to explain why the Tacoma Narrows Bridge had shaken itself apart in a moderate wind, a problem in what he called aeroelasticity.

Goals
  • Derive the equilibrium pressure of a solid rocket motor and its stability condition.
  • Understand why a burning-rate exponent below one is essential.
  • Distinguish vortex-induced vibration from flutter.

The solid-propellant motor (1940)

Propellant of density \(\rho_p\) burns on a surface of area \(A_b\), receding at a rate that depends on chamber pressure, empirically \(r=ap^n\). Gas leaves through a nozzle throat of area \(A_t\) at the choked mass flow \(pA_t/c^*\), where the characteristic velocity \(c^*\) depends on the propellant chemistry. Malina and Kármán's analysis, in modern notation, is a mass balance in the chamber.

Equilibrium chamber pressure

Neglecting the gas stored in the chamber, steady burning requires \(\rho_pA_b\,ap^n=pA_t/c^*\), so

\[ p_{eq}=\Big(\rho_p\,a\,c^*\,\frac{A_b}{A_t}\Big)^{1/(1-n)}. \]

This equilibrium is stable if and only if \(n<1\). A motor burning on a surface of constant area, as in an end-burning grain, then gives constant pressure and thrust for the whole burn.

Stability. Including the gas mass \(m=pV/(RT)\) in the chamber, \(\dot m=\rho_pA_bap^n-pA_t/c^*\). For a small perturbation \(p=p_{eq}(1+\epsilon)\), the generation term grows like \(1+n\epsilon\) and the discharge like \(1+\epsilon\), so \(\dot\epsilon\propto(n-1)\epsilon\). With \(n<1\) a pressure rise increases discharge faster than generation and dies away. With \(n>1\) it feeds on itself, and the motor explodes. The GALCIT group's first jet-assisted take-off, on 12 August 1941 at March Field, used solid units built on this principle. In 1942 Jack Parsons devised a castable asphalt and perchlorate propellant, a forerunner of modern composite propellants.

Aeroelasticity and Tacoma Narrows (1940)

On 7 November 1940 the Tacoma Narrows suspension bridge twisted itself to destruction in a wind of about 19 m/s. Kármán, with Othmar Ammann and Glenn Woodruff, formed the board of engineers that reported on the failure in March 1941. Vortex shedding (Lesson 05) was the first suspect. Later analysis showed that the decisive motion was torsional flutter: once the deck began to twist, the aerodynamic moment acquired a component in phase with the angular velocity, acting as negative damping. A minimal model is

\[ I\ddot\theta+(c_s-c_a(U))\,\dot\theta+k\theta=0, \]

where \(c_s\) is the structural damping and \(c_a(U)\) the aerodynamic contribution, increasing with wind speed. Oscillations grow without bound when \(c_a(U)>c_s\). After the collapse, wind-tunnel testing of bridge models became standard practice, and Kármán, with his student Ernest Sechler and others at Caltech, developed aeroelasticity into a discipline for aircraft wings as well as bridges.

Exercise 10.1A propellant has \(n=0.7\). By how much does the chamber pressure rise if a crack increases the burning area by 10%? Repeat for \(n=0.3\).

Show solution
\(p_{eq}\propto A_b^{1/(1-n)}\). For \(n=0.7\), \(1.1^{1/0.3}=1.1^{3.33}\approx1.37\): a 37% rise. For \(n=0.3\), \(1.1^{1/0.7}\approx1.15\). Low exponents make motors tolerant of defects, which is why propellant chemists aim for small \(n\).

Exercise 10.2The Tacoma deck was about 2.4 m deep. Taking \(\mathrm{St}\approx0.2\) as a rough value, estimate the vortex-shedding frequency at 19 m/s, and compare it with the torsional oscillation observed on the last day, roughly 0.2 Hz. What does the comparison suggest?

Show solution
\(f\approx0.2\times19/2.4\approx1.6\) Hz, about eight times the frequency of the destructive twisting. Simple resonance with periodic shedding cannot explain a large motion at 0.2 Hz. The motion was self-excited: the twisting deck organized its own wake, and the aerodynamic forces fed energy into the torsional mode, as in the negative-damping model.

Exercise 10.3In the model \(I\ddot\theta+(c_s-\beta U)\dot\theta+k\theta=0\), find the critical wind speed for flutter and the frequency of oscillation at onset.

Show solution
Solutions \(e^{st}\) satisfy \(Is^2+(c_s-\beta U)s+k=0\). The real parts of the roots are \(-(c_s-\beta U)/(2I)\), which become positive when \(U>U_c=c_s/\beta\). At onset the damping term vanishes and \(s=\pm i\sqrt{k/I}\): the structure oscillates at its natural torsional frequency. Doubling the structural damping doubles \(U_c\); stiffening the deck in torsion, as in the replacement bridge, raises both \(k\) and the aerodynamic stability.
Part III · Life and legacy
Lesson 11

Caltech, JPL and the Air Force, 1930–1949

At Caltech Kármán directed the Guggenheim Aeronautical Laboratory (GALCIT) and its ten-foot wind tunnel. He gave American aeronautics what Göttingen had given Germany: theory joined to measurement. His seminars drew aircraft designers from the Los Angeles industry, and his doctoral students formed a generation of leaders. They included Maurice Biot (1932), Frank Wattendorf (1933), Ernest Sechler (1934), Qian Xuesen (1939), Frank Malina (1940), Homer Joe Stewart (1940), Chia-Chiao Lin (1944), Wallace Hayes (1947) and Frank Marble (1947). In 1933 he helped found the Institute of the Aeronautical Sciences. He was elected to the National Academy of Sciences in 1938 and gave the Gibbs Lecture of the American Mathematical Society in 1939.

From student rockets to JPL

In 1936 Malina, with Jack Parsons, Edward Forman and other students, asked to test rocket motors as a GALCIT project. Kármán gave them permission and his support, and the group fired its first motor in the Arroyo Seco above Pasadena that autumn. The Army Air Corps funded jet-assisted take-off research from 1939, the first American JATO flight took place on 12 August 1941, and in March 1942 Kármán, Malina, Parsons, Forman, Martin Summerfield and the lawyer Andrew Haley founded the Aerojet Engineering Corporation, with Kármán as president. In 1943 Kármán and his colleagues analysed intelligence on the German long-range rocket programme and proposed a missile research programme. In 1944 the Arroyo Seco group became the Jet Propulsion Laboratory, with Kármán as its first director; when he left for Washington at the end of that year, Malina succeeded him.

Adviser to the Air Force

In 1944 General Henry Arnold asked Kármán to chair a Scientific Advisory Group for the Army Air Forces. In 1945 Kármán led a mission to Europe that interviewed German scientists, including his old teacher Prandtl in Göttingen, and inspected wind tunnels and rocket sites. The group's multi-volume report, Toward New Horizons, set out a programme of research in jet propulsion, supersonic flight, guided missiles and electronics, and is often credited with shaping the post-war US Air Force's commitment to science. Kármán chaired the group's successor, the Air Force Scientific Advisory Board. He received the US Medal for Merit in 1946 and resigned his Caltech directorships in 1949.

DiscussionKármán moved between pure research, company founding and military advice with little sense of conflict. Using this lesson and Lesson 10, discuss what this mixture made possible, what it put at risk, and how universities now try to manage similar roles.

Lesson 12

Statesman of aeronautics, 1949–1963

After 1949 Kármán spent much of his time in Europe and on international organizations. He persuaded NATO to create the Advisory Group for Aeronautical Research and Development (AGARD), established in 1952, and chaired it until his death. In 1956 he founded near Brussels a training centre for experimental aerodynamics, renamed the von Karman Institute for Fluid Dynamics after his death; in the same year he helped found the International Council of the Aeronautical Sciences, and in 1960 the International Academy of Astronautics. His book Aerodynamics: Selected Topics in the Light of Their Historical Development (1954) remains a lucid introduction.

The postwar years also brought a bitter episode. His student Qian Xuesen, a founder of JPL and a colonel in the US Air Force during the 1945 mission, lost his security clearance in 1950 amid accusations of communist sympathies, was held under restriction for five years and was deported to China in 1955. There he led the Chinese missile and space programmes. Kármán supported him throughout.

In February 1963 President John F. Kennedy presented Kármán with the first National Medal of Science. He died on a visit to Aachen on 6 May 1963, five days before his 82nd birthday. (Some sources give 7 May.) He never married. His body was returned to the United States and entombed in the Beth Olam Mausoleum at what is now Hollywood Forever Cemetery in Los Angeles. His autobiography, The Wind and Beyond, written with Lee Edson, appeared in 1967.

Legacy

  • Fluid mechanics: the vortex street, the momentum integral, the rotating-disk flow, the Kármán constant, the Kármán–Howarth equation and the von Kármán spectrum are standard textbook material.
  • Solids and physics: the Föppl–von Kármán equations govern thin elastic sheets, from aircraft panels to crumpled paper and graphene; Born–von Kármán boundary conditions underlie band theory and phonon calculations.
  • Aeronautics and space: compressibility corrections, transonic similarity and slender-body theory guided the design of high-speed aircraft; the Kármán line is the conventional edge of space.
  • Institutions: the Jet Propulsion Laboratory, Aerojet, AGARD, the von Karman Institute and the international congresses of applied mechanics all carry his imprint.
  • Honours: Foreign Member of the Royal Society (1946), Franklin Medal (1948), Daniel Guggenheim Medal (1955), first Ludwig-Prandtl-Ring (1957), first National Medal of Science (1963); craters on the Moon and Mars, and a 1992 US postage stamp.

Final projectChoose one model from Part II (the vortex street, the Born–von Kármán chain, the rotating disk, the similarity hypothesis or the rocket mass balance). Write a two-page account for a strong undergraduate that states the model, solves it, and explains which features of the real system it keeps and which it throws away. End with one experiment that tests it.

Reference

Connections

People linked to Theodore von Kármán across this collection. A name without a link has no page here yet. See Theodore von Kármán in the network of all 337 lives →

Teachers and advisors

Students

  • Chia-Chiao Lin 1944
  • Ernest Sechler 1934
  • Frank E. Marble 1947
  • Frank Malina 1940
  • Frank Wattendorf 1933
  • Homer Joe Stewart 1940
  • Maurice Anthony Biot 1932
  • Qian Xuesen 1939
  • Wallace D. Hayes 1947
  • Wolfgang Klemperer 1926

Influenced by

Influenced

Colleagues and collaborators

Family

  • Mór Kármán

Also mentioned in this life

Fields

Reference

Sources and further reading

Biography and history

  1. J. J. O'Connor and E. F. Robertson, Theodore von Kármán, MacTutor History of Mathematics.
  2. Theodore von Kármán, Wikipedia (doctoral students, honours).
  3. Kármán, Theodor von, Neue Deutsche Biographie.
  4. Toward New Horizons: Theodore von Kármán, in memoriam, Air Force Magazine, June 1963.
  5. U.S. Centennial of Flight Commission, Early rocketry: the GALCIT project.
  6. Kármán Tódor, Fizikai Szemle 2011, no. 7–8 (in Hungarian).
  7. T. von Kármán with L. Edson, The Wind and Beyond, Little, Brown, 1967.
  8. M. H. Gorn, The Universal Man: Theodore von Kármán's Life in Aeronautics, Smithsonian Institution Press, 1992.
  9. H. L. Dryden, Theodore von Kármán 1881–1963, Biographical Memoirs of the National Academy of Sciences 38 (1965).

Physics

  1. Th. von Kármán, Untersuchungen über Knickfestigkeit, Mitteilungen über Forschungsarbeiten des VDI 81 (1910).
  2. Th. von Kármán, Festigkeitsprobleme im Maschinenbau, Encyklopädie der mathematischen Wissenschaften IV.4 (1910).
  3. Th. von Kármán, Über den Mechanismus des Widerstandes, den ein bewegter Körper in einer Flüssigkeit erfährt, Nachr. Ges. Wiss. Göttingen (1911, 1912).
  4. M. Born and Th. von Kármán, Über Schwingungen in Raumgittern, Physikalische Zeitschrift 13 (1912).
  5. Th. von Kármán, Über laminare und turbulente Reibung, ZAMM 1 (1921).
  6. Th. von Kármán, Mechanische Ähnlichkeit und Turbulenz, Nachr. Ges. Wiss. Göttingen (1930).
  7. Th. von Kármán, E. E. Sechler and L. H. Donnell, The strength of thin plates in compression, Trans. ASME 54 (1932).
  8. Th. von Kármán and N. B. Moore, Resistance of slender bodies moving with supersonic velocities, Trans. ASME 54 (1932).
  9. Th. von Kármán and L. Howarth, On the statistical theory of isotropic turbulence, Proc. R. Soc. Lond. A 164 (1938).
  10. Th. von Kármán, Compressibility effects in aerodynamics, J. Aeronaut. Sci. 8 (1941).
  11. Th. von Kármán, The similarity law of transonic flow, J. Math. Phys. 26 (1947).
  12. Th. von Kármán, Aerodynamics: Selected Topics in the Light of Their Historical Development, Cornell University Press, 1954 (Dover reprint).
  13. H. Schlichting and K. Gersten, Boundary-Layer Theory, 9th ed., Springer, 2017 (momentum integral, rotating disk, turbulent friction laws).
  14. H. Lamb, Hydrodynamics, 6th ed., Cambridge, 1932 (vortex rows and the vortex street).
  15. N. W. Ashcroft and N. D. Mermin, Solid State Physics, 1976, chapters 8 and 22 (Born–von Kármán conditions, lattice dynamics).