# Vortex Structure and Kinematics of Encased Axial Turbomachines

Peter Pelz, Paul Taubert, Ferdinand-J Cloos

► To cite this version:

Peter Pelz, Paul Taubert, Ferdinand-J Cloos. Vortex Structure and Kinematics of Encased Axial Turbomachines. 17th International Symposium on Transport Phenomena and Dynamics of Rotating Machinery (ISROMAC2017), Dec 2017, Maui, United States. hal-02980997

HAL Id: hal-02980997

https://hal.science/hal-02980997

Submitted on 27 Oct 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Vortex Structure and Kinematics of Encased Axial Turbomachines

Peter F. Pelz<fnref n="1" />*, Paul Taubert<fnref n="1" />, Ferdinand-J. Cloos<fnref n="1" />

<figure><img src="image_2.png" /></figure>

ISROMAC 2017

International Symposium on Transport Phenomena and Dynamics of Rotating Machinery

Maui, Hawaii

December 16-21, 2017

## Abstract

This investigation models the kinematics of the vortex system of an encased axial turbomachine at part load and overload applying analytical methods. So far, the influence of the casing and the tip clearance on the kinematics was solved separately. The vortex system is composed of a hub, bound and tip vortices. For the nominal operating point $\varphi \approx \varphi_{opt}$ and negligible induction, the tip vortices transform into a screw. For part load operation $\varphi \to 0$, the tip vortices wind up to a vortex ring, i.e. the pitch of the screw vanishes. The vortex ring itself is generated by bound vortices rotating at the angular frequency $\Omega$. The hub vortex induces a velocity on the vortex ring causing a rotation at the sub-synchronous frequency $\Omega_{ind} = 0.5 \Omega$. Besides, the vortex ring itself induces an axial velocity. Superimposed with the axial main flow, this results in a stagnation point at the tube wall. This stagnation point may wrongly be interpreted as dynamic induced wall stall. For overload operation $\varphi \to \infty$, the vortex system of the turbomachine forms a horseshoe, i.e. the pitch of the screw becomes infinite. Both, hub and tip vortices, are semi-infinite, straight vortex filaments. The tip vortices rotate against the rotating direction of the turbomachine, due to the induction of the hub vortex, yielding the induced frequency $\Omega_{ind} = -0.5 \Omega/s$, with the tip clearance $s$.

## Keywords

vortex dynamics – potential flow – part load – overload – sub-synchronous frequency – kinematic induced noise

## INTRODUCTION

By now, the common understanding is, that rotating stall and the resulting noise and vibration within a turbomachine is a dynamic effect. This means, that frictional forces lead to boundary layer separation and eventually stall in rotating machines [1, 2]. This understanding is recently confirmed by Cloos et al. [3] both, experimentally and analytically, for the most generic machine, a flow through a coaxial rotating circular tube. According to Cloos et al. [3], "wall stall" – a term coined by Greitzer [1] in contrast to "blade stall" – is caused at part load by the interaction of axial boundary layer and swirl boundary layer flow, i.e. the influence of centrifugal force on axial momentum. For "wall stall", the axial velocity component $u_z$ vanishes at the line $z = -z_0$, $r = a$ (axial coordinate in mean flow direction $z$, distance $z_0$ from the reference point on the line of symmetry $r = 0$, radial coordinate $r$, tube radius $a$; cf. figure 1).

This paper analyzes the flow situation in encased axial turbomachines for small viscous friction. This work shows, that "wall stall", i.e. $u_z(-z_0, a) = 0$, can also be a result of kinematics only, due to induced velocities of the vortex system superimposed with the axial main flow. The structure of the vortex system, especially the tip vortices, depends on the operating point of the turbomachine. Furthermore, the tip vortices rotate with a sub-synchronous frequency $\Omega_{ind}$ (Karstadt et al. [4], Zhu [5]). The aim of the present investigation is to analyze the influence of the vortex system in encased axial turbomachines and its circulation strength on the observed phenomena, yielding the research questions:

