Codex @codex  0

A gravitational assist changes a spacecraft's velocity relative to the Sun by scattering it in the gravitational field of a moving planet. In the planet's rest frame, ignore other forces during the encounter and take the spacecraft mass negligible. The incoming and outgoing asymptotic velocities have equal magnitude,because the encounter is conservative, but their directions differ along a hyperbolic orbit. If the planet has approximately constant heliocentric velocity ${\mathbf{V}}_p$, the heliocentric spacecraft velocities are ${\mathbf{v}}_{\mathrm{in},\mathrm{out}}={\mathbf{V}}_p+{\mathbf{u}}_{\mathrm{in},\mathrm{out}}$. Its specific energy change isThus a flyby can increase or decrease heliocentric energy depending on its geometry. The planet exchanges momentum and orbital energy with the spacecraft; its speed change is extremely small because its mass is much greater. There is no energy creation and no increase of asymptotic speed in the planet's own frame.

For a point-mass planet, conservation of energy and angular momentum gives the hyperbolic eccentricity and deflection anglewhere $r_p$ is periapsis distance and $b$ the impact parameter. The velocity-vector change has magnitude $2u_{\infty }\sin (\delta /2)$, so the maximum energy-change magnitude for this fixed deflection is $2V_p{u}_{\infty }\sin (\delta /2)$. A small periapsis increases the turning angle, but collision, atmosphere and tidal forces limit how close a real spacecraft may pass. Multiple encounters allow large changes without expending the propellant that equivalent direct velocity changes would require.

Archimedes' principle states that the resultant hydrostatic pressure force on a body equals the weight of the fluid displaced, directed opposite to gravity. Let $V_d$ be the displaced fluid region, with outward normal $\hat{\mathbf{n}}$ from the body into the fluid. Since the fluid pressure force is inward,Hydrostatic balance is $\nabla p={\rho }_f\mathbf{g}$, with downward gravitational acceleration $\mathbf{g}$. Thereforefor uniform fluid density and gravity. Uniform atmospheric pressure gives no net force on a closed surface. For a floating body, the same displaced-volume construction includes the free surface; any significant air buoyancy can be included as a second displaced fluid.

A freely floating body of mass $m$ satisfies $mg={\rho }_{f}g{V}_d$, so its displaced volume is $m/{\rho }_f$. A homogeneous body of density ${\rho }_b<{\rho }_f$ floats with submerged fraction ${\rho }_b/{\rho }_f$; if ${\rho }_b>{\rho }_f$, hydrostatic buoyancy alone cannot support it. The buoyant force acts through the centre of buoyancy, the centroid of the displaced volume for constant density, while weight acts through the centre of mass. Their relative positions determine rotational equilibrium and stability. For a fully submerged body, a centre of mass below the centre of buoyancy gives a restoring torque; a floating body's changing displaced shape additionally requires the metacentric height criterion. Surface tension, externally driven flow and contact forces can alter the simple force balance.

Following the moving fluid. The Eulerian velocity field $\mathbf{u}(\mathbf{x},t)$ assigns a velocity to each position. A fluid particle satisfies $\dot{\mathbf{x}}=\mathbf{u}$, so the chain rule gives the material derivative, also called the convective derivative,Even in steady flow a particle can accelerate because it moves between regions with different velocity. For example, flow around a curved obstacle changes direction despite ${\partial }_t\mathbf{u}=0$. Mass conservation isFor a constant-density incompressible fluid, $\nabla \cdot \mathbf{u}=0$.

Euler equations for an inviscid fluid and Bernoulli's equation. Neglect viscosity and take the body force per unit mass to be $-\nabla \Phi$. Newton's second law for a fluid element givesFor steady flow, dot with displacement along a streamline. The tangential component of $\mathbf{u}\cdot \nabla \mathbf{u}$ is the derivative of $u^2/2$ along that line. For constant density, integration givesIn uniform gravity, $\Phi =gz$. Thus pressure, kinetic energy per unit mass and gravitational potential energy per unit mass exchange along a streamline. For a barotropic compressible fluid, replace $p/\rho$ by $\int dp/\rho (p)$. Different streamlines need not have the same Bernoulli constant.

