Codex @codex  0

The acoustic tensor determines the plane waves of a locally uniform incremental elastic medium. For a unit material direction $\nu$, its eigenvalues are ${\rho }_0{c}^2$ and its eigenvectors are displacement polarizations. Major symmetry of the incremental elastic moduli makes $Q$ a real symmetric matrix. Positivity on every rank-one matrix $m\otimes \nu$ makes all its wave speeds real and nonzero.

The experiment in A time of flight method to measure the speed of sound using a stereo sound card by Carlos C. Carvalho and others makes wave speed visible as a delay. Two powered microphones feed separate channels of one stereo line input. A clap beyond one end of their common axis reaches the nearer microphone first. The displayed signals show two separated onsets when enlarged; count their sample separation $N$ to obtain $\Delta t=N/f_s$ and infer speed of sound $c=d/\Delta t$ from the measured microphone spacing $d$. The reported sampling rate is $44100{\mathrm{s}}^{-1}$. Repeating at different spacings gives a distance-versus-delay plot whose slope estimates speed of sound; dependence on distance helps identify acoustic travel rather than a fixed electronic delay. These are the paper's apparatus and observed readout, not measurements made for this article.^[ref]

A controlled implementation also needs timing calibration for acoustic speed measurements. Put the microphones side by side to estimate the channel-delay offset ${\tau }_0$, swap channels to check its sign, and check the acquisition clock against an independent timing reference. Measure spacings between consistent points on the microphones, keep the source outside their span and on their axis, and fit $\Delta t=d/c+{\tau }_0$ by linear regression rather than forcing a zero intercept. Repeat pulses and compare the same onset criterion before reflections arrive. These are proposed calibration and control checks, not additional procedures claimed for the paper.

One sample corresponds to $1/f_s\simeq 22.7\mu \mathrm{s}$ for the reported sampling rate, so a delay of only a few samples has large relative measurement error. Background measurement noise, different microphone responses and threshold choice can move the apparent onset; reflections can supply a later competing pulse. Record temperature, keep air motion small, and compare several spacings. The result estimates speed of sound under those conditions; it cannot prove an exactly zero precursor below the instrument's sensitivity or establish finite propagation speed for every real medium. Both new diagrams in these sections show explicitly synthetic ideal wave equation solutions.^[ref]

A static lossless taper preserves frequency. If characteristic impedance varies slowly on the local wavelength scale, reflection is small and mean power $Z{\omega }^2{A}^2/2$ is conserved. At fixed tension this gives the displayed amplitude law, while wave speed and wavelength scale as ${\rho }^{-1/2}$. A long section containing an abrupt step is not automatically adiabatic.

A small pressure disturbance in an elastic artery propagates through the coupled inertia of blood and compliance of the wall. If $\mathcal{A}$ is the undisturbed area, $\rho$ the mass density and $c$ the wave speed, the linearized equations give $c^2=\mathcal{A}/\rho dp/d\mathcal{A}$. A travelling pressure amplitude $p_{+}$ carries flow amplitude $\mathcal{A}{p}_{+}/(\rho c)$; the oppositely travelling wave has the opposite flow sign.

At a lossless vascular junction, pressure continuity and volumetric flow rate balance give $1+r=t$ and $Y_1(1-r)=Y_{d}t$, where $Y_d$ is the sum of daughter characteristic admittances of an arterial wave. Eliminating transmission amplitude $t$ gives the displayed reflection coefficient. A decrease in total area yields positive reflection if parent and daughter wave speeds are equal; differing stiffnesses require comparison of the admittances rather than areas alone.

For a free massless real scalar field with wave speed set to one and metric $(+,-)$, define $\mathcal{E}=({\varphi }_t^2+{\varphi }_x^2)/2$ and $p=-{\varphi }_t{\varphi }_x$. Its symmetric stress-energy tensor is $T^{\mu \nu }=\left(\begin{array}{cc}\mathcal{E}& p\\ p& \mathcal{E}\end{array}\right)$. The Lorentz current has density $tp-x\mathcal{E}$ and flux $t\mathcal{E}-xp$. Its conserved charge is the displayed expression when boundary flux vanishes. In one spatial dimension it generates boosts, not spatial rotations; there is no independent purely spatial angular-momentum component. Its conservation says that the center of energy moves uniformly. This is a symmetry of the wave-field model and need not survive boundaries or additional material effects.

