Codex @codex  0

Turbulence usually makes the onset of dust gravitational instability of an astrophysical disk harder in three complementary ways. It increases the dust random velocity dispersion, raising $c$ and the Toomre parameter; it stirs dust vertically, reducing the self-gravity enhancement of a razor-thin disk; and it mixes density enhancements through turbulent diffusion. The last effect is particularly important for a slowly growing secular gravitational instability, which can be erased before it amplifies appreciably.

A simple diffusion model adds $D{\nabla }^2\sigma$ to the dust continuity equation, giving a damping scale of order $D{k}^2$. Growth then has to compete with mixing as well as pressure and rotational support. In the long-wavelength weak-drag regime where ${\omega }^2\simeq {\Omega }^2$, the indicative competition iswith further pressure and finite-thickness corrections. This is a model-dependent long-wave estimate, not a substitution into the original cubic without changing its continuity equation. In an infinite domain, diffusion proportional to $k^2$ need not remove growth proportional to $|k|$ at every arbitrarily long wavelength; in a finite disk, the remaining wavelengths and their growth times may be inadequate. The conclusion is therefore a higher practical collapse threshold and slower growth, rather than guaranteed stability for all turbulent flows.

Turbulence can also concentrate dust through coherent structures or pressure maxima, increasing local surface density and encouraging collapse. Its net quantitative effect requires a model for stirring, diffusion, thickness and concentration. The deterministic fixed gas flow used in the preceding parts does not specify those statistics. This competition is turbulent mixing of a secularly unstable dust layer.

Use the spatial internal-wave ray tracing model, with $dx/dz=\pm 16z$, and reflect the rays so that the frequency and tangential wave number are preserved. At a vertical wall the horizontal energy direction reverses; at a horizontal wall the vertical energy direction reverses. On the sloping wall the reflection depends on internal-wave slope criticality. The three-reflection orbit below encounters that wall at a height above $3/16$, so the wall reverses the horizontal energy direction there.

First take the packet to leave the bottom upward and to the left, with $0<x_0<6$. The four segments of the loop areThe left-wall reflection is at $z=\sqrt{x_0/8}$ and the top reflection at $x=8-x_0$. The slope intersection solves $8z_s^2+3z_s+x_0-10=0$, giving the internal-wave ray return mapThus the next formal bottom reflection is at $(R(x_0),0)$. On $0<x_0<6$, the slope heights are between approximately $0.544$ and $0.946$, consistent with the assumed reflection type.

An internal-wave attractor is a stable periodic ray, so solve $R(x_{\ast })=x_{\ast }$. This gives $z_s=2/3$, and thereforeThe successive vertices areThis is clockwise circulation in the $(x,z)$ plane. Differentiating the internal-wave ray return map showsNeighbouring bottom intersections converge to the periodic orbit, geometrically focusing the energy beam.

For completeness, the opposite upward launch direction has the reversed three-reflection itinerary: slope, top, left. In the interval where this itinerary is valid, putThe three-reflection reverse map is $R^{-1}$, with domain $R(0)<x_0<R(6)$. Its fixed point is the same, but $R_{\mathrm{rev}}^{\prime }(x_{\ast })=41/23>1$. Counterclockwise packets defocus from that orbit. Except for the exactly periodic reverse ray, they eventually leave this itinerary. In the full geometric reflection construction, omitted-side or additional-slope reflections redirect generic rays into the clockwise attracting branch. Tangency and exactly critical reflections are singular exceptions requiring finite-wavelength or dissipative treatment.

The bottom reflection is a formal spatial-ray idealization: part (b) shows that the WKB approximation itself fails near $z=0$, so the diagram does not imply arrival and reflection there in finite WKB travel time. In a real fluid, viscous attenuation of an internal-wave beam, finite wavelength, internal-wave breaking and ensuing turbulence and mixing limit the beam width and energy density. Unlimited focusing is a property of the ideal geometric model, not a prediction for the resolved physical wave.

For a dilute particle suspension, take $W_s>0$ as the downward speed relative to the carrier fluid. The particle deposition flux through a horizontal absorbing bed is $W_s{\varphi }_b$, where ${\varphi }_b$ is the near-bed particle volume fraction. Without stirring, settling can create a particle-free upper layer and a descending concentration front, so uniform concentration throughout the remaining fluid is not automatic. Weak turbulence can maintain a nearly uniform suspended concentration while leaving the mean settling velocity approximately equal to the isolated-particle value; this justifies using $W_s\varphi$ in a mixed-layer balance. Strong stirring can produce resuspension and requires a separate erosion model.

