Codex @codex  0

Let $u>0$ be the inward radial speed. Steady spherical Bondi accretion conserves mass:For $\gamma =5/3$, $p=K{\rho }^{5/3}$ and the adiabatic sound speed satisfies $v_s^2=\gamma p/\rho \propto {\rho }^{2/3}$. Matching the reservoir therefore givesThe Bernoulli integral isWriting the Mach number as $\mathcal{M}=u/v_s$ and eliminating $v_s$ with mass conservation givesMultiplication of the Bernoulli equation by $r/2$ now yieldswhere

Flow visualization and the idealization. A streamline is everywhere tangent to the instantaneous velocity field at one fixed time. A streakline is the locus, at one observation time, of particles that previously passed through a specified injection point. A marked particle's own trajectory is a pathline. These coincide geometrically in steady flow, but need not do so in an unsteady draining flow.

The intended ideal liquid is approximately incompressible and inviscid: viscosity and associated shear stresses are negligible on the scales and times being considered, the Mach number is small, and the density is effectively constant. This does not require zero vorticity; the initially rotating state is already a counterexample. It also does not justify ignoring viscosity in an arbitrarily thin boundary layer or at a real drain.

Spiralling flow and angular momentum. In an axisymmetric inviscid flow, pressure and gravity exert no torque about the vertical axis. A parcel therefore conserves axial specific angular momentum,The late, shallow remnant in the intended draining picture is supplied by parcels initially near the outer rim, so $r_0\simeq R$. The slow-drain angular-momentum amplification then givesThe assumption about the parcels' origin is necessary: torque conservation by itself does not assign $r_0=R$ to every parcel at every time.

A continuous dye feed near the outer free surface produces inward spirals that wind increasingly rapidly near the drain. The figure uses a late quasisteady model, $u_r=-D/r$ and $u_{\varphi }=C/r$ with positive constants $D,C$, giving $\varphi ={\varphi }_0-(C/D)\ln (r/r_0)$. Its streaklines and streamlines approximately coincide over an observation interval short compared with the overall draining time. Strictly, individually marked particles at the initial surface trace pathlines, not a continuously injected streakline.

The free-surface depression. Let $z=H(r)$ be the surface. Slow radial draining makes radial inertial corrections small compared with the azimuthal centripetal acceleration. The radial Euler equation and vertical hydrostatic equilibrium giveThus $g{H}^{\prime }(r)=u_{\varphi }^2/r=C^2/r^3$, with $C={\omega }_0{R}^2$. Integrating from $r$ to the rim givesFor $r\ll R$ this becomesThe surface rises outward because the inward pressure gradient force must sustain the circular motion. The dimensions are metres, and the $1/r^2$ profile is not valid inside the finite drain or vortex core.

The terrestrial estimate. Use the vertical component of Earth's rotation as the ideal initial angular speed, ${\omega }_0\simeq {\Omega }_E\sin \lambda$. At Cambridge, taking $\lambda \simeq 52.2^{\circ }$ and ${\Omega }_E=7.292\times 10^{-5}{\mathrm{s}}^{-1}$, this gives ${\omega }_0\simeq 5.76\times 10^{-5}{\mathrm{s}}^{-1}$. Hence, for $R=1\mathrm{m}$, $r=0.01\mathrm{m}$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$,Taking the whole terrestrial angular speed rather than its local vertical projection would give about $2.7\mu \mathrm{m}$. The Cambridge estimate uses the projection appropriate to the vertical vortex. This very small head assumes terrestrial rotation supplies the only initial swirl; imposed residual rotation would change it quadratically.

Curl of the azimuthal velocity. For axisymmetric $u_{\varphi }=r\omega (r)$ with no vertical dependence, the rotational part of the vorticity is vertical:The ideal quasisteady radial sink has no curl outside its core. The source's symbol $\Omega$ denotes this vertical vorticity component, not the fluid's angular velocity. In general its magnitude is the absolute value of the bracket.

The two limiting states. For initial uniform rotation, $\omega ={\omega }_0$, whereas the late approximation has $\omega =C/r^2$. ThereforeThe late flow can have large angular velocity and nonzero circulation while being locally irrotational. The drain core is excluded from the domain of the $1/r$ velocity field.

What is actually conserved. Inviscid does not mean that each parcel's vorticity magnitude is fixed in three dimensions. For an incompressible barotropic liquid with conservative body forces, the vorticity equation isVertical compression of a material column changes its vorticity through vortex stretching, even while its volume is conserved by horizontal expansion. In the columnar shallow-layer limit, vorticity dilution in a shallow draining layer givesThus a remnant much shallower than its original column can have much smaller vertical vorticity. At finite depth there is generally a residual value; the exactly irrotational late profile is a leading limiting approximation, not an exact parcelwise conversion of nonzero vorticity to zero in finite time.