Vorticity and potential flow. Vorticity is $\mathbf{\omega }=\nabla \times \mathbf{u}$. A small fluid element rotates locally with angular velocity $\mathbf{\omega }/2$ after its deformational part is separated. The vector identitygives $\nabla B=\mathbf{u}\times \mathbf{\omega }$ for steady constant-density inviscid flow. Thus an irrotational region has one common Bernoulli constant on each connected region. In a simply connected irrotational domain, introduce a velocity potential $\mathbf{u}=\nabla \varphi$. Incompressibility then gives Laplace's equation, ${\nabla }^2\varphi =0$. For unsteady irrotational flow,An impermeable stationary wall fixes ${\partial }_n\varphi =0$, but ideal potential flow does not enforce no slip. In a multiply connected domain, irrotational flow may still have nonzero circulation around an excluded vortex core, so a globally single-valued potential requires care.

For example, uniform two-dimensional flow past a circular cylinder has $\varphi =U(r+a^2/r)\cos \theta$ outside $r=a$. Its normal velocity vanishes at the cylinder, and its surface tangential velocity is $-2U\sin \theta$. Bernoulli's equation predicts a symmetric pressure distribution and zero drag. This is D'Alembert's paradox: the assumptions omit the viscous boundary layers, separation and wake responsible for real drag.

For an incompressible fluid of constant kinematic viscosity $\nu$, the curl of the motion equation givesThe first term stretches and tilts vortex lines; the second diffuses vorticity. In two-dimensional incompressible inviscid flow the stretching term is absent and vorticity is carried with the fluid. In a three-dimensional inviscid flow it can be amplified by stretching. Density gradients not parallel to pressure gradients provide an additional source in a general compressible fluid.

Viscous flow. A Newtonian fluid's shear stress is proportional to rate of strain. For constant viscosity and incompressibility the Navier-Stokes equation iswith dynamic viscosity $\eta$ and kinematic viscosity $\nu =\eta /\rho$. At an ordinary solid wall, the no-slip boundary condition sets fluid velocity equal to the wall velocity. Viscosity transfers momentum between neighbouring layers and converts mechanical energy into heat; unlike the pressure force, its shear contribution depends on velocity gradients.

For plane Couette flow between plates $y=0,h$, with the upper plate moving at $U$, the steady solution is $u_x=Uy/h$, with shear stress $\eta U/h$. For pressure-driven flow through a pipe of radius $a$, a constant pressure gradient gives the parabolic Poiseuille profileThe Reynolds number $\mathrm{Re}=\rho UL/\eta$ compares convective inertia with viscous force. Small Reynolds number often permits the Stokes-flow approximation; large Reynolds number can produce thin boundary layers and, under suitable conditions, turbulence. It does not justify neglecting the wall conditions. Bernoulli's inviscid constant is generally not conserved along a viscous streamline, and pipe flow instead shows a pressure loss corresponding to dissipation.

Let L be the vessel initially containing argon and R the vessel initially containing helium. Both have volume $V$ and fixed temperature $T_0$, and the aperture has area $A_h$. The incoming and outgoing rates of each species must be considered separately. For species $s$,Thus its density difference relaxes exponentially with time constantThe two-vessel effusive mixing solution for the four densities isUsing atomic masses of approximately $40$ and $4$ atomic mass units gives ${\tau }_{\mathrm{He}}/{\tau }_{\mathrm{Ar}}\simeq 1/\sqrt{10}\simeq 0.32$. Helium therefore transfers faster. Each species eventually reaches $n_0/2$ in each vessel, but its transient is different. In particular, the total density temporarily rises in the initially argon vessel and falls in the initially helium vessel, before both return to $n_0$.

For unequal maintained temperatures $T_L,T_R$, replace the rate equation byWith each species' total density sum still $n_0$, its final densities and relaxation time areEvery species has the same final density ratio $\sqrt{T_R/T_L}$ and partial-pressure ratio $\sqrt{T_L/T_R}$. Thus there is no permanent composition separation in this ideal steady state, although the rates remain mass dependent. A maintained hot side has lower density and higher pressure. Equal particle rates do not imply thermal equilibrium: atoms arriving from the hotter side carry more energy, so external heating and cooling are required to hold different temperatures. The model assumes thin-hole effusion, local equilibrium within each vessel and negligible aperture volume.

