0
Research Papers: Fundamental Issues and Canonical Flows

# Experiments on the Effects of Reynolds Number and Advance Ratio on the Unfolding of Disorganization in Low-Speed Underwater Propulsors With Vibrating BladesOPEN ACCESS

[+] Author and Article Information

Fellow ASME
Naval Undersea Warfare Center Division,
Newport 1176 Howell Street,
Newport, RI 02841

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received May 17, 2017; final manuscript received September 21, 2017; published online November 17, 2017. Assoc. Editor: Riccardo Mereu. This material is declared a work of the U.S. Government and is not subject to copyright protection in the United States. Approved for public release; distribution is unlimited. License CC-BY 4.0.

J. Fluids Eng 140(4), 041202 (Nov 17, 2017) (12 pages) Paper No: FE-17-1283; doi: 10.1115/1.4038241 History: Received May 17, 2017; Revised September 21, 2017

## Abstract

Ships and submarines are acoustic hazards to marine life. The rational control of acoustic radiation would be possible at least at low Reynolds numbers if the underlying organization buried in seeming randomness is revealed. We build a novel low-speed propulsor where all blades undergo small-amplitude pitch oscillation while spinning at large pitch angles at transitional chord Reynolds numbers (3.75 × 103 ≤ Rec ≤ 3.75 × 104) and advance ratios (0.51 ≤ J ≤ 4.89). We measure and model time-averaged and temporal thrust. The relationship between the time-averaged and the temporal thrust is observed when the latter is mapped as limit cycle oscillation (LCO), or departure from it. High-thrust coefficients occurring at large (30 deg and 45 deg) angles of amplitude of blade vibration are modeled assuming poststall lift enhancement due to flapping blades when a leading edge vortex (LEV) forms, while the lower thrust coefficients occurring at 20 deg are modeled by its absence. The disorganization in temporal thrust increases with J and Rec. An external orthogonal oscillator, perhaps a vibration, is modeled to couple with the thrust oscillator for temporal control of disorganization. The unfolding disorganization is seen as a departure from LCO, and it is attenuated by smooth-wall boundary-layer fencing, compared to unfenced smooth and rough surfaces. When the fencing properties of the leading edge tubercles of whale fins are recognized, the ratio of the spacing of the fences and chord is found to be similar (0.5–1.0) in both whale flippers and aircraft wings.

<>

## Introduction

The world is becoming acoustically noisier. For example, due to the rapid rise in shipping, the low-frequency ambient acoustic noise in the ocean is increasing at an average rate of 0.5 dB per year [1]. Such noise occurs at high Reynolds numbers whose hallmark is seeming disorganization that is randomness which is thought to be statistical and not deterministic.

Experiments in turbulent boundary layers in the momentum thickness Reynolds number range of 500–17,500 have shown that organized hairpin vortex structures populate the boundary layer when viewed in orthogonal cross-stream planes one of them being slanted at 45 deg to the flow direction [26]. Therefore, organizations can remain masked in seeming disorganization unless viewed in appropriate orientations. In a recent dynamical analysis of such turbulent boundary layers, the experimentally determined spatiotemporal distributions of several state variables have been reproduced with remarkable accuracy but only at transitional low Reynolds numbers so far [7]. Furthermore, it was possible to reproduce accurately the dermal denticles and riblets of sharks [8] as devices for control of lateral diffusion of vorticities (called μ-control). This work indicates that: (1) the spatiotemporal evolution of the flow at least at low transitional Reynolds numbers is deterministic and not random and (2) lateral control of vorticity diffusion organizes the boundary layer turbulence at least in the wall-layers.

Therefore, in the context of a propulsor before we march to higher Reynolds numbers, first we need to pinpoint exactly when and how disorganization in thrust emerges after laminar to turbulence transition has already taken place. To simplify, in a propulsor, we consider temporal growth rather than spatiotemporal growth of the flow. We also explore devices that control lateral (spatial) diffusion of vorticity in lifting surfaces in vogue in nature and engineering.

Assuming a propulsor to be a point source of acoustic radiation related to the time derivative ($˙$) of fluctuating ($′$) pressure ($p′˙$), we write $p′˙∝Fx′˙/A$, where $Fx′$ is the fluctuating component of thrust Fx, and $A$ is the cross-sectional surface area. To reduce error, instead of time derivative, we map temporal history ($Fx′(t)$), where t is the time. In this work, we look for the emergence of disorganization in $Fx′$ at low Reynolds numbers and advance ratios and the properties of their control.

Nonlinear time series analyses have led to chaos predictions [9]. Such approaches are useful when the underlying determinism cannot be extracted from the noise, and the initial condition is unknown. However, since chaotic phenomena are basically deterministic although complex, if the system is oscillatory and if we have the option of starting from the pristine (unchaotic) initial condition, then it may be possible to uncover the route to thrust disorganization. This is the approach we take in this work on an underwater rotating propulsor.

To understand the measurements, one could fully solve the Navier–Stokes equations at low Reynolds numbers, or carry out turbulence modeling. Instead, we have taken the dynamical approach, because the propulsor is oscillatory. Furthermore, the nonlinearly coupled blades and wakes tend to be self-regulating. These are key features of many fluid dynamic systems [7,10]. We hope that these features would offer preferential properties to the flow thereby simplifying our search for organization. However, in dynamical analysis considerable insight of the flow mechanism and knowledge of key parameter ratios and of the interactions between the main and control oscillators are needed. For these reasons in engineering, dynamical analysis is not routinely used. However, in view of ours and others recent successes in oscillatory propulsion [7,1012], here, we use the dynamical approach. We had success even in low Reynolds number turbulent boundary layers, where the oscillatory nature of production is masked in seeming randomness [7]. Although the following fluid dynamic example is distant, we are encouraged to notice that dynamical meteorological modeling is now supplanting statistical models in annual rainfall predictions in large geophysical land masses such as the U.S. and India indicating their maturity and acceptance for planning.

