Editing Drag equation

{{Short description|Equation for the force of drag}}
{{Use mdy dates |date=June 2022}}

In [[fluid dynamics]], the '''drag equation''' is a formula used to calculate the force of [[drag (physics)|drag]] experienced by an object due to movement through a fully enclosing [[fluid]]. The equation is:
<math display="block">F_{\rm d}\, =\, \tfrac12\, \rho\, u^2\, c_{\rm d}\, A</math>
where
*<math>F_{\rm d}</math> is the drag [[force]], which is by definition the force component in the direction of the flow velocity,
*<math>\rho</math> is the [[mass density]] of the fluid,<ref>For the [[Earth's atmosphere]], the air density can be found using the [[barometric formula]]. Air is {{cvt|1.293|kg/m3|lb/cuft}} at {{cvt|0|C|F|0}} and 1  [[atmosphere (unit)|atmosphere]]</ref>
*<math>u</math> is the [[flow velocity]] relative to the object,
*<math>A</math> is the reference [[area]], and
*<math>c_{\rm d}</math> is the [[drag coefficient]] – a [[dimensionless number|dimensionless]] [[physical coefficient|coefficient]] related to the object's geometry and taking into account both [[skin friction]] and [[form drag]]. If the fluid is a liquid, <math>c_{\rm d}</math> depends on the [[Reynolds number]]; if the fluid is a gas, <math>c_{\rm d}</math> depends on both the Reynolds number and the [[Mach number]].

The equation is attributed to [[Lord Rayleigh]], who originally used ''L''<sup>2</sup> in place of ''A'' (with ''L'' being some linear dimension).<ref>See Section 7 of Book 2 of Newton's [[Principia Mathematica]]; in particular Proposition 37.</ref>

The reference area ''A'' is typically defined as the area of the [[orthographic projection]] of the object on a plane perpendicular to the direction of motion. For non-hollow objects with simple shape, such as a sphere, this is exactly the same as the maximal [[cross section (geometry)|cross sectional]] area. For other objects (for instance, a rolling tube or the body of a cyclist), ''A'' may be significantly larger than the area of any cross section along any plane perpendicular to the direction of motion. [[Airfoils]] use the square of the [[chord (aircraft)|chord length]] as the reference area; since airfoil chords are usually defined with a length of 1, the reference area is also 1. Aircraft use the wing area (or rotor-blade area) as the reference area, which makes for an easy comparison to [[lift (force)|lift]]. [[Airship]]s and [[Solid of revolution|bodies of revolution]] use the volumetric coefficient of drag, in which the reference area is the square of the [[cube root]] of the airship's volume. Sometimes different reference areas are given for the same object in which case a drag coefficient corresponding to each of these different areas must be given.

For sharp-cornered [[bluff body|bluff bodies]], like square cylinders and plates held transverse to the flow direction, this equation is applicable with the drag coefficient as a constant value when the [[Reynolds number]] is greater than 1000.<ref>[http://www.ac.wwu.edu/~vawter/PhysicsNet/Topics/Dynamics/Forces/DragForce.html Drag Force<!-- Bot generated title -->] {{webarchive |url=https://web.archive.org/web/20080414225930/http://www.ac.wwu.edu/~vawter/PhysicsNet/Topics/Dynamics/Forces/DragForce.html |date=April 14, 2008 }}</ref> For smooth bodies, like a cylinder, the drag coefficient may vary significantly until Reynolds numbers up to 10<sup>7</sup> (ten million).<ref>See Batchelor (1967), p. 341.</ref>

== Discussion ==
The equation is easier understood for the idealized situation where all of the fluid impinges on the reference area and comes to a complete stop, building up [[Pressure#Stagnation pressure|stagnation pressure]] over the whole area. No real object exactly corresponds to this behavior. <math>c_{\rm d}</math> is the ratio of drag for any real object to that of the ideal object. In practice a rough un-streamlined body (a bluff body) will have a <math>c_{\rm d}</math> around 1, more or less. Smoother objects can have much lower values of <math>c_{\rm d}</math>. The equation is precise – it simply provides the definition of <math>c_{\rm d}</math> ([[drag coefficient]]), which varies with the [[Reynolds number]] and is found by experiment.

Of particular importance is the <math>u^2</math> dependence on flow velocity, meaning that fluid drag increases with the square of flow velocity. When flow velocity is doubled, for example, not only does the fluid strike with twice the flow velocity, but twice the [[mass]] of fluid strikes per second. Therefore, the change of [[momentum]] per time, i.e. the [[force]] experienced, is multiplied by four. This is in contrast with solid-on-solid [[Friction#Kinetic friction|dynamic friction]], which generally has very little velocity dependence.

== Relation with dynamic pressure ==
The drag force can also be specified as
<math display="block">F_{\rm d} \propto P_{\rm D} A</math>
where ''P''<sub>D</sub> is the pressure exerted by the fluid on area ''A''. Here the pressure ''P''<sub>D</sub> is referred to as [[dynamic pressure]] due to the [[kinetic energy]] of the fluid experiencing relative flow velocity ''u''. This is defined in similar form as the kinetic energy equation:
<math display="block">P_{\rm D} = \frac12 \rho u^2</math>

== Derivation == <!-- [[Drag (physics)]] links here -->
The '''drag equation''' may be derived to within a multiplicative constant by the method of [[dimensional analysis]]. If a moving fluid meets an object, it exerts a force on the object. Suppose that the fluid is a liquid, and the variables involved – under some conditions – are the:
* speed ''u'', 
* fluid density ''ρ'',
* [[Viscosity#Dynamic and kinematic viscosity|kinematic viscosity]] ''ν'' of the fluid, 
* size of the body, expressed in terms of its wetted area ''A'', and 
* drag force ''F''<sub>d</sub>.
Using the algorithm of the [[Buckingham π theorem]], these five variables can be reduced to two dimensionless groups: 
* [[drag coefficient]] c<sub>d</sub> and
* [[Reynolds number]] Re.

That this is so becomes apparent when the drag force ''F''<sub>d</sub> is expressed as part of a function of the other variables in the problem:

<math display="block">
f_a(F_{\rm d}, u, A, \rho, \nu) = 0.
</math>

This rather odd form of expression is used because it does not assume a one-to-one relationship.  Here, ''f<sub>a</sub>'' is some (as-yet-unknown) function that takes five arguments. Now the right-hand side is zero in any system of units; so it should  be possible to express the relationship described by ''f<sub>a</sub>'' in terms of only dimensionless groups.

There are many ways of combining the five arguments of ''f<sub>a</sub>'' to form dimensionless groups, but the Buckingham π theorem states that there will be two such groups. The most appropriate are the Reynolds number, given by

<math display="block">
  \mathrm{Re} = \frac{u\sqrt{A}}{\nu}
</math>

and the drag coefficient, given by

<math display="block">
  c_{\rm d} = \frac{F_{\rm d}}{\frac12 \rho A u^2}.
</math>

Thus the function of five variables may be replaced by another function of only two variables:

<math display="block">
  f_b\left(\frac{F_{\rm d}}{\frac12 \rho A u^2}, \frac{u \sqrt{A}}{\nu} \right) = 0.
</math>

where ''f<sub>b</sub>'' is some function of two arguments.
The original law is then reduced to a law involving only these two numbers.

Because the only unknown in the above equation is the drag force ''F''<sub>d</sub>, it is possible to express it as

<math display="block">\begin{align}
  \frac{F_{\rm d}}{\frac12 \rho A u^2} &= f_c\left(\frac{u \sqrt{A}}{\nu} \right) \\
  F_{\rm d} &= \tfrac12 \rho A u^2 f_c(\mathrm{Re}) \\
  c_{\rm d} &= f_c(\mathrm{Re})
\end{align}</math>

Thus the force is simply {{sfrac|1|2}} ''ρ'' ''A'' ''u<sup>2</sup>'' times some (as-yet-unknown) function ''f<sub>c</sub>'' of the Reynolds number Re – a considerably simpler system than the original five-argument function given above.

Dimensional analysis thus makes a very complex problem (trying to determine the behavior of a function of five variables) a much simpler one: the determination of the drag as a function of only one variable, the Reynolds number.

If the fluid is a gas, certain properties of the gas influence the drag and those properties must also be taken into account. Those properties are conventionally considered to be the absolute temperature of the gas, and the ratio of its specific heats. These two properties determine the speed of sound in the gas at its given temperature. The Buckingham pi theorem then leads to a third dimensionless group, the ratio of the [[relative velocity]] to the speed of sound, which is known as the [[Mach number]]. Consequently when a body is moving relative to a gas, the drag coefficient varies with the Mach number and the Reynolds number.

The analysis also gives other information for free, so to speak.  The analysis shows that, other things being equal, the drag force will be proportional to the density of the fluid.  This kind of information often proves to be extremely valuable, especially in the early stages of a research project.

==Air viscosity in a rotating sphere==
Air viscosity in a rotating sphere has a coefficient, similar to the drag coefficient in the drag equation.<ref>{{cite tech report
 | first1=C.
 | last1=Cena
 | first2=L.F.S.
 | last2=Guizado
 | first3=J.V.B.
 | last3=Ferreira
 | doi=10.1590/1806-9126-RBEF-2023-0351
 | title=Determination of skin friction on a rotating sphere in magnetic levitation
 | publisher=SciELO
 | location=US
 | date=2 February 2024
}}</ref>

==Experimental methods==
{{unreferenced section|date=January 2022}}
To empirically determine the Reynolds number dependence, instead of experimenting on a large body with fast-flowing fluids (such as real-size airplanes in [[wind tunnel]]s), one may just as well experiment using a small model in a flow of higher velocity because these two systems deliver [[Similitude (model)|similitude]] by having the same Reynolds number. If the same Reynolds number and Mach number cannot be achieved just by using a flow of higher velocity it may be advantageous to use a fluid of greater density or lower viscosity.<ref name="Zohuri2015">{{cite book |last=Zohuri |first=Bahman |title=Dimensional Analysis and Self-Similarity Methods for Engineers and Scientists |year=2015 |publisher=Springer, Cham |doi=10.1007/978-3-319-13476-5_2 |chapter=Similitude Theory and Applications|pages=93–193 |isbn=978-3-319-13475-8 }}</ref>

== See also ==
* [[Aerodynamic drag]]
* [[Angle of attack]]
* [[Morison equation]]
* [[Newton's sine-square law of air resistance]]
* [[Stall (flight)]]
* [[Terminal velocity]]

==References==
{{Reflist}}

== External links ==
*{{cite book
 | first=G.K.
 | last=Batchelor
 | authorlink=George Batchelor
 | title=An Introduction to Fluid Dynamics
 | year=1967
 | publisher=Cambridge University Press
 | isbn=0-521-66396-2 }} 
*{{cite book
 | last = Huntley | first = H. E.
 | year = 1967
 | title = Dimensional Analysis
 | publisher = Dover
 | id = LOC 67-17978
 }}
*{{cite web
 | first=Tom
 | last=Benson
 | url=https://www.grc.nasa.gov/www/k-12/VirtualAero/BottleRocket/airplane/falling.html
 | title=Falling Object with Air Resistance
 | publisher=NASA
 | location=US
 | access-date=2022-06-09
 }}
*{{cite web
 | first=Tom
 | last=Benson
 | url=https://www.grc.nasa.gov/www/k-12/rocket/drageq.html
 | title=The drag equation
 | publisher=NASA
 | location=US
 | access-date=2022-06-09
 }}

[[Category:Drag (physics)]]
[[Category:Equations of fluid dynamics]]
[[Category:Aircraft wing design]]