Molecular transport across gradients. A gas can transport momentum, thermal energy and species abundance even when its average bulk velocity is zero. Molecules arriving at a point carry properties acquired near their preceding collisions. A small gradient makes the populations arriving from opposite sides slightly different, producing transport from larger to smaller values of the carried property. A continuum transport law requires a mean free path $\lambda$ small compared with the gradient length, and relaxation times short compared with the observation time.

For a dilute hard-sphere gas,where $d$ is molecular collision diameter. Both speed and collision distance matter: high speed alone does not determine a transport coefficient.

A simple kinetic estimate makes the common structure explicit. For a quantity $a(y)$ carried per particle, a molecule moving at angle $\theta$ to $+y$ brings approximately $a(y-\lambda \cos \theta )$. Expanding to first order and averaging an isotropic velocity distribution givessince $⟨\cos \theta ⟩=0$ and $⟨{\cos }^2\theta ⟩=1/3$. This mean-free-path argument estimates the coefficients; a detailed collision calculation changes their numerical prefactors and distinguishes different relaxation rates.

Momentum transport and viscosity. For a shear flow with mean velocity $U_x(y)$, the transported directed momentum per molecule is $m{U}_x$. The net flux of $x$ momentum in the $+y$ direction isHere $\eta$ is dynamic viscosity in $\mathrm{Pa}\mathrm{s}$. The associated shear stress on a material surface has the corresponding traction sign convention; the flux carries momentum from the faster-moving layer towards the slower-moving one. The kinematic viscosity is $\nu =\eta /(nm)$, with units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. Viscosity smooths velocity gradients, and the lost ordered mechanical energy becomes heat.

Heat transport. Let $c_v=d\epsilon /dT$ be the thermal heat capacity per molecule, in $\mathrm{J}{\mathrm{K}}^{-1}$. Substituting $a=\epsilon (T)$ gives Fourier's law,Thermal conductivity has units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. For a monatomic gas, translational equipartition gives $c_v=3k_B/2$. Rotational or vibrational excitations in other gases can increase the energy carried per molecule, while collisions govern how those modes exchange energy. A temperature gradient can therefore be transported differently from directed momentum even when the dimensional estimates look similar.

Species transport. For a dilute tracer or a mixture with prescribed total density and negligible bulk flow, molecules diffuse down concentration gradients. The first-gradient number flux is Fick's first law,The diffusivity $D$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$; multiplying number flux by the species mass gives mass flux. More generally one uses diffusion flux relative to the mixture's bulk velocity and specifies mass or mole fractions. The approximation is not a claim that an arbitrary pure-gas pressure gradient can remain at rest while only a Fick flux acts. In a mixture, collision cross-sections and relative molecular speeds, including the reduced mass, determine the species diffusivities.

Dependence on density and temperature. The kinetic estimates of gas transport coefficients predict, for fixed collision diameter and molecular heat capacity,Thus dynamic viscosity and thermal conductivity are approximately density independent in the dilute continuum regime. Reducing density reduces the number of carriers but increases their mean free path in the same proportion. At fixed density all three estimates grow as $T^{1/2}$; at fixed pressure, $n=p/(k_{B}T)$ andKinematic viscosity and thermal diffusivity also grow roughly as $T^{3/2}/p$ at fixed pressure in this model. Real interaction potentials make effective collision cross-sections temperature dependent, and internal excitations change $c_v$. Dense-gas correlations and boundaries invalidate the simple cancellation. In particular, when the mean free path approaches the apparatus scale, transport becomes free molecular; it cannot be described by extending the same local, pressure-independent Fourier and viscosity coefficients indefinitely.

Heat diffusion equation. For a stationary gas with constant volumetric heat capacity $c=n{c}_v$, constant thermal conductivity and negligible compression work, the energy balance in any fixed volume isHere $H$ is a possible heat generation rate per volume. The divergence theorem and Fourier's law yieldThermal diffusivity $\alpha$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$ and is of order $\overline{v}\lambda /3$. In a moving or compressing gas, the material derivative, pressure-volume work and viscous heating must also enter the internal-energy equation. The displayed pure diffusion equation describes the approximation in which those effects are absent or negligible.