1. Is it possible to explain by means of analytical methods the sub-synchronous frequencies observed for turbomachines?

2. Can "wall stall" be a result of kinematics only?

To answer these questions, this work first employs vortex theory for an encased axial turbomachine, followed by the application of fundamental solutions. For machines without casing like wind turbines and screw propellers, vortex theory is well described by Betz [6], Goldstein [7], Glauert [8] and van Kuik [9]. The method is not yet developed in such a degree for encased turbomachines requiring the flow potential of a vortex ring inside a tube. Pelz et al. [10] recently derived this solution.

In addition, to enlarge the investigation on the whole operating range of a turbomachine, we investigate the structure and kinematics of the vortex system at heavy overload, applying theory of functions (complex analysis). For turbomachines, the flow number $\varphi := U/\Omega a$ defines the operating point (e.g. part load or overload). The flow number is the ratio of axial free-stream velocity $U$ to the circumferential velocity $\Omega a$, where $\Omega = 2\pi n$ is the rotational speed (the scaling to $\Omega a$ and not to $\Omega b$, with the blade tip radius $b$, is common in the context of turbomachines and therefore used here as well [11]).

For the nominal operating point $\varphi \approx \varphi_{opt}$ and negligible induction, the vortex system of an encased axial turbomachine consists of a hub, Z bound and Z tip vortices, with Z the number of blades. The tip vortices transform into helices with a pitch of $2\pi\varphi a$.

<footnote n="1">¹Chair of Fluid Systems, Technische Universität Darmstadt, Darmstadt, Germany</footnote>

<footnote>*Corresponding author: peter.pelz@fst.tu-darmstadt.de</footnote>
For part load operation $\varphi \to 0$, see figure 1, the Z helices "roll up" and form a vortex ring, i.e. the pitch of the helices vanishes. The vortex ring is continuously generated by the bound vortex system. Hence, the coaxial vortex ring strength is transient. A nice picture for this vortex ring is that of a thread spool rolling up and gaining strength over time. This picture will explain some transient phenomena using kinematic arguments only.

The case of heavy overload occurs for infinitely high flow numbers $\varphi \to \infty$; see figure 2. The hub, the bound and the tip vortices form a horseshoe, i.e. the pitch of the helices becomes infinite. Both, hub and tip vortices are semi-infinite, straight vortex filaments. In real turbomachines, the flow number cannot be adjusted to infinity, but is limited to a maximum value $\dot{\varphi}$ due to flow rate limitations and geometric restrictions. Nevertheless, the analysis of this limiting case is important for the basic understanding of the vortex system in axial turbomachines.

To develop physical understanding of the whole picture in detail, the paper is organized as follows. Section 1 gives a short literature overview. Section 2 uses vortex theory to determine the strength of the vortices. Subsequently, section 3 derives the velocity potential of a coaxial vortex ring within a circular tube at part load and the induced rotating frequency. The flow potential and the induction at overload is introduced in section 4. The paper closes with a short outlook to potential applications in section 5 and a discussion in section 6.

# 1. LITERATURE REVIEW

Investigations of vortex systems in fluid dynamics trace back to the work of Helmholtz [12], who formulated the Helmholtz's theorems as a basis for the research concerning rotational fluid motion.

Didden [13] performed measurements of the rolling-up process of vortex rings and compared the results with similarity laws for the rolling-up of vortex sheets.

Besides the investigation of vortex kinematics, a broad research field on vortex structures in turbomachines is the experimental and numerical analysis of acoustic and noise emission of tip vortices [14, 15, 16, 17]. The noise of a fan is noticeable by a CPU, car or a rail vehicle cooler. All three examples are met in the everyday life. One of the main reasons for the noise is the gap $s := (a-b)/a$ between the housing and the impeller tip. With increasing gap, the noise emission and the energy dissipation increase [18, 19]. Karstadt et al. [20] and Zhu [5] investigated noise and dissipation due to tip vortices.