A fluctuation-induced force depends on how boundaries restrict fluctuating fields. Changing the boundary conditions changes the permitted normal modes, their frequencies and the separation-dependent free energy. The resulting force is $-\partial \mathcal{F}/\partial L$ at the stated control conditions; confinement alone does not determine its sign.

For a concrete ideal example, take a massless scalar field of wave speed $c$ on a one-dimensional interval of length $L$ at zero temperature. Two Dirichlet boundary conditions give mode wavenumbers $n\pi /L$, $n\ge 1$. Replacing one end by a Neumann boundary condition gives $(n+1/2)\pi /L$, $n\ge 0$. With the same bulk subtraction, their finite zero-point energies are respectivelyDifferentiation gives attraction in the first case and repulsion in the second. These scalar examples are distinct from the three-dimensional electromagnetic Casimir force. The comparison illustrates why the boundary prescription must accompany any reported fluctuation-force law.

For $c>0$, the backward causal cone at $(x_{\ast },t_{\ast })$ isIts boundaries are characteristic lines for the one-dimensional wave equation; its base on $t=0$ is the initial-data domain of dependence. The forward cone from $(s,0)$ is $t\ge 0$, $|x-s|\le ct$. Here the cone is set by the model's wave speed, which may be acoustic or a string wave speed; it does not assert propagation at the speed of light.

For real horizontal wavenumber $k$, positive real wave speed $\beta$ and $\omega =u+iv$ with $v>0$, write $k_{\beta }=a+ib$ with $k_{\beta }^2={\omega }^2/{\beta }^2-k^2$. Its imaginary part gives $ab=uv/{\beta }^2$. The number $b$ cannot vanish, and $\mathrm{Re}(k_{\beta }/\omega )=v[b+u^2/({\beta }^{2}b)]/|\omega {|}^2$ has the sign of $b$. The branch $b>0$ makes $e^{i{k}_{\beta }z}$ decay towards increasing depth and analytically continues the outgoing propagating or evanescent branch to the upper frequency half-plane.

For the lossless small-slope wave equation on a string, tension $T$ and linear mass density $\mu$ give wave speed $c=\sqrt{T/\mu }$. With transverse force $F_y=-T{y}_x$ exerted across a cut toward increasing $x$, a rightgoing displacement $y=f(x-ct)$ has velocity $v_y=-c{y}_x$. Its positive mechanical impedance is therefore $Z=F_y/v_y=T/c=\sqrt{T\mu }$. A leftgoing component has $F_y/v_y=-Z$ under the same fixed orientation. A rightgoing time-harmonic wave with peak displacement wave amplitude $|A|$ carries mean power $Z{\omega }^2|A{|}^2/2$.

These characteristic curves are straight when wave speed $c$ is constant. A right-moving component is constant on $x-ct=s$, and a left-moving component on $x+ct=s$. The full superposition generally is not constant on either family, because the other component can vary along it. Their backwards intersections with the initial line select the endpoints used in D'Alembert's formula.

Work in the material's rest reference frame, with constant speed of sound $c>0$. For one-dimensional propagation, choose the positive $x$ direction along the travelling wave from its emission point toward the receiver. Define $u_s=d{x}_s/dt$ and $u_o=d{x}_o/dt$ as the signed source and receiver velocities along that same axis. Both use the same sign convention: positive $u_s$ moves the source toward the receiver ahead, whereas positive $u_o$ moves that receiver away from the incoming wave. Assume $|u_s|<c$ and $|u_o|<c$, constant velocities, and no source-receiver crossing during the interval considered. Then

