## 1. Introduction

Vortices are found in classical, relativistic and quantum fluids [1–7]. It is well known from classical theory of vortices that the dynamics of vortices is described by Kelvin-Helmholtz’s flux conservation theorem. There are three theorems of Helmholtz on vortex tubes in the classical theory of fluids [2]. The product of the magnitude of fluid’s vorticity vector and of vortex tube’s cross-sectional area is called the strength of a vortex tube. Helmholtz’s first theorem states that the strength of vortex tube remains constant along the vortex tube. This result is purely kinematical in nature in the sense that its derivation does not require Euler’s equations which govern the motion of a fluid. The second theorem tells us that the vortex lines are material lines in the case of perfect fluid flows. The third theorem says that the vortex tube’s strength multiplied by the fluid’s chemical potential remains constant along the flow lines (or streamlines) for perfect fluid flows. Since a two-dimensional surface across which fluid’s vorticity flux passes is usually bounded by a closed circuit, Kelvin’s circulation theorem is intimately connected with that of Helmholtz’s flux conservation theorem for perfect fluid flows. The circulation of a fluid motion around any closed curve on a vertex tube wall is equal to the flux of vorticity along a vertex tube. Kelvin’s circulation theorem asserts that the circulation of a perfect fluid motion around a closed material curve is conserved in time as the fluid evolves. This assertion in turn implies that the vortex tube is a material tube and moves with the fluid (i.e., vortex lines are material lines) [1, 2].

Since relativistic fluid dynamics differs from that of Newtonian fluid dynamics in many ways, Greenberg [8] developed a theory of a family of spacelike curves in order to study the kinematical behavior of vortex lines (i.e., spacelike curves) appearing in an isentropic perfect fluid motion. Greenberg [8] formulated a relativistic analogue of Helmboltz’s first theorem on the basis of Ehlers’ divergence identity of fluid’s vorticity vector [9] in the Newtonian limit (i.e., ignoring the fluid’s acceleration term). The relativistic divergence identity of fluid’s vorticity vector involves fluid’s acceleration vector which embodies the curvature of space-time, and thereby, it differs from classical divergence identity of vorticity vector. This is the reason that the exact relativistic version of Helmholtz’s first theorem is not obtainable kinematically. Greenberg [8] further extended his analysis to obtain a relativistic version of Helmholtz’s third theorem using Euler’s equations of motion for an isentropic perfect fluid. In this derivation, the concept of vortex lines to be material lines is automatically implemented because of his transport law (i.e., it is presently known as Greenberg’s transport law).

The necessary and sufficient conditions for the existence of conservation of flux associated with a vector field have been given firstly by Bekenstein and Oron [10] and then secondly by Carter [11]. Of two conditions, the first condition is satisfied if a vorticity 2-form is expressed as a curl of momentum covector associated with a fluid and the second condition when imposed on vorticity 2-form gives Euler’s equations that govern fluid’s motion. Vorticity 2-form satisfying these two conditions led to the formulation of a relativistic version of Kelvin’s circulation conservation theorem and Helmholtz’s flux conservation theorem [10, 11].

There is a link between Greenberg’s transport law [8] and Carter’s formulation [11] of vorticity flux conservation. Carter’s vorticity flux conservation is based on vorticity 2-form of rank 2 (i.e., simple vorticity bivector field). The matrix of vorticity 2-form admits eigen vectors with zero eigen values. Such eigen vectors are referred to as flux vectors. The flux vectors span two-dimensional tangent subspaces at each point of a four-dimensional space-time. These tangent subspaces mesh together to form a family of timelike 2-surfaces because simple vorticity bivector field satisfies Frobenius condition [12] of 2-surface forming. Such timelike 2-surfaces are called flux surfaces. It is to be noted that Greenberg’s transport law is an alternative version of Frobenius condition [12] due to which congruences of fluid flow lines and vortex lines form a family of timelike two-dimensional surfaces in the case of an isentropic perfect fluid flows. Furthermore, vortex lines are material lines. Flux surfaces are also material surfaces.