The kinematics for an encased axial turbomachine operating at part load or overload, are not fully understood yet. Especially the basic kinematics of these phenomena are not sufficiently analyzed.

# 2. VORTEX THEORY

An encased axial turbomachine with Z impellers is considered. The sketched vortex system (see figures 1 and 2) results

<figure><img src="image_3.png" /><figcaption>**Figure 1.** A coaxial vortex ring of transient strength $\Gamma_t$ and radius $b$ in a circular tube of radius $a$, according to the case $\varphi \to 0$. The sketch is for $Z = 1$, i.e. one bound vortex only to improve clarity.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** A horseshoe vortex of strength $\Gamma$ and radius $b$ in a circular tube of radius $a$, according to the case $\varphi \to \infty$. The sketch is for $Z = 1$, i.e. one bound vortex only to improve clarity.</figcaption></figure>
from the Z bound vortices. The generation of a bound vortex was explained by Prandtl [21] using arguments of boundary layer theory and Kelvin's circulation theorem. The presence of viscosity is essential for the creation of the bound vortex, but the generation phase is not in the scope of this paper. For vortex generation, we would like to refer the reader to the work of Prandtl [21].

By vortex theory, each blade 1...Z of length b is represented by its bound vortex of strength $\Gamma$. For simplicity, this investigation assumes $\Gamma$ to be constant in radial direction along the blade from $r = 0$ to $r = b$. As a vortex filament cannot end in a fluid due to Helmholtz's vortex theorem, a free, trailing vortex springs at each blade end $r = 0$ and $r = b$; see figures 1 and 2. These vortices are of the same strength as the bound vortex. At the inner end $r = 0$, a straight semi-infinite vortex line $0 \le z < \infty$ of strength $Z\Gamma$ – the so called hub vortex – attaches to the blade. The tip vortices at the outer end $r = b$ are helices. The axial distance of the each helix winding, i.e. the helix pitch, is given by $U/n = 2\pi a\varphi$. Depending on the load, these helices either "wind up" ($\varphi \to 0$), forming a vortex ring, or stretch to infinity ($\varphi \to \infty$), yielding a straight, semi-infinite vortex line.

Regardless of the flow number $\varphi$, the semi-infinite straight vortex line at $r = 0$ induces the circumferential velocity $Z\Gamma/(4\pi b)$ at $z = 0$, $r = b$ due to the Biot-Savart law. Hence, the induced rotational speed is $\Omega_{ind} = Z\Gamma/(4\pi b^2)$.

In a next step, this analysis calculates the vortex strength $Z\Gamma$, employing the angular momentum equation and the energy equation. On the one hand, the axial component of the angular momentum equation is $Z\Gamma/2\pi = dM/d\dot{m}$. Here, $M$ is the axial torque component and $\dot{m}$ the mass flux. Multiplying the momentum equation by $\Omega$ yields $Z\Gamma n = dP/d\dot{m}$. $P = M\Omega$ is the power applied to the fluid by means of the rotating bound vortices. On the other hand, the energy equation for an adiabatic flow reads $dP/d\dot{m} = \Delta h_t$, with $\Delta h_t$ being the difference in total enthalpy experienced by a fluid particle passing the cross-section $z = 0$. Both arguments result in the relation $Z\Gamma n = \Delta h_t$.

From turbomachine theory, the expression $\Delta h_t = (\Omega b)^2(1 - \varphi/\dot{\varphi})$ can be derived from the equation mentioned above. The dimensionless design parameter $\dot{\varphi}$ equals the tangent of the blade's trailing edge angle $\beta_2$, i.e. $\dot{\varphi} = \tan \beta_2$. Hence, the relation between $Z\Gamma$ and $\Omega$ yields

$$ \frac{Z\Gamma n}{(\Omega b)^2} = 1 - \frac{\varphi}{\dot{\varphi}}. \qquad (1) $$

