Hostname: page-component-cd9895bd7-fscjk Total loading time: 0 Render date: 2024-12-26T16:24:29.539Z Has data issue: false hasContentIssue false

Three-dimensional flow around and through a porous screen

Published online by Cambridge University Press:  16 May 2024

Olivier C. Marchand*
Affiliation:
Laboratoire d'Hydrodynamique (LadHyX), UMR7646, CNRS, École polytechnique, Institut Polytechnique de Paris, 91120 Palaiseau, France
Sophie Ramananarivo
Affiliation:
Laboratoire d'Hydrodynamique (LadHyX), UMR7646, CNRS, École polytechnique, Institut Polytechnique de Paris, 91120 Palaiseau, France
Camille Duprat
Affiliation:
Laboratoire d'Hydrodynamique (LadHyX), UMR7646, CNRS, École polytechnique, Institut Polytechnique de Paris, 91120 Palaiseau, France
Christophe Josserand
Affiliation:
Laboratoire d'Hydrodynamique (LadHyX), UMR7646, CNRS, École polytechnique, Institut Polytechnique de Paris, 91120 Palaiseau, France
*
Email address for correspondence: [email protected]

Abstract

We investigate the three-dimensional (3-D) flow around and through a porous screen for various porosities at high Reynolds number $Re = {O}(10^4)$. Historically, the study of this problem has been focused on two-dimensional cases and for screens spanning completely or partially a channel. Since many recent problems have involved a porous object in a 3-D free flow, we present a 3-D model initially based on Koo & James (J. Fluid Mech., vol. 60, 1973, pp. 513–538) and Steiros & Hultmark (J. Fluid Mech., vol. 853, 2018 pp. 1–11) for screens of arbitrary shapes. In addition, we include an empirical viscous correction factor accounting for viscous effects in the vicinity of the screen. We characterize experimentally the aerodynamic drag coefficient for a porous square screen composed of fibres, immersed in a laminar air flow with various solidities and different angles of attack. We test various fibre diameters to explore the effect of the space between the pores on the drag force. Using PIV and hot wire probe measurements, we visualize the flow around and through the screen, and in particular measure the proportion of fluid that is deviated around the screen. The predictions from the model for drag coefficient, flow velocities and streamlines are in good agreement with our experimental results. In particular, we show that local viscous effects are important: at the same solidity and with the same air flow, the drag coefficient and the flow deviations strongly depend on the Reynolds number based on the fibre diameter. The model, taking into account 3-D effects and the shape of the porous screen, and including an empirical viscous correction factor that is valid for fibrous screens may have many applications including the prediction of water collection efficiency for fog harvesters.

Type
JFM Papers
Copyright
© The Author(s), 2024. Published by Cambridge University Press

1. Introduction

The flow around porous structures has been largely investigated throughout the recent decades and has many engineering applications. It can be applied to parachute problems for the determination of drag and stability (Sarpkaya & Lindsey Reference Sarpkaya and Lindsey1990; Johari & Desabrais Reference Johari and Desabrais2005), to vertical axis wind turbines (Ayati et al. Reference Ayati, Steiros, Miller, Duvvuri and Hultmark2019) as well as to blockage correction in wind tunnels (Steiros, Bempedelis & Cicolin Reference Steiros, Bempedelis and Cicolin2022). Laws & Livesey (Reference Laws and Livesey1978) highlighted the possibility of using screens in flow to control the velocity distribution and change the flow direction. The understanding and improvement of the water collection of fog harvesters in arid regions require a quantitative description of the flow in the vicinity of the net (Regalado & Ritter Reference Regalado and Ritter2016; Moncuquet et al. Reference Moncuquet, Mitranescu, Marchand, Ramananarivo and Duprat2022). Furthermore, such a quantitative description may provide a first step in the physical understanding of respiratory flows in the presence of a face mask, as used to reduce the propagation of an airborne virus such as SARS-CoV-2 (Mittal, Ni & Seo Reference Mittal, Ni and Seo2020; Bourrianne et al. Reference Bourrianne, Xue, Nunes, Abkarian and Stone2021). In these cases, the flow can either pass through the porous net or mask, or is deviated around or through the leaks. More generally, there may be an interest in reducing the constraints exerted on high panels or masts exposed to a flow for safety reasons, leading to an increased interest in large porous structures as mentioned by Giannoulis et al. (Reference Giannoulis, Stathopoulos, Briassoulis and Mistriotis2012). In some buildings, a permeable layer is added at a certain distance from the façade for energy efficiency reasons or to block a part of the sun rays. Also, windbreak panels are usually used in industry and power plants to control wind and dust pollution. In these cases an estimation of the cladding wind load is useful information for architects and engineers (Pomaranzi et al. Reference Pomaranzi, Daniotti, Schito, Rosa and Zasso2020). Furthermore, the modelling of the interaction of flow with arrays of fibres and the prediction of the corresponding drag can be helpful to understand the mechanism of filter feeding for numerous marine organisms for which arrays of bristles move in water to capture food particles (see Cheer & Koehl Reference Cheer and Koehl1987; Hood, Jammalamadaka & Hosoi Reference Hood, Jammalamadaka and Hosoi2019).

The main physical characteristics of the screen involved in flow resistance is the porosity and the permeability. For very thin porous screens, the porosity can be represented by the solidity, which is the ratio between the solid surface area and the total surface area of the screen. The permeability is defined as the parameter relating the pressure gradient within a porous media to the local velocity of the flow, depending on the geometry of the pores. The porosity effects on pressure drop across porous screens and on drag force has largely been investigated, both theoretically and experimentally, whereas the permeability effect is much less understood, especially for very thin porous structures. However, some recent work at low Reynolds numbers has been conducted highlighting its influence: in particular, Ledda et al. (Reference Ledda, Siconolfi, Viola, Gallaire and Camarri2018) has shown that the permeability has a strong effect on the wake characteristics, and Pezzulla et al. (Reference Pezzulla, Strong, Gallaire and Reis2020) has demonstrated that the drag coefficient of the screen depends on the permeability. On the other hand, Steiros & Hultmark (Reference Steiros and Hultmark2018) have developed a model to predict the drag coefficient for two-dimensional (2-D) perforated plates as a function of the solidity only. Although they obtained a good agreement with experimental data, this approach does not take explicitly into account the Reynolds number based on the scale of the pores, while it is known that the behaviour of the flow in the vicinity of the screen depends on the specific geometry of the pores as well as the material and thickness of the screen (see § 4.1).

Several approaches have been adopted to model the flow. For instance, Carvajal et al. (Reference Carvajal, Silva-Llanca, Larraguibel and González2020) used three-dimensional (3-D) numerical simulations to access the aerodynamic characteristics during fog collection and model the net as a porous medium using Darcy's law, whereas de Dios Rivera (Reference de Dios Rivera2011) used the superposition principle applied to a flow passing around a solid plate and a flow forced to pass through the net to find an approximation of the velocity at the screen. The difficulty of the problem lies in the multiscale physical phenomena, from the characteristic scale of the flow around the screen that is of order $0.1$ to 1 m to that of the flow through the screen (i.e. the pore size) that is of order ${O}(10^{-3})$ to ${O}(10^{-6})$ m; the scale of the Reynolds number thus varies from low Reynolds numbers at the pore scale (of order ${O}(10^{-1})$) to large ones at the net scale (of order ${O}(10^{6})$). Therefore, the physical mechanisms of the whole system cannot be easily captured by a numerical simulation resolving all scales as noticed by Shklyar & Arbel (Reference Shklyar and Arbel2008), and the method generally used consists in modelling the porous surface as an imaginary interface where transfers of mass and momentum occur. The macroscopic jump laws for the velocities and pressure at the interface are deduced from a microscopic model at the vicinity of the screen where the fluid is generally governed by the steady Stokes equations. This can be obtained by periodic homogenization theory as done by Zampogna & Gallaire (Reference Zampogna and Gallaire2020) and has been recently used for permeable shells by Ledda et al. (Reference Ledda, Boujo, Camarri, Gallaire and Zampogna2021). However, these methods introduce some parameters like the permeability of the porous surface that are difficult to measure experimentally, although they could be obtained using pore-scale simulation.

We aim here at predicting the drag coefficient of elevated porous panels of arbitrary shape placed in a laminar flow, and the corresponding flow through the structure, in conditions relevant for many applications (e.g. fog harvesting). We focus on porous screens composed of fibres, but the model can be applied without major changes to other kinds of porous screens. Experimentally, we consider rectangular meshes of woven fibres, with fibre diameters between $6.0\,\mathrm {\mu }$m and 1.9 mm and with typical pore sizes of the order of $10\,\mathrm {\mu }$m to 1.0 cm, placed in a uniform laminar flow of velocity varying from $0.5$ to $13\,{\rm m}\,{\rm s}^{-1}$ with three orientation angles. For a given fibre diameter, we vary the porosity of the screen by changing the fibres arrangements and pore sizes; conversely, we study screens of the same porosity but different fibre diameters, in order to study the effect of the screen porosity and the flow around or through the screen, as well as the effect of the microstructure, i.e. of the local pore Reynolds number associated to the flow through the individual pores. We measure the drag coefficient and perform velocimetry measurements of the flow around the screen with particle image velocimetry (PIV) and a hot wire probe.

Our goal is to obtain a model of the flow for porous surfaces based only on the porosity, the large-scale geometry of the screen and the Reynolds number at the scale of the holes (instead of the permeability), which are easy to access. From this model, we wish to obtain values of the drag coefficient as well as the proportion of the flow that is deviated around the screen. We focus on the model first proposed by Taylor (Reference Taylor1944) and used by numerous authors (O'Neill Reference O'Neill2006), which consists in considering the screen as a distribution of sources. This approach has been adopted by Koo & James (Reference Koo and James1973) who proposed a 2-D mathematical model for a screen confined in a channel with two parallel boundaries. Recently, Steiros & Hultmark (Reference Steiros and Hultmark2018) derived the drag coefficient of a porous plate based partially on Koo & James (Reference Koo and James1973) and 2-D potential flow. Since their prediction showed a good agreement with their experimental data, here we propose to extend the model of Koo & James (Reference Koo and James1973) to a 3-D free flow, keeping the same main hypothesis but taking into account the shape of the screen, 3-D effects and the base suction effect considered by Steiros & Hultmark (Reference Steiros and Hultmark2018). We then derive the equations to predict the drag coefficient following the method proposed by Steiros & Hultmark (Reference Steiros and Hultmark2018), and discuss the limits of this approach. While the pressure jump derived following the method of Taylor (Reference Taylor1944) and Steiros & Hultmark (Reference Steiros and Hultmark2018) is based on clear physical assumptions, it fails to reproduce the viscous effects observed experimentally, i.e. finite pore Reynolds number. We thus incorporate in the model the empirical law of Brundrett (Reference Brundrett1993) for the pressure jump. We compare the prediction of the drag coefficient and of the proportion of deviated fluid for various porosities and Reynolds number based on the fibres’ diameter to our experimental results (both with the theoretical and empirical law of the pressure jump).

Furthermore, very few studies report measurements on the drag coefficient for elevated panels with different angles of attack (to the best of our knowledge only Prandtl & Flachsbart Reference Prandtl and Flachsbart1932; Letchford Reference Letchford2001). We therefore also study the flow for arbitrary angle of attacks and compare our theoretical prediction with experimental data for two angles of attack ($65^\circ$ and $43^\circ$). Finally, we discuss the asymptotic behaviour of the model at solidities close to 1.

2. Experiments

We measure experimentally the drag coefficient for a series of porous planar structures consisting of regularly woven nylon yarns in a square mesh and other types of meshes like rod screens (parallel fibres), of size $L=10\,{\rm cm}\times 10\,{\rm cm}$. The characteristics of the porous screens as well of their solidity are detailed in Appendix C, and we also present a synthesis in table 1. The solidity $s$ is defined as the ratio between the solid surface area of the screen to its total area $S_p$. We also characterize the drag coefficient for a classical surgical face mask for which the physical characteristics such as the fibre diameter and the solidity are taken from Monjezi & Jamaati (Reference Monjezi and Jamaati2021) and Du et al. (Reference Du, Iacoviello, Fernandez, Loureiro, Brett and Shearing2021) (screen P9 in table 1). The porous structure is held in a planar configuration by a square frame with a width of 0.5 cm (that is 2.1 % of the surface area of the porous structure), representing a small portion (7.6 %) of the total cross-section of the laminar flow generated by the wind tunnel. Following Letchford (Reference Letchford2001), the square frame is fixed on a 21.5 cm high mast to avoid boundary layer effects. The set-up is shown in figure 1.

Table 1. Porous screen characteristics. The interval of the Reynolds number $Re_d$ is calculated with the fibre diameter $d$ as characteristic size, with a velocity $v_0 = 5.0\,{\rm m}\,{\rm s}^{-1}$ for the minimum value, a velocity $v_0 = 13.0\,{\rm m}\,{\rm s}^{-1}$ for the maximum value and a kinematic viscosity $\nu = 15.6 \times 10^{-6}\,{\rm m}^2\,{\rm s}^{-1}$. For the screens P2, P12, P26, P27, P28 P29 and P30, we put the detailed values of the drag coefficient at each velocity (and calculated the corresponding Reynolds number $Re_d$). The uncertainties as well as the fibre diameters, material type and geometry can be found in the Appendix C.

Figure 1. Diagram of the experimental set-up for the drag coefficient measurement, and examples of porous screens essentially in nylon woven mesh.

A force balance (SIXAXES, FX2.6, No 1026, $\pm$5 N, sensitivity of $\pm$0.001 N) is used to measure the force applied to the whole system. The laminar airflow is generated by an open jet wind tunnel with a square test section of width 40 cm. We also perform experiments on inclined porous screens; in that case eight different velocities have been used. For both cases, the velocities vary from $0.5$ to $13\,{\rm m}\,{\rm s}^{-1}$. We can thus define the (global) Reynolds number of the problem as

(2.1)\begin{equation} Re=\frac{L_0v_0}{\nu}, \end{equation}

where $v_0$ is the uniform velocity of the flow far upstream of the screen, $L_0$ its typical size and $\nu$ the kinematic viscosity of the fluid (the air for the configurations considered here, so $\nu = 15.6 \times 10^{-6}\,{\rm m}^2\,{\rm s}^{-1}$). For a screen of the size of a few tens of centimetres with a velocity of the order of $1\,{\rm m}\,{\rm s}^{-1}$, we obtain $Re\approx {O}(10^4) \gg 1$. We perform measurements with screens of different solidities, as well as made of different fibre radii and pore sizes while keeping the solidity almost constant, which allows us to probe the local effect of viscosity. The details of the wind tunnel as well as the characteristics of the flow can be found in the thesis of Du Pontavice (Reference Pontavice2016). The system is placed at a distance of approximately 50 cm from the outlet of the wind tunnel in order to obtain a free flow. Temperature and humidity were taken either from meteorological data of the site or local instruments placed upstream to reduce some of the uncertainty in the air density value. All the uncertainty estimates can be found in Appendix C. To compute the drag coefficient at normal incidence, $12$ different velocities have been used as shown in figure 2(a) where we plotted the force as a function of velocity for several screens. From the force curves, we can then deduce a drag coefficient $C_D$ defined as

(2.2)\begin{equation} C_D=\frac{F_D}{\dfrac12 \rho v_0^2 S_p}, \end{equation}

where $\rho$ is the fluid (air) density and $S_p$ the surface area of the screen. We note that here the surface $S_p$ is the total surface area of the screen and not its projected area along the far-field stream direction $\boldsymbol {v_0}$. Figure 2(b) shows the measured drag coefficient $C_D$ as a function of the velocity $v_0$ for different screens with different solidities $s$ and fibre diameter $d$, after subtraction of the contribution of the frame and the mast detailed in Appendix C, and rescaled by the value at the largest velocity of the wind tunnel $C_{D}(v_{0_{max}})$. In most experiments, the drag coefficient does almost not vary with the incoming fluid velocity, i.e. the drag force is proportional to the square of the fluid velocity $v_0$ upstream from the screen. This is consistent with the fact that the Reynolds number is always high in our experiments ($32\,000< Re<83\,000$). However, we observe that for small fibre radii and low solidities, the drag coefficient decreases with increasing velocity, suggesting viscous effects. At higher velocity, the drag coefficient tends toward a constant. For screens that show almost constant value of the drag coefficient with respect to the velocity, the data can then be fitted with a quadratic law to obtain the drag coefficient $C_{D}$. We calculated the coefficient of determination $R^2$ associated to the fitting in table 1. For these screens, with almost constant drag coefficient, we always have $R^2\geq 0.998$. For screens that show more important variation of the drag coefficient, we keep the value of the drag of each velocity between 5.0 and $13.0\,{\rm m}\,{\rm s}^{-1}$. Those values are reported in table 1. For comparison, we also calculated the coefficient of determination $R^2$ as if they were fitted by a quadratic law.

Figure 2. (a) Drag force $(F_D)$ measured for different square porous screens at normal incidence. (b) Drag coefficient $C_D$ as a function of velocity $v_0$, renormalised by the drag coefficient at maximum velocity for each screen $C_{D}(v_{0_{max}})$, including screens showing the higher variation of the drag coefficient with respect to the velocity. Colour bars: fibre diameter $d$.

Figure 3 shows the characteristic drag coefficient of the square screens as a function of the solidity $s$. The drag coefficient increases with increasing solidity, until it reaches a constant value $C_D\simeq 1$ at high solidities (for $s\leq 0.7$). This evolution of the drag coefficient is qualitatively consistent with previous experiments in particular Prandtl & Flachsbart (Reference Prandtl and Flachsbart1932), as well as with the 2-D model derived by Steiros & Hultmark (Reference Steiros and Hultmark2018). Furthermore, we observe that for a given solidity, the drag coefficient increases for decreasing $d$; this effect is particularly important at small solidities. In order to take into account the effect of the fibre diameter $d$, we define a local Reynolds number for the flow around each fibre denoted $Re_d$, for each screen, following

(2.3)\begin{equation} Re_d = \frac{v_0 d}{\nu} \end{equation}

