Codex @codex  0

Net gravitational effect of thermal and compositional density contrasts between a slab and its surroundings, including phase changes.

A model of the Sun obtained from stellar structure equations with specified composition, opacity and energy generation. Uniform density is useful for scale estimates but cannot reproduce the strong central concentration.

For repulsive contact interactions in a three-dimensional isotropic harmonic trap, neglecting the kinetic term in the Gross–Pitaevskii equation gives $n(r)=(\mu -m{\omega }^2{r}^2/2)/U_0$ inside $R=\sqrt{2\mu /(m{\omega }^2)}$ and zero outside. Normalization yields $N=8\pi \mu R^3/(15U_0)$, hence $\mu =(\hslash \omega /2)(15Na/a_0{)}^{2/5}$ and $R=a_0(15Na/a_0{)}^{1/5}$, with oscillator length $a_0=\sqrt{\hslash /(m\omega )}$. The density is parabolic but the wavefunction is its square root. Neglected kinetic energy rounds the edge and supplies a tail.

A uniform overdensity evolves as a separate spherical pressureless region until shell crossing. In an Einstein-de Sitter background, formal collapse has linearly extrapolated threshold ${\delta }_c=(3/5)(3\pi /2{)}^{2/3}\simeq 1.686$. This is distinct from the actual nonlinear density. Virialization and nonspherical tides require further physics.

Take $v=-{\mathcal{X}}_z$, $w={\mathcal{X}}_y$, constant density gradient $S<0$, and meridional eddy density flux $F=\sin (ly)[-1_{z<0}]$. If the divergence of the eddy momentum flux is zero, bounded Eulerian circulation and decaying residual mean circulation give $\mathcal{X}=\sin (ly)\chi (z)$ with $\chi =-A_0(1-e^{\kappa z}/2)$ below zero and $\chi =-A_0{e}^{-\kappa z}/2$ above, where $A_0=-1/S$ and $\kappa =Nl/|f_0|$. The transformed streamfunction is ${\mathcal{X}}^{\ast }=\mathcal{X}+F/S$, yielding the displayed ${\chi }^{\ast }$ and two oppositely rotating cells. Its jump represents a horizontal delta-function transport at the discontinuous flux termination. A smooth flux cutoff replaces that sheet by a thin regular return flow. Density flux alone does not determine either circulation if the momentum-flux divergence is left unspecified.

For a steady incompressible flow of a Newtonian fluid with constant density $\rho$ and constant dynamic viscosity $\mu >0$, neglecting inertia in the Navier-Stokes equation givesThe first equation is local force balance; the second is mass conservation. The term $\mu {\nabla }^2\mathbf{u}$ is the divergence of the viscous part of the Newtonian fluid stress tensor. The body force $\mathbf{f}$ has units of force per unit volume, $\mathrm{N}{\mathrm{m}}^{-3}$, and equals $\rho \mathbf{g}$ for uniform gravity. Appropriate boundary data for velocity or traction, together with suitable distant conditions for an exterior flow, complete the problem; a rigid surface normally imposes the no-slip boundary condition.

With uniform gravity, write $p=p_h+p^{\prime }$ where the hydrostatic pressure satisfies $\nabla p_h=\rho \mathbf{g}$. Then the disturbance equation is $-\nabla p^{\prime }+\mu {\nabla }^2\mathbf{u}=0$. This removes gravity from the fluid disturbance equation, while the hydrostatic pressure still supplies buoyancy on an immersed body. Omitting it from the total force balance would give the wrong Stokes settling velocity.

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.