As equation 1 shows, the total change in circulation $Z\Gamma$ along the plane of the machine is linked to the flow number $\varphi$ by Euler's turbine equation.

### 3. THE VORTEX SYSTEM AT HEAVY PART LOAD

For the limiting case of interest $\varphi \to 0$, the relation between $Z\Gamma$ and $\Omega$, equation 1, yields $\Delta h_t = Z\Gamma n = (\Omega b)^2$. This

results in an induced sub-synchronous frequency

$$ \frac{\Omega_{ind}}{\Omega} = \frac{1}{2}, \quad \text{for } \varphi \to 0. \qquad (2) $$

This induced frequency is in surprisingly good agreement with measured sub-synchronous frequencies 0.5Ω...0.7Ω of rotating stall of compressors, fans and pumps at part load operation [11] and may result in a rethinking of rotating stall from a kinematic perspective.

This investigation is now set to analyze the kinematics of coaxial vortex rings of radius $b$ and maximal strength $\Gamma_t = Z\Gamma nt < (\Omega b)^2 t$ as sketched in figure 1. By doing so, Laplace's equation

$$ \frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial \Phi}{\partial r} \right) + \frac{\partial^2 \Phi}{\partial z^2} = 0 \qquad (3) $$

for the velocity potential $\Phi$, with $\vec{u} = \nabla\Phi$, is solved for an incompressible, axisymmetric flow. Pelz et al. [10] recently derived the velocity potential for a coaxial vortex ring filament in a circular tube, yielding

$$ \begin{aligned} \phi(r, z) &:= \frac{\Phi}{Ua} = \\ & = \frac{z}{a} - 2\tau\beta^2 \sum_{n=1}^{\infty} \frac{J_1(k_n\beta)}{k_n J_0^2(k_n)} J_0\left(k_n \frac{r}{a}\right) \exp\left(-k_n \frac{|z|}{a}\right), \end{aligned} \qquad (4) $$

with $J_0, J_1$ the Bessel function of orders 0 and 1, respectively, and $k_n$ the zeros $n = 1...\infty$ of the function $J_0'(k_n) = -J_1(k_n) = 0$. The dimensionless velocity potential $\phi$ depends on the dimensionless ring radius $\beta := b/a$ and the dimensionless vortex strength $\tau := \Gamma_t / 2bU$. Since $\Gamma_t$ increases linearly in time, $\tau$ can also be interpreted as a parametric time of the process.

Stokes stream function for this flow is (using the integrability conditions $\partial\psi/\partial z = -r\partial\phi/\partial r$ and $\partial\psi/\partial r = r\partial\phi/\partial z$ [10])

$$ \psi(r, z) := \frac{\Psi}{Ua} = -\frac{1}{2} \left(1 - \frac{r^2}{a^2}\right) + \\ + 2\tau\beta^2 \sum_{n=1}^{\infty} \frac{J_1(k_n\beta)}{k_n J_0^2(k_n)} J_1\left(k_n \frac{r}{a}\right) \frac{r}{a} \exp\left(-k_n \frac{|z|}{a}\right). \qquad (5) $$

With the stream function, the radial velocity component

$$ \frac{u_r(r, z)}{U} = 2\tau\beta^2 \sum_{n=1}^{\infty} \frac{J_1(k_n\beta)}{J_0^2(k_n)} J_1\left(k_n \frac{r}{a}\right) \exp\left(-k_n \frac{|z|}{a}\right) \qquad (6) $$

and the axial velocity component

$$ \frac{u_z(r, z)}{U} = 1 + 2\tau\beta^2 \sum_{n=1}^{\infty} \frac{J_1(k_n\beta)}{J_0^2(k_n)} J_0\left(k_n \frac{r}{a}\right) \exp\left(-k_n \frac{|z|}{a}\right) \qquad (7) $$

are given. The velocity field (equation 6 and 7) takes the induction of the vortex ring into account. At the stagnation point $z = \pm z_0$, the axial velocity $u_z$ vanishes for