In dynamical determination of the emergence of disorganization in thrust, one normally would plot ($Fx′$ − $Fx′˙$) where $Fx′˙$ represents the orthogonal force component, that is the relationship between roll (or spinning) and pitching motions. However, in the present propulsor, each blade has its own pitch axis and such graphing would be hard to interpret. Instead, we found that at the lowest values of $(Rec,J)$, limit cycle oscillations (LCOs) are present in axes ($Fx′-I$), where I is the total current into the propulsor which is proportional to rotational speed in linear motors (current consumption due to pitch vibrations and electronics is small). Here, $Fx′$ is the thrust fluctuation in the useful direction of motion, and I is the total current input, meaning ($Fx′-I$) is a relationship between input and output. We investigate the departures from such LCOs.

In underwater propulsion, the significance of viewing flows as oscillatory systems has relevance to control. The flows can be described using Stuart–Landau equations of wake circulation [7,10], or van der Pol equations [11,12]. Such equations denote the presence of self-regulation, and their properties are discussed in Refs. [1315]. It turns out that all animals, including swimming and flying animals also have motion controlling mechanisms that are self-regulating. For example, we have accurately calculated the motion of the cilia of a paramecium in water [16], and the planar motion of a bat in a room [10] using olivocerebellar control models, namely, the FitzHugh–Nagumo models. These wakes signify the similarity between temporal [10,17], or spatiotemporal [7] control and transitional turbulence mechanisms of hydrodynamics. The value of the recent advances in biohydrodynamics and biocontrol of propulsion [14,1821] in traditional engineering systems has not been explored.

The pristine initial condition in this work is set by very low advance ratios (J = 0.51) and Reynolds numbers (Rec = 3.75 × 103). The advance ratio (J) is defined as $U∞/(2fRavg)$, where $U∞$ is the tow velocity, $Ravg$ is the average blade radius $(rtip2+rroot2)/2$, $rtip$ and $rroot$ are the radii from the hub axis to the tip and root, respectively, f is the rate of rotation, and Reynolds number Rec is defined as $U∞c/ν$, where c is blade chord and $ν$ is kinematic viscosity of water. Disorganization and departure from the initial condition are measured by J and Rec, their ranges being varied by a factor of 10 (0.51 ≤ J ≤ 4.89; 3.75 × 103 ≤ Rec ≤ 3.75 × 104).

Our propulsor is operating at very low forward and rotational speeds where animals like penguins excel. In this work, the individual and isolated blades are called fins, because they are inspired by the pectoral fins of penguins [11,22]. In the context of ours and others novel lowspeed propulsors, we refer to them as blades [14,15,20,21].

In Fig. 2 of Ref. [19], the sound pressure levels versus frequency due to commercial ships and frigates are compared with those of fish and sea state 3. Figure 45 of Ref. [23] shows schematically the three unsteady force levels versus frequency in a typical propulsor attached behind a cylinder. They are: (1) turbulence ingestion from any upstream body which is a broad band source; (2) blade tonals which show up as discrete spikes in the spectrum due to wake cutting if there is any upstream stator; and (3) vibration of trailing edges which shows up as a narrow band at higher frequencies. In the present laboratory experimental work, there is a small hub, but there is no upstream body, and there is no upstream stator either. The speeds are low, and we do not expect any trailing edge vibration. Hence, the present propulsor experiences no external perturbations apart from mechanical vibration of carriage.

Below, we describe the propulsor and explain how the blade vibration is simulated. We also describe the propulsor blade surface and drag calculation. We give the measurements and models of time-averaged and temporal thrust. The relationship of the time-averaged thrust to the amplitude of blade pitch vibration angle $θo$ is examined and modeled. The relationship between the time-averaged and temporal thrust is explored. Temporal thrust is also modeled. We explore the deterministic nature of the unfolding disorganization in temporal thrust. We draw an analogy between the blade vibration in otherwise steady rotational blades and the unsteady propulsion mechanisms of the pectoral fins of penguins that undergo simultaneous orthogonal roll and pitch oscillations, where lift forces are enhanced due to the formation of leading edge vortices [11,22]. Finally, we show the common empirical relationship of the ratio of the boundary-layer control fence spacing to the fin/wing chord in unsteady whale flippers and in the steady wings of an aircraft.

## Description of the Propulsor

This is a large-diameter, low-speed propulsor (Fig. 1). The propulsor has six untwisted blades of small aspect ratio (span/chord = 3) [11,22]. The pitch angle is variable. The propulsor is hung from a six-component load cell on a carriage in a tow tank. Measurements of forces, power, current, blade position, rotational frequency, and tow speed are recorded. Current is normalized with a time-mean value. The propulsor blades were smooth, or roughened using a sand-grain roughness (60 grit), and a boundary layer fence (Fig. 1) was also used to make the flow chordwise in each otherwise smooth-walled blade. The fence is larger than the boundary layer thickness, but the roughness height is similar to the boundary layer thickness near the leading edges of the fins. The bare body drag versus speed was measured (D = 8.243 U1.88), where D(N) is the drag, and U (m/s) is the tow speed; the exponent is close to 2.0 (the theoretical value), and the bare body drag was removed from thrust measurements.

The blades have a National Advisory Committee for Aeronautics 0012-34 nominal section with a rounded leading edge and a rounded trailing edge (multiple layers of polyurethane coatings were applied to thicken and protect the trailing edge). The planform is an abstract penguin wing [11,22]. The blades are made of urethane of Shore hardness A20–A30, which is not as hard as metal. The chord, inner radius, outer radius, and planform area are 0.075 m, 0.0638 m, 0.22255 m, and 0.0104 m2/fin, respectively. The blades are hinged at a distance of c/3 from the leading edge [11,22] to allow efficient coupling between spinning (which has the same axis as rolling oscillation of pectoral fins of penguins) and pitch oscillations.

The thrust coefficient was calculated using the planform area of the blades. With a total of six blades, the total planform area is 0.0624 m2. A typical oscillating time trace of thrust is shown in Fig. 2. As the roll motor turns, the pitch motor applies a correction to maintain the amplitude of the pitch oscillation selected (see Ref. [14] on how roll and pitch motors can be made mechanically independent). This simulates blade vibration coupled to hub rotation.

## Time-Averaged Thrust Measurements and Modeling

