Codex @codex  0

Use $v=|u_r|>0$ for the radial speed and distinguish the central mass $M$ from the Mach number $\mathcal{M}=v/c$. In each smooth nonzero-flow region, mass conservation, the adiabatic equation of state, and the Bernoulli equation giveHere $K>0$ is constant on a smooth isentropic flow branch, $c$ is the adiabatic sound speed, and $\mathcal{E}$ is the Bernoulli function. Either sign of $u_r$ is possible. The positive-radius nondimensionalization in the question requires $\mathcal{E}>0$; flows with $\mathcal{E}\le 0$ do not have these positive $r_s,x,\lambda$ and require another scaling.

Indeed, the supplied adiabatic process equations give $D(p{\rho }^{-3/2})/Dt=0$, so $K$ is spatially constant along a steady nonzero radial fluid flow. The adiabatic sound speed is $c^2=(\partial p/\partial \rho {)}_s$, and the specific enthalpy is $h=c^2/(\gamma -1)=2c^2$. The radial Euler momentum equation, $vdv/dr=-{\rho }^{-1}dp/dr-GM/r^2$, integrates to $\mathcal{E}$ because $dh=dp/\rho$ on this isentropic flow branch.

Put $D=[j(3K/2{)}^2{]}^{2/5}$. Eliminating the mass density from the continuity equation givesWith $y={\mathcal{M}}^{2/5}$, the Bernoulli equation becomesTherefore the physical constants of the spherical polytropic flow with adiabatic exponent three halves are

For $t>0$, $F^{\prime }(t)=4(t^5-1)/t^2$, so $F$ decreases from infinity to its unique minimum $F(1)=5$, then increases to infinity. This completely determines the algebraic curves. For $0<\lambda <1$, there are two disconnected branches at every $x>0$, one subsonic flow with $y<1$ and one supersonic flow with $y>1$. For $\lambda =1$, the line $y=x$ and a decreasing branch meet at $(1,1)$. Their slopes are $+1$ and $-1$, since the leading expansion is $F(1+\eta )=5+10{\eta }^2+O({\eta }^3)$. For $\lambda >1$, solutions exist only where $F(x)\ge 5\lambda$: if $x_{-}<1<x_{+}$ solve $F(x_{\pm })=5\lambda$, there is a forbidden interval $x_{-}<x<x_{+}$. At either endpoint the two branches turn at $y=1$ and have an infinite slope.

A smooth sonic point needs $F^{\prime }(x)=F^{\prime }(y)=0$; hence a regular sonic transition requires $\lambda =1$ and $r=r_s$. The critical speed of a polytropic flow is then $c_s^2=GM/(2r_s)=2\mathcal{E}$, and the allowed mass flux is $j=(2\mathcal{E}{)}^{5/2}{r}_s^2/(3K/2{)}^2$. Here the subscript on $c_s$ denotes evaluation at the sonic point, rather than the isothermal sound speed convention used in the next solution.

The simple transonic branch is $y=x$. Direct substitution givesThus the radial velocity has constant magnitude, while the mass density and pressure scale as $r^{-2}$ and $r^{-3}$. For an outward fluid flow, this is a subsonic flow becoming a supersonic flow as its adiabatic sound speed decreases.

On the other transonic branch, $y$ decreases as $x$ increases. At large $r$, the terms $4/y$ and $x^4$ dominate, giving $y\sim 4x^{-4}$. At small $r$, the terms $y^4$ and $4/x$ dominate, giving $y\sim 4^{1/4}{x}^{-1/4}$. The resulting asymptotic expansions areFor inward fluid flow, these are the nearly static reservoir at infinity and the inner free fall region of Bondi accretion. The specific enthalpy is smaller than $GM/r$ in the inner limit, which explains free fall. Reversing the radial velocity gives a decelerating outflow with the same profiles. The formal small-$r$ limit describes the point-mass model: if the spherical gravitating body has a finite surface, the exterior solution stops there.

A stationary, nonradiative normal shock wave conserves the mass flux and the sum $v^2/2+h$, while the Newtonian gravitational potential is continuous across its thin layer. Hence $\mathcal{E}$ and $r_s$ are unchanged. The entropy production in a perfect-gas shock increases $K=p/{\rho }^{3/2}$, soMore explicitly, the Rankine-Hugoniot condition for $\gamma =3/2$ gives, with upstream $m={\mathcal{M}}_1^2>1$,Thus the shock wave jumps vertically, at the same $x$, from the supersonic flow to a subsonic flow on a larger-$\lambda$ curve. The figure illustrates an outward flow on $y=x$ with a shock wave at $x=1.4$; the downstream subsonic flow has $\lambda \simeq 1.224$ and continues to larger $x$. The shock wave radius is an additional boundary-data choice, rather than being determined by the smooth equations alone.

