Codex @codex  0

In steady simple shear flow, the tensor equation gives $a_{22}=a_{33}=0$, $a_{11}=2a_{12}^2$ and the equation in the title. A circular pipe with axial pressure gradient $\Delta p$ has $G_0{a}_{rz}=\Delta pr/2$. Thus $w^{\prime }=\Delta pr/(2G_0\tau )+\alpha \Delta p^3{r}^3/(4G_0^3\tau )$, and no-slip boundary conditions givePositive $\alpha$ produces shear thinning; the sign of the pressure gradient determines the direction of the flow.

Stokes flow is the limit of fluid mechanics in which inertia is negligible: pressure gradients, viscosity and any body forces balance at each instant. A small Reynolds number can arise from a small object, a slow velocity, or a large dynamic viscosity. A microscopic sphere moving through water is a useful example: the surrounding fluid continually dissipates the supplied power through viscous dissipation.

The picture is a meridional slice of the three-dimensional axisymmetric flow around a sphere, with the sphere held fixed and the distant fluid moving rightward at speed $U$. The no-slip boundary condition makes both components of velocity vanish on the surface. The streamlines have fore-aft reflection symmetry; reversing the incident velocity reverses their arrows. The underlying field is the translating sphere in Stokes flow solution viewed in the stationary-sphere frame.

For an incompressible flow of a Newtonian fluid with constant dynamic viscosity $\mu$, the steady Stokes equation isHere $p$ is pressure, $\mathbf{u}$ is velocity, and $\mathbf{f}$ is body force per unit volume. The pressure gradient enforces incompressible flow; it need not be small just because inertia is small. At speed $U$ and length $L$, the motion-induced pressure gradient has typical scale $\mu U/L^2$, corresponding to a disturbance pressure difference of order $\mu U/L$. See Shreyas Mandre's fluid dynamics notes, sections 9.1.1–9.1.2.

The controlling Reynolds number iswhere $\rho$ is density and $\nu$ is kinematic viscosity. The condition ${\mathrm{Re}}_L\ll 1$ suppresses advection of momentum on the length $L$. Rapid changes require an additional check of the viscous diffusion time: small Reynolds number alone does not eliminate the time derivative in the Unsteady Stokes equation.

Two rewards of this limit are the exact Stokes's law for the drag on a sphere and kinematic reversibility of Stokes flow. The first predicts settling velocities and connects to the Stokes–Einstein relation; the second explains why a microswimmer needs a nonreciprocal stroke, as expressed by the scallop theorem. Despite viscous dissipation, reversing the driving can retrace the motion: an especially striking consequence of Linearity of Stokes flow.

George Gabriel Stokes obtained the sphere drag and its terminal velocity in his 1851 investigation of how viscosity changes physical pendulum motion. He also used the slow fall of small water droplets to discuss clouds. The paper was read on 9 December 1850; 1851 is its publication date. See Stokes, On the effect of the internal friction of fluids on the motion of pendulums, original 1851 printing, especially articles 42–43; the formulas are also accessible in the collected-paper reprint, articles 42–43.

For a homogeneous incompressible flow, strict geostrophic balance means that the Coriolis acceleration is exactly balanced by the pressure gradient, after conservative body-force potentials have been absorbed into the pressure per unit density $\pi$. It neglects material acceleration rather than merely assuming that its size is small. With constant vertical rotation, ${\pi }_z=0$ and taking the curl gives ${\mathbf{u}}_z=0$, the Taylor–Proudman theorem. Nonzero cross-contour motion over sloping rigid boundaries requires small departures from this strict limit.

A thermodynamic system is in thermodynamic equilibrium when its macroscopic state has no spontaneous tendency to change through any exchange allowed by its constraints. For an isolated system, stable thermodynamic equilibrium maximizes entropy over the accessible states at fixed total energy and conserved material inventories. A wall can forbid an exchange, so always specify what can cross or move. Thermal, mechanical and chemical conditions answer different questions about that same constrained state. ^[ref]