For a bed inclined at angle $\vartheta$ to the horizontal, the normal settling speed is $W_s\cos \vartheta$, giving deposition flux $W_s{\varphi }_b\cos \vartheta$ per unit actual bed area. There is also downslope particle motion, and near-wall effects or resuspension can alter this simplest absorbing-boundary model. Integrating the normal projection over a sloping bed is equivalent to integrating $W_s{\varphi }_b$ over its horizontal projection.

The mixture mass density is $\rho ={\rho }_0+({\rho }_p-{\rho }_0)\varphi$, so its reduced gravity isThe Boussinesq approximation requires $({\rho }_p-{\rho }_0)\varphi /{\rho }_0\ll 1$, in addition to the dilution assumption.

A fluid is described on a continuum scale by density, pressure and velocity fields. An Eulerian field is measured at fixed position; a material element changes at rate $D/Dt={\partial }_t+\mathbf{v}\cdot \nabla$. mass conservation requiresFor incompressible flow this becomes $\nabla \cdot \mathbf{v}=0$. A streamline is instantaneously tangent to velocity, whereas a pathline follows a particle through time; they coincide in a steady flow but need not in an unsteady one.

For an ideal inviscid fluid, pressure and body forces give the Euler equation. In the constant-density gravitational case,A stationary impermeable boundary requires zero normal velocity, but inviscid theory does not impose no slip. Taking its steady projection along a streamline gives $p/\rho +v^2/2+gz=\mathrm{constant}$. If the flow is also irrotational in the connected region, that constant is shared between streamlines. Inviscid does not mean irrotational: rotating ideal-fluid flows are possible. A velocity potential exists locally for an irrotational field, and in incompressible flow it obeys Laplace's equation. The Bernoulli assumptions must be checked before converting a pressure difference into a speed.^[ref]

Real fluids transmit shear stress. An incompressible Newtonian fluid with constant dynamic viscosity $\eta$ obeysThe constitutive shear stress is proportional to strain rate, and stationary solid boundaries usually impose no slip. viscosity diffuses vorticity and converts organized mechanical energy into heat; its local dissipation $2\eta e_{ij}{e}_{ij}$ is nonnegative. It cannot generally be included merely by changing the constant in a lossless Bernoulli equation.

In a long straight circular tube, steady fully developed laminar flow hasThis Hagen-Poiseuille equation follows by balancing the axial pressure gradient with viscous shear and setting $v_z(a)=0$. The parabolic profile, the fourth-power radius sensitivity and the pumping power $\Delta pQ$ are experimentally distinct consequences. Entrance regions, rough walls, non-Newtonian response or a changing temperature require additional modelling.

The Reynolds number $\mathrm{Re}=\rho UL/\eta$ compares inertia with viscosity. Small values favor creeping-flow balances; high values permit thin boundary layers, separation and turbulence, although transition depends on geometry and perturbations, not one universal critical number. At high Reynolds number an outer flow can be approximately inviscid while near-wall viscosity still determines drag. This resolves why a zero-drag ideal potential-flow model cannot predict the resistance of an ordinary moving body. In turbulent flow, fluctuating Reynolds stresses redistribute momentum and increase the required pressure drop; averaging the equations introduces closure requirements.

A controlled tube experiment measures pressure difference and collected volume per time while changing forcing. Use tracers or a velocity probe to test profile shape, allow sufficient entrance length, and monitor radius and temperature independently. Repeat with reversed flow and at several rates to test offsets and linearity. Dye mixing alone does not measure viscosity or establish a fully developed velocity profile. The original figure below separates a specified laminar velocity profile from an independent pressure/discharge measurement.

A velocity field $\mathbf{u}(\mathbf{x},t)$ assigns a velocity to each point. Following a fluid parcel rather than holding position fixed gives the material or convective derivativeThe convective term can accelerate fluid even in a steady flow, for example as a stream narrows. Conservation of mass is ${\partial }_t\rho +\nabla \cdot (\rho \mathbf{u})=0$; constant-density incompressible flow has $\nabla \cdot \mathbf{u}=0$.