Take $v>0$ to mean the magnitude of the inward radial velocity. In the Gaussian units used for the magnetic pressure, approximate the surface density of a disk by $\Sigma \simeq 2H\rho$. The mass accretion rate and free-fall speed then giveAt the magnetospheric truncation radius, balance the magnetic pressure of the magnetic dipole field with the specified ram pressure:ThusDiscarding the convention-dependent numerical factor givesThe magnetic pressure grows inward as $r^{-6}$, faster than the ram pressure $\propto r^{-5/2}$, so the field disrupts the flow inside this radius. The numerical coefficient depends on the vertical averaging and on treating the inflow as radial free fall, as stipulated here; a detailed thin disk boundary model need not have the same coefficient.

The phrase equivalence principle covers several related claims. The weak equivalence principle concerns the universality of free fall for test particles of negligible self-gravity, with the same initial conditions and no relevant nongravitational forces. In a Newtonian notation,Composition-independent motion requires the ratio of passive gravitational mass to inertial mass to be universal. Its common value can then be absorbed into the normalization of $\mathbf{g}$; the empirical assertion is universality, not the numerical convention used to express the ratio.

The Einstein equivalence principle includes that universality and extends the claim to local nongravitational experiments through local Lorentz invariance and local position invariance. A sufficiently small freely falling laboratory should recover the nongravitational laws of special relativity. The strong equivalence principle extends the scope to local gravitational experiments and bodies with appreciable gravitational self-energy. An experiment establishing one part of the weak statement does not automatically establish these stronger claims.

The motivating elevator comparison makes both the power and the limitation clear. Release a small object in a freely falling cabin. The cabin and object share the local free fall, so the object initially floats relative to the cabin. In a cabin supported against gravity, a released object instead accelerates relative to the floor, as it does in a suitably accelerated cabin in flat spacetime. This suggested replacing a universal gravitational force by geometric inertial motion. It does not identify a finite planetary gravitational field with an accelerated frame everywhere: separated falling objects can have a measurable relative acceleration.

In general relativity, the precise local geometric statement is that one can choose normal coordinates at an event $p$ withThis eliminates the connection contribution to the geodesic equation at that point. It does not generally set the Riemann curvature tensor to zero. In normal coordinates the first metric corrections away from $p$ are quadratic in the coordinate separation, and geodesic deviation measures the associated tidal forces. Equivalence is a local approximation for finite experiments whose extent and duration make these corrections negligible; it is not an assertion that curvature has been transformed away.^[ref]

There are direct experimental consequences beyond the elevator thought experiment. An Eötvös experiment compares bodies of different composition on a torsion balance, using a differential signal to test the universality of gravitational acceleration. MICROSCOPE performed an analogous test in orbit using electrostatic control of titanium-alloy and platinum-alloy test masses. Its first published result in 2017 constrained their differential free fall; this was already available before the examination. It tested a specific weak-equivalence prediction, not every possible gravitational law.^[ref] Comparisons of different atomic clocks can instead test local position invariance, including the expected gravitational redshift. These measurements constrain related but distinct parts of the framework.

The principle remains important as a restriction on how matter couples to gravity and as an organizing guide for experimental tests. It does not, on its own, derive the Einstein field equations. Universally coupling matter to a metric tensor still leaves choices for that metric's action; the equivalence principle does not uniquely fix field equations. For example, a metric gravitational action containing a scalar curvature squared term can retain the same minimal matter coupling while changing the gravitational dynamics. Selecting the Einstein-Hilbert action requires further dynamical assumptions and empirical success. The present role of the principle is consequently precise and limited: it supports local relativistic physics and universal matter coupling, while curvature, gravitational field equations and their quantitative predictions remain additional content.