as done in Prandtl & Flachsbart (Reference Prandtl and Flachsbart1932). At moderately high $Re_d$ (${\sim }10^2$), where viscous effects should not dominate the flow through the screen, the 2-D model overestimates the drag coefficient obtained experimentally, in particular at high solidities. This highlights the importance to take into account 3-D effects. At low $Re_d$, the model underestimates the drag coefficient, i.e. viscous effects must be included. We note that the model assumes a steady wake, and thus, is not applicable in the presence of vortex shedding. However, our experiments suggest that vortex shedding is negligible for our range of parameters, as further discussed in § 5.3.

Figure 3. Drag coefficient $C_{D_{max}}$ as a function of the solidity for various square porous screens normal to the free flow and comparison with the 2-D model of Steiros & Hultmark (Reference Steiros and Hultmark2018). The colour of the points indicates the value of the fibre diameter $d$ following the colour scale on the right. The data from Prandtl & Flachsbart (Reference Prandtl and Flachsbart1932) and the value at solidity $s=1$ from Blevins (Reference Blevins1992) are also plotted.

To characterize the influence of the solidity on the flow deviation, we measure the flow field with a PIV method. We use a wind tunnel with square test section of width 22 cm at constant velocity $2.84 \pm 0.02\,{\rm m}\,{\rm s}^{-1}$. The fluid is seeded with micro-droplets of water of diameter $3.0 \pm 2.0\,\mathrm {\mu } {\rm m}$. The Stokes number is defined as

(2.4)\begin{equation} St=\frac{1}{18} \frac{\rho_p d_p^2 v_0}{\mu L_{0}}, \end{equation}

where $\rho _p$ and $v_p$ are the particles density and diameter, respectively, and $\mu = \rho \nu$ is the viscosity of the fluid. The Stokes number is of the order of $10^{-3} \ll 1$ and we can consider that the water droplets act as passive tracers of the flow. A 1 mm-thick laser sheet is used to highlight the particles in a plane parallel to the flow. The laser (Elforlight LTD.model FCHPG-3000) has a wavelength of 532 nm and maximum power of 6.0 W. A high-speed camera PHOTRON was used at a frame rate of 4000 fps to record successive images that have been analysed with PIVlab in MATLAB (version 2.62).

In figure 4 we show the trajectories of the particles in a plane orthogonal to the screen at mid-height (i.e. at the middle of the screen), obtained by the superposition of the maximum intensity of $2000$ successive images. The upstream region is at the bottom of the figures and the downstream region at the top. We observe that the flow deviation around the screen increases as the solidity increases. We can further observe that the velocity decreases as the solidity increases, as shown by the variation of the length of the bright lines, shorter in figure 4(b) than in 4(a). Furthermore, in figure 4(c) we barely observe any particles crossing the screen, while some particles appear to be mixed by the recirculation in the wake.

Figure 4. Experimental evidence of the deviation of the flow and streamlines around a square porous screen at different solidities $s$ under normal uniform upstream flow $v_0 = 2.84\pm 0.02\,{\rm m}\,{\rm s}^{-1}$. Top view of the screen. The white bar (110 mm) corresponds to the screen embedded in the frame, the screen alone is 100 mm long. Results are shown for (a) P14, $s=0.24$; (b) P6, $s=0.61$; (c) P23, $s=0.82$.

In the following, we aim at developing a model that takes into account both 3-D and viscous effects. We then use this model to describe the flow around rectangular screens and predict the corresponding drag coefficient. In order to model these effects, we start by extending the 2-D model first proposed by Koo & James (Reference Koo and James1973) to a 3-D free flow and arbitrary shape, taking into account the base suction effect considered by Steiros & Hultmark (Reference Steiros and Hultmark2018). We then derive the equations to predict the drag coefficient following the method proposed by Steiros & Hultmark (Reference Steiros and Hultmark2018), and discuss the limits of this approach. In particular, while this approach takes into account some viscous effects in the pressure jump, it fails to describe the dependency of our measurements on $Re_d$. We thus incorporate in the model an empirical law as proposed by Brundrett (Reference Brundrett1993) for the pressure jump. Our modelling approach is sketched in figure 6.

In addition to these 2-D deviations, we observe 3-D deviations. In order to quantify them, we look at a plane in the diagonal direction (figure 5). In both the transverse and diagonal plane, we measure the proportion of deviated streamlines (i.e. the position of the separation distances $D_{sd}$ and $D_{st}$). If the deviations were purely two dimensional, we would have $D_{sd}=\sqrt {2} D_{st}$, i.e. the streamtube formed by the separatrix would be square like the screens. If 3-D deviations occurs, the streamtube should be deformed, and, in particular, $D_{sd}>D_{st}$. In our experiments, we found that indeed $D_{sd}/D_{st} = 1.7 \pm 0.1 >\sqrt {2}$.

Figure 5. Experimental evidence of the 3-D deviation of the flow and streamlines around a square porous screen (transverse and diagonal directions) for a screen of solidity $s=0.61$ (P6) under normal uniform upstream flow $v_0 = 2.84\pm 0.02\,{\rm m}\,{\rm s}^{-1}$. Here $D_{sd}/2 = 55.8 \pm 1.5$ mm and $D_{st}/2 = 32.8 \pm 1.2$ mm; thus, $D_{sd}/D_{st} = 1.7 \pm 0.1 >\sqrt {2}$.

Figure 6. Diagram of the theoretical model derived in this study and the ingredients from previous models used. In blue, a complete 3-D potential model. In grey, the theoretical basis used for the application to rectangular mesh screens. In purple, the application to rectangular mesh screens. Here $\phi$ denotes the velocity potential, $p_{III}$ denotes the base pressure and $p_{I}$$p_{II}$ denotes the pressure jump across the screen.

3. Three-dimensional model

In this section we derive the equations describing the 3-D flow around and through a porous screen with arbitrary shape and solidity. First, we adopt the method used by Koo & James (Reference Koo and James1973) that showed a good agreement with experimental results for 2-D flow in a channel except at high solidity. This discrepancy at high solidity may come from the lack of base pressure and vortex shedding in their theory, as suggested by Steiros & Hultmark (Reference Steiros and Hultmark2018) (see §§ 3.4 and 3.5). We extend the model to the 3-D case for free flow, i.e. a case where there are no boundaries constraining the flow, which is one of the major differences with the model of Koo & James (Reference Koo and James1973). We also take into account the effect of the base pressure. We first obtain a general formulation of the equations for a porous screen of arbitrary shape. We then apply our equations to the case of a rectangular plate inclined in a laminar flow for which an analytical solution can be found.

3.1. General formulation

The flow around bluff bodies is complex. In order to obtain an analytical or semi-analytical description of the flow, a widely used approach consists in simplifying the governing equations using potential flow theory outside the wake while introducing free parameters such as base pressure to account for viscous and complex phenomena in the wake and near the solid structure (Parkinson & Jandali Reference Parkinson and Jandali1970). The model proceeds with the same idea. The system is separated into four regions delimiting four flow regimes, as shown in figure 7. The flow is assumed to be stationary, incompressible and inviscid everywhere except through the porous structure where viscous effects cannot be neglected.

Figure 7. Diagram of the model for a 3-D potential flow around and through a porous screen. The dashed lines are the separation streamlines used as a boundary between the regions. The dotted lines are the streamlines used in the model to calculate the velocities in regions I and II. The incoming flow is laminar and is extended over the entire height of the system. Here $\mathcal {C}_w$ denotes the section of the wake.

In region I we assume that the flow is potential and the velocity is denoted $\boldsymbol {v}_I(x,y,z)$. Region I is located upstream of the structure as well as downstream outside the wake zone contained by a well-defined streamtube attached to the contour of the porous surface as shown in figure 7. Therefore, the velocity derives from a velocity potential denoted $\phi _I(x,y,z)$ for region I. Using the method employed by Koo & James (Reference Koo and James1973), and initially suggested by Taylor (Reference Taylor1944), we calculate the flow by modelling the screen with a continuous source distribution with strength $\varOmega (x_s,y_s,z_s)$, where $(x_s,y_s,z_s)$ denotes a point on the surface. We obtain the potential flow in region I by superposing the resulting potential flow from the distribution of sources with the uniform laminar flow $\boldsymbol {v}_0$. The two streamfunctions $\psi _I$ and $\chi _I$ needed to describe general 3-D incompressible flows can then be deduced using the following relation (for a definition of 3-D streamfunctions, see Yih Reference Yih1957):

(3.1)\begin{equation} \boldsymbol{v}_I(x,y,z) = \boldsymbol{\nabla} \psi_I(x,y,z) \wedge \boldsymbol{\nabla} \chi_I(x,y,z). \end{equation}

Regions II and III are located downstream of the porous surface in the near wake. In these regions the flow can be rotational so that we cannot use anymore a velocity potential to describe the flow. In region II the pressure and the velocity are not constant since they are influenced by the surface. However, in region III the flow is sufficiently far from the screen so that the streamlines tend to be aligned with the uniform flow $\boldsymbol {v}_0$ as represented by the contour $\mathcal {C}_w$ in figure 7. Therefore, the pressure tends towards a base suction pressure $p_{III}$ that is a priori different (and lower) than the constant external pressure $p_0$. This region is mathematically at infinity (there are no finite separations between regions II and III), however, since the flow aligns rapidly with the uniform flow, we indicate a region III in the near wake in figure 7. In two dimensions, the approach of constant pressure along separating streamlines has been successfully used in free-streamline theory by Wu (Reference Wu1962), Parkinson & Jandali (Reference Parkinson and Jandali1970) and Roshko (Reference Roshko1954) to model the wake. In this model, since we consider three dimensions, we can not adopt free-streamline theory, but $p_{III}$ can be considered as having the same role as the constant pressure used in such theory. The flow in regions II and III is found with matching conditions as explained later.

In figure 7 we added a region IV that is located in the far wake where the mixing with the outer flow can not be ignored. In this region, the pressure should increase to reach again the pressure $p_0$ outside the wake. We assume that this region has little influence on the flow near the porous screen and on the aerodynamic forces, and therefore, it is not included in the model. Consequently, in our model the pressure in the far wake will remain equal to $p_{III}$.

3.2. Determination of the flow in region I

In region I the flow is potential, and we derive the velocity from the velocity potential. For simplicity, we set the reference frame so that the axis $(Oz)$ is aligned with the velocity $\boldsymbol {v}_0$ without loss of generality. Due to the linearity of the Laplacian, we first consider the potential flow $\phi (x,y,z)$ for the source distribution only, then we add the potential flow for the uniform flow. The velocity potential from the source distribution $\varOmega$ located on a general regular surface $\mathcal {S}_p$ is the solution of the equation

(3.2)\begin{equation} \Delta \phi(x,y,z) = \varOmega(x,y,z)\mathbb{1}_{\mathcal{S}_p}, \end{equation}

where $\mathds {1}_{\mathcal {S}_p}$ denotes the Dirac function associated to the surface $\mathcal {S}_p$. For a point source in three dimensions centred at the origin, the Green function of the Laplacian is

(3.3)\begin{equation} \varGamma(x,y,z) ={-} \frac{1}{4{\rm \pi}}\frac{1}{\sqrt{x^2+y^2+z^2}}. \end{equation}

Therefore, if we assume $\boldsymbol {\xi }: U \subset \mathbb {R}^2\to \mathbb {R}^3$ to be a surface patch of a general regular surface $\mathcal {S}_p$ with coordinates

(3.4)\begin{equation} \boldsymbol{\xi} =\begin{pmatrix} x_s(u,v)\\ y_s(u,v)\\ z_s(u,v) \end{pmatrix}, \end{equation}

parametrized by two parameters $u$ and $v$, with $(u,v)\in (U = [a,b]\times [c,d])$ with $(a,b,c,d)\in \mathbb {R}^4$, then the velocity potential $\phi (x,y,z)$ is expressed for all $(x,y,z)\in \mathbb {R}^3\setminus \mathcal {S}_p$ as (Pressley Reference Pressley2010)

(3.5)\begin{align} \phi(x,y,z) = \iint_{U} \varOmega(u,v)\varGamma(x-x_s(u,v),y-y_s(u,v), z-z_s(u,v)) \left\|\frac{\partial \boldsymbol \xi}{\partial u} \wedge \frac{\partial \boldsymbol \xi}{\partial v}\right\| \mathrm{d}u \,\mathrm{d}v, \end{align}

and the total velocity potential can be written as the sum

(3.6)\begin{equation} \phi_I(x,y,z) = {v_0} z + \phi(x,y,z). \end{equation}

We deduce the velocity in region I with $\boldsymbol {v_I}(x,y,z)=\boldsymbol {grad}(\phi _I(x,y,z))$. Note that the Green function can be changed without other modifications in the model to study the situation of a flow in a confined environment or near a wall.

3.3. Determination of the flow in regions II and III

In region II the flow can be rotational, and therefore, is not necessarily potential. The flow is obtained from the streamfunctions by considering, as done by Koo & James (Reference Koo and James1973), that the streamlines in region II have the same pattern as if they were obtained by the streamfunctions from the superposition of the distribution of sources and the uniform flow $\boldsymbol {v}_0$. This can be formulated in a general way by writing the two streamfunctions for the flow in region II as functions of the streamfunctions of region I. Let $\psi _{II}$ and $\chi _{II}$ be the streamfunctions in region II. As defined above, $\psi _I$ and $\chi _{I}$ are the streamfunctions deduced from the flow in region I, functions that can be considered in the whole space. Then, without loss of generality, we choose an analog formulation as the one proposed by Koo & James (Reference Koo and James1973) to describe the flow in region II, using the functions $f$ and $g$:

(3.7a,b)\begin{equation} \psi_{II}=f(\psi_I)\psi_I \quad {\rm and}\quad \chi_{II}=g(\chi_I)\chi_I . \end{equation}

This means that the velocity $\boldsymbol {v}_{II}$ in region II, and the velocity $\boldsymbol {v}_{I}$ obtained from the velocity potential $\phi _I(x,y,z)$, are co-linear at any point. We indeed obtain

(3.8)\begin{equation} {\boldsymbol{v}_{II}}(x,y,z) = \left(\frac{\mathrm{d}f}{\mathrm{d}\psi_I}+f(\psi_I)\right) \left(\frac{\mathrm{d}g}{\mathrm{d}\chi_I}+g(\chi_I)\right){\boldsymbol{v}_I}(x,y,z). \end{equation}

We define the attenuation function $E$ as

(3.9)\begin{equation} E(\psi_I,\chi_I)=\left(\frac{\mathrm{d}f}{\mathrm{d}\psi_I}+f(\psi_I)\right) \left(\frac{\mathrm{d}g}{\mathrm{d}\chi_I}+g(\chi_I)\right), \end{equation}

which is the crucial quantity to determine the flow in region II. Since $E$ is a function only of the streamfunctions, it is constant along a streamline and is therefore entirely defined by considering its value on the screen.

In addition to (3.8), the mass flow rate must be conserved when the fluid passes through the screen implying the continuity of the normal velocity at the screen between region I and region II:

(3.10)\begin{equation} {v_{I}}_n(x_s,y_s,z_s) = {v_{II}}_n(x_s,y_s,z_s). \end{equation}

At this point, our system thus contains two unknowns that are the attenuation function $E$ and the distribution of sources $\varOmega$. We have one equation (3.10), and another equation linking the velocities and pressures in the vicinity of the porous structure is required to close our system of equations. For this purpose, two streamlines are considered as shown in figure 7: $(AP)$ and $(PB)$ where the point $P$ is on the screen taken as a surface from a macroscopic point of view. Along each of these streamlines, Bernoulli's equation can be applied, and we thus obtain, for $(AP)$,

(3.11)\begin{equation} \tfrac{1}{2}\rho v_{I}^2(x_A,y_A,z_A)+p_{I}(x_A,y_A,z_A) = \tfrac{1}{2}\rho v_{I}^2(x_s,y_s,z_s)+p_{I}(x_s,y_s,z_s), \end{equation}

and, for $(PB)$,

(3.12)\begin{equation} \tfrac{1}{2}\rho v_{II}^2(x_s,y_s,z_s)+p_{II}(x_s,y_s,z_s) = \tfrac{1}{2}\rho v_{II}^2(x_B,y_B,z_B)+p_{II}(x_B,y_B,z_B). \end{equation}

We consider that the points $A$ and $B$ are far enough from the screen so that we can take constant values of the velocities and pressures (see Fail, Lawford & Eyre (Reference Fail, Lawford and Eyre1957) for flat plates normal to an air stream). Therefore, upstream we have $\boldsymbol {v}_{I}(x_A,y_A,z_A)= \boldsymbol {v}_0$ and $p_{I}(x_A,y_A,z_A)=p_0$; downstream we take the mean value of the velocity over a section of the wake orthogonal to the far-field stream direction ($\boldsymbol {v}_0$): $\boldsymbol {v}_{II}(x_B,y_B,z_B) = \overline {\lim _{z \rightarrow +\infty }\boldsymbol {v}_{II}(x,y,z)} = \overline {E(\psi _{I},\chi _{I})}\boldsymbol {v}_0$. In the rest of the paper, $\overline {E(\psi _{I},\chi _{I})}$ will be denoted as $\bar {E}$. Koo & James (Reference Koo and James1973) considered a far-downstream constant pressure $p_{0}$; however, it is known that the pressure in the wake is lower than the pressure outside the wake that contributes to aerodynamic forces. Steiros & Hultmark (Reference Steiros and Hultmark2018) therefore introduced a suction base pressure $p_{III}$ at the point $B$, which is assumed to be constant far enough from the screen in region III as explained above. Thus, we introduce a third free parameter $p_{III}$ that will be determined from conservation equations in § 3.5. Note that since we assume that the pressure $p_{III}$ is constant, the pressure will be discontinuous across the wake boundaries (as in the model of Koo & James (Reference Koo and James1973) in the 2-D case). By combining (3.11) and (3.12), decomposing the velocities according the tangential and normal components on the screen (respectively, ${v_I}_t$ and ${v_{I}}_n$) and by using (3.10), we obtain the pressure difference

(3.13)\begin{align} p_0 - p_{III} &= \tfrac{1}{2}\rho(1-E^2(\psi_{I},\chi_{I})){v_I^2}_t(x_s,y_s,z_s) + \tfrac{1}{2}\rho(\bar{E}^2-1) v_0^2 \nonumber\\ &\quad + p_I(x_s,y_s,z_s) - p_{II}(x_s,y_s,z_s). \end{align}

