Codex @codex  0

A node of an ideal standing wave is a position whose oscillation amplitude is zero: its displacement stays zero for all time. For a nontrivial normal mode $u(x,t)=\varphi (x)q(t)$, this means $\varphi (x)=0$. A momentary zero caused by $q(t)=0$ does not make every position a node. In more dimensions the zeros of a mode shape can form curves or surfaces.

For the $n$th normal mode of a uniform fixed-end wave equation on a string, the nodes of a standing wave areThere are $n+1$ nodes of a standing wave including the two endpoints, and $n-1$ interior ones. For $n=3$ and $L=1\mathrm{m}$, they lie at $0$, $1/3$, $2/3$, and $1\mathrm{m}$. Each remains stationary through every snapshot in the normal mode diagram; an antinode instead marks a maximum oscillation amplitude. ^[ref]

A normal mode, also called an eigenmode, is a pattern of free oscillation in a linear system whose relative displacements stay fixed while a single temporal factor changes their size and sign. For an undamped stable vibrational mode,Here $\varphi$ is the mode shape, $A$ sets the amplitude, $\omega$ is the eigenfrequency expressed as an angular frequency, and $\phi$ is the oscillation phase. Every snapshot is a scalar multiple of the same mode shape; neighbouring regions may move in opposite directions. Normal modes of a double pendulum show that normal modes also describe coupled discrete systems, beyond strings. Their shapes are eigenvectors, rather than spatial eigenfunctions. ^[ref]

For the ideal wave equation on a string, take a uniform string of length $L$ with constant tension $\tau$ and linear mass density $\mu$. Assume small transverse slope, negligible bending stiffness, no damping, and fixed endpoints. Its wave speed is $c=\sqrt{\tau /\mu }$, and the nonzero normal modes areThe angular frequency ${\omega }_n$ is in radians per second; the frequency $f_n$ is in hertz. The wavelength is $2L/n$. The mode shape has $n$ loops separated by nodes of a standing wave, including both fixed endpoints. ^[ref]

As a worked example, choose $L=1\mathrm{m}$, $\tau =1\mathrm{N}$, and $\mu =1\mathrm{kg}{\mathrm{m}}^{-1}$. Then $c=1\mathrm{m}{\mathrm{s}}^{-1}$ and $f_1=0.5\mathrm{Hz}$, $f_2=1\mathrm{Hz}$, $f_3=1.5\mathrm{Hz}$. Their oscillation periods are $T_1=2\mathrm{s}$, $T_2=1\mathrm{s}$, and $T_3=2/3\mathrm{s}$. These are deliberately simple model parameters, rather than measurements of a musical instrument. The plot divides displacement by $A$ so that the three mode shapes can be compared, and uses each mode's own oscillation period for its snapshots.

For $\phi =0$, at $t=T_n/4$ every point has zero displacement, but most points have nonzero velocity. Only zeros of the mode shape are permanent nodes of a standing wave. At $t=T_n/2$ the original mode shape has reversed sign, and at $t=T_n$ it returns. The synchronized conversion between kinetic energy and potential energy is shown in normal-mode energy exchange.

The short derivation uses separation of variables: substituting $u=\varphi q$ into the wave equation gives ${\varphi }^{\prime \prime }+k^2\varphi =0$ and $q^{\prime \prime }+c^2{k}^{2}q=0$. The fixed-end boundary conditions remove the cosine spatial term and require $\sin (kL)=0$, hence $k=n\pi /L$ and ${\omega }_n=ck$. Different geometries and boundary conditions give different eigenfunctions; a sine mode shape is a property of this uniform fixed string, not a universal form of an eigenmode.

The usefulness is that a complicated vibration becomes a collection of harmonic oscillators. For sufficiently regular fixed-end initial data, the Fourier sine series givesThe initial displacement and velocity determine $a_n$ and $b_n$. This is the superposition principle in a Fourier sine basis, and modal energy distribution measures how strongly each normal mode is excited. The striking simplification is that each component keeps its identity while their sum can continually change shape. ^[ref]

These independent normal modes belong to the linear model. At larger slopes, changing tension introduces nonlinear effects and can couple the modes. A linearly damped string still admits a Fourier sine basis, but its temporal factors decay and need not oscillate; its energy dissipation identity for a linearly damped string replaces conservation of energy. More general damping need not leave undamped mode shapes independent.

Write the perturbation Stokes streamfunction as $\psi =\varphi (r)e^{ikz-i\omega t}$, with $\omega =kc$ and nonzero $k$. Its radial and axial velocity amplitudes are $u=-ik\varphi /r$ and $w={\varphi }^{\prime }/r$. These automatically obey incompressibility. Linearizing the radial and axial Euler equations, with pressure amplitude $P$, givesThusDifferentiating the first expression, the terms involving $U^{\prime }{\varphi }^{\prime }$ cancel. Equating the two expressions for $P^{\prime }$ proves the axisymmetric inviscid pipe stability equation:At the axis, finite smooth axisymmetric velocity requires $\varphi =O(r^2)$, after choosing the irrelevant constant in the streamfunction. At the wall, impermeability gives $\varphi (a)=0$. There is no tangential no-slip condition in this inviscid problem. A decoupled azimuthal disturbance satisfies $ik(U-c)u_{\theta }=0$ and does not introduce a growing discrete mode when $U\ne c$.

