circulation and vorticity


0.1 Introduction

Vortex theory is essentially due to Sir W. Thomson, Lord Kelvin [1] and H. v. Helmholtz [3], although dโ€™Alembert, Euler, Cauchy, Lagrange, Hankel , Hadamard and Stokes also contibuted important ideas. A useful kinematic notion in many problems of hydrodynamics is that of circulation and vorticity. These concepts are usually related through the Stokesโ€™ theorem. There are also many purposes for which it is more convenient to think in terms of circulation and vorticity rather than in terms of the velocity field, despite the simpler physical character of the latter quantity. It also proves to be possible and useful, in many important cases of fluid flow, to separate the flow field into two regions with different properties, one of them being characterized by the vorticity and being approximately zero everywhere. It is well-known that such concepts are associated with viscous fluids 11It turns out that inside the boundary layer viscous effects take place, and therefore the flow is rotational (๐Ž=โˆ‡ร—๐ฏโ‰ ๐ŸŽ) whereas outside of it the flow is irrotational, i.e. the flow is potential(๐Ž=โˆ‡ร—๐ฏโ‰ก๐ŸŽโŸน๐ฏ=โˆ‡โกฯ•, where the scalar function ฯ•=ฯ•โข(๐ฑ,t) is the potential of the velocity field) but we will only treat the purely kinematical consequences.

0.2 Circulation

Let us consider a contractible closed curveMathworldPlanetmath embedded in the fluid flow that we shall call a circuitMathworldPlanetmath๐’ž and described in a counterclockwise sense. Each arc element can be considered as an infinitesimalMathworldPlanetmathPlanetmath vector dโข๐ฑ which is tangentPlanetmathPlanetmathPlanetmath to ๐’ž. As usual, ๐ฏ=๐ฏโข(๐ฑ,t) represents the instantaneous velocity field at each point. Then, if the scalar productMathworldPlanetmath of ๐ฏ and dโข๐ฑ is integrated around the circuit, the line integral

ฮ“=โˆฎ๐’ž๐ฏโ‹…๐‘‘๐ฑ=โˆฎ๐’žโˆฅ๐ฏโˆฅโขcosโก(๐ฏ,dโข๐ฑ)โข๐‘‘x, (1)

where dโขx=โˆฅdโข๐ฑโˆฅ, is called the circulation around ๐’ž. All values of ๐ฏ are to be taken for the same instant t. Although evident because the additiveness of integral (1) defining circulation, results illustrative to show that circulation is additive. Suppose we cut the circuit ๐’ž by some path AโขB and give the two new circuits AโขBโข๐’ž1โขA with circulation ฮ“1 and BโขAโข๐’ž2โขB with circulation ฮ“2, both the same sense as ๐’ž. Thus (1) shows that

ฮ“1+ฮ“2=โˆฎAโขBโข๐’ž1๐ฏโ‹…๐‘‘๐ฑ+โˆฎBโขAโข๐’ž2๐ฏโ‹…๐‘‘๐ฑ=โˆซAโขB๐ฏโ‹…๐‘‘๐ฑ+โˆซ๐’ž1๐ฏโ‹…๐‘‘๐ฑ+โˆซBโขA๐ฏโ‹…๐‘‘๐ฑ+โˆซ๐’ž2๐ฏโ‹…๐‘‘๐ฑ,

cancelling the integrals along AโขB and BโขA since ๐ฏ is the same and dโข๐ฑ is opposite on these two paths. Thus,

ฮ“1+ฮ“2=โˆซ๐’ž1๐ฏโ‹…๐‘‘๐ฑ+โˆซ๐’ž2๐ฏโ‹…๐‘‘๐ฑโ‰กโˆฎ๐’ž๐ฏโ‹…๐‘‘๐ฑ,

exactly the inegral around ๐’ž, and hence

ฮ“=ฮ“1+ฮ“2. (2)

We may generalize (2). Let ๐’œ be any open two-sided surface spanning ๐’ž and consider the side where the sense of description of ๐’ž appears counterclockwise, and normals to ๐’œ will always be drawn out from this side. On ๐’œ draw two sets of orthogonal curves forming a latticeMathworldPlanetmath and each generated mesh is a closed circuit, the sense of description being taken counterclockwise as viewed from the normal to the surface; its circulation denoted by ฮ“i,i=1,โ‹ฏ,m. By iterated application of (2) we obtain