If we are able to determine the pressure differences $p_0 - p_{III}$ and $p_I(x_s,y_s,z_s) - p_{II}(x_s,y_s,z_s)$ independently from these equations, then we can use (3.10) to find the attenuation function $E$ and (3.13) to find the distribution of sources $\varOmega$. In order to determine the suction pressure $p_{III}$, we follow the method of Steiros & Hultmark (Reference Steiros and Hultmark2018) in § 3.5, using the conservation law of the momentum in the control volume $V_1$ shown in figure 7, and at the vicinity of the screen. Before addressing this problem, we focus in the following § 3.4 on the pressure jump $p_I(x_s,y_s,z_s) - p_{II}(x_s,y_s,z_s)$.

3.4. Pressure jump across the screen

A summary of the models for the relation between the pressure jump and screen porosity can be found in Xu et al. (Reference Xu, Patruno, Lo and de Miranda2020). This problem has largely been discussed in many papers, raising the difficulties of a general formulation. In their model, Steiros & Hultmark (Reference Steiros and Hultmark2018) consider two streamlines passing through the screen in a hole where the velocity and the pressure are assumed to be uniform. Immediately upstream after the acceleration of the fluid, the characteristic pressure of the flow is denoted $p_h$ and is assumed to correspond physically to the mean pressure inside the hole. As noted by several authors including Taylor & Davies (Reference Taylor and Davies1944) and Wieghardt (Reference Wieghardt1953), the characteristic velocity immediately upstream should be regarded as the mean velocity after contraction of the flow within the holes, ${v_h}_1$, denoting the average velocity through the screen expressed as

(3.14)\begin{equation} {v_h}_1 = \frac{{v_I}_n(x_s, y_s, z_s)}{(1-s)}, \end{equation}

The velocity accelerates to ${v_h}_1$ in the hole and the pressure decreases to $p_h$ so that at this point there are no losses. Along this first streamline, Bernoulli's equation leads to

(3.15)\begin{equation} p_{I}(x_s, y_s,z_s) + \tfrac{1}{2}\rho {v_I^2}(x_s, y_s,z_s) = p_h + \tfrac{1}{2}\rho {v_{h}^2}_1. \end{equation}

Immediately downstream, considering also a homogenized velocity, the flow enlarges and the velocity reaches a value equal to

(3.16)\begin{equation} {v_h}_2 = {v_{II}}_n(x_s, y_s, z_s), \end{equation}

taken as the characteristic velocity just after the hole. As discussed by Taylor & Davies (Reference Taylor and Davies1944) and Steiros & Hultmark (Reference Steiros and Hultmark2018), there are pressure losses in the pores, i.e. $p_{II}< p_{I}$. Taylor & Davies (Reference Taylor and Davies1944) consider that a fraction $\lambda$ of the pressure lost in acquiring the velocity $v_h$ is regained when the stream becomes uniform again behind the sheet. Steiros & Hultmark (Reference Steiros and Hultmark2018) consider that the pressure loss during the fluid acceleration is not recovered, i.e. all surplus kinetic energy due to fluid acceleration is lost and not reconverted into pressure (i.e. $\lambda =0$). In that case, Bernoulli's equation along this second streamline leads to

(3.17)\begin{equation} p_h + \tfrac{1}{2}\rho {v_{h}^2}_2 = p_{II}(x_s, y_s, z_s) + \tfrac{1}{2}\rho {v_{II}^2}(x_s, y_s,z_s). \end{equation}

We combine (3.8), (3.10) and (3.14)–(3.17) to obtain

(3.18)\begin{align} {p_{II}}(x_s,y_s,z_s)-{p_{I}}(x_s,y_s,z_s) &= \tfrac{1}{2}\rho {v_I}_t^2(x_s, y_s,z_s) (1-E^2(\psi_I,\chi_I)) \nonumber\\ &\quad + \tfrac{1}{2}\rho{v_I}_n^2(x_s, y_s,z_s)\theta(s), \end{align}

with

(3.19)\begin{equation} \theta(s) = 1-\frac{1}{(1-s)^2}. \end{equation}

Injecting the pressure difference $p_{I}$$p_{II}$ obtained in (3.18) in (3.13), we have

(3.20)\begin{equation} p_0\unicode{x2013}p_{III} ={-}\tfrac{1}{2}\rho{v_I^2}_n(x_s,y_s,z_s)\theta(s) + \tfrac{1}{2}\rho(\bar{E}^2-1) v_0^2. \end{equation}

We assumed that $p_{III}$ is constant, therefore, the right-hand side of the above equation has to be constant, which leads to the condition

(3.21)\begin{equation} \boldsymbol{grad}_{\mathcal{S}_p}({v_I}_n^2(x_s,y_s,z_s)\theta(s))= \boldsymbol{0}, \end{equation}

with $\theta (s)$ that can vary for non-homogeneous porous surfaces. Therefore, under the assumptions made so far, we are looking for a source strength $\varOmega$ that satisfies this condition. It is possible to relax certain restrictions on the value of $\varOmega$ by considering a variable base suction pressure $p_{III}$ or a variable far wake velocity $\boldsymbol {v}_{II}(x_B,y_B,z_B) = E(\psi _{I},\chi _{I})\boldsymbol {v}_0$. In that case, the problem is far more difficult, and should be solved numerically. Such a resolution is beyond the scope of the present work.

3.5. Drag coefficient

At this point, our system of equations is not closed. Therefore, in this section we find another equation by adapting the model of Steiros & Hultmark (Reference Steiros and Hultmark2018) in three dimensions with a screen of arbitrary shape and orientation in the flow. It uses two expressions of the aerodynamic forces that can be combined to give a final equation, closing eventually our system of equations.

A first expression is given by the momentum balance around the surface of the screen, using the drag coefficient as defined in (2.2), giving

(3.22)\begin{equation} C_D = \frac{1}{\dfrac{1}{2}\rho S_{p} v_0^2} \left(\iint_{\mathcal{S}_p}(p_I\unicode{x2013}p_{II})\boldsymbol{e}_{z}\boldsymbol{\cdot} \boldsymbol{n}_s\,\mathrm{d}S+\rho\iint_{\mathcal{S}_p} {v_I}_n \boldsymbol{e}_{z}\boldsymbol{\cdot} (\boldsymbol{v}_{I}\unicode{x2013}\boldsymbol{v}_{II})\,\mathrm{d}S\right), \end{equation}

where the normal vector $\boldsymbol {n}_s$ points in the direction of region II. The pressure difference $p_I\unicode{x2013}p_{II}$ in (3.22) is given by (3.18).

The second expression of the drag coefficient is provided by the momentum balance in a control volume $V_1$ around the screen as shown in figure 7. We consider that this volume $V_1$ is large enough so that the velocities at the surfaces $\mathcal {S}_x$ and $\mathcal {S}_y$ on the sides of the block parallel to the $z$ axis are equal to $\boldsymbol {v}_0 + \boldsymbol {v}_{\epsilon }$ where $v_{\epsilon } \ll v_0$. With this approximation, the projection of momentum balance equation onto the far-field stream direction ($\boldsymbol {v}_0$) gives an expression of the drag $F_D$,

(3.23)\begin{equation} F_D =\rho S_w (1-\bar{E}^2){v^2_0} + S_w (p_0\unicode{x2013}p_{III}) - \rho v_0 \iint_{\mathcal{S}_x + \mathcal{S}_y} \boldsymbol{v}_I \boldsymbol{\cdot} \boldsymbol{n}\,\mathrm{d}S, \end{equation}

where $S_w$ is the area of the section of the wake (orthogonal to the $z$ axis), and the bar over $E$ denotes the mean over the considered surface. The normal vector $\boldsymbol {n}$ points outwards from the control volume. Then, a mass balance in the same control volume gives

(3.24)\begin{equation} \rho v_{0} S_w = \rho \iint_{\mathcal{S}_x + \mathcal{S}_y} \boldsymbol{v}_I \boldsymbol{\cdot} \boldsymbol{n}\,\mathrm{d}S + \rho \bar{E}v_{0}S_w. \end{equation}

This equation allows us to find the value of the last term of the conservation of momentum equation (3.23). Moreover, the section of the wake $S_w$ is determined with a mass balance through the screen

(3.25)\begin{equation} \rho \iint_{\mathcal{S}_p} \boldsymbol{v}_I \boldsymbol{\cdot} \boldsymbol{n}\,\mathrm{d}S = \rho \bar{E}v_{{0}} S_w. \end{equation}

Equations (3.23)–(3.25) are now combined to obtain a second expression of the drag coefficient, which reads

(3.26)\begin{equation} C_D = \frac{1}{\dfrac{1}{2}\rho S_p v_0^2}(p_0\unicode{x2013}p_{III} + \rho v_0^2 (\bar{E}-\bar{E}^2))\frac{1}{v_0 \bar{E}}\iint_{\mathcal{S}_p}{v}_{I_n}\,\mathrm{d}S. \end{equation}

Equation (3.20) giving an expression of $p_0\unicode{x2013}p_{III}$ is used so that the second expression of the drag (3.26) is calculated further:

(3.27)\begin{equation} C_D = \frac{1}{S_p}(-(1-\bar{E})^2 - {v}_{I_n}^2\theta(s)) \frac{1}{v_0 \bar{E}}\iint_{\mathcal{S}_p}{v}_{I_n} \,\mathrm{d}S. \end{equation}

The final equation of the problem is obtained by equalizing the two expressions of the drag coefficient (3.22) and (3.27), yielding

(3.28)\begin{align} & \bar{E}\iint_{\mathcal{S}_p} ({\tilde{v}_{I_t}^2}(1-E^2) + {\tilde{v}_{I_n}^2}\theta(s)) \boldsymbol{n}_s\boldsymbol{\cdot}\boldsymbol{e}_z + 2\tilde{v}_{I_n}(1-E){\boldsymbol{\tilde{v}}_{I_t}}\boldsymbol{\cdot} \boldsymbol{e}_z \,\mathrm{d}S \nonumber\\ &\quad =(-(1-\bar{E})^2 - \tilde{v}_{I_n}^2\theta(s)) \iint_{\mathcal{S}_p}\tilde{v}_{I_n}\,\mathrm{d}S, \end{align}

with dimensionless velocities $\tilde {v}={v}/{v_0}$.

Finally, we give the explicit expression for $v_{I_n}$ and $E$. Finding $v_{I_n}$ leads to a boundary surface potential problem. An explicit expression of $v_{I_n}$ is found from the gradient of the velocity potential (3.6) having therefore an integral term. This integral term is known as a harmonic double-layer potential with density $\varOmega$ (Gunter Reference Gunter1967), which is defined on a subdomain of $\mathbb {R}^3\setminus \mathcal {S}_p$. This integral term becomes singular if it is evaluated on the surface $\mathcal {S}_p$, however, it can be continuously extended on the surface for each side and the value depends on the side by which we approach the surface. From the definitions of regions I and II, for $\boldsymbol {m}=(x,y,z)\in \mathcal {S}_p$, we have

(3.29)\begin{equation} v_{I_n}(\boldsymbol{m}) = v_0\boldsymbol{e_z}\boldsymbol{\cdot}\boldsymbol{n}_s + v_n^-(\boldsymbol{m}), \end{equation}

and

(3.30)\begin{equation} v_{{II}_n}(\boldsymbol{m}) = E(\psi_I,\chi_I)(v_0\boldsymbol{e_z}\boldsymbol{\cdot}\boldsymbol{n}_s + v_n^+(\boldsymbol{m})), \end{equation}

with

(3.31)\begin{equation} v_n^{{\pm}}(\boldsymbol{m}) = \lim_{\epsilon \to 0^+}\frac{\partial \phi}{\partial \boldsymbol{n}_s} (\boldsymbol{m} \pm \epsilon \boldsymbol{n}_s), \end{equation}

with $\phi$ defined in (3.5) and the directional outward normal derivative

(3.32)\begin{equation} \frac{\partial \phi}{\partial \boldsymbol{n}_s}(\boldsymbol{m}) = \boldsymbol{grad}(\phi(\boldsymbol{m}))\boldsymbol{\cdot} \boldsymbol{n}_s. \end{equation}

From a theorem that can be found, for instance, in Kress (Reference Kress1999, p. 80), the value of the limits above can be expressed using an improper integral. For $\boldsymbol {m}=(x,y,z)\in \mathcal {S}_p$, this expression reads

(3.33)\begin{equation} v_{n}^{{\pm}}(\boldsymbol{m}) = \iint_{U} \varOmega(u,v) \frac{\partial \varGamma}{\partial \boldsymbol{n}_s} (\boldsymbol{m}-\boldsymbol{m}_s(u,v)) \left\|\frac{\partial \boldsymbol \xi}{\partial u} \wedge \frac{\partial \boldsymbol \xi}{\partial v}\right\| \mathrm{d}u\,\mathrm{d}v \pm \frac{1}{2}\varOmega(\boldsymbol{m}), \end{equation}

with $\boldsymbol {m}_s(u,v)=(x_s(u,v),y_s(u,v),z_s(u,v))$. Finally, $E$ is found by applying the continuity equation (3.10) for the normal velocities across the surface, with the expressions given in (3.29) and (3.30). For $\boldsymbol {m}=(x,y,z)\in \mathcal {S}_p$, we have

(3.34)\begin{equation} E(\boldsymbol{m}) = \frac{v_0\boldsymbol{e_z}\boldsymbol{\cdot}\boldsymbol{n}_s + v_n^-(\boldsymbol{m})}{v_0\boldsymbol{e_z}\boldsymbol{\cdot}\boldsymbol{n}_s + v_n^+(\boldsymbol{m})}. \end{equation}

The system of equations is now closed and enables us to obtain the velocities and pressure at any location in the flow by determining the source strength $\varOmega$, solution of (3.28). In the next section we solve the problem analytically in a basic but very common geometry.

3.6. Application to an inclined rectangular screen in free flow

For a rectangular geometry and homogeneous solidity, it is possible to obtain without major difficulty an analytical solution of the equations of our model. We therefore apply the 3-D model to the simple case of a rectangular screen (centred at $z=0$ as shown in the figure 8) in a free laminar flow in order to find the flow and the drag coefficient as a function of the solidity. In that case, the velocity potential (3.6) becomes

(3.35)\begin{align} \phi_I(x,y,z) = v_0z + c - \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{\varOmega(u,v)\,\mathrm{d}u\,\mathrm{d}v}{\sqrt{(x-v)^2+ (y-u\sin{(\beta)})^2+(z-u\cos{(\beta)})^2}}. \end{align}

Figure 8. Diagram of the inclined rectangular porous screen in the $(Oyz)$ plane in a 3-D free flow with an angle $\beta$. The dashed lines are the separation streamlines used as a boundary between the regions of the model. The dotted line between points A and B is the streamline used in the model to calculate the velocities in the regions I and III. The incoming flow is laminar.

We then calculate the velocities at the screen. The objective of the calculus is to obtain $\varOmega$ as a function of the solidity, which will determine all the other variables of the problem.

The normal component of the velocity in region I at the arbitrary position $(w,t)$ on the surface is

(3.36)\begin{equation} {v_{I}}_n^{{\pm}}(w,t) = v_0\sin{(\beta)} + \lim_{\epsilon \to 0^\pm}\frac{1}{4{\rm \pi}} \iint_{\mathcal{S}_p} f_{\epsilon}(u,v) \varOmega(u,v)\,\mathrm{d}u\,\mathrm{d}v, \end{equation}

with

(3.37)\begin{equation} f_{\epsilon}(u,v) = \frac{\epsilon\cos{(\beta)}\sin{(\beta)}}{((t-v)^2+(w-u)^2+(2(w-u) + \epsilon) \epsilon\cos{(\beta)})^{{3}/{2}}}. \end{equation}

The details of the calculations are given in Appendix B. At the screen, depending on the direction from which we approach the screen ($\epsilon \to 0^{\pm }$), the magnitude of the normal component of the velocity is constant and is equal to

(3.38)$$\begin{gather} {v_{I}}_n(w,t) = {v}_n^{-}(w,t) = v_0\sin{(\beta)} - \tfrac{1}{2} \varOmega(w,t), \end{gather}$$
(3.39)$$\begin{gather}{v_{II}}_n(w,t) = E(\psi_I,\chi_I){v}_n^+(w,t) = E(\psi_I,\chi_I)(v_0\sin{(\beta)} + \tfrac{1}{2} \varOmega(w,t)). \end{gather}$$

At this point, for an homogeneous screen ($s$ constant on the surface), (3.21) leads to $\boldsymbol {grad}(\varOmega (w,t))=\boldsymbol {0}$, i.e. the source strength is a constant.

Thus, (3.10) leads to a constant attenuation function

(3.40)\begin{equation} E = \frac{v_0\sin{(\beta)} - \dfrac{1}{2} \varOmega}{v_0\sin{(\beta)} + \dfrac{1}{2} \varOmega}. \end{equation}

If we found constant normal velocity, it is not the case of the tangential velocity for which the magnitude varies on the surface of the screen. The tangential component of the velocity at the surface is

(3.41)\begin{align} {v_I}_t(x,w) &= (\varOmega^2 \mathcal{I}_x^2(x,w) + (\varOmega(\sin{(\beta)}\mathcal{I}_y(x,w)+\cos{(\beta)}\mathcal{I}_z(x,w)) \nonumber\\ &\quad + \cos{(\beta)}v_0)^2)^{{1}/{2}}, \end{align}

with $\mathcal {I}_x$, $\mathcal {I}_y$ and $\mathcal {I}_z$ the following surface integrals that are calculated in Appendix A:

(3.42)$$\begin{gather} \mathcal{I}_x(x,w) = \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{x-v}{((x-v)^2+(w-u)^2)^{{3}/{2}}} \mathrm{d}u \,\mathrm{d}v , \end{gather}$$
(3.43)$$\begin{gather}\mathcal{I}_y(x,w) = \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{(w-u)\sin{(\beta)}}{((x-v)^2+(w-u)^2)^{{3}/{2}}} \mathrm{d}u \,\mathrm{d}v, \end{gather}$$
(3.44)$$\begin{gather}\mathcal{I}_z(x,w) = \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{(w-u)\cos{(\beta)}}{((x-v)^2+(w-u)^2)^{{3}/{2}}} \mathrm{d}u \,\mathrm{d}v. \end{gather}$$

For the sake of simplicity, we approximate the magnitude of the tangential component of the velocity as its root mean square. We note that since in (3.28), $v_{I_t}$ appears both in linear and quadratic form, this approximation becomes exact for a surface normal to the mean flow direction ($\beta ={{\rm \pi} }/{2}$). We thus take