Take downward as positive and let gravity have magnitude $g$. For a rigid sphere of density ${\rho }_p$ in otherwise stationary fluid of density ${\rho }_f$, its weight minus buoyancy is $\frac{4}{3}\pi a^3({\rho }_p-{\rho }_f)g$. At terminal velocity, this balances the Stokes drag law:This signed downward settling velocity is negative for a buoyantly rising sphere and zero for equal densities. For fixed density contrast and dynamic viscosity, the submerged weight grows as $a^3$ while linear drag at fixed speed grows as $a$; their balance gives the settling-velocity radius scaling $U_s\propto a^2$. Meanwhile,Thus doubling the radius multiplies the predicted Stokes settling speed by four and the particle Reynolds number by eight. This makes the formula useful for separating particles by size, while also providing its own check against extrapolation to large particles.

For a computed water example, take ${\rho }_f=1000\mathrm{kg}{\mathrm{m}}^{-3}$, $\mu =1.00\times 10^{-3}\mathrm{Pa}\mathrm{s}$, ${\rho }_p=2500\mathrm{kg}{\mathrm{m}}^{-3}$, $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, and $a=10.0\mu \mathrm{m}$. These are specified rounded inputs rather than a claim of exact water properties at every temperature. They giveThe diameter is $20.0\mu \mathrm{m}$, and both the net downward weight and the upward drag force are $6.16\times 10^{-11}\mathrm{N}$. At this terminal velocity, falling $1.00\mathrm{mm}$ takes $3.06\mathrm{s}$. At $a=20.0\mu \mathrm{m}$ the prediction is $1.308\mathrm{mm}{\mathrm{s}}^{-1}$ with ${\mathrm{Re}}_d=0.05232$; at $a=50.0\mu \mathrm{m}$ it is $8.175\mathrm{mm}{\mathrm{s}}^{-1}$ with ${\mathrm{Re}}_d=0.8175$, which is no longer a convincing small-Reynolds number case. The figure labels these as predictions of the Stokes settling velocity formula and shows where its inertia check fails.

The calculation describes the eventual terminal velocity, not the release transient governed by particle acceleration and possibly the Unsteady Stokes equation. Walls, hydrodynamic interactions, departures from a rigid sphere, or a different surface condition require a different resistance. For sufficiently small particles, Brownian motion can obscure individual settling trajectories; the formula still specifies the gravitational drift within its assumptions. The balancing argument already appears in article 43 of ^[ref].

For an isolated rigid sphere of radius $a$ moving at velocity $\mathbf{U}$ relative to a distant uniform Newtonian fluid, the drag force exerted by the fluid isThe assumptions are steady incompressible flow, constant dynamic viscosity $\mu$, an unbounded surrounding fluid, the no-slip boundary condition, and a small diameter-based particle Reynolds number, ${\mathrm{Re}}_d=2a{\rho }_{f}U/\mu \ll 1$. Integrating the surface traction in translating-sphere Stokes flow gives the factor $6\pi$: the disturbance pressure contributes $2\pi \mu aU$ to the drag magnitude and the viscous shear stress contributes $4\pi \mu aU$. See Richard Fitzpatrick, Axisymmetric Stokes Flow Around a Solid Sphere, equations 10.120–10.129.

The result is linear drag: doubling $U$ doubles the drag force, and doubling $a$ at fixed $U$ also doubles it. The required power is $6\pi \mu a{U}^2$, all supplied to viscous dissipation. Density enters the check of the particle Reynolds number, although it does not appear in the leading drag formula.

Nearby walls or other particles change the hydrodynamic resistance matrix and invalidate the isolated-sphere coefficient, even at very small Reynolds number. Finite Reynolds number introduces inertia corrections; the Oseen approximation resolves the distant region where the leading Stokes flow approximation ceases to be uniform. A clean gas bubble has tangential slip and instead obeys clean-bubble Stokes drag, with coefficient $4\pi \mu a$ in its spherical, negligible-internal-viscosity limit. See H. K. Moffatt, Six Lectures on General Fluid Dynamics, section 3.4 for the sphere solution and its small-Reynolds number limitation. The continuum mechanics and boundary assumptions must also remain valid as the particle becomes very small.

Balancing this Stokes drag law against gravity and buoyancy gives the Stokes settling velocity of an isolated rigid sphere.

