---
title: Rankine Oval Dynamics
url: https://www.emergentmind.com/topics/rankine-oval
type: topic
---

# Rankine Oval Dynamics

The Rankine oval is the closed stagnation streamline obtained by superimposing a uniform stream of speed \(U_\infty\) with a source–sink pair of equal and opposite strength \(m\), separated by a distance \(2a\). In classical potential flow it is an exact irrotational construction and is therefore “streamlined” in the inviscid sense; in viscous flow, however, the same geometry behaves as a bluff body, exhibiting boundary-layer separation, recirculation, wake instability, and eventually periodic vortex shedding. A systematic direct numerical simulation study of incompressible flow past Rankine ovals over \(Re=10\) to \(200\) and \(Ua/m=0\) to \(1\) establishes how the shape parameter \(Ua/m\) controls the transition from steady to periodic wakes, the scaling of drag and shedding frequency, and the extent to which potential-flow information remains useful in viscous drag prediction [2509.07192].

## 1. Classical definition and exact construction

In the standard construction, a source of strength \(+m\) is placed at \((x,y)=(-a,0)\), a sink of strength \(-m\) at \((x,y)=(a,0)\), and a uniform stream of velocity \(U_\infty\) is imposed in the \(x\)-direction. The body boundary is the streamline \(\psi=0\), with streamfunction
\[
\psi =\frac{m}{2\pi}{\theta _2}-\frac{m}{2\pi}{\theta_1}+{U_\infty}y.
\]

In Cartesian form, the same boundary is written as
\[
\frac{m}{2\pi}\tan^{-1}\left(\frac{y}{x+a}\right)-\frac{m}{2\pi}\tan^{-1}\left(\frac{y}{x-a}\right)+U_\infty y=0.
\]

The paper also gives an implicit dimensionless relation. The typesetting in the manuscript is corrupted, but the intended form is reported as
\[
\left(\frac{x}{a}\right)^2+\left(\frac{y}{a}\right)^2-1 = \frac{2(y/a)}{\tan\!\left(2\pi \frac{U_\infty a}{m}\frac{y}{a}\right)}.
\]

The governing geometric control parameter is the nondimensional quantity
\[
\frac{U_\infty a}{m},
\]
abbreviated as \(Ua/m\). Physically, it is the inverse of a nondimensional source strength. Small \(Ua/m\) corresponds to stronger source–sink influence relative to the freestream and produces fuller, more circular bodies; large \(Ua/m\) corresponds to stronger freestream influence and produces longer, more slender ovals [2509.07192].

The Rankine oval is therefore not merely a descriptive shape class but an exact streamline body generated by a specific potential-flow superposition. This exactness is central to its historical role in inviscid theory and to its later reinterpretation as a viscous bluff-body problem.

## 2. Shape parameterization, normalization, and limiting forms

The half-length \(l\) is the distance from the center to the fore and aft stagnation points \((\pm l,0)\), and the half-height \(h\) is the maximum half-width at \((0,\pm h)\). The manuscript gives the useful relation
\[
l = a\sqrt{1+\frac{1}{\pi C}}, \qquad C=\frac{Ua}{m}.
\]
It further states that \(h\) is determined from the \(x=0\) boundary condition and that \(a\), \(h\), and \(C\) are linked implicitly.

The aspect ratio is defined as
\[
AR=l/h.
\]

To compare shapes consistently, all ovals are rescaled so that
\[
h=0.5,\qquad D=2h=1.
\]
This normalization fixes the vertical diameter across the full family, so Reynolds number and force coefficients are evaluated on a common geometric basis.

Within this normalization, the geometry over \(0\le Ua/m\le 1\) follows a nearly linear trend. At \(Ua/m=0\), the body is a circle with \(l=0.500\), \(h=0.500\), and \(AR=1.000\). At \(Ua/m=1\), the body has \(l=1.497\), \(h=0.500\), and \(AR=2.994\), i.e. an oval about three diameters long. The approximation
\[
\frac{l}{h} \approx 1 + 2\frac{Ua}{m}
\]
is derived and used over the studied range [2509.07192].

| \(Ua/m\) | \(l/h\) trend | Geometric interpretation |
|---|---:|---|
| \(0\) | \(1.000\) | Circular cylinder |
| \(1\) | \(2.994\) | Slender oval, about three diameters long |