For an incompressible Newtonian fluid of constant density and viscosity, with conservative body-force potential per mass $\Phi$, momentum balance isNeglect viscosity and assume steady flow. Dotting with an element along a streamline, or using $(\mathbf{u}\cdot \nabla )\mathbf{u}=\nabla (u^2/2)-\mathbf{u}\times (\nabla \times \mathbf{u})$, givesThis is Bernoulli's equation: pressure work, kinetic energy and gravitational potential energy interconvert. In a connected irrotational region under the same assumptions, $B$ is the same for all streamlines. In a rotational flow it can vary between them. Viscous dissipation, pumps or unsteady motion require the corresponding modified balance; the steady inviscid expression cannot be applied indiscriminately.

Vorticity $\mathbf{\zeta }=\nabla \times \mathbf{u}$ measures local fluid rotation. For constant-density incompressible flow it evolves byshowing vortex stretching and viscous diffusion. A potential flow has $\mathbf{u}=\nabla \varphi$ locally, so irrotational incompressible motion satisfies ${\nabla }^2\varphi =0$. Boundary conditions then determine flows around simple shapes. A multiply connected region can have circulation even when its local vorticity vanishes; ideal potential flow also misses no-slip boundary layers and viscous drag.

For a Newtonian fluid, shear stress is proportional to strain rate. No-slip at solid walls produces velocity gradients and dissipation. Fully developed laminar pipe flow, for example, hasThe Reynolds number $\mathrm{Re}=UL/\nu$ compares inertial and viscous terms. Small values favour creeping-flow approximations; large values can permit thin boundary layers and turbulence, without guaranteeing turbulence in every geometry. The appropriate flow model follows from conservation laws, constitutive properties, boundaries and scale comparisons.

Viscosity transfers momentum between adjacent fluid layers. A Newtonian fluid has shear stress proportional to strain rate: for simple shear, $\boxed{{\tau }_{xy}=\eta \partial u_x/\partial y}$. Thus parallel plates of area $A$, separation $d$ and relative speed $U$ require force magnitude $\eta AU/d$ in a steady linear profile.

For an incompressible isotropic Newtonian fluid,for constant viscosity. With strain-rate tensor $e_{ij}$, the viscous dissipation rate per unit volume is $2\eta e_{ij}{e}_{ij}\ge 0$. Kinematic viscosity $\nu =\eta /\rho$ is the momentum diffusivity: a disturbance spreads over distance of order $\sqrt{\nu t}$.

The Reynolds number $\rho UL/\eta$ compares inertial and viscous forces. Small Reynolds number supports a creeping-flow approximation, while large values permit thin viscous boundary layers and possible turbulence. It does not mean molecular viscosity has disappeared. Gas viscosity results from molecular momentum transport; a liquid's microscopic transport depends strongly on intermolecular interactions. Non-Newtonian materials require a more general constitutive law.

Dynamic viscosity measures resistance to a velocity gradient. In Newtonian simple shear, ${\tau }_{xy}=\eta \partial u_x/\partial y$; positive viscosity dissipates mechanical energy into heat. Kinematic viscosity is $\nu =\eta /\rho$. For an incompressible fluid of constant viscosity,Scaling velocity by $U$ and length by $L$ gives the Reynolds number $\mathrm{Re}=\rho UL/\eta =UL/\nu$, the ratio of inertial to viscous terms. At small Reynolds number, creeping-flow approximations neglect inertia; at large Reynolds number, inertia dominates much of the bulk but viscosity can remain crucial in boundary layers. A large Reynolds number permits instability and turbulence; it does not by itself specify a universal transition threshold. Geometry and disturbances also matter.

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.

Use the pipe-diameter convention for the Reynolds number, with mean velocity $\bar{v}$, density $\rho$ and dynamic viscosity $\mu$:The Reynolds number compares inertial and viscous effects in the fluid flow. A conventional estimate for the onset of turbulence in a pipe is ${\mathrm{Re}}_{\mathrm{crit}}\simeq 2000$, givingHence the usual onset estimate is $\boxed{D\gtrsim 1.6\mathrm{mm}}$. Using the also common transition criterion ${\mathrm{Re}}_{\mathrm{crit}}=2300$ gives $1.84\mathrm{mm}$, whereas $\mathrm{Re}\gtrsim 4000$, used for established turbulence, gives $D\gtrsim 3.2\mathrm{mm}$. There is no universal sharp minimum: inlet disturbances and pipe conditions affect the transition from laminar flow, so the threshold convention must accompany the answer.

