{'fraction_non_alphanumeric': 0.03191653616251627, 'fraction_numerical': 0.015042900941835587, 'mean_word_length': 4.658055659553226, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 0, 'https://': 1, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Contrasting thrust generation mechanics and energetics of ﬂapping foil locomotory states characterized by a uniﬁed St - Re scaling

Self-propelled ﬂapping foils with distinct locomotion-enabling kinematic restraints exhibit a remarkably similar Strouhal number ( St )-Reynolds number ( Re ) dependence. This similarity has been hypothesized to pervade diverse forms of oscillatory self-propulsion and undulatory biolocomotion; however, its genesis and implications on the energetic cost of locomotion remain elusive. Here, using high-resolution simulations of translationally free and restrained foils that self-propel as they are pitched, we demonstrate that a generality in the St - Re relationship can emerge despite signiﬁcant disparities in thrust generation mechanics and locomotory performance. Speciﬁcally, owing to a recoil reaction induced passive heave, the ﬂuid’s inertial response to the prescribed rotational pitch, the principal source of thrust in unidirectionally free and towed conﬁgurations, ceases to produce thrust in a bidirectionally free conﬁguration. Rather, the thrust generated from the leading edge suction mechanics self-propels a bidirectionally free pitching foil. Owing to the foregoing distinction in the thrust generation mechanics, the St - Re relationships for the bidirectionally and unidirectionally free/towed foils are dissimilar and pitching amplitude dependent, but speciﬁcally for large reduced frequencies, converge to a previously reported uniﬁed power law. Importantly, to propel at a given mean forward speed, the bidirectionally free foil must counteract the out-of-phase passive heave through a more intense rotational pitch, resulting in an appreciably higher power consumption over the range 10 ≤ Re ≤ 10 3 . We highlight the critical role of thrust in introducing an offset in the St - Re relation, and through its ampliﬁcation, being ultimately responsible for the considerable disparity in the locomotory performance of differentially constrained foils.


Introduction
Rigid foils that undergo a periodic rotational pitch and/or transverse heave exhibit an exceptional tendency to self-propel unidirectionally (Alben & Shelley 2005;Spagnolie et al. 2010;Zhang, Liu & Lu 2010;Alben et al. 2012;Zhu, He & Zhang 2014;Deng & Caulfield 2016;Verma et al. 2017;Das, Shukla & Govardhan 2019;Lin, Wu & Zhang 2021). This tendency emerges as the thrust produced from a sufficiently intense prescribed harmonic excitation surmounts the hydrodynamic resistance to the foil's motion. At a mean self-propelling state, a precise balance between the mean thrust and average drag is established. This balance naturally leads to a relation between the mean self-propelling speed and the imposed kinematics and is commonly expressed in the form of a non-dimensional Strouhal number-Reynolds number dependence (Triantafyllou et al. 2005;Gazzola, Argentina & Mahadevan 2014;Smits 2019). For a typical self-propelled flapping foil state, the Strouhal number St = fA/ū p and the Reynolds number Re = u p c/ν represent the normalized forcing frequency and propulsive speed, respectively.
Hereū p , f , A, c and ν denote the mean forward speed (an overbar indicates time average), excitation frequency, peak-to-peak trailing edge excursion, foil's chord length and kinematic viscosity of the fluid medium, respectively.
The prescribed kinematics have a direct bearing on the trailing edge excursion A and, therefore, the Strouhal number. It is then reasonable to hypothesize that the St-Re relationship will exhibit significant variations across the diverse range of self-propelling states that flapping foils with distinct imposed kinematics display. Counterintuitively however, the St-Re relationship has been shown to possess a remarkable universal character. Specifically, early investigations have shown that the range of Strouhal numbers over which swimming fishes and flying birds cruise is exceptionally narrow (0.2 ≤ St ≤ 0.35, see Taylor, Nudds & Thomas (2003), Eloy (2012) and review articles by Triantafyllou et al. (2005) and Smits (2019). Similar observations have been made in the recent work of Gazzola et al. (2014) wherein a balance between the thrust, hypothesized to arise from the fluid's inertial reaction to lateral body motion, and the skin friction drag from a turbulent boundary layer has been shown to desirably yield a constant Strouhal number, in agreement with the observations spanning a wide range of natural swimmers. The narrow range 0.2 ≤ St ≤ 0.35 coincides with the range of Strouhal numbers over which high-Re oscillating foils produce thrust at a peak propulsive efficiency (Triantafyllou, Triantafyllou & Grosenbaugh 1993). The foregoing similarity in the range of St for undulatory locomotion and high efficiency thrust generation from oscillating foils is of substantial significance with wide-ranging implications. Notably, the similarity in the range of St justifies analysis of undulatory locomotion from a global unified viewpoint and lends support to the use of relatively simple foil-like flapping states as a model system to systematically investigate much more complex self-propelling states associated with undulating natural swimmers.
The St-Re relationship is fundamentally dynamical and arises naturally from a balance between a time-averaged thrust generated from body undulations and a viscous resistance to the forward motion. Its simple dynamical origin has an undeniable appeal. Nonetheless, a St-Re relation based dynamical viewpoint remains detached from the energetics of self-propulsion. The energetic cost of locomotion is often the single most important limiting constraint that determines viability of artificial locomotors and survivability of natural swimmers. A meaningful quantification of the energetic cost of locomotion is therefore central to any unified framework of self-propulsion hydrodynamics. Attempts have been made to link the narrow range of Strouhal numbers to an optimality in the energetics of the wake vorticity distribution across self-propelled organisms and thrust Contrasting dynamics and energetics of self-propelled foils producing flapping foils (Triantafyllou et al. 1993). However, the momentumless wakes associated with self-propelling bodies are distinct from the excess-momentum wakes associated with thrust-generating tethered foils. Moreover, recent works suggest that the optimal wake energetics are more likely a consequence rather than a cause for efficient propulsion (Arbie, Ehrenstein & Eloy 2016). The implication of a generality in the St-Re relation on the energetic cost of self-propulsion is therefore not quite apparent.
Here, we probe in detail the implications of a generality in the St-Re-relationship on the locomotory energetics and thrust generation mechanics of differentially constrained self-propelled pitching foils. Through detailed simulations and in-depth analysis of self-propelled pitching foils with distinct translational restraints, we demonstrate how a difference in the kinematic constraints leads to a profound disparity in the thrust generation mechanics. Notwithstanding this disparity, we find convergence to surprisingly similar St-Re power laws over an appropriate parametric space characterized by large reduced frequencies (reduced frequency k = πfc/ū p ). Specifically, we find that the St-Re relations for the dissimilarly constrained flapping foil states are generally quite distinct and prescribed kinematics dependent. However, for sufficiently large reduced frequencies, these distinct relations converge to similar power laws with a unified scaling exponent. Most importantly, we find that the input power requirement is extremely sensitive to the translational restraints. Hence, strikingly divergent, far from universal trends for the energetic cost of locomotion are obtained for self-propelled flapping foil states that are governed by a unified St-Re relationship. Our investigation is centred around self-propelled foils that pitch either in isolation or in combination with an induced passive heave. Configurations similar to ours but consisting of tethered foils undergoing pitch and heave in a uniform cross-flow have been extensively scrutinized in the recent past (see the review by Smits 2019). Much of the current understanding of a convergence in the St-Re dependence is based on the idea of a balance between the predominantly viscous resistance to locomotion from a laminar or turbulent boundary layer, and a counterbalancing thrust conjectured to arise principally from the fluid's inertial reaction to the periodic transverse motion. Our analysis supports a significant departure from this popular view in that a convergence to a unified St-Re power law occurs only in the large reduced frequency limit. It is only in this large reduced frequency limit that the thrusts generated from distinct leading edge suction and added mass related mechanisms assume a considerably simpler and familiar, forcing frequency and amplitude squared dependent form. The identification of this distinction in the origin of thrust and its non-trivial implications on the energetic cost of locomotion distinguish our work from the earlier ones. Our present analysis focuses exclusively on self-propelled pitching foils endowed with varying levels of translational freedom. Nonetheless, given the morphological and dynamical similarity between swimming fish and typical foils (Webb 1975;Triantafyllou et al. 2005;Lucas, Lauder & Tytell 2020), we anticipate that the distinction in thrust generation mechanics and the extreme sensitivity of energetics to the locomotion-enabling kinematics will be observed over a wide spectrum of rigid and flexible flapping foil self-propulsion and undulatory biolocomotion.
the frequency of the imposed rotational pitch and time, respectively. We define the inline x and the transverse y directions along and perpendicular to the foil's chord when the foil is at its mean position, respectively. We denote the coordinates of the pivot point located at the quarter chord by (x p , y p ).
The three configurations are distinguished by the translational degrees of freedom endowed to the self-propelling foil. In a bidirectionally free (BF) configuration we prescribe only the rotational pitch and allow the foil to translate freely in response to the unsteady thrust/drag and lift forces. In the BF configuration we therefore determine the foil's translational motion (forward surge and transverse heave) from the total force (thrust/drag and lift) exerted on it. The unidirectionally free (UF) configuration is similar to the BF configuration except that we prevent a transverse heave motion altogether by imposing a restraint y p (t) ≡ 0. Thus, in the UF configuration we prescribe a harmonic rotational pitch and impose a zero transverse motion while determining the foil's forward surge from the total thrust force exerted on it.
For completeness, we consider a third steered (S) configuration wherein we impose a harmonic rotational pitch about the pivot point while constraining it to undergo a uniform rectilinear motion in such a way that the cycle-averaged drag (or thrust) force exerted on the foil vanishes identically. In the S configuration therefore, the foil is towed at a constant forward speed so thatẋ p is thus a constant and y p (t) ≡ 0, while the foil is pitched such that the cycle-averaged drag/thrust force vanishes identically. Owing to Galilean invariance, the configuration S is equivalent to and reminiscent of the commonly investigated canonical configuration consisting of uniform flow past a foil that is pitched sinusoidally about its quarter chord (Das, Shukla & Govardhan 2016;Floryan et al. 2017). The arrangement is such that the foil's mean chordwise direction is aligned with the free stream. In this set-up, for fixed pitching amplitude θ 0 and a given uniform free stream u p , the hydrodynamic force exerted on the foil depends strongly on the pitching frequency. Specifically, for pitching frequencies below a critical frequency f , the foil experiences a net drag force. In contrast, a mean thrust is exerted on the foil when the forcing frequency exceeds f . Precisely when the frequency of the imposed pitch equals f , a mean self-propelling state is established. In this mean self-propelling state, the cycle-averaged drag force exerted on the foil vanishes identically (as does the cycle-averaged thrust force generated from the foil's pitching motion). The frequency f associated with this mean self-propelling state is unique in the sense that for a given combination of θ 0 and u p , the self-propelling state is achieved only when the frequency of the imposed pitch equals f . Therefore, in all the three S, UF and BF configurations, u p =ẋ p and f are interdependent.
In all the three configurations, we set an initial state that corresponds to a foil pitched impulsively in a quiescent fluid. The foil's pivot is located initially at the origin. We solve for the foil's position and velocity at subsequent times by computing its response to the hydrodynamic forces exerted on it. To this end, we rely on a high fidelity Lagrangian viscous vortex particle method (Cottet & Koumoutsakos 2000;Eldredge 2007), the details and validation tests can be found in our previous works (Das et al. 2016(Das et al. , 2019.
In both the BF and UF configurations, the foil accelerates initially and eventually attains a self-propelled state corresponding to a time-invariant cycle-averaged forward speed. Left frames of figure 1 illustrate this process for θ 0 = 5 • and a frequency that results in a forcing Reynolds number, Re f = fA f c/ν = 310, based on the maximum pitching speed at the trailing edge and the foil's chord length c as the characteristic velocity and length scales.
Here, A f = 1.5c sin θ 0 denotes the trailing edge excursion due to the imposed rotational pitch. Thus, A = A f in UF and S configurations. However, owing to a finite induced passive heave, A differs from A f in the BF configuration. The top left frame of figure 1 depicts the temporal evolution of the foil's forward and transverse velocity components u p =ẋ p and v p =ẏ p , where a dot denotes time derivative, in each of the three configurations. In the S configuration a mean self-propelling state is established at a towing speed that corresponds uniquely to the prescribed pitch. This unique towing speed overlaps perfectly with the mean forward speed of a UF foil. The moment M and the input power P per unit span (P = −Mθ − F xẋp − F yẏp , F x and F y being the instantaneous thrust and lift forces per unit span) in the S and UF configurations are remarkably similar as well, despite the presence of significant inline oscillations (fluctuations in u p ) in the latter. The insensitivity of the foil's dynamics to the fluctuations in u p is commensurate with the insensitivity of the propulsive attributes of thrust-generating foils to the streamwise oscillations introduced in an otherwise uniform incoming flow (Van Buren et al. 2018).
The foil's dynamics in the UF and BF configurations differ strikingly. A UF foil accelerates faster than a BF foil, attains a mean self-propelling state in far fewer cycles and exhibits pronounced inline fluctuations. Crucially, the foil's mean forward speed in the UF configuration is over three-fold higher than in the BF configuration. In the UF configuration v p ≡ 0, as the translational motion in the transverse direction is forbidden by the constraint y p (t) = 0. The BF foil exhibits prominent transverse oscillations, the corresponding fluctuations in v p are comparable toū p . Its transverse velocity component v p settles rather rapidly exhibiting no significant cycle-to-cycle variations right from the start. Notably, v p exhibits a cyclic variation that except for being completely out-of-phase is exactly analogous to the sinusoidal variation of the imposed rotational pitch. An instantaneous, direct correspondence between the induced passive heave and the imposed rotational pitch is indicative of the inviscid origin of the foil's transverse motion. The moment and power for BF and UF configurations differ strikingly, as evident from the middle frames of figure 1. Both M and P are however of similar functional form and 930 A27-5 exhibit a periodic temporal variation. The instantaneous cycle-to-cycle invariance of M and P points to their inviscid origin as well.
The out-of-phase passive heave of a BF foil counteracts the influence of imposed pitch, effectively reducing the trailing edge excursion and transverse speed. Consequently, the size and the intensity of the vortices shed into the wake of a BF foil are appreciably lower, as clearly evidenced from the right frames of figure 1. The vorticity distributions in the UF and S configurations are very nearly indistinguishable and have been depicted through a single top right frame in figure 1. Understandably, a similarity in foil dynamics in the UF and S configurations translates into a similarity in the vorticity distribution.
The disparities in the initial acceleration and the mean forward speed to which the UF and BF foils eventually settle are a direct consequence of a dissimilarity in the thrust generation mechanics. This dissimilarity is clearly evidenced from the striking contrast in the cycle-averaged pressure distribution on the BF and UF self-propelled foils (right frames of figure 1). The transverse heave of a BF foil results in a relatively large effective angle between the foil and the direction of locomotion (x direction). The periodic variation in the effective angle of attack alters the pressure distribution in such a way that a sharp drop in the mean pressure is encountered at the foil's leading edge. The thrust generated from this reduction in mean pressure (the so-called leading edge suction mechanism, Garrick (1937)) facilitates propulsion of a BF foil. Owing to a significantly diminished leading edge excursion, no such drop is evidenced in a UF configuration. Consequently, the inertial reaction of the fluid is principally responsible for self-propulsion of a UF foil. This contrast in the mechanics of thrust production is quantitatively established through a minimal model for self-propelled foils in the forthcoming § 3.

The St-Re relationship and the drag-thrust balance
To link the mean self-propelling speed to the prescribed rotational pitch, we perform runs over a wide range of forcing Reynolds numbers (10 ≤ Re f ≤ 1000) and angular amplitudes 5 • ≤ θ 0 ≤ 16 • for each of the three configurations. Our principal observations deduced in the last section for the specific case of θ 0 = 5 • and Re f = 310 generalize to other Re f and θ 0 as well. Specifically, at each Re f and θ 0 , the mean self-propelling speeds attained in the UF and S configurations overlap and consistently exceed the mean self-propelling speed attained in the BF configuration by over two folds.
We next express the dependency ofū p on the key kinematic parameters through the St-Re relationship. For fair comparison, we define the Strouhal number St = fA/ū p based on the actual trailing edge excursion as the characteristic length scale. As noted earlier, Re =ū p c/ν. Thus, owing to an explicit dependence onū p , both Re and St are in fact output quantities, the forcing Reynolds number Re f being the key input quantity. In figure 2 we depict the St-Re dependence for all our runs spanning the entire parametric space of Re f and θ 0 in the BF, UF and S configurations. We observe an expected near-perfect overlap in the St-Re relation for the UF and S configurations. Crucially, we find a convergence in the Strouhal number's dependence on the Reynolds number in that St ∼ Re n , where the scaling exponent n = −0.375 is invariant across the UF, BF and S configurations. The emergence of a power-law scaling was reported in our earlier, comparatively limited parametric space (θ 0 = 5 • ) investigation of the S (Das et al. 2016) and BF (Das et al. 2019) configurations.
We note here that our scaling exponent n = −0.375 differs from the exponent of −0.25 reported in the work of Gazzola et al. (2014). As mentioned in our previous work (Das et al. 2016), a disparity in the set-ups analysed is a likely reason for the aforementioned difference between the scaling exponents. We investigate pure rigid body motions that 930 A27-6 Downloaded from https://www.cambridge.org/core.
StRe −0.56 , is shown using thin solid lines. Thick solid lines depict θ 0 -independent power law that is obtained in the limit of large reduced frequencies (k 1). The drag-determining constant γ = 24 over the entire parametric space and across the three configurations.
arise from a combination of the free translation and prescribed time-periodic rotation. In contrast, Gazzola et al. (2014) investigate biolocomotion wherein a lateral motion induced reduction in boundary layer thickness could potentially be suppressed through travelling-wave-like undulatory motion of flexible tails. Nevertheless, our scaling exponent n = −0.375 is in excellent agreement with the power-law fit of −0.4 reported in the recent independent work on BF foils by Lin et al. (2021). Notably, using StRe = Re f , the St-Re power law can be recast into an equivalent form Re ∼ Re 1.6 f , where the scaling exponent of 1.6 deduced from our runs is in excellent agreement with the exponent of 5/3 reported in Lin et al. (2021). Furthermore, Lin et al. (2021) observe that the power-law relationship of the form St ∼ Re −0.4 (equivalent to Re ∼ Re 5/3 f ) fits the locomotory characteristics of freely translating heaving foils investigated by Alben & Shelley (2005) and Hu & Xiao (2014), as well as those of several natural aquatic swimmers (see figure 6 in Lin et al. 2021). Therefore, the power law St ∼ Re −0.375 very likely possesses a much wider range of applicability that extends to other self-propelled flapping foil configurations and natural aquatic swimmers as well.
The trends illustrated in figure 2 unequivocally point to a generality in the scaling exponent n across the three configurations. The precise value of the scaling exponent n that fits the discrete data set in figure 2 so closely has been shown to arise naturally from a balance between the mean drag coefficientC D ∼ √ StRe −0.56 and the cycle-averaged thrust coefficientC T ∼ St 2 in the configuration S (Das et al. 2019) as where ρū 2 p c/2 has been used for non-dimensionalization of thrust/drag. The specific formC T ∼ St 2 arises from a ( fA) 2 dependence of the inertial thrust conjectured to originate from the fluid's reaction to the transverse foil motion (Gazzola et al. 2014). The expression for the cycle-averaged drag coefficient used in (3.1) follows directly from the Bone-Lighthill boundary layer thinning hypothesis (Lighthill 1971) for enhanced viscous resistance experienced by a swimming body that undergoes significant lateral movement perpendicular to itself (Das et al. 2016). Specifically, the cycle-averaged drag coefficient Das, R.K. Shukla and R.N. Govardhan in (3.1) is deduced from the general expressionC D ∼ √ StC D0 , where a √ St dependence accounts for the boundary layer thinning due to transverse body motion (Ehrenstein, Marquillie & Eloy 2014). The term C D0 represents the drag coefficient associated with a similar body that is held stationary in a uniform cross-flow. Specifically for a NACA0012 foil, C D0 = 5.91Re −0.56 (Das et al. 2016). Computational results confirming the Bone-Lighthill boundary layer thinning hypothesis for uniform flow past two-and three-dimensional finite-span plates undergoing significant time-oscillatory wall-normal motion have appeared in Ehrenstein & Eloy (2013) and Ehrenstein et al. (2014). Clearly, the scaling exponent n conforms to the St-Re discrete data set for the UF and BF configurations equally well. Hence, a heuristic extension of the foregoing arguments to UF and BF self-propelled foils is particularly tantalizing.
A generalization of the above arguments to a UF configuration is justified given the insensitivity of UF foil's dynamics to inline oscillations. Generalization to a BF configuration is however dubious owing to the stark disparity in the mean pressure distribution, the principal source of thrust, over the BF and UF/S foils (see § 2 and figure 1). Furthermore, the St-Re relationship for a BF configuration is discernibly offset from the St-Re relationship for the UF and S configurations (figure 2). Unlike the scaling exponent, this offset can not be explained from a straightforward generalization of the scaling arguments from the foregoing paragraph. We emphasize here that the offset is not a consequence of the use of actual trailing edge excursion in the definition of St. The offset is in fact enhanced with a switch to the forcing amplitude A f based definition of the Strouhal number (A/A f ≈ 0.67 as shown in the forthcoming analysis).
To unravel the reasons underlying the offset in the St-Re relation, we analyse the mechanics of thrust production in UF and BF configurations by developing a minimal model for self-propelled foils wherein the pitch induced passive heave is either prevented (UF configuration) or determined by coupling the BF foil's transverse motion to the lift force deduced from the linear theory (Garrick 1937;von Kármán & Sears 1938;Fernandez-Feria 2016). The complete details of our rigorous linear theory based modelling approach for prediction of transverse foil dynamics and thrust generation mechanics associated with BF/UF configurations are provided in Appendix A. The analysis yields the following theoretical estimate for the transverse dynamics of the foil's pivot point (see Appendix A.1): Hereỹ p and φ h denote the amplitude and phase with respect to the prescribed pitch (y p = Here reduced frequency k = πfc/ū p and Theodorsen function C(k) = F(k) + iG(k).
In figure 3 we compare the theoretical predictions from (3.2a,b) with the heave amplitude and phase obtained from the simulations. The k-dependence in the heave amplitude and phase given by (3.2a,b) arises principally from a Coriolis-like term that is significant only at low reduced frequencies and is 90 • phase offset with respect to θ . Without this Coriolis term, the foil's transverse motion is k-independent and arises solely from the added mass effects. Thus, without the Coriolis term, or equivalently when k 1,ỹ p ≈ 0.16A f , φ h = 180 • (shown by dotted lines in figure 3) and A ≈ 0.67A f so that the induced heave and imposed pitch are completely out-of-phase. The dominance of the Coriolis effect at progressively lower reduced frequencies leads to an increasingly prominent passive heave that lags the imposed pitch by a monotonically decreasing phase φ h . The agreement between the simulation results and predictions from the linear theory  based model that fully accounts for the added mass and Coriolis effects is quite remarkable over the entire parametric space of θ 0 and Re. The agreement indicates that the BF foil's passive heave arises principally from a recoil reaction to the imposed pitch. With a reasonably accurate predictive model for transverse motion in place, we again appeal to the linear theory for an estimate of the mean thrust coefficient. For a foil undergoing simultaneous pitch and heave motions, as shown in Appendix A.2, where the reactive (added mass related), quasi-steady, wake and leading edge suction contributions to the mean thrust coefficient are given bȳ Clearly β = 0 for UF/S foils. An evaluation of the above contributions in the large reduced frequency limit yields C T ≈C R T = π 3 St 2 /9 in the UF and S configurations. Here the quasi-steady, wake and leading edge suction terms contribute negligibly and the mean thrust is generated principally from the fluid's reaction to the imposed pitch (see Appendix A.3 for details). Contrastingly, in a BF configurationC T ≈C LS T = π 3 St 2 /16 with the rest of the terms contributing negligibly to the mean thrust coefficient (see Appendix A.4). Thus, the reaction from the fluid ceases to be a significant source of thrust and rather, consistent with the pressure distribution contrast depicted in figure 1, the thrust produced via leading edge suction mechanics self-propels a BF foil.
The reason behind the equivalence in the scaling exponent and the offset between the St-Re power laws for the BF and UF configurations is now unravelled. The balance between the mean thrust coefficient from (3.3) and a cycle-averaged drag coefficient that 930 A27-9 A. Das, R.K. Shukla and R.N. Govardhan is consistent with the Bone-Lighthill boundary layer thinning hypothesis, γ being a drag-determining parameter, allows us to deduce a general k-dependent St-Re relation. This general relation is θ 0 specific, but for sufficiently large k (or equivalently large St), simplifies considerably to a unified θ 0 -independent power law across UF, BF and S configurations as Re −0.375 for a UF/S foil, (3.7) in the large reduced frequency limit (k 1). The similarity in the scaling exponent across the three configurations is thus a direct consequence of an equivalence in the high reduced frequency scaling of the thrust coefficient (C T ∼ St 2 for k 1 in UF, BF and S configurations). Importantly, the difference in the magnitude of mean thrusts (proportionality constant of π 3 /16 vs π 3 /9) produced in the BF and UF configurations is directly responsible for the appreciable offset in the St-Re relation. As evident from figure 2, the general k-dependent St-Re relation with just a single adjustable parameter γ fits the discrete data set from our simulations reasonably well and crucially explains both the offset between the St-Re relations for BF and UF/S foils, and the divergence from the unified power law at low k. Evidently, such low k, or St, is achieved at the highest Reynolds numbers (Re 500 as inferred from the inset of figure 2).

Energetic cost of self-propulsion
The impact of the passive heave on the energetic cost of self-propulsion is significantly more pronounced. In figure 4 we compare the mean power coefficient (C P = 2P/(ρū 3 p c)) for all our runs over the three configurations. Owing to the form of normalization employed, the power coefficient can also be viewed as a scaled cost of transport. At low Re (Re 200), the mean power coefficient for the BF configuration consistently exceeds power coefficients for UF and S configurations by a factor that varies from two to three (left frame of figure 4). The power coefficients for the UF and S configurations are expectedly in mutual agreement.
To linkC P and imposed kinematics, we illustrate the dependence ofC P on St in the right frame of figure 4. We observe an effective convergence to a unified power lawC P ∼ St 3 across the UF, BF and S configurations. The Re-dependent disparity in the power requirement, as evidenced from the left frame of figure 4, is therefore a direct consequence of a cubic amplification of the comparatively small offset in the St-Re relation (figure 2). To self-propel at a givenū p , a BF foil must overcome the counterproductive out-of-phase passive heave through a more intense rotational pitch (higher θ 0 and/or f ). The increase in pitching intensity significantly augments the mean energetic cost of self-propulsion.
Foregoing arguments hold reasonably well in the low Re regime wherein k 1. In this case, the imposed pitch and the passive heave remain nearly out-of-phase (φ h ≈ 180 • ).  highest θ 0 and/or f . At high θ 0 and/or f , with further increase in the pitching intensity, Re rises whereas St and k diminish considerably. At sufficiently high Re, or equivalently low k, owing to the dominance of the Coriolis-like term,ỹ p grows while φ h diminishes monotonically towards 90 • (see figure 3). In this scenario, a BF foil's excursions from its instantaneous direction of motion are reduced and its cost of transport is diminished. At sufficiently low k,C P for a BF foil reduces to such an extent that it is overwhelmed by theC P for a UF foil. The induced passive heave thus plays a contrasting role in that it enhances locomotory performance of a BF foil at high Re (low k), while reducing it at low Re (high k). This contrast in the role of passive heave at low and high k is consistent with the enhancement of thrust generated from tethered foils undergoing combined pitch and heave motions at φ h = 90 • , and its deterioration at φ h = 180 • (Anderson et al. 1998;Smits 2019). Our results point to the possibility of achieving an exceptionally efficient fish-like locomotory state in a relatively simple self-propelled system. Remarkably, the induced heave in this BF self-propelled system is passive and uncontrolled and, thus, not tuned to minimize the energetic cost of locomotion.

Conclusions
In summary, through detailed simulations, we demonstrated that bidirectionally and unidirectionally free self-propelled pitching foils exhibit a stark distinction in the thrust generation mechanics and locomotory performance, and yet are governed by a remarkably similar St-Re scaling. Our drag-thrust balance based analysis revealed a similarity in the high reduced frequency scaling of the thrust generated from the distinct propulsion-enabling mechanisms in bidirectionally and unidirectionally free configurations (leading edge suction in BF vs recoil reaction in UF) to be the principal cause of the similarity in the scaling exponent. In general, the St-Re dependence is therefore prescribed kinematics dependent, and a reduction to a previously reported, considerably simpler power-law form with a unified scaling exponent can only be expected in the large reduced frequency limit. Crucially, we established that an amplification of the offset in the St-Re relation, caused by the St 3 dependence of the scaled cost of transport, is directly responsible for the contrast between the locomotory performance of UF and BF foils.

A27-11
Being the principal determinant of the scaling exponent n in the St ∼ Re n power law, drag plays an important dynamical role in the differentially constrained configurations investigated in our work. This conclusion is consistent with the previous ones that have noted the importance of drag in dictating the optimal Strouhal number for flexible locomotion (Godoy-Diana & Thiria 2018) and maximization of thrust generation efficiency from tethered flapping foils (Floryan, Van Buren & Smits 2018). Here we stress on the critical role of thrust in introducing an offset in the St-Re dependence, and through its amplification, being ultimately responsible for the large disparity in the locomotory performance of differentially constrained foils. We anticipate that this key conclusion from our work will generalize to other forms of oscillatory flapping foil propulsion and undulatory biolocomotion. Specifically, the extreme sensitivity to variations in thrust will prevail in artificial as well as natural undulatory swimmers so that any deviation from a general St-Re trend in the form of an offset or general scatter will have enormous implications for the energetic cost of locomotion.
Funding. We acknowledge support received from NVIDIA Corporation (hardware donation program) and National Supercomputing Mission of the Department of Science and Technology (runtime on cray XC-40 housed in Supercomputing Education and Research Center-Indian Institute of Science).

Declaration of interests. The authors report no conflict of interest.
Author ORCIDs.

Appendix A
Here we provide details of the linear theory based model for predicting the thrust generated from unidirectionally and bidirectionally free foils that are pitched about their quarter chord. A BF foil can undergo lateral heave motion in addition to its forward surge. To begin with therefore, we analyse the imposed pitch induced passive heave in a BF foil.
A.1. Model for prediction of imposed pitch induced passive heave of a BF foil Apart from being almost out-of-phase, the transverse heave of a BF foil very nearly follows the imposed harmonic pitch. This combined with the fact that the transverse heave of a BF foil exhibits no significant cycle-to-cycle variation right from the start (as shown in figure 1) points to the inviscid origin of the induced lateral heave motion. Linear theory (Theodorsen 1935;Garrick 1937;von Kármán & Sears 1938;Fernandez-Feria 2016, 2017Mackowski & Williamson 2017) provides a natural framework for an inviscid analysis of thin foils undergoing combined pitch and heave motions. We therefore appeal to the inviscid linear theory for prediction of transverse dynamics of a BF foil.
We consider an upward linear motion and counterclockwise rotational motion of the foil to be positive. Utilizing this sign convention, the lift force per unit span, exerted on the foil, as deduced from linear theory is given by where and C(k) = F(k) + iG(k) is the Theodorsen function. The reduced frequency k = πfc/ū p , whereū p denotes the mean forward velocity of the foil with a as the pivot location normalized with respect to the foil's half-chord length c/2. For pitching motion about the quarter chord therefore, a = −1/2. The instantaneous heave is characterized using the vertical displacement of the pivot point y p =ỹ p exp(i(2πft + φ h )), where f andỹ p represent the frequency and amplitude of the induced passive heave. The instantaneous angular position of the foil is given by θ = θ 0 exp(i2πft), where θ 0 denotes the amplitude of the imposed rotational pitch. Given the insensitivity of the foil's dynamics to the inline oscillations (see § 2), we employ only the mean forward speedū p in our estimate (A1) of the lift force.
The transverse heave of a BF foil is a direct consequence of the instantaneous lift force exerted on it. We therefore obtain where m and y cm denote the foil's mass per unit span and the y-coordinate of its geometric centre, respectively. For small pitching angles, an assumption inherent in the linear theory, we can relate y cm to y p as y cm = y p + r cm/p θ andÿ cm =ÿ p + r cm/pθ , (A4a,b) where r cm/p denotes the distance between the geometric centre and the pivot point.
A rearrangement yields Substitution of the specific forms for pitch and heave allows us to deduce that A normalization using A f = 1.5θ 0 c, the trailing edge excursion due to the imposed pitch, yields the following for the non-dimensional amplitude of the transverse harmonic motion of the pivot point: Here r p c = 0.75c denotes the distance between the trailing edge and the pivot point. Substituting parameters specific to a NACA0012 foil, ignoring the influence of the foil's mass m (m/ρc 2 = 0.08 1), and utilizing r cm/p = 0.17c for the configuration under investigation, we obtain a 1 = −k 2 − 2kG(k), a 2 = 2kF(k), b 1 = 0.5k 2 + 2kG(k) − 2F(k) 930 A27-13 and b 2 = −2kF(k) − k − 2G(k). An expression for the trailing edge excursion of a BF foil can now be derived as Furthermore, We note here that our present model, developed within the rigorous framework of linear theory, is far more generic than the model developed in our previous work (Das et al. 2019). Our present model fully accounts for the leading edge suction mechanics, as well as the added mass effects, and also incorporates contributions from the quasi-steady terms and the wake. Unlike the model proposed in Das et al. (2019), our present linear theory based model remains applicable over a wide range of reduced frequencies and accurately predicts the strong k-dependence of the amplitude and phase of the induced passive heave in the BF configuration (see figure 3). It is worth noting that in the limit k 1, the expression (A10) simplifies to A ≈ 2A f /3.

A.2. Thrust generation mechanics of UF and BF self-propelled foils
Having characterized the passive heave of a BF foil, we next extend our analysis to predict the propulsion-enabling thrust force exerted on the BF and UF foils. Here, unlike several previous works on similar configurations (Gazzola et al. 2014;Das et al. 2016Das et al. , 2019, for instance), we do not make any intrinsic assumption about the reactive origin of the thrust. Rather, the determination of the mean thrust is tightly coupled to the theoretical prediction of the lateral foil motion from the previous subsection.
The instantaneous thrust force experienced by a foil undergoing simultaneous pitch (θ = θ 0 exp(i2πft)) and heave (y p =ỹ p exp(i(2πft + φ h ))) is given by Normal force contribution to thrust Leading edge suction contribution to thrust , where L denotes the normal force component oriented in a direction perpendicular to the foil's chord with S as the leading edge suction term. The normal force (L) and leading edge suction term (S) are given by

Contrasting dynamics and energetics of self-propelled foils
contribution from the wake L W . Thus, where Defining the thrust coefficient as C T = 2T/(ρū 2 p c) and using the split form of the lift force (A13), we deduce the following for the thrust coefficient: Here C R T , C QS T , C W T and C LS T denote the contribution to the thrust from the reactive normal force, quasi-steady circulation, wake and leading edge suction terms, respectively.
The reactive normal force contribution to the mean thrust coefficient is Substitution of the expression for L R from (A14a-c) yields To make a further simplification, we make use of the Strouhal number St = fA/ū p , where A represents the trailing edge excursion of the foil. For similar pitching amplitude and frequency, the trailing edge excursions differ for UF and BF free foils. Setting A f = 3αA and β = 2ỹ p /(θ 0 c), the reactive contribution to the thrust coefficient assumes the following general form: The three other constituents of the thrust coefficient are given bȳ  A.3. Thrust coefficients in a UF configuration In a UF configuration β ≡ 0 and α ≡ 1/3. Substituting these parameters in (A18) and (A19), we deduce the following for the constituents of the mean thrust coefficientC T : In the large reduced frequency limit, the above expressions forC T in a UF configuration simplify considerably. Specifically, making use of G(k) ≈ 0 and F(k) ≈ 1/2 for k 1, we obtainC R T = π 3 St 2 /9 withC QS T ≈ 0,C W T ≈ 0 andC S T ≈ 0 for a UF foil. In the top left frame of figure 5 we depict the dependence of thrust coefficient for a UF configurationC T on St, as predicted from (A15) and (A20), for θ 0 = 5 • , 8 • , 12 • and 16 • . We find an effective convergence to the asymptotic limit ofC T = π 3 St 2 /9 for St 0.3. For St 0.3, k ∼ O(1), and the general expression (A15) and (A20) exhibits significant θ 0 -dependent deviations from the k 1 asymptotic limit ofC T = π 3 St 2 /9. The bottom left frame of figure 5 depicts a comparison of the individual contributions from the reactive, quasi-steady, wake and leading edge suction terms to the net thrust coefficient. We find the reactive termC R T to be the dominant source of thrust and more so at 930 A27-16 large reduced frequencies, in accordance with the foregoing prediction from the theoretical analysis. Note that the ratiosC R T /C T ,C QS T /C T ,C W T /C T andC LS T /C T are all independent of θ 0 and, therefore, their individual dependence on k has been illustrated through a single trend line in the bottom left frame of figure 5.
A.4. Thrust coefficients in a BF configuration Substitution of (A7a,b) in (A18) and (A19) yields the following for the thrust components in a BF configuration: Making use of G(k) ≈ 0 and F(k) ≈ 1/2 for k 1, we obtainC R T ≈ 0,C QS T ≈ 0,C W T ≈ 0 andC LS T ≈ π 3 St 2 /16 for a BF foil. The top right frame of figure 5 depicts the dependence of thrust coefficient for a BF configuration on St, as predicted from (A15) and (A21), for θ 0 = 5 • , 8 • , 12 • and 16 • . We find an effective convergence to the asymptotic limit ofC T = π 3 St 2 /16 over the entire range 0.04 < St < 10. In the BF configuration the thrust does exhibit θ 0 dependence. The convergence observed in the top right frame of figure 5 is simply due to the normalization with respect to the actual trailing edge excursion A.
The bottom right frame of figure 5 depicts a comparison of the individual contributions from the reactive, quasi-steady, wake and leading edge suction terms to the net thrust coefficient. We find the leading edge suction termC LS T to be the dominant source of thrust and more so at large reduced frequencies. This observation is again consistent with the foregoing theoretical analysis. As in the case of the UF configuration, the ratiosC R T /C T , C QS T /C T ,C W T /C T andC LS T /C T are all independent of θ 0 and, therefore, their individual dependence on k has been illustrated through a single trend line in the bottom right frame of figure 5.
Clearly, the rigorous framework of linear theory utilized in our present model allows for unambiguous discrimination between the distinct origins of thrust. Notably, our present linear theory based model identifies leading edge suction and the fluid's inertial response (added mass) as distinct sources of thrust production in BF and UF/S configurations, respectively. This distinction in the source of thrust in BF and UF/S configurations is consistent with the predictions from our detailed simulations (see figure 1 and