Two limiting statements organize the family. First, \(Ua/m=0\) corresponds to a circular cylinder. Second, increasing \(Ua/m\) produces progressively more elongated Rankine ovals. This makes the family particularly useful for studying shape effects continuously, without changing the underlying analytic construction.

## 3. Inviscidly streamlined, viscously bluff

In inviscid potential flow, the Rankine oval is “streamlined” because it is literally defined by a smooth streamline of an irrotational velocity field. In that setting there is no boundary layer, no separation, and no wake in the viscous sense. The pressure field follows directly from the potential solution, and the body is an analytically convenient idealization.

In viscous flow, by contrast, the same shape behaves as a bluff body. The no-slip boundary condition generates boundary layers, separation, recirculating wake bubbles, and eventually periodic vortex shedding. The central conceptual point is that streamline geometry in potential flow does not prevent wake instability once viscosity is included [2509.07192].

This establishes a bridge between two classical viewpoints. On one side is potential-flow theory, in which Rankine ovals are canonical exact shapes defined by source–sink superposition. On the other side is practical wake dynamics, in which the same bodies exhibit the standard separated-flow sequence seen behind cylinders, ellipses, rectangles, and other bluff bodies: steady recirculation at low \(Re\), followed by a Hopf bifurcation to periodic shedding.

A common misconception is that a body derived from an inviscid streamline should remain dynamically streamlined in real flow. The reported results contradict that inference. A plausible implication is that analytic streamline constructions and viscous wake stability belong to different problem classes: the former constrains the outer inviscid geometry, whereas the latter is governed by no-slip-induced boundary layers and their instability.

## 4. Governing equations, simulation framework, and onset criterion