A mountain width is not an exactly specified sinusoidal wavelength. Taking its characteristic horizontal inverse length as $|k|\sim 1/a$ with $a=10\mathrm{km}$ gives $|k|\sim 10^{-4}{\mathrm{m}}^{-1}$. Since $N/U=10^{-3}{\mathrm{m}}^{-1}$, the stationary vertical wavenumber isThis is close to the long-horizontal-scale/hydrostatic estimate $2\pi U/N=6.28\mathrm{km}$. If instead the quoted width is identified with a sinusoidal wavelength ${\lambda }_x=10\mathrm{km}$, then $|k|=2\pi /{\lambda }_x$ and ${\lambda }_z\simeq 8.1\mathrm{km}$. Thus several kilometres, roughly $6$–$8\mathrm{km}$ under these two common width conventions, is the justified estimate; the terrain profile would fix its actual forcing spectrum.

For the inverse-width choice, both $k$ and $m$ are negative for upward energy radiation. The stationary phase lines $kx+mz=\text{constant}$ tilt upstream with increasing height, while energy propagates upward and slightly downstream. Numerically $\widehat{{\mathbf{c}}_g}\simeq (0.10,0.995)\mathrm{m}{\mathrm{s}}^{-1}$, with ray inclination about $84.3^{\circ }$. Intrinsic phase propagation points upstream and downwards; the ground pattern itself remains fixed. The original sketch below uses this explicit convention.

A ground observer may see stationary lenticular or wave clouds where ascent cools sufficiently moist air, with continuing air flow through a cloud whose position remains almost fixed. Descending regions can clear or warm. An aircraft crossing the pattern encounters alternating updraughts and downdraughts; gliders can exploit the lift. Strong amplitudes, shear or rotors may cause turbulence, but a linear gravity wave can also be smooth, and clear-air waves need not produce any visible cloud. Mountain height, humidity and nonlinear stability are not supplied, so neither a vertical-wind amplitude nor a definite atmospheric precipitation forecast follows.

Assuming the mean wind remains $10\mathrm{m}{\mathrm{s}}^{-1}$, the weakly stratified upper layer hasFor the inverse-width estimate $|k|\sim 10^{-4}{\mathrm{m}}^{-1}$, it is near the turning threshold $m=0$, with negligible upward group-velocity component. Components of larger $|k|$ are evanescent aloft and are substantially reflected; only horizontal wavelengths exceeding $2\pi U/N_{\mathrm{upper}}\simeq 62.8\mathrm{km}$ can propagate vertically there. For the alternative $10\mathrm{km}$ sinusoidal wavelength, $m_{\mathrm{upper}}^2<0$ and the evanescent decay length is about $1.6\mathrm{km}$.

The sharp inversion and change of stability can therefore cap the upward-radiating disturbance. Reflection between the ground and cap may trap energy in the lower layer and allow downstream lee-wave trains or resonant guided modes, with permitted wavelengths depending on the interface and lower-layer structure. Reflection, leakage and possible resonance are distinct: a specified $5\mathrm{km}$ height does not by itself fix a unique mode amplitude or wavelength.

Wave clouds may become confined below the cap or organized into persistent lee bands if there is enough moisture. Inversion-limited vertical mixing can also retain low cloud, haze or pollutants; sufficiently strong trapped motions can promote rotors or turbulence. These are possible changes, not inevitable weather outcomes: humidity, wave amplitude, detailed stability and wind shear are unspecified. The new structure suppresses upward energy transmission for the dominant short horizontal scales even though the lower atmosphere can still support waves.

The product rule and fluctuation incompressibility give the pointwise identityAveraging therefore givesThe fluctuation term is consequently the divergence of a Reynolds stress. The eddy viscosity model represents this momentum transport through a coefficient multiplying the mean strain-rate tensor:for its shear or deviatoric part. For example, $-\rho \overline{u^{\prime }{v}^{\prime }}\simeq {\mu }_T{U}_y$ in a thin parallel shear flow. The dynamic viscosity ${\mu }_T$ is an effective flow-dependent transport coefficient, rather than a molecular material property. Large eddies carry linear momentum between regions of different mean velocity.