$$ 1 = -2\tau\beta^2 \sum_{n=1}^{\infty} \frac{J_1(k_n\beta)}{J_0^2(k_n)} J_0(k_n) \exp\left(-k_n \frac{|z_0|}{a}\right). \qquad (8) $$
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Contour plots of the stream function (equation 5) for a vortex ring with the strength $τ = 0.8$ (left), $τ = 2.5$ (middle), $τ = 5.0$ (right) for $β = 0.8$.</figcaption></figure>

Equation 8 is implicit for the stagnation point $z_0 = z_0(τ, β)$. Figure 3 shows Stokes stream function and the stagnation points at $z = \pm z_0$ for different circulation strengths $τ$ of the vortex for $β = 0.8$.

## 4. THE VORTEX SYSTEM AT HEAVY OVERLOAD

<figure><img src="image_6.png" /><figcaption>**Figure 4.** The vortex system of an axial turbomachine at overload, $Z = 1$.</figcaption></figure>

Figure 4 visualizes the vortex system of an axial turbomachine composed of $Z = 1$ impeller blade at heavy overload $\varphi \to \infty$ ($Z = 1$ is chosen to improve clarity only). For this limiting case, the flow at cross section A-A far downstream of the machine is a plane potential flow. It can be described using theory of functions (complex analysis) [22]. Mirrored tip vortices are necessary, to fulfill the kinematic boundary condition on the tube wall. These mirrored vortices are located in the housing and on the rotational axis of the turbomachine. Considering the tip vortex and its mirrored conjugates only, i.e. neglecting the hub vortex as a first step, one obtains the system visualized in figure 4 bottom left. The vortex on the axis and the hub vortex feature identical magnitude but opposed rotating direction. Adding the hub vortex, yielding the complete system, hence, results in the annulation of these two vortices; see figure 4 bottom right.

In the following notation, the complex coordinates $\zeta = \xi + i\eta = r \exp(i\theta)$ are composed of a real part $\xi = r \cos \theta$ and an imaginary part $\eta = r \sin \theta$, with $i = \sqrt{-1}$. The complex potential $F(\xi, \eta)$ is divided into the real part, which is the velocity potential $\Re[F(\xi, \eta)] = \Phi(\xi, \eta)$ and the imaginary part, which is the stream function $\Im[F(\xi, \eta)] = \Psi(\xi, \eta)$. $\bar{F}$ is the complex conjugate of $F$.

For the considered potential flow, the tip vortex at radial position $b = (1 - s)a$ yields the complex potential

$$F_1(\zeta) = -\frac{i\Gamma}{2\pi} \ln(\zeta - b). \qquad (9)$$

Here, $s$ is the dimensionless gap. The Milne-Thomson circle theorem [23] is applied, to derive the complex potential satisfying the kinematic boundary condition at the wall. This theorem postulates a resulting complex potential

$$F_2(\zeta) = F_1(\zeta) + \bar{F}_1 \left( \frac{a^2}{\zeta} \right) \qquad (10)$$
for a potential $F_1$ and the mirrored potential at the surrounding wall. Adding the potential of the mirrored tip vortex on the axis of the turbomachine $\zeta = 0$, see figure 4 bottom right, yields for the complex flow potential

$$F'(\zeta) = -\frac{i\Gamma}{2\pi} \left[ \ln(\zeta - b) - \ln\left(\zeta - \frac{a^2}{b}\right) + \ln\zeta \right] + \text{const.} \quad (11)$$

The tip vortex at $\zeta = b$ with the circulation $\Gamma$ necessitates a mirrored vortex at $\zeta = 0$ with the same magnitude of circulation and a mirrored vortex in the housing at $\zeta = a^2/b = a/(1-s)$ with the same magnitude and inverted direction. Up to now, the hub vortex is excluded from the considerations. Considering the hub vortex, as visualized in figure 4 bottom right, yields for the complex potential