The variations in the measurements of the thrust coefficient with the effective pitch angle (EPA) when the fin surfaces are smooth, sand-grain roughened, or contains boundary layer fences (Fig. 1) are shown in Figs. 3(a)3(c). They are overlaid in Fig. 3(d). The laminar-to-turbulence transition locations on the fins are influenced by sand-grain roughening and fencing, and they could affect the location of boundary layer separation if the boundary layer does separate. These variations are accounted for by deducting 5 deg from the effective pitch angle in the sand-grain roughened data and deducting 10 deg from the fence data. The angles 5 deg and 10 deg are determined by trial for best fit. The overlaid data provided in Fig. 3(e) show better collapse than in Fig. 3(d), meaning that all variables are being correctly scaled between different blade surface flows.

Figure 3(e) shows a branching of the trends into modes of a lower and a higher thrust coefficient past the EPA of about 10 deg. The coefficient of thrust is weakly and strongly dependent on changes in the effective pitch angle, in the lower and upper modes, respectively.

Note that the (lift-pitch) graph of a single steady (nonflapping) fin of similar dimensions as the propulsor blade also departs from linearity at 10 deg and maintains an approximately flat stall up to 50 deg [14]. Remarkably, the single fin measurements in the stall range of 25–30 deg show a slight dip in lift and that the lower mode measurements of thrust in the propulsor in Figs. 3 and 4 have even such a reduction in thrust. Hence, in the lower mode, the blades act like that they are in a steady flow and no leading edge vortex (LEV) forms enhancing lift force unlike that in the flapping case [14] and in the higher thrust mode.

We suggest that insight into the effects of blade vibration in propulsors may come from lift trends of isolated blades of such propulsors subjected to flapping motions.

The two modes of time-averaged thrust (Fig. 3(e)) have been modeled (Fig. 4) assuming that the blades vibrate in pitch and that, in the mode of higher thrust coefficient, the blade motion due to spin and pitch oscillation couples but that it does not couple in the other mode of lower thrust coefficient (Fig. 2). This is justified because the mode correlates accurately with the amplitude of the blade pitch vibration angle; the coupled mode is produced when θo is increased from 20 deg to 30 deg and 45 deg.

It is assumed that, when the pitch angle is large making the blade vibration large, an LEV forms on each blade as in flapping fins [10,11]. In the coupled mode modeling, we use the quasi-steady assumption found to be applicable in flapping fin propulsion where the spatiotemporal angles, rather than their rate of change, determine thrust and where stall is delayed up to large angles of attack [10,11,14,22]. The LEV is accounted for by the extrapolated steady lift curve beyond steady-state stall.

We make use of our bio-inspired single-fin lift and drag measurements [14] with and without flapping for modeling the coupled and uncoupled modes [11] (Fig. 4). It is assumed that when the spinning and pitch vibration effects in the blade boundary layer do not couple, the single-fin lift and drag behavior (versus angle of attack) in the (nonflapping) steady case can be used to resolve the thrust produced once the effect of hub rotation on the pitch angle is accounted for. Conversely, when the propulsor blade boundary layer does couple with the pitch oscillation, the hydrodynamic behavior can be modeled by the flapping-fin lift and drag characteristics. In both cases, viscous and pressure drags of the blades are included.

The thrust generated in the propulsor is given by $CTj=K(CLj1 cos α-CDj1 sin α)$, where CTj is the coefficient of thrust of the propulsor, whereas CLj1 and CDj1 are the lift and drag coefficients of the single fin, j = 0 or 1 (j = 0 when the blade boundary layer is not coupled to the pitch oscillation, and j = 1 when the fin boundary layer is coupled to the pitch oscillation), and α is the effective blade pitch angle. Here, K (<1) is a factor that accounts for the interference effects of the hub and the fins.

Define rotational velocity $vr=2πRavg(N/60)$, where N is the blade revolutions per minute, effective pitch angle = $αp-αo$, $αo$ is the geometric pitch angle, and blade pitch angle $αp=tan-1$ (vr/U). The total thrust coefficient $CT=2Fx/(ρUtotal2Aplanform)$, where Fx is the measured force (baseline hub-propulsor drag is removed) in the tow direction; and $Utotal=(U∞2+vr2)$.

Fifth-order polynomials were fitted to the single-fin lift and drag versus angle of attack (αo) measurements for both flapping and nonflapping cases. Their slopes are not much different, but the stall behavior is different. The modeled thrust is included in Fig. 4. Recall that, in the roughened and smooth-wall fenced blade cases, the pitch angle in the model is decreased by 5 deg and 10 deg, respectively, to account for the transition delay. A value of K = 0.3 is used.

Some justification for the choice of K = 0.3 can be found in our measurement of hydrodynamic efficiency in single fins and in a cylinder to which six fins are appended—three at each end [13]. The former has an efficiency of 0.60, while the latter has an efficiency of 0.42 at 0.4 m/s, which drops to about 0.15 at 0.1 m/s. In addition to blade–blade and blade-hub interference, there is a low-Reynolds-number effect; the efficiency is independent of speed above 0.46 m/s. The parameter K acts like hydrodynamic efficiency.

## Relationship Between Time-Averaged and Temporal Thrust

In the unfolding of temporal thrust disorganization, Fig. 5(a) shows the first departure from a limit cycle (LCO), namely, the appearance of phase jitter—an effect of the coupling of the main oscillator with an external oscillator. Below the values of J and Rec of Fig. 5(a), the amplitude of the external oscillator was likely insufficient for phase jitter of the coupling to occur, so the map was a pure limit cycle. As J and Rec are increased (Fig. 5(c)), the disorganization takes hold. The thrust attractors of patterns in Figs. 5(a) and 5(c) are marked in Figs. 5(b) and 5(d).

Figures 68 show three collages of the thrust-current maps for amplitude of pitch angle vibration of 30 deg for smooth (Fig. 6), roughened (Fig. 7), and smooth-walled fenced blades (Fig. 8). Lines of constant Rec connect maps with increasing J. For example, at the bottom of Fig. 6, for Rec = 7470.1, as J increases to 1.08, 1.22, 1.41, 1.67, 2.04, and 2.62, the limit cycles shift bodily along the thrust axis, increasing the disorganization. The shift occurs when the current reaches a cycle-minimum value. The trend is similar at other values of Rec and surface condition (and at θo = 20 deg, 45 deg).

###### Synthesis of How Disorganization in Thrust Unfolds.