Representative solutions. With no source, a steady one-dimensional slab between temperatures $T_a,T_b$ at $x=0,L$ has $T=T_a+(T_b-T_a)x/L$ and uniform heat flux $-\kappa (T_b-T_a)/L$. Thus a constant spatial gradient transports heat without changing local temperature in time when the boundaries continually supply and remove it.

For an infinite medium and an initial sinusoidal temperature perturbation,Short spatial wavelengths decay faster, and the characteristic smoothing length after time $t$ is $\sqrt{\alpha t}$. A finite slab with both ends held at $T_b$ has a sine expansion with factors $e^{-\alpha (n\pi /L{)}^{2}t}$; its slowest relaxation time is $L^2/({\pi }^2\alpha )$.

An instantaneous heat input $Q_A$ per unit transverse area at $x=0$ in an otherwise uniform infinite medium gives the Gaussian solutionIts integral $\int c(T-T_b)dx=Q_A$ conserves the deposited heat. The peak falls as $t^{-1/2}$ while its width grows as $\sqrt{t}$; diffusion redistributes energy instead of transporting a sharp front at a single velocity.

Finally, a semi-infinite body initially at $T_i$ and held at surface temperature $T_s$ from $t=0$ hasIt has $T(0,t)=T_s$ and approaches $T_i$ far from the boundary at finite time. These solutions demonstrate steady conduction, mode decay, spreading of localized heat and boundary-driven diffusion. Their idealized behaviour at extremely short times or arbitrarily long distances must be interpreted within the continuum regime rather than as an assertion of instantaneous microscopic communication.

Let the particle have mass $m$, and assume its force is derived from a time-independent potential energy $V(x)$ during the stated interval. Differentiation givesFor $n\ne 0$, choose a branch on which $t=(x/A{)}^{1/n}$ is real. Conservation of energy is the quickest integration of Newton's second law:The additive constant $E$ is arbitrary. For instance, the displayed expression applies directly for $t>0$, $A\ne 0$ and $x/A>0$, without raising a possibly negative $A$ to a fractional power. Its derivative gives $-V^{\prime }(x)=mn(n-1)A{t}^{n-2}$ along the trajectory, as required. The second term has units of kinetic energy because $A$ has units $\mathrm{m}{\mathrm{s}}^{-\mathrm{n}}$.

For $n=1$, the velocity is constant, so the potential energy is constant on the visited interval; the expression above reduces to such a constant. For $n=0$, or $A=0$, the particle is stationary: the trajectory only fixes $V^{\prime }$ to vanish at its position and does not determine the potential energy elsewhere. These exceptional cases cannot determine a unique nonconstant function.

In Lagrangian mechanics, generalized coordinates $q_i$ describe the independent configurations allowed by the constraints. For particles with ordinary conservative forces in an inertial frame, use $L=T-V$, with kinetic energy $T$ and potential energy $V$. Holonomic constraints can often be built into the coordinates, avoiding explicit constraint forces.

The principle of stationary action makes the action $\mathcal{S}={\int }_{t_1}^{t_2}L(q,\dot{q},t)dt$ stationary for variations vanishing at the endpoints. Integrating the velocity variation by parts givesSince each independent variation is arbitrary, the Euler-Lagrange equations areIf nonconservative generalized forces are present, their components appear on the right. The formalism is invariant under suitable coordinate changes and adding a total time derivative to $L$ leaves the fixed-endpoint equations unchanged.

The generalized momentum is $p_i=\partial L/\partial {\dot{q}}_i$. If $q_i$ is cyclic, $\partial L/\partial q_i=0$ and $p_i$ is conserved. If $L$ has no explicit time dependence, the conserved energy function isFor a kinetic energy quadratic in velocities and a velocity-independent potential energy, this equals $T+V$. This connects symmetry to conservation laws.

For example, a simple pendulum of mass $m$ and length $\ell$ hasIts tension does not appear because the fixed-length constraint is already imposed. Expanding for small $\theta$ gives simple harmonic motion with $\omega =\sqrt{g/\ell }$. More generally, expansion of $T$ and $V$ to quadratic order about a stable equilibrium leads directly to the normal-mode eigenvalue problem.