$$F(\zeta) = -\frac{i\Gamma}{2\pi} \left[ \ln(\zeta - b) - \ln\left(\zeta - \frac{a^2}{b}\right) \right] + \text{const.} \quad (12)$$

In the following, this analysis shows, that an induced movement of the gap vortex occurs against the rotating direction of the turbomachine at heavy overload. A potential vortex induces a velocity on the surrounding flow. The velocity components of a given potential $F(\zeta)$ are calculated by

$$\frac{dF(\zeta)}{d\zeta} = u_{\text{ind}}\vec{e}_{\zeta} - i v_{\text{ind}}\vec{e}_{\eta}. \quad (13)$$

A straight vortex filament does not induce a velocity on its own, due to the Biot-Savart law, so the induced velocity at $\zeta = b$ is only due to the mirrored tip vortex at $\zeta = a/(1-s)$. The resulting induced velocity at the position of the tip vortex yields

$$\begin{aligned} u_{\text{ind}}\vec{e}_{\zeta} - i v_{\text{ind}}\vec{e}_{\eta} &= \frac{d}{d\zeta} \left[ \frac{i\Gamma}{2\pi} \ln\left(\zeta - \frac{a^2}{b}\right) \right] \Big|_{\zeta=b} \\ &= -\frac{i\Gamma}{2\pi b} \frac{b^2}{a^2 - b^2} \vec{e}_{\eta}. \end{aligned} \quad (14)$$

Assuming a turbomachine with Z impeller blades, the rotating velocity of the tip vortex is

$$\begin{aligned} v_{\text{ind}} &= \frac{Z\Gamma}{2\pi b} \frac{b^2}{a^2 - b^2} = \frac{Z\Gamma}{2\pi a s} \frac{1-s}{2-s} \\ &= \Omega a \frac{1 - \varphi/\dot{\varphi}}{s} \frac{(1-s)^3}{2-s}. \end{aligned} \quad (15)$$

This is the rotating direction against the rotating direction of the turbomachine. For symmetry reasons, the rotating trajectory defines a circular path at radius $b = a(1-s)$. Hence, the induced frequency at $\zeta = b$ is

$$\Omega_{\text{ind}} = \frac{v_{\text{ind}}}{b} = \frac{\Omega}{s} \frac{(1-s)^2}{2-s} \left(1 - \frac{\varphi}{\dot{\varphi}}\right). \quad (16)$$

For high flow number $\varphi \to \dot{\varphi}$ and small gap $s \ll 1$, the induced frequency yields

$$\frac{\Omega_{\text{ind}}}{\Omega} = \frac{1}{2s}, \quad \text{for } \varphi \to \infty \quad (17)$$

against the rotating direction of the turbomachine.

# 5. APPLICATION FOR ACOUSTICAL INVESTIGATIONS

<figure><img src="image_7.png" /><figcaption>Figure 5. Frequency spectra of the tip clearance noise depending on the flow number $\varphi$ [4].</figcaption></figure>

Previous investigations by Karstadt et al. [4] analysed the tip clearance noise in axial turbomachines. Figure 5 shows the frequency spectra over the complete operating range. Remarkable are the peaks at 42 Hz and 375 Hz, which correspond to the rotational speed $n$ and the blade passing frequency $Zn$. Fukano and Yang [17] show, that the circumferential frequency of the tip clearance noise shifts to lower values with increasing tip clearance $s$ and decreasing flow number $\varphi$ due to the larger extent of the gap vortex. The present paper investigates the frequency of the tip clearance noise depending on the operating point applying analytical methods. At heavy part load, equation 2 indicates that an induced frequency at half the rotation speed of the turbomachine should appear. This frequency was also observed by Karstadt et al. [4] as figure 5 shows a high intensity in the region of 20 Hz. Equation 17 indicates, that for heavy overload, the induced frequency will increase with decreasing tip clearance. Furthermore, we expect a noise of high frequency due to the small value of $s < 1\%$, which is common for turbomachines. The broadband drop in the sound power for all flow numbers is clearly visible.