Measuring disorganization by the standard deviation of $CT$ ($CT̃$), Fig. 9 shows the trend for all boundary conditions (BC) and all Rec. In Fig. 9, some of the trends have large amplification suggesting the presence of an underlying instability after laminar to turbulence transition has occurred. Figure 10 gives the individual trends. In Figs. 10(a)10(c), the distributions of $(J-CT̃$) are separated based on wall BCs—smooth (a), rough (b), and fenced (c). The unfolding of the disorganization occurs earlier at lower Reynolds numbers. In the present data range (J < 4), the $CT̃$ values are not high at the higher Rec. This means that the unfolding process being observed is viscosity dominated, and the experimental data cover a regime where the viscous domination is giving way to inertia effects.

When the blades are fenced, the tipward flow is restricted, and the flow becomes chordwise (see schematic in Fig. 12(d)). The fencing attenuates disorganization (Fig. 10(c)). The rough wall BC both allows tipward flow and energizes the near-wall flow, the former being similar to the smooth BC case. Control of the boundary layer streamline curvature controls disorganization. From a dynamical consideration (see below), smooth and rough walls act as natural flows, where there is no control of lateral diffusion of vorticity. In contrast, the fence acts as a device for controlling the lateral diffusion of vorticity. Fencing is an external impetus of body scale that is countering the natural body force of curvature.

The trends in Figs. 10(a)10(c) are synthesized in Fig. 10(d). Before the map becomes a blurred LCO (Fig. 5(a)), at lower Rec and J, a closed LCO without jitter must have occurred. For a blurred LCO (see below), $CT̃$ < 0.03, approximately. The jump in $CT̃$ that occurs in some instances is followed by a reduction in $CT̃$. The jump in $CT̃$ is an indication of the onset of a change in the pattern. The later reduction in $CT̃$ is an indication of bistable behavior (hysteresis).

The emergence of blurred LCO in phase maps, that is departure from LCO, is the first indication that the oscillation is no longer locked-in, and the condition for phase equal to 0 deg does not return any more. In other words, there is excess energy in the system. If destabilizing forces such as inertial forces increase further, such departures would also increase. Eventually, the nature of coupling of oscillations in orthogonal axes would also change and new input–output patterns may emerge.

The blurring is also an indication of mismatch of all higher derivatives of the state variable at the end of the cycle of oscillation. New maps of acoustic radiation can then be expected. In flapping fin oscillations, for example, that is common in large swimming animals, initially, the blurring is systematic, because the system indicates lock-in after eight cycles of oscillation [11]. This discussion suggests a path how disorganization in the propulsor may unfold.

Recall that in flapping fin propulsion in the range of $3558≤Rec≤70,895$, bistable thrust fluctuations appear over $8902≤Rec≤13,154$ and $31,345≤Rec≤35,686$ [11,13]. For theoretical modeling, see Ref. [10]. They indicate that the vortices shed from the pressure, and the suction sides of the vibrating blades have unequal absolute values of circulation. The order of inequality is switching in these two ranges of bistability between the pressure and the suction sides of the blades. Therefore, the root-mean-square thrust fluctuation graphs are measures of vortex–vortex interactions in the blade wakes. In other words, the thrust patterns in Figs. 68 are the histories of interactions of the wake vortices and the process appears to be organized. The dynamical modeling below suggests that this organization is temporally deterministic.

###### Modeling of Temporal Thrust.

Two models are given in Fig. 11. The first is systemic (Eq. (1); Figs. 11(a) and 11(b)), and the second (Eq. (2); Fig. 11(f)) is fluid-structural. Both models include effects that are internal and external to the main oscillator.

An oscillator of state (x(t)) under external periodic forcing is modeled as Display Formula

(1)$d2xdt2-1-xndxdt+Bsinkt+x=τo-Bkcos(kt)$

where the diffusion within the oscillator is a power relationship that is also a periodic function of the state and is given by $1-xndx/dt+Bsinkt$. The external forcing that is orthogonal to the main (roll) oscillator is also oscillatory and has some offset; it is given by $τo-Bkcos(kt)$, where t is the time. In Figs. 11(a) and 11(b), which shows the results of Eq. (1), the parameters are: B = 0.15, $τo$ = −0.75, k = 1.8, and n = 2 in (a) and 8 in (b).

In Fig. 11, two sample measurements of thrust trajectories (Figs. 11(c) and 11(d)) for smooth fins at low values of J and Rec are compared with the systemic model (Figs. 11(a) and 11(b): Eq. (1)). The model in Eq. (1) has two effects—a nonlinear diffusion internal to the thrust oscillator state x(t) (left-hand terms) and an external forcing (right-hand terms, where the state variable (x(t)) does not appear). The oscillator diffusion exponent n becomes more dominant with increasing J and Rec. The cosinusoidal external oscillator, which does not change with n, is orthogonal to the main sinusoidal oscillator, which indicates its weaker strength but decisive role in control. The disorganization unfolds with increasing values of the power law exponent n of the diffusion (n from 2 in Fig. 11(a) to 8 in Fig. 11(b)). Hence, n is a measure of a change in patterns due to amplification of some yet unknown instability triggered by increasing values of J and Rec (Figs. 68).

The trajectories for higher values of J and Rec (Fig. 11(f)) were modeled using the damped oscillatory equation for a structure as follows: Display Formula

(2)$Mϕ¨+cϕ˙+kϕ-Ωnomt=τϕ$

where M is the moment of inertia, c is the damping, k is the stiffness, $Ωnom$ is the attractor frequency (rad), $τϕ$ is the roll torque, $ϕ˙(t)$ is the roll rate, and $τϕ(ϕ˙,U)$ comes from quasi-steady fin modeling [10,11], where the tow speed of 0.2 units is assumed to have 10% noise. In (f) which shows the results of Eq. (2), M = 0.218, c = 0.1, $Ωnom=1/4π$, $ϕ˙$ is abscissa, and ordinate is thrust.

Figure 11(f) shows the fluid-structure model (Eq. (2)) compared with thrust measurements (Figs. 11(e) and 11(g)) at higher values of J and Rec for fenced and roughened blades. The system is modeled as a damped oscillator. The cyclic bodily shift in thrust is reproduced. In the unfolding of disorganization, this model indicates the increasing roles: (1) the coupling of the fin boundary-layer with the oscillatory wake, (2) the transition location of the fin boundary layer, and (3) the oscillations in forward motion with increasing J and Rec. Recall that similarly, in bio-inspired flapping fin propulsion the shed wake vortices alter the incoming angle of attack, whereby the wake and the fin boundary layers become coupled [1012].