(3.45)\begin{equation} {v_I}_t = (\varOmega^2 \gamma_0 + v_0^2\cos^2{(\beta)})^{{1}/{2}}, \end{equation}

where $\gamma _0$ can be considered as a shape factor (see Appendix A). Equation (3.18) becomes

(3.46)\begin{equation} {p_{II}}-{p_{I}} = \tfrac{1}{2}\rho ({v_I}_t^2(1-E^2)+{v_I}_n^2\theta(s)). \end{equation}

Equation (3.13) leads to

(3.47)\begin{align} {p_{III}}-{p_{0}} &= \tfrac{1}{2}\rho (1-E^2)(v^2_0-{v_I}_t^2) + p_{II} - p_{I} \nonumber\\ &= \tfrac{1}{2}\rho((1-E^2)v^2_0 + {v_I}_n^2\theta(s)). \end{align}

Using (3.23)–(3.25), we obtain a first expression of the magnitude of the drag force $F_D$:

(3.48)\begin{equation} F_D = \rho v_0(1-E){v_n} {S_p} + \frac{1}{v_0}(p_0 - p_{III}) \frac{{v_n}}{E} {S_p}. \end{equation}

The second expression for the drag force $F_D$ is obtained with (3.22),

(3.49)\begin{equation} F_D = ({p_I}\unicode{x2013}{p_{II}})\sin{(\beta)}{S_p} + \rho {v_n} v_0\cos^2(\beta)(1-E){S_p}. \end{equation}

Note that, for a rectangular screen orthogonal to the free flow, the second term of (3.49) vanishes since $\beta ={{\rm \pi} }/{2}$ and we obtain a drag coefficient proportional to the pressure difference ${p_I}\unicode{x2013}{p_{II}}$. Now, by denoting $\omega = {\varOmega }/{v_0}$ and combining (3.48) and (3.49), we obtain the following equation that we have to solve to find the value of the source strength:

(3.50)\begin{equation} -\tfrac{1}{8} \omega^4\theta(s) + \omega^2 \sin^2{(\beta)}(8\gamma_0 + \theta(s) - 2) -4\omega\sin{(\beta)} - 2\sin^4{(\beta)}\theta(s) = 0. \end{equation}

This equation has been solved using Python 3.9.5 and the method fsolve from scipy.optimize. The solution for a square porous plate at normal incidence with $\beta ={{\rm \pi} }/{2}$ is plotted in figure 9(a). At low solidity, our model for a square plate is close to the prediction of Steiros & Hultmark (Reference Steiros and Hultmark2018) and Taylor & Davies (Reference Taylor and Davies1944). Above a solidity $s=0.4$ the drag coefficient becomes slightly different from the prediction of Steiros & Hultmark (Reference Steiros and Hultmark2018) and the difference increases with increasing solidity. Three-dimensional effects are therefore important at high solidity.

Figure 9. Theoretical prediction of the drag coefficient and velocities on the surface as a function of the solidity obtained with the 3-D model applied to a square homogeneous screen normal to the flow ($\beta = {{\rm \pi} }/{2}$), and comparison with different 2-D models. The full lines are obtained with the 3-D model developed in this paper using the same pressure jump law across the screen as Steiros & Hultmark (Reference Steiros and Hultmark2018) (but taking account of geometric 3-D effects). (a) Drag coefficient. (b) Normal velocity and quadratic mean of the tangential velocity as defined respectively in (3.38) and (3.45). (a) Theoretical drag coefficient. (b) Theoretical velocity at the screen.

At high solidity ($s\lesssim 1$), the curve converges towards the drag coefficient of a flat solid plate. This value is predicted to be $1.2$ in our model, which is close to the experimental measurement at the global Reynolds number $Re \approx 10^4$, giving approximately $1.05$ according to Blevins (Reference Blevins1992), $1.17$ according to the synthesis of Hoerner (Reference Hoerner1965) on various drag measurements and $0.939$ in pour experiments.

Our model takes into account the shape of the porous surface through the parameter $\gamma _0$. For a rectangular screen, the aspect ratio, taken into account in $\gamma _0$, has an influence on the result as detailed in Appendix A. According to Hoerner (Reference Hoerner1965), for a flat plate normal to the flow, the drag coefficient increases very slowly when the aspect ratio is reduced until a ratio of approximately $0.1$. Beyond this point the increase becomes more pronounced, until it reaches $C_D=2.0$ for an infinitely thin plate ($1.90$ according to Blevins Reference Blevins1992). In our model, we indeed observe an increase of the drag when the aspect ratio decreases (see Appendix A) from $C_D=1.2$ for a square plate (with $\gamma _0=0.0998$), 1.29 for a rectangle with aspect ratio 1/10 (with $\gamma _0=0884$) to $C_D=1.33$ for an infinitely thin plate. This value for an infinitely thin plate is lower than the experimental value, that is, $C_D \approx 2.0$ according to Blevins (Reference Blevins1992) and Hoerner (Reference Hoerner1965). This difference may be due, in this case of very high aspect ratio, to the vortex shedding that is not taken into account in the present wake model, as noted by Steiros & Hultmark (Reference Steiros and Hultmark2018).

In figure 9(b) we compare the normal and tangential velocities in the 2-D and 3-D cases. The tangential velocity is taken in both cases as the quadratic mean over the whole surface, in three dimensions it is defined in (3.45), the normal velocity is constant on the surface in both cases also, in three dimensions it is defined in (3.38). While the normal velocity is the same in two dimensions and three dimensions, we see an influence of 3-D effects on the tangential velocity.

While taking into account the 3-D effects improve the prediction of $C_D$ at high solidities, it does not improve the prediction at moderate and low solidities. The limitations of potential flow theory to describe the complex flow around bluff bodies are arising at high solidity, but should be negligible at low and moderate solidity. Here, we hypothesize that the main reason for the discrepancy between our experimental results (figure 3) and the theoretical prediction (figure 9a) lies in the pressure loss at the pore scale. Indeed, the relation used to estimate the pressure jump $p_{I}\unicode{x2013}p_{II}$ (3.18) is a strong assumption, and does not consider the pressure losses by viscous friction. In particular, it does not explicitly account for the viscosity, and can not account for the variations of the drag coefficient with the local Reynolds number observed in the experiments.

4. Local viscous effects

As we have seen, the pressure losses are accounted for in (3.18) with the parameter $\theta (s)$ given by (3.19). We have adopted the assumption of Steiros & Hultmark (Reference Steiros and Hultmark2018), i.e. that all energy due to fluid acceleration is lost and not reconverted into pressure, which is a crucial assumption. The approach of Taylor & Davies (Reference Taylor and Davies1944) is to assume that a fraction $\lambda$ of energy is regained, which would modify $\theta (s)$ to read

(4.1)\begin{equation} \theta^{*}(s) = 1-\frac{1-\lambda}{(1-s)^2}. \end{equation}

In their experiments, Taylor & Davies (Reference Taylor and Davies1944) estimate $\lambda \simeq 0.4$. Taking into account this effect would decrease the drag coefficient, i.e. bring it closer to our experimental values at high $Re_d$. However, $\lambda$ can only be obtained empirically. Furthermore, in our experiments, $Re_d$ is smaller than the critical value of ${\approx }1000$ above which viscous effects can be neglected (Hoerner Reference Hoerner1952). We thus want to account for the pressure losses due to viscous friction as well. One possible approach is to tune the pressure loss term in (3.18) with a single prefactor $\theta (s) f$ that should depend on $Re_d$; it should increase as $Re_d$ decreases, accounting for increased viscous dissipation. This prefactor $\theta (s) f(Re_d)$ has been determined experimentally for wire screens by Brundrett (Reference Brundrett1993) and Bailey et al. (Reference Bailey, Montero, Parra, Robertson, Baeza and Kamaruddin2003). In this section we follow their approach to include viscous effects and compare our predictions with experimental results.

4.1. Pressure jump dependency on the local geometry of the pores and viscous effects

As explained above, in the equation of the pressure jump (3.18) with (3.19), as formulated by Steiros & Hultmark (Reference Steiros and Hultmark2018), we do not take into account the dependency of the pressure drop on the local Reynolds number $Re_d$, computed at the scale of the screen pores (and thus, much smaller than $Re$), the geometry of the holes and other possible dependency like the energy transfer between the material of the screen and the fluid. For instance, Ando et al. (Reference Ando, Nishikawa, Kaneda and Suga2022) showed that a layer of flexible fibres can have a higher permeability than the same layer of rigid fibres due to a flow-induced deformation. Moreover, as shown by Schubauer, Spangenberg & Klebanoff (Reference Schubauer, Spangenberg and Klebanoff1950), the angle of the screen relative to the laminar upstream free flow has an impact on the pressure drop. Kalugin et al. (Reference Kalugin, Epikhin, Chernukha and Kalugina2021) explained also that for inclined perforated plates, the structure of the flow in the holes depends on (1) the distance of the hole on the plate from the leading edge, and (2) the angle of attack of the plate. For a low angle of attack, the effective hole area can be significantly reduced due to the difficulty of the flow to deflect from its original direction mostly parallel to the surface of the plate. All these studies underline the current difficulty to obtain a general formulation of the pressure loss for an arbitrary porous screen. Therefore, in what follows, we adopt another method based on empirical laws in order to test whether this would be sufficient to estimate the drag accurately.

Numerous experimental investigations have shown that the pressure drop can be reasonably considered proportional to the square of the velocity normal to the screen at the vicinity of it through the resistance coefficient $k$, especially Taylor & Davies (Reference Taylor and Davies1944). More recently, Ito & Garry (Reference Ito and Garry1998) studied this problem in the 2-D case of a flow around and through a gauze for low resistance coefficient, while Eckert & Pflüger (Reference Eckert and Pflüger1942) studied the resistance coefficient for the case of a gauze spanning the entire section of a channel. In these cases, the pressure drop can be written as

(4.2)\begin{equation} \Delta p = \tfrac{1}{2} k \rho {v_{I}^2}_n(x_s,y_s,z_s). \end{equation}

The resistance coefficient $k$ depends on the geometry of the holes, the material of the screen and the Reynolds number based on the scale of the holes (sometimes $k$ is directly related to what is called the loss factor or friction factor). This has been mostly studied when the screen spans entirely a channel with normal incidence, and oblique incidence (Schubauer et al. Reference Schubauer, Spangenberg and Klebanoff1950; Reynolds Reference Reynolds1969). Note that in (3.18), $\theta (s)$ can be interpreted as a resistance coefficient, the intervention of the tangential velocity in (3.18) comes from the fact that in our case the fluid can pass around the porous structure.

To formulate the resistance coefficient dependency for porous screens, Pinker & Herbert (Reference Pinker and Herbert1967) has shown that the resistance coefficient $k$ can reasonably be considered as a product of a function of the solidity, $G(s)$, and a function $f$ of the local Reynolds number in a pore. In our definition of $Re_d$, we use the upstream velocity of the fluid far from the mesh; in order to take into account the changes in velocity close to the mesh, we define the local approach Reynolds number,

(4.3)\begin{equation} Re_n = Re_d \frac{v_{I_n}}{v_0}, \end{equation}

based on the scale of holes and the approach velocity that in our formulation corresponds to $v_{I_n}$. Among several fitted expressions for $G$ with respect to the solidity, Pinker & Herbert (Reference Pinker and Herbert1967) found that $G(s) = -\theta (s)$, which exhibits the best agreement with their data.

For the pressure drop, depending also on the inclination of the surface, we have to consider in fact the function $f(Re_n,\beta )$, as first proposed by Schubauer et al. (Reference Schubauer, Spangenberg and Klebanoff1950). Since the geometry of our porous screen is arbitrary, $\beta$ should be considered as a local characteristic of the inclination of the surface of the screen relative to the direction of the laminar free flow in the far upstream. Here $\beta = 0$ means that the surface is parallel to the flow $\boldsymbol {v}_0$. We thus propose the following relation for the pressure jump:

(4.4)\begin{align} {p_{II}}(x_s,y_s,z_s)-{p_{I}}(x_s,y_s,z_s) &= \tfrac{1}{2}\rho {v_I}_t^2(x_s, y_s,z_s) (1-E^2(\psi_I,\chi_I)) \nonumber\\ &\quad + \tfrac{1}{2}\rho{v_I}_n^2(x_s, y_s,z_s)\theta(s)f(Re_n,\beta). \end{align}

The expression of $f$ is thus expected to be found experimentally. As far as we know, there is no general physical formulation of the pressure drop through the holes at the microscopic level that can cover all types of screens, and there is no general demonstration of an analytical expression of $f$ for an arbitrary porous screen shape. Therefore, it is expected that some modifications are required for particular porous structures taking into account, for instance, the geometry of the holes. For screens composed of fibres, with an angle of attack $\beta$ with the upstream flow, and for $10^{-4} < Re_n < 10^4$, Brundrett (Reference Brundrett1993) gives the following empirical expression of $f$:

(4.5)\begin{align} f(Re_n,\beta) = \sin^2{(\beta)}\left(\frac{c_1}{Re_n\sin{(\beta)}}+ \frac{c_2}{\ln{(Re_n\sin{(\beta)}+1.25)}}+c_3\ln{(Re_n\sin{(\beta)})}\right). \end{align}

Here $c_1$, $c_2$ and $c_3$ are real constant. In this model (4.5), the first term corresponds to the laminar contribution, the second term to the turbulent friction and the last term is valid at large Reynolds numbers (Brundrett Reference Brundrett1993). Brundrett (Reference Brundrett1993) obtain a good fitting with his data for wire mesh screens by taking $c_1\approx 7.125$, $c_2\approx 0.88$, $c_3\approx 0.055$. Bailey et al. (Reference Bailey, Montero, Parra, Robertson, Baeza and Kamaruddin2003) found also a good fitting with their data by taking $c_1\approx 18$, $c_2\approx 0.75$ and $c_3\approx 0.055$. For the following sections of the paper, we take the geometric mean of the different values $c_1 \approx 11$, $c_2 \approx 0.8$ and $c_3 \approx 0.055$. For high $Re_n$, the function $f$ behaves as a logarithmic function of the Reynolds number $Re_n$, while for low $Re_n$, the variation of $f$ is much more pronounced (inverse function of the Reynolds number $Re_n$). We thus expect the viscous effects to be important for the typical flow speeds and mesh sizes considered in this work, and in applications such as fog harvesting and face masks.

We thus use this empirical expression, which is in principle only valid for fibrous screens, to account for local viscous effects in our experiments. The last term in the pressure jump (3.18) is thus multiplied by $f(\beta, Re_n)$. This is carried through all the calculation. In particular, our equation to determine the source strength (3.50) is modified such that

(4.6)\begin{align} & -\tfrac{1}{8} \omega^4\theta(s)f(Re_n,\beta) + \omega^2 \sin^2{(\beta)} (8\gamma_0 + \theta(s)f(Re_n,\beta) - 2) \nonumber\\ &\quad -4\omega\sin{(\beta)} - 2\sin^4{(\beta)}\theta(s)f(Re_n,\beta) = 0. \end{align}

This equation has been solved using Python 3.9.5 and the method fsolve from scipy.optimize. In the following, we keep and compare both formulations (3.50) and (4.6). For high solidity, and for a certain range of the two Reynolds numbers $Re$ and $Re_d$, the term $\theta (s)f(Re_n,\beta )$ should behave as the inertial term of the Darcy–Forchheimer equation for porous media. For porous screens composed of square fibre meshes, the inertial term of the Darcy–Forchheimer equation calculated with the method of Wang et al. (Reference Wang, Yang, Huang, Huang, Zhuan and Wu2021) has a value reasonably close to $\theta (s)$ when $s \approx 0.9$. For a very thin porous surface, as discussed by Teitel (Reference Teitel2010), the concept of permeability for porous media involved in the equation of Darcy, and Darcy–Forshheimer, may not always hold for the pressure loss through screens depending on the regime of the flow. According to Brundrett (Reference Brundrett1993) and Bailey et al. (Reference Bailey, Montero, Parra, Robertson, Baeza and Kamaruddin2003), the expression (4.5) seems to be valid over a larger flow regime.

4.2. Three-dimensional and viscous effects on the drag coefficient

As seen previously, our experiments suggest a strong effect of the local Reynolds number on the drag coefficient. We now compare in figure 10 our experimental results with the prediction of the 3-D model, including local viscous effects using the function $f(Re_n,\beta )$. For all mesh Reynolds numbers $Re_d$, the trend of the curve remains globally the same: there is an approximately linear increase of the drag coefficient at low solidity before the curve flattens and reaches a plateau, which is well represented by the 3-D model derived using the method proposed by Steiros & Hultmark (Reference Steiros and Hultmark2018), i.e. without $f(Re_n,\beta )$ (3.50). However, the slope of the initial linear part strongly depends on $Re_d$: it decreases with increasing mesh Reynolds number $Re_d$ for $Re_d<10^3$.

Figure 10. Drag coefficient as a function of the solidity for various square porous screens normal to the free flow with different $Re_d$, and comparison between the 3-D model. The full lines are obtained with the 3-D model developed in this paper, solution of (4.6). The black dash-dotted curve is the solution of (3.50), thus, without taking into account $Re_d$.

Using our model, (4.4) with the empirical formulation of Brundrett (Reference Brundrett1993), we obtain a good fitting with our data for solidity $s\leq 0.6$, for different Reynolds numbers. This shows the importance of both the 3-D effects and the viscous effects through the local Reynolds number $Re_d$. In particular, for a given solidity, the drag coefficient strongly decreases for increasing $Re_d$, as are plotted in figure 11.

Figure 11. Drag coefficient as a function of the Reynolds number $Re_d$ for three different narrow ranges of solidities: $s=0.1 \pm 0.02$, $s=0.28 \pm 0.04$, $s=0.585 \pm 0.025$. The bullets represent the experimental measurements, and the plain lines the results of our model ((4.6) taking into account the effect of mesh Reynolds number $Re_d$). The data points for the drag as well as the Reynolds number values are obtained for the velocity $5.0\leq v_0\leq 13\,{\rm m}\,{\rm s}^{-1}$. The dash-dotted lines represent the value of the drag coefficient for a solid plate (three dimensional) from our experiments and from Blevins (Reference Blevins1992).

