Marie Rastello
email: marie.rastello@legi.cnrs.fr

Jean-Louis Marié

Wake behind contaminated bubbles in a solid-body rotating flow

Keywords:

Bubbles are present in many environmental and industrial systems. Their surface in contact with the liquid enables transfers (of heat and mass) with the bulk (Clift, Grace & Weber 1978), which finds its application in numerous processes (waste water treatment, chemical reactors, etc.). Besides, by their dynamics, they may improve such

A bubble injected into a flow rotating about a horizontal axis comes to an equilibrium location. The drag and lift exerted on the bubble can be measured and the bubble wake visualized (Rastello et al., J. Fluid Mech., vol. 624, 2009, pp. 159-178;Rastello et al., J. Fluid Mech., vol. 682, 2011, pp. 434-459;Rastello et al., J. Fluid Mech., vol. 831, 2017, pp. 529-617). For a contaminated bubble, interface deformation remains limited. The bubble is freely rotating, which results in a complex separated wake, influenced by rotational flow and bubble spinning. The wake is described by analysing the near-and far-wake geometry and behaviour from laser-sheet visualizations, as a function of the relevant non-dimensional numbers: bubble Reynolds number Re, Rossby number Ro, and non-dimensional spinning rate Ω * . As the far-wake length increases with Re, it deflects towards the rotation axis of the flow, the deflection angle increasing with Re and being twice the angle that would occur without deflection. Deflection is stronger for bubbles located close to the rotation axis of the flow (small Ro). The far wake is more curved than the incoming streamlines. The near wake exhibits three distinct regimes as a function of Re. For Re 140, the near wake is structured by the bubble spinning. Its size is related to Ω * and grows faster with Re than for a stationary sphere in a uniform flow. As the bubble spinning is saturating (140 < Re 240), these differences vanish and Re dependences for the two situations becomes comparable. For Re > 240, wake instability generates a bubble precession that makes the near wake decrease rapidly for higher Re. These regimes coincide with the changes in the lift coefficients that we have noted in our studies. transfers and mixing in the liquid through various mechanisms [START_REF] Risso | Agitation, mixing, and transfers induced by bubbles[END_REF]. In most situations, bubbles are millimetric (Reynolds number Re 200) and the agitation caused by these rising bubbles is to a large extent related to the potential flow disturbance created around the surface and the liquid on average dragged by the wake. Depending on the size, shape and degree of contamination of the bubbles, the wake may be more or less significant and stable, which has an influence on the liquid disturbance generated and the hydrodynamic interactions between bubbles. Knowledge of the behaviour of bubble wakes, more specifically of the near-wake behaviour where most transfers occur [START_REF] Fan | Bubble Wake Dynamics in Liquids and Liquid-Solid Suspensions[END_REF], is therefore of some importance regarding these problems. Also, the development of the far wake is interesting to investigate, since it may affect the bubbles or particles located behind.

We recently performed an experiment on the forces acting on bubbles immersed in a solid-body rotating flow around a horizontal axis (SBRH). In this flow, the bubble comes to an equilibrium location, the coordinates of which can give access to the coefficients of the drag and lift forces experienced by the bubble. Two situations were investigated: the case where the interface is clean (silicone oils [START_REF] Rastello | Drag and lift forces on clean spherical and ellipsoidal bubbles in a solid-body rotating flow[END_REF]), and the case where the interface is partially or fully covered by impurities (water [START_REF] Rastello | Drag and lift forces on interfacecontaminated bubbles spinning in a rotating flow[END_REF][START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF]). The fact that the bubble is stationary is convenient to image the flow around it and in particular the wake behind as seen in figure 1. When the surface is contaminated, and only in that case, it was shown to rotate with characteristics that are very similar to those of solid spheres immersed in that kind of flow [START_REF] Bluemink | A sphere in a uniformly rotating or shearing flow[END_REF][START_REF] Bluemink | Drag and lift forces on particles in a rotating flow[END_REF]. This results in an extra 'Magnus'-like lift force and a separated wake behind the bubble, whose separation angle (the angle from the bubble rear at which the wake detaches from the bubble) is higher than that observed at the same Reynolds number on a solid non-rotating sphere in a uniform flow [START_REF] Johnson | Flow past a sphere up to a Reynolds number of 300[END_REF]. In the literature, the effect of a transverse rotation on the wake structure of a sphere was investigated in uniform flow [START_REF] Giacobello | Wake structure of a transversely rotating sphere at moderate Reynolds numbers[END_REF][START_REF] Poon | Flow past a transversely rotating sphere at Reynolds numbers above the laminar regime[END_REF]) and more scarcely in linear shear flow [START_REF] Kurose | Drag and lift forces on a rotating sphere in a linear shear flow[END_REF][START_REF] Bagchi | Effect of free rotation on the motion of a solid sphere in linear shear flow at moderate Re[END_REF]. Information about this question in solid rotating flows is even more scarce. The available data, ours [START_REF] Rastello | Drag and lift forces on interfacecontaminated bubbles spinning in a rotating flow[END_REF]) and those of [START_REF] Bluemink | Drag and lift forces on particles in a rotating flow[END_REF], give an estimate of the flow region that is affected by the wake and show that it may deflect strongly when the bubble is close to the cylinder axis. However, these results obtained by using particle image velocimetry (PIV) visualization techniques with seeding particles remain limited.

The objective of this paper is thus to provide a more detailed experimental description of the wake characteristics, as a function mainly of the Reynolds number, itself linked to the rotation rate of the tank. This mostly concerns bubbles in water whose surface is partially covered by impurities. A few results on bubbles with 'fully' contaminated surface are given for comparison. For an easy reading, we first briefly present the facility, the experimental procedure and the range of parameters investigated. The details can be found in our previous papers. We then define the wake characteristics on which we focus our study and make precise the way they were determined. The results obtained are analysed in terms of the non-dimensional numbers relevant for the effect we look at. The far wake is first addressed followed by the near wake. The paper ends with a conclusion.

Experimental procedure