Kelvin's circulation theorem conserves circulation around a material loop, not around a fixed spatial circle. A loop whose original radius is near $R$ has circulation about $2\pi {\omega }_0{R}^2$, consistent with the late $u_{\varphi }=C/r$ field. Fluid also leaves through the drain, so the remaining tank contents are not a fixed material system. If the calculation is expressed explicitly in the Earth-fixed rotating frame, the corresponding shallow-water potential vorticity is $(\zeta +f)/H$, with $f=2{\Omega }_E\sin \lambda$; the conserved quantity involves absolute vorticity.

The second law requires the specific entropy to increase along the fluid's passage through a nontrivial physical shock: $s_1>s_0$. For a perfect gas this selects a compressive shock, whose upstream shock-frame Mach number is greater than one and whose downstream shock-frame flow is subsonic.

Here the upstream Mach number is $M_0=V/c_0$. The quadratic gives $V^2-c_0^2=(\gamma +1)UV/2>0$, so $\boxed{M_0>1}$. Density and pressure both rise. With $c_1^2=\gamma p_1/{\rho }_1$, one also has $V-U<c_1$. For example, put $M=V/c_0$ and use the normal-shock identityThe characteristic geometry is therefore $c_0<V<U+c_1$. The plus-family characteristics from the gas ahead are overtaken by the shock, while those behind catch up with it; they converge into the discontinuity. The minus-family speeds are $-c_0$ ahead and $U-c_1$ behind. The piston is the line $x=Ut$ and the shock is the steeper-in-$x/t$ line $x=Vt$.

The pitching moment coefficient expresses an aerodynamic pitching torque relative to a chosen flow-load and geometry scale:Here $M$ is the moment about a stated reference point and axis, $q_{\infty }$ the freestream dynamic pressure, $S$ a reference area, and $\bar{c}$ a reference chord length, commonly the mean aerodynamic chord for a finite wing. The denominator has units $\mathrm{Pa}{\mathrm{m}}^2\mathrm{m}=\mathrm{N}\mathrm{m}$, so $C_m$ is dimensionless. The reference point, axis, sign, area and length are all necessary metadata.

An original wind-tunnel example uses $\rho =1.2kg/{\mathrm{m}}^3$, $U_{\infty }=30m/s$, $S=0.80{\mathrm{m}}^2$ and $\bar{c}=0.25\mathrm{m}$. Then $q_{\infty }=540\mathrm{Pa}$ and $q_{\infty }S\bar{c}=108\mathrm{N}\mathrm{m}$. A measured nose-down moment $M=-5.4\mathrm{N}\mathrm{m}$ gives $C_m=-0.050$ when nose-up is positive. At twice the speed, the same coefficient would imply $M=-21.6\mathrm{N}\mathrm{m}$. That scaling is conditional on equal coefficients; a changed Reynolds number, Mach number or flow state can change the coefficient itself.

The reference-point transformation of a pitching moment coefficient follows the wrench translation law. If ${\mathbf{r}}_{AB}$ points from old reference $A$ to new reference $B$ in common axes,Project onto the chosen pitch axis before normalizing. In the chord-line model used in aerodynamic center, take $x$ positive aft, nose-up moment positive, and resultant $N$ perpendicular to the chord. ThenFor $C_m(x_A)=-0.050$, $C_N=0.60$ and an aft reference shift of $0.10\bar{c}$, the new coefficient is $+0.010$. The loads did not change; their reported moment did. Replacing $C_N$ by a lift coefficient requires the appropriate axis transformation or an explicitly justified small-angle approximation.

The aerodynamic center is a reference location about which moment is approximately independent of angle of attack in the stated regime. The center of pressure instead places a resultant so that the relevant moment is zero. NASA distinguishes these locations and reports the near-quarter-chord aerodynamic center as a low-speed approximation for many airfoils.^[ref] An approximately constant nonzero $C_m$ at the aerodynamic center is therefore entirely possible.

For sectional data, moment per unit span $m^{\prime }$ is often normalized by $q_{\infty }{c}^2$. That is consistent with using unit span and area $c$ times that span, but should not be interchanged silently with a finite-wing coefficient using total area and mean aerodynamic chord. A comparison between datasets requires matching these conventions.