ฮ“=โˆ‘i=1mฮ“i, (3)

where m is the number of meshes in the lattice. Increasing the number of orthogonal curves in such a way that the lattice becomes more dense and all meshes become smaller, increasing so the number of terms in (3). Let a function ฮณ=ฮณโข(๐ฑ,t) be defined at each point P of ๐’œ as the limit of the quotient between the circulation along the contour of a mesh around P and the area of the mesh, becoming all of these one steadily smaller in all directions. So that for an arbitrary mesh the circulation is ฮ“iโ‰ˆฮณiโขdโข๐”ži, where ฮณi is the value of ฮณ at some point in the mesh, and dโข๐”ži its respective area. Thus, as the number of terms increases indefinitely, the right-hand member of (3) yields the surface integral of ฮณ over ๐’œ, i.e.

ฮ“=โˆซ๐’œฮณโข๐‘‘๐”ž. (4)

From the definition of ฮณ, we conclude that the value of this function at any point ๐ฑ of ๐’œ depends upon the distribution of the velocity ๐ฏ=๐ฏโข(๐ฑ,t) in the neighborhoodMathworldPlanetmath of this point and that ฮณ=โˆ‡ร—๐ฏโ‹…๐ง, 22This can be easily seen if we choose a mesh centered in a neighborhood of ๐ฑ equipped with Cartesian local coordinates {xi}, and impose the sufficient surface regularity at that point. Thus the circulation, around that mesh located in the local plane {x1,x2}, is given by (comma denotes partial differentiation respect to the indicated componentPlanetmathPlanetmathPlanetmath) dโขฮ“=(v2,1-v1,2)โขdโขx1โขdโขx2. Since dโขx1โขdโขx2 is the area of this mesh, the function ฮณ must have the value (v2,1-v1,2) at ๐ฑ. But this quantity is exactly the x3-component of โˆ‡ร—๐ฏ, defined as โˆ‡ร—๐ฏ=(v3,2-v2,3,v1,3-v3,1,v1,2-v2,1). This definition is valid in any rectangular right-handed coordinate system. Since the x3-direction is here that of the normal at ๐ฑ to the surface ๐’œ, the result in question follows. where ๐ง is the outward normal to the surface ๐’œ at an arbitrary point ๐ฑ. From this fact and by using (1) and (4) we have

โˆฎ๐’ž๐ฏโ‹…๐‘‘๐ฑ=โˆซ๐’œโˆ‡ร—๐ฏโ‹…๐งโข๐‘‘๐”ž, (5)

which is known as Stokesโ€™ theorem. It states that the circulation along any circuit is given by the surface integral of โˆ‡ร—๐ฏ over any surface spanning the circuit. It is clear that (5) can be applied only if it is possible to find some surface that has the given circuit as rim and on which โˆ‡ร—๐ฏ is defined everywhere. Thus, in the case of a fluid flow circulating around an infiniteMathworldPlanetmath cylindrical obstacle, no such surface can be found for any circuit surrounds the cylinderMathworldPlanetmath. However, the Stokesโ€™ theorem can be yet applied if we choose two circuits ๐’ž1 and ๐’ž2 about the obstacle (no intersecting themselves) and by making a cut AโขB, it can be combined into a single circuit for which a suitable spanning surface exists. Then the left side of (5) becomes (ฮ“1-ฮ“2). 33The contributions from AโขB and BโขA cancel. In particular, if the flow is irrotational (a term introduced by Lord Kelvin) in its domain, i.e. โˆ‡ร—๐ฏโ‰ก๐ŸŽ, then ฮ“1=ฮ“2. So that the circulation is equal for every circuit surrounding the obstacle.

0.3 Vorticity