2.1. Brief description of the experiment A more detailed description can be found in our aforementioned papers. Readers who would like further information on the device that is used are invited to refer to appendix A, where two technical schematics showing the mechanical components are given. The experiment is conducted in a Plexiglas cylindrical tank (Ø = 10 cm, L = 10 cm) rotating around a horizontal axis (figure 2a). The tank is mounted in a cylindrical counterbore and the contact is made using ball bearings (appendix A). An 80 W brushless electric motor is coupled to the tank using a toothed belt. The motor is driven through a Labview interface so as to obtain the best torque/speed control. This ensures that the assigned rotation speed is well reproduced and constant. This was checked by an optical tachometer. The facility is fixed on a damped heavy metallic casing so that no perturbing vibration is experienced by the flow. The tank is filled with demineralized water characterized by a resistivity of 0.3 M cm. This resistivity is between that of ultra-purified water, 18 M cm (Duineveld 1995), and that of tap water, 3 k cm. Because of operating constraints (bubble injection using a syringe, plunging of a thermometer for temperature measurements), it was extremely difficult to keep this water clean.

F D F L FIGURE 3.
The forces applied to a bubble at equilibrium at moderate to high Reynolds numbers [START_REF] Rastello | Drag and lift forces on interfacecontaminated bubbles spinning in a rotating flow[END_REF]. Here F B is the buoyancy force, F D the drag force, F L the lift force, and F A is the sum of the added-mass and inertial forces. In this diagram, the angular velocity ω and equilibrium angle θ e are negative.

Practically, it contains contaminants that are mainly solid impurities and/or traces of tensio-actives.

When injecting the bubble with the syringe, a first-order estimate of the volume and diameter of the bubble is obtained. The precise determination will then be achieved through the interface detection on the collected images. When the flow is in solid-body rotation at a constant speed, all the forces acting on the bubble balance and the bubble comes to an equilibrium point (figure 3). For the bubbles considered in this study, which have moderate to large Reynolds numbers, this force balance can be written as

F B + F A + F D + F L = 0.
(2.1)

Here F B is the buoyancy; F A is the sum of the added-mass and the inertial forces, the inertial force corresponding to the pressure gradient in the fluid and also called the pressure gradient force; F D is the drag force; and F L the lift. As F B and F A are known analytically, the lift C L and drag C D coefficients can be readily deduced from the measured polar coordinates (r e , θ e ) of the centre of the bubble at equilibrium, the rotation speed of the tank ω and the equivalent diameter d of the bubble. The diameter d is directly the diameter of the bubble for a spherical bubble, while for slightly ellipsoidal ones, d is the diameter of a sphere of identical volume accounting for their oblate spheroidal shape [START_REF] Rastello | Drag and lift forces on clean spherical and ellipsoidal bubbles in a solid-body rotating flow[END_REF]. The aspect ratio χ of the deformed bubbles investigated in this paper (Re 500) never exceeds 1.1 [START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF], making the resulting uncertainty negligible. The bubble centroid is always the barycentric position. The bubble at equilibrium is illuminated from behind and imaged with a Phantom V4.3 high-speed camera positioned in front of the tank (figure 2a).

The non-dimensional numbers used to describe the flow are the Reynolds number Re = |v -U|d/ν and the Rossby number Ro = |v -U|/ωd (equal to the inverse of the local shear rate S) for spherical bubbles, plus the Weber number We = ρ|v -U| 2 d/σ for deformed bubbles. The Morton number Mo = gν 4 ρ 3 /σ 3 for experiments in water is constant, equal to 2.5 × 10 -11 . When the bubble is at equilibrium, its velocity is v = 0 and the liquid velocity it sees is U = ωr e , so that the non-dimensional numbers simplify as Re = ωr e d/ν, Ro = r e /d and We = ρω 2 r 2 e d/σ . The ranges of 1. The rotation rate for the tank is constant, with a value that varies from 6 rad s -1 to 21 rad s -1 , while the equivalent bubble diameter ranges from 0.5 mm to 3.3 mm. Rossby numbers Ro 5 were not investigated, because they correspond to situations where bubbles are too close to the centre of rotation (less than five diameters) and destabilize the rotational flow [START_REF] Rastello | Drag and lift forces on interfacecontaminated bubbles spinning in a rotating flow[END_REF]).

Wake visualization -determination of boundaries

As described in [START_REF] Rastello | Drag and lift forces on interfacecontaminated bubbles spinning in a rotating flow[END_REF][START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF], the flow around the bubble has been visualized by seeding the fluid with micrometric fluorescent polymer particles (Duke Scientific Corporation, size = 3 µm, density = 1.05 g cm -3 ). Particles are illuminated with a laser sheet and imaged with the high-speed camera (figure 2b). Most of them tag the flow around the surface and the wake, while a small fraction collects at the surface where they materialize the rotation of the surface. The frame rate is selected such as to follow the particles rotating with the surface with a good accuracy, for example so that a minimum of 50 images are acquired per full spin. As a consequence, its precise value (between 100 and 200 frames per second) depends on the rotation rate of the tank. A typical example of instantaneous images thus recorded is shown in figure 1(b). Two regions can be distinguished behind the bubble: the near wake located close to the spinning surface, where particle velocity is low and the streamlines not well visible for that reason; and the far wake, where particles having a higher velocity look like small segments and thus better mark the streamlines. This is precisely the criterion that we use to separate the two regions. The visualization of the streamlines in the near wake was made possible by averaging between five and 10 successive images, depending on the sharpness and the contrast of the sequence processed.

The result is shown in figure 4(a). The image has been inverted for better visualizing the streamlines. We see that the averaging enables to follow particles along their trajectories a longer time, hence to better visualize the streamlines and the wake boundaries. As a proof, the average image is compared with the streamlines obtained on the same sequence by using a standard LaVision PIV software (figure 4b). Except for the white region, the agreement between the two panels is particularly good. On the right-hand side of the bubble, the white region corresponds to the shadow visible in figure 4(a). This shadow is caused by the light sheet being intercepted by the bubble's surface (figure 2b). The white region located in the upper right corner is due to a loss in the PIV velocity correlations resulting from the refraction of the laser sheet by the tank's cylindrical shape. Neither of the two regions intersect the wake, and no specific effort was made to analyse them. Average images were therefore used in addition to instantaneous ones to visualize structures in the wake and identify its boundaries. It was not possible to do the same with PIV, because in many sequences concentration in seeding particles was not sufficient to obtain two-frame correlations of good enough quality for all regions of the flow. For each bubble, five representative instantaneous images are extracted from the full recording.