For a light test particle of mass $m$ outside a spherical source of mass $M$, Newton's law of universal gravitation gives $\mathbf{F}=-GMm\hat{\mathbf{r}}/r^2$. It is a central force, so its torque about the source vanishes and angular momentum $\mathbf{L}=m\mathbf{r}\times \mathbf{v}$ is conserved. The motion therefore lies in a fixed plane when $L\ne 0$. Its gravitational potential energy is $-GMm/r$, and conservation of mechanical energy givesHere $ϵ$ is specific orbital energy and $h$ is specific angular momentum. A Kepler orbit has polar equationThe angle is measured from periapsis. The conic section is an ellipse for $ϵ<0$ and $0\le e<1$, a parabola for $ϵ=0$ and $e=1$, and a hyperbola for $ϵ>0$ and $e>1$. The circle is the $e=0$ elliptical case. The special case $h=0$ is radial motion and must be handled as a degenerate limit. Bound ellipses have the source at one focus, which is Kepler's first law.

For a bound orbit of semi-major axis $a$, $ϵ=-GM/(2a)$. Combining this with the energy equation gives the vis-viva equation:The particle therefore moves faster near periapsis and slower near apoapsis. Conservation of angular momentum also gives constant areal velocity, $dA/dt=h/2$, expressing Kepler's second law. Using the ellipse's area and $h^2=GMa(1-e^2)$ gives the orbital periodwhich is Kepler's third law.

A circular orbit follows from centripetal acceleration $v_c^2/r=GM/r^2$, giving $v_c=\sqrt{GM/r}$ and $T=2\pi \sqrt{r^3/(GM)}$. The minimum escape velocity from radius $r$ follows from the marginally unbound condition $ϵ=0$:These ideal orbits neglect drag force, other gravitating bodies and changes in the central source. For two comparable masses, the relative two-body problem uses $G(M+m)$ instead of $GM$. An orbiting object is continuously in free fall; the inward gravitational force changes its direction even when its speed is constant. These inverse-square results assume the physical variation of gravity with radius, rather than extending a near-surface constant-$g$ approximation to an entire orbit.

Introduce the dimensionless quantities $\tau =\mu t/r_0$ and $w=\mu v/(g{r}_0)$. ThenThe Taylor series at $t=0$ givesThus the initial slope is $g$, as for free fall. At late times,The oblique asymptote is $v=gt/4+g{r}_0/(4\mu )$, approached from below. Direct differentiation givesThe speed increases monotonically with decreasing slope, from $g$ towards $g/4$. It does not reach a terminal velocity in this drag-free model. Without mass accumulation, and still without drag, a drop released from rest has $v=gt$: its acceleration remains $g$, whereas the growing drop continually shares its momentum with newly acquired material.

The weak equivalence principle is universality of free fall for sufficiently small test bodies, independent of composition, when non-gravitational forces and self-gravity are negligible. Equality of inertial and passive gravitational mass makes this possible in the Newtonian limit. The Einstein equivalence principle adds local Lorentz invariance and local position invariance for nongravitational experiments: a small freely falling laboratory has the local physics of special relativity.

This motivates a universal spacetime metric. Freely falling unforced test particles follow geodesics, $u^{\nu }{\nabla }_{\nu }{u}^{\mu }=0$. At a point one can choose local inertial coordinates with vanishing connection, but curvature and tidal relative acceleration remain measurable over a finite region. Gravity cannot be transformed away globally. A uniformly accelerating small laboratory and a local gravitational field can have the same nongravitational experimental behaviour without having identical global geometry.

Clock comparisons predict a fractional rate difference of order $g\Delta h/c^2$ in a weak static field. The Pound-Rebka experiment tested the associated photon redshift through resonant gamma-ray emission and absorption at different heights, using controlled Doppler shifts to match the tiny energy displacement. Such a result tests the metric redshift prediction, not uniquely every component of the gravitational field equations.^[ref]

general relativity adds the dynamical Einstein field equations,The equivalence principle alone does not derive this unique equation; assumptions about covariance, dynamics and the Newtonian limit are also needed. The strong equivalence principle extends the universality claim to self-gravitating systems and gravitational experiments. Alternative metric theories can satisfy the weaker local principles while differing in gravitational dynamics, so composition tests, clock/redshift tests and self-gravity tests probe distinct claims.

Universality of free fall equates the passive gravitational response to inertial mass up to a universal choice of units. The Einstein equivalence principle further states that local nongravitational experiments in a freely falling laboratory obey special relativity independently of the laboratory's location and velocity. A supported laboratory instead has proper acceleration, just as an accelerating laboratory in flat spacetime does.