The analysisMathworldPlanetmath of the relative motion near a point of the fluid outcomes that the force exerted by one portion of fluid on an adjacentPlanetmathPlanetmathPlanetmath portion depends on the way in which the fluid is being deformed by the motion, and it is necessary as a preliminary to dynamical considerations, to make an analysis of the character of the motion in the neighborhood of any point. That analysis it has to do with the study of local rate of strain and rate of rotationMathworldPlanetmath. Thus, the velocity field of the fluid at the place ๐ฑ and time t is given by ๐ฏ=๐ฏโข(๐ฑ,t) and the simultaneous velocity at a neighboring position ๐ฑ+ฮดโข๐ฑ is ๐ฏ+ฮดโข๐ฏ. Thus, for Cartesian coordinatesMathworldPlanetmath,

ฮดโขvi=โˆ‚โกviโˆ‚โกxjโขฮดโขxj, (6)

where the usual summation index convention applies here and that equation is correct to the first order in the small distanceMathworldPlanetmath ฮดโข๐ฑ between the two points located in the cited neighborhood. The kinematical character of the relative velocity ฮดโข๐ฏ, considered as a linear functionMathworldPlanetmath of ฮดโข๐ฑ, can be recognized by decomposing the velocity gradientMathworldPlanetmath โˆ‚โกvi/โˆ‚โกxj, which is a tensor of second rank, into parts which are symmetrical and skew-symmetrical in the indices i and j. That is,

ฮดโขvi=ฮดโขvi(d)+ฮดโขvi(w),

where

ฮดโขvi(d)=diโขjโขฮดโขxj,ฮดโขvi(w)=wiโขjโขฮดโขxj,

and

diโขj=12โข(โˆ‚โกviโˆ‚โกxj+โˆ‚โกvjโˆ‚โกxi),wiโขj=12โข(โˆ‚โกviโˆ‚โกxj-โˆ‚โกvjโˆ‚โกxi).

The first above equation corresponds to the well-known rate of deformation tensor, but we are here interested in the second one so-called the vorticity tensor, i.e.

wiโขj:=12โข(xห™i,j-xห™j,i)=12โข(vi,j-vj,i), (7)

where xห™i,jโ‰กvi,j are the Cartesian components of the velocity gradient in connectionMathworldPlanetmath to the velocity field ๐ฏโข(๐ฑ,t)โ‰ก๐ฑห™โข(๐ฑ,t). Since (7) is skew-symmetric, the associate axial vectorMathworldPlanetmath in the Euclidean space โ„3 is given by (ฯตiโขjโขk is the Levi-Civitaโ€™s isotropic Cartesian tensor)

2โขฮฉkโ‰ก-ฯตiโขjโขkโขvi,j=ฯตiโขjโขkโขvj,i,

which was introduced by Lagrange and Cauchy [4] and was shown by Cauchy and Stokes [6] to represent a local instantaneous rate of rotation in a neighborhood of some point ๐ฑ in the fluid media, that in the time being we call it the local angular velocity in such neighborhood (a dynamical cause is generally due to the fluid viscosity). It is usually called vortex vector or simply vorticity, and is defined by

๐Ž:=โˆ‡ร—๐ฏ. (8)

Notice that vorticity is, by definition, the twice of the local angular velocity, that is, ๐Žโ‰ก2โข๐›€.

0.4 The vorticity distribution

One consequence from the definition of vortex vector is the identityPlanetmathPlanetmath

โˆ‡โ‹…๐Žโ‰ก0. (9)

A line in the fluid whose tangent is everywhere parallelMathworldPlanetmathPlanetmathPlanetmath to the local vortex vector is termed a vortex line. The family of such lines at any instant is defined by an equation analogous to the streamlines. The surface in the fluid flow formed for all the vortex lines passing through a given contractible closed curve drawn in the fluid is said to be a vortex tube. The flux of the vortex vector across an open surface bounded by this same closed curve and lying entirely in the fluid flow is

โˆซ๐’œ๐Žโ‹…๐งโข๐‘‘๐”ž,

i.e. the Stokesโ€™ theorem right-hand side. We can use (9) to prove that this integral has the same value for any open surface lying in the fluid flow and bounded by any contractible closed curve which lies in the vortex tube and passes round it once. For if ๐งโขdโข๐”ž and ๐งโ€ฒโขdโข๐”žโ€ฒ are vector elements of area of two such open surfaces, with ๐ง,โ€‰ ๐งโ€ฒ having the same sense relative to the vortex tube, the Gauss-Green divergence theoremMathworldPlanetmathPlanetmath applied to the control volume enclosed by these two surfaces and the connecting (lateral) portion of the vortex tube, shows that