On each image, wake boundaries are determined manually by delimiting both the instantaneous near wake and far wake. An average over the five instantaneous wakes is then performed for both the near and far wakes, giving average boundaries for each wake. The resulting wakes are then superimposed on the average images and PIV when available to confirm the accuracy of the boundaries. The characteristics Mean standard deviation of the surface of the far wake (column 1) and near wake (column 3) compared with the surface of the bubble. Mean standard deviation of the length of the chord for the far wake (column 2) and near wake (column 4) compared with the diameter of the bubble.

that we investigated are presented on average images in the schematic in figure 5.

They are obviously two-dimensional. No three-dimensional information on the wake structure can be inferred from this kind of visualization. An uncertainty analysis was undertaken to quantify the degree of accuracy of the method. The results are listed in table 2.

No important bias in the results is to be expected from the procedure of determination of the wake properties. In contrast, we were not able to estimate properly the separation angle at which the wake detaches from the bubble surface. The reason is that, in the range of Re investigated (mostly Re 50), the points of the surface where it detaches are located very close to the horizontal equatorial plane of the bubble and almost diametrically opposed. The resolution of the visualization is not sufficient to identify their location at the surface precisely. This makes the accurate estimation of the separation angle from the centre of the bubble rather difficult and imprecise. None of the measured wake surfaces suffer from this lack of precision given that only a negligible surface is at stake. The separation angle defined from the rear point where the bubble spinning is maximum (figure 7 in [START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF] typically ranges between 75 • and 90 • .

Results

The spinning of the bubbles partially covered by impurities is non-uniform, with a maximum in the rear region where the wake is separated [START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF]. Their mean spinning rate Ωb is higher than ω and increases with Re, without any clear Ro dependence being visible over the range investigated. Its Re variation is well approximated by

Ωb ω = 1.6 -0.6 exp(-0.011Re), (3.1)
while the maximum spinning rate varies as

Ω m ω = 1 + 0.15Re 0.55 (3.2)
and can be locally much higher than Ωb . The effect of this spinning on the wake characteristics can be described through the non-dimensional rate Ω * , defined as ( Ωb d/2)/U, with U = ωr e . This non-dimensional rate quantifies the importance of the bubble spinning velocity relative to the upstream flow velocity. The bubble spinning is induced by the vorticity of the flow at the equilibrium distance r e to which the bubble stabilizes. This distance depends on the rotation rate of the tank ω and the particle diameter d. As a result, Re and Ro cannot be vary independently but according to a phase diagram (Re, Ro) given in figure 6. We see

Ω * = 1 2 Ωb ω 1 Ro . (3.3)
As shown in figure 7(a), this non-dimensional spinning rate increases from 0.05 up to 0.15 with Re, a range that is more restricted than the ones ([0-1]) generally covered by simulations. However, replacing Ωb by Ω m in (3.3) yields an Ω * that may be locally higher: from 0.07 up to 0.55. The mean dependence of Ω * with Re resembles that of Ωb /ω, with the existence of a plateau at Re > 250. The smallest values of Ω * are obtained when Ro > 8.5 for the few moderate Re numbers with large equilibrium distance r e compared to the bubble diameter. Note that it is more difficult to estimate r e and d from the laser-sheet visualizations than from the shadow technique used for force measurements ( § 2.1). In some cases, the flow centre is not visible and r e then must be inferred from the streamlines in this region like in figure 9(a). On the other hand, the bubble surface looks more or less thick (figure 9d), which renders the accurate estimation of d difficult. The uncertainty on the interface position is typically ±2 pixels for the small bubbles (<1 mm) and ±6 pixels for the big ones (>1.5 mm). This gives an uncertainty for d of the order of 5 % for Re 100 and of 10 % for Re > 100.

All these difficulties influence the estimation of Re and Ro, which contributes to the spreading of the data. However, we tried to minimize the errors on these nondimensional numbers by estimating them otherwise. When a bubble is at equilibrium,

C D = 4 3 gd cos θ e r 2 e ω 2
(3.4) yielding an uncertainty of the order of 20 % for the spherical bubbles and 30 % for the deformed ones. This explains why the data scattering is more accentuated in Ro than in Re.

Examples of bubble wakes for increasing Re at two given Ro values (Ro ≈ 7.3 and Ro ≈ 11) are displayed in figures 9 and 10, respectively.

Wake behaviour -general trends

In a uniform flow, the wake behind a stationary non-rotating sphere separates for Re ≈ 24 [START_REF] Taneda | Experimental investigation of the wake behind a sphere at low Reynolds[END_REF]), becomes non-axisymmetric above Re ≈ 212 and unsteady for Re > 274 [START_REF] Johnson | Flow past a sphere up to a Reynolds number of 300[END_REF]. Here separation also occurs at small Re (figure 10a), but the near and far wake in the steady regime are not axisymmetric. The high-velocity flow arriving on the right-hand side of the bubble rotates with the spinning surface around the bubble rear and moves upwards on the left side (figure 10d) till it is pushed downwards by the downstream flow on the low-velocity side, there forming a high-shear layer. This causes the stagnation point observed behind the bubbles in the clean situation (figure 1a) to displace upwards on the low-velocity side of the bubble. This trend clearly appears on the PIV data in figure 4(b) and is consistent with those reported in a uniform flow [START_REF] Giacobello | Wake structure of a transversely rotating sphere at moderate Reynolds numbers[END_REF][START_REF] Poon | Flow past a transversely rotating sphere at Reynolds numbers above the laminar regime[END_REF]. On the other hand, separation angles are much higher than those observed at the same Re on solid spheres in uniform flows [START_REF] Johnson | Flow past a sphere up to a Reynolds number of 300[END_REF]. As explained in [START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF], this may result from the fact that bubbles are partially contaminated and spinning. It is worth noting that it is difficult to accurately delimit the regions or recirculation in the near wake, unlike in simulations. Whatever Ro, we see that the increase of Ω * with Re increases the curvature of the streamlines in the near wake. At the same time, the near wake is reduced and even vanishes at high Re and Ω * . The far wake, whose length increases with Re, deflects more or less towards the rotation axis. For a given Re (compare figure 9 10(d) or 11 for a direct comparison), the deflection is more pronounced for the small Ro (small distances from the centre). This can be attributed to two effects. One is the bubble spinning, since Ω * is then higher (figure 7a). The other is the radial pressure gradient ∂P/∂r, which is also higher at small distances (figure 7b). In a rotating flow, ∂P/∂r is equal to ρu 2 θ e /r e , with u θ e the angular velocity. This means that, for a given pressure gradient, the retarded fluid in the wake adopts a smaller radius [START_REF] Bluemink | Drag and lift forces on particles in a rotating flow[END_REF]) and, as a result, the wake may deflect towards the cylinder axis. Our few wake visualizations in silicone oils (non-spinning bubbles show that this effect is indeed visible when a bubble is not spinning. However, deflections in this case remain weak even at small distance from the centre. The reason is that the Reynolds numbers investigated are low (Re 20) and the wakes are thin and purely laminar. From such data, it not easy to estimate the importance of a radial pressure gradient on the wake deflection, when the wake is separated and the bubble is spinning.

The far-wake deflection and the bubble spinning seem to delay the transition at which vortex shedding generally occurs for a stationary sphere in uniform flow (Re > 274). This appears in figure 9(d) where Ω * is significant (≈ 0.11) and where, for Re ≈ 300 and Ro ≈ 7.4, no clear vortex shedding signature can be detected. This trend is to some extent reminiscent of the findings of [START_REF] Giacobello | Wake structure of a transversely rotating sphere at moderate Reynolds numbers[END_REF] and [START_REF] Kim | Laminar flow past a sphere rotating in the transverse direction[END_REF], who showed that the rotation of a sphere may trigger or suppress vortex shedding depending on the Re and Ω * considered. For large Re (figure 12), the low-velocity side of the wake begins to oscillate. This behaviour has some common features with the wake instability of Kelvin-Helmholtz type noted by Giacobello drastically increases, and we can reasonably think that it destabilizes the layer. In comparison, the high-velocity side where the shear is much smaller, since the surface turns with the flow, looks constantly more steady. These high Re numbers correspond to a range where bubbles are slightly deformed (aspect ratio ≈ 1.1-1.2) and undergo small-amplitude circular precession motions around their equilibrium position, due to the wake instability. The near wake is then reducing towards disappearing. As it reduces to a thin layer of fluid behind the bubble, the precession induces a periodic pinching of the near wake that splits it into two rather symmetrical parts attached to each side of the wake. A small recirculation can be observed inside the two regions of the near wake thus formed (figures 9(d) and 12). For this very specific situation, FIGURE 12. Large bubble: Re = 560. The high-velocity side of the wake remains steady (lower yellow line). The low-velocity side of the wake exhibits an oscillation (upper yellow line). The near-wake layer that is thin for these high Re numbers is here pinched due to the precession of bubble motion, leaving only two small recirculation regions (dotted yellow line). it is not possible to estimate a chord length l NW . The surface of this split near wake S NW (blue symbols in the near-wake figures) can still be defined but takes on a special meaning.

The next subsection presents the characteristics of the far and near wake as a function of Re, Ro and Ω * .

3.2. Far-wake characteristics 3.2.1. Far-wake shape and geometry

The variation of the curvilinear length of the far wake (s FW -d/2) normalized by d is given in figure 13 as a function of Re. As expected, this normalized length increases with meaning that the region behind the bubble that is affected by the wake is increasing. A least-squares (LS) fit of the data (including all Ro) yields an increase as Re 1/4 . When maximum, the curvilinear length of the far wake is approximately equal to 0.3 times the full rotation distance 2πr e at the bubble equilibrium location. We deduce that the bubbles in these runs never feel their own wake as if it were the wake of an upstream bubble. This is all the more true given that the wake tends to deflect towards the centre. The data scattering in figure 13 is amplified by the difficulty in estimating r e accurately, as has already been pointed out. The few data available with fully contaminated bubbles (triangles) has been added for comparison. We see that their far-wake length follows the same general trend as the partially contaminated ones, although their spinning is uniform like spheres [START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF].

The mean radius of curvature of the wake r FW was estimated by approximating the curvilinear wake length s FW by an arc of a circle. e versus Re number. Far wake: ×, Ro < 6; E , 6 Ro 8.5; @ , Ro > 8.5. Streamlines: +, Ro < 6; u , 6 Ro 8.5; 6 , Ro > 8.5. Solid line, κ + e = (125/3)Re -2/5 . and the reciprocal of the value thus calculated gives the mean far-wake curvature κ FW . This curvature normalized by the curvature of the tank is denoted by κ + FW . Its Ro and Re dependence is displayed in figure 14, together with the non-dimensional curvature of the flow streamline at the bubble equilibrium position κ + e (κ e = 1/r e ). An LS fit of all the data shows that, on average, κ e decreases as Re -2/5 .

As expected, for a given Re, the curvature κ + e is greater for smaller Ro, the reason being that bubbles located closer to the centre regarding their diameter see an incoming rotational flow that is more curved. This is particularly visible for Re in the range [100,200]. Whatever Re in that range, κ + FW is clearly higher than κ + e when Ro 8.5, the difference increasing as Ro decreases. The far wake is therefore more curved than the incoming flow streamline and this trend amplifies as the relative bubble distance from the centre decreases. This is what can be noted when comparing figures 9(b) and 10(d). In contrast, no significant difference in curvature is on average .45 . Fully contaminated bubbles, q . (b) Non-dimensional far-wake area versus the non-dimensional area that would be obtained with an axisymmetric straight far wake: ×, Re < 150; E , 150 Re 250; @ , Re > 250. Dashed line, y = x; solid line, 0.3x 1.6 . Fully contaminated bubbles, q .

observed for Ro > 8.5. At intermediate distances (6 Ro 8.5), the two curvatures are very close to each other up to Re = 150, while above, κ + FW is systematically greater than κ + e . The data for the fully contaminated bubbles exhibit the same trends and so have not been reported.

The far-wake surface normalized by the bubble surface S B is plotted in figure 15(a). An LS fit of all the data shows that it evolves as Re 0.45 , thus increasing with Re. The influence of the curvature on the wake area S FW was investigated by comparing S FW with a reference area S FWref defined as the product of the bubble diameter and the curvilinear length. It corresponds to the surface that would be obtained with an axisymmetric straight far wake. The results are shown in figure 15(b) in non-dimensional form, with all areas normalized by the bubble surface S B .

For Re up to 150-200, the exact value depending on the Ro number, S FW is slightly smaller than S FWref , while above it becomes rather rapidly higher. Referring to figure 14, the low-Reynolds-number behaviour corresponds to situations where both curvatures κ + FW and κ + e are important. In that case, the difference of velocity seen by both sides of the bubble together with the curvature of the incoming flow streamlines reduces the area of the far wake, which can explain that S FW is slightly lower than S FWref . In contrast, for high Re (400-500), the flow and the bubble wake are weakly curved and there is no obvious link between this smaller curvature and the fact that S FW > S FWref . A more plausible explanation in that case seems to be that, as the far wake grows in length (figure 13), its cross-section also thickens (figure 10), making S FW > S FWref . The cross-section growing can be attributed to various conjugated effects. The spinning of the surface is close to its maximum (3.1), the bubble deformation is no longer negligible (aspect ratio 1.2) and, as explained before, small-amplitude circular precessing motions appear. It follows that the wake cross-section is larger than the mean equivalent diameter of the bubble d and explains the trends observed.

Far-wake deflection

The wake deflection is quantified by the angle γ 2 between the chord l FW and the vertical line passing through the centre of the bubble (figure 5a). Clearly, this angle Partially contaminated bubbles: ×, Ro < 6; E , 6 Ro 8.5; @ , Ro > 8.5. Fully contaminated bubbles, q . Solid line, γ 2 = 5.5Re 0.45 ; dotted line, 'mean deflection angle' γ 2 -γ 1 , same equation as γ 2 but with factor 5.5/2.

increases with Re (figure 16) up to 90 • for Re numbers of the order of 400-500. This physically occurs when the far wake extends far downstream, until approximately 180 • away from the bubble around the axis of rotation, impacting a large region of the tank. The trends are very similar to those exhibited by the curvature. For Re in the range [100, 200], deflection is higher when Ro 8.5, the difference increasing as Ro decreases, that is, as the bubble is closer to the centre. At intermediate distances (6 Ro 8.5), no Ro significant effect is observed, the smaller Re are less deflected, while the large ones are more deflected (see examples in figure 9a,d). The behaviour of the 'fully contaminated data' is rather similar. An LS fit of all the data shows that the deflection angle on average varies as Re 0.45 . To go further in the analysis, figure 17 compares the angle γ 2 with the angle γ 1 that would be measured if the wake had followed the incoming flow streamlines at the bubble equilibrium position. The interesting result is that γ 2 is, on a first-order basis, twice the angle γ 1 that would exist in this case. The far wake therefore always deflects from the incoming flow streamlines with a mean 'deflection angle' γ 2 -γ 1 that also varies as Re 0.45 , but with half the factor of γ 2 . This mean angle is represented by the dotted line in figure 16. As discussed in § 3.1, it is difficult to infer from the measurements which effect, between the flow pressure gradient and the bubble spinning, contributes more to the far-wake deflection. We may conjecture that the deflection is driven by both bubble spinning and pressure gradient at small distances behind the bubble and that pressure gradient (centripetal effect) has more influence at larger distance.

Near-wake characteristics

We focus in this section on the region located close behind the bubble where the flow structure is directly impacted by the surface spinning. As a rule, knowledge of the near wake is important for the forces acting on the bubble and for other aspects like the heat and mass transfers [START_REF] Clift | The rise velocity and shape of bubbles in pure water at high Reynolds number[END_REF][START_REF] Fan | Bubble Wake Dynamics in Liquids and Liquid-Solid Suspensions[END_REF]. For a steady sphere in a uniform flow, this region is well described in the literature. [START_REF] Taneda | Experimental investigation of the wake behind a sphere at low Reynolds[END_REF] was among the first to visualize and characterize the recirculations that structure this near wake. Later on, [START_REF] Johnson | Flow past a sphere up to a Reynolds number of 300[END_REF] provided a full description of the structure of this near wake at Re up to 300, both experimentally and numerically. More recently, [START_REF] Kim | Laminar flow past a sphere rotating in the transverse direction[END_REF], [START_REF] Giacobello | Wake structure of a transversely rotating sphere at moderate Reynolds numbers[END_REF] and [START_REF] Poon | Laminar flow structures from a rotating sphere: effect of rotating axis angle[END_REF][START_REF] Poon | Hydrodynamic forces on a rotating sphere[END_REF] numerically studied this near wake when the sphere is rotating transversely and identified various behaviours as a function of Ω * and Re. In our experiment, this near-wake region can be rather well identified on the instantaneous images as a region where the velocity of seeding particles is lower than in the main flow stream ( § 2.2).

Near-wake shape and geometry

The geometrical parameters that characterize the near wake are its chord length (l NW -d/2) and its surface S NW (figure 5b). Figure 18 compares the Re dependence of this chord length normalized by the bubble diameter with that of the non-dimensional near-wake length L NW(uniform) behind a steady sphere in a uniform flow, as simulated by [START_REF] Giacobello | Wake structure of a transversely rotating sphere at moderate Reynolds numbers[END_REF] and [START_REF] Tomboulides | Numerical investigation of transition and weak turbulent flow past a sphere[END_REF], and measured by [START_REF] Nakamura | Steady wake behind a sphere[END_REF]. It shows that the dissymmetry and the deflection of the near wake that are induced by the bubble spinning and the flow rotation result in a chord length that increases with Re, but is systematically larger than L NW(uniform) . This holds for Re up to 150, where both recover the same order of magnitude. This corresponds to the range of Re where bubble spinning has reached its maximum. Regarding the surface (figure 19), we see that S NW is smaller than the surface S NWref that would be obtained for the same bubble diameter with an axisymmetric straight near wake. The difference between the two increases with Re, that is, when bubble spinning and deflection are much more important. By curving the high-velocity-side streamline more or less towards the low-velocity-side streamline, these two effects control the point where they meet and the boundaries of the near wake. This gives this typical shape in figures 5(b) and 10(d) with a sharp tip at the point where they meet and an external streamline more or less curved. The trends are globally the same for a fully contaminated interface. . Non-dimensional near-wake chord length versus Re. Partially contaminated bubbles, ×, Ro < 6; E , 6 Ro 8.5; @ , Ro > 8.5. Fully contaminated bubbles, q . Dashed line, [START_REF] Giacobello | Wake structure of a transversely rotating sphere at moderate Reynolds numbers[END_REF] and [START_REF] Tomboulides | Numerical investigation of transition and weak turbulent flow past a sphere[END_REF]; dash-dotted line, [START_REF] Nakamura | Steady wake behind a sphere[END_REF].

Near-wake surface dependences

Since the near wake is the region behind the bubble that is a priori the most directly influenced by the surface spinning, we have plotted its surface evolution as a function of Ω * . This yields figure 20, where two asymptotic behaviours appear. For moderate Re in the range [20,140] where Ω * (Re) is given by the equation in figure 7(a). This highlights that bubble spinning strongly influences the flow in the near wake for Re ∈ [20,140]. When normalized by the bubble surface, this dependence becomes The impact of the bubble interface spinning on the surface S NW for Re ∈ [20, 140] can explain why S NW is higher than what can be observed for inclusions without an interface rotation for the same Re and d. This is illustrated in figure 21, where S NW /S B is plotted as a function of Re, together with the normalized near-wake surface for spheres in a uniform flow. The normalized surface for non-rotating spheres in a uniform flow was estimated in two ways that agree with one another: from the correlation proposed by [START_REF] Fan | Bubble Wake Dynamics in Liquids and Liquid-Solid Suspensions[END_REF], and from the separation angles and the wake length calculated by [START_REF] Giacobello | Wake structure of a transversely rotating sphere at moderate Reynolds numbers[END_REF]. Up to 140, S NW /S B is clearly higher than the correlations for non-rotating spheres in uniform flows, while it is well reproduced by (3.9). In contrast, from 140 up to 240, the trend changes and S NW /S B then follows this correlation rather well, even above its limit of validity (Re = 200). The increase of S NW /S B is no longer linked to the surface spinning, which has reached its maximum, but to the increase of the flow velocity U = ωr e seen by the bubble, with, as a result, a near wake that extends farther downstream (figure 18).

The behaviour at high Re (Re 240), where S NW /S B suddenly decreases and saturates, corresponds to situations where the wake is unstable. As pointed out in § 3.1, this instability manifests by the occurrence of a small-amplitude circular precessing bubble motion that induces a periodic pinching of the near wake. The pinching splits ×, Ro < 6; black E , 6 Ro 8.5; @ , Ro > 8.5; filled E , split near wake. Fully contaminated bubbles: q . Thick line, (S NW /S B )(20 Re 140) = 25(exp(-0.014 × 20)exp(-0.014Re))(Re/20) -1.10 ; solid line, (S NW /S B )(Re 240) = 2(Re/240) -5/2 ; dashed line, surface variation for a bubble in a uniform flow (Re -20) 1.12 /200 (20 Re 200), [START_REF] Fan | Bubble Wake Dynamics in Liquids and Liquid-Solid Suspensions[END_REF], extrapolated up to Re = 250; dash-dotted line, [START_REF] Giacobello | Wake structure of a transversely rotating sphere at moderate Reynolds numbers[END_REF] and [START_REF] Tomboulides | Numerical investigation of transition and weak turbulent flow past a sphere[END_REF].

the near wake into two almost symmetrical parts and contributes to decrease its area. Also, the far wake extends far downstream and is strongly deflected. The angle γ 2 -γ 1 is typically of the order of 30 • to 45 • (figure 16) and the high-velocity-side streamline is deflected inwards so much (figure 9d) that it prevents the downstream development of the near wake. These two effects (instability plus deflection) explain why the near wake then reduces to a very thin layer of fluid behind the bubble that progressively vanishes at Re above 300. The surface of this layer normalized by the bubble surface clearly decreases with Re as (figure 21)

S NW S B (Re 240) = 2 Re 240 -5/2 . (3.10)
These changes in the behaviour of the near wake coincide with the changes in the lift coefficients noted by [START_REF] Rastello | Drag and lift forces on interfacecontaminated bubbles spinning in a rotating flow[END_REF][START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF] between Re < 200 and Re > 250. It was shown that the total lift coefficient C L was the sum of the lift coefficient C L (SBRH) induced by the solid-body rotating flow and an excess lift C L (BS) induced by the bubble spinning. Subtracting the excess lift from the total lift yields the lift contribution of the solid-body rotation flow (see figure 20 in [START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF]). This lift contribution is close to that obtained by [START_REF] Bluemink | A sphere in a uniformly rotating or shearing flow[END_REF][START_REF] Bluemink | Drag and lift forces on particles in a rotating flow[END_REF] for non-spinning spheres and behaves differently depending on Re. It increases for Re up to approximately 150 while it decreases with Re above 250. According to the latter authors, these different regimes can be attributed to changes in the structure of the wake behind the sphere. The data presented here bring experimental evidence of these changes through the near-wake surface evolution and confirm this statement. They suggest that vorticity around the bubble is changing with Re and so is the lift. Coefficient C L (SBRH) was also found to increase with Ro, but this effect is more difficult to identify. The only visible trend in figure 21 is that, for 150 Re 200, the near-wake surface is on average smaller for Ro < 6 (×) than for large Ro. This is precisely the Re range where the flow-induced lift is most influenced by Ro.

Conclusion

This work provides new results that complete our previous studies on the drag and lift forces on a bubble immersed in a solid-body rotating flow around a horizontal axis [START_REF] Rastello | Drag and lift forces on interfacecontaminated bubbles spinning in a rotating flow[END_REF][START_REF] Rastello | Drag and lift forces on clean spherical and ellipsoidal bubbles in a solid-body rotating flow[END_REF][START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF]. We focus on the situation where the bubble surface is contaminated and the bubble is freely spinning, thus creating an excess lift force in comparison with a clean surface. The objective was to describe the complex separated wake that forms behind the bubble in that case. The two-dimensional geometry and behaviour of this wake have been inferred from laser-sheet visualizations with fluorescent seeding microparticles. They have been analysed as a function of the local parameters that characterize the rotational flow at the point where the bubble is at equilibrium (Re and Ro numbers), and the non-dimensional spinning rate Ω * that quantifies the importance of the bubble spinning velocity relative to the upstream flow velocity.

We find that the flow rotation and the bubble spinning cause a significant dissymmetry of the wake. Under the influence of the surface spinning, the highvelocity flow arriving on one side of the bubble turns around the bubble rear and moves upwards on the other side till it meets the downward low-velocity flow, forming there a high-shear layer. This causes the stagnation point observed behind the bubbles in the clean situation to displace upwards on the low-velocity side of the bubble. This effect increases as Re and thus Ω * increase.

With increasing Re, the wake extends farther downstream and deflects towards the flow rotation axis. The angle of deflection relative to the incoming streamline at the bubble equilibrium position increases with Re concomitantly with the length of the far wake. It is typically of the order of 50 • at Re = 450. For a given Re, this angle is larger for bubbles located close to the centre of the flow (small Ro) reflecting a stronger deviation. The far wake for these bubbles is also more curved than the incoming streamlines. The wake deflection can be attributed to two main effects: the bubble spinning, which generates the wake dissymmetry close to the surface, and the centripetal action of the pressure gradient. However, it is difficult to evaluate which of these effects contributes more to the deflection.

The study of the near wake has yielded the following conclusions.

(i) For Re 140, the surface spinning strongly influences the structure of the near wake, so that its size increases faster with Re than for a stationary sphere in a uniform flow. An expression giving the normalized surface of the near wake as a function of Ω * and Re has been proposed. (ii) For 140 < Re 240, the influence of the bubble spinning reaches its maximum and is no longer dominant. The size of the near wake is then linked to the increase of the free-stream velocity U = ωr e with Re and comparable with that of a stationary sphere in a uniform flow. (iii) For Re > 240, the near wake decreases rapidly with Re until it disappears for higher Re. This decrease is caused by the wake instability and the small circular precession motions of the bubble that result. The precession contributes to cyclically split the near wake whose area vanishes at high Re into two small recirculating zones. 240) turns to be the range where the coefficient of the lift induced on the bubble by the solid-body rotation flow [START_REF] Rastello | Drag and lift forces on interfacecontaminated bubbles spinning in a rotating flow[END_REF][START_REF] Rastello | Clean versus contaminated bubbles in a solid-body rotating flow[END_REF] stops increasing with Re and reaches a maximum before decreasing at higher Re values. We can therefore reasonably conclude that the changes of this lift coefficient with Re are closely linked to the changes in the near-wake development around the spinning bubbles and thus to the changes in vorticity near the surface that its implies. that can rotate around a horizontal axis thanks to ball bearings (3). The ball bearings are installed in two stainless-steel rigid disks (4), fixed together by three rods (5) to ensure a good rigidity of the system. The cylindrical counterbore is entrained in rotation by a toothed belt (6) fixed at its back and coupled to an 80 W Crouzey brushless electric motor (7), through a transmission system (8). The ensemble of the device is screwed onto a thick iron casing (9), equipped with four adjustable legs (10), whose tips are in a damping material to remove vibrations.