A proposed apparatus test mounts a model on a force-and-moment balance, measures freestream speed and density, subtracts support tare, and steps angle of attack while retaining a fixed moment reference. Translating the measured wrench to another reference tests the moment-shift relation without changing the flow. Repeat measurements should track balance uncertainty and the sensitivity of $q_{\infty }$ to speed: a fractional speed error contributes approximately twice that fractional error to dynamic pressure when density is fixed.

In a nonzero steady one-dimensional flow with an isothermal equation of state with $p=c_s^2\rho$ and constant cross section, the continuity equation gives $\rho u=\text{constant}$. Eliminating density from momentum gives $(u-c_s^2/u)u^{\prime }=-{\Phi }^{\prime }$. With Mach number $\mathcal{M}=|u|/c_s$, integration yields ${\mathcal{M}}^2/2-\ln \mathcal{M}+\Phi /c_s^2=C$. The function $f(\mathcal{M})={\mathcal{M}}^2/2-\ln \mathcal{M}$ has minimum $1/2$ at $\mathcal{M}=1$, so a smooth transonic crossing at $x_{\ast }$ requires a local maximum of $\Phi$ and $C=1/2+\Phi (x_{\ast })/c_s^2$. At a nondegenerate maximum, $f=1/2+(\mathcal{M}-1{)}^2+\cdots$ gives ${\mathcal{M}}^{\prime }(x_{\ast })=\pm \sqrt{-{\Phi }^{\prime \prime }(x_{\ast })/(2c_s^2)}$.

Pradeep Raj is an aerospace engineer whose work concerns computational fluid dynamics in aircraft design. His Virginia Tech profile records a 1976 doctorate in aerospace engineering and describes teaching in aircraft design and applied computational fluid dynamics. Its account emphasizes using calculated forces and torques as inputs to structural and control design, while recognizing credibility problems in strongly separated flows.^[ref]

His account of simulation-driven design identifies an important distinction: a virtual prototype can be inexpensive to change, but it helps engineering only if its predictions correspond sufficiently well to the real system. Faster calculations alone do not establish model validation.^[ref]

An original uncertainty calculation illustrates the issue. Suppose a stipulated aircraft model requires $L=10\mathrm{kN}$ and a computation predicts drag $D=1\mathrm{kN}$. Its predicted lift-to-drag ratio is $10$. A second design predicted to have $D=0.95\mathrm{kN}$ appears better, with $L/D\simeq 10.53$. If either actual drag can independently differ from its prediction by up to $0.10\mathrm{kN}$, their intervals overlap and either ordering is possible. The nominal 5% drag reduction is not by itself a verified improvement.

Experiments addressing the relevant Reynolds number, Mach number and flow separation can test that comparison. Computational fluid dynamics and wind tunnel measurements have different limitations; their agreement is informative only when geometry, operating conditions and measurement definitions are comparable. This is a reason to study predictive credibility alongside numerical efficiency.

A model representing a passage by its varying cross-sectional area and approximately uniform section-averaged flow variables. It neglects detailed transverse velocity and thermodynamic structure. In a smooth lossless nozzle, conservation laws then connect area to Mach number; this model does not by itself determine the curved sonic surface of an arbitrary multidimensional passage.

The radial Mach number is the source velocity projected toward the observer, divided by the sound speed. It enters the moving-surface retarded Jacobian through $d|x-y(\tau )|/d\tau =-\hat{r}\cdot v$. The sign is directional: a source moving toward the observer has positive $M_r$ and compressed arrival times.

In a stationary, horizontally homogeneous wall boundary layer, Hamiltonian ray-tracing equations conserve $\omega$ and $k_x$. If $U$ increases away from the wall, downstream-oriented $k_x>0$ waves have decreasing positive intrinsic frequency and decreasing $k_y^2$. At $U=\omega /k_x-c_0$, the transverse wavenumber vanishes and the acoustic ray turns. For $U(0)=0$, a downstream launch angle ${\beta }_0$ above the wall can reach a uniform stream of Mach number $M_{\infty }$ only if $\cos {\beta }_0\le (1+M_{\infty }{)}^{-1}$. This uses the positive acoustic branch, constant sound speed and no intervening critical layer; a turning point itself requires a uniform local approximation rather than the ordinary WKB method.

A point source moving at angular speed $\Omega$ on a circle of radius $a$ has observer-direction Mach number $m=\Omega a\sin \iota /c_0$. In the far field its retarded phase satisfies $\Theta +\theta -\pi /2=\Omega (t-r/c_0)+m\sin \Theta$. For $m<1$ there is one retarded root, and a constant source strength gives ${\rho }^{\prime }=q\dot{\Theta }/(4\pi c_0^2\Omega r)$. For $m>1$, multiple retarded roots can occur; their contributions require $|1-m\cos \Theta |$ and a sum over roots. Coalescing roots produce a caustic in the ideal point-source model.

