Hicks equation explained

In fluid dynamics, Hicks equation, sometimes also referred as Bragg–Hawthorne equation or Squire–Long equation, is a partial differential equation that describes the distribution of stream function for axisymmetric inviscid fluid, named after William Mitchinson Hicks, who derived it first in 1898.^{[1] [2] [3]} The equation was also re-derived by Stephen Bragg and William Hawthorne in 1950 and by Robert R. Long in 1953 and by Herbert Squire in 1956.^{[4] [5] [6]} The Hicks equation without swirl was first introduced by George Gabriel Stokes in 1842.^{[7] [8]} The Grad–Shafranov equation appearing in plasma physics also takes the same form as the Hicks equation.

Representing

(r,\theta,z) as coordinates in the sense of cylindrical coordinate system with corresponding flow velocity components denoted by (v _{r,v} _{\theta,v} _{z)} , the stream function \psi that defines the meridional motion can be defined asrv _{r} = - (\partial \psi)/(\partial z) , rv _{z} = (\partial \psi)/(\partial r)

that satisfies the continuity equation for axisymmetric flows automatically. The Hicks equation is then given by ^[9]

(\partial ^{2} \psi)/(\partial r ^{2}) - (1)/(r) (\partial \psi)/(\partial r) + (\partial ^{2} \psi)/(\partial z ^{2}) = r ^{2} (d H)/(d \psi) - \Gamma (d \Gamma)/(d \psi)

where

H(\psi) = (p)/(\rho) + (1)/(2) (v ^{2 ),} _{z} \Gamma(\psi) = rv _{\theta}

where

H(\psi) is the total head, c.f. Bernoulli's Principle. and 2 \pi\Gamma is the circulation, both of them being conserved along streamlines. Here, p is the pressure and \rho is the fluid density. The functions H(\psi) and \Gamma(\psi) are known functions, usually prescribed at one of the boundary; see the example below. If there are closed streamlines in the interior of the fluid domain, say, a recirculation region, then the functions H(\psi) and \Gamma(\psi) are typically unknown and therefore in those regions, Hicks equation is not useful; Prandtl–Batchelor theorem provides details about the closed streamline regions.

Derivation

Consider the axisymmetric flow in cylindrical coordinate system

(r,\theta,z) with velocity components (v _{r,v} _{\theta,v} _{z)} and vorticity components (\omega _{r,\omega} _{\theta,\omega} _{z)} . Since \partial/\partial \theta =0 in axisymmetric flows, the vorticity components are \omega _{r} =- (\partial v _{\theta})/(\partial z) , \omega _{\theta=} (\partial v _{r})/(\partial z) - (\partial v _{z})/(\partial r) , \omega _{z} = (1)/(r) (\partial (rv _{\theta)})/(\partial r) .

Continuity equation allows to define a stream function

\psi(r,z) such thatv _{r =- (1)/(r)} (\partial \psi)/(\partial z) , v _{z} = (1)/(r) (\partial \psi)/(\partial r)

(Note that the vorticity components

\omega _{r} and \omega _{z} are related to rv _{\theta} in exactly the same way that v _{r} and v _{z} are related to \psi ). Therefore the azimuthal component of vorticity becomes\omega _{\theta} = - (1)/(\left( r) (\partial ^{2\psi})/(\partial r ^{2}) - (1)/(r) (\partial \psi)/(\partial r) + (\partial ^{2\psi})/(\partial z ^{2}) \right).

The inviscid momentum equations

\partial\boldsymbol{v}/\partial t - \boldsymbol{v} x \boldsymbol{\omega} = - \nabla H , where H = (1)/(2) (v ^{2 )} _{z} + (p)/(\rho) is the Bernoulli constant, p is the fluid pressure and \rho is the fluid density, when written for the axisymmetric flow field, becomes\begin{align} v _{\theta} \omega _{z} - v _{z\omega} _{\theta} - (\partial v _{r})/(\partial t) & = (\partial H)/(\partial r) ,\\ v _{z\omega} _{r} - v _{r} \omega _{z} - (\partial v _{\theta})/(\partial t) & =0 ,\\ v _{r\omega} _{\theta} - v _{\theta} \omega _{r} - (\partial v _{z})/(\partial t) & = (\partial H)/(\partial z) \end{align}

in which the second equation may also be written as

D(rv _{\theta)/Dt=0} , where D/Dt is the material derivative. This implies that the circulation 2 \pi rv _{\theta} round a material curve in the form of a circle centered on z -axis is constant.

If the fluid motion is steady, the fluid particle moves along a streamline, in other words, it moves on the surface given by

\psi = constant. It follows then that H = H(\psi) and \Gamma = \Gamma(\psi) , where \Gamma = rv _{\theta} . Therefore the radial and the azimuthal component of vorticity are\omega _{r} = v _{r (d\Gamma)/(d\psi)} , \omega _{z} = v _{z (d\Gamma)/(d\psi)} .

The components of

\boldsymbol{v} and \boldsymbol{\omega} are locally parallel. The above expressions can be substituted into either the radial or axial momentum equations (after removing the time derivative term) to solve for \omega _{\theta} . For instance, substituting the above expression for \omega _{r} into the axial momentum equation leads to^[9] \begin{align} (\omega _{\theta})/(r) & = (v _{\theta} \omega _{r})/(r v _{r}) + (1)/(rv _{r}) (dH)/(d\psi) (\partial \psi)/(\partial z) \\ & = (\Gamma)/(r ^{2}) (d\Gamma)/(- d\psi) (dH)/(d\psi) . \end{align}