Around an event, freely falling coordinates can be chosen so that $g_{\mu \nu }={\eta }_{\mu \nu }$ and ${\Gamma }^{\alpha }{}_{\mu \nu }=0$ there. The first metric tensor derivatives vanish at that event, and nearby corrections begin at order curvature times squared separation. The metric tensor's second derivatives retain the Riemann tensor. Thus a local inertial frame removes Christoffel symbols at an event, while tidal gravity remains measurable over a finite region. Relative motion of separated freely falling test particles detects this curvature through geodesic deviation. Freely falling axes carried along an observer's world line extend the local construction, but do not make an extended curved region globally Minkowskian.

Local equivalence and tides. In Newtonian mechanics, inertial mass $m_{\mathrm{i}}$ measures the response to a force, whereas passive gravitational mass $m_{\mathrm{g}}$ couples to an imposed Newtonian gravitational field. A freely falling body obeysThe observed composition independence of free fall says that this ratio is universal. With a common choice of units it is one. This is the Newtonian content of the weak equivalence principle for test particles; it does not by itself specify all nongravitational physics.

The Einstein equivalence principle further says that local nongravitational experiments in a sufficiently small freely falling laboratory have the form predicted by special relativity, independently of the laboratory's location and motion. In such a local inertial frame, the metric tensor is locally the Minkowski metric and the Christoffel symbols vanish at the selected event. The strong equivalence principle extends the universality claim to local gravitational experiments and bodies with significant gravitational self-energy. None of these statements removes measurable tidal forces across a finite apparatus.

In an ideal small freely falling lift, a released ball initially floats relative to the cabin, because both undergo the same free fall. Two separated balls, however, need not retain their separation: they sample different parts of the gravitational field. For a point source, linearizing $\mathbf{g}=-GM\mathbf{r}/r^3$ across a small separation $\mathbf{\xi }$ gives the tidal tensor predictionA radial separation stretches with coefficient $2GM/r^3$, whereas a transverse separation is compressed with coefficient $-GM/r^3$. Initially parallel velocities therefore develop relative motion even though the individual bodies experience zero proper acceleration.

The relativistic statement is geodesic deviation. With the Riemann curvature tensor conventiona connecting vector between neighboring freely falling worldlines satisfiesHere $D/d\tau$ is the covariant derivative along the reference geodesic, and $U^{\mu }$ its four-velocity. Curvature can be nonzero where the Christoffel symbols have been made zero. In the weak, static, slow-motion limit this equation reduces to ${\ddot{\xi }}^i=-{\partial }_i{\partial }_j\Phi {\xi }^j$, recovering the tidal tensor above. A finite falling lift can detect gravity through these relative accelerations; local equivalence does not imply a globally uniform field.

Weak-field connection and potential. Take $x^0=ct$ and the $+---$ metric signature. There is an index error in the printed formula for the Christoffel symbols: its first differentiated metric component must be $g_{\nu 0}$, not $g_{\mu 0}$. The correctly contracted expression isThe offending term vanishes under the stipulated time independence, so this correction does not change the intended calculation. Write $g_{\mu \nu }={\eta }_{\mu \nu }+h_{\mu \nu }$ and treat both $h$ and its derivatives as first order. To that order $g^{00}=1-h_{00}$, $g^{0i}=O(h)$, and $g^{ij}=-{\delta }^{ij}+O(h)$. A stationary metric tensor givesFor the temporal component the inverse-metric coefficient is already first order, soFor a spatial component the zeroth-order inverse metric contributes:The time component of the approximate geodesic equation says that $dt/d\tau$ is constant to the retained order. In its spatial component, convert the proper time derivatives to coordinate-time derivatives and cancel $(dt/d\tau {)}^2$:Comparison with the Newtonian equation $d^2\mathbf{x}/d{t}^2=-\nabla \Phi$ identifies $h_{00}=2\Phi /c^2+C$. For an isolated source, choose the gravitational potential to tend to zero at infinity and normalize the time coordinate there, giving $C=0$. ThusThe quantity $\Phi$ is the Newtonian gravitational potential, or gravitational potential energy per unit mass, in $\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}={\mathrm{m}}^2{\mathrm{s}}^{-2}$. For a point source, $\Phi =-GM/r$. A static clock has $d\tau =\sqrt{g_{00}}dt\approx (1+\Phi /c^2)dt$, connecting the Newtonian limit with gravitational time dilation.

Frequency of the orbiting source. Treat the gas blob as a test particle on a timelike circular Schwarzschild geodesic, with pressure and other nongravitational forces negligible. Treat our very distant location as infinity, so the observer's own gravitational potential is negligible. These are the assumptions that make a unique answer in terms of $R_S/r$ possible.