We indeed observe a strong effect of the mesh Reynolds number, that also depends on the solidity. We observe a decrease of drag coefficient with increasing $Re_d$ that is well captured by the empirical formulation (4.6). This effect has also been observed for perforated plates by de Bray (Reference de Bray1957). At low $Re_d$, all curves tend towards the value of the drag coefficient of a flat solid plate, i.e. no fluid is passing through the screen and most of it is deviated around. We then note rapid variations for intermediate values of $Re_d$ ($5< Re_d<50$ depending on the solidity), to finally converge towards a constant value at high $Re_d$. This reduction of drag decreases with increasing solidity. As expected, for high solidity ($s=0.9$), there are no effects of $Re_d$, as the flow through the screen is weak. We obtain a similar drag coefficient using the values for the constants $c_1$, $c_2$ and $c_3$ from either Brundrett (Reference Brundrett1993) or Bailey et al. (Reference Bailey, Montero, Parra, Robertson, Baeza and Kamaruddin2003). Indeed, we see that for all solidity and moderate Reynolds numbers $Re_d={O}(200)$ that correspond to most of the screens we tested, the difference is small ($1\leq {C_D({\rm {Bailey}})}/{C_D({\rm {Brundrett}})}\leq 1.05$), and the difference is even lower with increasing Reynolds number $Re_d$. For low Reynolds numbers $Re_d={O}(5)$, the difference is however higher.

To further check the validity of our model, we made screens that have the same solidity but different hole size and number, keeping however $Re_d$ constant (see table 1 for $s=0.24$ (P14, P15 and P16) and $s=0.7$ (P8 and P25)), but still with a periodic distribution. As expected from our model, we do not observe any difference in the drag coefficient within the bounds of the measurement uncertainty, demonstrating that for these regular screens, the friction coefficient depends only on $s$ and $Re_d$.

4.3. Flow visualization

We can gain some insights into those behaviours by plotting the streamlines and velocity magnitude of the flow around and through the screen with our model (figure 12). For a screen normal to the flow, we observe that as the solidity increases, a larger part of the flow is deviated around the screen and that the flow is strongly slowed down in front of the screen, which is consistent with the measured increase of drag coefficient (figure 12ac).

Figure 12. Streamlines and velocity magnitude obtained with the 3-D model for a square screen normal to the flow at $v_0=2.0\,{\rm m}\,{\rm s}^{-1}$ for different solidities and the pressure jump defined in (3.18) (thus, without dependency on the mesh Reynolds number $Re_d$). Results are shown for (a) $s=0.3$ and $\beta = 90^{\circ }$; (b) $s=0.6$ and $\beta = 90^{\circ }$; (c) $s=0.9$ and $\beta = 90^{\circ }$.

This analysis is taken further by plotting in figure 13 both the PIV measurements (as described in § 2) in the symmetry plane of the screen, as well as the theoretical predictions using the pressure jump equation (3.18) for the sake of simplicity (thus, not accounting for the mesh Reynolds number $Re_d$ effect). As the solidity increases, we can observe a stronger attenuation of the velocity behind the screen, a larger deviation of the flow around the screen, as well as the apparition of a slower region upstream of the screen. We note that the theoretical prediction of the velocity magnitudes are in good agreement with the measured velocities even for the velocity attenuation downstream in the wake. We also observe that the attenuation of the velocity upstream is well captured by the model. However, for high solidity, the width of the wake appears larger in experiments than in the model, which might also be due to the presence of a thicker frame around the mesh in the experiments. Indeed, as observed in the streamlines in figure 4, there are vortices attached to the edges of the frame that may impact both the normal velocity and the shape of the wake.

Figure 13. Comparison between the 3-D theoretical model and our experiments for the prediction of the velocity field in the $(z,x)$ plane with $y=0$ for three different solidities (low, moderate and high). The theoretical velocity field has been obtained using the theoretical pressure jump law (3.18) and (3.19). Uncertainties on the velocities are estimated to be around $0.1\,{\rm m}\,{\rm s}^{-1}$. Results are shown for (a) P14, $s=0.24$, experimental; (b) P14, $s=0.24$, theoretical; (c) P6, $s=0.61$, experimental; (d) P6, $s=0.61$, theoretical; (e) P23, $s=0.82$, experimental; (f) P23, $s=0.82$, theoretical.

By mass conservation and since the normal velocity is assumed constant on the surface in this study, the proportion of the incoming fluid that goes through the screen is directly given by the dimensionless normal velocity ${v_n}/{v_0}$. Experimentally, we measure the velocity using a constant temperature anemometer from Dantec Dynamics (MiniCTA, with probe 55P11, tungsten wire with diameter $5\,\mathrm {\mu }$m and length 1.25 mm, precision of $0.01\,{\rm m}\,{\rm s}^{-1}$, minimum velocity of $0.20\,{\rm m}\,{\rm s}^{-1}$). In figure 14 we plot this ratio ${v_{I_n}}/{v_0}$ as a function of the solidity $s$ for different meshes, i.e. different $Re_d$ as presented in table 2. For a high mesh Reynolds number and low solidity, the model exhibits good agreements with the data. However, we observe a strong effect of $Re_d$ on the normal velocity: notably, for a small mesh Reynolds number $Re_d$, the normal velocity drops more rapidly than predicted by the model. The model also overpredicts the normal velocity for high solidity. For screens commonly used for fog harvesting ($Re_d = {O}(100)$), the normal velocity in the case where $Re_d$ is taken into account can be up to $21\,\%$ greater than the normal velocity in the case where $Re_d$ is not taken into account (according to our theoretical results, this is true for $0.25 < s < 0.75$), and is thus not negligible. These data are coherent with the PIV measurements.

Figure 14. Dimensionless normal velocity ${v_{I_n}}/{v_0}$ obtained with (4.6) for the solid lines, and with (3.50) for the dotted black line. The bullets represent the experimental measurements obtained using a hot wire anemometer.

Table 2. Reynolds number $Re_d$ calculated with the fibre diameter $d$ as characteristic size and with a velocity $v_0 = 2.84\,{\rm m}\,{\rm s}^{-1}$ and a kinematic viscosity $\nu = 15.6 \times 10^{-6}\,{\rm m}^2\,{\rm s}^{-1}$.

In figure 15 we plot the dimensionless tangential velocity distribution on the porous surface for different solidities $s$ and Reynolds number $Re_d$. We see that the tangential velocity is much lower than $v_0$ for a large part of the screen. We thus expect important differences in the streamlines that are deviated around the screen when the solidity or $Re_d$ varies.

Figure 15. Dimensionless tangential velocity magnitude ${|v_{I_t}|}/{v_0}$ on the screen surface obtained with the 3-D model for a square screen normal to a flow for different solidities and Reynolds numbers ((4.6) taking into account the effect of mesh Reynolds number $Re_d$). Results are shown for (ac) $Re_d = 10$; (df) $Re_d = 100$; (a) $s=0.1$ and $Re_d = 10$; (b) $s=0.5$ and $Re_d = 10$; (c) $s=0.8$ and $Re_d = 10$; (d) $s=0.1$ and $Re_d = 100$; (e) $s=0.5$ and $Re_d = 100$; (f) $s=0.8$ and $Re_d = 100$.

In figure 16 we give a representation of the proportion of the fluid passing through the porous surface (and thus, the proportion of the fluid going around it) for different solidities and Reynolds numbers ($Re_d = 10$ and $Re_d = 100$). The contours of these sections $S_a$ represent the separation between the fluid going through the surface and the fluid going around the surface at $z=-\infty$, i.e. the section of the streamtube or separation surface. If we place a passive tracer inside these sections at $z=-\infty$ then the passive tracer will flow across the surface. We observe a clear effect of the 3-D nature of the flow on these separation surfaces. As measured in our experiments, the separation distance along a diagonal is much higher than along the transverse direction, and this effect is amplified as $s$ increases and $Re_d$ decreases. The fraction ${S_a}/{S_0}$ is always equal to ${v_n}/{v_0}$. However, in applications such as filtration or fog harvesting, where the main drop capture mechanism is inertial impact (Moncuquet et al. Reference Moncuquet, Mitranescu, Marchand, Ramananarivo and Duprat2022), the tangential velocity and the shape of the streamlines are expected to impact the amount of droplets that leave the streamlines and are collected on the fibres.

Figure 16. Upstream section $S_a$ of the streamtube containing all the fluid passing through the surface, at $z=-\infty$, seen from the front, for $Re_d = 10$ and $Re_d = 100$.

4.4. Drag coefficient for low angle of attack

We then vary the orientation angle of the screen in the flow $\beta$. When the orientation angle is increased away from the normal incidence, the deviation of the streamlines is less important and the velocity is slightly higher (figure 17a,b). For the extreme case of a solid plate in figure 17(c), two streamlines are deviated along the plate, the flow slows down first and increases again with a peak value above the plate as the fluid particle leaves it and is re-entrained in the surrounding flow. However, if the asymptotic behaviour of the model may provide some indications about the global flow and the aerodynamic forces, it is expected to be outside of the assumptions of the model as discussed in the next section.

Figure 17. Streamlines and velocity magnitude obtained with the 3-D model for a square screen inclined to a flow at $v_0=2.0\,{\rm m}\,{\rm s}^{-1}$ for different solidities and the pressure jump defined in (3.18) (thus, without dependency on the mesh Reynolds number $Re_d$). Results are shown for (a) $s=0.8$ and $\beta = 60^{\circ }$; (b) $s=0.8$ and $\beta = 30^{\circ }$; (c) $s=1$ and $\beta = 30^{\circ }$.

The experimental results are compared with our model for the two angles $\beta =65^\circ$ and $\beta =43^\circ$ on figure 18. As for panels at a normal incidence, the drag coefficient increases with increasing solidity. The drag coefficient decreases with decreasing angle of attack, i.e. as we get further away from the normal incidence. Our model is in good agreement with the experimental data at $\beta =65^\circ$ (figure 18a). However, it underestimates the drag at orientations further from the normal incidence ($\beta =43^\circ$ figure 18b). In the case of an inclined porous screen composed of fibres, the effective solidity should increase as the angle of attack decreases. The use of this effective solidity would result in shifting our data closer to the yellow curve ($87 \leq Re_d \leq 225$) in figure 18(b). We finally observe that both the slope of the linear part and the final value at $s=1$ depend on the angle and that the drag coefficient decreases with decreasing angle.

Figure 18. Comparison between the 3-D model and our experiments for the prediction of the drag coefficient at two angles of inclination for various porous square screens: (a) $\beta = 65^\circ$ (b) $\beta = 43^\circ$ using (3.50). The black dash-dotted line is the theoretical result obtained with (4.6).

In addition, our model predicts a maximum of the drag coefficient at high solidity, as shown in figure 3 (or 10) and more pronounced in figure 18. Such a non-monotonic behaviour of the drag coefficient of permeable shells with solidity has been shown and demonstrated in the viscous regime for Stokes flow by Ledda et al. (Reference Ledda, Boujo, Camarri, Gallaire and Zampogna2021). In our case, at high $Re_L$ and a square screen normal to the flow, several of the porous screens have indeed (on average) a higher drag coefficient than a solid screen (e.g. for P3, P5, P8, P9, P22, P23, P24 and P25). However, due to the uncertainty and the interference drag with the frame used for the support of the screens (see Appendix C), the difference is not significant enough to draw a clear conclusion. We do not observe such a maximum in our measurements for inclined plates at an angle $43^{\circ }$. However, for an angle of $65^{\circ }$, the drag coefficient of the screen $P23$ is on average higher than the solid screen. It is also possible that in our experiments we were outside the flow regime for such a non-monotonic behaviour in the drag coefficient.

5. Model extensions

5.1. Vena contracta

In our model, as done also by Steiros & Hultmark (Reference Steiros and Hultmark2018), we neglected the possible contraction of the flow within the holes. Considering the vena contracta within the pores, as shown in figure 19, the characteristic velocity immediately upstream given by (3.14), i.e. the mean velocity after contraction of the flow within the holes ${v_h}_1$, can be expressed as

(5.1)\begin{equation} {v_h}_1 = \frac{{v_I}_n(x_s, y_s, z_s)}{C(1-s)}, \end{equation}

where $C$ is the contraction coefficient. This implies that the coefficient $\theta (s)$ given in (3.19) is

(5.2)\begin{equation} \theta(s) = 1-\frac{1}{C^2(1-s)^2}. \end{equation}

As far as we know, there are no measurements of vena contracta for screens constituted of fibres and especially in the 3-D case of free flow. Simmons & Cowdrey (Reference Simmons and Cowdrey1945) performed measurements of the velocity profile behind a porous screen made of a square mesh of woven material spanning a section of a channel and suggest an estimate of the vena contracta assuming a uniform velocity ${v_h}_1 = v_0$ in the holes. In these experiments the screens are made of circular rods with diameters ranging from 0.112 to 0.373 mm arranged in a square mesh. The velocity used was from $2.44$ to $10.36\,{\rm m}\,{\rm s}^{-1}$. The local Reynolds number $Re_d$ thus varies from $18$ to $248$, which is the same order of most of the porous screens used in our experiments. For a solidity $s\approx 0.5$, the coefficient $C$ should be between $0.9$ and $1.0$. Note that with the formulation (5.2), if $C \neq 1$, the limit $s=0$ leads to a non-zero pressure difference, suggesting that $C$ may depend on the solidity $s$ at least for low solidities. We note that if the contraction of the flow is not neglected in $\theta (s)$, i.e. $C\neq 1$, there is an increase of the drag coefficient compared with the curve plotted in figure 9(a). We would thus still overestimate the drag compared with our data in figure 3.

Figure 19. Diagram of the flow at the scale of the pores in case of parallel fibres (left), and corresponding imaginary surface as used in the model (right).

5.2. Asymptotic behaviour at large solidities

Although the model is built for porous surfaces, it is interesting to explore its asymptotic behaviour when the solidity tends to $1$. In this limit, one can question whether a porous description of the surface is still valid.

When the screen tends to a solidity equal to one, we can obtain an expression of the drag coefficient. Taking (4.6), dividing by $\theta (s)$, which tends to $+\infty$ when $s \to 1$, we have

(5.3)\begin{equation} (\omega^2-4\sin^2{(\beta)})^2=0. \end{equation}

Therefore (excluding the case when $\omega$ is negative, which would not correspond to the type of flow we study in this paper), we obtain

(5.4)\begin{equation} \omega = 2\sin{(\beta)}, \end{equation}

which gives ${v_I}_n=0$ as expected. Introducing this value into the expression of the drag coefficient gives

(5.5)\begin{equation} C_D = 2\sin^3{(\beta)}(1-4\gamma_0). \end{equation}

Note that the solid drag coefficient does not depend on the assumption of the pressure loss across the screen through the expression of $\theta (s)f(Re_n,\beta )$ (in (4.4)) as expected.

We measured the drag coefficient as a function of the angle of attack for a square plate at solidity $1$. We plotted the result in figure 20, and as we can see, there is a gap between the prediction (5.5) and the experimental value of the drag coefficient especially at a low angle of attack.

Figure 20. Comparison between the drag coefficient prediction of the 3-D model and the experimental measurement for a square plate in a free flow for different angles of inclination. Data are added from Blevins (Reference Blevins1992), Okamoto & Azuma (Reference Okamoto and Azuma2011) and Torres & Mueller (Reference Torres and Mueller2004).

Several explanations for this gap can be listed. First, for the solid plate, assuming a no-slip boundary condition, the tangential component of the velocity on the surface is equal to zero and increases gradually in the boundary layer. Here, this component of the velocity is not equal to zero but takes a value that corresponds to the conservation of the pressure head along the streamlines in region I. Moreover, it is clearly possible that the pressure difference $p_I\unicode{x2013}p_{II}$ is not always well determined through the assumptions of the model, especially for low angle of attacks (or particular shapes) where detachment of the boundary layer and the formation of separation bubbles can occur that strongly affects the pressure distribution around the surface (Crompton & Barrett Reference Crompton and Barrett2000).

5.3. Vortex shedding

Our model is only valid for a steady wake, i.e. without vortex shedding. This could be taken into account, but is beyond the scope of the current study. However, we believe vortex shedding has little influence on our experimental results. For instance, it has been shown that vortex shedding was suppressed for the flow past a porous cylinder (Ledda et al. Reference Ledda, Siconolfi, Viola, Gallaire and Camarri2018). In fact, in our flow field measurements, we do not see any vortex shedding close to the plate. To evaluate the influence of the vortex shedding on our measurements, we used two splitter plates of different length located in the wake as done by Steiros & Hultmark (Reference Steiros and Hultmark2018). We note that this method was proposed for a 2-D flow, and might not apply directly to our 3-D configuration. In fact, we are not aware of experiments to suppress vortex shedding in three dimensions but we believe that splitter plates can still give pertinent information on vortex shedding. Following Apelt & West (Reference Apelt and West1974), for Reynolds numbers in the range $10^4< Re<5\times 10^4$, which corresponds to our case, the use of a splitter plate three times longer than the plate already suppresses the vortex shedding. In our experiments the splitter plate has the same height as the solid plate on which drag is measured but with length 11 and 30 cm and thickness of 3 mm. Figure 21 shows the measured drag force for a solid plate (solidity equal to $1$, square of $11.0\times 11.0\,{\rm cm}^2$, 3 mm thickness) with and without a splitter plate. We see a reduction of the drag coefficient of $0.07$ with a splitter plate three times longer than the plate. This difference is expected to be less and less significant with decreasing solidity.

Figure 21. Comparison of the drag force $F_D$ with and without a splitter plate for a solid plate (without frame support) of dimensions $11.0\times 11.0\,{\rm cm}^2$ and thickness 3 mm. The fitting curves are obtained using a quadratic law.

While these results are not sufficient to rule out completely the effect of vortex shedding, we note that vortex shedding increases the drag significantly and reduces the base pressure, as seen in drag measurements in two dimensions by Steiros & Hultmark (Reference Steiros and Hultmark2018); these effects cannot explain the difference between the theory and the experimental data since suppressing vortex shedding could lead to a lower value of the drag coefficient compared with our measurements, thus, the difference between the theoretical model and the data would be more important.

6. Conclusion