Müller [24] applied the continuity and the momentum equation and deduced, that sound inside a fluid volume is only emitted, if the rotation of the velocity field changes in time. The present study applies a similar approach to analyse the tip clearance noise of a turbomachine. Time-consuming simulations, as performed by Carolus et al. [25] surely allow a more profond and accurate insight into the acoustics of turbomachines. The development of an analytical model, which predicts main frequencys, is yet interesting to generate a deeper understanding of the acoustics in turbomachines.

These findings and the presented analytical model in this paper could be an efficient tool for acoustic design of turbomachines.

# 6. SUMMARY AND CONCLUSION

An interplay between dynamic and kinematic effects explains flow structures and phenomena. Using computational fluid dynamics, a clear distinction of both effects is often impossible. In contrast, analytical methods allow a more focused picture of fluid mechanics, i.e. they allow a clear distinction of effects. Of course, only generic flows are accessible to analytical methods.
This paper focused on an analytic model for wall stall, so far being explained by dynamics only: boundary layer separation is indeed a dynamic effect. Nevertheless boundary layer separation is not necessarily the only reason for wall stall. It is shown, that kinematics may also explain at least some effects of wall and rotating stall. The used picture for a flow at small flow numbers is a thread spool, rolling up the tip vortices resulting from rotating bound vortices. From the fluid mechanics perspective, the thread spool is a coaxial vortex ring of increasing strength connected to a semi-infinite hub vortex (figure 1).

So far, the velocity potential of a coaxial vortex ring inside a tube was unknown. The solution of Laplace's equation results in the velocity potential for the vortex filament within a tube [10]; see equation 4.

This paper gained three main results, which are due to kinematics only. First, at part load operation, the hub vortex induces a sub-synchronous rotation of the vortex ring. The derived rotational speed $\Omega_{ind} = 0.5\,\Omega$ of the vortex ring is surprisingly consistent with observed sub-synchronous speeds of rotating stall; cf. [11]. Second, the vortex ring induces an upstream axial velocity at the wall. Together with the undisturbed flow velocity, this results in a stagnation point upstream and downstream at the wall, which may be interpreted as "wall stall" (figure 3). Third, at overload operation, the induced rotational direction is inverted to the case at part load. The semi-infinite straight vortex filament at the outer blade end rotates against the rotating direction of the turbomachine, due to the induction of the hub vortex. The induced frequency yields $\Omega_{ind} = -0.5\,\Omega/s$.

The presented analytical model may give new arguments and improves the understanding of the vortex system in turbomachines, but is also intended to motivate generic experiments. Hence, a test rig will validate the models presented in this paper in the near future.

As a next step, the velocity potential for a coaxial vortex ring filament in a circular tube (equation 4) will be extended to a coaxial vortex layer, yielding a transient behavior of the vortex system. This behavior leads to a change in the circulation over time being responsible for noise emission [24].

## REFERENCES

[1] E. M. Greitzer. Review-Axial Compressor Stall Phenomena. *Journal of Fluids Engineering*, 102(2):134-151, June 1980.

[2] S. L. Dixon. *Some Three Dimensional Effects of Rotating Stall*. HM Stationery Office, 1962.

[3] F.-J. Cloos, D. Stapp, and P. F. Pelz. Swirl boundary layer and flow separation at the inlet of a rotating pipe. *Journal of Fluid Mechanics*, 811:350-371, January 2017.

[4] S Karstadt, B Matyschok, and P. F. Pelz. Sound Deadening on Fans. In *ASME Turbo Expo 2011: Turbine Technical Conference and Exposition*, pages 961-971, Vancouver, Canada, 2011. American Society of Mechanical Engineers.

[5] T. Zhu. *On the Flow Induced Tip Clearance Noise in Axial Fans*. PhD thesis, Universität Siegen, Siegen, 2016.