FIGURE 1 .

 1 FIGURE 1. Two examples of wake imaging behind a bubble at equilibrium location (Rastello et al. 2009; Rastello, Marié & Lance 2017): (a) clean surface (no spinning), Re ≈ 15; and (b) contaminated surface (spinning), Re ≈ 120.

FIGURE 2 .

 2 FIGURE 2. Experimental optical set-up. (a) Force measurements. The position, size and shape of the bubble are determined from images recorded with a backlight on a fast camera located in front of the tank. (b) Surrounding flow investigation. The surface spinning and the wake behind the bubble at its equilibrium position are visualized with a laser sheet as shown. With this lighting, a shadow appears behind the bubble on the bubble side that is opposed to the laser source. The fast camera (not visible in the figure) is still positioned in front of the tank.

FIGURE 4 .

 4 FIGURE 4. Bubble wake for Re = 100, Ro = 7.7, Ω * = 0.09, equation (3.3): (a) streamlines obtained by averaging 10 successive frames; and (b) streamlines from averaging 20 PIV velocity vector fields.

FIGURE 5 .

 5 FIGURE 5. Schematic of the wake region and its parameters. (a) Far wake (FW): l FW = chord length, s FW = curvilinear length, S FW = area, γ 2 = deflection angle, γ 1 = streamline angle with the same l FW . (b) Near wake (NW): l NW = chord length, S NW = area.