The viscous-flow study is based on two-dimensional direct numerical simulations of incompressible flow past Rankine ovals at zero angle of attack. The governing equations are the incompressible Navier–Stokes equations,
\[
\frac{\partial{u}_i}{\partial {x_i} = 0,
\]
\[
\frac{\partial {u}_i}{\partial {t} + {u}_j\frac{\partial {u}_i}{\partial {x_j} = - \frac{\partial p}{\partial x_i} + \nu \frac{\partial ^2 {u}_i}{\partial x_j\partial x_j}.
\]
The simulations are 2D, so \(i=1,2\), corresponding to \(x\) and \(y\).

The solver uses a cell-centered finite volume method on hybrid unstructured meshes, a second-order Roe scheme for convection, a reconstructed central scheme for viscous flux, and second-order implicit dual-time stepping for unsteady simulations. The Reynolds number is defined by
\[
Re = \frac{U_\infty D}{\nu},
\]
with \(D=1\) under the chosen normalization. The study also defines
\[
St = \frac{fD}{U_\infty}, \qquad Cl = \frac{F_l}{0.5\rho U_\infty^2D}, \qquad Cd = \frac{F_d}{0.5\rho U_\infty^2D},
\]
\[
Cp = \frac{p - p_\infty}{0.5\rho U_\infty^2},
\]
and decomposes drag as
\[
Cd = Cd_p + Cd_f.
\]

The simulations cover \(Re=10\) to \(200\) and \(Ua/m=0\) to \(1\), with finer sampling near shedding onset for each shape. The computational domain is
\[
[-100D,100D]\times[-100D,100D],
\]
with maximum blockage ratio \(0.005\). Boundary conditions are free-slip at the far field and no-slip on the Rankine oval surface. The first near-wall grid spacing is less than \(0.01D\), with growth rate \(1.05\), and total cell count is about \(7.8\times10^4\) to \(1.2\times10^5\), depending on geometry [2509.07192].

Validation is performed against standard cylinder data at \(Re=100\) and \(200\). For the circular-cylinder onset, the study obtains
\[
Re_c = 46.9,\qquad St_c=0.115,
\]
which is reported to match the literature well.

The onset of vortex shedding is identified using Stuart–Landau/Hopf-bifurcation analysis. The amplitude equation is
\[
\frac{dA}{dt} = \sigma A - L|A|^2 A,
\]
where \(A(t)\) is taken from the lift-coefficient envelope and \(\sigma=\sigma_r+i\sigma_i\). The onset criterion is based on the real growth rate \(\sigma_r\), using
\[
\frac{1}{|A|}\frac{d|A|}{dt} = \sigma_r\left(1 - \frac{|A|^2}{|A|_{max}^2}\right).
\]
The critical Reynolds number is then determined by linearly fitting \(\sigma_r\) versus \(Re\) and locating \(\sigma_r=0\). This procedure directly ties the instability threshold to Hopf-bifurcation dynamics rather than to a purely visual wake criterion.

## 5. Wake instability, critical Reynolds number, and wake organization

The principal instability result is that the critical Reynolds number increases approximately linearly with \(Ua/m\). The reported empirical fit is
\[
Re_c = 60\,Ua/m + 45,
\]
stated to be reliable over
\[
0.1 \le Ua/m \le 1.0,
\]
with minor deviation near \(Ua/m=0\). Numerically, \(Re_c \approx 46.9\) at \(Ua/m=0\) and \(Re_c \approx 104.0\) at \(Ua/m=1.0\) [2509.07192].

The physical interpretation given is that increasing \(Ua/m\) makes the body more slender and more streamlined geometrically, which tends to lengthen the recirculation region and stabilize the wake, so a higher \(Re\) is required to trigger the Hopf bifurcation. Fuller, more circle-like shapes therefore shed earlier.

Near \(Re_c\), the wake already exhibits the canonical anti-symmetric development associated with vortex shedding. At the phase of maximum lift, a single distinct vortex bubble is observed; when lift is zero, two wake vortex bubbles are visible. The asymmetry grows gradually as the wake approaches periodic shedding.

The study compares this trend with other two-dimensional bluff bodies. Rankine ovals and rectangles show roughly linear increases of \(Re_c\) with aspect ratio, whereas ellipses and diamonds show more nonlinear behavior. This suggests that Rankine ovals form an intermediate class: they are analytically generated by potential flow, yet their viscous instability thresholds vary with shape in a manner characteristic of separated bluff-body wakes.

The wake patterns are grouped into five classes based on vorticity intensity and shedding frequency. Increasing \(Re\) and decreasing \(Ua/m\) both intensify vorticity and increase shedding frequency. Using rms velocity maps, the wake width \(W\) is defined as the distance between the two rms peaks, and the vortex formation length \(L_f\) as the streamwise distance from the body center to the midpoint between those peaks. At fixed \(Re\), increasing \(Ua/m\) increases both \(L_f\) and \(W\); at fixed \(Ua/m\), increasing \(Re\) decreases both. At \(Re=100\), for example,
\[
L_f = 2.73 \text{ for } Ua/m=0.1,\qquad L_f = 15.62 \text{ for } Ua/m=0.9.
\]
At \(Re=200\), \(W\) becomes nearly shape-independent, around \(0.86\)–\(0.90\).

## 6. Forces, Strouhal scaling, and reduced-order drag prediction

At the critical point, increasing \(Ua/m\) decreases the critical mean drag coefficient from about \(1.34\) to \(1.01\) and increases the critical Strouhal number slightly from about \(0.12\) to \(0.13\). The approximate invariance
\[
\overline{Cd}_c\, St_c \approx 0.14 \pm 0.01
\]
is reported as consistent with previously observed bluff-body trends [2509.07192].

Over \(Re\in[50,200]\), the mean drag decreases with increasing Reynolds number for every shape, and at fixed \(Re\) it decreases with increasing \(Ua/m\). The fit
\[
\overline{Cd} = a Re^{-1/2} + b,
\]
with
\[
b \approx \overline{Cd}_p,
\]
implies
\[
\overline{Cd} = a Re^{-1/2} + \overline{Cd}_p,
\]
so the friction contribution scales approximately like \(Re^{-1/2}\). The decomposition further shows that \(\overline{Cd}_p\) is nearly constant with \(Re\), whereas \(\overline{Cd}_f\) decreases as \(Re\) increases. At \(Ua/m=0.5\), for example, \(\overline{Cd}_p\approx 0.60\) over the studied range, while \(\overline{Cd}_f\) drops from \(0.52\) to \(0.35\). For sufficiently elongated shapes, friction drag can become comparable to or exceed pressure drag; at \(Ua/m=0.9\) and \(Re=100\),
\[
\overline{Cd}_f = 0.53 > \overline{Cd}_p = 0.50.
\]

The instantaneous drag fluctuation is dominated by pressure-drag fluctuation, while friction drag remains comparatively steady. At \(Ua/m=0.3\), the amplitude of \(Cd_p\) is \(0.007\), compared with \(0.001\) for \(Cd_f\). As \(Ua/m\) increases, the amplitude of pressure-drag oscillation drops sharply, from \(0.042\) at \(Ua/m=0.1\) to \(0.001\) at \(Ua/m=0.9\). At fixed \(Re=200\), increasing \(Ua/m\) reduces the amplitude of both lift and drag oscillations: the maximum instantaneous \(Cl\) falls from \(0.51\) at \(Ua/m=0.1\) to \(0.13\) at \(Ua/m=0.9\), and drag oscillation amplitude falls from \(0.0235\) to \(0.0004\).

The Strouhal number shows two distinct descriptions. On the one hand, in the raw geometric parameterization, \(St\) increases with \(Re\) for a given shape and decreases with \(Ua/m\) at fixed \(Re\); by \(Re=200\), \(St\approx 0.18\) for all shapes, indicating weaker geometry dependence in that regime. At \(Re=100\),
\[
St \approx 0.16 \text{ at } Ua/m=0.1,\qquad St \approx 0.13 \text{ at } Ua/m=0.9.
\]
The study also notes that once \(Re\) exceeds \(Re_c\), \(St\) appears near \(\approx 0.13\) rather than growing gradually from zero.

On the other hand, the data-driven dimensional analysis identifies reduced variables in which shape dependence largely collapses. For friction drag, the selected variables are \(\rho,\mu,l,a,U,m\), with half-length-based coefficient
\[
\overline{Cd}_{fl}=\frac{F_d^{fric}{0.5\rho U_\infty^2l}.
\]
The dimensional-analysis form is
\[
\overline{Cd}_{fl}=f \left(\frac{\rho l U}{\mu},\frac{Ua}{m},\frac{l}{a} \right),
\]
and the dominant controlling group is essentially
\[
Re_l=\frac{\rho l U}{\mu}.
\]
The reported conclusion is that \(\overline{Cd}_{fl}\) is determined by \(Re_l\) independently of \(Ua/m\), with fit
\[
\overline{Cd}_{fl} = 4.0 Re_l^{-0.54}.
\]

For vortex shedding, the selected wake variables are the vortex formation length \(L_f\) and the largest reverse-flow speed \(U_r\), yielding
\[
Re_r = \frac{\rho U_r L_f}{\mu}.
\]
The paper states that \(St\) is independent of \(Ua/m\) and controlled only by \(Re_r\), with fitted relation
\[
St = -0.54 Re_r^{0.54} + 0.27.
\]
The sign in the printed formula is identified in the source text as likely typographically problematic, because the plotted trend and accompanying discussion indicate a meaningful positive correlation over the sampled range. The explicit claim retained by the study is that the dependence collapses onto a single law in terms of \(Re_r\), essentially independent of the shape parameter.

A further reduced-order result concerns pressure drag estimation. At \(Re=200\), the time-averaged surface pressure distributions from DNS approach inviscid potential-flow predictions as \(Ua/m\) increases. This motivates the piecewise approximation
\[
C_{p, \text{visc}(x)} = \begin{cases} C_{p, \text{pot}(x)}, & \text{if } x < 0 \\ C_{p, \text{pot}(0)}, & \text{if } x \geq 0 \end{cases}
\]
for which the intended expression is explicitly reconstructed in the source text. The approximation works poorly for short, full bodies, with pressure-drag error exceeding \(50\%\) at \(Ua/m=0.1\) and still \(32.5\%\) at \(Ua/m=0.5\), but becomes accurate for sufficiently slender ovals; at \(Ua/m=0.9\), the estimated pressure drag is \(0.56\), within \(5.6\%\) of DNS [2509.07192].

Combining this pressure estimate with
\[
\overline{Cd}_{fl} = 4.0 Re_l^{-0.54},
\]
the total drag for sufficiently large \(Ua/m\) can be predicted without full numerical simulation by obtaining \(Cd_p\) from potential flow plus rear correction, computing \(l\) from \(Ua/m\), evaluating \(Re_l\), obtaining \(Cd_{fl}\), converting to \(Cd_f\), and then summing \(Cd_p + Cd_f\). This suggests a limited but practically useful persistence of potential-flow structure within a viscous separated-flow regime.

The Rankine oval thus occupies a distinctive position among canonical bodies. It is an exact streamline body in inviscid theory, yet in viscous flow it falls squarely within the bluff-body instability framework. The contrast between its elegant potential-flow construction and its separated, shape-sensitive wake dynamics is the defining insight of the recent viscous-flow study.

Source: https://www.emergentmind.com/topics/rankine-oval