We have studied both experimentally and theoretically the flow around and through porous screens. In particular, we have derived a model with a universal core based on clear physical assumptions that allows us to probe 3-D and shape effects. The main issue lies in the modelling of the pressure jump law across the screen. At low Reynolds numbers, Zampogna & Gallaire (Reference Zampogna and Gallaire2020) derived an analytical expression using homogenization tools. However, it is still an issue at moderate or high Reynolds numbers where no explicit relation has been derived directly from Navier–Stokes equations (Wang et al. Reference Wang, Yang, Huang, Huang, Zhuan and Wu2021). We propose here to use an empirical relation obtained for fibrous screens in the literature. The semi-empirical model accurately predicts our experiments, with the only fitting parameters in the pressure jump; we note that these parameters are taken identical to those obtained in previous studies for other screens, suggesting that they are universal for any wire screens (or textiles) that are widely used in applications. For other types of screens, including perforated plates, once the pressure jump law across is known, it can be directly implemented in the model.

In addition, the empirical function includes both laminar and turbulent viscous friction, and is thus applicable to a wide range of local Reynolds numbers. The interest of the semi-empirical model lies in the fact that it uses only three types of information that are readily available: the solidity, the macroscopic geometry of the screen and the Reynolds number based on the fibre diameter. Both viscous effects (i.e. the local Reynolds number $Re_d$) and 3-D effects have a significant impact on the flow and aerodynamic forces. We performed experiments on more than $30$ porous screens composed of fibres to measure the drag force and the flow deviation around the screen. We show both theoretically and experimentally that for screens of identical solidity, the lower the local Reynolds number based on the fibre diameter is, the higher the drag coefficient is, and the stronger the deviations are. The drag coefficient for square porous screens, either normal or with a high angle of attack, show a good agreement with the model except at very high solidity for which however the prediction is improved compared with previous models. Furthermore, the model can predict the proportion of the flow deviated around the screen, as well as the detailed separation surface and the tangential velocity, which are important quantities in many applications such as filtration or fog harvesting.

Our model might also be useful for a non-homogeneous porous screen. Indeed, for high solidity perforated square plates, de Bray (Reference de Bray1957) found that the drag coefficient depends on the distribution of the perforations, i.e. is slightly higher with outer holes than inner holes. It is worth emphasizing that it should be possible in our model to implement such a surface with non-homogeneous solidity and, therefore, try to reproduce such a difference in the drag coefficient. However, this consideration goes beyond the scope of this paper and will be explored in a further study.

We note that this model may find an application to the wind tunnel blockage correction (Steiros et al. Reference Steiros, Bempedelis and Cicolin2022) and turbines modelling (Ayati et al. Reference Ayati, Steiros, Miller, Duvvuri and Hultmark2019). According to Steiros et al. (Reference Steiros, Bempedelis and Cicolin2022), the use of a porous plate potential model as the one we use in the present paper improves the blockage correction accuracy for moderate and high solidity compared with other models. However, these models are in two dimensions (Ayati et al. Reference Ayati, Steiros, Miller, Duvvuri and Hultmark2019; Steiros et al. Reference Steiros, Bempedelis and Cicolin2022), while the turbine can have a circular 3-D structure. Our theoretical and experimental results show that the use of a 3-D model may improve the accuracy at moderate and high solidity for the drag compared with the 2-D models (see figures 9(a) and 3). Moreover, our experimental results show that the normal velocity $v_{I_n}$ on the screen may be significantly overestimated at high solidity by both 2-D and 3-D models. We stress again that our model assumes a steady wake, and thus, is not applicable in the presence of vortex shedding that occurs at high solidity for a certain range of Reynolds numbers.

Further study should focus on more complex shapes with curvature; a straightforward formulation would be to use several small rectangular plates like the well-known panel method (adapting therefore the velocity potential $\phi$). Furthermore, they should focus on the pressure and velocity distribution around the porous screen, as there are very few, if any, data for free 3-D flow. Indeed, we only found data for the drag in 3-D free flow in Prandtl & Flachsbart (Reference Prandtl and Flachsbart1932) and more recently Letchford (Reference Letchford2001).

Acknowledgements

We are grateful to the reviewers for their careful reading of the paper. Their detailed comments have led to an improvement of the present study.

Declaration of interests

The authors report no conflict of interest.

Appendix A

This appendix contains the details of the calculation of the velocities for a rectangular screen inclined with an angle $\beta$ (with respect to the $z$ axis) in a laminar flow.

Using the previous notations in § 2, the velocity potential in region I is

(A1)\begin{align} \phi_I(x,y,z) &= v_0z + c + \phi(x,y,z) \nonumber\\ &= v_0z + c - \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{\varOmega(u,v)\,\mathrm{d}u\,\mathrm{d}v}{\sqrt{(x-v)^2+(y-u\sin{(\beta)})^2+(z-u\cos{(\beta)})^2}}. \end{align}

A.1. Normal component of the velocity

From the velocity potential we deduce the velocity in region I:

(A2)\begin{equation} \boldsymbol{v_I}(x,y,z) = \left(v_0 + \frac{\partial \phi}{\partial z}(x,y,z)\right)\boldsymbol{e_z} + \frac{\partial \phi}{\partial y}(x,y,z)\boldsymbol{e_y} + \frac{\partial \phi}{\partial x}(x,y,z) \boldsymbol{e_x}. \end{equation}

The vector normal to the surface with an angle $\beta$ is

(A3)\begin{equation} \boldsymbol{e_n} ={-}\cos{(\beta)}\boldsymbol{e_y}+\sin{(\beta)}\boldsymbol{e_z}. \end{equation}

Therefore, the component of the velocity normal to the surface is

(A4)\begin{equation} \boldsymbol{{v_I}_n}(x,y,z) = \left(\sin{(\beta)}\left(v_0 + \frac{\partial \phi}{\partial z}(x,y,z)\right) -\cos{(\beta)} \frac{\partial \phi}{\partial y}(x,y,z)\right)\boldsymbol{e_n}, \end{equation}

with

(A5)$$\begin{gather} \frac{\partial \phi}{\partial z}(x,y,z) = \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{(z-u\cos{(\beta)})\varOmega(u,v)}{((x-v)^2+(y-u\sin{(\beta)})^2+ (z-u\cos{(\beta)})^2)^{{3}/{2}}}\mathrm{d}u\,\mathrm{d}v, \end{gather}$$
(A6)$$\begin{gather}\frac{\partial \phi}{\partial y}(x,y,z) = \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{(y-u\sin{(\beta)})\varOmega(u,v)}{((x-v)^2+(y-u\sin{(\beta)})^2+ (z-u\cos{(\beta)})^2)^{{3}/{2}}}\mathrm{d}u\,\mathrm{d}v. \end{gather}$$

The normal velocity magnitude is therefore

(A7)\begin{equation} {v_I}_n(x,y,z) = \sin{(\beta)}v_0 + \frac{1}{4{\rm \pi}}I(x,y,z), \end{equation}

with

(A8)\begin{equation} I(x,y,z) = \iint_{\mathcal{S}_p} \frac{(\sin{(\beta)}z - \cos{(\beta)}y) \varOmega(u,v)}{((x-v)^2+(y-u\sin{(\beta)})^2+(z-u\cos{(\beta)})^2)^{{3}/{2}}}\mathrm{d}u\,\mathrm{d}v. \end{equation}

When this component is evaluated on the surface, the integral becomes singular at the position of this evaluation. We have to calculate the value of this singularity. To simplify the calculation, we introduce two parameters:

(A9)\begin{equation} \left.\begin{gathered} h = y\sin{(\beta)} + z\cos{(\beta)}, \\ t = y\cos{(\beta)} - z\sin{(\beta)}. \end{gathered}\right\} \end{equation}

Thus, we have inversely

(A10)\begin{equation} \left.\begin{gathered} y = h\sin{(\beta)} + t\cos{(\beta)}, \\ z = h\cos{(\beta)} - t\sin{(\beta)}. \end{gathered}\right\} \end{equation}

Denoting $\tilde {I}$ the integral $I$ with the new parameters and without the source strength $\varOmega$ assumed to be continuous, we obtain

(A11)\begin{align} \tilde{I}(x,h,t) &={-}\int_{v_a}^{v_b}\int_{u_a}^{u_b} \frac{t}{((x-v)^2+t^2+(h-u)^2)^{{3}/{2}}}\mathrm{d}u\,\mathrm{d}v \nonumber\\ &= F_{u_a,v_a}(x,h,t) - F_{u_b,v_a}(x,h,t) - F_{u_a,v_b}(x,h,t) + F_{u_b,v_b}(x,h,t), \end{align}

with

(A12)\begin{equation} F_{u_a,v_a}(x,h,t) = \arctan{\left(\frac{(h-u_a)(v_a-x)}{t\sqrt{(x-v_a)^2+t^2+(h-u_a)^2}}\right)}. \end{equation}

We can now introduce

(A13)\begin{equation} \left.\begin{gathered} y = w\sin{(\beta)}, \\ z = (w+\epsilon)\cos{(\beta)}. \end{gathered}\right\} \end{equation}

The parameters become

(A14)\begin{equation} \left.\begin{gathered} h = w, \\ t ={-}\epsilon\cos{(\beta)}\sin{(\beta)}. \end{gathered}\right\} \end{equation}

At the surface $\epsilon \to 0^{\pm }$ ($\pm$ depending on the direction from which we approach the surface, upstream or downstream), denoting $\tilde {\tilde {I}}$ the integral $I$ with the new parameters $w$, we have

(A15)$$\begin{gather} {v_I}_n(x,w) = v_n^-, \end{gather}$$
(A16)$$\begin{gather}{v_{II}}_n(x,w) = E v_n^+, \end{gather}$$

with

(A17)\begin{equation} v_n^{{\pm}} = \sin{(\beta)}v_0 + \frac{1}{4{\rm \pi}} \lim_{\epsilon \to 0^{{\pm}}} \tilde{\tilde{I}}(x,w). \end{equation}

Whatever the constant integral bounds, we have the following limit:

(A18)\begin{equation} \lim_{\epsilon \to 0^{{\pm}}} \tilde{\tilde{I}}(x,w) ={\pm} 2{\rm \pi} \varOmega(x,w). \end{equation}

Therefore, the normal component of the velocity at the position $(x,w)$ on the porous surface is

(A19)\begin{equation} {v_I}_n(x,w) = \sin{(\beta)}v_0 - \tfrac{1}{2}\varOmega(x,w) \end{equation}

and

(A20)\begin{equation} {v_{II}}_n(x,w) = E\left(\sin{(\beta)}v_0 + \tfrac{1}{2}\varOmega(x,w)\right). \end{equation}

A.2. Tangential component of the velocity

We define the tangential vector on the porous surface as

(A21)\begin{equation} \boldsymbol{t} = \boldsymbol{{t_1}} + \boldsymbol{{t_2}}, \end{equation}

with

(A22)$$\begin{gather} \boldsymbol{{t_1}} = \boldsymbol{e_x}, \end{gather}$$
(A23)$$\begin{gather}\boldsymbol{{t_2}} = \sin{(\beta)}\boldsymbol{e_y}+ \cos{(\beta)}\boldsymbol{e_z}. \end{gather}$$

Thus, the magnitude of the tangential component of the velocity in region I is

(A24)\begin{align} {v_I}_t(x,y,z) &= \|(\boldsymbol{v_I}\boldsymbol{\cdot}\boldsymbol{{t_1}})\boldsymbol{{t_1}} + (\boldsymbol{v_I}\boldsymbol{\cdot}\boldsymbol{{t_2}})\boldsymbol{{t_2}}\| \nonumber\\ &= \left(\left(\frac{\partial \phi}{\partial x}\right)^2 + \left(\sin{(\beta)}\frac{\partial \phi}{\partial y} + \cos{(\beta)}\left(v_0 + \frac{\partial \phi}{\partial z}\right)\right)^2\right)^{{1}/{2}}. \end{align}

If we consider a constant source strength $\varOmega$ then we can write the tangential velocity at the point $(x,w)$ on the surface with $w$ defined in (A13) as

(A25)\begin{align} {v_I}_t(x,w) = (\varOmega^2 \mathcal{I}_x^2(x,w) + (\varOmega(\sin{(\beta)}\mathcal{I}_y(x,w)+\cos{(\beta)}\mathcal{I}_z(x,w)) + \cos{(\beta)}v_0)^2)^{{1}/{2}}. \end{align}

For simplicity, we use the root mean square of the magnitude of the tangential component of the velocity on the surface. We have

(A26)\begin{align} {v_I}_t &= \frac{1}{\sqrt{S_p}}\left(\iint_{\mathcal{S}_p} {v_I}_t^2(x,w)\,\mathrm{d}\kern 0.06em x\,\mathrm{d}w\right)^{{1}/{2}} \nonumber\\ &= (\varOmega^2 \gamma(\beta) + \varOmega \alpha(\beta)v_0 + v_0^2\cos^2{(\beta)})^{{1}/{2}}, \end{align}

with

(A27)\begin{equation} \gamma(\beta) = \frac{1}{S_p}\iint_{\mathcal{S}_p} \mathcal{I}_x^2(x,w) + (\sin{(\beta)} \mathcal{I}_y(x,w) + \cos{(\beta)} \mathcal{I}_z(x,w))^2\,\mathrm{d}\kern 0.06em x\,\mathrm{d}w, \end{equation}

and

(A28)\begin{equation} \alpha(\beta) = \frac{1}{S_p}\iint_{\mathcal{S}_p} 2\cos{(\beta)} (\sin{(\beta)} \mathcal{I}_y(x,w) + \cos{(\beta)} \mathcal{I}_z(x,w))\,\mathrm{d}\kern 0.06em x\,\mathrm{d}w, \end{equation}

with the following integrals:

(A29)$$\begin{gather} \mathcal{I}_x(x,w) = \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{x-v}{((x-v)^2+(w-u)^2)^{{3}/{2}}}\mathrm{d}u \,\mathrm{d}v, \end{gather}$$
(A30)$$\begin{gather}\mathcal{I}_y(x,w) = \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{(w-u)\sin{(\beta)}}{((x-v)^2+(w-u)^2)^{{3}/{2}}}\mathrm{d}u\,\mathrm{d}v, \end{gather}$$
(A31)$$\begin{gather}\mathcal{I}_z(x,w) = \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{(w-u)\cos{(\beta)}}{((x-v)^2+(w-u)^2)^{{3}/{2}}} \mathrm{d}u\,\mathrm{d}v. \end{gather}$$

There is no particular difficulty to compute the integrals, if we define a function $F_{a,b}: (x,w)\mapsto F_{a,b}(x,w)$ as well as $G_{a,b}: (x,w)\mapsto G_{a,b}(x,w)$ where $\vert a \vert \geq \vert x\vert$, $\vert b \vert \geq \vert w\vert$,

(A32)\begin{equation} F_{a,b}(x,w) = \frac{1}{4{\rm \pi}} \ln{(\sqrt{(x-a)^2+(w-b)^2} + w - b)} \end{equation}

and

(A33)\begin{equation} G_{a,b}(x,w) = \frac{1}{4{\rm \pi}} \ln{(\sqrt{(x-a)^2+(w-b)^2} + x - a)}, \end{equation}

then if we integrate over the rectangular domain $[v_a,v_b]\times [u_a,u_b]$, we have

(A34)\begin{align} \mathcal{I}_x(x,w) &= F_{v_a,u_a}(x,w) - F_{v_a,u_b}(x,w) - F_{v_b,u_a}(x,w) + F_{v_b,u_b}(x,w), \end{align}
(A35)\begin{align} \mathcal{I}_y(x,w) &= \sin{(\beta)}(G_{v_a,u_a}(x,w) - G_{v_a,u_b}(x,w) \nonumber\\ &\quad - G_{v_b,u_a}(x,w) + G_{v_b,u_b}(x,w)), \end{align}
(A36)\begin{align} \mathcal{I}_z(x,w) &= \cos{(\beta)}(G_{v_a,u_a}(x,w) - G_{v_a,u_b}(x,w) \nonumber\\ &\quad - G_{v_b,u_a}(x,w) + G_{v_b,u_b}(x,w)), \end{align}

and those expressions are integrable again over the rectangular domain $[v_a,v_b]\times [u_a,u_b]$.

For a screen normal to the free flow $\boldsymbol {v_0}$, $\beta = {{\rm \pi} }/{2}$; thus, the mean tangential velocity is reduced to

(A37)\begin{equation} {v_I}_t = \varOmega\sqrt{\gamma\left(\frac{\rm \pi}{2}\right)}. \end{equation}

A.3. Value of the shape parameters for a rectangular screen with different aspect ratio

Starting from (A27) and (A28) we can rewrite the expression of $\gamma (\beta )$ and $\alpha (\beta )$ in order to separate the shape terms and the inclination terms, we see for this particular case of a rectangular plate that these expressions can be simplified as

(A38)\begin{equation} \gamma(\beta) = \gamma_0 \end{equation}

and

(A39)\begin{equation} \alpha(\beta) = 2\alpha_0\cos{(\beta)}, \end{equation}

with

(A40)\begin{align} \mathcal{J}(x,w) &= \frac{1}{4{\rm \pi}}\iint_{\mathcal{S}_p} \frac{w-u}{((x-v)^2+(w-u)^2)^{{3}/{2}}}\mathrm{d}u \,\mathrm{d}v \nonumber\\ &= G_{v_a,u_a}(x,w) - G_{v_a,u_b}(x,w)- G_{v_b,u_a}(x,w) + G_{v_b,u_b}(x,w), \end{align}
(A41)\begin{align} &\qquad\quad \left.\begin{gathered} \gamma_0 = \frac{1}{S_p}\iint_{\mathcal{S}_p} \mathcal{I}^2_x(x,w) + \mathcal{J}^2(x,w)\,\mathrm{d}\kern 0.06em x\,\mathrm{d}w,\\ \alpha_0 = \frac{1}{S_p}\iint_{\mathcal{S}_p} \mathcal{J}(x,w)\,\mathrm{d}\kern 0.06em x\,\mathrm{d}w. \end{gathered}\right\} \end{align}

For symmetry reasons, $\alpha _0 = 0$. Thus, the mean tangential component of the velocity at the surface for a rectangular plate inclined in the flow is actually

(A42)\begin{equation} {v_I}_t = (\varOmega^2 \gamma_0 + v_0^2\cos^2{(\beta)})^{{1}/{2}}. \end{equation}

We computed the values of the parameter $\gamma _0$ in table 3 for different aspect ratios ($v_b = - v_a = {l_1}/{2}$ and $u_b = -u_a = {l_2}/{2}$). We used Python 3.9.5 with the method nquad from scipy.integrate.