FIGURE 8 .FIGURE 9 .

 89 FIGURE 8. The Re dependence of d obtained in two ways. (1) By measuring d and r e directly from the images taken with the backlight (force measurements): E , Ro < 10; ×, Ro 10.(2) With d measured from the laser sheet (bubble rotation and wake measurements) and Re estimated by solving (3.5): @ , Ro < 10; +, Ro 10; solid line, 0.0006 + 3.85 × 10 -6 Re; dotted line, 9.5 × 10 -5 Re 0.5 .

FIGURE 10 .

 10 FIGURE 10. Bubble wake for increasing Re and Ro ≈ 11: (a) Re = 34, Ro = 11.9, Ω * = 0.05; (b) Re = 47, Ro = 10.8, Ω * = 0.07; (c) Re = 68, Ro = 10.8, Ω * = 0.06; (d) Re = 109, Ro = 10.7, Ω * = 0.07.

  (a) with figure 10(b), and figure 9(b) with figure

FIGURE 11 .

 11 FIGURE 11. Wake deflection as a function of Ro at comparable Re: (a) Re = 49, Ro ≈ 7.3; (b) Re = 47, Ro ≈ 11; (c) Re = 102, Ro ≈ 7.3; (d) Re = 109, Ro ≈ 11.

  et al. (2009) and Poon et al. (2014) on a rotating sphere at high Re and Ω * . The maximum rotation rate Ω * being locally very high, the shear in the layer developing on the low-velocity side of the bubble where the surface turns against the flow (the advancing side of the sphere in Giacobello et al. (2009) and Poon et al. (2014))

  FIGURE 15. (a) Non-dimensional far-wake area versus Re number. Partially contaminated bubbles: ×, Ro < 6; E , 6 Ro 8.5; @ , Ro > 8.5. Solid line, S FW /S B = 1.1Re 0.45 . Fully contaminated bubbles, q . (b) Non-dimensional far-wake area versus the non-dimensional area that would be obtained with an axisymmetric straight far wake: ×, Re < 150; E , 150 Re 250; @ , Re > 250. Dashed line, y = x; solid line, 0.3x 1.6 . Fully contaminated bubbles, q .