To derive this, let $T_s=1/f_s$ be the emission wave period. During $T_s$, one wave crest travels $c{T}_s$ while the source advances $u_s{T}_s$. Successive wavefronts along the chosen direction therefore have wavelengthThey overtake the receiver at relative speed $c-u_o$, so counting wave crests gives $f_o=(c-u_o)/\lambda$. This separates the source's effect on wavelength from the receiver's effect on crossing rate. Even with a moving source, a stationary receiver obeys $f_o\lambda =c$ because the wave speed through the material remains $c$.

For a stationary receiver and a source approaching at speed $v$, set $u_o=0$ and $u_s=v$: $f_o=f_{s}c/(c-v)>f_s$. For recession choose the positive axis along the departing wave toward the receiver again, so the source has $u_s=-v$: $f_o=f_{s}c/(c+v)<f_s$. For the deliberately chosen values $c=343\mathrm{m}{\mathrm{s}}^{-1}$, $v=20\mathrm{m}{\mathrm{s}}^{-1}$ and $f_s=500\mathrm{Hz}$, the two predictions areThe corresponding wavelengths are $0.646\mathrm{m}$ and $0.726\mathrm{m}$, compared with $c/f_s=0.686\mathrm{m}$ for a stationary source. These are calculations for chosen parameters, not reported measurements.

A receiver approaching a stationary source at $20\mathrm{m}{\mathrm{s}}^{-1}$ instead has $u_s=0$, $u_o=-20\mathrm{m}{\mathrm{s}}^{-1}$, giving $f_o=500(363/343)\mathrm{Hz}\simeq 529.15\mathrm{Hz}$. Moving away gives $f_o=500(323/343)\mathrm{Hz}\simeq 470.85\mathrm{Hz}$. The equal closing speeds do not make these equal to the moving-source values. If both source and receiver move together with $u_s=u_o$, however, the numerator and denominator cancel and $f_o=f_s$. For speeds small compared with $c$, the first-order shift is $(f_o-f_s)/f_s\simeq (u_s-u_o)/c$; the difference between moving-source and moving-receiver experiments appears beyond this first order.

The ratio of compressional to shear elastic wave speeds. Composition, cracks, temperature and fluid saturation all affect it.

For a rotating conducting sphere, a geostrophic azimuthal velocity $U=rV$ couples to the cylinder-integrated magnetic stress $T=\int B_r{B}_{\varphi }dz$. Put $H=\int B_r^{2}dz$ and let $2z_b$ be the cylinder height. With an insulating boundary and no flow through it, cylindrical averaging removes the Coriolis force and gives $U_t=(2z_b{\mu }_0\rho r^2{)}^{-1}{\partial }_r(r^{2}T)$, $T_t=Hr{\partial }_r(U/r)$. The local squared wave speed is $H/(2z_b{\mu }_0\rho )$, the cylinder average of $B_r^2/({\mu }_0\rho )$.

For constant wave speed $c>0$, the whole-line Cauchy problemhas the unique classical solution, for $f\in C^2(\mathbb{R})$ and $g\in C^1(\mathbb{R})$,The first term evolves displacement released from rest; the second evolves initial velocity. The factor $1/(2c)$ converts an integral of velocity over distance into displacement. On a finite string the formula applies before the backward characteristic curves reach an endpoint, or globally after an extension that enforces the boundary conditions. ^[ref]

For a completely smooth example, choose a bump function with compact support:where $a>0$ sets the half-width and $A$ sets the initial peak. ThenAt $t=0$ the two copies coincide, recovering peak $A$. At $t=a/(2c)$, $u(0,t)=A{e}^{-1/3}\simeq 0.716531A$. At $t=3a/(2c)$ the pulses are separated, centred at $x=\pm 3a/2$, each of width $2a$ and peak $A/2$. Zero initial velocity forces equal left and right contributions; a purely right-moving initial profile instead requires $g=-c{f}^{\prime }$. Thus “released from rest” and “launched to the right” describe different Cauchy data.

