Characterization of coupling between inertial particles and turbulent wakes from porous disk generators

<jats:p>This study presents the findings of a wind tunnel experiment investigating the behaviour of micrometric inertial particles with Stokes numbers around unity in the turbulent wake of a stationary porous disk. Various concentrations <jats:inline-formula>        <jats:alternatives>   <jats:inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="png" xlink:href="S0022112021010958_inline2.png" />   <jats:tex-math>$\varPhi _{v}\in ([6-19] \times 10^{-6})$</jats:tex-math>        </jats:alternatives>      </jats:inline-formula> of poly-disperse water droplets (average diameter 40–50 <jats:inline-formula>        <jats:alternatives>   <jats:inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="png" xlink:href="S0022112021010958_inline3.png" />   <jats:tex-math>$\mathrm {\mu }$</jats:tex-math>        </jats:alternatives>      </jats:inline-formula>m) are compared with sub-inertial tracer particles. Hot-wire anemometry, phase Doppler interferometry and particle image velocimetry were implemented in the near- and far-wake regions to study the complex dynamics of such particles. Quadrant analysis is used to explore the shear effects of the particle wake interaction. Turbulence statistics and particle size distributions reveal distinct differences in the structure of the wake when inertial particles are present in the flow. Additionally, there are different structures in the near and far wake regions and structures change with particle volume fraction.</jats:p>

section_content: Many flows of interest in the natural environment and engineering applications combine turbulent wakes and inertial particles. Turbulent wakes and inertial particles are of relevance where re-circulation, entrainment, energy deficits and pressure fluctuations affect the design and efficiency of wake-creating bodies such as wind turbines, aircraft, buildings and bridges (Cheynet et al. 2017; Frej Vittale & Davenport 2017; Ali et al. 2018; Hertwig et al. 2019; de Jong Helvig et al. 2021; Smith et al. 2021). Cloud formation, industrial spray applications, combustion and air pollution are also of interest (Chandrakar et al. 2017; Stodt, Kiefer & Fritsching 2019; Wei et al. 2020; Gai et al. 2021). This study explores the effects of inertial particles in the wake behind a stationary porous disk in homogeneous isotropic turbulence (HIT). 

For some time, the wakes behind porous disks have been studied as analogues to more complicated wakes such as parachutes and rotating wind turbines in single-phase flow applications. Roberts (1980) studied drag coefficients of solid and porous disks with open area ratios of 2-33 per cent for parachute applications, and Aubrun et al. (2013) compared wake properties of a porous disk with a rotating wind turbine model in wind tunnel experiments. Lignarolo et al. (2016) presented an experimental analysis comparing the near wakes of a wind turbine model and porous disk. Their results establish a good match between the turbine and the disk for velocity, pressure and enthalpy fields, but show differences in turbulence intensity and turbulent mixing. Aloui et al. (2013) compared particle image velocimetry measurements of steady and unsteady wakes behind porous disks. Their proper orthogonal decomposition analysis showed that alternating vortices form in the unsteady wake. 

In the case of wind energy, various porous disk designs have been studied in single-phase flows to approximate wind turbine wakes, as stationary disks are easier to simulate than rotating blades. (Martinez et al. 2012; Naderi & Torabi 2017) A measurement campaign by Aubrun et al. (2019) compared uniform with non-uniform porous disks in nine different wind tunnels to investigate mean streamwise velocity and turbulence intensity profiles four diameters downstream of the disk. They found the results collapse reasonably well across facilities and that the non-uniform disk created a smaller velocity deficit and greater turbulence intensity in the centre of the wake. Camp & Cal (2016) developed a non-uniform porous disk, and compared the disk array with an array of rotating turbines and quantified the differences in the mean kinetic energy transport within the wakes. They ascertained the primary difference was in the spanwise mean velocity component in the near-wake region. 