FIGURE 16 .

 16 FIGURE 16. Deflection angle γ 2 as a function of Re. Partially contaminated bubbles: ×, Ro < 6; E , 6 Ro 8.5; @ , Ro > 8.5. Fully contaminated bubbles, q . Solid line, γ 2 = 5.5Re 0.45 ; dotted line, 'mean deflection angle' γ 2 -γ 1 , same equation as γ 2 but with factor 5.5/2.

FIGURE 17 .

 17 FIGURE 17. Deflection angle γ 2 in degrees as a function of γ 1 in degrees. Partially contaminated bubbles: ×, Ro < 6; E , 6 Ro 8.5; @ , Ro > 8.5. Fully contaminated bubbles, q . Dashed line, y = x; solid line, y = 2x.

  FIGURE 18. Non-dimensional near-wake chord length versus Re. Partially contaminated bubbles, ×, Ro < 6; E , 6 Ro 8.5; @ , Ro > 8.5. Fully contaminated bubbles, q . Dashed line,[START_REF] Giacobello | Wake structure of a transversely rotating sphere at moderate Reynolds numbers[END_REF] and[START_REF] Tomboulides | Numerical investigation of transition and weak turbulent flow past a sphere[END_REF]; dash-dotted line,[START_REF] Nakamura | Steady wake behind a sphere[END_REF].