For smooth blades in (c), J = 0.51, Rec = 3735, effective angle of attack = 36 deg. For smooth blades in (d), J = 3.52, Rec = 18,675, effective angle of attack = 25 deg. For fenced blades (e), J = 2.04, Rec = 7470.1, and effective angle of attack = 30 deg. For roughened blades (g), J = 1.67, Rec = 7470.1, and effective angle of attack = 32 deg.

Figures 68 show three effects on the unfolding of disorganization. In these three figures, the LCO patterns are “coat-hanger-like” when the blade is smooth but are less so when the blade is fenced or roughened. If two propulsor cases are compared, it is seen that their blade transition location also affects the unfolding of the disorganization, in addition to J and Rec.

## Discussion

For a three-dimensional state variable $w(t)$, the generalized spatiotemporal equation of motion is $ρDw/Dt=F+P$, where F is the body force, and P is the surface force. In vector notation, $P=-gradp+μ∇2w$, where μ is absolute viscosity of the fluid, $∇2$ is the Laplace operator $∇2=∂2/∂x2+∂2/∂y2+∂2/∂z2$, and $μ∇2w$ is the viscous term.

The existence of many niche actuators in low Reynolds numbers in aquatic swimming animals suggests a relationship with the ratio of the destabilizing and stabilizing forces, and that the relationship is like a staircase with the stairs being hysteretic [10,13]. Here, we are looking for an analogous behavior in the case of conventional propulsors.

In this work, 3735 < Rec < 37,500 covers the transition from the domination of viscous forces $μ∇2w$ over the inertial forces ($ρDw/Dt$) to the domination of ($ρDw/Dt$) over ($μ∇2w$); here, ρ is the fluid density. For Rec ≥ 14,940, the disorganization ($CT̃$) does not show amplification as at lower values of Rec. As inertial forces increase, structural vibration modes may act as external perturbations to thrust. Such effects are modeled as external oscillators in Eq. (1). Blade vibration is such an external oscillator.

The amplification of $CT̃$ in the baseline case is large over 3735 ≤ Rec ≤ 11,205, but, if F = 0 (fenced case when streamlines are not curved), the amplification of $CT̃$ is attenuated (Fig. 10(c)). The useful result is that, at least at very low Reynolds numbers, the growth of disorganization in thrust can be delayed by the reduction in streamline curvature.

The leading edges of whale flippers have several tubercles that are similar to a boundary layer fence. The dermal denticles of shark skins have numerous tiny, longitudinal ridges (also known as riblets) deeply buried in their turbulent boundary layers; these ridges are akin to wall-layer scale fences. When the ridges [19] are modeled [7,18] as devices for the control of the lateral diffusion of vorticity fluctuation (called μ-control), the wanderings of the wall-layer streaks are remarkably stopped (Fig. B-4 in Ref. [18]).

Similarly, the propulsor fence, or the laterally arrayed protrusions of the whale flipper leading edge that produces streamwise vortices effectively acting as fences, can act as lateral diffusion control devices (boundary-layer scale μ-control in Figs. 1 and 12). Such control is explicitly treated in the Stuart–Landau equation approach. The critical aspect in the μ-control is the ratio (Lc) of the wavelength of the instability due to spanwise streamline curvature (fence spacing λ) to the mean flow wavelength (or the chord c). See schematic in Fig. 12(d). In the propulsor, the hub acts as a fence. The fence height is controlled by the mechanism of lateral instability. Note that quasi-steady assumptions apply to unsteady flapping propulsion [22]. The ratio λ/c is 2 in the propulsor (arbitrarily selected) and approximately 0.5–0.75 in the whale pictures in Fig. 12. It is 0.5 and 1.0 in the MiG-17 wings. All of these examples of wings have large aspect ratios (AR = span/chord ≥ 3). Hence, the control of lateral diffusion of vorticity fluctuation of body scale (curvature of streamlines) is the mechanism of the reduction in disorganization due to fencing. In an engineering context—steady or unsteady, if the aspect ratio is large, then, to accommodate several fences, the ratio Lc should be within 0.50–1.0.

To reduce the exposure of aquatic animals to ship radiated noise, vibration is sometimes isolated (as opposed to being introduced externally and minimally) to prevent its transmission from the propulsor to the hull (see Fig. 2 in Ref. [19]). But this remedy is not designed to organize the temporal propulsor thrust. This work, on the other hand, shows the disorganization in temporal thrust in a practical context as a deterministic process, albeit at low transitional Reynolds numbers, thus opening the door to rational control when transition cannot be delayed.

## Conclusions

We have built a lowspeed underwater propulsor where the blades undergo synchronous pitching oscillations simulating blade vibration, while the hub is spinning. Three types of blades are used: smooth, or roughened, or smooth-walled and fenced. The advance ratios and Reynolds numbers are transitional: 0.51 ≤ J ≤ 4.89 and 3.75 × 103 ≤ Rec ≤ 3.75 × 104.

We find that large amplitude pitch vibrations act similar to flapping fin propulsion of penguins. We have examined the emergence of flow disorganization by examining the departure from LCO in a novel attempt at flow control of the oscillatory flow in this propulsor in order to gain some understanding of the general principles of unsteady flows.

Animals having a diverse array of propulsive actuators swim in the above range of Reynolds numbers. Theoretical oscillatory wake modeling and olivocerebellar control have similar self-regulating mechanisms in swimming animals [1014]. For cruising, animals maintain LCO for optimal swimming. They maintain the Strouhal number of oscillation in a narrow range indicating resonant tuning [24,25]. Conversely, departure from LCO is costly and animals in their niche of low Reynolds numbers have evolved strategies for control of disorganization.

In this work, we have taken a similar approach to arrive at an internally consistent new method of analysis of temporal thrust measurements in a propulsor. Boundary layer fences are found to reduce the departure from LCO. The leading edge tubercles in whale fins and at least some aircraft wings of similar aspect ratios have fence spacing to chord ratios in the range 0.5–1.0. The suggestion is that in both nature and engineering, similar principles of control of disorganization at least at low Reynolds numbers are found to be beneficial when unsteadiness is present.