Put $f=1-R_S/r$ and use a local static observer at the emission point as an intermediate measuring frame. The remote observer lies on the symmetry axis of the orbit. A ray arriving on that axis has zero conserved axial angular momentum, so its local direction is perpendicular to the emitter's azimuthal velocity. Bending within a meridional plane does not alter this orthogonality. The transverse relativistic Doppler effect therefore giveswhere $v$ is the orbital speed measured by that local static observer. The gravitational redshift between static Schwarzschild observers then givesA gravitational circular orbit has $v^2=GM/r+O((GM/r{)}^2/c^2)$ in the weak-field limit, so $v^2/c^2=R_S/(2r)+O((R_S/r{)}^2)$. Expanding both factors,Solving for the emitted frequency givesThe gravitational part is $R_S/(2r)$ and the transverse kinematic time dilation contributes another $R_S/(4r)$; the absence of line-of-sight motion does not remove this second effect.

As an exact check, the local speed on a timelike circular Schwarzschild geodesic satisfies $v^2/c^2=\mu /(r-2\mu )$, where $\mu =GM/c^2=R_S/2$. Therefore the face-on frequency shift of a circular Schwarzschild source iswhose expansion agrees with the boxed answer. If the gas is supported or accelerated by pressure or magnetic stresses rather than following a geodesic, its orbital speed is not fixed by the mass and radius alone; the appropriate first-order result is instead ${\nu }_{\mathrm{em}}\approx {\nu }_{\mathrm{ob}}[1+GM/(r{c}^2)+v^2/(2c^2)]$, with $v$ supplied independently.

For a source point mass $M$ at the origin, let $\mathbf{R}=R\mathbf{n}$, where $\mathbf{n}$ is an outward unit vector, and let $\mu =GM$ be the source's standard gravitational parameter. Differentiating $g_i=-\mu R_i/R^3$ gives the tidal tensorwhere ${\delta }_{ij}$ is the Kronecker delta and $I$ is the identity matrix. In an orthonormal basis with first axis radial, the tidal tensor isThe radial eigenvalue $2\mu /R^3$ separates nearby freely falling test particles; the two transverse eigenvalues $-\mu /R^3$ bring them together. The farther radial test particle falls inward less rapidly, so its acceleration relative to the nearer one points outward. Transverse gravitational accelerations converge toward the source. Thus an initially resting small cloud stretches radially and squeezes in both transverse directions. These are instantaneous relative accelerations; the subsequent shape also depends on initial velocities and the changing tidal tensor along the cloud's free fall.

Outside the source, the matrix trace is $2\mu /R^3-\mu /R^3-\mu /R^3=0$. This matches $\mathrm{tr}T=-{\nabla }^2\Phi =0$ in a region with zero mass density, by the Poisson equation for Newtonian gravity. Zero matrix trace means no leading change of a small initially resting cloud's volume at order $t^2$, not that each separation stays fixed. Inside matter, the same sign convention gives $\mathrm{tr}T=-4\pi G\rho$ instead.

In a fixed central Newtonian gravitational potential, pressureless free fall has speed $u=a{R}^{-1/2}$, and randomly moving stars have typical gas-relative speed $V=b{R}^{-1/2}$. A steady, undepleted spherical mass flux $A$ gives $\rho =A{R}^{-3/2}/(4\pi a)$. The Bondi--Hoyle--Lyttleton accretion rate therefore becomes $\dot{M}=\lambda M^2$ with $\lambda =G^{2}A/(a{b}^3)$, independent of radius. Integrating $d(1/M)/dt=-\lambda$ gives the displayed solution before $\lambda M_{0}t=1$. Its finite-time divergence signals failure of the fixed potential, undepleted gas supply and negligible feedback assumptions; it is not a physically unlimited mass supply.

A dipole's magnetic pressure scales as $B_{\ast }^2{R}_{\ast }^6/({\mu }_0{r}^6)$. Spherical free fall with $\rho =\dot{M}/(4\pi r^{2}v)$ and $v\propto (GM/r{)}^{1/2}$ has ram pressure proportional to $\dot{M}(GM{)}^{1/2}{r}^{-5/2}$. Equating them gives the seventh-root scaling. Disc geometry changes the coefficient; the estimate ceases to describe an exterior magnetosphere if the resulting radius lies inside the star.