For real $k$, suppose $c=c_r+i{c}_i$ with $c_i\ne 0$, so division by $U-c$ is legitimate. Multiply the ordinary differential equation by ${\varphi }^{\ast }/[r(U-c)]$, integrate, and apply integration by parts. The endpoint terms vanish by the preceding regularity and wall conditions, leavingTaking its imaginary part yieldsA growing temporal mode therefore requires a sign change of $(U^{\prime }/r{)}^{\prime }$, rather than a sign change of $U^{\prime \prime }$ alone. This is the cylindrical counterpart of Rayleigh's inflection-point theorem. If this derivative vanishes identically, the positive first integral rules out a nontrivial growing real-wavenumber mode as well. Isolated zeros without a sign change do not make the second integral vanish for a nonzero regular mode.

For the parabolic Poiseuille flow, $U^{\prime }/r=-2U_0/a^2$ is constant. Under the specified exclusion of critical levels, division givesSetting $\varphi =rf$, or equivalently choosing $n=1$ in $r^{-n}\varphi$, turns this intoIt is the Modified Bessel differential equation of order one. Hence $\varphi =r[A{I}_1(kr)+B{K}_1(kr)]$. Since $r{K}_1(kr)$ tends to a nonzero constant at the axis, it gives a singular radial velocity and must be excluded. The admissible mode shape is thereforeFor $k=0$, direct solution gives $\varphi =C{r}^2$ after axis regularity, and the wall condition again forces the trivial solution. For nonzero $k$, the nonzero roots of $I_1$ are imaginary. Writing $j_{1,l}>0$ for the zeros of the ordinary Bessel function $J_1$, the evanescent potential modes of inviscid Poiseuille flow areTheir perturbation vorticity vanishes: ${\partial }_{z}u-{\partial }_{r}w=k^2\varphi /r-({\varphi }^{\prime }/r{)}^{\prime }=0$. These modes grow exponentially in one axial direction and decay in the other. They can describe end-forced or spatially localized potential fields, but neither sign gives a bounded normal mode on a whole infinite pipe. Moreover, the reduced radial equation contains no $c$: it supplies no temporal eigenvalue condition. An arbitrary imaginary part assigned to $\omega =kc$ here is consequently not evidence for a physical temporal instability.

The real-wavenumber calculation does exclude growing discrete axisymmetric modes of this inviscid profile. It does not establish every form of stability: singular neutral disturbances with $c$ in the range of $U$, transient dynamics, non-axisymmetric disturbances, and finite-amplitude transition require separate treatment. The exclusion $U\ne c$ specifically removes the inviscid continuous spectrum.

With viscosity, linearize the Navier-Stokes equations instead and eliminate pressure to obtain a fourth-order cylindrical counterpart of the Orr-Sommerfeld equation. Impose both $\varphi (a)=0$ and ${\varphi }^{\prime }(a)=0$, with appropriate regularity at the axis. The resulting eigenvalue problem involves the Reynolds number and determines complex frequencies, rather than leaving $c$ arbitrary. Viscosity supplies diffusion, no-slip boundary layers and critical-layer regularization, so the inviscid sign argument alone no longer decides its spectrum. That observation does not imply that viscous pipe flow must have a growing linear eigenmode; finite-amplitude pipe transition is a different question.

The wave equation is a second-order linear partial differential equation describing waves with finite propagation speed. In one spatial dimension its homogeneous, constant-coefficient form isHere $x$ is position, $t$ is time, $u(x,t)$ is the disturbance, and subscripts denote partial derivatives. For a wave equation on a string, $u$ is transverse displacement; in linear acoustics, it can be a pressure perturbation. In several spatial dimensions, $u_{xx}$ becomes the Laplacian $\Delta u$. The same mathematical structure connects vibrating strings and acoustic waves, although their physical assumptions and wave speeds differ. ^[ref]

For the whole-line Cauchy problem, prescribe both displacement $u(x,0)=f(x)$ and velocity $u_t(x,0)=g(x)$. A classical solution has continuous partial derivatives through order two; $f\in C^2(\mathbb{R})$ and $g\in C^1(\mathbb{R})$ suffice for D'Alembert's formula to produce one. A finite string additionally needs boundary conditions. These statements concern the unforced model with constant $c$; damping, spatially varying coefficients and nonlinear terms change the solution formula.

A travelling wave $F(x-ct)$ carries its entire profile to the right at wave speed $c$; $G(x+ct)$ travels left. Their sum is again a solution by the superposition principle. A standing wave instead has a fixed mode shape whose size and sign oscillate: in $\sin (kx)\cos (ckt)$, each node of a standing wave stays stationary. The distinction is especially clear in the pulse and string diagrams below. For this ideal model, the dispersion relation is ${\omega }^2=c^2{k}^2$: each propagation branch is nondispersive, so a localized profile translates without spreading. This makes the equation a useful reference against which to understand wave dispersion and damping.

The appealing combination is local interaction with global predictability: neighbouring displacements determine acceleration through $u_{xx}$, yet D'Alembert's formula explicitly identifies exactly which initial locations can influence a later observation. Finite propagation speed turns that observation into a domain of dependence, rather than requiring knowledge of the entire initial string.