FIGURE 19 .

 19 FIGURE 19. Normalized near-wake area versus the normalized one that would be obtained with an axisymmetric straight near wake. Partially contaminated bubbles: ×, Re < 150; E , 150 Re 250; @ , Re > 250. Fully contaminated bubbles, q . Dashed line, y = x; solid line, y = 0.7x1.35 .

  Re 140) = 25(exp(-0.014 × 20) -exp(-0.014Re)) above Re = 150, S NW / SB and Ω * (Re) are uncorrelated. While Ω * (Re) saturates because Ωb /ω reaches a constant value, S NW / SB still increases before decreasing at high Re.

FIGURE 21 .

 21 FIGURE 21. Normalized near-wake surface versus Re. Partially contaminated bubbles: ×, Ro < 6; black E , 6 Ro 8.5; @ , Ro > 8.5; filled E , split near wake. Fully contaminated bubbles: q . Thick line, (S NW /S B )(20 Re 140) = 25(exp(-0.014 × 20)exp(-0.014Re))(Re/20) -1.10 ; solid line, (S NW /S B )(Re 240) = 2(Re/240) -5/2 ; dashed line, surface variation for a bubble in a uniform flow (Re -20) 1.12 /200 (20 Re 200),[START_REF] Fan | Bubble Wake Dynamics in Liquids and Liquid-Solid Suspensions[END_REF], extrapolated up to Re = 250; dash-dotted line,[START_REF] Giacobello | Wake structure of a transversely rotating sphere at moderate Reynolds numbers[END_REF] and[START_REF] Tomboulides | Numerical investigation of transition and weak turbulent flow past a sphere[END_REF].

  FIGURE 22. Three-dimensional general view of the device.