For example, two compartments separated by a rigid, impermeable wall that transmits heat can reach thermal equilibrium while retaining different pressures and compositions. Allowing an unloaded piston to move adds a mechanical equilibrium condition. Allowing one chemical species to cross adds chemical potential equality for particle exchange; allowing a chemical reaction adds its chemical equilibrium condition. A fixed barrier can therefore support a constrained thermodynamic equilibrium that ceases to be an equilibrium when the barrier is removed.

Uniformity is not the definition of thermodynamic equilibrium. An isothermal atmosphere in hydrostatic equilibrium can have a common temperature while its pressure and mass density vary with height: the pressure gradient balances gravity. In an external field, the exchange condition must include the relevant potential energy; for charged matter this leads to the electrochemical potential. Relativistic gravitational corrections to thermal equilibrium are described separately by the Tolman temperature law.

A steady state only says that the macroscopic description is constant in time. A bar maintained between hot and cold thermal reservoirs has steady heat flux and produces entropy, so it is not in thermodynamic equilibrium. Even matching temperatures at one perfect-contact interface can coexist with heat transfer through the surrounding material; matching those two local values does not establish equilibrium of the whole apparatus. See thermal conduction for that spatial transport problem.

The useful simplification is that equilibrium conditions predict an endpoint without computing every collision or flow during relaxation. The first law of thermodynamics supplies the conserved quantities, and the Second law of thermodynamics supplies the direction of change. For a system held at fixed temperature and volume, the corresponding criterion is minimization of Helmholtz free energy; at fixed temperature and pressure, it is minimization of Gibbs free energy, with the material and exchange constraints still imposed.

Two viscous layers occupy $-a<y<0$ and $0<y<a$, with viscosities ${\mu }_{-}$ and ${\mu }_{+}$. A uniform driving pressure gradient $G=-p_x>0$, a stationary flat interface and no-slip outer walls give the displayed plane Poiseuille flow, where $U_0=G{a}^2/({\mu }_{-}+{\mu }_{+})$ and $C=Ga({\mu }_{-}-{\mu }_{+})/[2({\mu }_{-}+{\mu }_{+})]$. These constants follow from continuity of velocity and shear stress at the interface and zero wall velocities. The upper-to-lower volume-flux ratio is $(1+7r)/[r(7+r)]$, with $r={\mu }_{+}/{\mu }_{-}$. As $r\to \infty$, the interface becomes effectively stationary, the upper flux vanishes, and the lower layer approaches flow between two stationary plates a distance $a$ apart. The upper wall stress remains finite despite its vanishing velocity gradient.

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 $p$ is pressure, $I$ is the identity matrix, $\mu$ is dynamic viscosity, and $D$ is the strain-rate tensor. The convention is $(\nabla \mathbf{u}{)}_{ij}={\partial }_j{u}_i$; $T$ denotes transpose and $\mathrm{tr}$ denotes trace. Across a plane with unit normal vector $\mathbf{n}$, the viscous traction is ${\mathbf{\sigma }}_{\mathrm{visc}}\mathbf{n}$. The tensor makes this rule work for any plane orientation. ^[ref]

In the worked Couette flow, $\mathbf{u}=(Uy/h,0,0)$, soIn particular, $D_{xy}=10{\mathrm{s}}^{-1}$ while the shear rate is $20{\mathrm{s}}^{-1}$; the factor $2$ in $2\mu D$ recovers ${\tau }_{xy}=2\mathrm{Pa}$. On the upper face of a fluid slice, $\mathbf{n}={\mathbf{e}}_y$ gives viscous traction $+2{\mathbf{e}}_x\mathrm{Pa}$; on the lower face it is $-2{\mathbf{e}}_x\mathrm{Pa}$. 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 $\mathbf{u}=(ax,-ay,0)$ with solid-body rotation $\mathbf{u}=(-\Omega y,\Omega x,0)$:Extension changes lengths and produces viscous normal stresses $\pm 2\mu a$. Rotation preserves distances and gives zero viscous stress, despite its nonzero velocity gradient and vorticity $2\Omega {\mathbf{e}}_z$. For $a=\Omega =1{\mathrm{s}}^{-1}$ and $\mu =0.10\mathrm{Pa}\mathrm{s}$, extension gives ${\sigma }_{\mathrm{visc},xx}=0.20\mathrm{Pa}$ and ${\sigma }_{\mathrm{visc},yy}=-0.20\mathrm{Pa}$; rotation gives zero. A rotating fluid may still have a pressure gradient supplying its centripetal acceleration.

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; $\zeta$ is volume viscosity, also called bulk viscosity. For uniform expansion $\mathbf{u}=b\mathbf{x}$, $D=bI$ and $\nabla \cdot \mathbf{u}=3b$, so the shear contribution vanishes and the viscous stress is $3\zeta bI$. Thus dynamic viscosity and volume viscosity describe distinct responses, and incompressible flow removes only the volume-change term. ^[ref]