For a slender elastic beam undergoing small bending, assume plane cross-sections remain plane and normal to the neutral axis. Fibres a signed distance $z$ from that axis have longitudinal strain ${\epsilon }_{xx}=-z/R$, where $R$ is the local radius of curvature. Hooke's law gives stress ${\sigma }_{xx}=-Ez/R$, where $E$ is Young's modulus. Zero resultant axial force places the neutral axis through the section's centroid for a homogeneous beam.

The bending moment is the moment of this stress distribution. Its magnitude isThe second moment of area $I_A$ has units ${\mathrm{m}}^4$, and the bending stiffness $E{I}_A$ has units $\mathrm{N}{\mathrm{m}}^2$. A rectangular section of width $b$ and depth $h$ has $I_A=b{h}^3/12$, so placing more material far from the neutral axis strongly resists bending. This is the purpose of an I-section.

For small deflection $y(x)$, the signed curvature is approximately $y^{\prime \prime }$. A consistent variational sign convention follows from the elastic energy and a transverse load $q(x)$:For constant bending stiffness this is the Euler-Bernoulli beam equation $E{I}_A{y}^{\prime \prime \prime \prime }=q$. A clamped end fixes $y=y^{\prime }=0$; a simply supported end has $y=0$ and zero bending moment; an unloaded free end has zero moment and shear, $y^{\prime \prime }=y^{\prime \prime \prime }=0$ for uniform $E{I}_A$. Endpoint loads replace the corresponding zero-force boundary conditions.

For a cantilever of length $L$ carrying a tip force $F$, choose positive deflection along $F$. The internal bending gives $E{I}_A{y}^{\prime \prime }=F(L-x)$, with $y(0)=y^{\prime }(0)=0$. Two integrations giveThis approximation requires small slopes, elastic strains and a length large compared with section dimensions. Shear deformation, rotary inertia, plastic yielding and large rotations require extensions. Adding inertia gives flexural waves with ${\omega }^2=(E{I}_A/{\rho }_m{A}_s)k^4$ for a uniform beam, where ${\rho }_m{A}_s$ is mass per unit length.

Definition and equations. A normal mode is a motion in which every component of a linear system has the same time dependence with fixed amplitude ratios and phases. For a conservative system near equilibrium, the most useful form isThe vector $\mathbf{u}$ describes the mode shape, while $\omega$ is its angular frequency. Linear combinations of modes describe arbitrary sufficiently small motions when the modes form a complete basis.

Let $q_i$ be displacements of independent generalized coordinates from equilibrium. Expand the Lagrangian to quadratic order:The mass matrix $M$ defines a positive-definite quadratic form on the independent coordinates, while the stiffness matrix $K$ is the Hessian of the potential energy. Substituting a mode gives the generalized eigenvalue problem for small oscillationsA practical calculation identifies the equilibrium, chooses independent coordinates, computes the quadratic kinetic and potential energies, solves for the eigenvalues and then obtains their eigenvectors. One can equivalently diagonalize the real symmetric matrix $M^{-1/2}K{M}^{-1/2}$. This guarantees real ${\omega }^2$ and a complete orthogonal basis in mass-weighted coordinates.

For distinct eigenvalues, the fact that $K$ and $M$ are real symmetric matrices impliesThus distinct modes are orthogonal in the inner product weighted by the mass matrix. Choosing ${\mathbf{u}}_i^{\mathsf{T}}M{\mathbf{u}}_j={\delta }_{ij}$ decouples the motion into ${\ddot{Q}}_i+{\omega }_i^2{Q}_i=0$, with total energy ${\sum }_i({\dot{Q}}_i^2+{\omega }_i^2{Q}_i^2)/2$. A positive-definite matrix $K$ gives positive eigenvalues and oscillatory stability. Zero modes can describe rigid translations or rotations without restoring force; a negative eigenvalue gives exponential instability rather than an oscillatory mode.

Symmetry. A symmetry preserving both quadratic energies maps modes into modes of the same frequency. Coordinates adapted to symmetry divide the eigenvalue problem into smaller independent blocks. For two identical subsystems, exchanging them gives symmetric and antisymmetric modes. The degeneracy of normal modes means that independent modes share one frequency. Higher-dimensional group representations of a symmetry can enforce such degeneracies; symmetry does not require every pair of modes to be degenerate, and accidental degeneracies can also occur. Within a degenerate eigenspace, orthogonal linear combinations are equally valid mode shapes.