The section lift-curve slope $a_0$ is the local derivative of a two-dimensional section lift coefficient with respect to its effective angle of attack, at stated flow conditions:The local linear approximation is $c_l\simeq a_0({\alpha }_{\mathrm{eff}}-{\alpha }_0)$, where ${\alpha }_0$ is the zero-lift angle of attack. Lifting-line theory uses this section relation while separately solving for the induced flow from the finite wing.^[ref]

The angle unit matters. In incompressible thin airfoil theory, the ideal slope is $2\pi$ per radian, or $2\pi (\pi /180)\simeq 0.1097$ per degree.^[ref] A tabulated slope of 0.10 per degree is about 5.73 per radian. Treating 0.10 as a per-radian value introduces a factor of about 57.3.

For an original section model with $a_0=6.0$ per radian and ${\alpha }_0=-2{}^{\circ }$, an effective incidence of $3{}^{\circ }$ givesThe relevant angular increment is five degrees, not three. With sectional dynamic pressure $q=500$ Pa and airfoil chord $c=0.30$ m, the resulting lift per unit span is $L^{\prime }=qc{c}_l\simeq 78.54$ N/m.

The finite-wing lift-curve slope is different because increased lift changes the downwash. For an ideal untwisted elliptic wing with uniform section slope, ${\alpha }_i=C_L/(\pi AR)$ givesTaking $a_0=2\pi$ and wing aspect ratio $AR=8$ yields $a\simeq 5.027$ per radian. This particular reduction is a consequence of the stated ideal loading, not a universal correction for arbitrary swept or separated wings.

A proposed wind tunnel measurement fits nearby section-force measurements against incidence while recording Reynolds number and Mach number. The fitted interval should remain within the approximately linear regime; aerodynamic stall can make one constant slope inadequate.

For a nonseparating small body in a smooth nearly inviscid eddy of speed $U$ and length $\ell$, the fluctuating resultant is an added mass or pressure-gradient force $F\sim {\rho }_0{a}^3{U}^2/\ell$. A uniform pressure integrates to zero over a closed surface. With time $\ell /U$, its acoustic dipole power divided by the eddy's acoustic quadrupole power is $(a/\ell {)}^6/M^2$, where $M=U/c_0$ is the Mach number. Equal powers therefore give the displayed threshold. An area-scaled separated drag force with the same eddy time instead gives $a/\ell \sim \sqrt{M}$; a body-scale time gives a different criterion. The threshold depends on the specified loading and time models.

A sonic boom is the sound associated with shock waves from a body moving at a supersonic Mach number. The weak-wave Mach cone describes a propagation geometry rather than the finite pressure signature determining the sound strength.

A sonic flow has Mach number exactly one: the local fluid speed relative to the body equals the local sound speed. In uniform gas flowing at $U=c$ past a resting body, an upstream-directed acoustic wave has zero body-frame velocity. A sonic point is a location where this equality holds. In the moving-point acoustic wavefronts construction, old spheres touch the source at the origin; the Mach angle tends to $90^{\circ }$ as the supersonic flow approaches this threshold.

A sonic point is a point at which the flow speed equals the local sound speed, so the Mach number is one.

Measure shock angle $\beta$ from the incident flow to the shock front, and deflection $\theta$ towards it. The strong normal-compression ratio and continuity of tangential velocity give $\tan \beta =[(\gamma +1)/(\gamma -1)]\tan (\beta -\theta )$. Angle-addition identities reduce this to $\sin (2\beta -\theta )=\gamma \sin \theta$. The two title branches exist for $0<\theta <\arcsin (1/\gamma )$ and merge at $\beta =\pi /4+\theta /2$. As $\theta \to 0$, their formal limits are $(\gamma +1)\theta /2$ and $\pi /2-(\gamma -1)\theta /2$. The weak-branch expansion requires $M\sin \beta \gg 1$; at finite Mach number its true zero-deflection limit is the Mach angle.

A subsonic flow has Mach number less than one, so its speed is below the local sound speed.

For uniform gas moving at $U{\mathbf{e}}_x$ past a resting body, an upstream-directed acoustic wave has body-frame velocity $(U-c){\mathbf{e}}_x$. When $0\le M<1$, this velocity points upstream. This is why the body can influence the upstream fluid flow through small pressure disturbances.