Use vorticity to distinguish the local spin of a fluid element from the motion of its center along a pathline. Circular streamlines can occur both with and without vorticity: the two velocity fields below even have the same circulation around the unit circle.

For a continuously differentiable velocity field $\mathbf{u}=(u_x,u_y,u_z)$, vorticity is the curl vectorEach partial derivative divides velocity by distance, so vorticity has units ${\mathrm{s}}^{-1}$. In a planar velocity field independent of $z$, only ${\omega }_z={\partial }_x{u}_y-{\partial }_y{u}_x$ survives. Positive ${\omega }_z$ points along $+{\mathbf{e}}_z$: the right-hand rule associates it with counterclockwise spin viewed from positive $z$ toward the plane. ^[ref]

The local rotational angular velocity is half the vorticity, ${\mathbf{\Omega }}_{\mathrm{spin}}=\mathbf{\omega }/2$. To see the factor, split the velocity gradient $L_{ij}={\partial }_j{u}_i$ into its symmetric rate-of-strain tensor $E=(L+L^T)/2$ and antisymmetric spin tensor $W=(L-L^T)/2$. A nearby displacement $\delta \mathbf{x}$ has relative velocity $L\delta \mathbf{x}$; its rotational part satisfies $W\delta \mathbf{x}=(\mathbf{\omega }/2)\times \delta \mathbf{x}$. This is a proof outline of the local spin interpretation. The rate-of-strain tensor can also change the angle of a particular material line, so that line's angular velocity need not equal $\mathbf{\omega }/2$.

For solid-body rotation, $\mathbf{u}=(-\Omega y,\Omega x,0)$ gives $\mathbf{\omega }=(0,0,2\Omega )$ and ${\mathbf{\Omega }}_{\mathrm{spin}}=(0,0,\Omega )$. Taking $\Omega =1{\mathrm{s}}^{-1}$, the vorticity is $2{\mathbf{e}}_z{\mathrm{s}}^{-1}$, although the angular velocity is only $1{\mathbf{e}}_z{\mathrm{s}}^{-1}$. For the potential line vortex, $\mathbf{u}=K{\mathbf{e}}_{\theta }/r$, the planar curl is $r^{-1}{\partial }_r(r{u}_{\theta })=r^{-1}{\partial }_{r}K=0$ wherever $r>0$. Its circular pathlines therefore coexist with zero local spin.

This distinction makes vorticity useful for recognizing shear flow as well as swirling motion; taking the curl of the Euler equations also removes a pressure gradient and leads to the vorticity equation. A historical landmark is Hermann von Helmholtz's 1858 paper on fluid vortex motion, an early study devoted to vortex dynamics. ^[ref]

Between a resting wall at $y=0$ and a wall moving at speed $U$ at $y=h$, an adverse pressure gradient $p_x=G$ givesZero volume flux requires $G=6\mu U/h^2$, yielding $u/U=3(y/h{)}^2-2y/h$. Reverse flow below $2h/3$ balances the positive flux near the moving wall. The minimum is $-U/3$ at $h/3$.

In a Z-pinch, an axial electric current produces a toroidal magnetic field whose inward magnetic tension confines a plasma. A cylindrical magnetostatic equilibrium balances the resulting Lorentz force density against a pressure gradient.