โˆซ๐’œ๐Žโ‹…๐งโข๐‘‘๐”ž-โˆซ๐’œโ€ฒ๐Žโ€ฒโ‹…๐งโ€ฒโข๐‘‘๐”žโ€ฒ=โˆซ๐”ณโˆ‡โ‹…๐Žโข๐‘‘๐”ณ=0,

where ๐”ณ is the Eulerian description of the control volume in question. Note that there is no contribution to the surface integral from the portion of the vortex tube. The flux of vorticity along a vortex tube is thus independent of the choice of the open surface used to measure it, and is termed the strength of the vortex tube. In the case of a vortex tube of infinitesimal cross section (usually called filament-tube), such strength is equal to the productPlanetmathPlanetmath of cross-sectional area and the magnitude of the local vortex vector, being the same at all stations along the vortex tube. It is very important to mention that a vortex tube cannot begin or end in the interior of the fluid flow, but must either be a closed tube, like a torus, or else (provided it does not meet a boundary) must extend ad infinitum in either direction. For at an end, if were one, a continuousMathworldPlanetmathPlanetmath transition would be possible, along the mantle of the vortex tube, between contractible closed curves ๐’ž1 located there and cross-sectional contractible one ๐’ž2, which is inconsistentPlanetmathPlanetmath with the fact that ฮ“1=0 while ฮ“2=cโขoโขnโขsโขtโขaโขnโขtโ‰ 0.
An extensive and detailed bibliographical data is given in [7].

0.5 Acknowledgement

To Cameron McLeman=mathcam, for his clear explanation about the distinction between reducible and contractible curves and its โ€˜closednessโ€™. Until about the mid-past century, mechanicists used the word โ€˜reducibleโ€™ like synonymous of โ€˜contractibleโ€™. See, for example, [8]

References

  • 1 W. Thomson, On vortex motion, Trans. Roy. Soc. Edinburgh, 25, 1869.
  • 2 W. Thomson, P.G. Tait, Treatise on Natural Philosophy, Part I (1879), Part II (1883), Cambridge University Press, 1912.
  • 3 H. v. Helmholtz, รœber Integrale der hydrodynamischen Gleishungen, welche den Wirbelbewegugen entsprechen, J. reine angew. Math. 55, pp. 25-55, 1858.
  • 4 J. L. Lagrange, Mรƒยƒร‚ยฉmoire sur la thรƒยƒร‚ยฉorie du mouvement des fluides, Nouv. Mรƒยƒร‚ยฉm. Acad. Berlin, pp. 151-198, 1781 = Ouevresโ€‰ 4, pp. 695-748, 1783.
  • 5 A. L. Cauchy, Memoire sur les dilatations, les condensation, et les rotations produites par un changement de forme dans un systรƒยƒร‚ยจme de points matรƒยƒร‚ยฉriels, Ouvres CompletesPlanetmathPlanetmathPlanetmathPlanetmathPlanetmathPlanetmath, Ser. 2, Vol. 12, pp. 343-377, Paris: Gauthier-Villars, 1916.
  • 6 G. G. Stokes, On the theories of the internal friction of fluids in motion, and of the equilibrium and motion of elastic solids, Trans. Cambridge Phil. Soc.โ€‰ 8, p. 287 ff., 1845.
  • 7 C. Truesdell, The Kinematics of Vorticity, Bloomington: Indiana Univ. Press, 1954.
  • 8 G. K. Batchelor, An Introduction to Fluid Dynamics,โ€‰ 2, p. 92, Cambridge University Press, 1967.
Title circulation and vorticity
Canonical name CirculationAndVorticity
Date of creation 2016-05-24 23:17:48
Last modified on 2016-05-24 23:17:48
Owner perucho (2192)
Last modified by perucho (2192)
Numerical id 13
Author perucho (2192)
Entry type Definition
Classification msc 53A45
Related topic SourcesAndSinksOfVectorField