But

\omega _{\theta} can be expressed in terms of \psi as shown at the beginning of this derivation. When \omega _{\theta} is expressed in terms of \psi , we get(\partial ^{2} \psi)/(\partial r ^{2}) - (1)/(r) (\partial \psi)/(\partial r) + (\partial ^{2} \psi)/(\partial z ^{2}) = r ^{2} (d H)/(d \psi) - \Gamma (d \Gamma)/(d \psi) .

This completes the required derivation.

Example: Fluid with uniform axial velocity and rigid body rotation in far upstream

Consider the problem where the fluid in the far stream exhibit uniform axial velocity

U and rotates with angular velocity \Omega . This upstream motion corresponds to \psi = (1)/(2) Ur ^{2,} \Gamma = \Omega r ^{2,} H = (1)/(2) U ^{2} + \Omega ^{2} r ^{2.}

From these, we obtain

H(\psi) = (1)/(2) U ^{2} + (2 \Omega ^{2})/(U) \psi, \Gamma(\psi) = (2 \Omega)/(U) \psi

indicating that in this case,

H and \Gamma are simple linear functions of \psi . The Hicks equation itself becomes(\partial ^{2} \psi)/(\partial r ^{2}) - (1)/(r) (\partial \psi)/(\partial r) + (\partial ^{2} \psi)/(\partial z ^{2}) = (2 \Omega ^{2})/(U) r ^{2} - (4 \Omega ^{2})/(U ^{2}) \psi

which upon introducing

\psi(r,z) = Ur ^{2/2} + r f(r,z) becomes(\partial ^{2} f)/(\partial r ^{2}) + (1)/(r) (\partial f)/(\partial r) + (\partial ^{2} f)/(\partial z ^{2}) + \left(k ^{2- (1)/(r ^{2})} \right) f = 0

where

k =2 \Omega/U .

Yih equation

For an incompressible flow

D\rho/Dt =0 , but with variable density, Chia-Shun Yih derived the necessary equation. The velocity field is first transformed using Yih transformation(v _{r',v} _{\theta',v} _{z')} = \sqrt{ (\rho)/(\rho _{0}) }(v_r,v_\theta,v_z)

where

\rho _{0} is some reference density, with corresponding Stokes streamfunction \psi' defined such thatrv _{r'} = - (\partial \psi')/(\partial z) , rv _{z'} = (\partial \psi')/(\partial r) .

Let us include the gravitational force acting in the negative

z direction. The Yih equation is then given by^{[10] [11]} (\partial ^{2} \psi')/(\partial r ^{2}) - (1)/(r) (\partial \psi')/(\partial r) + (\partial ^{2} \psi')/(\partial z ^{2}) = r ^{2} (d H)/(d \psi') - r ^{2} (d \rho)/(d \psi') (g)/(\rho _{0}) z - \Gamma (d \Gamma)/(d \psi')

where

H(\psi') = (p)/(\rho _{0}) + (\rho)/(2 \rho _{0}) (v ^{2 )} _{z'} + (\rho)/(\rho _{0}) g z, \Gamma(\psi') = rv _{\theta'}

Notes and References

1. Hicks, W. M. (1898). Researches in vortex motion. Part III. On spiral or gyrostatic vortex aggregates. Proceedings of the Royal Society of London, 62(379–387), 332–338. https://royalsocietypublishing.org/doi/pdf/10.1098/rspl.1897.0119
2. Hicks, W. M. (1899). II. Researches in vortex motion.—Part III. On spiral or gyrostatic vortex aggregates. Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character, (192), 33–99. https://royalsocietypublishing.org/doi/pdf/10.1098/rsta.1899.0002
3. Smith, S. G. L., & Hattori, Y. (2012). Axisymmetric magnetic vortices with swirl. Communications in Nonlinear Science and Numerical Simulation, 17(5), 2101–2107.
4. Bragg, S. L. & Hawthorne, W. R. (1950). Some exact solutions of the flow through annular cascade actuator discs. Journal of the Aeronautical Sciences, 17(4), 243–249
5. Long, R. R. (1953). Steady motion around a symmetrical obstacle moving along the axis of a rotating liquid. Journal of Meteorology, 10(3), 197–203.
6. Squire, H. B. (1956). Rotating fluids. Surveys in Mechanics. A collection of Surveys of the present position of Research in some branches of Mechanics, written in Commemoration of the 70th Birthday of Geoffrey Ingram Taylor, Eds. G. K. Batchelor and R. M. Davies. 139–169
7. Stokes, G. (1842). On the steady motion of incompressible fluids Trans. Camb. Phil. Soc. VII, 349.
8. Lamb, H. (1993). Hydrodynamics. Cambridge university press.
9. Batchelor, G. K. (1967). An introduction to fluid dynamics. Section 7.5. Cambridge university press. section 7.5, p. 543-545
10. Yih, C. S. (2012). Stratified flows. Elsevier.
11. Yih, C. S. (1991). On stratified flows in a gravitational field. In Selected Papers By Chia-Shun Yih: (In 2 Volumes) (pp. 13-21).