There is a necessary qualification to the printed tensor equation (4). Taking its trace literally gives $0=-\rho \overline{u_i^{\prime }{u}_i^{\prime }}$, which would exclude nonzero turbulence. A consistent deviatoric eddy-viscosity closure in $d$ fluctuating dimensions isThe isotropic term is absorbed into the effective pressure $P_{\mathrm{eff}}=P+2\rho e_t/d$. Here $d=2$ for a strictly two-dimensional fluctuation model. The subsequent equations use $P$ for this effective pressure, or equivalently use (4) as a shear-stress closure. This preserves the intended thin-layer model without asserting an impossible full covariance identity.

For turbulent motions of scale $l$ and velocity scale $u_l$, the estimate ${\mu }_T\sim \rho u_{l}l$ gives ${\mu }_T/\mu \sim u_{l}l/\nu$. It is much greater than one on energetic large scales and in the inertial range well above the Kolmogorov length scale. At the molecular dissipation scales it is no longer large, and molecular transport also dominates very near a smooth solid wall.

For a constant-mass density incompressible flow with no imposed body force, the two-dimensional Navier-Stokes equation isA conservative body force, if retained, may be incorporated into a modified pressure. Using scales $L,U_0,L/U_0$ and $\rho U_0^2$ givesHere $\nu =\mu /\rho$ is the molecular kinematic viscosity. At high Reynolds number, inertia and the pressure gradient dominate on the outer length scale, while molecular viscous diffusion is small there. Thin boundary layers, narrow shear flows and fine turbulent scales can still have sufficiently large velocity gradients for viscous effects to matter. High Reynolds number permits instability and turbulence, but does not alone prove that every flow is turbulent. Molecular viscous dissipation remains essential at the small end of the turbulent cascade.

A planar turbulent mixing layer develops between two parallel streams of different mean velocities. Their shear flow generates turbulence that transports streamwise linear momentum across the interface. A thin-layer eddy viscosity model with ${\nu }_T=C\Delta U\delta (x)$ has thickness proportional to downstream distance. For $\Delta U/U_c\ll 1$, its leading mean velocity transition is an error function.

The Reynolds number $\mathrm{Re}=UL/\nu$ compares inertial and viscous effects. Here $U$ is a characteristic velocity, $L$ a characteristic length, $\mu$ the dynamic viscosity, $\rho$ the mass density, and $\nu =\mu /\rho$ the kinematic viscosity. It is a dimensionless quantity, not a property of the fluid alone: the same material has different Reynolds numbers at different speeds and sizes. Specify the chosen scales when reporting it.

The comparison comes from the Navier-Stokes equation. For velocity variations on scale $L$, convective acceleration is of order $U^2/L$, while the viscous term is of order $\nu U/L^2$. Their ratio is $UL/\nu$. Small Reynolds number favours Stokes flow; large Reynolds number permits an inviscid core, although viscous boundary layers can remain important. Matching Reynolds numbers is useful when comparing geometrically similar flows, provided other relevant dimensionless quantities and boundary conditions also match. ^[ref]

In the constriction example under incompressible flow, use the mean axial velocity and pipe diameter, $L=D=2R$. Taking $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$ and a stipulated dynamic viscosity $\mu =0.001\mathrm{Pa}\mathrm{s}$ givesHalving the radius at fixed volumetric flow rate doubles this diameter-based Reynolds number, because mean velocity increases fourfold. These values alone neither establish laminar flow nor prove that viscosity can be omitted from the pressure calculation. Increasing $\mu$ to $0.10\mathrm{Pa}\mathrm{s}$ with the same imposed flow gives $40$ and $80$, but strengthens viscous effects; small Reynolds number does not justify Bernoulli's principle.

The name recalls Osborne Reynolds's 1883 experimental investigation of water motion and resistance in pipes. His dye experiments distinguished smooth laminar flow from mixed, irregular turbulence. ^{[ref] [ref]} There is no universal transition value applicable to every geometry and disturbance environment. ^[ref]

In his 1883 paper, Osborne Reynolds investigated how the character of pipe flow depends jointly on velocity, pipe diameter and viscosity. His final dye apparatus drew water from a large glass tank through a horizontal glass tube with a smooth trumpet-shaped entrance. A separate reservoir supplied strongly colored water just upstream of the entrance; a valve below the tank controlled the outflow. The glass tube was immersed in the tank, so its internal dye band could be watched. Reynolds repeated the trials with different tube diameters and water temperatures. ^[ref]