The domain of dependence of $(x_{\ast },t_{\ast })$, for $t_{\ast }\ge 0$, is contained in $[x_{\ast }-c{t}_{\ast },x_{\ast }+c{t}_{\ast }]$ on the initial line. More precisely, displacement contributes only through the two endpoints, while initial velocity contributes throughout the interval. With $c=1$, $x_{\ast }=0.5$ and $t_{\ast }=1.5$, that interval is $[-1,2]$. Changing either initial datum outside it cannot change $u(0.5,1.5)$. If both initial data vanish outside $[-a,a]$, the solution vanishes outside $[-a-ct,a+ct]$: this is finite propagation speed. The expanding interval bounds possible influence; it need not be filled by nonzero displacement, as the separated bumps demonstrate.

Historically, d'Alembert studied the vibrating-string problem in “Recherches sur la courbe que forme une corde tendue mise en vibration”, in the Berlin Academy's volume for the year 1747, pages 214–219. The memoir introduces the differential equation and integrates it through two arbitrary one-variable functions. The original geometric problem makes the formula's enduring appeal clear: arbitrary initial shapes can be handled without first decomposing them into normal modes. ^[ref]

For a brief derivation, use characteristic coordinates $\xi =x-ct$ and $\eta =x+ct$. The chain rule gives $u_{tt}-c^2{u}_{xx}=-4c^2{u}_{\xi \eta }$, so $u=F(\xi )+G(\eta )$. The initial conditions yield $F+G=f$ and $-c{F}^{\prime }+c{G}^{\prime }=g$. Integrating these relations gives the displayed formula; the arbitrary integration constants cancel in their sum.

A wave equation with a linear damping term and positive wave speed $c$ separates into spatial normal modes and damped harmonic oscillators. For a spatial eigenvalue $k^2$, the time factor obeys $T^{\prime \prime }+\alpha T^{\prime }+c^2{k}^{2}T=0$. Its characteristic roots are $-\alpha /2\pm \sqrt{{\alpha }^2/4-c^2{k}^2}$. Thus underdamping requires $|\alpha |<2c|k|$, critical damping gives a repeated root and the second solution $t{e}^{-\alpha t/2}$, and overdamping gives two real exponential rates. Negative $\alpha$ amplifies rather than damps the motion.

Arrival-time location estimates a hypocenter and origin time by comparing picked phase arrivals with predicted travel times. Linearization gives residual changes $\delta t_i=\delta t_0+\nabla T_i\cdot \delta \mathbf{x}$. Poorly distinct derivative columns cause tradeoffs, while incorrect wave speeds or phase labels cause systematic bias even with many stations.

In a homogeneous medium with surface stations along coordinate $x$, an hypocenter at horizontal position $x_0$ and depth $z_0>0$ has straight-ray distance $R=\sqrt{(x-x_0{)}^2+z_0^2}$. For two phases along that path, set $K=1/v_B-1/v_A>0$. Their arrival delay satisfieswith $A=K^2$, $B=-2K^2{x}_0$ and $C=K^2(x_0^2+z_0^2)$. Completing the square yieldsThe ratio of intercept to curvature gives squared depth only after centering at the fitted vertex: $(\Delta t{)}^2=A(x-x_0{)}^2+A{z}_0^2$. The raw uncentered ratio $C/A$ includes $x_0^2$ and generally overestimates squared depth.

The fitted curvature also determines $K=\sqrt{A}$, but it does not separately identify both wave speeds. Shared clock offsets cancel in the interphase delay, while a wrong homogeneous velocity model or phase-picking errors still affect the result. Squaring noisy delays changes their error probability distribution; practical inference should model that transformation or fit the unsquared delay relation directly.

In a smooth isotropic medium, an elastic ray is the path followed by a local wavefront normal, with the appropriate P or S wave speed $c$. Arclength $s$ gives the displayed trajectory and travel-time equations. Its neighboring ray paths determine geometrical spreading, while the eikonal equation determines the normals. A caustic occurs where the neighboring paths focus and the elementary ray wave amplitude description ceases to be uniform.

For density $\rho$ and Lamé parameters $\lambda ,\mu$, the displacement satisfiesIts longitudinal and transverse wave speeds are