A streakline at observation time $T$ is the set of present positions of all fluid elements that passed through one fixed release point $\mathbf{a}$ during a specified history of injection. For continuous release over ${\tau }_0\le \tau \le T$, solve a different initial value problem for each release time $\tau$:Thus a pathline varies $t$ for one particle; a streakline varies $\tau$ while holding $T$ fixed. Interrupted release gives only the corresponding subsets of the streakline. A moving nozzle requires a time-dependent release position instead, so the standard fixed-point definition should not silently be applied to it. Continuous release makes streaklines easy to observe and useful for showing where material from one source has reached. ^[ref]

In the analytical velocity field $(U,at)$ of pathline, a particle released at the origin at time $\tau$ followsObserve all particles together at time $T$, after continuous release from $\tau =0$:For $U=1\mathrm{m}{\mathrm{s}}^{-1}$, $a=1\mathrm{m}{\mathrm{s}}^{-2}$ and $T=2\mathrm{s}$, the streakline is the downward-curving parabola $y=2x-x^2/2$ in metre coordinates. Releases at $\tau =0,1,2\mathrm{s}$ occupy $(2,2)$, $(1,1.5)$ and $(0,0)\mathrm{m}$, respectively. At $x=1\mathrm{m}$, the streakline lies above the pathline $y=x^2/2$ and below the origin's streamline $y=2x$. All particles nevertheless have the same present velocity $(1,2)\mathrm{m}{\mathrm{s}}^{-1}$. The streakline slope $dy/dx=a\tau /U$ generally differs from the present velocity slope $aT/U$; they agree at the nozzle. The oldest particle is the downstream endpoint, which explains why the pathline and streakline meet there without being the same curve.

Each panel fixes one observation time. Coloured dots are the present positions of particles with the indicated release times; dashed pathlines show those particles' earlier histories. The black streakline connects present positions of the continuously released family, including release times between the sampled dots. After $T=1\mathrm{s}$, particles already present keep moving and new ones enter at the origin, so the later streakline is not a frozen record of the earlier one. Injection beginning only at $\tau =1\mathrm{s}$ would show just $0\le x\le 1\mathrm{m}$ of the $T=2\mathrm{s}$ curve. Both panels are diffusion-free analytical models, not experimental photographs.

Kline's contracting-channel demonstration switches from marking an element to continuous bubble production at a marking wire: the resulting trace is a streakline. The apparatus therefore changes the curve being displayed by changing the marking history, even with the same channel. A wire produces a family of traces; conceptually each fixed emitting location supplies its own streakline. Individual bubbles are markers rather than perfect fluid elements: Jean Hertzberg's teaching notes describe hydrogen production by electrolysis at a thin cathode wire and warn that larger bubbles have more buoyancy. ^{[ref] [ref]}

A dye version uses a small fixed injection port upstream in a transparent water channel, a dilute premix of dye with the channel water, illumination and a fixed optical camera. Start continuous injection at a recorded time ${\tau }_0$ and take a short exposure at $T$: the visible filament approximates a streakline, not the accumulated pathline of one particle. Repeat with the same port and imaging geometry after changing the flow with time; compare the filament against local directions estimated from short displacements of sparse markers at $T$. Agreement in a steady flow and possible disagreement when the flow changes are the predictions being tested; no measured curves are claimed here. The protocol follows Hertzberg's advice to use small ports, little injected volume and a dilute premix with matched density and viscosity. Check for nozzle disturbance by reducing injection rate, for buoyancy by watching the premix in still water, and for diffusion by comparing filament widths at different ages. Keep exposure short enough to suppress motion blur. Dye concentration approximates a passive scalar only when it has negligible feedback on the velocity field; finite port width, diffusion and out-of-plane motion limit how literally a visible filament can be read as a geometrical streakline. ^[ref]

A stress glut tensor describes a seismic source through an equivalent body-force density $\mathbf{f}=-\nabla \cdot \mathbf{\Gamma }$. For a symmetric stress glut tensor supported away from the boundary, integration by parts makes its coupling to a normal mode equal to $\int \mathbf{\Gamma }:{\mathbf{e}}^{\ast }dV$, where $\mathbf{e}$ is that mode's strain.

For an elastic filament with reference arclength $s$, curvature $c$, axial strain $ϵ$ and excess twist rate $\Omega$, a harmonic energy density is $(K_b{c}^2+S{ϵ}^2+K_t{\Omega }^2)/2$. Bending and torsional rigidities have units of energy times length; the stretch modulus has units of force. Factoring out $k_{B}T$ therefore makes $A,C$ lengths but $B$ an inverse length. Axial strain under tension is approximately $F/S$. Twist is controlled by torque and rotational boundary conditions, not force alone. A three-dimensional thermally fluctuating worm-like chain has persistence length $K_b/(k_{B}T)$.

For $|a|\ll \mu$, the equal-density ion-neutral cubic has a rapidly damped mode ${\omega }_d=-2i\mu +i{a}^2/(4\mu )+O(a^4/{\mu }^3)$ and two propagating modes ${\omega }_{\pm }=\pm a/\sqrt{2}-i{a}^2/(8\mu )+O(a^3/{\mu }^2)$. The reduced phase speed uses the total density ${\rho }_i+{\rho }_n=2{\rho }_i$. The exact displacement ratio is ${\xi }_n/{\xi }_i=\mu /(\mu -i\omega )$. It is nearly $-1$ for the rapidly damped counterflow and nearly $1\pm ia/(\sqrt{2}\mu )$ for the propagating modes. Thus the latter have small ion-neutral slippage of relative size $|a|/\mu$, responsible for ambipolar damping.

The buoyancy of the mantle beneath a descending slab. Thermal, compositional or melt-related density anomalies can modify pressure and traction on the slab, potentially changing its contact with the overriding plate. It differs from the slab's own slab buoyancy. Inferring it from low seismic velocities requires independent constraints because temperature, fluids, melt and composition can produce overlapping signatures.

For relativistic dust with nowhere-zero density, stress-energy conservation implies zero four-acceleration. If its normalized four-velocity $u$ has hypersurface orthogonality, then $u^{♭}\wedge d{u}^{♭}=0$ and ${\iota }_{u}d{u}^{♭}=0$. The interior product of the first identity gives $d{u}^{♭}=0$. The Poincaré lemma supplies $u^{♭}=-dt$ locally; spatial coordinate functions constant along the smooth flow of $u$ give $u={\partial }_t$, $g_{tt}=-1$ and $g_{ti}=0$. The resulting synchronous coordinates describe the dust until the local coordinate construction fails.

For uniform density between concentric radii $R_1,R_2$, subtract the missing solid core's central moment of inertia and then shift the axis by $R_2$. Here $M=4\pi \rho (R_2^3-R_1^3)/3$. The thin-shell limit gives $I=5M{R}_2^2/3$, while a full solid sphere gives $7M{R}_2^2/5$.

A correction for the difference between real surrounding relief and the assumed simple topographic slab. Its computation needs terrain geometry and an assumed density, and can matter in rough surface topography.

A unit of mass per length, equal to $1g/km$. Force divided by a linear density stated in tex gives tensile tenacity in newtons per tex.

For two ideal gas layers at a common interfacial pressure, the lower layer has smaller density precisely when $M_{-}/T_{-}<M_{+}/T_{+}$, where $M$ is mean molar mass. That ordering is unstable to an infinitesimal interchange. A heavy condensable constituent can therefore inhibit thermal convection in a light background gas even when the lower layer is hotter. This local criterion assumes no phase change during the initial displacement.

Space Science Reviews 134, pp. 155–190. The article discusses plasma emission, absorption, ionization and practical X-ray diagnostics. Its institutional NED version includes a cooling-function section with explicit density normalization.