## Acknowledgements

W. H. Nedderman is thanked for the industrial design and fabrication of the propulsor. Daniel Thivierge built the motor controller. Albert Fredette assisted with tow tank data acquisition. Aren Hellum assisted with the self-regulation modeling and figure preparation.

## Funding Data

• The U.S. Office of Naval Research, Biology-Inspired Autonomous Systems Program (ONR 341) to PRB (Grant No. N0001414AF00002).

## Nomenclature

• c =

damping

• c =

fin chord

• CDj1 =

drag coefficient of single fin

• CLj1 =

lift coefficient of single fin

• CT =

total thrust coefficient

• CTj =

coefficient of thrust of propulsor

• $CT̃$ =

standard deviation of CT

• D(N) =

drag

• deg =

degree

• f =

frequency of oscillation

• F =

body force

• Fx =

thrust, measured force in the tow direction

• J =

• k =

stiffness

• K =

factor accounting for interference effects of hub and blades

• Lc =

ratio of fence spacing λ to chord c

• M =

moment of inertia

• n =

oscillator diffusion exponent

• N =

fin revolutions per minute

• P =

surface force

• p′ =

fluctuating pressure

• rroot =

• rtip =

• Ravg =

• Rec =

Reynolds number based on chord length

• t =

time

• U =

tow speed

• $U∞$ =

tow velocity

• α =

effective fin pitch angle

• αo =

angle of attack

• θo =

amplitude of blade pitch vibration angle

• λ =

fence spacing

• μ =

absolute viscosity

• νr =

rotational velocity

• ρ =

fluid density

• τϕ =

roll torque

• υ =

kinematic viscosity of water

• $ϕ˙(t)$ =

roll rate

• Ωnom =

## References

Ross, D. , 2005, “ Ship Sources of Ambient Noise,” IEEE J. Oceanic Eng., 30(2), pp. 257–261.
Bandyopadhyay, P. R. , 1980, “ Large Structure With a Characteristic Upstream Interface in Turbulent Boundary Layers,” Phys. Fluids, 23(11), pp. 2326–2327.
Head, M. R. , and Bandyopadhyay, P. R. , 1981, “ New Aspects of Turbulent Boundary-Layer Structure,” J. Fluid Mech., 107, pp. 297–338.
Perry, A. E. , and Chong, M. S. , 1982, “ On the Mechanism of Wall Turbulence,” J. Fluid Mech., 119, pp. 173–217.
Adrian, R. J. , Meinhart, C. D. , and Tomkins, C. D. , 2000, “ Vortex Organization in the Outer Region of the Turbulent Boundary Layer,” J. Fluid Mech., 422, pp. 1–54.
Wu, X. , and Moin, P. , 2009, “ Direct Numerical Simulation of Turbulence in a Nominally Zero-Pressure-Gradient Flat Plate Boundary Layer,” J. Fluid Mech., 630, pp. 5–41.
Bandyopadhyay, P. R. , and Hellum, A. M. , 2014, “ Modeling How Shark and Dolphin Skin Patterns Control Transitional Wall-Turbulence Vorticity Patterns Using Spatiotemporal Phase Reset Mechanisms,” Sci. Rep., 4, p. 6650. [PubMed]
Bandyopadhyay, P. R. , 1996, “ REVIEW—Mean Flow in Turbulent Boundary Layers Disturbed to Alter Skin Friction,” ASME J. Fluids Eng., 108(2), pp. 127–140.
Golestani, A. , and Gras, R. , 2014, “ Can We Predict the Unpredictable?,” Sci. Rep., 4, p. 6834. [PubMed]
Bandyopadhyay, P. R. , Leinhos, H. A. , and Hellum, A. M. , 2013, “ Handedness Helps Homing in Swimming and Flying Animals,” Sci. Rep., 3, p. 1128. [PubMed]
Bandyopadhyay, P. R. , Beal, D. N. , Hrubes, J. D. , and Mangalam, A. , 2012, “ Relationship of Roll and Pitch Oscillations in a Fin Flapping at Transitional to High Reynolds Numbers,” J. Fluid Mech., 702, pp. 298–331.
Von Ellenrieder, K. D. , Parker, K. , and Soria, J. , 2008, “ Fluid Mechanics of Flapping Wings,” Exp. Therm. Fluid Sci., 32(8), pp. 1578–1589.
Bandyopadhyay, P. R. , and Leinhos, H. A. , 2013, “ Propulsion Efficiency of Bodies Appended With Multiple Flapping Fins: When More Is Less,” Phys. Fluids, 25(4), p. 041902.
Bandyopadhyay, P. R. , 2015, “ A Novel Large Slosh-or-Spin Low-Speed Underwater Propulsor Bridges the Unsteady and Steady Propulsion Mechanisms of Nature and Engineering,” IEEE J. Oceanic Eng., 41(4), pp. 868–881.
Bandyopadhyay, P. R. , 2016, “ Highly Maneuverable Biorobotic Underwater Vehicles—Part B: 11.1 Autonomous Ocean Vehicles, Subsystems and Control,” Springer Handbook of Ocean Engineering, M. R. Dhanak and X. I. Nikolas , eds., Springer, New York, pp. 281–300.
Bandyopadhyay, P. R. , and Hansen, J. C. , 2013, “ Breakup and Then Makeup: A Predictive Model of How Cilia Self-Regulate Hardness for Posture Control,” Sci. Rep., 3, p. 1956. [PubMed]
Huerre, P. , and Monkewitz, P. A. , 1990, “ Local and Global Instabilities in Spatially Developing Flows,” Annu. Rev. Fluid Mech., 37(22), pp. 1145–1151.
Bandyopadhyay, P. R. , and Hellum, A. M. , 2014, “Llinas' Phase Reset Mechanism Delays the Onset of Chaos in Shark and Dolphin Wall Turbulence,” Naval Undersea Warfare Center, NUWC-NPT, Newport, RI, Technical Report No. 12,154.
Bandyopadhyay, P. R. , 2009, “ Swimming and Flying in Nature—The Route Toward Applications: The Freeman Scholar Lecture,” ASME J. Fluids Eng., 131(3), p. 031801.
Usab, W. J. , Hardin, J. , and Bilanin, A. J. , 2004, “ Bioinspired Delayed Stall Propulsor,” IEEE J. Oceanic Eng., 29(3), pp. 756–765.
Bandyopadhyay, P. R. , Singh, S. N. , Thivierge, D. P. , Annaswamy, A. M. , Leinhos, H. A. , Fredette, A. R. , and Beal, D. N. , 2008, “ Synchronization of Animal-Inspired Multiple Fins in an Underwater Vehicle Using Olivo-Cerebellar Dynamics,” IEEE J. Oceanic Eng., 33(4), pp. 563–578.
Bandyopadhyay, P. R. , Beal, D. N. , and Menozzi, A. , 2008, “ Biorobotic Insights Into How Animals Swim,” J. Exp. Biol., 211, pp. 206–214. [PubMed]
Bandyopadhyay, P. R. , 2005, “ Trends in Biorobotic Autonomous Undersea Vehicles,” IEEE J. Oceanic Eng., 30(1), pp. 109–139.
Triantafyllou, M. S. , and Triantafyllou, G. S. , 1995, “ An Efficient Swimming Machine,” Sci. Am., 272, pp. 64–70.
Taylor, G. K. , Nudds, R. L. , and Thomas, A. L. R. , 2003, “ Flying and Swimming Animals Cruise at a Strouhal Number Tuned for High Power Efficiency,” Nature, 425, pp. 707–710. [PubMed]
View article in PDF format.