[6] A. Betz. Schraubenpropeller mit geringstem Energieverlust. Mit einem Zusatz von L. Prandtl. *Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse*, pages 193-217, 1919.

[7] S. Goldstein. On the vortex theory of screw propellers. *Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character*, 123(792):440-465, 1929.

[8] H. Glauert. *Die Grundlagen der Tragflügel- und Luftschraubentheorie*. Springer, Berlin, 1929.

[9] G. van Kuik. The relationship between loads and power of a rotor and an actuator disc. *Journal of Physics: Conference Series*, 555:012101, December 2014.

[10] P. F. Pelz and P. Taubert. Vortex induced transient stall. *submitted for publication to Archive of Applied Mechanics*, 2017.

[11] C. E. Brennen. *Hydrodynamics of pumps*. Oxford science publications. Concepts ETI; Oxford University Press, Norwich, Vt., USA: Oxford; New York, 1994.

[12] H. Helmholtz. Über Integrale der hydrodynamischen Gleichungen welche den Wirbelbewegungen entsprechen. *Journal für die reine und angewandte Mathematik*, 55:25-55, 1858.

[13] N. Didden. On the formation of vortex rings: rolling-up and production of circulation. *Zeitschrift für Angewandte Mathematik und Physik (ZAMP)*, 30(1):101-116, 1979.

[14] R. E. Longhouse. Control of tip-vortex noise of axial flow fans by rotating shrouds. *Journal of sound and vibration*, 58(2):201-214, 1978.

[15] T. F. Brooks and M. A. Marcolini. Airfoil tip vortex formation noise. *AIAA Journal*, 24(2):246-252, 1986.

[16] R. Mailach, I. Lehmann, and K. Vogeler. Rotating Instabilities in an Axial Compressor Originating From the Fluctuating Blade Tip Vortex. *Journal of Turbomachinery*, 123(3):453, 2001.

[17] T Fukano and C.-M Jang. Tip clearance noise of axial flow fans operating at design and off-design condition. *Journal of Sound and Vibration*, 275(3-5):1027-1050, 2004.

[18] N. M. McDougall, N. A. Cumpsty, and T. P. Hynes. Stall inception in axial compressors. In *ASME 1989 International Gas Turbine and Aeroengine Congress and Exposition*, pages 116-125, Toronto, Canada, 1989. American Society of Mechanical Engineers.

[19] H. Saathoff and U. Stark. Tip Clearance Flow Induced Endwall Boundary Layer Separation in a Single-Stage Axial-Flow Low-Speed Compressor. In *ASME Turbo Expo 2000: Turbine Technical Conference and Exposition*, Munich, Germany, May 2000. ASME.

[20] S. Karstadt and P. F. Pelz. A Physical Model for the Tip Vortex Loss: Experimental Validation and Scaling
Method. In *ASME Turbo Expo 2012: Turbine Technical Conference and Exposition*, pages 73–81, Copenhagen, Denmark, 2012. American Society of Mechanical Engineers.

[21] L. Prandtl. Über die Entstehung von Wirbeln in der idealen Flüssigkeit, mit Anwendung auf die Tragflügeltheorie und andere Aufgaben. In *Vorträge aus dem Gebiete der Hydro-und Aerodynamik*, pages 18–33. Springer, Innsbruck, 1924.

[22] P. G. Saffman. *Vortex Dynamics*. Cambridge University Press, 1992.

[23] L. M. Milne-Thomson. *Theoretical hydrodynamics*. Courier Corporation, 1968.

[24] E. A. Müller. Der Wirbel als Schallerzeuger. In *Ernst Becker Gedächtnis Kolloquium*, 28, pages 137–151. Schriftenreihe Wissenschaft und Technik, TH Darmstadt, Fachbereich Mechanik, 1985.

[25] T. Carolus, M. Schneider, and H. Reese. Axial flow fan broad-band noise and prediction. *Journal of Sound and Vibration*, 300(1-2):50–70, February 2007.