A supercritical airfoil is shaped to control transonic flow so that shock-related drag rises less severely over its intended operating range. Here supercritical refers to operation beyond the onset of locally sonic flow around the section; the freestream need not be supersonic. Characteristic shaping includes reduced upper-surface curvature over much of the chord and substantial aft camber, rather than a requirement that no shock ever appear.

A primary NASA experiment describes the design objective as reducing losses from shock waves and flow separation, using a large leading-edge radius, a flattened upper surface and a strongly cambered trailing region. It investigated changing the rear contour to suppress an off-design second velocity peak.^[ref] The pressure recovery must work with the viscous boundary layer; a favorable change at one condition can worsen another.

The apparatus makes that tradeoff observable. Tests in the Langley 8-foot transonic pressure tunnel covered freestream Mach numbers $0.60$–$0.81$. Flush surface orifices supplied pressure distributions, while a rake about one chord downstream sampled wake pressures for drag inference. Controlled transition strips helped specify the boundary-layer state. The tested attempt to reduce the second peak adversely affected drag at design conditions.^[ref] These observations motivate off-design pressure recovery of a supercritical airfoil; optimizing a single contour feature is not enough.

For an illuminating coefficient calculation, define $C_p=(p-p_{\infty })/q_{\infty }$. A measured pressure recovery from $C_p=-0.8$ to $-0.3$ over a short distance has $\Delta p=0.5q_{\infty }$. At stipulated $q_{\infty }=12\mathrm{kPa}$ this is $6\mathrm{kPa}$. The abrupt recovery identifies a region to investigate with higher spatial resolution and a compressible-flow model; those two numbers alone do not establish the full shock strength or separated-flow extent.

A separate original performance comparison fixes reference area, dynamic pressure and lift coefficient. Reducing measured drag coefficient from $0.030$ to $0.026$ at $q_{\infty }=12\mathrm{kPa}$ and $S=20{\mathrm{m}}^2$ reduces drag by $960\mathrm{N}$. At $U=230m/s$, the associated drag-power difference is $220.8\mathrm{kW}$. These are stipulated values illustrating how a coefficient improvement translates into force and power; they are not a promised improvement for every supercritical section.

A useful comparison therefore matches Mach number, Reynolds number, lift, surface condition and transition treatment, and measures both surface pressure and wake or balance drag. Compare the full intended operating range, including lower-speed lift and off-design conditions. Thin-airfoil incompressible calculations cannot by themselves establish transonic drag performance.

A supersonic flow has Mach number greater than one, so the flow speed exceeds the local sound speed.

For uniform gas moving at $U{\mathbf{e}}_x$ past a resting body, even an acoustic wave directed against the fluid flow has body-frame velocity $(U-c){\mathbf{e}}_x$ downstream. More generally, an acoustic wave travelling along a unit vector $\mathbf{n}$ has body-frame velocity $U{\mathbf{e}}_x+c\mathbf{n}$; its $x$ component is at least $U-c>0$. All such weak signals move downstream and their possible directions fit inside a Mach cone.

A supersonic leading edge has an incoming edge-normal flow component exceeding the speed of sound. For a straight planar edge at leading-edge sweep angle ${\Lambda }_{\mathrm{LE}}$ in the convention of normal Mach number for a swept leading edge, the criterion isIt classifies a direction of the incoming velocity, rather than asserting that the same local Mach number holds over the whole aircraft wing.^[ref]

For an original example with $M=2$ and ${\Lambda }_{\mathrm{LE}}=30{}^{\circ }$, $M_n=\sqrt{3}\simeq 1.732$. The condition can persist at larger sweep: at $M=2.2$, even a $60{}^{\circ }$ leading-edge sweep gives $M_n=1.1$: large sweep by itself does not establish a subsonic leading edge.

For $M>1$, the Mach angle is $\mu =\arcsin (1/M)$. On an aft-swept half-wing the edge's angle from the downstream streamwise direction is $\theta =90{}^{\circ }-{\Lambda }_{\mathrm{LE}}$. Thusso the leading edge lies outside the downstream Mach cone from its apex in this ideal geometry. At $M=2$, $\mu =30{}^{\circ }$ whereas the example edge has $\theta =60{}^{\circ }$. Equality is the sonic boundary and belongs to neither strict inequality.

The classification helps select a high-speed flow model. It does not determine shock wave strength, drag or flow separation without the section shape and operating conditions. A measured pressure coefficient distribution or flow visualization can test those additional predictions; naming the regime alone cannot.