FIGURE 23 .

 23 FIGURE 23. Lateral view of the device.

TABLE 2 .

 2 Statistical uncertainties.

	σ SFW	σ lFW	σ SNW	σ lNW
	S B	d	S B	d
	0.04	0.05	0.04	0.07

Downloaded from https://www.cambridge.org/core. SID UGA INP, on 10 Dec 2019 at 11:49:24, subject to the Cambridge Core terms of use, available at https://www.cambridge.org/core/terms. https://doi.org/10.1017/jfm.2019.936

  exp(-0.014Re)). Fully contaminated bubbles, q . The dotted lines represent 0.05Ω * increments. (b) Pressure gradient at bubble location r e ; same symbols.

	https://doi.org/10.1017/jfm.2019.936 https://www.cambridge.org/core/terms.	0.25 0.20 0.15 0.10 0.05 0 FIGURE 7. (a) Non-dimensional rotation rate Ω * as a function of Re. Partially conta-) 6 5 (a) (b) Ø * 100 200 300 400 4 3 ™P/™r (kPa m -1 500 Re 0 100 200 300 400 500 Re 2 1 minated bubbles: ×, Ro < 6; E , 6 Ro = 0.12 8.5; @ , Ro > 8.5; solid line, Ω * (1
	at
	available
	use,
	of
	terms
	Core
	Cambridge
	the
	to
	subject
	11:49:24,
	at
	2019
	Dec
	10
	on
	INP,
	UGA
	SID
	https://www.cambridge.org/core.
	from
	Downloaded

  Then r FW is the solution of FIGURE 13. Normalized far-wake curvilinear chord length versus Re. Partially contaminated bubbles: ×, Ro < 6; E , 6 Ro 8.5; @ , Ro > 8.5. Fully contaminated bubbles, q . Solid line, (s FW -d/2)/d = 2Re 1/4 .

			20
			15
	https://www.cambridge.org/core/terms.	(s FW -d/2)/d	0 5 10	100	200	Re	300	400	500
	at
	available		20
	use,
	of
	terms		15
	Cambridge Core	k + FW , k +	10
	the		5
	to
	subject
	11:49:24,		0	100	200	Re	300	400	500
	at
	2019
	Dec
	10
	on
	INP,
	UGA
	SID
	https://www.cambridge.org/core.
	Downloaded from		l FW -2r FW sin	s FW 2r FW	= 0,		(3.7)

https://doi.org/10.1017/jfm.2019.936 e FIGURE 14. Non-dimensional curvature of the far wake κ + FW and of the streamlines at bubble equilibrium position κ +

  , we note that S NW normalized by the average bubble surface SB increases quasi-proportionally to the non-dimensional rotation rate as S FIGURE 20. Normalized near-wake surface versus non-dimensional rotation rate Ω * . Partially contaminated bubbles: ×, Ro < 6; black E , 6 Ro 8.5; @ , Ro > 8.5; filled E , split near wake. Fully contaminated bubbles, q . Solid line: (S NW / SB )(20 Re 140) = 31(Ω * (Re) -Ω * (Re = 20)).

	https://www.cambridge.org/core/terms. https://doi.org/10.1017/jfm.2019.936	S NW /S	0 2 1 6 5 4 3	0.05	Ø *	0.10	0.15
	at
	available
	use,
	of
	terms
	Core
	Cambridge
	the
	to
	subject
	11:49:24,
	at
	2019
	Dec
	10
	on
	INP,
	UGA
	SID
	https://www.cambridge.org/core.
	from
	Downloaded

NW SB (20 Re 140) = 31(Ω * (Re) -Ω * (Re = 20)), (3.8) B

Acknowledgements

The authors would like to thank Nathalie Grosjean for her technical assistance in the attempt to process the visualizations using the LaVision PIV software.