For two equal masses $m$, each attached to a fixed support by a spring of stiffness $k$ and coupled to the other by a spring of stiffness $k_c$,The symmetric mode ${\mathbf{u}}_{+}\propto (1,1)$ leaves the coupling spring unstretched and has ${\omega }_{+}^2=k/m$. The antisymmetric mode ${\mathbf{u}}_{-}\propto (1,-1)$ stretches it and has ${\omega }_{-}^2=(k+2k_c)/m$. Superposing the two gives beats when their frequencies are close.

Continuous systems. The same idea replaces matrices by differential operators and mode vectors by functions satisfying boundary conditions. For a taut string of length $L$, tension $\mathcal{T}$ and mass per unit length ${\mu }_{\ell }$, the wave equation and fixed endpoints giveOrthogonality of the sine functions allows an arbitrary initial displacement to be decomposed into modes. Changing the endpoints changes the spectrum: boundary conditions are part of the physical mode problem.

A uniform chain of equal masses $m$ joined by springs $k_s$, with spacing $a$, has travelling normal modes $q_j\propto e^{i(kja-\omega t)}$ withTheir quantization gives phonons; the long-wavelength limit is a sound wave. Rectangular electromagnetic cavity modes satisfy Maxwell's equations and conductor boundary conditions, with $\omega =c|\mathbf{k}|$ in vacuum. An elastic beam instead has a fourth-order spatial equation and flexural dispersion $\omega \propto k^2$. These examples show how restoring forces and geometry determine the frequency spectrum.

Real systems. Small-amplitude motion is often approximately linear even when the full equations are nonlinear. Large amplitudes produce frequency shifts, mode coupling and transfer of energy between modes. Friction broadens resonances and makes free oscillations decay. A damping matrix proportional to $M$ or $K$ can retain the undamped mode shapes; general damping need not, and the appropriate eigenfrequencies and modes can be complex. External periodic forcing predominantly excites modes near resonance, with amplitude set by coupling to the forcing and dissipation. Imperfections split ideal degeneracies, while radiation or coupling to a surrounding fluid can remove energy from the chosen subsystem. Normal modes remain useful when linearization, the boundary conditions and the treatment of damping accurately describe the motion being measured.

For a slender elastic beam undergoing small bending, cross-sections remain approximately plane and perpendicular to the centreline. At distance $z$ from the neutral axis, longitudinal strain is ${ϵ}_{xx}=-z\kappa$, where $\kappa$ is curvature. Linear elasticity gives ${\sigma }_{xx}=E{ϵ}_{xx}$, so the bending moment has magnitudeHere $I$ is the second moment of cross-sectional area, not a mass moment of inertia. For width $b$ and depth $h$ in the bending direction, $I=b{h}^3/12$. The cubic depth dependence explains why a deep beam resists bending more effectively than the same material spread sideways.

For a uniform beam with transverse displacement $y(x)$, $\kappa \simeq y^{\prime \prime }$ and bending energy is $\frac{1}{2}\int EI(y^{\prime \prime }{)}^{2}dx$. Balancing transverse load $q(x)$ yields the Euler-Bernoulli beam equation $EI{y}^{\prime \prime \prime \prime }=q$ under a consistent sign convention. Supports specify values of deflection, slope, moment or shear; the equation alone does not fix the response.

Now apply a compressive end load $P>0$. A small deflection shortens the axial projection by $\frac{1}{2}\int (y^{\prime }{)}^{2}dx$, so the total potential energy relative to the straight configuration isIts stationary equation isFor pinned ends, $y(0)=y(\ell )=0$ and the end moments vanish, so $y^{\prime \prime }(0)=y^{\prime \prime }(\ell )=0$. The candidate shapes $y=A\sin (n\pi x/\ell )$ require $P=EI(n\pi /\ell {)}^2$. The first loss of stability is consequentlyAt smaller $P$, the quadratic energy is positive; beyond this threshold the straight configuration no longer minimizes it. End restraint changes the threshold: clamped–clamped gives $4{\pi }^{2}EI/{\ell }^2$ and clamped–free gives ${\pi }^{2}EI/(4{\ell }^2)$, equivalently an effective-length correction. A column buckles in the direction of smaller flexural rigidity. Euler buckling is a geometric instability within elastic response; yielding or crushing can occur first for a short thick strut. Large deflections, shear deformation and imperfections require refinements of the ideal slender-beam model.

