Huggins Vorticity Flux Tensor Study

• Huggins vorticity flux tensor is an anti-symmetric quantity that describes vorticity transport through advection, tilting, stretching, and viscous diffusion.
• Its dual axial-vector representation directly connects the tensor to the Josephson–Anderson relation, linking vorticity dynamics with drag production in various flows.
• Applications in laminar, turbulent, and periodic-channel flows illustrate how advective, viscous, and body-force fluxes contribute to drag and pressure recovery.

Searching arXiv for papers on the Huggins vorticity-flux tensor and Josephson–Anderson relation. arxiv_search(query="Huggins vorticity flux tensor Josephson Anderson drag", max_results=10, sort_by="relevance") Attempting to retrieve relevant arXiv records. The Huggins vorticity flux tensor, usually denoted $\Sigma$, is the anti-symmetric tensor that appears in the Helmholtz equation for vorticity,

$\partial_t \omega + \nabla\cdot\Sigma = 0,$

and encodes the transport of vorticity by advection, vortex tilting and stretching, viscosity, and, in generalized formulations, non-conservative body-force. In recent analyses of viscous bluff-body flows, it serves as the local quantity that links vorticity dynamics to drag through the detailed Josephson–Anderson relation. That role is developed for flow over a sphere at $Re=\{200,3700\}$ and over a prolate spheroid at $Re=3000$ and $20^\circ$ incidence, and is also extended to streamwise-periodic channel flows (Du et al., 12 Jul 2025, Kumar et al., 2024).

1. Definition and equivalent representations

For the bluff-body flows over a sphere and a prolate spheroid, $\Sigma$ is introduced in index form as

$$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$

or, equivalently,

$$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$

The tensor is anti-symmetric, $\Sigma_{ij}=-\Sigma_{ji}$, so it contains three independent components rather than nine (Du et al., 12 Jul 2025).

The same study also introduces the dual axial-vector form

$$\eta \equiv u\times\omega - \nu\nabla\times\omega,\qquad \Sigma_{ij}=\epsilon_{ijk}\eta_k,\qquad \eta_i=\tfrac12\epsilon_{ijk}\Sigma_{jk}.$$

This representation is central in the drag formula because the Josephson–Anderson relation is written directly in terms of $\partial_t \omega + \nabla\cdot\Sigma = 0,$0 as well as $\partial_t \omega + \nabla\cdot\Sigma = 0,$1 (Du et al., 12 Jul 2025).

In the streamwise-periodic channel formulation, the tensor is written more generally as

$\partial_t \omega + \nabla\cdot\Sigma = 0,$2

with $\partial_t \omega + \nabla\cdot\Sigma = 0,$3 any non-conservative body-force. In that formulation, $\partial_t \omega + \nabla\cdot\Sigma = 0,$4 is sometimes denoted $\partial_t \omega + \nabla\cdot\Sigma = 0,$5 in other works (Kumar et al., 2024).

2. Constituent fluxes and planar interpretation

A standard decomposition is

$\partial_t \omega + \nabla\cdot\Sigma = 0,$6

with

$\partial_t \omega + \nabla\cdot\Sigma = 0,$7

The first term is the advective flux, which physically transports each component $\partial_t \omega + \nabla\cdot\Sigma = 0,$8 of vorticity along the velocity $\partial_t \omega + \nabla\cdot\Sigma = 0,$9 and re-orients it via tilting and stretching. The second term is the viscous flux, the anti-symmetric gradient of vorticity, and embodies down-gradient diffusion of vorticity by viscosity (Du et al., 12 Jul 2025).

In the channel generalization, a third contribution appears when non-conservative body-force is present. The paper describes the convective term as a “Reynolds-stress” contribution and the body-force term as a Magnus-type contribution (Kumar et al., 2024).

Contribution	Expression	Interpretation
Advective flux	$Re=\{200,3700\}$0	Transport plus tilting and stretching
Viscous flux	$Re=\{200,3700\}$1	Down-gradient diffusion
Body-force flux	$Re=\{200,3700\}$2	Magnus-type contribution

For any fixed direction $Re=\{200,3700\}$3, the planar conservation law is

$Re=\{200,3700\}$4