An intimate connection between Kelvin’s circulation theorem and the fluid helicity (i.e., helicity of vorticity vector field) conservation was long ago demonstrated by Moffat [13]. The geometric structure of streamlines and vortex lines (forming a vortex tube) inherent in the helicity conservation was discovered by expressing the conservation of the sum of writhe and twist of vortex tubes undergoing continuous deformation in a fluid flow [14, 15]. The fluid helicity is conserved when vortex lines are frozen lines (i.e., material lines) in a fluid flow [14]. Such topological consequences have been extensively investigated by several authors [16, 17] in classical fluid dynamics.

Although relativistic version of fluid helicity conservation has been formulated by Carter [11] and Bekenstein [18], yet its topological consequences remained unknown. In a relativistic frame work, frozen-in property of spacelike vortex lines is describable by Greenberg’s transport law [8], which in turn leads to the fact that an isolated vortex is a timelike two-dimensional connected manifold called vortex world sheet. The deformation of timelike fluid flow lines and spacelike vortex lines caused by gravitation forbids the application of Biot-Savart-like law and Frenet-like transport frame along a spacelike vortex line. Furthermore, the kinematic deformation of a vortex tube is not independent of the kinematic deformation of fluid flow lines forming a stream tube. The deformation of a vortex tube is determined in terms of the expansion, shear and acceleration associated with fluid flow lines and the magnetic part of Weyl curvature tensor representing the free gravitational field [19]. In particular, the variation in cross-sectional area of a vortex tube along the tube is controlled by free gravitational field even if the acceleration vector of the fluid motion and the shear are ignored [19].

An inherent connection between helicity conservation and streamline invariant has already been pointed out by Bekenstein [18] in order to extract relevant information of astrophysical significance from nonlinear equations of relativistic fluid (or magnetofluid) dynamics. But such relationship between fluid helicity conservation and inherently connected streamline invariant is still to be investigated. This is the idea which motivates to work further in order to understand the role of vorticity in the rotational evolution of relativistic fluid, since most, if not all, compact stars composed of fluids rotate. It may be conceived that the gravitational effect on vorticity can be significant in the understanding of the internal structures of compact stars. It is known that differential rotation of fluid elements of which stars are composed arises from the gravitational collapse of massive stellar cores [20, 21]. Differential rotation produces vorticity which twists streamlines in fluid’s motion. It is expected that the winding up of vortex lines caused by differential rotation will change the angular velocity profile in a similar way as has been found in magnetized stars due to the presence of magnetic fields [22].

In the present work, we confine our attention to the derivation of streamline invariants which adhere to the fluid helicity conservation and a formulation of vorticity winding due to meridional circulation via differential rotation along poloidal components of vorticity. Furthermore, we derive new conservation law for vorticity flux.

The present paper is organized as follows. In Section 2, relativistic versions of Helmholtz’s three theorems on vorticity flux conservation are derived using an alternative approach based on Greenberg’s theory of spacelike congruence generated by vortex lines and (1 + 1) + 2 decomposition of gradient of fluid’s 4-velocity. Section 3 is devoted to the derivation of streamline invariants and the formulation of vorticity winding. In Section 4, we obtain new conservation law for the vorticity flux.

*Convention*: Space-time metric is of signature +2. Semi-colon (;) and comma (,) are, respectively, used to denote covariant derivative and partial derivative. Speed of light c is assumed to be unity.

## 2. Helmholtz’s theorems

In this section, we obtain an alternative derivation of exact relativistic version of Helmholtz’s theorems. We demonstrate that a relativistic version of Helmholtz’s first theorem describes vorticity flux conservation inside a vortex tube when the variation of flux is measured by a comoving observer with the fluid along a vortex tube. Such conservation of flux is independent of proper time. We will also show that a relativistic version of Helmholtz’s third theorem describes vorticity flux conservation in a streamline tube whose cross-sectional area lies in a spacelike two-dimensional subspace orthogonal to both stream and vortex lines. Such flux conservation is in proper time as the perfect fluid evolves and vortex lines are material lines. In order to prove these two results, we begin with Euler’s equation that governs the motion of a perfect fluid. Euler’s equation is as follows [23]:

where *n* is conserved

The first law of thermodynamics is [23]:

where *T* and *s* denote, respectively, the local temperature and the entropy per baryon. The relativistic enthalpy per baryon called chemical potential is expressed as:

If we assume that the entropy per baryon

Substituting Eqs. (2.4) and (2.5) in Eq. (2.1), we can reduce Euler’s equation in the following form:

It is known from relativistic version of Kelvin’s circulation theorem for a perfect fluid that the vorticity flux equals the closed contour line integral *C*. The conservation of vorticity flux means that the value of this integral is the same for all times when each point in the contour *C* is dragged along the fluid flow (mathematically, it is Lie transported along

Contraction of Eq. (2.7) with

Dualizing Eq. (2.7) and contracting the resulting equation with

where an overhead star (*) is used for Hodge dualization.

Inverting Eq. (2.9), we get

which satisfies

It is evident from Eqs. (2.8) and (2.11) that

Since

Because

Substituting Eq. (2.12) in Eq. (2.13) and contracting the resulting equation with

and

where

A spacelike congruence of vortex lines is generated by unit spacelike vector field *V* denotes the magnitude of the vorticity flux vector

Since the vorticity flux vector *u*^{a} and *m*^{a}, respectively. Greenberg’s expansion parameter

and

where *σ* and

which is an exact relativistic version of Helmholtz’s first theorem. It is evident from Eq. (2.18) that the variation in vorticity flux of a vortex tube is measured along the tube and the flux remains constant inside a tube. The proper time *τ* plays no role in such variation of flux. Since an observer employed to measure flux is comoving along vortex line, a vortex tube is a material tube. In relativistic framework, a vortex tube to be a material tube is automatically associated with Helmholtz’s first theorem because its derivation originates from Euler’s equation. An exact relativistic analogue of Helmholtz’s first theorem is not obtainable kinematically.

In order to obtain Helmholtz’s third theorem, we use (1 + 1) + 2 decomposition of the gradient of 4-velocity

and

where *τ* and

which is a relativistic version of Helmholtz’s third theorem. It is observed from Eq. (2.21) that the variation in vorticity flux is measured along a streamline and the flux inside a stream tube remains constant in proper time *τ* as the fluid evolves according to Euler’s equation of motion. This result can be understood in a sense that the volume occupied by a stream tube is the product of cross-sectional area of a stream tube and its length measured along a vortex line (being a material line). The vorticity flux passing through such cross-sectional area remains constant in proper time *τ*.

On account of Eq. (2.12), Eq. (2.13) can be put in two different forms which are given as below:

where

In the derivation of Eq. (2.23), baryon conservation law Eq. (2.2) is used.

Projecting Eq. (2.22) orthogonal to both

which is Greenberg’s transport law [8]. Eq. (2.24) tells us that the timelike congruence generated by

Eq. (2.23) is in a form as has been derived by Bekenstein and Oron [10] in the case of the magnetic field vector in order to show the frozen-in property of the magnetic field in a perfectly conducting magnetofluid dynamics. Similar conclusion holds for the vorticity flux vector

## 3. Vorticity winding in axisymmetric stationary fluid configuration

This section is devoted to the study of vorticity winding caused by meridional circulation via differential rotation of a perfect fluid which is assumed to be axisymmetric and stationary. In this case, a space time admits two linearly independent commuting Killing vectors: *t* and *∅* are usually called toroidal coordinates. We choose poloidal coordinates *t* and *∅*.

From Eq. (2.8), we get

Since

where the square bracket around indices indicates skew-symmetrization. From Eq. (3.2), we get

It follows from Eqs. (3.1a) and (3.3a) that

Similarly, substitution of Eq. (3.1b) in Eq. (3.3b) gives

From Eqs. (3.4a) and (3.4b), we get

which on integration along the streamline generated by

where *A* is constant along the streamline. Eq. (3.4b) can be rewritten as:

From Eq. (2.2), we get

Substitution of Eq. (3.8) in Eq. (3.7) gives

which on integration along the streamline gives

where *B* is constant along the streamline.

Substituting Eq. (3.10) in Eq. (3.1b), we get

Substituting Eqs. (3.6) and (3.10) in Eq. (3.1c), we get

where

It follows from Eqs. (3.6), (3.10) and (3.11) that

Eq. (2.9) gives

Substitution of Eqs. (3.10)–(3.13) in the coordinate expansion of Eq. (3.14) gives

which is the required expression for the vorticity flux vector. We now use Eq. (3.15) to find conserved quantity along a streamline from the conservation of the fluid helicity. The fluid helicity vector is given by Carter [11, 24] and Bekenstein [18] as follows

whose divergence vanishes,

This is fluid helicity conservation. Eq. (3.17) can be cast in form

or equivalently

Substituting Eq. (3.15) in Eq. (3.19), we get

which on integration along a streamline gives

where the constant *B* is absorbed in a constant *C* along a streamline. The energy and angular momentum scalars in axisymmetric stationary fluid configuration are of the forms [11]

From Eqs. (3.21) and (3.22), it can be seen that

which shows that the linear combination of energy and angular momentum is constant along a streamline but varies from streamline to streamline. Thus, *C* is a streamline invariant.

We now consider in view of Eq. (3.15) as follows

Setting *μ*

and using Eq. (3.21) in the resulting equation, we get

In order to provide a meaningful interpretation of Eq. (3.25), we proceed as follows. In the absence of meridional circulation (i.e., *A* is the angular velocity of vorticity flux surfaces commoving with the fluid in the absence of meridional circulation. But in the presence of meridional circulation (i.e.,

where

It is evident from Eq. (3.26) that the toroidal components

The left-hand side of Eq. (3.28) represents the differential rotation along the poloidal components of the vorticity vector, whereas the right-hand side gives the proper time rate of change of a combination of toroidal components of the vorticity vector. This implies that the differential rotation along the poloidal vorticity cannot vanish until the vanishing of toroidal vorticity. Since the presence of meridional circulation ensures the existence of toroidal vorticity, meridional circulation causes the kind of effect that mimics like vorticity winding by stretching frozen-in vortex lines via differential rotation.

Because

which is the law of gravitational isorotation as is pointed out by Glass [26]. Thus, it seems that the meridional circulation causes vorticity winding due to the breakdown of the gravitational isorotation.

Further investigation is needed on the lines of recent work by Birkl et al. [27] in order to understand the role of meridional circulation in relation to vorticity winding via differential rotation.

## 4. Vorticity energy conservation in axisymmetric stationary case

This section is devoted to the derivation of vorticity energy conservation assuming that the perfect fluid configuration is axisymmetric and stationary. Substituting Eq. (2.12) in Eq. (2.13) and contracting the resulting equation with *V*_{a}, we get

Eliminating

which can be converted to the form

The last term on the right-hand side of Eq. (4.3) with the help of Eq. (3.15) can be simplified in the following form

Substitution of Eq. (4.4) in Eq. (4.3) gives

Since *B* and *A* are streamline invariant, we obtain from Eq. (4.5) that

which on integration along a streamline gives that

Since *V*^{a} = *μω*^{a}, *V*^{2} = *μ*^{2}*ω*^{2}, *V*_{t }= *μω*_{t}, and *V*_{∅ }= *μω*_{∅}. Thus, Eq. (4.7) takes the form

The first term of Eq. (4.8) represents the energy of the vorticity flux vector per baryon. The second term indicates the presence of covariant toroidal components of the vorticity vector *ω*^{a}. The quantity *D* is streamline invariant of the motion of frozen-in vortex lines and can be thought of as a vorticity energy conservation arising from fluid helicity conservation.