Others have studied various aspects of inertial particles in turbulent flows without wakes. Specifically, heavy particles at Stokes numbers of the order of unity, where a maximum coupling between particles and turbulence occurs, as described in Toschi & Bodenschatz (2009). The mechanisms behind particle settling velocity modification in HIT are still being explored, and different studies have found both and increase and decrease in particle settling velocity, depending on various parameters and experimental and numerical conditions. For example, Aliseda et al. (2002) found the particle settling velocity to be larger than the quiescent fluid (an effect generally associated with the preferential sweeping mechanisms (Maxey 1987 ) and settling velocity enhancement increases with volume fraction, while Good et al. (2014) found that nonlinear drag could induce settling hindering (associated with the loitering effect Nielsen 1993) where regions of high velocity root-mean-square (r. m. s. ) anisotropy generally coincide with regions of settling velocity reductions. Numerical simulations by Rosa et al. (2016) found, however, that hindering by loitering rarely occurs unless particle horizontal motion is artificially blocked. 

Settling velocity is also related to particle preferential concentration or clustering in turbulent flow structures. Sumbekova et al. (2016) and Obligado et al. (2019) studied settling velocity and preferential concentration of inertial particles in HIT. Sumbekova et al. (2016) determined particles within clusters settle faster than particles in voids. This 933 A42-2 confirms observations by Aliseda et al. (2002). The same observation was also made by Huck et al. (2018), with a model to interpret the enhancement of settling as a collective 'effective rheology' effect. Obligado et al. (2019) found cluster settling velocity has a strong dependence on the Taylor scale Reynolds number Re λ or turbulence acceleration ratio, which is the ratio of the standard deviation of fluid acceleration fluctuations to gravity. Both Sumbekova et al. (2017) and Obligado et al. (2019) show that clustering increases with volume fraction (Φ v ) and Re λ, and the mean size of clusters increases with Re λ. However, Obligado et al. (2019) (working at a fixed value of the free-stream velocity U ∞ ) finds cluster size decreases with particle volume fraction Φ v, while Sumbekova et al. (2017) found cluster size increases with Φ v. Gustavsson, Vajedi & Mehlig (2014) studied clustering of particles falling in turbulent flow and found the clustering is dependent on Stokes number and may be very anisotropic. Wang, Gu & Zheng (2020) compared large and very large scale motions of single-phase flows with sand-laden flows with particles less than 10 μm and particle volume fractions < 10 -6. They found the streamwise turbulent kinetic energy is enhanced across all scales and that particles can affect turbulence at large scale structures. 

Fewer studies have explored inertial particles in turbulent wakes. For instance, Fessler & Eaton (1999) found turbulence modification in a particle-laden flow over a backward facing step was a function of particle Stokes number and Reynolds number. More recently, Homann & Bec (2015) performed direct numerical simulation of inertial particles in the wake of a sphere and found particles form preferential concentrations in the wake. They found a region downstream of the sphere where particle concentration exceeds the inflow and suggest this is due to ejections of particles from detached and advected vortices. 

This study aims to characterize the behaviour of inertial particles in the axisymmetric turbulent wake of a porous disk. This is a first step toward answering the question if and how stationery disks can be used as analogues to rotating turbines in the presence of inertial particles. As this is a first study on the subject, the focus is on the HIT incoming flow case. The persistence of the wake, how particles are entrained, particle discrimination by size, and particle settling velocity are studied. To the authors' knowledge, this is the first experimental study on the coupling of inertial particles with a self-similar, large Re λ flow behind a porous disk generator. The experimental set-up and data collection techniques are presented in § 2, and results are presented in § 3. The application of the analysis techniques and discussion follow in § 4, and concluding remarks on the implications are given in § 5.
section_name: Introduction

section_content: Experiments were conducted at the Université Grenoble Alpes at Laboratoire des Ecoulements Géophysiques et Industriels (LEGI) in the Lespinard Wind Tunnel. The closed-circuit wind tunnel test section is 4 m long with a cross-sectional area of 0. 75 m × 0. 75 m as shown in figure 3. Figure 1 shows the passive spray grid with 36 water misting nozzles (0. 4 mm diameter) at the tunnel inlet that inject water droplets into the tunnel flow. A stationary open grid is located 15 cm upstream of the spray grid to override and mix the turbulence added by the injectors, and the combined grids produce turbulence intensity of the order of 2. 4 %-3. 0 %. A 150 bar pump supplied water to the spray grid, producing a uniform spray of poly-disperse water droplets with most probable diameters between 40-50 μm. Water volumetric flow rate was controlled with a manual variable regulator and data were collected after the water flow rate, humidity and tunnel velocity reached steady state. 

A non-uniform porous disk attached to a 12. 7 mm diameter aluminium tube was mounted in the wind tunnel and attached at the tunnel floor as shown in figures 3 and 1. The disk has a diameter of D = 120 mm, a thickness of 3. 175 mm, a porosity of approximately 56 % and a tunnel blockage ratio of 1. 57 % including the mounting tube. A small blockage ratio allows unimpeded expansion of the wakes within the tunnel. Holes in the disk range from 0. 05-8. 5 mm in size (refer to figure 2 for disk geometry). The disk was designed by Camp & Cal (2016) and the porosity varies radially to mimic the design of a wind turbine rotor. 

Hot-wire anemometry (HWA), phase Doppler interferometry (PDI) and particle image velocimetry (PIV) measurements were taken of the background flow (i. e. a passive grid flow with no wake) and compared with locations at one disk diameter (1D) and 9. 6 disk diameters (9. 6D) downstream of the disk to analyse characteristics of different flow regimes in the wake. The disk was positioned such that the measurement location Figure 3. Schematic of wind tunnel experimental set-up, side view. (Schematic is not to scale. ) Note that the measurement location is fixed at 3 m downstream of the grid and the disk model was moved for near-and far-wake measurements. Grey boxes represent PIV measurement regions and green points represent PDI and hot-wire measurement locations; D is the diameter of the disk. (Smith et al. 2021). 

was always 3 m downstream of the grid, thus assuring the same amount of background turbulence intensity. The disk was positioned in the centre of the tunnel cross-section for the HWA and PDI measurements, and 20 cm from the tunnel wall for the PIV measurements. The off-centre location of the disk for PIV measurements was due to hardware mounted to the top outer side of the tunnel preventing a central placement of the laser optic.
section_name: Experimental set-up

section_content: The inflow conditions without inertial particles were characterized with HWA measurements using a Dantec Streamline constant temperature anemometer with a tungsten wire probe. The Pt-W 55P01 type probe had a sensing length of 1. 25 mm and a wire diameter of 5 μm. Data acquisition time was 300 s at either 20 kHz or 35 kHz at constant temperature and pressure, and velocities were calculated using King's law and Taylor's hypothesis. Adequate resolution was achieved to resolve the Kolmogorov length scale η, as κ max η ≥ 1, where κ max is the highest resolvable wavenumber (rad m -1 ), estimated as 2πF s /U ∞ with F s the sampling frequency (Hz) and U ∞ the mean velocity (m s -1 ) (assuming Taylor's hypothesis). To quantify turbulence in the absence of particles, the hot-wire was aligned vertically with the centre of the disk at 36. 5 cm from the tunnel wall and floor. Measurements were taken in the tunnel at 3 m downstream of the grid with no wake (no disk model), and in the near and far wakes at 1D and 9. 6D, respectively, for tunnel speeds U ∞ = 2. 6, 4. 9, 8. 4, 10. 6, 12. 0 and 15. 8 m s -1. 

Table 1 presents the following turbulence statistics calculated for the background flow (grid-generated turbulence but with no wake):

(2. 1) Figure 4 shows the normalized power spectral density (PSD) normalized with the velocity variance (σ u ) 2 as a function of wavenumber κ for the measured instantaneous velocities for the background flow and disk wake locations at 1D and 9. 6D downstream. The presence of the wake increases the dissipation wavenumber cutoff, with the largest cutoff in the near wake at 1D. This indicates the wake is increasing energy dissipation and producing smaller eddies than the background flow. The inertial range is globally much clearer in the far wake at 9. 6D where turbulence is enhanced. The low frequency peak at 1D is of the order of 2 cm, which is the size of the hub in the centre of the disk. 

HWA calculations are plotted for tunnel speeds U ∞ = 4. 9, 8. 3, 10. 6, 12. 0 and 15. 8 m s -1. Figure 5 shows turbulence intensity vs the local mean streamwise velocity Ū for the background flow and wake locations at 1D and 9. 6D downstream of the disk. The local mean velocity Ū is reduced by approximately 54 %-66 % in the near wake at 1D and 7 %-19 % in the far wake at 9. 6D. Turbulence intensity remains fairly constant for the background flow and far wake, and varies with Ū in the near wake, indicating the near wake is a highly turbulent and non-homogeneous region. Figures 4(b ) and 5 confirm that that near-wake flow is within the recirculation bubble and therefore cannot be characterized via HWA in figures 6-11, as measurements are contaminated by velocity components other than the streamwise one. 

Figures 6-11 compare turbulence statistics (see table 1) for the far wake at 9. 6D with the background flow conditions. Figure 6 shows the Taylor scale Reynolds number Re λ as a function of local mean velocity Ū. For both cases Re λ increases with increasing Ū. 

Figure 7 shows ε, where the dissipation rate at 1D is an order of magnitude greater than the background flow and the far wake. Figures 7 and 8 show ε also increases with increasing Re λ. Figure 9 shows the transverse Taylor microscale λ as a function of Re λ. λ represents the intermediate length scale at which viscosity significantly affects turbulent eddies, and values of λ are smaller in the wake indicating the wake is generating more turbulent dissipation. The value of Re λ is greater in the wake due to greater σ u of the turbulent signal. In both cases λ decreases with increasing Re λ as higher velocities are also producing more turbulence. 

The Kolmogorov length scale η is shown in figure 10 as a function of Re λ. The value of η decreases with increasing Re λ for both cases, as η and λ are linked by the exact relation: λ/η = 15 1/4 Re 1/2 λ. For each Re λ, η is smaller in the wake as the disk injects energy to smaller scales where dissipation occurs. 

Figure 11 compares the dissipation equation constant C ε and L int /λ vs Re λ. The value of C ε is larger for the background flow where it ranges from 0. 83 to 1. 4, these values are in agreement with values presented in Sreenivasan (1998) and Vassilicos (2015) C ε is reduced by approximately half in the wake where it ranges from 0. 33-0. 41. The value of L int /λ steeply increases with increasing Re λ for both cases, due to the decrease in λ at higher Re λ. The linear slopes of L int /λ agree with Kolmogorov's 1941 theory that L int /λ ∼ C ε Re λ, and that the background flow and far wake have HIT and the Taylor hypothesis holds. HWA measurements presented in figures 6-11 show that in the absence of inertial particles, classical HIT scaling works reasonably well. The disk produces more turbulence, smaller eddies and greater dissipation of turbulent kinetic energy than the background flow.
section_name: Characterization of the turbulence in the wake

section_content: Table 2 lists case parameters for the PDI measurements. Water was delivered to the spray grid at 1. 7 and 2. 0 l min -1 and at three different free-stream tunnel velocities resulting in various global water volume fractions calculated as

, where Q L and Q G are the liquid and gas volumetric flow rates, respectively. The Stokes number based on the most probable particle diameter of 40. 5 μm (that remains approximately stable for all U ∞ and Φ v ) was calculated for the background flow (no wake) as

) and the minimum and maximum water volume fractions are listed for each tunnel velocity. This range of Stokes numbers (St ∼ 1) provides maximum coupling between turbulence and heavier inertial particles that can be centrifuged out of vortex cores and form clusters influenced by turbulent structures. 

The PDI measurement location is denoted with circles in figure 3 and was located 3 m downstream of the grid, 36. 5 cm from the floor and receiver side of tunnel and 38. 5 cm from the laser transmitter side. The PDI system remained stationary while the disk was positioned in the tunnel such that the PDI measurement location was centred behind the disk at 1D and 9. 6D downstream. PDI data were collected with an Artium Technologies PDI-200MD system, capable of detecting particle diameter ranges from 0. 5-2000 μm with a velocity range of -100 to +300 m s -1, with ±0. 2 m s -1 accuracy. Two diode pumped solid state lasers with wavelengths of 532 and 473 nm were split into two beams of equal intensity with the 532 nm beam oriented to measure vertical velocity and particle diameter, and the 473 nm beam oriented to measure horizontal velocity. The transmitter and receiver were positioned on opposite sides of the tunnel and had focal lengths of 1000 and 500 mm respectively. The receiver aperture was 200 μm. Signals were analysed with the Artium advanced signal analyser with a maximum sampling frequency of 320 MHz and resolution of 0. 01 % of sampling frequency; 500 × 10 3 signals were collected per case with the Automated Instrument Management System (AIMS) 5. 2 software to assure good statistical convergence. Particle diameters and velocities were measured for water flow-rates of 1. 7 and 2. 0 l min -1 at the free-stream tunnel speeds listed in table 2, leading to particle volume fraction varying in the range

2. 3. Two-phase flow experiments -PIV methodology PIV measurements were taken with water flow rates of 1. 2, 1. 7 and 2. 0 l min -1 and non-inertial tracer particles (Φ v = 0) at a constant tunnel fan speed as shown in table 3. Measurements were taken of the combined air-particle velocity by using the inertial particles as tracers in the two-phase flow cases. Tracer particles were used for the single-phase flow cases. Figures 14-20 compare tracer particle behaviour in the single-phase flow with the inertial particle behaviour in the two-phase flow. To prevent laser interference from water droplets adhering to the top inside tunnel window, a 12 mm angle was attached on the inside upper window upstream of the laser. The vertical laser sheet was aligned with the free stream and was centred with respect to the disk giving a velocity field in the x-y plane. Note that the wake is axisymmetric in the x-z plane, but not axisymmetric in the x-y plane due to the influence of the tube holding the disk. The model was painted black to prevent reflections, however, some reflections were still present resulting in a cropped region of interest for the near wake at 1D. Regions of interest for the near and far wakes are shown as rectangles in figure 3. PIV data were collected with the use of a Litron LD30-527, double pulse, Nd:YLF laser with a frequency of 3. 0 kHz, a wavelength of 527 nm and a t of 295-310 μs. For non-inertial tracer particles, the flow was seeded with two Antari alpha F-80Z smoke machines with Antari high density smoke Z-Fluid. For the inertial particle cases, the water droplets were used as seeding. A Phantom V2511 high speed camera with a 50 mm Nikon lens and resolution of 1280 × 800 pixels was set perpendicular to the laser sheet, and the camera captured images through an opening in the tunnel wall. The PIV grid size was x/D = y/D = 21 × 10 -4 and 8000 pairs of images were analysed to generate converged turbulence statistics. It was found from approximately 10 × 10 3 snapshots that the number of images needed for convergence of mean velocity was below 2000. Non-inertial particle measurements were analysed
section_name: Two-phase flow experiments -PDI methodology

section_content: Table 3. PIV Case Parameters for figures 14-20. 

-4 -3 -2 -1 0 1

Streamwise velocity u -(m s -1 )

Φ 1. 7 9. 6D wake Φ 1. 7 9. 6D background Φ 2 9. 6D wake
section_name: Porous disk turbulent wakes and inertial coupling

section_name: Results

section_content: The following results were obtained using PDI within the wake of the disk and in the background flow. Figure 12 shows the probability density function (p. d. f. ) of particle counts as a function of particle velocity for near and far wakes at a constant wind tunnel fan speed resulting in an average free-stream speed of 4. 9 m s -1. Symbols Φ 1. 7 and Φ 2 represent water volume fractions of 10. 1 × 10 -6 and 12. 3 × 10 -6 respectively, where the subscripts refer to 1. 7 and 2. 0 l min -1 water volumetric flow rates. 

Particle velocities in figure 12 are reduced in the near wake by approximately 144 %-141 %, and are negative, indicating recirculation. This flow reversal confirms why the hot-wire measurements cannot accurately characterize the near wake at 1D. Particles in the far wake have a velocity reduction of 7. 4 %-12. 2 % confirming that the wake is still Diameter (µm) present at 9. 6D downstream. Increasing the water volume fraction for each tunnel speed reduces the particle streamwise velocities as the bulk flow gains mass. 

Figure 13 shows p. d. f. s of particle counts as a function of particle diameter for two different water flow rates (columns) and 1D and 9. 6D wake locations (rows) for three different wind tunnel speeds. Most probable particle diameters range from 13 μm in the near wake to 41 μm in the background flow. For all cases, diameters are 2-10 μm larger for the higher flow rate, with smaller diameters at 1D (3-21 μm difference). At 9. 6D particle diameter distributions are similar between the wake and the background flow, with the most probable particle diameter between 30 and 41 μm. In contrast, the wake at 1D produces a much higher probability of particles around 17 μm, implying preferential trapping of smaller particles within the near-wake region. Local Stokes numbers in the near wake calculated with a particle diameter of 17 μm range from 0. 049 for U ∞ ≈ 3 m s -1 to 0. 17 for U ∞ ≈ 7. 8 m s -1, which are smaller than the Stokes numbers of the background flow. Particle trapping may be occurring in the wake at 1D as this was also seen at 1D in the wake of a model turbine by Smith et al. (2021). Particle trapping could be confirmed in future experiments with particle tracking velocimetry. Figures 12 and 13 illustrate recirculation of small particles at 1D, and that the wake still influences the particles at 9. 6D downstream.
section_name: Particle size and velocity distribution along the centreline of the wake

section_content: Turbulence statistics calculated from PIV measurements are compared for inertial particles (water droplets) and sub-inertial particles (smoke) in near and far regions of the disk wake. (p-t ) Φ 1. 2, Φ 1. 7 and Φ 2 represent two-phase particle velocity fields with increasing volume fractions with water flow rates of 1. 2, 1. 7 and 2 l min -1 and volume fractions of 1. 77 × 10 -6, 10. 3 × 10 -6 and 12. 3 × 10 -6 respectively. Note that the colour scales vary so that spatial features can be identified and the flow is from the left. 

Inertial particles were used as the seeding for PIV measurements to observe particle behaviour. Figures 14 and 15 show contour plots of various time-averaged normalized statistics at the same wind tunnel speed, which corresponds to Re λ = 88. 7 in the single-phase flow. Areas of interest span from x/D = 0. 75D-1. 75D for the near wake, and 9. 5D-11D for the far wake. Row Φ 0 shows contour plots of the single-phase flow wake; and rows Φ 1. 2, Φ 1. 7, and Φ 2 show contours of inertial particle behaviour in the wake with water flow rates of 1. 2, 1. 7, and 2 l min -1 and corresponding particle volume fractions of 7. 11 × 10 -6, 10. 1 × 10 -6 and 12. 3 × 10 -6, respectively. The centre of the disk is at y/D = 0 and the flow is from the left. Each column in figure 14 represents the following normalized statistics: mean streamwise velocity ū/U ∞, mean vertical velocity v/U ∞, and Reynolds stresses

Figure 15 shows the same quantities, but in the far wake. 

In figure 14, the near wake Φ 0 velocity remains positive at the core of the wake, with a distinct steep velocity gradient around y/D ∼ 0. 5, above which is the free stream where ū/U ∞ is unity. The lower half of the wake does not have a distinct edge due to the influence of the tube that attaches the disk to the tunnel floor. In contrast to the single-phase flow, the inertial particles reverse direction where they are trapped and recirculated behind the disk. (p-t ) Φ 1. 2, Φ 1. 7 and Φ 2 represent two-phase particle velocity fields with increasing volume fractions with water flow rates of 1. 2, 1. 7 and 2 l min -1 and volume fractions of 4. 26 × 10 -6, 6. 03 × 10 -6 and 7. 09 × 10 -6 respectively. Note that the colour scales vary so that spatial features can be identified and the flow is from the left. 

Inertial particles shift the top of the wake to y/D > 0. 7. Particle ū/U ∞ PIV measurements confirm particle recirculation in the near wake with regimes of negative ū, and reveal these regions grow in size with increasing volume fraction. 

As shown in the second column of figure 14, the near-wake normalized mean vertical velocity v/U ∞ remains negative (downward) for the Φ 0 case, with the maximum downward flow at y/D ∼ -0. 5. Inertial particles reveal recirculation with a band of negative v/U ∞ at the top of the wake at y/D ∼ 0. 5 and a band of positive v/U ∞ (upward) at the bottom of the wake at y/D ∼ -0. 5. 

The third column of figure 14 displays the normal Reynolds stress u u /U ∞, which represents the contributions of streamwise turbulence fluctuations to momentum within the single-phase flow wake. In the near wake, the single-phase flow normal stress is near zero at y/D ∼ 0. 3, and increases downward and upstream (left). The particle velocity fields show bands of normal stress at the top and bottom edges of the wake, at y/D ∼ ±0. 5. Interestingly, the middle volume fraction Φ 1. 7 differs from the other particle cases, with greater stress in the lower band of the wake, and less stress in the upper band. The fourth column of figure 14 shows shearing Reynolds stress u v /U ∞, which represents vertical fluctuation influence on streamwise turbulent momentum for single-phase flows. In the near-wake single-phase flow, the shearing stress is negative at y/D ∼ 0. 1 and three bands of positive shear stress appear at y/D ∼ -0. 2, -0. 5 and -0. 9. The inertial particle velocity fields show pronounced bands of low and high u v /U ∞ with a negative band at y/D ∼ 0. 5 and a wider positive band at y/D ∼ -0. 5. 

The column on the right in figure 14 shows v v /U ∞, which is the vertical Reynolds normal stress and represents energy contributions within the single-phase flow wake from the vertical turbulent fluctuations. The single-phase flow v v /U ∞ in the near wake, is approximately zero above y/D ∼ 0. 25, and increases as you move downward and to the left. The near-wake particle velocity field v v /U ∞ forms bands at the top and bottom of the wake, but is less distinct than u v /U ∞, and increases with increasing volume fraction. In the near wake Φ 1. 2 case, two distinct bands of higher v v /U ∞ are centred behind the disc at y/D ∼ 0. 3 and -0. 3. 

In figure 15, the far-wake single-phase flow (Φ 0 ), the core tapers downward so that the centre of the wake core is around y/D ∼ -0. 3. The top edge of the wake is clearly defined and the bottom of the wake has moved lower than the region of interest, indicating the wake grows wider as it moves downstream. In the inertial particle cases, the far wake ū/U ∞ is lower in magnitude and positive as particles move downstream. All three volume fractions have similar contours that show a reduction in ū/U ∞ toward the bottom of the region of interest, suggesting the particle far wakes have shifted downward out of the region of interest. Figure 12 obtained with the PDI confirms that the wake is still present at 9. 6D, and the differences in particle velocity are small as one varies the volume fraction, which is in agreement with figure 14. 

The far wake v/U ∞ in figure 15, decreases with height for the Φ 0 case and a local minimum occurs below y/D ∼ -0. 5 and below x/D ∼ 10. For the inertial particle cases, v/U ∞ is more negative at the top edge of the wake, and approaches but does not reach zero below y/D ∼ -0. 5. 

In figure 15, the behaviour of u u /U ∞ in the far wake single-phase flow is different than the near wake, where a band of higher u u /U ∞ is clearly visible at y/D ∼ 0. 25, and it increases as it moves downstream. Note that the vertical line of slightly higher values at x/D ∼ 10. 2 is due to a laser reflection. Generally, the far-wake inertial particle cases have lower magnitudes for the middle volume fraction Φ 1. 7, and show larger bands of greater normal stress in the lower part of the region of interest below y/D ∼ -0. 5. All inertial particle cases have a thinner but distinct band of higher u u /U ∞ at y/D ∼ 0. 5. Inertial particle velocity fields have different structures than the single-phase flow, where particles tend to have greater u u /U ∞ at the top edge of the wake. 

In the far wake in figure 15, the magnitude of the single-phase v v /U ∞ is an order of magnitude lower than the near wake. In the single-phase flow, bands of positive and negative shear stress are wider with less defined edges as the wake expands downstream. The Φ 0 case has an upper band of negative shear stress centred at y/D ∼ 0. 2 and a lower band of positive shear stress centred at y/D ∼ -0. 7, with bands of zero shear stress at y/D > 0. 5 and ∼ -0. 25. The far-wake particle velocity fields also have bands of negative shear stress at the top of the wake, and positive shear stress below, but these bands have moved away from the centre of the wake, with the upper band centred around y/D ∼ 0. 53 and less defined lower bands of positive shear stress starting below y/D ∼ -0. 45. Both the Φ 1. 2 and Φ 1. 7 cases have a band of zero shear stress at y/D ∼ 0, but this band moves downward to y/D ∼ -0. 3 for the Φ 2 case. Travis and others -0. 5 0 0. 5 1. 0 1. 5

Figure 16. Normalized mean turbulence statistics as a function of height at 1D (solid lines) and 9. 6D (dashed lines). Here, Φ 0 represents single-phase flow, while Φ 1. 2, Φ 1. 7 and Φ 2 represent two-phase flow with water flow rates of 1. 2, 1. 7 and 2 l min -1 and volume fractions of 4. 26 × 10 -6, 6. 03 × 10 -6 and 7. 09 × 10 -6, respectively. Flow direction is from the left, and y/D = 0 is the centre of the disk. 

In figure 15 far wake Φ 0 case, v v /U ∞ occurs below y/D ∼ 0. 5. In the two-phase flow cases, the greatest values of v v /U ∞ occur in the bottom half of the wake. The largest volume fraction has a different structure than the others with the greatest intensity of v v /U ∞ in a thin band just above y/D ∼ 0. 

Vertical profiles of quantities displayed in figures 14 and 15 are represented in figures 16 and 17, and are averaged spatially over 2 mm. In figures 16 and 17, single-phase flows (grey and black lines) are compared with the three different volume fractions of inertial particles (coloured lines) at 1D in the near wake (solid lines) and 9. 6D in the far wake (dashed lines). Figure 16 compares the near and far wakes and figure 17 compares different volume fractions in the far wake. 

In the second column of figure 16, the inertial particle 1D vertical profiles have the opposite sign of the Φ 0 single-phase flow grey line. The absence of recirculation in the far wake allows for the measurement of the particle settling velocity. The second column in figure 17 reveals the settling velocity is highest for the smallest volume fraction. All far-wake cases below y/D ∼ -0. 6 have similar v/U ∞ magnitude. The Φ 1. 2 and Φ 1. 7 cases have nearly identical profiles and magnitudes, while the greatest volume fraction (Φ 2 ) profile diverges from this pattern and has the lowest negative velocity. This may imply a dependence of settling velocity on volume fraction. We nevertheless refrain from doing a systematic study on this parameter, as it has a non-trivial dependency with the finite size effect generated by the mass of water moving through the wind tunnel (Mora et al. 2021). Results may present different biases with U ∞ and φ v and we therefore leave such study for future works. 

Figure 17. Normalized mean turbulence statistics as a function of height at 9. 6D. Here, Φ 0 represents single-phase flow, while Φ 1. 2, Φ 1. 7 and Φ 2 represent two-phase flow with water flow rates of 1. 2, 1. 7 and 2 l min -1 and volume fractions of 4. 26 × 10 -6, 6. 03 × 10 -6 and 7. 09 × 10 -6, respectively. Flow direction is from the left, and y/D = 0 is the centre of the disk. 

The Reynolds shear stresses in figure 16 for the near wake reveal particle cases have distinct peaks near y/D ∼ ±0. 5. While figure 17 shows the far wake particle Reynolds shear stress is greater in magnitude and has a different profile shape for the Φ 2 case. The particle velocity fields show different patterns of u v /U ∞ than the single-phase flow, with bands of negative and positive particle field shear stress concentrating at the top and bottom of the near wake, respectively. 

Figure 17 shows that, although the changes in far-wake particle magnitudes are small, the Φ 2 case has greater v v /U ∞ in the top half of the wake, and lower v v /U ∞ below. Figure 17 v/U ∞ and v v /U ∞ profiles show that the far-wake inertial particle cases are similar in the lower part of the wake, but above y/D = 0 the Φ 2 case diverges from the other two with a lower settling velocity and higher v v /U ∞.
section_name: Velocity statistics of the global wake modulation through particle injection

section_content: Quadrant analysis is presented here to characterize the relationship of u and v in the wake. Quadrant analysis was first presented by Wallace, Eckelmann & Brodkey (1972) and applied to a turbulent boundary layer near the wall. The method decomposes the fluctuating velocity signal into four events with conditional averaging and is described in works by Yue et al. (2007)

Figure 18. Vertical profiles of quadrant analysis at wake location 1D. Dashed lines represent each quadrant for the single-phase flow case (Φ 0 ), while solid lines represent each quadrant for Φ 1. 2 (a), Φ 1. 7 (b) and Φ 2 (c) which represent two-phase flow with water flow rates of 1. 2, 1. 7 and 2 l min -1 and volume fractions of 7. 11 × 10 -6, 10. 3 × 10 -6 and 12. 3 × 10 -6, respectively. Flow direction is from the left, and y/D = 0 is the centre of the disk. separated by the sign of the fluctuating velocities into positive and negative event quadrants shown in table 4. 

These events are interpreted as directionality of the fluctuating flow compared with the mean flow and are presented in figures 18 and 19 for wake locations at 1D and 9. 6D, respectively. Dashed lines represent each quadrant for the Φ 0 case, and solid lines represent quadrants of the two-phase flow cases for each particle volume fraction. 

Figure 18 shows that for the near-wake Φ 0 case, Q2 and Q4 (ejections and sweeps) have a local maximum at just above y/D = 0 at the centre of the wake. In the lower half of the wake, Q1 and Q3 (outward and inward interactions) are greatest and increase in

Figure 19. Vertical profiles of quadrant analysis at wake location 9. 6D. Dashed lines represent each quadrant for the single-phase flow case (Φ 0 ), while solid lines represent each quadrant for Φ 1. 2 (a), Φ 1. 7 (b) and Φ 2 (c) which represent two-phase flow with water flow rates of 1. 2, 1. 7 and 2 l min -1 and volume fractions of 7. 11 × 10 -6, 10. 3 × 10 -6 and 12. 3 × 10 -6, respectively. Flow direction is from the left, and y/D = 0 is the centre of the disk. magnitude with decreasing height. For the two-phase flow cases, Q2 (ejections) have the greatest magnitude except at y/D = -0. 5 for the Φ 2 case. Generally, two-phase flow Q2 and Q4 have a local maximum at y/D = 0. 5, with Q2 and Q4 increasing with increasing volume fraction. The lower half of the two-phase flow wakes are dominated by Q2 for the two lower volume fractions. The Φ 2 case has a local maximum just below y/D = -0. 5 for Q1 and Q3 which are greater than Q2 and Q4. 

Quadrant analysis for the far wake is shown in figure 19. Here, the magnitudes of the two-phase flow cases are larger than the single-phase flow case compared with the near wake. The single-phase flow case has a peak of Q2 and Q4 in the top half of the wake, and Q1 and Q3 in the bottom half. Otherwise, patterns are similar to the near wake with Q2 dominating the Φ 1. 2 and Φ 1. 7 cases and the greatest magnitude of Q2 occurs at y/D = 0. 5 for the largest volume fraction. The magnitudes of the two-phase flow cases increase with increasing volume fraction. Overall, ejections (Q2) are the dominant event for the two-phase flow cases, and represent fluctuating events that push upward and backward compared with the mean flow. Q2 magnitudes are greatest at the top edge of the wake at y/D = 0. 5. 

Further comparison of ejections (Q2) and sweeps (Q4) is obtained with S 0 = S 4 -S 2 shown in figure 20. Introduced by Raupach (1981) The dashed lines represent the single-phase flow case (Φ 0 ), while solid lines represent two-phase flow cases Φ 1. 2, Φ 1. 7 and Φ 2, with water flow rates of 1. 2, 1. 7 and 2 l min -1 and volume fractions of 7. 11 × 10 -6, 10. 3 × 10 -6 and 12. 3 × 10 -6, respectively. Flow direction is from the left, and y/D = 0 is the centre of the disk. 

ejection events and values less than zero represent dominant sweeps. For both the near and far wakes, the top edge of the wake at y/D = 0. 5 is dominated by ejection events which can also be observed in figures 18 and 19. In the near wake (figure 20a ), the two-phase flow cases have greater magnitudes and different behaviour than the single-phase flow case. The centre of the wake is dominated by sweeps, and then transitions back to ejections at the bottom edge of the wake at y/D = -0. 5. In the far wake (b), the single-phase flow case has ejections at the wake edges and sweeps in the centre, while the two-phase flow cases are all ejection events. The largest volume fraction ejection events are three times greater in magnitude than the lower volume fraction cases at the top edge of the wake, and smaller than the lower volume fraction cases at the bottom edge of the wake.
section_name: Quadrant analysis of global wake modulation through particle injection

section_content: Hot-wire results (figures 4-11) confirm HIT behaviour in the far wake and background flow (no wake), with anisotropic behaviour in the near wake. As expected, the single-phase wake reduces Ū, increases turbulence intensity, and generates more turbulent dissipation than the background flow. 

Poly-disperse inertial particles were introduced to the flow in particle volume fractions from Φ v = 6. 31 × 10 -6 -1. 86 × 10 -5. Particles have negative horizontal velocity at 1D, indicating recirculation in the near wake, and reduced particle velocities compared with the background flow show the wake is still present at 9. 6D downstream. The PDI measurements also show that increasing Φ v increases the probability of larger particle diameters. Background flow particle diameters were measured on average to be 9 μm bigger in the larger Φ v case, and may be due to droplet coalescence. Particles are smaller in the wake than the background flow, and smallest in the near wake. This indicates the near-wake recirculation may be trapping smaller particles. Interestingly, the near-wake particle diameters are smallest at U ∞ = 4. 8 m s -1, and larger at U ∞ = 3. 2 and 5. 6 m s -1, suggesting a Reynolds number influence on particle diameter entrainment. 

Both PIV measurements and PDI measurements (figures 12 and 14) show particle ū/U ∞ is reversed in the near wake. Interestingly, negative particle ū/U ∞ does not occur in the single-phase flow PIV measurements. This is probably due to the fact that the water droplets adhered to the disk and blocked the smallest holes near the centre of the disk, Porous disk turbulent wakes and inertial coupling creating a larger recirculation zone. PIV measurements reveal that the region of negative particle ū/U ∞ increases in size in the y/D direction with increasing Φ v, and the blockage may have increased as more particles impacted the disk at higher volume fractions. The recirculation zone behind the disk in the single-phase flow was most likely smaller and to the left of the measurement window. 

Particle recirculation in the near wake is also apparent in the figure 14 v/U ∞ contour plots, where particles in the top half of the wake have negative vertical velocity (downward), and particles in the lower half of the wake have positive vertical velocity (upward). These regions of particle v/U ∞ widen and intensify with increasing volume fraction. 

In the far wake, where turbulence may be starting to approach HIT conditions, particle vertical velocity (figures 15 and 17) is smallest for the largest volume fraction, which may be in agreement with findings in Sumbekova et al. (2016) and Good et al. (2014). As stated above, a more detailed study is needed to disentangle the dependency of the settling velocity with U ∞ and Φ v with the finite-size effects from the wind tunnel. This lower vertical velocity result is interesting, and occurs in the top half of the wake where u v /U 2 ∞ is negative, and u u /U 2 ∞ and v v /U 2 ∞ are greatest for the greatest volume fraction. Possible explanations for this behaviour are nonlinear drag on particles, the loitering effect where falling particles spend more time in regions with upward flow, or regions with low preferential sweeping (Stout, Arya & Genikhovich 1995; Rosa et al. 2016). Preferential sweeping is enhancement of settling velocity where inertial particles accumulate on the downward side of eddies. Rosa et al. (2016) found in DNS simulations that when preferential sweeping is disabled, the settling velocity of droplets is reduced and droplets smaller than 55 μm are sensitive to ε. 

In the single-phase flow wake in figures 14 and 15, u u /U ∞, u v /U ∞ and v v /U ∞ are greater in magnitude in the near wake and smaller in the far wake, indicating streamwise fluctuations and Reynolds shear stress decrease as the wake moves downstream. The particle u u /U ∞, u v /U ∞ and v v /U ∞ are also greater in the near wake and concentrated where ū/U ∞ has a steep gradient near y/D ∼ ±0. 5. Regions of stress at the top and bottom edges of the single-phase flow wake are where the greatest particle field stresses are located. Figure 17 shows particle u v /U 2 ∞ and v v /U 2 ∞ are similar in magnitude for all volume fractions in the lower part of the wake, and deviate for the greatest volume fraction in the top half of the wake, suggesting a Reynolds stress influence on the settling velocity. 

Quadrant analysis also shows a difference in behaviour when particles are added to the wake. In the near wake, while inward and outward interactions dominate the lower half of the single-phase flow, ejections are most prominent in the two-phase flow cases, especially at the top edge of the wake. The parameter S 0 shows that ejection events are most prevalent at the top and bottom edges of the near wake and ejection events are most prevalent everywhere for the two-phase flow cases in the far wake. Ejection events (Q2) are turbulent bursts upward at velocities slower than the mean, and may be contributing to the slower settling velocity of the greatest volume fraction, as this case has the greatest magnitude of Q2. 

Comparison of turbulence statistics between single-phase and two-phase flow gives insight into how the turbulent wake and particles interact. There are distinct differences between the single-phase flow and the inertial particle velocity fields, as well as differences in the particle fields with changing volume fraction. The lack of symmetry between values at y/D ∼ ±0. 5 may be due to the presence of the tube that attaches the disk to the tunnel floor and the fact that particles are falling as they move downstream. Although the particle diameters are smaller than η for all PIV results presented, the Stokes numbers increase with increasing ε, and ε increases in the presence of the wake. The calculated respective near-and far-wake Stokes numbers are 0. 77 and 1. 7, which are larger compared with the background flow (no wake) of 0. 58. This suggests that the Stokes number influences the particle/turbulence interaction, and particles could be affecting the turbulence in the wake. Preferential concentration is a phenomenon well known to happen in turbulent flows seeded with inertial particles (Squires & Eaton 1991; Monchaux, Bourgoin & Cartellier 2012), and could lead to concentrated regions of higher Φ v. This effect may onset two-way coupling effects and explain why the particle fields have distinctly different behaviour than the background flow. However, our experimental set-up does not allow us to quantify two-way coupling with the measurements presented here.
section_name: Discussion

section_content: This study presents the findings of a wind tunnel experiment investigating the behaviour of inertial particles in the turbulent wake of a porous disk. The incoming flow was characterized as homogeneous, isotropic turbulence, and the wake of the disk was measured in regions around 1D and 9. 6D downstream. Poly-dispersed particles (water droplets) were injected uniformly upstream of the disk at varying volume fractions. Both single-phase and two-phase flow measurements were taken with HWA, PDI and PIV to investigate the wake particle interaction. 

Hot-wire measurements showed the local mean streamwise velocity is reduced by 54 %-66 % at 1D and 7 %-19 % at 9. 6D. In single-phase flow, the wake produces smaller turbulent eddies and increases energy dissipation. Essentially, the disk acts as a turbulence generator, with anisotropic flow in the near wake. 

Analysis of PDI measurements found the most probable particle diameter at 1D was 17 μm, while the most probable particle diameters in the far wake ranged from 30-41 μm, suggesting smaller particles are trapped in the near wake. Particle velocities were reversed at 1D and reduced by 41 %-44 %, revealing that particles are recirculating in the near wake. At 9. 6D particle velocities were reduced by 7. 4 %-12. 2 % compared with particles in the background two-phase flow (no wake). These measurements confirm there is still a wake/particle interaction at 9. 6D downstream. 

PIV contours also confirm a reversal of particle ū/U ∞ at 1D, and a reduced streamwise velocity at 9. 6D. Particle field v/U ∞ also shows particle recirculation at 1D with horizontal bands of downward particle velocity in the top half of the wake and upward particle velocity in the lower half. Particle Reynolds stress magnitudes tend to be greatest at y/D ∼ ±0. 5 at the edges of the wake where the gradient of ū/U ∞ is steep and shearing occurs. This supports the idea that particles cluster in regions of high strain and low vorticity. 

It is also shown that at 9. 6D, particle downward velocity was smallest for the largest volume fraction in the upper portion of the wake. This occurs where shear stress is negative and vertical stress increases with increasing volume fraction. Quadrant analysis suggests that ejection events may be contributing to this phenomenon, as these are upward fluctuations at the top edge of the wake. Whether particle settling velocity is enhanced or reduced is greatly debated in the literature, and depends on many parameters such as drag nonlinearity, vortex trapping/centrifuging, loitering, preferential sweeping, particle size, Stokes number, Re λ, dissipation rate ε and volume fraction. These parameters vary throughout the wake and are a starting point for further exploration. fractions capable of modifying turbulence (two-way coupling) and give more insight into the differences observed between the single-phase flow and inertial particle velocity fields within the wake. 

These results show that at Re D = 6. 50 × 10 4, there are distinct differences between turbulence statistics of a single-phase wake behind a porous disk, and the behaviour of the combined bulk flow containing inertial particles at volume fractions of Φ v ∈ ([6 -19] × 10 -6 ). We see that particles in near and far wake have a very different dynamics. These findings imply that there will be many challenges to using porous disks as analogues to rotating turbines if inertial particles are present in the turbulent flow. Furthermore, both the wake and particles seem to be affected by volume fraction. We therefore open the door for the study of inertial particles in axisymmetric wakes and a fundamental interest in terms of self-similar flows.
section_name: Conclusions

section_content: Acknowledgements. We would like to acknowledge V. Govart for machining experimental fixtures. Declaration of interests. The authors report no conflict of interest. Author ORCIDs. Kristin N. Travis https://orcid. org/ 0000-0002-7776-2586 ; Mickaël Bourgoin https://orcid. org/ 0000-0001-9442-7694 ; Martín Obligado https://orcid. org/ 0000-0003-3834-3941.