Table 3. Shape parameter for rectangular plate with different aspect ratios. Here $l_1$ and $l_2$ are the lengths of the sides of the rectangular plate. The aspect ratio is ${l_2}/{l_1}$.

Appendix B

This appendix contains the details of the calculation of the source strength $\varOmega$ for the case of a rectangular screen inclined to the laminar free flow with an angle $\beta$ with respect to the $z$ axis. We start from the following equations.

The pressure difference is expressed as

(B1)$$\begin{gather} {p_{II}}-{p_{I}} = \tfrac{1}{2}\rho ({v_I}_t^2(1-E^2)+{v_I}_n^2\theta(s)f(Re_n,\beta)), \end{gather}$$
(B2)$$\begin{gather}{p_{III}}-{p_{0}}= \tfrac{1}{2}\rho((1-E^2)v^2_0 + {v_I}_n^2\theta(s)f(Re_n,\beta)). \end{gather}$$

The drag forces are expressed as

(B3)$$\begin{gather} F_D = \rho v_0(1-E){v_n} {S_p} + \frac{1}{v_0}(p_0\unicode{x2013}p_{III})\frac{{v_n}}{E} {S_p}, \end{gather}$$
(B4)$$\begin{gather}F_D = ({p_I}\unicode{x2013}{p_{II}})\sin{(\beta)}{S_p} + \rho {v_n} v_0\cos^2(\beta)(1-E){S_p}. \end{gather}$$

In these expressions the normal and tangential components of the velocity at the surface of the screen are

(B5)$$\begin{gather} {v_I}_n = \sin{(\beta)}v_0 - \tfrac{1}{2}\varOmega, \end{gather}$$
(B6)$$\begin{gather}{v_I}_t = (\varOmega^2 \gamma_0 + v_0^2\cos^2{(\beta)})^{{1}/{2}}. \end{gather}$$

Then, by denoting $\omega = {\varOmega }/{v_0}$ we have, for the velocities,

(B7)$$\begin{gather} \frac{{v_I}_n}{v_0} = \sin{(\beta)} - \frac{1}{2}\omega, \end{gather}$$
(B8)$$\begin{gather}\frac{{v_I}_t}{v_0} = (\omega^2 \gamma_0 + \cos^2{(\beta)})^{{1}/{2}}. \end{gather}$$

The attenuation coefficient $E$ is

(B9)\begin{equation} E = \frac{\sin{(\beta)} - \dfrac{1}{2}\omega}{\sin{(\beta)} + \dfrac{1}{2}\omega}. \end{equation}

For the pressure differences, we obtain

(B10)$$\begin{gather} \frac{{p_{II}}-{p_{I}}}{\dfrac{1}{2}\rho v_0^2} = (\omega^2 \gamma_0 + \cos^2{(\beta)})\frac{2\omega\sin{(\beta)}}{\left(\sin{(\beta)} + \dfrac{1}{2}\omega\right)^2}+\left(\sin{(\beta)}-\frac{1}{2}\omega\right)^2\theta(s)f(Re_n,\beta), \end{gather}$$
(B11)$$\begin{gather}\frac{{p_{III}}-{p_{0}}}{\dfrac{1}{2}\rho v_0^2} = \frac{2\omega \sin{(\beta)}}{\left(\sin{(\beta)} + \dfrac{1}{2}\omega \right)^2}+ \left(\sin{(\beta)} - \frac{1}{2}\omega \right)^2\theta(s)f(Re,\beta). \end{gather}$$

The first expression of the drag force is

(B12)\begin{equation} \frac{F_D}{\dfrac{1}{2}\rho v_0^2 S_p} = \frac{2\omega\left(\sin{(\beta)}-\dfrac{1}{2}\omega\right)}{\sin{(\beta)} + \dfrac{1}{2}\omega} + \left(\sin{(\beta)} + \frac{1}{2}\omega\right)\frac{{p_{0}}\unicode{x2013}{p_{III}}}{\dfrac{1}{2}\rho v_0^2}. \end{equation}

The second expression of the drag force is

(B13)\begin{equation} \frac{F_D}{\dfrac{1}{2}\rho v_0^2 S_p} = \frac{{p_{I}}\unicode{x2013}{p_{II}}}{\dfrac{1}{2} \rho v_0^2}\sin{(\beta)} + \cos^2(\beta)\frac{2\omega\left(\sin{(\beta)} - \dfrac{1}{2}\omega\right)}{\sin{(\beta)}+ \dfrac{1}{2}\omega}. \end{equation}

The two expressions of the drag coefficient are

(B14)\begin{align} C_D &= \frac{-\omega^2}{\sin{(\beta)} + \dfrac{1}{2}\omega} - \left(\sin{(\beta)} + \frac{1}{2}\omega \right)\left(\sin{(\beta)} - \frac{1}{2}\omega \right)^2\theta(s)f(Re_n,\beta), \end{align}
(B15)\begin{align} C_D &= \cos^2(\beta)\frac{2x\left(\sin{(\beta)} - \dfrac{1}{2}\omega\right)}{\sin{(\beta)}+ \dfrac{1}{2}\omega } - \sin{(\beta)}\left(\sin{(\beta)} - \frac{1}{2}\omega\right)^2\theta(s)f(Re_n,\beta) \nonumber\\ &\quad - (\omega^2 \gamma_0 + \cos^2{(\beta)} ) \frac{2\omega\sin^2{(\beta)}}{\left(\sin{(\beta)} + \dfrac{1}{2}\omega\right)^2}. \end{align}

By combining the expression (B14) and (B15), we obtain the equation

(B16)\begin{align} & -\tfrac{1}{8} \omega^4\theta(s)f(Re_n,\beta) + \omega^2 \sin^2{(\beta)} (8\gamma_0 + \theta(s)f(Re_n,\beta) - 2) \nonumber\\ &\quad -4\omega\sin{(\beta)} - 2\sin^4{(\beta)}\theta(s)f(Re_n,\beta) = 0. \end{align}

Note that, for a rectangular screen normal to the free flow, this equation is reduced to

(B17)\begin{align} -\tfrac{1}{8} \omega^4\theta(s)f(Re_n,\beta) + \omega^2 (8\gamma_0 + \theta(s)f(Re_n,\beta) - 2) - 4\omega - 2\theta(s)f(Re_n,\beta) = 0. \end{align}

If the length $L$ (or the width $D$) of the rectangular screen is infinitely long then we can show that

(B18)\begin{equation} \lim_{L\to\infty} \gamma_0 = \frac{1}{12}, \end{equation}

and we tend to the 2-D case studied by Steiros & Hultmark (Reference Steiros and Hultmark2018).

Appendix C

This appendix contains the details of the experimental data and their processing.

C.1. Screen samples

The characteristics of the different porous screens used for the measurement are presented in table 4. They are mainly square wire mesh screens, as represented in figure 22. In addition, five other types of screen were used to test the robustness of the model.

Table 4. Porous screen characteristics. The Reynolds number $Re_d$ is calculated with a velocity $v_0 = 10\,\textrm {m}\,\textrm {s}^{-1}$ and a kinematic viscosity $\nu = 15.6 \times 10^{-6}\,\textrm {m}^2\,\textrm {s}^{-1}$, except for the screens P27, P28 and P29 where the velocity is $v_0 = 5\,\textrm {m}\,\textrm {s}^{-1}$.

Table 5. Summary of the values of the coupling drag force for a flat plate (solidity $s=1$) for different velocities. All the values have an uncertainty of approximately 0.020 N. For the measurement of $F_p$ and $F_r$, a solid flat plate with 4.0 mm thickness has been used, corresponding to the thickness of the frame used for measuring $F_c$ and $F_t$. The drag coefficients are obtained using the dimensions shown in figure 23 with fluid density estimated from the measurement of temperature, pressure and humidity during the different experiments.

Figure 22. Diagram of the fibre mesh.

C.2. Correction of the coupling drag

Due to the mast and frame supporting the porous structure, to obtain the drag coefficient of the porous structure from the raw data composed of the measured forces denoted $F_{t+m}$, the contributions of each part must be decoupled.

In our analysis of the data, we consider the most simple assumption that the coupling is negligible. Therefore, the drag force of the porous screen is the drag force of the total system (mast, frame and porous screen) minus the drag force of the mast and the frame (without the porous screen, measured before the series of measurements). In order to justify this assumption, we estimated the coupling drag force (the additional term due to the interference between the frame and the screen) for what we assumed to be the worst case, that is, for the solidity equal to $1$.

As far as we know, this coupling is nonlinear and there is no general method. We adopt the approach we detail here, based on different measurements with and without the frame illustrated in figure 23. Due to the elongated shape of the mast and the way it is connected to the screen, it is reasonable to assume that the drag force of the mast $F_m$ and the rest of the system $F_t$ add up (giving what we denote $F_{t+m}$). For the frame, the coupling with the porous structure is expected to be more important. We measured the drag force of the mast alone $F_m$ and subtracted the value of $F_{t+m}$. We measured the drag force of the mast alone $F_m$ and substracted it from the total force $F_{m+t}$. The drag force of the porous screen is denoted $F_p$ and the coupling term is denoted $\Delta F$. The total force measured of the frame and the porous screen $F_t$ can be written as

(C1)\begin{equation} F_t = F_c + F_p + \Delta F. \end{equation}

Figure 23. Diagram of the screens (in dark gold) and the frame support (in green).

Figure 24. Drag force for the plates and the frame represented in figure 23. Values and drag coefficients are in table 5. Data are fitted using a quadratic law with respect to the velocity $v_0$. The surface is orthogonal to the mean far-field flow direction $\boldsymbol {v}_0$.

To determine the coupling (or interference) term $\Delta F$, we know exactly the value of the drag force of the porous structure for two points: the first at solidity $s=0$ where $F_p= 0$, and the second at solidity $s=1$ for which we can determine easily the drag force without a frame since we can use a solid plate with the same surface area and thickness.

Note that the value of the coupling is not independent of the solidity: indeed, if we assume a monotone dependency, then the lower the solidity, the lower the coupling will be, until it reaches a zero value at zero solidity. These two points allow us to estimate the coupling.

The frame has a thickness of 0.2 mm each side of the system. Therefore, to quantify the influence of this extrusion on the drag force, we measured for the solidity $s=1$, the drag force $F_r$ of a plate in figure 23(d) and the drag force $F_t$ of the system composed of the frame and the porous screen in figure 23(c), with the same width at the border. The drag forces are shown in table 5 and in figure 24.

We notice that correcting our data with a constant coupling (interference) drag term would not change the curve shape. Doing so with a linear coupling would have a minor effect. Indeed, if the use of a frame to stretch the textiles seems to underestimate the drag coefficient, the difference for the worst case in the drag coefficient between the screen measured directly ($0.993$) and the screen measured with the frame after subtraction of its drag force ($0.939$) is closed to the order of the uncertainty calculated in the next section. Even if we added the coupling drag force estimated for solidity $1$ in all our data, this would not change our conclusions. Therefore, it is reasonable to neglect the interference drag, and simply proceed with the subtraction of the frame drag force from the total drag force.

C.3. Determination of the drag coefficient

To determine the drag coefficient, we calculate a nonlinear regression of the corrected data with the method of the least squares. The model function is

(C2)\begin{equation} {F_p}_i = f(v_i,C_D) = \tfrac{1}{2}\rho S_p C_D v_i^2, \end{equation}

where $C_D$ is the adjustable parameter. We minimize the sum of the square residuals $S$,

(C3)\begin{equation} S = \sum_{i = 1}^n (y_i - f(v_i,C_D))^2. \end{equation}

This leads to

(C4)\begin{equation} C_D = \frac{2}{\rho S}\frac{\displaystyle\sum_{i=1}^n {F_p}_i v_i^2}{\displaystyle\sum_{i=1}^n v_i^4}. \end{equation}

C.4. Measurement uncertainty

In what follows, we estimate the measurement uncertainty $u$ of the physical quantities. We assume, for simplicity, that the parameters $\rho$, ${F_p}_i$, $v_i$ and $S$, for $i \in [\kern-1pt[ 1,n ]\kern-1pt]$ are mutually independent, and that their respective uncertainty is small compared with their value. We neglect the uncertainty on the velocity, then the uncertainty can be calculated with

(C5)\begin{equation} u^2(C_D) = \frac{4}{\rho^2 S^2 \left(\displaystyle\sum_{i=1}^n v_i^4\right)^2} \left(\left(\displaystyle\sum_{i=1}^n {F_p}_i v_i^2\right)^2\left(\frac{u^2(\rho)}{\rho^2}+ \frac{u^2(S)}{S^2}\right)+\displaystyle\sum_{i=1}^n v_i^4 u^2({F_p}_i)\right). \end{equation}

We assume that the force uncertainty is the same for all the data, this assumption is supported by the different repeated measurements we performed for several porous screens, we take the mean value of the deviation we obtained for the force uncertainty.

C.4.1. Solidity uncertainty

The solidity of the porous screens with a regular nylon woven mesh are determined using an image of the screen at the scale of a hundred meshes, taken with a microscope. Several image analyses are used to estimate the solidity and the associated uncertainty. The uncertainty is estimated to ${\pm }0.02$.

C.4.2. Parameter and drag force uncertainties

The uncertainty of the drag force arises from the error of the force balance. By repeating the measurement several times, for different porous screens, we estimated the drag force uncertainty to approximately ${\pm }0.02$ N. The uncertainty of the air density is ${\pm }0.005\,\textrm {kg}\,\textrm {m}^{-3}$. Finally, the uncertainty on the surface of the screen is estimated to ${\pm }4\times 10^{-6}\,\textrm {m}^2$.

References