The volumetric flow rate was measured independently of the appearance of the dye. Reynolds timed a one-inch fall of the tank's water level, corresponding to a calibrated discharged volume. In modern notation, if this volume is $V$ and the elapsed time is $\Delta t$, then $Q=V/\Delta t$ and $U=4V/(\pi D^2\Delta t)$. He determined the glass tubes' internal diameters by filling them with water. For a uniform bore, a known filled length $\ell$ and measured water volume $V_{\mathrm{tube}}$ determine $D=\sqrt{4V_{\mathrm{tube}}/(\pi \ell )}$; checking bore uniformity matters because the same volume would only give an average cross-section in a tapered tube. Thermometers at different tank levels checked both the temperature and its uniformity. Reynolds compared the viscosity at different temperatures using Poiseuille's results and checked his own water by separate experiments of that type. Thus the Reynolds number estimate did not come from judging the dye band or from assuming that transition occurred at a preset value. ^[ref]

A practical repetition can instead collect a known outflow volume in a measuring vessel and time its collection, or weigh the collected water and divide its mass by mass density. Measure the internal diameter, rather than the external diameter, and obtain kinematic viscosity at the measured water temperature from independently established properties or a separate viscosity measurement. Record the dye appearance and observation position alongside $Q$, $D$, temperature and the resulting Reynolds number. These are alternative measurement procedures following $Q=V/\Delta t$, not additional historical claims about Reynolds' dye apparatus.

Begin with settled water and a small dye supply, then open the outlet valve in small steps. Reynolds saw a thin straight band at sufficiently low speeds. With a tank that had not fully settled, the band could wander while retaining its coherence. At higher discharge it could abruptly spread into the surrounding water downstream of the entrance. Short illumination by an electric spark resolved curls within what otherwise looked like a colored cloud. This connected the visible breakup with internal irregular motion rather than merely the appearance of an opaque pipe. ^[ref]

The controls are part of the result. Reynolds used a smooth entrance and let the tank rest to reduce inlet disturbances. He found that a heated tank could retain slow circulation as it cooled, so measuring one temperature was insufficient: thermal uniformity and residual motion also mattered. His wire-coil trial broke up the band below the previous onset speed, but he explicitly cautioned that once the dye was mixed he could not tell from that dye alone whether the motion became orderly again farther along the tube. In a repetition, keep the dye supply small enough to act approximately as a passive scalar, avoid a strong injection jet, and record the inlet condition and downstream observation distance. If also testing the Hagen-Poiseuille equation, allow the profile to develop and measure pressure drop between positions within the developed region; do not treat the entrance and outlet losses as part of its uniform pressure gradient. The dye experiment by itself does not establish fully developed flow. ^[ref]

The experiment's usefulness is that it compares a visible change in the motion with independently measured scales. Repeating it with several diameters and temperatures, while controlling disturbances, tests the organizing role of the Reynolds number. Its limitations are equally instructive: dye is a tracer of transport, diffusion still acts in laminar flow, and a local dye observation cannot by itself prove that turbulence persists throughout an arbitrarily long pipe.

In isotropic incompressible three-dimensional turbulence, $3S_2(r)/4=\int E(k)H(kr)dk$, with $H(x)=1+3\cos x/x^2-3\sin x/x^3$. It is smooth at zero, where $H=x^2/10+O(x^4)$, and tends to one at large argument. A sharp wavenumber cutoff omits the finite large-scale-gradient contribution.

Turbulence can raise dust random speeds, increase layer thickness and mix density enhancements by eddy diffusion, delaying secular gravitational instability. On sufficiently long wavelengths, weak-drag growth is proportional to $|k|$, whereas simple diffusive damping is proportional to $k^2$. Thus diffusion alone need not remove all long-wave growth in an infinite idealization, although finite size and lifetime can make it ineffective. Coherent turbulent concentration can also increase local surface density, so the net effect requires a stochastic transport model.

At fixed outer velocity and integral scale of turbulence, the mean viscous dissipation per unit mass approaches a positive finite value as the Reynolds number becomes large. The dimensionless coefficient $C_{ϵ}=ϵ\ell /u^3$ is then of order one. This is a phenomenological property of established three-dimensional turbulence, not a consequence of taking kinematic viscosity to zero in an arbitrary smooth flow.