Accordingly, $Re=\{200,3700\}$5 is the instantaneous flux of the vorticity component $Re=\{200,3700\}$6 within the plane normal to $Re=\{200,3700\}$7. This formulation is used explicitly for azimuthal vorticity in the axisymmetric sphere and for streamwise vorticity in the prolate spheroid (Du et al., 12 Jul 2025).

3. Detailed Josephson–Anderson relation

For an isolated bluff body held fixed in a uniform stream $Re=\{200,3700\}$8, the detailed Josephson–Anderson relation is written as

$Re=\{200,3700\}$9

Here $Re=3000$0 is the background potential-flow velocity, $Re=3000$1 labels streamtubes of the potential flow, and $Re=3000$2 is a line element along $Re=3000$3 (Du et al., 12 Jul 2025).

Because $Re=3000$4 selects only those $Re=3000$5 whose flux index $Re=3000$6 and vorticity index $Re=3000$7 are both orthogonal to the streamline direction $Re=3000$8, the relation equates the instantaneous drag-power $Re=3000$9 to the net rate at which streamwise mass currents cut through vorticity. In that interpretation, outward crossing of negative vorticity gives positive drag-power, whereas inward crossing gives anti-drag (Du et al., 12 Jul 2025).

The periodic-channel extension preserves the same structure. For fixed mass flux $20^\circ$0, the channel paper obtains

$20^\circ$1

and also rewrites $20^\circ$2 as a surface integral of $20^\circ$3 and, in a canonical channel with $20^\circ$4, as

$20^\circ$5

In that setting, the total drag is the integrated wall-normal flux of spanwise vorticity (Kumar et al., 2024).

4. Axisymmetric sphere at $20^\circ$6

For the laminar sphere, the flow is axisymmetric about $20^\circ$7, and only the azimuthal component $20^\circ$8 is nonzero. Using cylindrical coordinates $20^\circ$9, the study defines local bases $\Sigma$0 aligned with the streamfunction, its normal, and the azimuthal direction. The tensor is then decomposed as

$\Sigma$1

Projecting onto $\Sigma$2 gives

$\Sigma$3

The corresponding scalar Josephson–Anderson integrands are

$\Sigma$4

(Du et al., 12 Jul 2025).

The analysis identifies three distinct roles for the tensor components. First, the viscous flux near the wall drives $\Sigma$5 into the fluid and is responsible for the local drag. Second, the advective flux in the detached shear layer moves $\Sigma$6 outward across streamlines, producing drag, and later inward near the wake centerline, producing anti-drag. Third, viscous annihilation of azimuthal vorticity at the wake centerline contributes to pressure recovery. The paper states that the Josephson–Anderson relation is first demonstrated for the sphere at $\Sigma$7 (Du et al., 12 Jul 2025).

5. Impulsively started and turbulent sphere at $\Sigma$8

For the sphere at $\Sigma$9, the tensor is evaluated both instantaneously during impulsive start and in the stationary turbulent regime. In the averaged turbulent wake, the decomposition is extended to

$$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$0

where $$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$1 is mean-flow advection of mean vorticity, $$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$2 is the Reynolds-stress-like turbulent flux of vorticity, and $$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$3 is mean viscous diffusion. The corresponding projected Josephson–Anderson contributions are $$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$4, $$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$5, and $$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$6, and the averaged drag relation is

$$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$7

(Du et al., 12 Jul 2025).

The impulsive-start analysis assigns all of the early singular drag to wall diffusion via $$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$8, rather than to a conventional pressure-plus-shear partition. As the boundary layer detaches and the wake turns turbulent, $$\Sigma_{ij} = u_i\omega_j - \omega_i u_j - \nu(\partial_i\omega_j-\partial_j\omega_i),$$9 grows and eventually dominates drag production, while $$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$0 decays (Du et al., 12 Jul 2025).

In the fully turbulent wake, mean advective flux in the shear layer dominates near the body, whereas the turbulent flux $$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$1 takes over near the centerline. The turbulent flux enhances the transport of mean azimuthal vorticity toward the wake centerline and is identified as the driver of enthalpy recovery downstream. The same mechanism is described as rapid pressure recovery, or anti-drag, farther downstream. Mean viscous flux is negligible except immediately at the wall and in the early wake (Du et al., 12 Jul 2025).

