Unsteady Stokes equation 2026-10-14
For variations on a time scale , the ratio of local inertia to viscous resistance is . Steady Stokes flow at each instant requires both and . If only the first is small, retain local acceleration:
For an oscillating sphere with angular frequency , the second condition becomes on the radius scale . The Unsteady Stokes equation is still linear, but needs initial data and has memory; its velocity field is not determined by instantaneous boundary velocities alone. See E. J. Hinch, Slow Viscous Flow, example sheet 1, question 1 for the distinction between steady and unsteady limits of an oscillating sphere.
A U-tube differential manometer contains a connected column of liquid with one gas connection above each liquid meniscus. Write for the left and right pressures, and for the meniscus heights on an upward vertical coordinate. Assume constant mass density and downward gravitational acceleration of magnitude ; neglect the gas column's mass density and differences in capillary pressure. The signed vertical length difference then gives
Larger left pressure means a lower left liquid meniscus: implies . Equal pressures give equal meniscus heights.
Figure 1. . Original schematic of two same-height pressure tappings connected to a resting water column. The constructed example uses and , so the vertical difference is . Dashed lines mark both meniscus heights and the equal-pressure datum; the dotted bottom line identifies two equal-height points within the connected liquid. The arrow in the pipe indicates fluid flow; the manometer liquid has no continuing flow at equilibrium. Details.
The balance follows directly from hydrostatic pressure. At a common height below both surfaces, the two points in the connected homogeneous liquid have equal pressure:
Subtracting cancels the common bottom depth and leaves . The same-height rule follows from hydrostatic equilibrium, : a horizontal displacement changes no pressure, and for constant the quantity is constant throughout the connected liquid. Equal height alone is insufficient for unrelated vessels or for fluid flow with significant acceleration.
For an original constructed example, choose water of stipulated mass density , gravitational acceleration , and vertical length difference . Then
The SI pressure unit is the pascal: the product has units . Swapping the two connections gives and reverses which liquid meniscus is lower. In equal-bore arms, volume conservation makes each surface move from its initial equal-height position; using only one surface's displacement would halve the inferred pressure. With unequal arm areas , a left fall and right rise instead obey , while the vertical length difference is still .
If the right port is open to the local atmosphere, is the left gauge pressure. Only with that reference and known atmospheric pressure can one obtain absolute pressure from
For an explicitly stipulated atmospheric pressure of , the example gives absolute pressure. If the right port is connected to another vessel or a pressure tapping, the same reading is a differential pressure and supplies no absolute pressure without knowing that reference.[ref]
When a vacuum vessel's window shatters, a nearly atmospheric initial pressure difference gives inward pressure force on frontal area . Dividing by the glass mass estimates its initial collective acceleration. An ideal atmospheric pressure-work budget for vessel volume would give a collective kinetic energy estimate if assigned entirely to the fragments. This is a model-dependent energy budget, not a universal bound on every fragment's speed: air motion, fragmentation, pressure equalization and wall impacts alter the transfer, and unequal sharing can accelerate a small fragment faster than the collective estimate.
Velocity Verlet method 2026-10-18
For an ordinary differential equation with acceleration independent of velocity, the velocity Verlet method advances position and velocity by
Here is the time step and . For sufficiently smooth , the global error on a fixed time interval is , so halving reduces the leading global error by a factor of four once the asymptotic regime is reached.
Vertical ascent under linear drag Created 2026-10-05 Updated 2026-10-06
For constant downward gravitational acceleration and linear drag , the upward velocity solves with . Its first zero occurs at the displayed time. The dimensionless control variable is , and the limit gives .
Viscous stress tensor Created 2026-09-28 Updated 2026-10-14
The viscous stress tensor is the part of the Cauchy stress tensor associated with deformation rates. For an isotropic Newtonian fluid in incompressible flow,
Here is pressure, is the identity matrix, is dynamic viscosity, and is the strain-rate tensor. The convention is ; denotes transpose and denotes trace. Across a plane with unit normal vector , the viscous traction is . The tensor makes this rule work for any plane orientation. [ref]
In the worked Couette flow, , so
In particular, while the shear rate is ; the factor in recovers . On the upper face of a fluid slice, gives viscous traction ; on the lower face it is . Constant stress can therefore transmit a force between the plates even though its divergence, and the net viscous force on each slice, vanish.
Compare planar extension with solid-body rotation :
Extension changes lengths and produces viscous normal stresses . Rotation preserves distances and gives zero viscous stress, despite its nonzero velocity gradient and vorticity . For and , extension gives and ; rotation gives zero. A rotating fluid may still have a pressure gradient supplying its centripetal acceleration.
Figure 1.
Deformation produces Newtonian viscous stress; rigid rotation does not
. Blue arrows show instantaneous velocity. The dashed square is a material patch at ; its orange outline at is obtained from the exact trajectories of planar extension and solid-body rotation, with . The marked corner makes rotation visible. Both motions preserve area, but only extension changes the patch's shape. Details.
For an isotropic Newtonian fluid in three-dimensional compressible flow, separate shape change from volume change:
The first term is deviatoric stress, with zero trace; is volume viscosity, also called bulk viscosity. For uniform expansion , and , so the shear contribution vanishes and the viscous stress is . Thus dynamic viscosity and volume viscosity describe distinct responses, and incompressible flow removes only the volume-change term. [ref]
Wave equation Created 2026-09-24 Updated 2026-10-14
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 is
Here is position, is time, is the disturbance, and subscripts denote partial derivatives. For a wave equation on a string, is transverse displacement; in linear acoustics, it can be a pressure perturbation. In several spatial dimensions, becomes the Laplacian . 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 and velocity . A classical solution has continuous partial derivatives through order two; and suffice for D'Alembert's formula to produce one. A finite string additionally needs boundary conditions. These statements concern the unforced model with constant ; damping, spatially varying coefficients and nonlinear terms change the solution formula.
A travelling wave carries its entire profile to the right at wave speed ; 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 , 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 : 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 , 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.