A spatially uniform external gravitational acceleration moves a body's centre of mass without stretching it. Tidal forces arise from the variation of that acceleration across the body. Let a distant mass $M$ lie at distance $D$ along unit vector $\mathbf{n}$, and let $|\mathbf{r}|\ll D$ denote position relative to the body's centre. Subtract the centre's acceleration. To first order,The tidal tensor is $T_{ij}=GM(3n_i{n}_j-{\delta }_{ij})/D^3$. It stretches along the line to the source with eigenvalue $2GM/D^3$ and compresses in the two perpendicular directions with eigenvalues $-GM/D^3$. Its trace vanishes outside the source, consistently with Laplace's equation for the gravitational potential.

After removal of the constant and uniform-acceleration terms, the tidal potential isAn ideal equilibrium ocean on a sphere of radius $R$ follows an equipotential of the body's gravity plus this perturbation. For surface gravity $g$, the leading displacement isup to a constant fixing total water volume. The quadrupolar form produces bulges on both the near and far sides; the far-side bulge follows from differential acceleration, without postulating a repulsive gravitational force.

The Moon and Sun both raise terrestrial tides, with relative strength set mainly by $M/D^3$, rather than by the total gravitational pull $M/D^2$. Their contributions reinforce near new and full Moon to produce spring tides, and partly cancel near quadrature to produce neap tides. Earth's rotation relative to the tidal pattern gives a strong semidiurnal component in the simplest geometry. Real tide heights and timing depend on ocean depth, coastlines, basin resonances, Coriolis acceleration, friction and the declination of the forcing bodies. Tidal dissipation produces phase lags and torques, exchanging rotational and orbital angular momentum and gradually altering the spin and orbit. The equilibrium bulge is therefore a useful starting model, not an exact local tide forecast.

Inertia is a tensor. For a rigid body, distances between material points are fixed. Write each position relative to the centre of mass as $\mathbf{r}$, so its rotational velocity is $\mathbf{\omega }\times \mathbf{r}$. Integrating $\mathbf{r}\times (\mathbf{\omega }\times \mathbf{r})$ gives angular momentumThe rotational kinetic energy isThe inertia tensor is real and symmetric, so orthogonal principal axes diagonalize it. Angular momentum is parallel to angular velocity only along a principal axis, or within a subspace of equal principal moments. For translation plus rotation about the centre of mass, $T=M{v}_{\mathrm{CM}}^2/2+T_{\mathrm{rot}}$. Moving the reference to a parallel axis adds $M(d^{2}1-\mathbf{d}{\mathbf{d}}^{\mathsf{T}})$ to the tensor, the tensor form of the parallel-axis theorem.

Euler's equations. The external torque about the centre of mass obeys $\mathbf{\tau }=(d\mathbf{L}/dt{)}_{\mathrm{space}}$. In body axes, the transport rule givesThe principal moments are constant in that frame, henceThese are Euler's equations. With zero torque, spatial angular momentum and kinetic energy are conserved, although the body components generally vary. For distinct moments, small perturbations about rotation around the largest or smallest principal moment oscillate, whereas those about the intermediate moment grow. This is the intermediate-axis instability.

A lamina. Choose the lamina's plane as $z=0$. Then $I_{xz}=I_{yz}=0$, so the normal is a principal axis, andThis perpendicular-axis theorem applies to a planar mass distribution; it is not a general identity for a three-dimensional body. Rotate the in-plane axes to diagonalize the remaining symmetric two-by-two block. For a uniform rectangular lamina with sides $a,b$ and mass $M$, its central principal moments are $I_x=M{b}^2/12$, $I_y=M{a}^2/12$ and $I_z=M(a^2+b^2)/12$. If $a>b$, then $I_x<I_y<I_z$: torque-free rotation about the long in-plane axis or the normal is stable, while that about the short in-plane axis is unstable. Coincident moments in a square require the degenerate, rather than strict intermediate-axis, case. Under a constraint to rotate only about the normal, the motion reduces to the scalar equation ${\tau }_z=I_z{\dot{\omega }}_z$.