## References

Ross, D. , 2005, “ Ship Sources of Ambient Noise,” IEEE J. Oceanic Eng., 30(2), pp. 257–261.
Bandyopadhyay, P. R. , 1980, “ Large Structure With a Characteristic Upstream Interface in Turbulent Boundary Layers,” Phys. Fluids, 23(11), pp. 2326–2327.
Head, M. R. , and Bandyopadhyay, P. R. , 1981, “ New Aspects of Turbulent Boundary-Layer Structure,” J. Fluid Mech., 107, pp. 297–338.
Perry, A. E. , and Chong, M. S. , 1982, “ On the Mechanism of Wall Turbulence,” J. Fluid Mech., 119, pp. 173–217.
Adrian, R. J. , Meinhart, C. D. , and Tomkins, C. D. , 2000, “ Vortex Organization in the Outer Region of the Turbulent Boundary Layer,” J. Fluid Mech., 422, pp. 1–54.
Wu, X. , and Moin, P. , 2009, “ Direct Numerical Simulation of Turbulence in a Nominally Zero-Pressure-Gradient Flat Plate Boundary Layer,” J. Fluid Mech., 630, pp. 5–41.
Bandyopadhyay, P. R. , and Hellum, A. M. , 2014, “ Modeling How Shark and Dolphin Skin Patterns Control Transitional Wall-Turbulence Vorticity Patterns Using Spatiotemporal Phase Reset Mechanisms,” Sci. Rep., 4, p. 6650. [PubMed]
Bandyopadhyay, P. R. , 1996, “ REVIEW—Mean Flow in Turbulent Boundary Layers Disturbed to Alter Skin Friction,” ASME J. Fluids Eng., 108(2), pp. 127–140.
Golestani, A. , and Gras, R. , 2014, “ Can We Predict the Unpredictable?,” Sci. Rep., 4, p. 6834. [PubMed]
Bandyopadhyay, P. R. , Leinhos, H. A. , and Hellum, A. M. , 2013, “ Handedness Helps Homing in Swimming and Flying Animals,” Sci. Rep., 3, p. 1128. [PubMed]
Bandyopadhyay, P. R. , Beal, D. N. , Hrubes, J. D. , and Mangalam, A. , 2012, “ Relationship of Roll and Pitch Oscillations in a Fin Flapping at Transitional to High Reynolds Numbers,” J. Fluid Mech., 702, pp. 298–331.
Von Ellenrieder, K. D. , Parker, K. , and Soria, J. , 2008, “ Fluid Mechanics of Flapping Wings,” Exp. Therm. Fluid Sci., 32(8), pp. 1578–1589.
Bandyopadhyay, P. R. , and Leinhos, H. A. , 2013, “ Propulsion Efficiency of Bodies Appended With Multiple Flapping Fins: When More Is Less,” Phys. Fluids, 25(4), p. 041902.
Bandyopadhyay, P. R. , 2015, “ A Novel Large Slosh-or-Spin Low-Speed Underwater Propulsor Bridges the Unsteady and Steady Propulsion Mechanisms of Nature and Engineering,” IEEE J. Oceanic Eng., 41(4), pp. 868–881.
Bandyopadhyay, P. R. , 2016, “ Highly Maneuverable Biorobotic Underwater Vehicles—Part B: 11.1 Autonomous Ocean Vehicles, Subsystems and Control,” Springer Handbook of Ocean Engineering, M. R. Dhanak and X. I. Nikolas , eds., Springer, New York, pp. 281–300.
Bandyopadhyay, P. R. , and Hansen, J. C. , 2013, “ Breakup and Then Makeup: A Predictive Model of How Cilia Self-Regulate Hardness for Posture Control,” Sci. Rep., 3, p. 1956. [PubMed]
Huerre, P. , and Monkewitz, P. A. , 1990, “ Local and Global Instabilities in Spatially Developing Flows,” Annu. Rev. Fluid Mech., 37(22), pp. 1145–1151.
Bandyopadhyay, P. R. , and Hellum, A. M. , 2014, “Llinas' Phase Reset Mechanism Delays the Onset of Chaos in Shark and Dolphin Wall Turbulence,” Naval Undersea Warfare Center, NUWC-NPT, Newport, RI, Technical Report No. 12,154.
Bandyopadhyay, P. R. , 2009, “ Swimming and Flying in Nature—The Route Toward Applications: The Freeman Scholar Lecture,” ASME J. Fluids Eng., 131(3), p. 031801.
Usab, W. J. , Hardin, J. , and Bilanin, A. J. , 2004, “ Bioinspired Delayed Stall Propulsor,” IEEE J. Oceanic Eng., 29(3), pp. 756–765.
Bandyopadhyay, P. R. , Singh, S. N. , Thivierge, D. P. , Annaswamy, A. M. , Leinhos, H. A. , Fredette, A. R. , and Beal, D. N. , 2008, “ Synchronization of Animal-Inspired Multiple Fins in an Underwater Vehicle Using Olivo-Cerebellar Dynamics,” IEEE J. Oceanic Eng., 33(4), pp. 563–578.
Bandyopadhyay, P. R. , Beal, D. N. , and Menozzi, A. , 2008, “ Biorobotic Insights Into How Animals Swim,” J. Exp. Biol., 211, pp. 206–214. [PubMed]
Bandyopadhyay, P. R. , 2005, “ Trends in Biorobotic Autonomous Undersea Vehicles,” IEEE J. Oceanic Eng., 30(1), pp. 109–139.
Triantafyllou, M. S. , and Triantafyllou, G. S. , 1995, “ An Efficient Swimming Machine,” Sci. Am., 272, pp. 64–70.
Taylor, G. K. , Nudds, R. L. , and Thomas, A. L. R. , 2003, “ Flying and Swimming Animals Cruise at a Strouhal Number Tuned for High Power Efficiency,” Nature, 425, pp. 707–710. [PubMed]

