Flow Physics

Reynolds Number and Flow Transition, Through Visualization

A cylinder wake sits near Strouhal 0.2, pipe flow turns unsteady from Re 2300 to 4000, and a tenfold rise in Re multiplies a DNS mesh by about 178.

2026.03.25

Poiseuille Flow: Analytical Solution Derivation and Finite Difference Validation

At h = G = μ = 1 the centerline speed is 1/2 and the flow rate is 2/3. FD L2 starts at 2.78e-16 and never leaves roundoff, so ~2.0 is not a measured order.

2026.03.28

Bugs Swim in Honey — The World of Flow Changed by Reynolds Number

A whale at Re 2.5e8 and a bacterium at 3e-5 share the same water. The factor 1/Re on the viscous term decides which balance actually survives.

2026.04.21

Why Water Streams Break into Droplets — Plateau–Rayleigh Instability and Breakup Regimes

The fastest Plateau–Rayleigh wave sits near 9.02 radii. Weber and Ohnesorge then decide whether that jet stays in fat drops or turns to mist.

2026.04.23

Two Faces of a Whirl — Free Vortex, Forced Vortex, and Helmholtz

A free vortex revolves with zero vorticity except at the center. The Rankine core spins at 2U/R and meets the 1/r tail smoothly at the core radius r = R.

2026.04.28

Water Boils at 20°C Too — The Physics of Cavitation

At 20°C vapor pressure is 2.34 kPa. If a venturi throat drives static pressure under that line, water boils and the cavitation number turns negative.

2026.04.30

When the Stream Breaks — Reynolds' Critical Number and the Layered Turbulent Boundary

Pipe transition is usually quoted near Re 2300, yet a quiet pipe can stay laminar to 50,000. The wall then splits into sublayer, buffer, and log layer.

2026.05.05

Why a Spinning Ball Curves — Circulation, Magnus, and the Kutta–Joukowski Theorem

At U=10 and R=0.5 the three circulations send top-surface Cp from -3 to -15. At critical spin the lower stagnation point sits at Cp=1.

2026.05.07

Above Mach 1, You Must Widen to Speed Up — The Laval Nozzle Paradox

Area ratio 2.5 exits near Mach 2.44 on the supersonic branch, static pressure at 6.3 percent of stagnation. Raise back pressure and a shock stands inside.

2026.05.12

The Paramecium Swims Through Honey — Where Viscosity Comes From and Why Fourier and Fick Wrote the Same Line

At 1 atm the mean free path is about 68 nm and the mean speed about 470 m/s. One flux, μ ≃ ρ c̄ ℓ / 3, is the line Fourier and Fick used for heat and mass.

2026.05.14

Three Faces of Porous Flow — Darcy, Forchheimer, and Brinkman

On gravel-bed constants, viscous and Forchheimer losses meet at a Darcy velocity of 0.1 m/s, and a straight pump curve under-reads the drop after that.

2026.05.19

Two Functions Added, and a Cylinder Appeared — Complex Potential and d'Alembert's Paradox

Adding a doublet to a uniform stream draws a cylinder without sketching one. With zero circulation the pressure integral leaves only roundoff for drag.

2026.05.21

Why Smokestack Plumes Curl Up Instead of Lying Flat — The Jet in Crossflow and Its Counter-Rotating Vortex Pair

With A=1.6 and B=0.33, the centerline at x=10d sits at 5.44d, 8.66d, and 13.78d for r=2, 4, and 8. Doubling r raises height by about 1.59, not by 2.

2026.05.26

Three Portraits of a Flow — Why Streamlines, Pathlines, and Streaklines Diverge

Drop the contraction to 0.25 and throat speed jumps 4x while the local term stays at 0. With V0=0.6 and omega=2, streamline y never leaves ±0.6.

2026.05.28

Turbulence Runs Three Rulers at Once — The Kolmogorov Cascade and the -5/3 Law

A tenfold rise in Reynolds shrinks η/L to about 0.18 of its value, and a lab tunnel at Re=10^4 already asks a 3D DNS for about 10^9 cells.

2026.06.02

Fluid That Warms Up While Standing Still — Natural Convection and Rayleigh 1708

A 0.30 m wall at 60 C in 20 C air prints Ra 7.03e7, Nu 54.9, and a heat flux of 197.6 W/m^2, still on the laminar side of 1e9.

2026.06.04

Why a Rougher Ball Flies Farther — The Drag Crisis and Boundary-Layer Separation

At 70 m/s and Re = 1.99e5 the same ball gives 2.049 N smooth and 0.488 N dimpled. Separation at 82 degrees versus 120 opens that gap.

2026.06.09

Heating It Hotter Made It Cooler — The Boiling Curve and Critical Heat Flux

Zuber's CHF for water at 1 atm is 1.259 MW/m². Rohsenow's nucleate flux reads 1.094 at 20°C and 3.692 at 30°C, already past that ceiling.

2026.06.11

To Go Supersonic, Widen the Pipe — The De Laval Nozzle and Choking

Area ratio 4 designs to Mach 2.940 and pe/p0 of 0.0298. Gamma 1.4 locks the choked mass-flow factor near 0.5787, and back pressure stays out of it.

2026.06.18

Measure the Drag, Get the Heat Transfer — The Reynolds Analogy

At Re_x = 1e5 on a flat plate, Cf/2 stays at 1.050e-3, and St times Pr to the 2/3 matches that value at Pr = 0.7, 1, and 7.

2026.06.23

Why Gasoline Engines Stop at a Compression Ratio of 12 — Otto vs Diesel

Same compression ratio: Otto is 56.5% at r=8 and 63.0% at r=12, Diesel at cutoff 2 is 49.0% and 56.7%. Knock, not the heat input, stops gasoline near 12.

2026.06.25

Where the Ideal Gas Starts to Lie — Compressibility Factor Z and the van der Waals Equation

Nitrogen at 100 bar and 200 K has Z = 0.7905. The ideal-gas density, 168.4 kg/m3, is about 26 percent lighter than the real 213.0.

2026.06.30

The Shock Wave in Your Kitchen Sink — Froude Number and the Hydraulic Jump

At Fr 2.86 a 0.2 m supercritical sheet jumps to Fr 0.42 and a depth of 0.714 m, scattering 0.238 m of specific energy into the roller.

2026.07.02

Chaos That Begins at the Bottom of a Pot — Rayleigh–Bénard Convection and the Lorenz Equations

Rolls switch on at Rayleigh number 1708 with rigid walls. In the Lorenz truncation that heating, written as rho, stops settling once rho passes 24.74.

2026.07.09

Waves Hanging in the Sky — The Kelvin-Helmholtz Instability

On air and water the first wave shows up near 6.6 m/s, about 1.69 cm long, growing at 4.89 per second. Below that wind the surface stays calm.

2026.07.14

Why a Spinning Ball Curves — the Magnus Effect and Circulation

Random spin drifts a musket ball sideways. Circulation Γ=2πa²ω sets the side force, and spin 2 misses by 6.90 m, not twice the 3.71 m printed at spin 1.

2026.07.16

A Single Ripple Outruns Its Own Group — Wave Dispersion and Group Velocity

Deep-water gravity waves move energy at half the phase speed. The slowest ripple is near 1.73 cm and 0.232 m/s, where gravity balances surface tension.

2026.07.21

How Far Must Water Flow Before It Becomes a Parabola — Entrance Length and Hagen–Poiseuille

In laminar pipe flow a flat inlet needs about 0.05 Re diameters to become a parabola, and 64/Re is not valid until that station.

2026.07.23

Where Turbulence Stops Speeding Up a Flame — Damköhler's Wrinkling Hypothesis and Bending

At u'/S_L = 12 the large-scale wrinkling factor stays near 4.03, while the small scale drops to 2.00 once Ka passes 29.

2026.07.28

Four Shocks Beat One — The Oblique Shock Train in a Supersonic Inlet

At Mach 3 one normal shock recovers 0.3283 of stagnation pressure. Four 8° ramps plus a terminal shock lift that figure to 0.8235.

2026.08.06

Viscosity Is Not What Kills a Vortex — Newton's 1/r and the Lamb–Oseen Core

At r=0.5 the speed is unchanged to four decimals after 50 diffusion times. At n=1 the viscous force is zero; only a wall removes circulation.

2026.08.11

Sending the smoke down 1 m cost 7 Pa of draft — Franklin's stove and the buoyancy budget

One metre of down-leg removes 6.16 Pa and the longer chimney gives back less, 38.32 Pa instead of 39.15. Net draft falls from 44.72 Pa to 37.72 Pa.

2026.08.20

The Coriolis Force Did Zero Work, Yet the Impeller Added 1800 J/kg — Energy in a Rotating Frame

A rotating frame conserves rothalpy, not total enthalpy. Coriolis work stays at 0 J/kg while impeller head still rises by 1800 J/kg.

2026.08.27

Where Total Pressure Dropped 1,080 Pa, Dissipation Was Zero — When the Bernoulli Constant Crosses Streamlines

Inside a spinning vortex core the total pressure falls 1,080 Pa while viscous dissipation stays exactly zero. A Rankine vortex marks where a head deficit stops meaning loss.

2026.09.03

49.9% of the Dynamic Coefficients Came Out Negative — the Germano Identity and the Averaging Operator

The Germano identity puts the Smagorinsky coefficient negative at half the points. Where the average is taken decides whether the run stays stable.

2026.09.15

2.0000 in laminar flow, 1.0322 in turbulent — the coefficient in front of the velocity head

A mean-velocity head is not the kinetic energy a pipe actually carries. Alpha is 2.0000 for laminar Poiseuille flow and 1.0322 at Re = 10^6.

2026.09.17

Four times the speed, still 19.47 degrees — where water depth blows the Kelvin wedge open to 90

The Kelvin half-angle in deep water stays at 19.47° even if speed quadruples. A shallow bottom opens it toward 90° once the depth Froude number reaches 1.

2026.09.24

Refining the Grid 8x Never Moved 2.4674 — Two Constants Fourier Left Behind in 1807

The convergence objection that got Fourier's 1807 memoir rejected comes back on today's grids as two stubborn constants: 8.9490% and pi^2/4.

2026.09.29

Deriving μ and k from One Formula Put Pr Off by 1.98× — the 9/4 Behind a 123.5 K Wall Temperature

Deriving μ and k from one mean free path pins air's Pr at γ = 1.400 against a measured 0.707. Eucken's 9/4 drops the Couette wall gap from 123.5 K.

2026.10.01

One χ=50, two temperatures: 1759 K and 818 K — where a flamelet library folds

One χ=50 returns 1759 K and 818 K because the S-curve folds. A 61-point flamelet shows why a progress variable C makes that table single-valued.

2026.10.06

A 1 μm Aluminum Particle Heated 4.6× Slower Than the Continuum Formula — Sherman Interpolation in the Knudsen Transition Regime

At 1500 K and 1 atm a 1 μm aluminum particle sits at Kn = 0.67. Nu = 2 shortens heat-up by 4.58; Sherman's harmonic mean puts that time back.

2026.10.15

Dividing the Reverse Rate by Kp Split 99.6% of AlO at 2500 K — The (p°/RT)^Δn Hidden in the Equilibrium Constant

Plugging Kp into kr = kf/K skews the equilibrium of any reaction with Δn≠0 by (p°/RT)^Δn. An AlO dissociation reactor reproduces the SI and CGS versions of the mistake.

2026.10.22

A Small Wobble in Mixture Fraction Drops the Mean Flame Temperature by 964 K — Presumed β-PDF and the Infinity at Z=0

At stoichiometric mean Z, the β-PDF mean temperature sits 964 K below T(Z̄). When a<1, a 100-cell uniform-grid integral also overshoots by 115 K. Here is why, with code.

2026.10.27