6. Prolate spheroid at $$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$2 and $$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$3 incidence

For the prolate spheroid, the computed quantity is $$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$4, the instantaneous planar flux of streamwise vorticity $$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$5 in vertical cuts. In this setting, $$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$6 shows advection of $$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$7 along primary and secondary separated layers into the free stream, while $$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$8 shows diffusion of wall-generated $$\Sigma = u\,\omega^T - \omega\,u^T - \nu\bigl(\nabla\omega-(\nabla\omega)^T\bigr).$$9. The integrated Josephson–Anderson contributions $\Sigma_{ij}=-\Sigma_{ji}$0 and $\Sigma_{ij}=-\Sigma_{ji}$1 are plotted on vertical slices at $\Sigma_{ij}=-\Sigma_{ji}$2 (Du et al., 12 Jul 2025).

The principal drag mechanism is transport of vorticity along the separated boundary layers. Advection of $\Sigma_{ij}=-\Sigma_{ji}$3 in the primary separation, denoted $\Sigma_{ij}=-\Sigma_{ji}$4 in the paper’s figure discussion, is the main driver of drag. The secondary layer, denoted $\Sigma_{ij}=-\Sigma_{ji}$5, and the large-scale vortices redistribute vorticity and can locally reduce drag. The abstract states the same result more compactly: primary and secondary separation contribute oppositely to the drag force, while the large-scale vortices only re-distribute vorticity (Du et al., 12 Jul 2025).

The viscous flux remains confined to the attached and immediately separated boundary layers, where it establishes the local wall-pressure gradients that drive the initial detachment. The study also proposes a mechanism for secondary separation based on the theory of vortex-induced separation (Du et al., 12 Jul 2025).

7. Generalizations, scope, and boundary-condition correction

The periodic-channel study places the Huggins tensor in a broader setting. It states that the detailed Josephson–Anderson relation was first derived by Huggins for quantum superfluids, but also holds for internal flows of classical fluids and for external flows around solid bodies, corresponding there to relations of Burgers, Lighthill, Kambe, Howe and others. In all of these results, a background potential Euler flow with the same inflow and outflow as the physical flow is used, in the same spirit as Kelvin’s minimum energy theorem (Kumar et al., 2024).

A technical issue arises in streamwise-periodic channels. The paper argues that Huggins’s original use of pure Neumann boundary conditions for the background potential creates an unphysical vortex sheet in a periodic channel, because only normal-velocity continuity is enforced and a tangential-velocity jump remains across the periodic join. The proposed correction is to impose Dirichlet conditions on $\Sigma_{ij}=-\Sigma_{ji}$6 at inflow and outflow, together with Neumann conditions on the sidewalls. The resulting $\Sigma_{ij}=-\Sigma_{ji}$7 is $\Sigma_{ij}=-\Sigma_{ji}$8-periodic and $\Sigma_{ij}=-\Sigma_{ji}$9, with no spurious vortex sheet; the orthogonality $$\eta \equiv u\times\omega - \nu\nabla\times\omega,\qquad \Sigma_{ij}=\epsilon_{ijk}\eta_k,\qquad \eta_i=\tfrac12\epsilon_{ijk}\Sigma_{jk}.$$0 and the Josephson–Anderson derivation are preserved (Kumar et al., 2024).

This extension reinforces the tensor’s role as a unified field quantity. In the bluff-body study, $$\eta \equiv u\times\omega - \nu\nabla\times\omega,\qquad \Sigma_{ij}=\epsilon_{ijk}\eta_k,\qquad \eta_i=\tfrac12\epsilon_{ijk}\Sigma_{jk}.$$1 and its projections $$\eta \equiv u\times\omega - \nu\nabla\times\omega,\qquad \Sigma_{ij}=\epsilon_{ijk}\eta_k,\qquad \eta_i=\tfrac12\epsilon_{ijk}\Sigma_{jk}.$$2 identify where wall diffusion, shear-layer advection, turbulent mixing, and wake-centerline annihilation contribute positively or negatively to drag-power. In the periodic-channel study, the same tensor equates work by pressure drop to integrated spanwise-vorticity flux and reveals how vortex separation from a smooth bump creates drag at each time instant (Du et al., 12 Jul 2025, Kumar et al., 2024).