Ando, S., Nishikawa, M., Kaneda, M. & Suga, K. 2022 Numerical simulation of filtration processes in the flow-induced deformation of fibrous porous media by a three-dimensional two-way fluid–structure interaction scheme. Chem. Engng Sci. 252, 117500.CrossRefGoogle Scholar
Apelt, C.J. & West, G.S. 1974 The effects of wake splitter plates on bluff-body flow in the range $10^{4}< Re<5 \times 10^4$. Part 2. J. Fluid Mech. 71, 145160.CrossRefGoogle Scholar
Ayati, A.A., Steiros, K., Miller, M.A., Duvvuri, S. & Hultmark, M. 2019 A double-multiple streamtube model for vertical axis wind turbines of arbitrary rotor loading. Wind Energy Sci. 4 (4), 653662.CrossRefGoogle Scholar
Bailey, B.J., Montero, J.I., Parra, J.P., Robertson, A.P., Baeza, E. & Kamaruddin, R. 2003 Airflow resistance of greenhouse ventilators with and without insect screens. Biosyst. Engng 86 (2), 217229.CrossRefGoogle Scholar
Blevins, R.D. 1992 Applied Fluid Dynamics Handbook. Krieger Publishing.Google Scholar
Bourrianne, P., Xue, N., Nunes, J., Abkarian, M. & Stone, H.A. 2021 Quantifying the effect of a mask on expiratory flows. Phys. Rev. Fluids 6 (11), 110511.CrossRefGoogle Scholar
de Bray, B.G. 1957 Low speed wind tunnel tests on perforated square flat plates normal to the airstream: drag and velocity fluctuation measurements. Current Paper 323. Ministry of Supply, Aeronautical Research Council.Google Scholar
Brundrett, E. 1993 Prediction of pressure drop for incompressible flow through screens. J. Fluids Engng 115 (2), 239–242.CrossRefGoogle Scholar
Carvajal, D., Silva-Llanca, L., Larraguibel, D. & González, B. 2020 On the aerodynamic fog collection efficiency of fog water collectors via three-dimensional numerical simulations. Atmos. Res. 245, 105123.CrossRefGoogle Scholar
Cheer, A.Y.L. & Koehl, M.A.R. 1987 Paddles and rakes: fluid flow through bristled appendages of small organisms. J. Theor. Biol. 129, 1739.CrossRefGoogle Scholar
Crompton, M.J. & Barrett, R.V. 2000 Investigation of the separation bubble formed behind the sharp leading edge of a flat plate at incidence. Proc. Inst. Mech. Engrs G 214 (3), 157176.CrossRefGoogle Scholar
de Dios Rivera, J. 2011 Aerodynamic collection efficiency of flog collectors. Atmos. Res. 102, 335342.CrossRefGoogle Scholar
Du, W., Iacoviello, F., Fernandez, T., Loureiro, R., Brett, D.J.L. & Shearing, P.R. 2021 Microstructure analysis and image-based modelling of face masks for Covid-19 virus protection. Commun. Mater. 2 (69), 217229.CrossRefGoogle Scholar
Eckert, B. & Pflüger, F. 1942 The resistance coefficient of commercial round wire grids. Technical Memorandum 1003. National Advisory Committee for Aeronautics.Google Scholar
Fail, R., Lawford, J.A. & Eyre, R.C.W. 1957 Low-speed-experiments on the wake characteristics of flat plates normal to an air stream. Rep. and Memoranda 3120. Aeronautical Research Council.Google Scholar
Giannoulis, A., Stathopoulos, T., Briassoulis, D. & Mistriotis, A. 2012 Wind loading on vertical panels with different permeabilities. J. Wind Engng Ind. Aerodyn. 107–108, 116.CrossRefGoogle Scholar
Gunter, N.M. 1967 Potential Theory and its Application to Basic Problems of Mathematical Physics. Frederick Ungar.Google Scholar
Hoerner, S.F. 1952 Aerodynamic properties of screens and fabrics. Textile Res. J. 22 (4), 274280.CrossRefGoogle Scholar
Hoerner, S.F. 1965 Fluid Dynamic Drag. Published by the Author.Google Scholar
Hood, K., Jammalamadaka, M.S.S. & Hosoi, A.E. 2019 Marine crustaceans with hairy appendages: role of hydrodynamic boundary layers in sensing and feeding. Phys. Rev. Fluids 4, 114102.CrossRefGoogle Scholar
Ito, A. & Garry, K.P. 1998 Pressure measurements around a two-dimensional gauze at incidence. J. Fluids Struct. 12, 171181.CrossRefGoogle Scholar
Johari, H. & Desabrais, K.J. 2005 Vortex shedding in the near wake of a parachute canopy. J. Fluid Mech. 536, 185207.CrossRefGoogle Scholar
Kalugin, V.T., Epikhin, A.S., Chernukha, P.A. & Kalugina, M.D. 2021 The effect of perforation on aerodynamic characteristics and the vortex flow field around a flat plate. IOP Conf. Ser. 1191, 1–7.CrossRefGoogle Scholar
Koo, J.-K. & James, D.F. 1973 Fluid flow around and through a screen. J. Fluid Mech. 60, 513538.CrossRefGoogle Scholar
Kress, R. 1999 Linear Integral Equations, Applied Mathematical Sciences, vol. 82. Springer.CrossRefGoogle Scholar
Laws, E.M. & Livesey, J.L. 1978 Flow through screens. Ann. Rev. Fluid Mech. 10, 247266.CrossRefGoogle Scholar
Ledda, P.G., Boujo, E., Camarri, S., Gallaire, F. & Zampogna, G.A. 2021 Homogenization-based design of microstructured membranes: wake flow past permeable shells. J. Fluid Mech. 927, A31.CrossRefGoogle Scholar
Ledda, P.G., Siconolfi, L., Viola, F., Gallaire, F. & Camarri, S. 2018 Suppression of von Kármán vortex streets past porous rectangular cylinders. Phys. Rev. Fluids 3, 103901.CrossRefGoogle Scholar
Letchford, C.W. 2001 Wind loads on rectangular signboards and hoardings. J. Wind Engng Ind. Aerodyn. 89, 135151.CrossRefGoogle Scholar
Mittal, R., Ni, R. & Seo, J.-H. 2020 The flow physics of Covid-19. J. Fluid Mech. 894, F2.CrossRefGoogle Scholar
Moncuquet, A., Mitranescu, A., Marchand, O.C., Ramananarivo, S. & Duprat, C. 2022 Collecting fog with vertical fibres: combined laboratory and in-situ study. Atmos. Res. 277, 106312.CrossRefGoogle Scholar
Monjezi, M. & Jamaati, H. 2021 The effects of face mask specifications on work of breathing and particle filtration efficiency. Med. Engng Phys. 98, 3643.CrossRefGoogle ScholarPubMed
Okamoto, M. & Azuma, A. 2011 Aerodynamic characteristics at low Reynolds numbers for wings of various planforms. AIAA J. 49 (6), 1135–1150.CrossRefGoogle Scholar
O'Neill, F.G. 2006 Source models of flow through and around screens and gauzes. Ocean Engng 33, 18841895.CrossRefGoogle Scholar
Parkinson, G.V. & Jandali, T. 1970 A wake source model for bluff body potential flow. J. Fluid Mech. 40, 577594.CrossRefGoogle Scholar
Pezzulla, M., Strong, E.F., Gallaire, F. & Reis, P.M. 2020 Deformation of porous flexible strip in low and moderate Reynolds number flows. Phys. Rev. Fluids 5, 084103.CrossRefGoogle Scholar
Pinker, R.A. & Herbert, M.V. 1967 Pressure loss associated with compressible flow through square-mesh wire gauzes. J. Mech. Engng Sci. 9 (1), 1123.CrossRefGoogle Scholar
Pomaranzi, G., Daniotti, N., Schito, P., Rosa, L. & Zasso, A. 2020 Experimental assessment of the effects of a porous double skin façade system on cladding loads. J. Wind Engng Ind. Aerodyn. 196, 104019.CrossRefGoogle Scholar
Pontavice, E.D. 2016 Propulsion par cerf-volant: envol et pérégrinations. PhD thesis, Université Paris-Saclay.Google Scholar
Prandtl, L. & Flachsbart, O. 1932 Widerstand von seidengazefiltern, runddraht- und bledistreifensieben mit quadratischen maschen. Ergebnisse der Aerodynamischen Versuchsanstalt zu Göttingen.Google Scholar
Pressley, A. 2010 Elementary Differential Geometry. Springer.CrossRefGoogle Scholar
Regalado, C.M. & Ritter, A. 2016 The design of an optimal fog water collector: a theoretical analysis. Atmos. Res. 178–179, 4554.CrossRefGoogle Scholar
Reynolds, A.J. 1969 Flow deflection by gauze screen. J. Mech. Engng Sci. 11 (65), 290–294.CrossRefGoogle Scholar
Roshko, A. 1954 A new hodograph for free-streamline theory. Tech. Note 3168. National Advisory Committee for Aeronautics.Google Scholar
Sarpkaya, T. & Lindsey, P.J. 1990 Unsteady flow about porous cambered shells. J. Aircraft 28 (8), 502–508.Google Scholar
Schubauer, G.B., Spangenberg, W.G. & Klebanoff, P.S. 1950 Aerodynamic characteristics of damping screens. National Advisory Committee for Aeronautics Collection 2001.Google Scholar
Shklyar, A. & Arbel, A. 2008 Numerical simulations of turbulent flow through screen mesh by transient SST and $k-\epsilon$ turbulent models. Acta Hortic. 802, 6166.CrossRefGoogle Scholar
Simmons, M.A. & Cowdrey, C.F. 1945 Measurement of the aerodynamic forces acting on pours screens. Rep. and Memoranda 2276. National Advisory Committee for Aeronautics.Google Scholar
Steiros, K., Bempedelis, N. & Cicolin, M.M. 2022 An analytical blockage correction model for high-solidity turbines. J. Fluid Mech. 948, A57.CrossRefGoogle Scholar
Steiros, K. & Hultmark, M. 2018 Drag on flat plates of arbitrary porosity. J. Fluid Mech. 853, 111.CrossRefGoogle Scholar
Taylor, G.I. 1944 Air resistance of a flat plate of very porous material. Rep. and Memoranda 2236. Aeronautical Research Council.Google Scholar
Taylor, G.I. & Davies, R.M. 1944 The aerodynamics of porous sheets. Rep. and Memoranda 2237. Aeronautical Research Council.Google Scholar
Teitel, M. 2010 Using computational fluid dynamics simulations to determine pressure drops on woven screens. Biosyst. Engng 105, 172179.CrossRefGoogle Scholar
Torres, G.E. & Mueller, T.J. 2004 Low-aspect-ratio wing aerodynamics at low Reynolds numbers. AIAA J. 42 (5), 865–873.CrossRefGoogle Scholar
Wang, Y., Yang, G., Huang, Y., Huang, Y., Zhuan, R. & Wu, J. 2021 Analytical model of flow-through-screen pressure drop for metal wire screens considering the effects of pore structures. Chem. Engng Sci. 229, 116037.CrossRefGoogle Scholar
Wieghardt, K.E.G. 1953 On the resistance of screens. Aeronaut. Q. 4, 186192.CrossRefGoogle Scholar
Wu, T.Y.-T. 1962 A wake model for free-streamline flow theory. Part 1. Fully and partially developed wake flows and cavity flows past an oblique flat plate. J. Fluid Mech. 13, 161192.CrossRefGoogle Scholar
Xu, M., Patruno, L., Lo, Y.-L. & de Miranda, S. 2020 On the use of the pressure jump approach for the simulation of separated external flows around porous structures:a forward facing step. J. Wind Engng Ind. Aerodyn. 207, 104377.CrossRefGoogle Scholar
Yih, C.-S. 1957 Stream functions in three-dimensional flows. La Houille Blanche 43, 439450.CrossRefGoogle Scholar
Zampogna, G.A. & Gallaire, F. 2020 Effective stress jump across membranes. J. Fluid Mech. 892, A9.CrossRefGoogle Scholar
Figure 0

Table 1. Porous screen characteristics. The interval of the Reynolds number $Re_d$ is calculated with the fibre diameter $d$ as characteristic size, with a velocity $v_0 = 5.0\,{\rm m}\,{\rm s}^{-1}$ for the minimum value, a velocity $v_0 = 13.0\,{\rm m}\,{\rm s}^{-1}$ for the maximum value and a kinematic viscosity $\nu = 15.6 \times 10^{-6}\,{\rm m}^2\,{\rm s}^{-1}$. For the screens P2, P12, P26, P27, P28 P29 and P30, we put the detailed values of the drag coefficient at each velocity (and calculated the corresponding Reynolds number $Re_d$). The uncertainties as well as the fibre diameters, material type and geometry can be found in the Appendix C.

Figure 1

Figure 1. Diagram of the experimental set-up for the drag coefficient measurement, and examples of porous screens essentially in nylon woven mesh.

Figure 2

Figure 2. (a) Drag force $(F_D)$ measured for different square porous screens at normal incidence. (b) Drag coefficient $C_D$ as a function of velocity $v_0$, renormalised by the drag coefficient at maximum velocity for each screen $C_{D}(v_{0_{max}})$, including screens showing the higher variation of the drag coefficient with respect to the velocity. Colour bars: fibre diameter $d$.

Figure 3

Figure 3. Drag coefficient $C_{D_{max}}$ as a function of the solidity for various square porous screens normal to the free flow and comparison with the 2-D model of Steiros & Hultmark (2018). The colour of the points indicates the value of the fibre diameter $d$ following the colour scale on the right. The data from Prandtl & Flachsbart (1932) and the value at solidity $s=1$ from Blevins (1992) are also plotted.

Figure 4

Figure 4. Experimental evidence of the deviation of the flow and streamlines around a square porous screen at different solidities $s$ under normal uniform upstream flow $v_0 = 2.84\pm 0.02\,{\rm m}\,{\rm s}^{-1}$. Top view of the screen. The white bar (110 mm) corresponds to the screen embedded in the frame, the screen alone is 100 mm long. Results are shown for (a) P14, $s=0.24$; (b) P6, $s=0.61$; (c) P23, $s=0.82$.

Figure 5

Figure 5. Experimental evidence of the 3-D deviation of the flow and streamlines around a square porous screen (transverse and diagonal directions) for a screen of solidity $s=0.61$ (P6) under normal uniform upstream flow $v_0 = 2.84\pm 0.02\,{\rm m}\,{\rm s}^{-1}$. Here $D_{sd}/2 = 55.8 \pm 1.5$ mm and $D_{st}/2 = 32.8 \pm 1.2$ mm; thus, $D_{sd}/D_{st} = 1.7 \pm 0.1 >\sqrt {2}$.

Figure 6

Figure 6. Diagram of the theoretical model derived in this study and the ingredients from previous models used. In blue, a complete 3-D potential model. In grey, the theoretical basis used for the application to rectangular mesh screens. In purple, the application to rectangular mesh screens. Here $\phi$ denotes the velocity potential, $p_{III}$ denotes the base pressure and $p_{I}$$p_{II}$ denotes the pressure jump across the screen.

Figure 7

Figure 7. Diagram of the model for a 3-D potential flow around and through a porous screen. The dashed lines are the separation streamlines used as a boundary between the regions. The dotted lines are the streamlines used in the model to calculate the velocities in regions I and II. The incoming flow is laminar and is extended over the entire height of the system. Here $\mathcal {C}_w$ denotes the section of the wake.

Figure 8

Figure 8. Diagram of the inclined rectangular porous screen in the $(Oyz)$ plane in a 3-D free flow with an angle $\beta$. The dashed lines are the separation streamlines used as a boundary between the regions of the model. The dotted line between points A and B is the streamline used in the model to calculate the velocities in the regions I and III. The incoming flow is laminar.

Figure 9

Figure 9. Theoretical prediction of the drag coefficient and velocities on the surface as a function of the solidity obtained with the 3-D model applied to a square homogeneous screen normal to the flow ($\beta = {{\rm \pi} }/{2}$), and comparison with different 2-D models. The full lines are obtained with the 3-D model developed in this paper using the same pressure jump law across the screen as Steiros & Hultmark (2018) (but taking account of geometric 3-D effects). (a) Drag coefficient. (b) Normal velocity and quadratic mean of the tangential velocity as defined respectively in (3.38) and (3.45). (a) Theoretical drag coefficient. (b) Theoretical velocity at the screen.

Figure 10

Figure 10. Drag coefficient as a function of the solidity for various square porous screens normal to the free flow with different $Re_d$, and comparison between the 3-D model. The full lines are obtained with the 3-D model developed in this paper, solution of (4.6). The black dash-dotted curve is the solution of (3.50), thus, without taking into account $Re_d$.

Figure 11

Figure 11. Drag coefficient as a function of the Reynolds number $Re_d$ for three different narrow ranges of solidities: $s=0.1 \pm 0.02$, $s=0.28 \pm 0.04$, $s=0.585 \pm 0.025$. The bullets represent the experimental measurements, and the plain lines the results of our model ((4.6) taking into account the effect of mesh Reynolds number $Re_d$). The data points for the drag as well as the Reynolds number values are obtained for the velocity $5.0\leq v_0\leq 13\,{\rm m}\,{\rm s}^{-1}$. The dash-dotted lines represent the value of the drag coefficient for a solid plate (three dimensional) from our experiments and from Blevins (1992).

Figure 12

Figure 12. Streamlines and velocity magnitude obtained with the 3-D model for a square screen normal to the flow at $v_0=2.0\,{\rm m}\,{\rm s}^{-1}$ for different solidities and the pressure jump defined in (3.18) (thus, without dependency on the mesh Reynolds number $Re_d$). Results are shown for (a) $s=0.3$ and $\beta = 90^{\circ }$; (b) $s=0.6$ and $\beta = 90^{\circ }$; (c) $s=0.9$ and $\beta = 90^{\circ }$.

Figure 13

Figure 13. Comparison between the 3-D theoretical model and our experiments for the prediction of the velocity field in the $(z,x)$ plane with $y=0$ for three different solidities (low, moderate and high). The theoretical velocity field has been obtained using the theoretical pressure jump law (3.18) and (3.19). Uncertainties on the velocities are estimated to be around $0.1\,{\rm m}\,{\rm s}^{-1}$. Results are shown for (a) P14, $s=0.24$, experimental; (b) P14, $s=0.24$, theoretical; (c) P6, $s=0.61$, experimental; (d) P6, $s=0.61$, theoretical; (e) P23, $s=0.82$, experimental; (f) P23, $s=0.82$, theoretical.

Figure 14

Figure 14. Dimensionless normal velocity ${v_{I_n}}/{v_0}$ obtained with (4.6) for the solid lines, and with (3.50) for the dotted black line. The bullets represent the experimental measurements obtained using a hot wire anemometer.

Figure 15

Table 2. Reynolds number $Re_d$ calculated with the fibre diameter $d$ as characteristic size and with a velocity $v_0 = 2.84\,{\rm m}\,{\rm s}^{-1}$ and a kinematic viscosity $\nu = 15.6 \times 10^{-6}\,{\rm m}^2\,{\rm s}^{-1}$.

Figure 16

Figure 15. Dimensionless tangential velocity magnitude ${|v_{I_t}|}/{v_0}$ on the screen surface obtained with the 3-D model for a square screen normal to a flow for different solidities and Reynolds numbers ((4.6) taking into account the effect of mesh Reynolds number $Re_d$). Results are shown for (ac) $Re_d = 10$; (df) $Re_d = 100$; (a) $s=0.1$ and $Re_d = 10$; (b) $s=0.5$ and $Re_d = 10$; (c) $s=0.8$ and $Re_d = 10$; (d) $s=0.1$ and $Re_d = 100$; (e) $s=0.5$ and $Re_d = 100$; (f) $s=0.8$ and $Re_d = 100$.

Figure 17

Figure 16. Upstream section $S_a$ of the streamtube containing all the fluid passing through the surface, at $z=-\infty$, seen from the front, for $Re_d = 10$ and $Re_d = 100$.

Figure 18

Figure 17. Streamlines and velocity magnitude obtained with the 3-D model for a square screen inclined to a flow at $v_0=2.0\,{\rm m}\,{\rm s}^{-1}$ for different solidities and the pressure jump defined in (3.18) (thus, without dependency on the mesh Reynolds number $Re_d$). Results are shown for (a) $s=0.8$ and $\beta = 60^{\circ }$; (b) $s=0.8$ and $\beta = 30^{\circ }$; (c) $s=1$ and $\beta = 30^{\circ }$.

Figure 19

Figure 18. Comparison between the 3-D model and our experiments for the prediction of the drag coefficient at two angles of inclination for various porous square screens: (a) $\beta = 65^\circ$ (b) $\beta = 43^\circ$ using (3.50). The black dash-dotted line is the theoretical result obtained with (4.6).

Figure 20

Figure 19. Diagram of the flow at the scale of the pores in case of parallel fibres (left), and corresponding imaginary surface as used in the model (right).

Figure 21

Figure 20. Comparison between the drag coefficient prediction of the 3-D model and the experimental measurement for a square plate in a free flow for different angles of inclination. Data are added from Blevins (1992), Okamoto & Azuma (2011) and Torres & Mueller (2004).

Figure 22

Figure 21. Comparison of the drag force $F_D$ with and without a splitter plate for a solid plate (without frame support) of dimensions $11.0\times 11.0\,{\rm cm}^2$ and thickness 3 mm. The fitting curves are obtained using a quadratic law.

Figure 23

Table 3. Shape parameter for rectangular plate with different aspect ratios. Here $l_1$ and $l_2$ are the lengths of the sides of the rectangular plate. The aspect ratio is ${l_2}/{l_1}$.

Figure 24

Table 4. Porous screen characteristics. The Reynolds number $Re_d$ is calculated with a velocity $v_0 = 10\,\textrm {m}\,\textrm {s}^{-1}$ and a kinematic viscosity $\nu = 15.6 \times 10^{-6}\,\textrm {m}^2\,\textrm {s}^{-1}$, except for the screens P27, P28 and P29 where the velocity is $v_0 = 5\,\textrm {m}\,\textrm {s}^{-1}$.

Figure 25

Table 5. Summary of the values of the coupling drag force for a flat plate (solidity $s=1$) for different velocities. All the values have an uncertainty of approximately 0.020 N. For the measurement of $F_p$ and $F_r$, a solid flat plate with 4.0 mm thickness has been used, corresponding to the thickness of the frame used for measuring $F_c$ and $F_t$. The drag coefficients are obtained using the dimensions shown in figure 23 with fluid density estimated from the measurement of temperature, pressure and humidity during the different experiments.

Figure 26

Figure 22. Diagram of the fibre mesh.

Figure 27

Figure 23. Diagram of the screens (in dark gold) and the frame support (in green).

Figure 28

Figure 24. Drag force for the plates and the frame represented in figure 23. Values and drag coefficients are in table 5. Data are fitted using a quadratic law with respect to the velocity $v_0$. The surface is orthogonal to the mean far-field flow direction $\boldsymbol {v}_0$.