## Figures

Fig. 1

Six-bladed propulsor with boundary layer fence in each blade, where c is the chord, and λ is the spacing from the hub

Fig. 2

Example of thrust oscillations in the spinning propulsor. The solid line denotes the long-time-averaged thrust of 0.62 N (the bare body drag has been removed). Tow speed = 0.2 m/s; rotational speed = 15 rpm.

Fig. 3

Measurements of transitional thrust. Data from (a)–(c) are overlaid in (d) and (e). In the rough and fence data in (e), the effective pitch angles are adjusted: angles of 5 deg and 10 deg determined by trial are deducted from the rough blade data and smooth blade fence data, respectively. Filled triangles represent the condition of thrust equal to drag within ±10%. Symbols in (d) and (e) represent smooth (O), rough (X), fence (◻). In this figure and in Figs. 4 and 5, color bar: Rec, blue: 0; red: 30,000; and size of measurement symbol (numerals in Fig. 4): J.

Fig. 4

Comparison of models of time-averaged transitional-thrust (lines) with measurements from Fig. 3(e) indicated by numerals. The blue line (monotonic up line) represents the modeled coupled, unsteady condition; the red line (bent line) represents the modeled uncoupled, steady condition. Instead of symbols as in Fig. 3(e), the plotted measurements in Fig. 4 are marked by the numbers 20, 30, and 45 that represent the amplitudes of the blade pitch vibration angles in deg (θo). (Color bar: Rec, blue: 0; red: 29,880.)

Fig. 5

Paired plots showing the effects of J and Rec on the temporal thrust trajectories (a) and (c) of measured time-averaged transitional-thrust (b) and (d). The “clock hands” in (b) and (d) show the perspective of the time-averaged values of the temporal thrust in (a) and (c), respectively. (a) J = 0.51, Rec = 3735, effective angle of attack = 36 deg; (c): J = 2.04, Rec = 7470, effective angle of attack = 37 deg. Color bars (chord Reynolds number): blue (0), red (30,000).

Fig. 6

Collage of transitional-thrust trajectories in smooth blades showing the effects of advance ratio and Reynolds number. θo: 30 deg. The colored lines connecting the plots show that J decreases (left to right and down) for fixed values of Rec. The pair of small orthogonal lines at the left margin of the top left of plot indicate 0.05 of unit coefficient of thrust (ordinate) and normalized current (abscissa).

Fig. 7

Collage of transitional-thrust trajectories in sand-roughened blades showing the effects of advance ratio and Reynolds number. θo: 30 deg. The colored lines connecting the plots show J decreasing, left to right and down, for fixed values of Rec. The pair of orthogonal lines in top left plot indicate 0.05 of unit coefficient of thrust (ordinate) and normalized current (abscissa).

Fig. 8

Collage of transitional-thrust trajectories in fenced smooth-walled blades showing the effects of advance ratio and Reynolds number. θo: 30 deg. The colored lines connecting the plots show J decreasing, left to right and down, for fixed values of Rec. The pair of orthogonal lines in top left plot indicate 0.05 of unit coefficient of thrust (ordinate) and normalized current (abscissa).

Fig. 9

Root-mean-square of thrust trace (CT̃) versus J for all data (all Rec). Re is Rec. The CT̃ values do not amplify much in the fenced case (dotted line, black fill) (see Fig. 10(c)).

Fig. 10

Root-mean-squares of disorganization in thrust patterns. (a)–(c): Rec = 3735 (magenta), 7470 (blue), 11,205 (green), 14,940 (red), 22,410 (cyan); θo = 20 deg (O), 30 deg (▵), and 45 deg (◻) (as in Fig. 9). (d) Schematic shifts in LCO patterns in (a)–(c); the pattern index (CT̃) is scaled to the boundary condition; solid line is data; broken line is hypothesized; the jumps in patterns are hysteretic when the pattern is experiencing a bimodal instability [10,11,13]; the colored boxes represent the zones of J where the maps have characteristic patterns; and LCO is limit cycle oscillation.

Fig. 11

Dynamical modeling of temporal transitional-thrust compared with sample measurements. Systemic models (a) and (b) are compared with sample measurements of thrust trajectories (c) and (d) for smooth blades at low values of J and Rec. Fluid-structural model (f) is compared with thrust measurements (e) and (g) at higher values of J and Rec for fenced and roughened blades, respectively.

Fig. 12

The lateral diffusion control ratio Lc = λ/c (ratio controlling lateral diffusion of vorticity [7,18]). (a) Humpbacked whale flipper (photograph courtesy of Giles Breton). (b) Humpbacked whale in Atlantic; aspect ratio = 15 (photograph courtesy of Giles Breton). (c): MiG-17; aspect ratio = 3.7 (photograph courtesy of Balon Greyjoy). In both steady and unsteady wings, 0.50 ≤ λ/c ≤ 1.00. (d) Schematic of a streamline on a large aspect ratio wing and hydrofoil showing amplified boundary-layer instability of wavelengths λ and c of spanwise and chordwise flows in the unfenced smooth-wall case.

## Discussions

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections