Airfoil 2026-10-13
An airfoil (British spelling: aerofoil) is the cross-sectional profile of a lifting surface, such as a wing or blade. Its chord is the straight segment from leading edge to trailing edge, with length ; its mean camber line runs midway between its upper and lower surfaces, and its thickness measures their separation. Camber is the departure of that mean line from the chord. The angle of attack describes the profile's orientation relative to the incoming flow. These independent choices let a designer alter the distribution of pressure, lift and drag, rather than simply enlarging the surface. [ref]
https://www1.grc.nasa.gov/beginners-guide-to-aeronautics/wing-geometry/. NASA Glenn, Wing Geometry, Airfoil Geometry subsection
The distinction between an airfoil and a finite wing matters. A two-dimensional section model assumes no variation along span and predicts lift per unit span . A wing of span has a spanwise distribution and total lift . Its tip flow and trailing vortices generate a downward induced velocity, called downwash, which changes the local flow seen by the sections. Thus one cannot generally multiply the isolated two-dimensional by . The wing aspect ratio measures its span relative to its reference area and helps characterize these three-dimensional effects. [ref]
https://pressbooks.lib.vt.edu/aerodynamics/chapter/chapter-3-additional-aerodynamic-tools/. James F. Marchman III, Aerodynamics and Aircraft Performance, chapter 3.4: finite wings
For a thin section with small camber slopes and small angle of attack, thin airfoil theory predictswhere both angles are in radians and is the zero-lift angle of attack. The model assumes two-dimensional, steady flow, incompressible flow, inviscid flow, an attached outer flow, and smooth departure at a sharp trailing edge (the Kutta condition). This last condition selects the circulation used by the Kutta–Joukowski theorem; with counterclockwise circulation positive and , that theorem gives for positive lift towards . The beauty is that a detailed velocity field collapses to a simple force law, while the geometry enters through . [ref]
https://ocw.mit.edu/courses/16-01-unified-engineering-i-ii-iii-iv-fall-2005-spring-2006/d721171c42af48a056aaccec784a6d10_f03_0304.pdf. MIT Unified Engineering, F3 lecture notes, pp. 2–4: camber and lift slope
[ref]https://ocw.mit.edu/courses/16-01-unified-engineering-i-ii-iii-iv-fall-2005-spring-2006/cc3f20c40c9451bde6a595fdd96a7e9e_f02_mud_0304.pdf. MIT Unified Engineering, F2 questions 7 and 9: symmetric sections and Kutta condition
A symmetric airfoil has in this model: at the predicted lift coefficient is , and reversing reverses lift. A stipulated cambered section with instead predicts at zero angle of attack. These are ideal calculations, not measurements or guarantees for a particular profile. Air parcels on the upper and lower surfaces have no obligation to arrive together at the trailing edge; the pressure and velocity field must satisfy the flow problem. A symmetric lifting section already makes an unequal-path-length explanation inadequate. [ref]
https://www1.grc.nasa.gov/beginners-guide-to-aeronautics/bernoulli-and-newton/. NASA Glenn, Bernoulli and Newton: the equal-transit explanation is incorrect
Real airfoils have a viscous boundary layer, even when inviscid flow approximates the outer motion. At sufficiently demanding incidence, flow separation can become extensive and the lift coefficient departs from its approximately linear law, reaches a maximum and may fall: this is aerodynamic stall. Lift does not generally vanish at aerodynamic stall. The onset depends on profile, Reynolds number, surface condition and motion history; predicts neither a universal stall angle nor post-stall behavior. The pressure/balance experiment in lift makes this limit observable by measuring changes in the surface pressure distribution alongside the total force. [ref]
https://pressbooks.lib.vt.edu/aerodynamics/chapter/chapter-3-additional-aerodynamic-tools/. James F. Marchman, chapter 3.3: separation invalidates thin-airfoil assumptions
Bernoulli's principle 2026-10-14
The Bernoulli principle relates pressure, kinetic energy per unit volume and Newtonian gravitational potential along a streamline of steady flow with inviscid flow. For constant mass density, the Bernoulli equation gives constant along that streamline. Constancy across different streamlines requires further hypotheses, such as irrotational flow.
Precisely, suppose a fluid of constant positive mass density obeys the Euler equations for an inviscid fluid, its velocity field and pressure are steady, and the body force per unit mass is . ThenEvery term has units of energy per unit mass. For a uniform gravitational field, Newtonian gravitational potential is , up to an additive constant. The pressure term represents work by surrounding fluid, so acceleration can exchange pressure work with kinetic energy. Dotting the steady Euler equations with gives the displayed directional derivative: this is why constancy follows a streamline. [ref]
https://www.damtp.cam.ac.uk/user/mem2/FLUIDS-IB/dyn.pdf. Cambridge, Fluid Dynamics IB, Dynamics, section 2.4
Continue the horizontal water constriction example from incompressible flow. Assume negligible viscous losses, steady flow and approximately uniform flow at both measurement sections, so the streamline speeds there equal the mean speeds and . Equal height cancels the Newtonian gravitational potential terms. The ideal pressure drop isThis is the Venturi effect in its ideal limit: the gain in kinetic energy density is supplied by pressure work. Cambridge's treatment uses this relation to connect a measured pressure difference to volumetric flow rate. [ref]
https://www.damtp.cam.ac.uk/user/mem2/FLUIDS-IB/dyn.pdf. Cambridge, Fluid Dynamics IB, Dynamics, section 2.5
The result is conditional, rather than a prediction of the actual viscous pipe. Mass conservation fixes and the ratio of mean speeds even with viscosity, but nonuniform velocity profiles require an integrated energy balance, and wall friction dissipates mechanical energy. A pump would supply additional work. In particular, constant mean speed in a straight pipe does not imply constant pressure: Poiseuille flow needs a continuing pressure gradient to overcome viscosity. [ref]
https://uclnatsci.github.io/2021/Electromagnetism-Fluids-and-Waves/fluids/bernoulli.html. UCL, Bernoulli's equation, assumptions and boundary limitations
Choking by wall thickness in an elastic tube Created 2026-10-14 Updated 2026-11-05
For steady uniform-density inviscid flow with pressure law and a prescribed thickness , mass conservation and the Bernoulli equation give , where and . The function has a unique maximum at , producing the displayed choking bound. Below it there are two positive roots; continuity from the upstream state selects the lower-speed root for and higher-speed root for . At , the maximum is zero and any positive prescribed thickness is incompatible with this model.
D'Alembert's paradox 2026-10-13
Ideal steady flow with incompressible flow, inviscid flow and irrotational flow can predict zero drag on a body that has nonzero drag in a real fluid. The zero-circulation cylinder pressure coefficient distribution gives an explicit cancellation example; viscosity and flow separation are absent from that model.
Displacement thickness 2026-10-14
For constant-mass density incompressible flow, the displacement thickness measures the lost volume flux per unit span divided by the outer velocity:An outer inviscid flow encountering the same volume flux deficit behaves, at leading order, as though the wall were displaced outward by . In a Blasius boundary layer, numerical integration gives . This equivalent displacement is smaller than the boundary layer thickness and is not the position of a physical surface.[ref]
Dynamic pressure 2026-10-14
Dynamic pressure is the bulk kinetic energy per unit volume associated with fluid motion at speed in a specified reference frame:Here is mass density. The word "pressure" identifies its units and its use as a pressure scale; is not an extra isotropic stress to add to the local static pressure in every surface traction. A stationary fluid can have large static pressure with . Conversely, knowing and fixes without fixing static pressure. The definition remains useful in viscous fluid flow and compressible flow, even when the simplest Bernoulli equation is inapplicable. [ref]
For stipulated mass density and reference speed , . At with the same mass density, . This quadratic scaling makes dynamic pressure a useful normalization for the lift coefficient, drag coefficient and pressure coefficient; it predicts four times a particular force only if its corresponding coefficient remains unchanged.
The connection with stagnation pressure requires hypotheses. For steady flow, constant mass density, inviscid flow, and gravity as the only body force, the Bernoulli equation gives, along one streamline,Here is static pressure, is gravitational acceleration, and is height. Bring that flow to rest without losses at the same height and the limiting static pressure is the stagnation pressure :There must be no intervening pump supplying work or other energy input. Constancy along one streamline does not assert equal on all streamlines; irrotational flow in a connected region supplies that additional conclusion. Changing height requires retaining . In compressible flow the stopping pressure difference generally differs from , and viscous dissipation can reduce recovered stagnation pressure. [ref]
A Pitot-static tube turns this relation into a measurement: its upstream-facing opening senses stagnation pressure, while separate side openings sample static pressure. Connect these to the positive and negative ports of a differential pressure transducer or manometer. Under the hypotheses above, the measured difference givesThus the stipulated and imply , not an independently measured speed. Align the probe with the incoming flow and keep its passages clear: blockage and pinched tubing invalidate the comparison. At very low speed, the difference can be smaller than instrument error; at high Mach number, compressibility and possibly a shock wave require a different reduction. [ref]
A useful control is differential-pressure calibration. For a digital manometer, zeroing establishes the reading at equal port pressures; it does not establish the gain away from zero. The manufacturer's Series 475 Handheld Digital Manometer: Specifications - Installation and Operating Instructions by Dwyer Instruments, Inc. directs the user to vent both ports to the atmosphere before zeroing. Its stated accuracy is a fraction of full scale, not of the reading, and depends on instrument temperature. Select a suitable range and record the instrument variant and units; the manual covers several ranges.[ref]
Series 475 Handheld Digital Manometer: Specifications - Installation and Operating Instructions by Dwyer Instruments, Inc.. p. 1: Specifications, Model Chart and ZERO instructions
An independent physical check can use a liquid-column manometer: apply a stable pressure difference to both instruments and compare the digital reading with the hydrostatic pressure inferred from the vertical level difference. For equal-height connections carrying a gas of mass density above a manometric liquid of mass density ,This follows by balancing hydrostatic pressure down and up the two connected columns. Use vertical height, even if a tube is inclined. Known mass densities, measured height and gravitational acceleration set the comparison scale. Check several differences across the working range and reverse the port connections to check the sign. This is a proposed independent check, not a claim that the cited teaching experiment performed a full instrument calibration. Recheck zero after reconnecting the sealed tubing; a zero adjustment alone cannot detect a blocked passage or an incorrect scale.
Gravity-capillary wave Created 2026-10-14 Updated 2026-10-20
A gravity-capillary wave has restoring contributions from both gravity and surface tension. Consider an initially resting, horizontally unbounded liquid layer of constant depth , with uniform mass density , incompressible flow, inviscid flow and irrotational flow. The upper gas exerts constant pressure and has negligible inertia. Assume a clean interface with constant surface tension : no surfactant stresses, imposed current or wind. The wave amplitude must satisfy and , so linearized boundary conditions apply. Deep water means , rather than a particular depth in metres. [ref]
https://www.ocw.mit.edu/courses/18-354j-nonlinear-dynamics-ii-continuum-systems-spring-2015/6163837c5efa5d6ed969e535ec8f890c_MIT18_354JS15_lectureNotes.pdf. Brenner, Peacock, Stocker, Reis and Dunkel, MIT Nonlinear Dynamics II, §23.1, pp. 105–106
Let be the undisturbed free surface, with liquid below it. Potential flow gives a velocity potential with , and incompressible flow gives the Laplace equation. The linearized kinematic boundary condition for a free-surface graph and the Unsteady Bernoulli equation, including gravity and the Young–Laplace equation, giveHere is the interface displacement; the signed small-slope pressure jump is . A wave crest has , so its capillary pressure is positive. Choose a right-travelling Fourier mode of positive wavenumber and angular frequency :The factor decays into the deep liquid, satisfies the Laplace equation, and makes at . Substituting in the dynamic boundary condition leavesThis derivation displays the two restoring contributions with the same sign. Dimensional analysis also checks the result: has units , so both and have units . [ref]
https://www.ocw.mit.edu/courses/18-354j-nonlinear-dynamics-ii-continuum-systems-spring-2015/6163837c5efa5d6ed969e535ec8f890c_MIT18_354JS15_lectureNotes.pdf. MIT Nonlinear Dynamics II, §23.1, equations (547)–(551): boundary conditions and dispersion
At finite depth, a flat impermeable bottom at requires . Replacing by the normalized vertical factor makes the surface ratio . Eliminating the mode amplitudes then givesThus the same model supplies its own depth check: use the deep-water formula only when is sufficiently close to one for the required accuracy. For example, reduces by about relative to deep water, while reduces it by about . These percentages are calculations from this factor. [ref]
In deep water the phase velocity magnitude follows a wave crest, while the group velocity magnitude follows the wave envelope of a narrow wave packet. Differentiating the positive-frequency dispersion relation givesFor , the surface gravity wave limit has . For , the deep-water capillary-wave dispersion relation has : a wave crest slips backward relative to the faster wave envelope, even while both advance. Neither speed is the oscillating fluid element's own velocity. [ref]
There is a particularly neat coincidence at the crossover. Differentiate : its unique minimum is at , exactly where the restoring terms are equal. ThereforeFor the stipulated parameters in capillary length, the minimum is . At , direct substitution gives , and . This short wavelength is capillary dominated, but its group velocity is still slightly below the pure-capillary ratio because gravity has not vanished.
Incompressible flow Created 2026-09-24 Updated 2026-10-13
An incompressible flow preserves the volume of each moving fluid element. Its local condition is zero velocity divergence, . Combined with mass conservation, it says that each parcel's mass density stays unchanged as it moves. This does not require every parcel to have the same mass density: an initially nonuniform distribution can be transported without compression. Nor does incompressibility require a constant velocity or constant pressure.
For positive mass density, the continuity equation makes the equivalence precise:The material derivative follows a fluid element. Thus if and only if the velocity field is a divergence-free vector field. The distinction between following a parcel and observing a fixed position prevents confusing incompressibility with an unchanging density field. [ref]
https://uclnatsci.github.io/2021/Electromagnetism-Fluids-and-Waves/fluids/conservation.html. UCL, Conservation law derivations, sections 30.1–30.2
For a rigid pipe with impermeable walls, no leaks and steady flow, mass conservation at constant mass density implies a constant volumetric flow rate:Here points downstream, is the section's area and is the mean normal velocity, not necessarily the centreline speed. No assumption of inviscid flow or a uniform velocity profile is needed. If unequal amounts entered and left a fixed segment, its mass would accumulate, contradicting steady flow.
For a worked example, take horizontal water flow with stipulated mass density , inlet radius , outlet radius and inlet mean velocity . The circular areas giveThe radius halves, the area falls to one quarter, and the mean velocity quadruples. The mass flow rate is . The same amount of water passes both sections each second; narrowing changes its speed rather than compressing it. The Reynolds number and Bernoulli's principle discussions use this same example to examine viscous effects and the additional assumptions needed to calculate pressure.
Inviscid flow Created 2026-10-06 Updated 2026-11-05
An inviscid flow neglects viscosity in the fluid flow. Its momentum balance is described by the Euler equations for an inviscid fluid.
Irrotational circulation in an exterior domain Created 2026-10-14 Updated 2026-11-05
In smooth barotropic inviscid flow with conservative body force, Kelvin's circulation theorem preserves circulation of material curves. Initially zero vorticity therefore gives zero circulation on contractible loops for all times. A loop winding around an excluded obstacle can have a nonzero constant circulation: the planar field is irrotational outside the origin but has circulation on a positively oriented enclosing loop.
Kutta–Joukowski theorem Created 2026-09-29 Updated 2026-11-05
For steady two-dimensional inviscid flow of density , far-field velocity , and circulation , the force per unit span is perpendicular to the flow. With counterclockwise circulation taken positive and , it is .
Following a fluid element. Let be a pathline, satisfying . The material derivative, also called the Lagrangian time derivative, is the rate observed while following that element. The multivariable chain rule givesThe partial derivative holds position fixed; the advection term accounts for motion through a spatially varying field. For a vector field, the same formula applies component by component in fixed Cartesian axes.
Momentum balance. For an inviscid flow, pressure is the only surface stress. A small material volume has net pressure force , obtained by the divergence theorem. Its mass is and its material acceleration is . Applying Newton's second law gives the Euler equations for an inviscid fluid,This is the local momentum conservation equation when viscosity and physical body forces are absent. A prescribed body force per unit mass would add on the right.
Changing to rotating coordinates. In a rotating frame with constant angular velocity , the lab acceleration is the rotating-frame material acceleration plus . ThereforeFor uniform rotation,With constant , this gives the Euler equation in a uniformly rotating frame,The term subtracted from is the equilibrium centrifugal contribution. Dividing the bracket by defines a reduced pressure , with units . The Coriolis acceleration remains a velocity-dependent term.
Linear disturbances. Linearize about rest in the rotating coordinates, with the equilibrium part of removed. The advective acceleration is quadratic in disturbance amplitude, so the linearized Euler equations and incompressible flow condition areTaking the curl removes the pressure gradient. The identitygives the inertial-wave vorticity equation,
Plane-wave frequencies. Let and use the given plane wave convention. Incompressible flow requires , and the vorticity equation givesFor a nonzero frequency, apply again, using . It follows thatThis is the inertial-wave dispersion relation. If is the angle of the wavevector to the rotation axis, the positive angular frequency is . It is largest, , for a parallel wavevector, decreases to zero for a perpendicular one, and is independent of the magnitude of the wavenumber at fixed direction. The two signs describe the two oscillatory polarization branches. When , the zero-frequency result is obtained directly from the linear equations rather than by dividing by ; it describes the nonoscillating balanced limit.
The first argument mixes a material fluid element with a fixed spatial volume. For a material volume element following the fluid, its volume is generally time dependent:Therefore the material derivative cannot be taken outside the product with as in the proposed cancellation. Instead, mass conservation statesThe momentum balance for a material fluid element, using Newton's second law, is consequentlyFor a fixed control volume, fluid crosses its boundary and carries momentum; a momentum flux term must be included. Its local conservative form is , equivalent to the displayed Euler equations after using the continuity equation. The erroneous equation omits the term needed to convert a material derivative of momentum density into this conservative form.
For the second task, take a homogeneous fluid with constant mass density , prescribed gravitational acceleration , incompressible flow and inviscid flow. Its equations areUnder the replacement , the same velocity solves these equations if , up to a spatially constant pressure offset. This is density scaling invariance of incompressible Euler flow:The comparison assumes the same geometry and compatible velocity or gravitational free surface conditions; an independently specified nonuniform boundary pressure would supply additional information. The gravitational field is taken as prescribed. If the fluid's own self-gravity must be included, its field changes with mass density, so the invariance need not hold. A single freely moving particle's density-independent gravitational acceleration does not justify neglecting that distinction.
Take positive circulation to mean counterclockwise in the diagram, with right and up. The upper surface is then traversed right to left. To first order in ,Thus the bound circulation is clockwise. Any exterior contour enclosing the hydrofoil has the same circulation in an irrotational flow; a globally completed two-dimensional flow must contain a circulation component decaying as at large distance. The simple local upper/lower velocity prescription cannot describe the entire exterior. For a start from rest, Kelvin's circulation theorem requires compensating shed vorticity, represented by a starting vortex in a physical wake.
Numerically,If the question's fluid force includes buoyancy, and the total upward fluid force is about , using . Weight of the hydrofoil is not specified and is not part of this fluid force. The motion-induced lift also agrees with the Kutta–Joukowski theorem, in this orientation.
The chord-based Reynolds number isAt this large Reynolds number, smooth laminar flow over the whole body is unlikely in ordinary conditions; wall layers may undergo transition and separate. A thin-layer inviscid flow approximation outside them can remain useful, but does not justify globally laminar flow.
Past exam of the fluids and waves course of the University of Oxford 2017 2 a ii Solution 2026-10-31
For smooth inviscid flow of constant mass density under a conservative body force, Kelvin's circulation theorem states that the circulation around any closed material curve is conserved:The constant-mass density assumption is the usual homogeneous incompressible flow setting here. Merely requiring in a fluid of varying mass density would not by itself make the pressure term below a gradient.
Parameterize the material curve by a periodic particle label , writing . Differentiate its line integral, keeping fixed:The second term vanishes by periodicity. If the conservative body force per unit mass is , the Euler momentum equation givesIts line integral around a closed curve is zero, proving the claim. The velocity field must remain smooth on the moving curve, which cannot pass through a line vortex singularity during this argument. No assumption of irrotational flow was needed.
The denominator of the travelling response vanishes atFor the positive speeds in the paper, a finite travelling response exists whenever . At , the dynamic condition would reduce to , which is impossible. The finite-depth dispersion relation of a surface gravity wave with wavenumber is . At , its phase speed equals : the moving pressure pattern drives a natural surface gravity wave at its own speed, producing resonance.
The restriction concerns a bounded travelling response. To distinguish this from existence of a solution of an initial value problem, set , , and . Ifthe linearized free-surface boundary conditions reduce to the forced harmonic oscillatorFor example, if the forcing is switched on at from a flat surface with zero initial velocity, then, for ,This satisfies and contains both free oscillations and the forced component. At a solution of the initial value problem still exists, but its particular solution containswhich grows linearly with time. Thus the ideal inviscid flow model predicts resonant growth rather than a finite constant amplitude. Additional free oscillations can set the initial conditions but cannot remove that secular growth.
For one-dimensional inviscid flow, conservation of mass and the Euler equation areSubstitute and the stated background-plus-perturbation expressions. Dropping products of perturbations gives the linearizationThe Taylor expansion of the equation of state givesSince the background coefficients are constant, commutes with . Applying to the linearized continuity equation and using the linearized Euler equation therefore yieldsThis is the convected acoustic wave equation: in the coordinate , it becomes the ordinary one-dimensional wave equation. Positive is required for propagating sound waves.
Take the wind to blow toward positive , so . In subsonic flow, both stationary observers record the source angular frequency and, for the imposed density boundary data, the same acoustic density perturbation amplitude . Their spatial records distinguish the wavelengths and directions:The downstream pattern moves at and the upstream pattern at signed velocity . Thus the downstream observer sees longer wavelengths while the upstream observer sees shorter ones. The interval containing multiple wavelengths lets each measure the spatial spacing and its direction of travel. The arrival delays at are and respectively; after arrival their temporal traces differ by these propagation phases, not by temporal frequency.
In supersonic flow, both characteristic curves have positive slopes . With an initially undisturbed upstream region, the observer at sees no source-generated disturbance. The observer at can see two downstream travelling waves, withThe faster signal arrives first, at ; the slower one arrives at . Once both are present, their different wavenumbers can produce spatial modulation, while the time dependence at a fixed position remains harmonic with angular frequency .
There is a qualification to the source data: the single prescribed value of fixes only . It does not determine the two complex amplitudes separately. A second source condition, such as its downstream acoustic particle velocity with amplitude , or a specified mass/force source model, is needed. Consequently two downstream branches are allowed; the stated density data alone do not require both to have nonzero amplitude. Extending a formally upstream-going subsonic formula into when would impose a signal arriving from upstream infinity and violate the causal source interpretation.
As , one laboratory phase velocity tends to zero. Its wavelength tends to zero, its wavenumber diverges, and its arrival time diverges; a finite-time experiment never attains the assumed established upstream field at fixed . Finite source size, viscosity, thermal losses and eventually molecular scales invalidate this singular inviscid flow limit. With effective acoustic diffusivity , a loss estimate is , which diverges. Moreover the finite-amplitude change in the characteristic speed is of order , so uniform linearization requires . At exactly sonic background speed, the zero-speed branch cannot support a nonzero- harmonic solution. Nonlinear acoustics and losses must resolve the near-sonic regime; an unlimited spatial compression of the ideal wave is not physical.
The coefficient of the growing term is fixed by the far-field complex potential: . On , the no-penetration boundary condition makes the cylinder a streamline, so the imaginary part of must be constant. Choosing the two terms to be complex conjugates makes that constant zero, givingThis is the zero-circulation solution, also obtained by the Milne-Thomson circle theorem. Impermeability and the leading far field alone permit an additional circulation term ; the requested algebraic form implicitly selects . No no-slip boundary condition is imposed on an inviscid flow.
The streamfunction in polar coordinates isIt gives andThe liquid still exerts normal pressure on each surface element even though tangential viscous traction is absent. For this steady flow with zero vorticity, the Bernoulli equation has one constant throughout the connected fluid region:There is no finite uniform far-field pressure to substitute for , because this straining velocity field grows with radius. Taking outward from the cylinder, the fluid force per unit axial width isIndeed while : diametrically opposite elements cancel exactly. Zero resultant force does not mean zero local pressure or zero deformation tendency.
Choose the stream function convention , , so that this is an incoming stagnation point flow rather than the reversed flow. The velocity field isAt the plate , satisfying the no-penetration boundary condition. The nonzero tangential velocity is permitted for an inviscid flow; the no-slip boundary condition is not imposed on the Euler equations for an inviscid fluid.
The steady convective acceleration is . Thus the Euler equations for an inviscid fluid requireIntegrating givesThe printed formula omits from its quadratic term. The source uses physical pressure, not pressure divided by mass density, as is clear from its term and its Navier-Stokes equation. The omitted factor is also required by dimensional analysis. For a concrete substitution check, the printed expression gives rather than the required when in the chosen numerical units.