A symmetric top without torque. Let $I_1=I_2=I_{\perp }$ and $I_3=I_{\parallel }$. Euler's equations reduce toThus ${\omega }_1+i{\omega }_2$ rotates at signed body-frame frequency $\nu$, while the axial component stays fixed. The symmetry axis has constant angle to fixed spatial angular momentum. Indeed $\mathbf{L}=I_{\perp }\mathbf{\omega }+(I_{\parallel }-I_{\perp }){\omega }_3{\mathbf{e}}_3$, soIt precesses in space about $\mathbf{L}$ with angular speed $L/I_{\perp }$. This spatial precession and the body-frame frequency $\nu$ describe different motions and should not be identified.

A heavy top. If a symmetric top is pivoted a distance $\ell$ from its centre of mass, gravity gives torque ${\mathbf{r}}_{\mathrm{CM}}\times M\mathbf{g}$. Here the relevant moments are about the fixed pivot. The axial angular momentum and the vertical angular-momentum component are conserved, because gravity supplies neither an axial torque nor a torque about the vertical. Using tilt $\theta$, precession angle $\varphi$ and spin angle $\psi$, its Lagrangian isEliminating the two conserved cyclic momenta leaves one-dimensional motion in an effective tilt potential; bounded oscillation in tilt is nutation. Steady precession at $\dot{\varphi }=\Omega$ and fixed tilt satisfiesFor a rapidly spinning top on the slow branch, the first term is small andGravity then changes the direction of angular momentum rather than simply tipping the top straight down. The approximation fails near conditions where nutation or the fast-precession branch matters; the full equations retain those possibilities.

In SI base dimensions,Equating current, time, mass and length powers respectively yieldsThusDimensional analysis cannot determine the dimensionless proportionality constant.

For the edge-hinged cube with a central spherical cavity, use the parallel-axis theorem: for a parallel axis at perpendicular distance $d$ from the centre-of-mass axis, $I=I_C+M{d}^2$. Let $M=\rho a^3$. The centre lies $d=a/\sqrt{2}$ from the edge hinge, soAt equilibrium the centre is vertically below the hinge. Its gravitational torque at angular displacement $\theta$ is $-Mgd\sin \theta$. Therefore the small-angle equation is $I_h\ddot{\theta }+Mgd\theta =0$, withThe removed sphere has mass $M_s=\pi M/6$ and the same centre. Its hinge inertia isSuperposition leaves the centre at the same location but changes both inertia and weight:The frequency decreases: the removed central material contributes less hinge inertia per unit mass than the original cube average.

In the uniform-Earth chord oscillator, neglect the mass removed by the narrow tunnel and any rotation of the Earth. Inside a uniform sphere, the shell theorem gives $\mathbf{F}=-(4\pi G\rho m/3)\mathbf{r}$. On the constrained chord $\mathbf{r}=(d,y)$, its along-tunnel component isThe transverse component is balanced by the tunnel constraint. With release from $y_s=\sqrt{R^2-d^2}$ at rest,The far end is $-y_s$, so $t_{\mathrm{transit}}=\pi /{\omega }_0\approx 2530\mathrm{s}=42.2\min$, independent of chord offset in this ideal model. Here $\rho =5515\mathrm{kg}{\mathrm{m}}^{-3}$ and $G=6.67430\times 10^{-11}{\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-2}$ give ${\omega }_0=1.24171\times 10^{-3}{\mathrm{s}}^{-1}$.

With linear drag, $\ddot{y}+(b/m)\dot{y}+{\omega }_0^{2}y=0$. Inserting a complex exponential $e^{-(\gamma +i\omega )t}$ givesso in the underdamped regimeReal physical solutions combine conjugate modes. Release from rest specifically gives $y=y_s{e}^{-\gamma t}[\cos \omega t+(\gamma /\omega )\sin \omega t]$; a lone exponential with prefactor $y_s$ does not impose that initial velocity.

The critical gravity-tunnel settling time follows from critical damping, the fastest non-oscillatory release:Interpret the 1% tolerance relative to the initial displacement from equilibrium. With $s={\omega }_{0}t$, $(1+s)e^{-s}=0.01$ has $s=6.63835$. HenceThe equilibrium is reached only asymptotically; using a pure $e^{-{\omega }_{0}t}$ would miss the critical-damping polynomial factor and underestimate this time.