For steady spherical isentropic flow with $p=K{\rho }^{3/2}$ around a point mass, put $j=r^2\rho v>0$ and $\mathcal{E}=v^2/2+2c^2-GM/r>0$. Eliminating the mass density gives $c^2=[j(3K/2{)}^2{]}^{2/5}{r}^{-4/5}{\mathcal{M}}^{-2/5}$. With $r_s=GM/(4\mathcal{E})$, $x=(r/r_s{)}^{1/5}$, $y={\mathcal{M}}^{2/5}$ and $\lambda =[j(3K/2{)}^2{]}^{2/5}/(2\mathcal{E}{r}_s^{4/5})$, the Bernoulli equation reduces to $\lambda F(y)=F(x)$, where $F(t)=t^4+4/t$. Since $F$ has its unique minimum at $t=1$, a regular sonic point requires $\lambda =1$ and $r=r_s$. One transonic branch is $y=x$, which has $v=\sqrt{2\mathcal{E}}$, $\rho \propto r^{-2}$ and $p\propto r^{-3}$. The decreasing branch has inner free fall and an outer nearly static reservoir, giving Bondi accretion for inward fluid flow. A nonradiative stationary shock wave keeps $j,\mathcal{E},r_s$ fixed but raises $K$, so $\lambda$ increases by the factor $(K_2/K_1{)}^{4/5}$. The inward transonic solution is the spherical $n=2$, $\gamma =3/2$ specialization of transonic accretion in a power-law tube.

First published composition-dependent free fall comparison from MICROSCOPE, using electrostatically controlled orbital test masses. The abstract gives the measured differential-acceleration parameter and its statistical and systematic uncertainties.

A timelike geodesic has timelike tangent and vanishing covariant acceleration in an affine parameter. Choosing proper time gives $g(\dot{\gamma },\dot{\gamma })=-1$.

Physically, an ideal massive test particle in free fall follows a timelike geodesic when nongravitational forces and appreciable self-gravity can be neglected. The normalization above uses $c=1$; in SI units it is $g(u,u)=-c^2$. Zero covariant acceleration does not mean constant coordinate velocity: a dropped object accelerates relative to the ground. Conversely, a person standing on the ground is pushed away from free fall by the supporting surface. A photon follows a null geodesic, whose parameter cannot be proper time. See Tong's geodesics in spacetime.

An ideal fluid issuing from a small hole a vertical depth $h$ below a free surface has speed

Here $h$ is the vertical distance from the free surface to the centre of the orifice, and $g$ is the positive magnitude of gravitational acceleration. The ideal prediction assumes incompressible flow, the same ambient pressure at the free surface and emerging jet, negligible velocity of the reservoir's free surface, and negligible viscous dissipation and surface tension effects. For a falling free surface, a quasi-steady approximation additionally assumes that the local outflow adjusts much faster than $h$ changes. It applies the Bernoulli equation at each instant, rather than claiming that the whole drainage is steady flow.^[ref]

A short proof uses a streamline connecting the reservoir's free surface at elevation $z_s$ to the emerging jet at elevation $z_o$. With constant mass density $\rho$, negligible surface speed and no mechanical energy losses, the Bernoulli equation givesEqual pressures cancel, leaving $v^2/2=gh$: a unit mass gains the same kinetic energy as a body in free fall through height $h$. This makes Torricelli's law a useful first estimate before the detailed velocity field is known. At one quarter of the original head, ideal jet speed is one half as large; it is not proportional to head.

Far from the accelerating outlet region, the reservoir has approximately hydrostatic pressure $p_{\mathrm{ambient}}+\rho gh$ at the outlet's elevation. The emerging jet instead has approximately ambient pressure when surface tension is negligible. Do not assign that bulk hydrostatic pressure to the emerging jet as well as adding the height difference: that would count the available mechanical energy twice. A sealed reservoir, a submerged outlet or significant surface tension changes the pressure conditions and requires a different balance.

Actual volumetric flow rate is commonly represented by a discharge coefficient:where $a$ is the geometric area of the orifice. This definition compares volumetric flow rates, not just speeds. If a nearly uniform jet at its vena contracta has area $a_j$ and speed $v_j$, define the two separate ratios $C_c=a_j/a$ and $C_v=v_j/\sqrt{2gh}$. Then $Q=a_j{v}_j$ gives $C_d=C_c{C}_v$. A contracted jet can consequently have a speed close to $\sqrt{2gh}$ while $C_d$ is appreciably below one. No particular $C_d$ is universal: calibrate it for the opening and operating range.